Comparative transcriptome analysis of eggplant (Solanum melongena L.) and turkey berry (Solanum torvum Sw.): phylogenomics and disease resistance analysis
BMC Genomics volume 15, Article number: 412 (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.
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
To maximize the range of transcript diversity and completeness, mixed RNA samples from three tissues of each plant were prepared for Illumina sequencing. We obtained 2.24 Gb and 3.94 Gb of sequence from eggplant and turkey berry respectively (Table 1), and the raw paired-end data were deposited in the NCBI Sequence Read Archive. The cleaned reads were aligned to the genomes of the closely related Solanum species tomato and potato to assess sequencing completeness. As shown in Figure 1 (rings A1–A3), the depth distribution of eggplant and turkey berry fit well to the tomato gene distribution. Similarly, the eggplant and turkey berry reads fit well with the potato gene distribution (Additional file 1: Figure S1, rings A1–A3). These results indicate that the sequencing reads obtained from eggplant and turkey berry covered the majority of genes in these species.
Clean reads from the two Solanum species were then separately assembled into contigs and clustered into transcripts using the de novo transcriptome assembler Trinity, which can efficiently reconstruct full-length transcripts across a broad range of expression levels and sequencing depths . The clustering step substantially improved the assembly quality, as indicated by elevated N50 values and decreased total length, by eliminating redundant contigs (Table 1 and Figure 2A). Similar transcripts in the same cluster are thought to be isoforms (splice variants) at the gene locus . To further eliminate redundant transcripts and to obtain the primary representative of each gene locus, only the longest transcript in each cluster was regarded as the final assembled unigene. This process identified 34,174 unigenes for eggplant and 38,185 unigenes for turkey berry (Table 1), which included 9,743 (28.51%) and 10,762 (28.18%) unigenes longer than 1 kb respectively. We observed a decrease in N50 values of unigenes compared with transcripts, suggesting that longer genes may tend to generate more isoforms. This hypothesis was confirmed by plotting unigene length against the average number of isoforms in each bin and performing a Pearson's correlation coefficient test (Figure 2B), which showed a significant positive correlation for both eggplant and turkey berry.
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 provides important information on gene function and structure. We were able to annotate 81.98% (28,016) of the eggplant unigenes and 78.16% (29,845) of the turkey berry unigenes with a threshold of 1e–5 by performing a BLASTX search against diverse protein databases. When we extracted and aligned the putative CDSs, 86.96% (29,717) of eggplant unigenes and 84.03% (32,086) of turkey berry were annotated (Table 2). These results further confirmed the high quality of the de novo assembly.
In a BLASTX homolog search against the NCBI non-redundant (NR) protein database, 27,393 eggplant unigenes and 29,072 turkey berry unigenes had matches (Table 2), 78.0% and 75.4% respectively, of which showed >80% identity (Figure 3A), indicating the high accuracy of the assembly. For both species, the top hit species was tomato, followed by potato and then grape (Vitis vinifera) (Figure 3B). Interestingly, only 2.1% of the top hits were assigned to potato, which is much less than the 86.6% of eggplant and 84.3% of turkey berry hits that were assigned to tomato. A similar result was observed in an EST-based comparative analysis of eggplant , suggesting that these two species are more closely related to tomato than potato.
Comparative analysis of gene sets between plants
A total of 427,731 proteins from eggplant (29,717), turkey berry (32,086) and 11 other plant species, including tomato, potato, Arabidopsis thaliana, Carica papaya, V. vinifera, Prunus persica, Citrus sinensis, Medicago truncatula, Zea mays and Oryza sativa japonica, were binned into 36,627 orthologous groups (gene families) using OrthoMCL v2.0.9  following self-self-comparison with the BLASTP program. The average number of genes in each gene family (Table 3), the number of unique gene families (Figure 4A), and number of genes in the unique gene families (Figure 4B) of eggplant and turkey berry were less than those of tomato, potato and other plants. This suggests that either eggplant and turkey berry have distinct gene family features or that our gene sets are incomplete. Although our RNA libraries were derived from mixed tissue samples, it is likely that not all genes in the genome are represented in our transcriptomes.
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.
The CDS alignments from the 276 refined single-copy orthologous groups were first concatenated to form one supergene for each species, each of which was then subjected to phylogenetic analyses with the maximum likelihood method in PhyML3.1 . Unexpectedly, the phylogenies obtained (Additional file 1: Figure S4A) were incongruent with the well-recognized Angiosperm Phylogeny Group III (APG III) system . Notably, the branch lengths (indicating substitutions per site) varied considerably in our tree, indicating relatively variable evolution rates among species. Quite different substitution rates are commonly observed for the three positions within codons, with the third position being especially variable as a result of the degeneracy of the genetic code. Third-position substitutions are likely to be saturated and may accumulate mutational bias, which may influence the accuracy of phylogeny estimations . Therefore, the CDS alignments of each of the 276 gene families were separated into three datasets corresponding to each of the three codon positions in the CDS, and another three supergenes were assembled and used to estimate phylogeny. As predicted, the three maximum likelihood trees were identical (Figure 5 and Additional file 1: Figure S4B–D) and placed the monocot, Asterids, grape and Eurosids clades in accordance with the APG III system. Notably, all the clades leading to Asterid species had 100% bootstrap support values, even in the uncorrected tree (Additional file 1: Figure S4), implying that the RNA-Seq assemblies may not be responsible for the incongruence of phylogenies that we observed when using full-length CDS sequences and also providing robust support for the known relationships in Asterid species. As shown in Figure 5, eggplant was most closely related to its Leptostemonum clade relative turkey berry, and further separated from the members of the Potatoe clade, tomato and potato [14, 15].
Estimation of divergence time
The three codon position–based supergene sets from the 276 single-copy orthologous groups were used for combination analysis of multi-partitions in the MCMCTree program (PAML4.7 package) . The same substitution model was used, but different parameters were assigned and estimated for each set. Moreover, because of the variable evolution rate among species we observed, the clock model with independent rates among lineages specified by a log-normal probability distribution was adopted . To check the robustness of results, we ran the MCMCTree analysis twice and obtained similar results, and a chronogram (Figure 6) was produced using FigTree v1.4.0 (http://tree.bio.ed.ac.uk/) from the first run. Another dataset containing only the first two supergene sets (after removing the fast-evolving third position) was subjected to MCMCTree analysis, and a similar chronogram was obtained (Additional file 1: Figure S5).
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.
Amino acid sequences for the 112 reference resistance genes were downloaded from the PRGdb  and used to identify and classify putative resistance genes in Arabidopsis, eggplant, turkey berry, tomato and potato (Table 4), and the resistance gene distributions were plotted (Figure 1 and Additional file 1: Figure S1). This conservative approach revealed 336 resistance genes in Arabidopsis, including 44 CNL and 100 TNL class genes, which is comparable to results from domain prediction–based methods  in which 48 CNL and 89 TNL class genes were identified.
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
Two datasets were generated from the CDS alignments used for divergence time estimation: (1) a dataset containing the first two partitions, the first and second codon positions of the sequences; and (2) a set containing all the three partitions corresponding to all the three codon positions in the sequences. Divergence times were estimated under a relaxed clock model in the MCMCTree program in the PAML4.7 package , with “Independent rates model (clock = 2)” and “JC69 model” selected for our calculations. The MCMC process performs 40,000 iterations after a burn-in of 15,000 iterations. Other parameters were the default settings of MCMCTree. We ran the program twice for each dataset to confirm that the results were similar between runs. The following constraints were used for time calibrations:
140–150 Mya, monocot–dicot split 
94 Mya, lower boundary for Vitis–Eurosid split 
68–76 Mya, Caricaceae–Brassicaceae split 
44 Mya, upper boundary for the Solaneae 
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).
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.
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.
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.
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 Press
Arumuganathan K, Earle E: Nuclear DNA content of some important plant species. Plant Mol Biol Rep. 1991, 9 (3): 208-218. 10.1007/BF02672069.
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.
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.
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.
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-
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.
Guri A, Sink K: Interspecific somatic hybrid plants between eggplant (Solanum melongena) and Solanum torvum. Theor Appl Genet. 1988, 76 (4): 490-496.
Van Eck J, Snyder A: Eggplant (Solanum melongena L.). Agrobacterium Protocols. 2006, New Jersey: Humana Press, 439-448.
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.
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.
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.
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-
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24 (8): 1586-1591. 10.1093/molbev/msm088.
Rannala B, Yang Z: Inferring speciation times under an episodic molecular clock. Syst Biol. 2007, 56 (3): 453-466. 10.1080/10635150701420643.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Birney E, Clamp M, Durbin R: GeneWise and genomewise. Genome Res. 2004, 14 (5): 988-995. 10.1101/gr.1865504.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.
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.
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.
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.
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).
The authors have declared that no competing interests exist.
Conceived and designed the experiments: XY. Performed the experiments: YFC, YM. Analyzed the data: XY, YFC, CD, ZWW, YM. Contributed reagents/materials/analysis tools: XY, CD, XHC, LBX. Deposited the sequences in databases: YFC, CD. Wrote the manuscript: XY. All authors read and approved the final manuscript.
Xu Yang, Yu-Fu Cheng contributed equally to this work.
Electronic supplementary material
Additional file 1: Figure S1: Distributions of genomic elements of potato eggplant and turkey berry on potato genome. Figure S2: Distributions of depth of reads and densities of genes on tomato genome. Figure S3: Distributions of depth of reads and densities of genes on potato genome. Figure S4: Maximum likelihood trees based on 276 single-copy genes. Figure S5: Estimation of divergence time using the first and second codon positions. (DOCX 2 MB)
Authors’ original submitted files for images
About this article
Cite this article
Yang, X., Cheng, YF., Deng, C. et al. Comparative transcriptome analysis of eggplant (Solanum melongena L.) and turkey berry (Solanum torvum Sw.): phylogenomics and disease resistance analysis. BMC Genomics 15, 412 (2014). https://doi.org/10.1186/1471-2164-15-412