Comparative transcriptome analysis of eggplant (Solanum melongena L.) and turkey berry (Solanum torvum Sw.): phylogenomics and disease resistance analysis
- Xu Yang†1Email author,
- Yu-Fu Cheng†1,
- Cao Deng2,
- Yan Ma1,
- Zhi-Wen Wang2,
- Xue-Hao Chen1 and
- Lin-Bao Xue1
© Yang et al.; licensee BioMed Central Ltd. 2014
Received: 17 February 2014
Accepted: 22 May 2014
Published: 31 May 2014
Eggplant (Solanum melongena L.) and turkey berry (S. torvum Sw.), a wild ally of eggplant with promising multi-disease resistance traits, are of great economic, medicinal and genetic importance, but genomic resources for these species are lacking. In the present study, we sequenced the transcriptomes of eggplant and turkey berry to accelerate research on these two non-model species.
We built comprehensive, high-quality de novo transcriptome assemblies of the two Leptostemonum clade Solanum species from short-read RNA-Sequencing data. We obtained 34,174 unigenes for eggplant and 38,185 unigenes for turkey berry. Functional annotations based on sequence similarity to known plant datasets revealed a distribution of functional categories for both species very similar to that of tomato. Comparison of eggplant, turkey berry and another 11 plant proteomes resulted in 276 high-confidence single-copy orthologous groups, reasonable phylogenetic tree inferences and reliable divergence time estimations. From these data, it appears that eggplant and its wild Leptostemonum clade relative turkey berry split from each other in the late Miocene, ~6.66 million years ago, and that Leptostemonum split from the Potatoe clade in the middle Miocene, ~15.75 million years ago. Furthermore, 621 and 815 plant resistance genes were identified in eggplant and turkey berry respectively, indicating the variation of disease resistance genes between them.
This study provides a comprehensive transcriptome resource for two Leptostemonum clade Solanum species and insight into their evolutionary history and biological characteristics. These resources establish a foundation for further investigations of eggplant biology and for agricultural improvement of this important vegetable. More generally, we show that RNA-Seq is a fast, reliable and cost-effective method for assessing genome evolution in non-model species.
KeywordsSolanum torvum Sw. Solanum melongena L. Comparative transcriptomics Evolution Plant resistance genes
Eggplant (Solanum melongena L.) is the third most agriculturally important crop from the genus Solanum after potato (S. tuberosum)  and tomato (S. lycopersicum) . This large and diverse genus of flowering plants comprises >1400 species having a wide range of genetic and phenotypic variation . In 2011, 46.8 million tons of eggplant was produced in the top four producing countries, namely China (27.7 million tons), India (11.8 million tons), Egypt (1.1 million tons) and Turkey (8.2 million tons), according to the Food and Agriculture Organization of the United Nations (http://faostat.fao.org). There are three closely related cultivated species of eggplant, all of Old World origin: S. aethiopicum L. (scarlet eggplant), S. macrocarpon L. (gboma eggplant) and S. melongena L. (brinjal or aubergine eggplant) . The brinjal or aubergine eggplant, hereafter referred to as eggplant, is cultivated worldwide and is an autogamous diploid with 12 chromosomes (2n = 2x = 24) . Eggplant is susceptible to many bacterial and fungal pathogens and insects, such as the Verticillium dahlia fungus and nematodes , which cause significant yield losses. As such, improving resistance to biotic and abiotic stresses is one of the main objectives of eggplant breeding programs.
Solanum torvum Sw., commonly known as turkey berry, is a wild relative of eggplant and is found in tropical Africa, Asia and South America. Turkey berry is widely consumed and is an important folk medicinal plant in tropical and subtropical countries . More importantly, turkey berry is resistant to root-knot nematodes and the most serious soil-borne diseases, such as those caused by Ralstonia solanacearum, V. dahlia Klebahn and Fusarium oxysporum f. sp. Melongenae , providing promising genetic resources for improvement of eggplant. Traditional grafting techniques are now used worldwide in eggplant cultivation, in which eggplant tissues are grafted onto disease-resistant rootstock of turkey berry [8–10]. Also, attempts have been made to introduce turkey berry resistance into eggplant through conventional breeding and biotechnological techniques, however, progress is limited. Owing to sexual incompatibilities, however, attempts at crossing eggplant with turkey berry have had limited success , and sterile hybrids were obtained, with difficulty, only when eggplant was used as the female parent . Other biotechnological techniques, such as embryo rescue, somatic hybridization and Agrobacterium-mediated transformation, have been difficult to apply to eggplant [12, 13] because of the limited genetic information available for this species.
Solanum crops that belong to the Potatoe clade, which includes potato and tomato, have been targets for comprehensive genomic studies [1, 2]. However, genomic resources are lacking for the Leptostemonum clade (the “spiny solanums”), which comprises almost one-third of the genus distributed worldwide  and includes eggplant and turkey berry. For eggplant, 98,861 nucleotide sequences have been deposited in the National Center for Biotechnology Information (NCBI) GenBank database (as of December 18, 2013), and the vast majority of them (98,086) were provided recently by a comparative analysis of ESTs . In that analysis, however, only 16,245 unigenes were constructed, which is approximately half the number of genes identified in the closely related potato (39,031)  and tomato (34,727) , implying that these unigenes represent only a limited portion of the whole eggplant transcriptome. In addition, large numbers of short-read sequences have been generated from turkey berry in attempts to identify single nucleotide polymorphisms and simple sequence repeats using restriction site–associated DNA tag sequencing strategies; however, this approach provides only limited information on full-length genes, and such information is vital for identifying trait-related genes and for quantitative gene expression analysis. Recent studies reported 6,296 unigenes from S. torvum cultivar Torubamubiga  and 36,797 unigenes from S. torvum Sw. accession TG1 transcriptome assemblies . In the latter study, however, sequencing was confined to the 3′ end of the transcripts, resulting in fragmentary assembled transcripts as revealed by an N50 value (the 50% of the entire assembly is contained in sequences equal to or larger than this value) of only 514 bp and an N10 value of only 715 bp. Therefore, there is an urgent need to obtain more high-quality genomic information about eggplant and turkey berry, and a promising technology to accomplish this is RNA sequencing (RNA-Seq).
High-quality transcriptome data would not only facilitate genetic and molecular breeding approaches in eggplant and allow genomic resource mining in turkey berry but also be valuable for comparative biology studies, such as phylogenomics. For example, RNA-Seq data have been used to explore the evolution of paleopolyploidy in plants [17, 18] and to reconstruct deep phylogenies in flowering plants of the grape family (Vitaceae) . These studies suggest that transcriptome data can be very useful and practical in the reconstruction of phylogenies in flowering plants.
The specific goals of this study were to (1) generate high-quality transcripts and unigenes of eggplant and turkey berry using RNA-Seq, which will provide reference transcriptomes for further analysis, such as trait-related gene mining and quantitative expression analysis; (2) produce a dated phylogeny of the Potatoe and Leptostemonum clades and of the Leptostemonum-nested eggplant (Old Word clade) and turkey berry (Torva clade), which will deepen our understanding of phylogenetic relationships and ultimately assist crop improvement; and (3) identify and compare disease resistance genes in eggplant and turkey berry to take a first glance at the variation of resistance genes among them using RNA-Seq data.
Results and discussion
De novotranscriptome assembly and annotation captures high-quality transcripts and unigenes
Summary of the eggplant and turkey berry transcriptome assemblies
Total raw reads
27,387,245 × 2
15,576,018 × 2
72 + 72
72 + 72
Total raw reads data size (bp)
To evaluate the completeness of our assemblies, the transcripts and unigenes were aligned with the tomato and potato sequences to obtain the corresponding reference genes, and then the unigene and transcript distributions were plotted against the tomato and potato reference genomes. The unigene and transcript distribution patterns were similar to the gene distribution patterns of both the tomato (Additional file 1: Figure S2) and potato (Additional file 1: Figure S3), indicating the completeness of the unigene assemblies.
Our assemblies were of substantially higher quality than those generated in previous studies [15, 16]. In a comparative analysis of eggplant ESTs , only 16,245 unigenes were constructed, which is less than half of our 34,174 unigenes and of the genes identified in the closely related potato (39,031)  and tomato (34,727) . Global transcriptome profiling aimed at gaining insight into the mechanisms underpinning turkey berry resistance against Meloidogyne incognita produced 36,797 unigenes from S. torvum Sw. accession TG1. Although this number is comparable to our results, to improve coverage and conserve specificity, sequencing in that study was confined to the 3′ end of the transcripts, resulting in a fragmented assembly, as indicated by low N50 (514 bp) and N10 (715 bp) values. Without introduced bias, our N50 value was 1,349 bp, which is similar to the N50 of the non-redundant coding sequences (CDS) from tomato (1,467 bp) and potato (1,257 bp). Taken together, these results suggest that the quality and completeness of our sequencing and assembly were high enough for annotation and further analyses.
Annotation results of the eggplant and turkey berry unigenes
Comparative analysis of gene sets between plants
Summary of orthologous groups between 13 species
Number of genes
Genes in families
Number of families
Average genes per family
S. melongena L.
S. torvum Sw.
O. sativa japonica
Nevertheless, 4,900 orthologous groups were shared by all 13 species (Figure 4A), which is comparable to previous studies. Wang et al. found 9,525 shared core orthologous groups between Gossypium raimondii, Theobroma cacao, A. thaliana and Z. mays, D’Hont et al. found 7,674 shared gene families between Musa acuminata, Phoenix dactylifera, A. thaliana, O. sativa, Sorghum bicolor and Brachypodium distachyon, and Peng et al.  found 9,451 shared gene families among five grass genomes. The numbers of orthologous groups that we observed were smaller, but the groups included more species, which may indicate that our analysis was more stringent and therefore may represent only highly conserved orthologous groups among dicotyledonous and monocotyledonous plants. Among the 4,900 core orthologous groups, 559 contained only one ortholog in each species (single copy, Figure 4B). These groups were suitable for inferring phylogenetic relationships and for estimating divergence time.
Inferring phylogenetic relationships
To maximize the information content of our sequences and minimize the impact of missing data, the 559 single-copy orthologous groups were further filtered with stricter constraints on length (minimum 200 amino acids) and sequence alignment (maximum missing data 50% in the CDS alignments), and the resultant 276 groups were used for phylogenetic tree reconstruction.
Estimation of divergence time
All of the geological times estimated for nodes leading to non-Asterid species were well matched to data deposited in TimeTree , a public knowledge-base of divergence times among organisms, demonstrating the high reliability of this molecular clock dating strategy. As shown in Figure 6, the divergence between eggplant and turkey berry appears to have occurred ~6.66 (4.9–8.8) million years ago (Mya), during the late Miocene. The Leptostemonum and Potatoe clades shared a common ancestor during the middle Miocene and appear to have diverged ~15.75 (12.7–18.8) Mya, which is in agreement with the 11.60–16.00 Mya estimated by Wang et al.. A whole-genome triplication in tomato  and potato  has been estimated at 71 (±19.4) Mya on the basis of synonymous substitutions of paralogous genes, which is much earlier than the splitting of Leptostemonum and Potatoe clades. This timeline implies, therefore, that both eggplant and turkey berry underwent genome triplication, but this remains to be verified by complete genome sequences.
Disease resistance genes
A fundamental strategy for controlling diseases in agriculturally important plants is the isolation of resistance genes from their less susceptible relatives to be used in conventional breeding, genetic engineering and biotechnological approaches [12, 13]. Because of limited genetic resources for eggplant and turkey berry, however, only one resistance gene, a Ve-like gene (StVe), has been identified in these species, to our knowledge . Moreover, a large number of plant resistance genes have been identified and deposited in the Plant Resistance Genes database (PRGdb, http://prgdb.crg.eu/wiki/Main_Page) . Of these entries, 112 were manually curated to confirm that they were described in the literature to confer resistance to pathogens, and they are grouped into seven distinct classes based on the presence of specific domains or partial domains [34, 35]: N-terminal coiled coil–nucleotide-binding site–leucine-rich repeat (CNL), Toll interleukin1 receptor–nucleotide-binding site–leucine-rich repeat (TNL), receptor-like kinase (RLK), receptor-like protein (RLP), three truncated classes (Kinase, NL and TN) and ‘Other’ which has no typical resistance related domains. Of the 112 entries, 36 (32.14%) are from Solanaceae, 37 (33.04%) are from Poaceae, 25 (22.32%) are from Brassicaceae, and only 14 (12.50%) are from other families. The high percentage of closely related sequences (from Solanaceae) and outgroup sequences (from monocot, Poaceae) made it possible to identify and classify both recently arisen and ancient orthologous resistance genes through homology-based approaches.
Summary of plant resistance genes in Solanum species and Arabidopsis
Compared with Arabidopsis, each of the four Solanum species contained approximately twice the number of resistance genes, with 621 in eggplant, 815 in turkey berry, 505 in tomato, and 774 in potato. The wide intraspecific variation in number of resistance genes may underlie the species-specific differences in resistance to different types and quantities of pathogens and differences in the degree of responses to the same pathogen. The different resistance capability between eggplant and turkey berry may partly result from variation in the number of resistance genes, as turkey berry carries nearly 200 more resistance genes than eggplant. Resistance genes are frequently clustered in the genome—the result of both segmental and tandem duplications [36, 37]—and this was also observed in tomato (Figure 1, B1 ring) and potato (Additional file 1: Figure S1, B1 ring). Resistance genes also appeared to be clustered in eggplant (Figure 1, B2 ring and Additional file 1: Figure S1, B2 ring) and turkey berry (Figure 1, B3 ring and Additional file 1: Figure S1, B3 ring), but this observation needs verification with genome data.
Another difference between the Solanum species and Arabidopsis was the composition of resistance gene classes. TNL genes outnumbered CNL genes in the four Solanum species, which is similar to what has been observed in both grape and poplar (P. trichocarpa) but in contrast to what has been found in apple (Malus domestica), soybean (Glycine max) and Arabidopsis . The CNL and TNL classes are the two major NL proteins, which are believed to act intracellularly , and the RLK and RLP classes are the two major membrane-localized receptor proteins that sense various pathogens and transduce signals to downstream intra- and intercellular networks . The numbers of genes of all of these four classes were larger in turkey berry than in eggplant (Table 4). This may reflect amplification of the entire disease resistance pathway in turkey berry rather than duplication of a particular gene or class of genes to enhance pathogen defense and consequently improve fitness. The variation in the number of resistance genes was also evidenced by plotting the distribution of eggplant and turkey berry resistance genes against the tomato genome (Figure 1 B2 and B3 rings). As shown in Figure 1, the distribution patterns were similar (presence or absence) overall, but numbers of genes varied.
Our results deepen our understanding of phylogenetic relationships, which will ultimately assist in eggplant improvement efforts. Furthermore, these high-quality unigenes will be useful in trait-related gene mining, as we demonstrated with the identification of plant resistance genes and comparison of these genes between species. Results from resistance genes identification indicated the high variation of resistance genes between them. In addition, these datasets can serve as reference transcriptomes for further analyses, such as quantitative gene expression profiling, to broaden our understanding of eggplant biology and to improve this agriculturally important vegetable.
None of the species used in this study are endangered or protected, and all plants were grown in greenhouses, which complies with all relevant regulations. Therefore, no specific permits were required for the collection of samples.
Plant materials and transcriptome sequencing
All samples of eggplant and turkey berry were collected from the experimental farm of the Department of Horticulture in Yangzhou University, Jiangsu Province, and were grown in pots containing peat, vermiculite and perlite (3:1:1, v/v) in a greenhouse at 28/18°C (12/12 h) day/night temperature with relative humidity ranging 70%–85%. For each species, the following tissues were sampled from seedling at the four true leaves stage: root, stem and young leaves. All samples were immediately frozen in liquid nitrogen and stored at −70°C for later use. The RNA extraction, library construction and RNA-Seq were performed at Beijing BioMarker Technologies (Beijing, China) following the protocol of Han et al..
Sequence data analysis and assembly
To obtain high-quality clean reads for transcript de novo assembly, the raw reads from transcriptome sequencing were filtered with the following criteria: (1) reads with adaptor contamination were removed, (2) low-quality reads were designated with “N” and (3) reads in which >10% of the bases had a Q-value < 20 were discarded. The clean reads were then assembled into contigs using Trinity  (http://trinityrnaseq.sourceforge.net/) with an optimized k-mer length of 31 for de novo assembly. Based on the paired-end information, the contigs (longer than 47 bp) were linked into transcripts. Finally, to eliminate redundant sequences, transcripts longer than 200 bp were clustered based on sequence similarities, and the longest transcript in each cluster represented the final assembled unigene that was subjected to functional and structural annotation.
Evaluation of sequence and assembly completeness
Using TopHat2  with default parameters, the clean sequencing reads from eggplant and turkey berry were aligned to the tomato and potato genomes. Tomato (ITAG2.3 release) and potato (PGSC DM 3.4 release) data were obtained from Sol Genomics Network (http://solgenomics.net/). The resultant accepted bam files were assessed for call depth at each nucleotide site using SAMtools , and the depth distribution was plotted for eggplant and turkey berry relative to the tomato and potato genomes.
The corresponding tomato and potato homologs of transcripts and unigenes of the eggplant and turkey berry were identified using BLASTX. Transcripts and unigenes were aligned with the parameters: −evalue 1e-5 -outfmt 6 -max_target_seqs 1 -seg no, and then the alignments were filtered for minimum alignment length of 50 amino acids and identity value of ≥30%. The distributions of eggplant and turkey berry unigenes and transcripts relative to the tomato and potato genomes were then plotted.
Functional and structural annotation
To determine the functional categories of the unigenes, a BLASTX search with a cut-off E-value ≤ 105 was performed against public protein databases, including the NCBI NR, SwissProt  and KEGG  databases and the potato (PGSC DM 3.4) and tomato (ITAG 2.3) protein sets. KEGG pathways were retrieved from the KEGG web server (http://www.genome.jp/kegg/) . The output of the KEGG analysis includes orthology assignments and pathways that are populated with the orthology assignments. Domain-based alignments were carried out against the NCBI COG database  (http://www.ncbi.nlm.nih.gov/COG/) with a cut-off E-value of ≤ 1e−5. The resulting NR BLASTX hits were processed with BLAST2GO software  to retrieve the associated gene ontology terms with E-values ≤ 10−5 describing biological processes, molecular functions and cellular components .
The CDSs of each putative unigene were extracted according to the BLASTX results (homologous approach), with a minimum 150-bp cutoff value and the priority order of SwissProt, Solanum (tomato and potato) protein datasets and NR database if conflicting results were obtained. ESTSCAN software  was also used to determine the direction of sequences that did not align to any of the databases, and CDSs shorter than 150 bp were removed. To avoid missing potential coding transcripts, the unigenes for which CDSs were not predicted by either homologous or ESTSCAN approaches were subjected to an in-house script, which, like most gene prediction programs, uses fifth-order hidden Markov chains to model coding regions . Again, the CDSs shorter than 150 bp were removed. The resultant CDSs extracted from the eggplant and turkey berry unigenes were translated into amino acid sequences with the standard codon table.
Identification of gene orthologous groups
The translated eggplant and turkey berry amino acid sequences were pooled into a protein database with sequences (>50 amino acids) from another 11 plant species: S. lycopersicum (Sol Genomics Network ITAG2.3), S. tuberosum (Sol Genomics Network PGSC DM 3.4), A. thaliana (TAIR release 10), C. papaya (http://www.life.illinois.edu/plantbio/People/Faculty/Ming), V. vinifera (http://www.genoscope.cns.fr/externe/GenomeBrowser/Vitis/), P. trichocarpa (JGI release v2.0 annotation v2.2), P. persica (Phytozome v9.0), C. sinensis (http://citrus.hzau.edu.cn/orange/download/), M. truncatula (Medicago Genome Sequence Consortium release Mt 3.0), Z. mays (Maize Genome Project 5b.60 B73) and O. sativa japonica (MSU Release 7.0).
Self-to-self BLASTP was conducted for all amino acid sequences with a cut-off E-value of 1e−5, and hits with identity < 30% and coverage < 30% were removed. Orthologous groups were constructed from the BLASTP results with OrthoMCL v2.0.9  using default settings.
Phylogenetic tree reconstruction
Single-copy gene families were retrieved from OrthoMCL as described above and used for the following phylogenetic tree reconstruction steps. The families containing any sequences shorter than 200 amino acids were removed, the amino acid sequences in each family were aligned using MUSCLE v3.8.31  with default parameters, and the corresponding CDS alignments were back-translated from the corresponding amino acid sequence alignments. The families were further filtered if the CDS alignment contained any taxon for which >50% of the data was missing. The remaining CDS alignments of each family were separated into three sets corresponding to each of the three codon positions. The four supermatrices (all codon positions and each codon position) were then separately assembled into supergenes using an in-house Perl script. The refined supergene data were then subjected to maximum likelihood phylogenetic analyses using PhyML3.1 . The HKY85 + gamma substitution model was selected, and bootstrap values were calculated using the aLRT model (parameters: −d nt –m HKY85 –b −4 –a e -c 4). TreeBeST (version 1.9.2, http://treesoft.sourceforge.net/) was used to root the trees if necessary.
Estimation of divergence time
Identification of plant resistance genes
Amino acid sequences for 112 reference resistance genes were downloaded from the Plant Resistance Genes database (PRGdb; http://prgdb.crg.eu/wiki/Main_Page) . BLASTP was used to identify and classify putative resistance genes in eggplant, turkey berry, tomato potato and Arabidopsis (parameters: −evalue 1e-5 -outfmt 6 -max_target_seqs 1). By parsing tabular outputs using in-house PERL scripts, results were filtered with a threshold cut-off of 40% identity and 50% coverage, and then homologous sequences were extracted and classified.
The sequences reported in this paper have been deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) and Transcriptome Shotgun Assembly (TSA). Raw paired-end reads are available through the NCBI SRA under accession numbers [SRA: SRR1104129] (eggplant) and [SRA: SRR1104128] (turkey berry). Transcripts are available through the NCBI TSA under accession number GBEF00000000 (eggplant) and GBEG00000000 (turkey berry).
The work presented here was supported by the National Natural Science Foundation of China (Grant No. NSFC31171954), the Natural Science Foundation of Jiangsu Province (Grant No. BK2010320), China Postdoctoral Science Foundation (Grant No. 20110491463) and Jiangsu Postdoctoral Science Foundation (Grant No. 1002001C).
- Xu X, Pan S, Cheng S, Zhang B, Mu D, Ni P, Zhang G, Yang S, Li R, Wang J: Genome sequence and analysis of the tuber crop potato. Nature. 2011, 475 (7355): 189-195. 10.1038/nature10158.PubMedView ArticleGoogle Scholar
- Zouine M, Latché A, Rousseau C, Regad F, Pech J-C, Philippot M, Bouzayen M, Delalande C, Frasse P, Schiex T: The tomato genome sequence provides insights into fleshy fruit evolution. Nature. 2012, 485: 635-641. 10.1038/nature11119.View ArticleGoogle Scholar
- Knapp S, Vorontsova MS, Prohens J: Wild relatives of the eggplant (Solanum melongena L.: Solanaceae): new understanding of species names in a complex group. PLoS One. 2013, 8 (2): e57039-10.1371/journal.pone.0057039.PubMed CentralPubMedView ArticleGoogle Scholar
- Daunay M, Lester R, Gebhardt C, Hennart J, Jahn M, Frary A, Doganlar S: Genetic Resources of Eggplant (Solanum melongena L.) and Allied Species: a New Challenge for Molecular Geneticists and Eggplant Breeders. 2001, Nijmegan, The Netherlands: Nijmegen University PressGoogle Scholar
- Arumuganathan K, Earle E: Nuclear DNA content of some important plant species. Plant Mol Biol Rep. 1991, 9 (3): 208-218. 10.1007/BF02672069.View ArticleGoogle Scholar
- Collonnier C, Fock I, Kashyap V, Rotino G, Daunay M, Lian Y, Mariska I, Rajam M, Servaes A, Ducreux G: Applications of biotechnology in eggplant. Plant Cell Tiss Org Cult. 2001, 65 (2): 91-107. 10.1023/A:1010674425536.View ArticleGoogle Scholar
- Jaiswal B: Solanum torvum: a review of its traditional uses, phytochemistry and pharmacology. International Journal of Pharma and Bio Sciences. 2012, 3 (4):
- Yamaguchi H, Fukuoka H, Arao T, Ohyama A, Nunome T, Miyatake K, Negoro S: Gene expression analysis in cadmium-stressed roots of a low cadmium-accumulating solanaceous plant, Solanum torvum. J Exp Bot. 2010, 61 (2): 423-437. 10.1093/jxb/erp313.PubMed CentralPubMedView ArticleGoogle Scholar
- Gousset C, Collonnier C, Mulya K, Mariska I, Rotino GL, Besse P, Servaes A, Sihachakr D: Solanum torvum, as a useful source of resistance against bacterial and fungal diseases for improvement of eggplant (S. melongena L.). Plant Sci. 2005, 168 (2): 319-327. 10.1016/j.plantsci.2004.07.034.View ArticleGoogle Scholar
- Garibaldi A, Minuto A, Gullino M: Verticillium wilt incited by Verticillium dahliae in eggplant grafted on Solanum torvum in Italy. Plant Dis. 2005, 89 (7): 777-Google Scholar
- Daunay M, Lester R, Laterrot H: The use of wild species for the genetic improvement of Brinjal eggplant (Solanum melongena) and tomato (Lycopersicon esculentum). Solanaceae III: Taxonomy, Chemistry, Evolution. 1991, 27: 389-413.Google Scholar
- Guri A, Sink K: Interspecific somatic hybrid plants between eggplant (Solanum melongena) and Solanum torvum. Theor Appl Genet. 1988, 76 (4): 490-496.PubMedView ArticleGoogle Scholar
- Van Eck J, Snyder A: Eggplant (Solanum melongena L.). Agrobacterium Protocols. 2006, New Jersey: Humana Press, 439-448.View ArticleGoogle Scholar
- Levin RA, Myers NR, Bohs L: Phylogenetic relationships among the “spiny solanums” (Solanum subgenus Leptostemonum, Solanaceae). Am J Bot. 2006, 93 (1): 157-169. 10.3732/ajb.93.1.157.View ArticleGoogle Scholar
- Fukuoka H, Yamaguchi H, Nunome T, Negoro S, Miyatake K, Ohyama A: Accumulation, functional annotation, and comparative analysis of expressed sequence tags in eggplant (Solanum melongena L.), the third pole of the genus Solanum species after tomato and potato. Gene. 2010, 450 (1): 76-84.PubMedView ArticleGoogle Scholar
- Bagnaresi P, Sala T, Irdani T, Scotto C, Lamontanara A, Beretta M, Rotino G, Sestili S, Cattivelli L, Sabatini E: Solanum torvum responses to the root-knot nematode Meloidogyne incognita. BMC Genomics. 2013, 14 (1): 540-10.1186/1471-2164-14-540.PubMed CentralPubMedView ArticleGoogle Scholar
- Barker MS, Vogel H, Schranz ME: Paleopolyploidy in the Brassicales: analyses of the Cleome transcriptome elucidate the history of genome duplications in Arabidopsis and other Brassicales. Genome Biol Evol. 2009, 1: 391-PubMed CentralPubMedView ArticleGoogle Scholar
- McKain MR, Wickett N, Zhang Y, Ayyampalayam S, McCombie WR, Chase MW, Pires JC, Leebens-Mack J: Phylogenomic analysis of transcriptome data elucidates co-occurrence of a paleopolyploid event and the origin of bimodal karyotypes in Agavoideae (Asparagaceae). Am J Bot. 2012, 99 (2): 397-406. 10.3732/ajb.1100537.PubMedView ArticleGoogle Scholar
- Wen J, Xiong Z, Nie Z-L, Mao L, Zhu Y, Kan X-Z, Ickert-Bond SM, Gerrath J, Zimmer EA, Fang X-D: Transcriptome sequences resolve deep relationships of the grape family. PLoS One. 2013, 8 (9): e74394-10.1371/journal.pone.0074394.PubMed CentralPubMedView ArticleGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011, 29 (7): 644-652. 10.1038/nbt.1883.PubMed CentralPubMedView ArticleGoogle Scholar
- Fischer S, Brunk BP, Chen F, Gao X, Harb OS, Iodice JB, Shanmugam D, Roos DS, Stoeckert CJ: Using OrthoMCL to Assign Proteins to OrthoMCL‒DB Groups or to Cluster Proteomes Into New Ortholog Groups. Curr Protoc Bioinformatics. 2011, 35: 6.12. 11-16.12. 19.Google Scholar
- Wang K, Wang Z, Li F, Ye W, Wang J, Song G, Yue Z, Cong L, Shang H, Zhu S: The draft genome of a diploid cotton Gossypium raimondii. Nat Genet. 2012, 44 (10): 1098-1103. 10.1038/ng.2371.PubMedView ArticleGoogle Scholar
- D’Hont A, Denoeud F, Aury J-M, Baurens F-C, Carreel F, Garsmeur O, Noel B, Bocs S, Droc G, Rouard M: The banana (Musa acuminata) genome and the evolution of monocotyledonous plants. Nature. 2012, 488 (7410): 213-217. 10.1038/nature11241.PubMedView ArticleGoogle Scholar
- Peng Z, Lu Y, Li L, Zhao Q, Feng Q, Gao Z, Lu H, Hu T, Yao N, Liu K: The draft genome of the fast-growing non-timber forest species moso bamboo (Phyllostachys heterocycla). Nat Genet. 2013, 45 (4): 456-461. 10.1038/ng.2569.PubMedView ArticleGoogle Scholar
- Guindon S, Dufayard J-F, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010, 59 (3): 307-321. 10.1093/sysbio/syq010.PubMedView ArticleGoogle Scholar
- Bremer B, Bremer K, Chase M, Fay M, Reveal J, Soltis D, Soltis P, Stevens P: An update of the Angiosperm Phylogeny Group classification for the orders and families of flowering plants: APG III. Botanical Journal of the Linnean Society. 2009, 161 (2): 105-121.View ArticleGoogle Scholar
- Jeffroy O, Brinkmann H, Delsuc F, Philippe H: Phylogenomics: the beginning of incongruence?. Trends Genet. 2006, 22 (4): 225-231. 10.1016/j.tig.2006.02.003.PubMedView ArticleGoogle Scholar
- Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24 (8): 1586-1591. 10.1093/molbev/msm088.PubMedView ArticleGoogle Scholar
- Rannala B, Yang Z: Inferring speciation times under an episodic molecular clock. Syst Biol. 2007, 56 (3): 453-466. 10.1080/10635150701420643.PubMedView ArticleGoogle Scholar
- Hedges SB, Dudley J, Kumar S: TimeTree: a public knowledge-base of divergence times among organisms. Bioinformatics. 2006, 22 (23): 2971-2972. 10.1093/bioinformatics/btl505.PubMedView ArticleGoogle Scholar
- Wang Y, Diehl A, Wu F, Vrebalov J, Giovannoni J, Siepel A, Tanksley SD: Sequencing and comparative analysis of a conserved syntenic segment in the Solanaceae. Genetics. 2008, 180 (1): 391-408. 10.1534/genetics.108.087981.PubMed CentralPubMedView ArticleGoogle Scholar
- Fei J, Chai Y, Wang J, Lin J, Sun X, Sun C, Zuo K, Tang K: cDNA cloning and characterization of the Ve homologue gene StVe from Solanum torvum Swartz. Mitochondrial DNA. 2004, 15 (2): 88-95.Google Scholar
- Sanseverino W, Hermoso A, D’Alessandro R, Vlasova A, Andolfo G, Frusciante L, Lowy E, Roma G, Ercolano MR: PRGdb 2.0: towards a community-based database model for the analysis of R-genes in plants. Nucleic Acids Res. 2013, 41 (D1): D1167-D1171. 10.1093/nar/gks1183.PubMed CentralPubMedView ArticleGoogle Scholar
- van Ooijen G, van den Burg HA, Cornelissen BJ, Takken FL: Structure and function of resistance proteins in solanaceous plants. Annu Rev Phytopathol. 2007, 45: 43-72. 10.1146/annurev.phyto.45.062806.094430.PubMedView ArticleGoogle Scholar
- Sanseverino W, Roma G, De Simone M, Faino L, Melito S, Stupka E, Frusciante L, Ercolano MR: PRGdb: a bioinformatics platform for plant resistance gene analysis. Nucleic Acids Res. 2010, 38 (suppl 1): D814-D821.PubMed CentralPubMedView ArticleGoogle Scholar
- Meyers BC, Kozik A, Griego A, Kuang H, Michelmore RW: Genome-wide analysis of NBS-LRR–encoding genes in Arabidopsis. Plant Cell Online. 2003, 15 (4): 809-834. 10.1105/tpc.009308.View ArticleGoogle Scholar
- Leister D: Tandem and segmental gene duplication and recombination in the evolution of plant disease resistance genes. Trends Genet. 2004, 20 (3): 116-122. 10.1016/j.tig.2004.01.007.PubMedView ArticleGoogle Scholar
- Wu J, Wang Z, Shi Z, Zhang S, Ming R, Zhu S, Khan MA, Tao S, Korban SS, Wang H: The genome of the pear (Pyrus bretschneideri Rehd.). Genome Res. 2013, 23 (2): 396-408. 10.1101/gr.144311.112.PubMed CentralPubMedView ArticleGoogle Scholar
- Han X-J, Wang Y-D, Chen Y-C, Lin L-Y, Wu Q-K: Transcriptome sequencing and expression analysis of terpenoid biosynthesis genes in Litsea cubeba. PLoS One. 2013, 8 (10): e76890-10.1371/journal.pone.0076890.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL: TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14 (4): R36-10.1186/gb-2013-14-4-r36.PubMed CentralPubMedView ArticleGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The sequence alignment/map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-2079. 10.1093/bioinformatics/btp352.PubMed CentralPubMedView ArticleGoogle Scholar
- Boeckmann B, Bairoch A, Apweiler R, Blatter M-C, Estreicher A, Gasteiger E, Martin MJ, Michoud K, O'Donovan C, Phan I: The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 2003, 31 (1): 365-370. 10.1093/nar/gkg095.PubMed CentralPubMedView ArticleGoogle Scholar
- Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M: KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012, 40 (D1): D109-D114. 10.1093/nar/gkr988.PubMed CentralPubMedView ArticleGoogle Scholar
- Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T: KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36 (suppl 1): D480-D484.PubMed CentralPubMedGoogle Scholar
- Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, Krylov DM, Mazumder R, Mekhedov SL, Nikolskaya AN: The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003, 4 (1): 41-10.1186/1471-2105-4-41.PubMed CentralPubMedView ArticleGoogle Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.PubMedView ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT: Gene Ontology: tool for the unification of biology. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.PubMed CentralPubMedView ArticleGoogle Scholar
- Iseli C, Jongeneel CV, Bucher P: ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. ISMB. 1999, 1999: 138-148.Google Scholar
- Birney E, Clamp M, Durbin R: GeneWise and genomewise. Genome Res. 2004, 14 (5): 988-995. 10.1101/gr.1865504.PubMed CentralPubMedView ArticleGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.PubMed CentralPubMedView ArticleGoogle Scholar
- Chaw S-M, Chang C-C, Chen H-L, Li W-H: Dating the monocot–dicot divergence and the origin of core eudicots using whole chloroplast genomes. J Mol Evol. 2004, 58 (4): 424-441. 10.1007/s00239-003-2564-9.PubMedView ArticleGoogle Scholar
- Crepet WL, Nixon KC, Gandolfo MA: Fossil evidence and phylogeny: the age of major angiosperm clades based on mesofossil and macrofossil evidence from Cretaceous deposits. Am J Bot. 2004, 91 (10): 1666-1682. 10.3732/ajb.91.10.1666.PubMedView ArticleGoogle Scholar
- de Sa MM, Drouin G: Phylogeny and substitution rates of angiosperm actin genes. Mol Biol Evol. 1996, 13 (9): 1198-1212. 10.1093/oxfordjournals.molbev.a025685.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.