Solanum aculeatissimum and Solanum torvum chloroplast genome sequences: a comparative analysis with other Solanum chloroplast genomes.

BACKGROUND: Solanum aculeatissimum and Solanum torvum belong to the Solanum species, and they are essential plants known for their high resistance to diseases and adverse conditions. They are frequently used as rootstocks for grafting and are often crossbred with other Solanum species to leverage their resistance traits. However, the phylogenetic relationship between S. aculeatissimum and S. torvum within the Solanum genus remains unclear. Therefore, this paper aims to sequence the complete chloroplast genomes of S. aculeatissimum and S. torvum and analyze them in comparison with 29 other previously published chloroplast genomes of Solanum species. RESULTS: We observed that the chloroplast genomes of S. aculeatissimum and S. torvum possess typical tetrameric structures, consisting of one Large Single Copy (LSC) region, two reverse-symmetric Inverted Repeats (IRs), and one Small Single Copy (SSC) region. The total length of these chloroplast genomes ranged from 154,942 to 156,004 bp, with minimal variation. The highest GC content was found in the IR region, while the lowest was in the SSC region. Regarding gene content, the total number of chloroplast genes and CDS genes remained relatively consistent, ranging from 128 to 134 and 83 to 91, respectively. Nevertheless, there was notable variability in the number of tRNA genes and rRNAs. Relative synonymous codon usage (RSCU) analysis revealed that both S. aculeatissimum and S. torvum preferred codons that utilized A and U bases. Analysis of the IR boundary regions indicated that contraction and expansion primarily occurred at the junction between SSC and IR regions. Nucleotide polymorphism analysis and structural variation analysis demonstrated that chloroplast variation in Solanum species mainly occurred in the LSC and SSC regions. Repeat sequence analysis revealed that A/T was the most frequent base pair in simple repeat sequences (SSR), while Palindromic and Forward repeats were more common in long sequence repeats (LSR), with Reverse and Complement repeats being less frequent. Phylogenetic analysis indicated that S. aculeatissimum and S. torvum belonged to the same meristem and were more closely related to Cultivated Eggplant. CONCLUSION: These findings enhance our comprehension of chloroplast genomes within the Solanum genus, offering valuable insights for plant classification, evolutionary studies, and potential molecular markers for species identification.

CAGEE: Computational Analysis of Gene Expression Evolution.

Despite the increasing abundance of whole transcriptome data, few methods are available to analyze global gene expression across phylogenies. Here, we present a new software package (Computational Analysis of Gene Expression Evolution [CAGEE]) for inferring patterns of increases and decreases in gene expression across a phylogenetic tree, as well as the rate at which these changes occur. In contrast to previous methods that treat each gene independently, CAGEE can calculate genome-wide rates of gene expression, along with ancestral states for each gene. The statistical approach developed here makes it possible to infer lineage-specific shifts in rates of evolution across the genome, in addition to possible differences in rates among multiple tissues sampled from the same species. We demonstrate the accuracy and robustness of our method on simulated data and apply it to a data set of ovule gene expression collected from multiple self-compatible and self-incompatible species in the genus Solanum to test hypotheses about the evolutionary forces acting during mating system shifts. These comparisons allow us to highlight the power of CAGEE, demonstrating its utility for use in any empirical system and for the analysis of most morphological traits. Our software is available at https://github.com/hahnlab/CAGEE/.

Abundance of the bat Sturnira lilium (Phyllostomidae) in relation Solanum mauritianum (Solanaceae) diaspores in an Atlantic Forest fragment of southern Brazil.

Frugivore bats are important seed dispersers in forests and their abundance are associated with the presence of zoochoric plants. In this context, the aim of our study was to investigate the association of the frugivore bat S. lilium with the diaspores of the zoochoric plant S. mauritianum, a common arboreal species present in forest fragments of southern Brazil. We also investigated the diet of the species based on seed content present in feces of individuals. Bats were mist-netted from November 2017 to April 2018 in a fragment of Atlantic Forest. The proportion of immature and mature diaspores of S. mauritianum was estimated in the same area where bats were sampled, and feces were sampled from captured individuals. In total, 61 individuals of S. lilium were captured, and 795 seeds were sampled from their feces. The abundance of S. lilium was significantly associated with the proportion of immature diaspores of S. mauritianum. We identified seeds of two botanical families: Solanaceae (89%) and Moraceae (11%) in the fecal samples. Our findings support the view that S. lilium is a legitimate disperser of S. mauritianum, and that its ecological function is probably a result of co-adaptation.

The evolution of inflorescence diversity in the nightshades and heterochrony during meristem maturation.

One of the most remarkable manifestations of plant evolution is the diversity for floral branching systems. These "inflorescences" arise from stem cell populations in shoot meristems that mature gradually to reproductive states in response to environmental and endogenous signals. The morphology of the shoot meristem maturation process is conserved across distantly related plants, raising the question of how diverse inflorescence architectures arise from seemingly common maturation programs. In tomato and related nightshades (Solanaceae), inflorescences range from solitary flowers to highly branched structures bearing hundreds of flowers. Since reproductive barriers between even closely related Solanaceae have precluded a genetic dissection, we captured and compared meristem maturation transcriptomes from five domesticated and wild species reflecting the evolutionary continuum of inflorescence complexity. We find these divergent species share hundreds of dynamically expressed genes, enriched for transcription factors. Meristem stages are defined by distinct molecular states and point to modified maturation schedules underlying architectural variation. These modified schedules are marked by a peak of transcriptome expression divergence during the reproductive transition, driven by heterochronic shifts of dynamic genes, including transcriptional regulators with known roles in flowering. Thus, evolutionary diversity in Solanaceae inflorescence complexity is determined by subtle modifications of transcriptional programs during a critical transitional window of meristem maturation, which we propose underlies similar cases of plant architectural variation. More broadly, our findings parallel the recently described transcriptome "inverse hourglass" model for animal embryogenesis, suggesting both plant and animal morphological variation is guided by a mid-development period of transcriptome divergence.

The Solanum commersonii Genome Sequence Provides Insights into Adaptation to Stress Conditions and Genome Evolution of Wild Potato Relatives.

Here, we report the draft genome sequence of Solanum commersonii, which consists of ∼830 megabases with an N50 of 44,303 bp anchored to 12 chromosomes, using the potato (Solanum tuberosum) genome sequence as a reference. Compared with potato, S. commersonii shows a striking reduction in heterozygosity (1.5% versus 53 to 59%), and differences in genome sizes were mainly due to variations in intergenic sequence length. Gene annotation by ab initio prediction supported by RNA-seq data produced a catalog of 1703 predicted microRNAs, 18,882 long noncoding RNAs of which 20% are shown to target cold-responsive genes, and 39,290 protein-coding genes with a significant repertoire of nonredundant nucleotide binding site-encoding genes and 126 cold-related genes that are lacking in S. tuberosum. Phylogenetic analyses indicate that domesticated potato and S. commersonii lineages diverged ∼2.3 million years ago. Three duplication periods corresponding to genome enrichment for particular gene families related to response to salt stress, water transport, growth, and defense response were discovered. The draft genome sequence of S. commersonii substantially increases our understanding of the domesticated germplasm, facilitating translation of acquired knowledge into advances in crop stability in light of global climate and environmental changes.

Greatly reduced phylogenetic structure in the cultivated potato clade (Solanum section Petota pro parte).

PREMISE OF THE STUDY: The species boundaries of wild and cultivated potatoes are controversial, with most of the taxonomic problems in the cultivated potato clade. We here provide the first in-depth phylogenetic study of the cultivated potato clade to explore possible causes of these problems. METHODS: We examined 131 diploid accessions, using 12 nuclear orthologs, producing an aligned data set of 14,072 DNA characters, 2171 of which are parsimony-informative. We analyzed the data to produce phylogenies and perform concordance analysis and goodness-of-fit tests. KEY RESULTS: There is good phylogenetic structure in clades traditionally referred to as clade 1+2 (North and Central American diploid potatoes exclusive of Solanum verrucosum), clade 3, and a newly discovered basal clade, but drastically reduced phylogenetic structure in clade 4, the cultivated potato clade. The results highlight a clade of species in South America not shown before, 'neocardenasii', sister to clade 1+2, that possesses key morphological traits typical of diploids in Mexico and Central America. Goodness-of-fit tests suggest potential hybridization between some species of the cultivated potato clade. However, we do not have enough phylogenetic signal with the data at hand to explicitly estimate such hybridization events with species networks methods. CONCLUSIONS: We document the close relationships of many of the species in the cultivated potato clade, provide insight into the cause of their taxonomic problems, and support the recent reduction of species in this clade. The discovery of the neocardenasii clade forces a reevaluation of a hypothesis that section Petota originated in Mexico and Central America.

Detecting Macroevolutionary Self-Destruction from Phylogenies.

Phylogenetic analyses have lent support to the concept of lineage selection: that biological lineages can have heritable traits that influence their capacity to persist and diversify, and thereby affect their representation in biodiversity. While many discussions have focused on "positive" lineage selection, where stably heritable properties of lineages enhance their diversification rate, there are also intriguing examples that seem to represent "negative" lineage selection, where traits reduce the likelihood that a lineage will persist or speciate. In this article, we test whether a particular pattern of negative lineage selection is detectable from the distributions of the trait on a phylogeny. "Self-destructive" traits are those that arise often but then disappear again because they confer either a raised extinction rate or they are prone to a high rate of trait loss. For such a trait, the reconstructed origins will tend to be dispersed across the tips of the phylogeny, rather than defining large clades of related lineages that all share the trait. We examine the utility of four possible measures of "tippiness" as potential indicators of macroevolutionary self-destruction, applying them to phylogenies on which trait evolution has been simulated under different combinations of parameters for speciation, extinction, trait gain, and trait loss. We use an efficient simulation approach that starts with the required number of tips with and without the trait and uses a model to work "backwards" to construct different possible trees that result in that set of tips. We then apply these methods to a number of case studies: salt tolerance in grasses, color polymorphism in birds of prey, and selfing in nightshades. We find that the relative age of species, measured from tip length, can indicate a reduced speciation rate but does not identify traits that increase the extinction rate or the trait loss rate. We show that it is possible to detect cases of macroevolutionary self-destruction by considering the number of tips with the trait that arise from each inferred origin, and the degree to which the trait is scattered across the phylogeny. These metrics, and the methods we present, may be useful for testing macroevolutionary hypotheses from phylogenetic patterns.

Dissimilar evolutionary histories of two resistance gene families in the genus Solanum.

Genomic analyses have shown that most genes in eukaryotic lineages belong to families. Gene families vary in terms of number of members, nucleotide similarity, gene integrity, expression, and function. Often, the members of gene families are arranged in clusters, which contribute to maintaining similarity among gene copies and also to generate duplicates through replication errors. Gene families offer us an opportunity to examine the forces involved in the evolution of the genomes and to study recombination events and genomic rearrangements. In this work, we focused on the evolution of two plant resistance gene families, Sw5 and Mi-1, and analyzed the completely sequenced nuclear genomes of potato and tomato. We first noticed that the potato genome carries larger resistance gene families than tomato, but all gene copies are pseudogenes. Second, phylogenetic analyses indicated that Sw5 and Mi-1 gene families had dissimilar evolutionary histories. In contrast to Sw5, Mi-1 homologues suffered repeated gene conversion events among the gene copies, particularly in the tomato genome.

Addition of straw from hyperaccumulator plants to cadmium-contaminated soil increases cadmium uptake by loquat seedlings.

The straw from three different cadmium (Cd) hyperaccumulators (Galinsoga parviflora, Youngia erythrocarpa, and Solanum photeinocarpum) was added to Cd-contaminated soil, and its effects on plant growth and Cd accumulation in loquat seedlings were evaluated. Straw from each of G. parviflora, Y. erythrocarpa, and S. photeinocarpum was added to Cd-contaminated soil before planting seedlings of two varieties of loquat (Dawuxing and Chuanzao). Addition of straw from G. parviflora and S. photeinocarpum increased the root and shoot biomasses of both loquat varieties, compared with that in the control. Addition of straw also increased Cd uptake by loquat seedlings. The treatments could be ranked, from highest Cd contents in roots and shoots of loquat seedlings to lowest, as follows: S. photeinocarpum straw > Y. erythrocarpa straw > G. parviflora straw > control. All three types of hyperaccumulator straw increased the amount of Cd extracted by shoots of two loquat seedlings, with the maximum effect in the S. photeinocarpum straw treatment. Addition of hyperaccumulator straw also resulted in increased soil invertase, urease, and catalase activities to varying degrees. Among the three types of hyperaccumulator straw, S. photeinocarpum straw was the most effective to increase Cd accumulation in loquat seedlings. Therefore, this material has the potential to increase the phytoremediation capacity of loquat seedlings in Cd-contaminated orchards.

Transcriptome analysis and molecular marker discovery in Solanum incanum and S. aethiopicum, two close relatives of the common eggplant (Solanum melongena) with interest for breeding.

BACKGROUND: Solanum incanum is a close wild relative of S. melongena with high contents of bioactive phenolics and drought tolerance. S. aethiopicum is a cultivated African eggplant cross-compatible with S. melongena. Despite their great interest in S. melongena breeding programs, the genomic resources for these species are scarce. RESULTS: RNA-Seq was performed with NGS from pooled RNA of young leaf, floral bud and young fruit tissues, generating more than one hundred millions raw reads per species. The transcriptomes were assembled in 83,905 unigenes for S. incanum and in 87,084 unigenes for S. aethiopicum with an average length of 696 and 722 bp, respectively. The unigenes were structurally and functionally annotated based on comparison with public databases by using bioinformatic tools. The single nucleotide variant calling analysis (SNPs and INDELs) was performed by mapping our S. incanum and S. aethiopicum reads, as well as reads from S. melongena and S. torvum available on NCBI database (National Center for Biotechnology Information), against the eggplant genome. Both intraspecific and interspecific polymorphisms were identified and subsets of molecular markers were created for all species combinations. 36 SNVs were selected for validation in the S. incanum and S. aethiopicum accessions and 96 % were correctly amplified confirming the polymorphisms. In addition, 976 and 1,278 SSRs were identified in S. incanum and S. aethiopicum transcriptomes respectively, and a set of them were validated. CONCLUSIONS: This work provides a broad insight into gene sequences and allelic variation in S. incanum and S. aethiopicum. This work is a first step toward better understanding of target genes involved in metabolic pathways relevant for eggplant breeding. The molecular markers detected in this study could be used across all the eggplant genepool, which is of interest for breeding programs as well as to perform marker-trait association and QTL analysis studies.

The first de novo transcriptome of pepino (Solanum muricatum): assembly, comprehensive analysis and comparison with the closely related species S. caripense, potato and tomato.

BACKGROUND: Solanum sect. Basarthrum is phylogenetically very close to potatoes (Solanum sect. Petota) and tomatoes (Solanum sect. Lycopersicon), two groups with great economic importance, and for which Solanum sect. Basarthrum represents a tertiary gene pool for breeding. This section includes the important regional cultigen, the pepino (Solanum muricatum), and several wild species. Among the wild species, S. caripense is prominent due to its major involvement in the origin of pepino and its wide geographical distribution. Despite the value of the pepino as an emerging crop, and the potential for gene transfer from both the pepino and S. caripense to potatoes and tomatoes, there has been virtually no genomic study of these species. RESULTS: Using Illumina HiSeq 2000, RNA-Seq was performed with a pool of three tissues (young leaf, flowers in pre-anthesis and mature fruits) from S. muricatum and S. caripense, generating almost 111,000,000 reads among the two species. A high quality de novo transcriptome was assembled from S. muricatum clean reads resulting in 75,832 unigenes with an average length of 704 bp. These unigenes were functionally annotated based on similarity of public databases. We used Blast2GO, to conduct an exhaustive study of the gene ontology, including GO terms, EC numbers and KEGG pathways. Pepino unigenes were compared to both potato and tomato genomes in order to determine their estimated relative position, and to infer gene prediction models. Candidate genes related to traits of interest in other Solanaceae were evaluated by presence or absence and compared with S. caripense transcripts. In addition, by studying five genes, the phylogeny of pepino and five other members of the family, Solanaceae, were studied. The comparison of S. caripense reads against S. muricatum assembled transcripts resulted in thousands of intra- and interspecific nucleotide-level variants. In addition, more than 1000 SSRs were identified in the pepino transcriptome. CONCLUSIONS: This study represents the first genomic resource for the pepino. We suggest that the data will be useful not only for improvement of the pepino, but also for potato and tomato breeding and gene transfer. The high quality of the transcriptome presented here also facilitates comparative studies in the genus Solanum. The accurate transcript annotation will enable us to figure out the gene function of particular traits of interest. The high number of markers (SSR and nucleotide-level variants) obtained will be useful for breeding programs, as well as studies of synteny, diversity evolution, and phylogeny.

Evolution of a complex locus for terpene biosynthesis in solanum.

Functional gene clusters, containing two or more genes encoding different enzymes for the same pathway, are sometimes observed in plant genomes, most often when the genes specify the synthesis of specialized defensive metabolites. Here, we show that a cluster of genes in tomato (Solanum lycopersicum; Solanaceae) contains genes for terpene synthases (TPSs) that specify the synthesis of monoterpenes and diterpenes from cis-prenyl diphosphates, substrates that are synthesized by enzymes encoded by cis-prenyl transferase (CPT) genes also located within the same cluster. The monoterpene synthase genes in the cluster likely evolved from a diterpene synthase gene in the cluster by duplication and divergence. In the orthologous cluster in Solanum habrochaites, a new sesquiterpene synthase gene was created by a duplication event of a monoterpene synthase followed by a localized gene conversion event directed by a diterpene synthase gene. The TPS genes in the Solanum cluster encoding cis-prenyl diphosphate-utilizing enzymes are closely related to a tobacco (Nicotiana tabacum; Solanaceae) diterpene synthase encoding Z-abienol synthase (Nt-ABS). Nt-ABS uses the substrate copal-8-ol diphosphate, which is made from the all-trans geranylgeranyl diphosphate by copal-8-ol diphosphate synthase (Nt-CPS2). The Solanum gene cluster also contains an ortholog of Nt-CPS2, but it appears to encode a nonfunctional protein. Thus, the Solanum functional gene cluster evolved by duplication and divergence of TPS genes, together with alterations in substrate specificity to utilize cis-prenyl diphosphates and through the acquisition of CPT genes.

The effects of abscisic acid (ABA) addition on cadmium accumulation of two ecotypes of Solanum photeinocarpum.

The study of the effects of exogenous abscisic acid (ABA) addition on cadmium (Cd) accumulation of two ecotypes (mining and farmland) of Solanum photeinocarpum was operated through a pot experiment. The results showed that the biomass and chlorophyll content of the two ecotypes of S. photeinocarpum increased with increasing ABA concentration. Applying exogenous ABA increased Cd content in the two ecotypes of S. photeinocarpum. The maximum Cd contents in shoots of the two ecotypes of S. photeinocarpum were obtained at 20 µmol/L ABA; shoot Cd contents respectively for the mining and farmland ecotypes were 33.92 and 24.71% higher than those for the control. Applying exogenous ABA also increased Cd extraction by the two ecotypes of S. photeinocarpum, and the highest Cd extraction was obtained at 20 µmol/L ABA with 569.42 µg/plant in shoots of the mining ecotype and 520.51 µg/plant in shoots of the farmland ecotype respectively. Therefore, exogenous ABA can be used for enhancing the Cd extraction ability of S. photeinocarpum, and 20 µmol/L ABA was the optimal dose.

Salinity-induced expression of HKT may be crucial for Na(+) exclusion in the leaf blade of huckleberry (Solanum scabrum Mill.), but not of eggplant (Solanum melongena L.).

Reduced Na(+) accumulation in the leaf blade is an important aspect of salinity tolerance and high affinity K(+) transporters (HKTs) are known to play a significant role in the process. Huckleberry and eggplant have previously been shown to display 'excluder' and 'includer' characteristics, respectively, under salt stress, but the underlying mechanisms have not been investigated. Here, we isolated the cDNA of the HKT homologs, Solanum scabrum HKT (SsHKT) from huckleberry and Solanum melongena HKT (SmHKT) from eggplant, and analyzed their expressions in different tissues under salt stress. SsHKT expression was markedly induced in the root (28-fold) and stem (7-fold), with a corresponding increase in Na(+) accumulation of 52% and 29%, respectively. Conversely, eggplant accumulated 60% total Na(+) in the leaf blade, with a lower SmHKT expression level in the root (3-fold). Huckleberry also maintained a higher K(+)/Na(+) ratio in the leaf blade compared to eggplant, due to the reduction of its Na(+) concentration and unaltered K(+) concentration. Functional analysis demonstrated that SsHKT-mediated Na(+) influx inhibited yeast growth under Na(+) stress, and that SsHKT did not complement the growth of the K(+) uptake-deficient CY162 strain under K(+)-limiting conditions. These results suggest that the Na(+) accumulation characteristics of both plants are caused by the differential expression of HKT genes, with SsHKT exerting a greater control over the ability of Na(+) to reach the leaf blade in huckleberry, than SmHKT does in eggplant.

The development of type VI glandular trichomes in the cultivated tomato Solanum lycopersicum and a related wild species S. habrochaites.

BACKGROUND: Type VI glandular trichomes represent the most abundant trichome type on leaves and stems of tomato plants and significantly contribute to herbivore resistance, particularly in the wild species. Despite this, their development has been poorly studied so far. The goal of this study is to fill this gap. Using a variety of cell imaging techniques, a detailed record of the anatomy and developmental stages of type VI trichomes in the cultivated tomato (Solanum lycopersicum) and in a related wild species (S. habrochaites) is provided. RESULTS: In both species, the development of these structures follows a highly reproducible cell division pattern. The two species differ in the shape of the trichome head which is round in S. habrochaites and like a four-leaf clover in S. lycopersicum, correlating with the presence of a large intercellular cavity in S. habrochaites where the produced metabolites accumulate. In both species, the junction between the intermediate cell and the four glandular cells constitute a breaking point facilitating the decapitation of the trichome and thereby the quick release of the metabolites. A strongly auto-fluorescent compound transiently accumulates in the early stages of development suggesting a potential role in the differentiation process. Finally, immuno-labelling with antibodies recognizing specific cell wall components indicate a key role of pectin and arabinogalactan components in the differentiation of type VI trichomes. CONCLUSIONS: Our observations explain the adaptive morphologies of type VI trichomes for metabolite storage and release and provide a framework for further studies of these important metabolic cellular factories. This is required to better exploit their potential, in particular for the breeding of pest resistance in tomato.

Diversification and distinctive structural features of S-RNase alleles in the genus Solanum.

The multigenic and multiallelic S-locus in plants is responsible for the gametophytic self-incompatibility system, which is important to prevent the detrimental effects of self-fertilization and inbreeding depression. Several studies have discussed the importance of punctual mutations, recombination, and natural selection in the generation of allelic diversity in the S-locus. However, there has been no wide-ranging study correlating the molecular evolution and structural aspects of the corresponding proteins in Solanum. Therefore, we evaluated the molecular evolution of one gene in this locus and generated a statistically well-supported phylogenetic tree, as well as evidence of positive selection, helping us to understand the diversification of S alleles in Solanum. The three-dimensional structures of some of the proteins corresponding to the major clusters of the phylogenetic tree were constructed and subsequently submitted to molecular dynamics to stabilize the folding and obtain the native structure. The positively selected amino acid residues were predominantly located in the hyper variable regions and on the surface of the protein, which appears to be fundamental for allele specificity. One of the positively selected residues was identified adjacent to a conserved strand that is crucial for enzymatic catalysis. Additionally, we have shown significant differences in the electrostatic potential among the predicted molecular surfaces in S-RNases. The structural results indicate that local changes in the three-dimensional structure are present in some regions of the molecule, although the general structure seems to be conserved. No previous study has described such structural variations in S-RNases.

Testing Taxonomic Predictivity of Foliar and Tuber Resistance to Phytophthora infestans in Wild Relatives of Potato.

Potato late blight, caused by the oomycete phytopathogen Phytophthora infestans, is a devastating disease found in potato-growing regions worldwide. Long-term management strategies to control late blight include the incorporation of host resistance to predominant strains. However, due to rapid genetic changes within pathogen populations, rapid and recurring identification and integration of novel host resistance traits is necessary. Wild relatives of potato offer a rich source of desirable traits, including late blight resistance, but screening methods can be time intensive. We tested the ability of taxonomy, ploidy, crossing group, breeding system, and geography to predict the presence of foliar and tuber late blight resistance in wild Solanum spp. Significant variation for resistance to both tuber and foliar late blight was found within and among species but there was no discernable predictive power based on taxonomic series, clade, ploidy, breeding system, elevation, or geographic location. We observed a moderate but significant correlation between tuber and foliar resistance within species. Although previously uncharacterized sources of both foliar and tuber resistance were identified, our study does not support an assumption that taxonomic or geographic data can be used to predict sources of late blight resistance in wild Solanum spp.

The WRKY transcription factor genes in eggplant (Solanum melongena L.) and Turkey Berry (Solanum torvum Sw.).

WRKY transcription factors, which play critical roles in stress responses, have not been characterized in eggplant or its wild relative, turkey berry. The recent availability of RNA-sequencing data provides the opportunity to examine WRKY genes from a global perspective. We identified 50 and 62 WRKY genes in eggplant (SmelWRKYs) and turkey berry (StorWRKYs), respectively, all of which could be classified into three groups (I-III) based on the WRKY protein structure. The SmelWRKYs and StorWRKYs contain ~76% and ~95% of the number of WRKYs found in other sequenced asterid species, respectively. Positive selection analysis revealed that different selection constraints could have affected the evolution of these groups. Positively-selected sites were found in Groups IIc and III. Branch-specific selection pressure analysis indicated that most WRKY domains from SmelWRKYs and StorWRKYs are conserved and have evolved at low rates since their divergence. Comparison to homologous WRKY genes in Arabidopsis revealed several potential pathogen resistance-related SmelWRKYs and StorWRKYs, providing possible candidate genetic resources for improving stress tolerance in eggplant and probably other Solanaceae plants. To our knowledge, this is the first report of a genome-wide analyses of the SmelWRKYs and StorWRKYs.

Comparative transcriptome analysis of eggplant (Solanum melongena L.) and turkey berry (Solanum torvum Sw.): phylogenomics and disease resistance analysis.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

Genetic structure of the four wild tomato species in the Solanum peruvianum s.l. species complex.

The most diverse wild tomato species Solanum peruvianum sensu lato (s.l.) has been reclassified into four separate species: Solanum peruvianum sensu stricto (s.s.), Solanum corneliomuelleri, Solanum huaylasense, and Solanum arcanum. However, reproductive barriers among the species are incomplete and this can lead to discrepancies regarding genetic identity of germplasm. We used genotyping by sequencing (GBS) of S. peruvianum s.l., Solanum neorickii, and Solanum chmielewskii to develop tens of thousands of mapped single nucleotide polymorphisms (SNPs) to analyze genetic relationships within and among species. The data set was condensed to 14,043 SNPs with no missing data across 46 sampled plants. Origins of accessions were mapped using geographical information systems (GIS). Isolation by distance, pairwise genetic distances, and number of clusters were estimated using population genetics approaches. Isolation by distance was strongly supported, especially between interspecific pairs. Eriopersicon (S. peruvianum s.s., S. corneliomuelleri, S. huaylasense) and Arcanum (S. arcanum, S. neorickii, S. chmielewskii) species groups were genetically distinct, except for S. huaylasense which showed 50% membership proportions in each group. Solanum peruvianum and S. corneliomuelleri were not significantly differentiated from each other. Many thousands of SNP markers were identified that could potentially be used to distinguish pairs of species, including S. peruvianum versus S. corneliomuelleri, if they are verified on larger numbers of samples. Diagnostic markers will be valuable for delimiting morphologically similar and interfertile species in germplasm management. Approximately 12% of the SNPs rejected a genome-wide test of selective neutrality based on differentiation among species of S. peruvianum s.l. These are candidates for more comprehensive studies of microevolutionary processes within this species complex.