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

PE (Prickly Eggplant) encoding a cytokinin-activating enzyme responsible for the formation of prickles in eggplant

Eggplant is one of the most important vegetables worldwide, with some varieties displaying prickles. These prickles, present on the leaves, stems, and fruit calyxes, posing challenges during cultivation, harvesting, and transportation, making them an undesirable agronomic trait. However, the genetic mechanisms underlying prickle morphogenesis in eggplant remain poorly understood, impeding genetic improvements. In this study, genetic analyses revealed that prickle morphogenesis is governed by a single dominant nuclear gene, termed PE (Prickly Eggplant). Subsequent bulk segregant RNA-sequencing (BSR-seq) and linkage analysis preliminarily mapped PE to chromosome 6. This locus was then fine mapped to a 9233 bp interval in a segregating population of 1109 plants, harboring only one candidate gene, SmLOG1, which encodes a LONELY GUY (LOG)-family cytokinin biosynthetic enzyme. Expression analyses via transcriptome and qRT-PCR demonstrate that SmLOG1 is predominantly expressed in immature prickles. CRISPR-Cas9 knockout experiments targeting SmLOG1 in prickly parental line 'PI 381159' abolished prickles across all tissues, confirming its critical role in prickle morphogenesis. Sequence analysis of SmLOG1 pinpointed variations solely within the non-coding region. We developed a cleaved amplified polymorphic sequences (CAPS) marker from a distinct SNP located at -735-bp within the SmLOG1 promoter, finding significant association with prickle variation in 190 eggplant germplasms. These findings enhance our understanding of the molecular mechanisms governing prickle development in eggplant and facilitate the use of marker-assisted selection (MAS) for breeding prickleless cultivars.

Carboxypeptidase inhibitors from Solanaceae as a new subclass of pathogenesis related peptide aiming biotechnological targets for plant defense

Background: Plant protease inhibitors play a crucial role in inhibiting proteases produced by phytopathogens and exhibiting inhibitory effects on nematodes, fungi, and insects, making them promising candidates for crop protection. Specifically, carboxypeptidase inhibitors, a subset of proteinase inhibitors, have been extensively studied in potato and tomato of Solanaceae plant family. However, further research is needed to fully understand the functions and biotechnological potential of those inhibitors in plants. This work aimed to in silico characterize carboxypeptidase inhibitors from Solanaceae as potential antimicrobial and defense agents focused on biotechnological targets. Methods: The methodology employed involved search in UniProt, PDB, KNOTTIN, NCBI, and MEROPS databases for solanaceous carboxypeptidase inhibitors, phylogenetic relationships and conservation patterns analyzes using MEGA-X software and Clustal Omega/MView tools, physicochemical properties and antimicrobial potential prediction using ProtParam, ToxinPred, iAMPred, and APD3 tools, and structural features prediction using PSIPRED. Results and discussion: A systematic literature search was conducted to identify relevant studies on Solanaceae carboxypeptidase inhibitors and their activities against pathogens. The selected studies were reviewed and the main findings compiled. The characterization of Solanaceae carboxypeptidase inhibitors proposed for the first time the global sequence consensus motif CXXXCXXXXDCXXXXXCXXC, shedding light on carboxypeptidase inhibitors distribution, sequence variability, and conservation patterns. Phylogenetic analysis showed evolutionary relationships within the Solanaceae family, particularly in Capsicum, Nicotiana, and Solanum genera. Physicochemical characteristics of those peptides indicated their similarity to antimicrobial peptides. Predicted secondary structures exhibited variations, suggesting a broad spectrum of action, and studies had been demonstrated their activities against various pathogens. Conclusion: Carboxypeptidase inhibitors are being proposed here as a new subclass of PR-6 pathogenesis-related proteins, which will aid in a focused understanding of their functional roles in plant defense mechanisms. These findings confirm the Solanaceae carboxypeptidase inhibitors potential as defense agents and highlight opportunities for their biotechnological applications in pathogen control.

Comparative Proteomic Analysis of Roots from a Wild Eggplant Species Solanum sisymbriifolium in Defense Response to Verticillium dahliae Inoculation

Eggplant verticillium wilt, caused by Verticillium spp., is a severe eggplant vascular disease. Solanum sisymbriifolium, a wild species of eggplant that is resistant to verticillium wilt, will be beneficial for genetically modifying eggplants. To better reveal the response of wild eggplant to verticillium wilt, proteomic analysis by iTRAQ technique was performed on roots of S. sisymbriifolium after exposure to Verticillium dahliae, and some selected proteins were also validated using parallel reaction monitoring (PRM). After inoculation with V. dahliae, the phenylalanine ammonia lyase (PAL) and superoxide dismutase (SOD) enzymes and the malondialdehyde (MDA) and soluble protein (SP) of S. sisymbriifolium roots all exhibited an increase in activity or content compared with that of the mock-inoculated plants, especially at 12 and 24 h post-inoculation (hpi). A total of 4890 proteins (47.04% of the proteins were from S. tuberosum and 25.56% were from S. lycopersicum according to the species annotation) were identified through iTRAQ and LC-MS/MS analysis. A total of 369 differentially expressed proteins (DEPs) (195 downregulated and 174 upregulated) were obtained by comparison of the control and treatment groups at 12 hpi, and 550 DEPs (466 downregulated and 84 upregulated) were obtained by comparison of the groups at 24 hpi. The most significant Gene Ontology (GO) enrichment terms at 12 hpi were regulation of translational initiation, oxidation-reduction, and single-organism metabolic process in the biological process group; cytoplasm and eukaryotic preinitiation complex in the cellular component group; and catalytic activity, oxidoreductase activity, and protein binding in the molecular function group. Small molecule metabolic, organophosphate metabolic, and coenzyme metabolic processes in the biological process group; the cytoplasm in the cellular component group; and catalytic activity and GTPase binding in the molecular function group were significant at 24 hpi. Then, KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis was performed, and 82 and 99 pathways (15 and 17, p-value < 0.05) were found to be enriched at 12 and 24 hpi, respectively. Selenocompound metabolism, ubiquinone, and other terpenoid-quinone biosyntheses, fatty acid biosynthesis, lysine biosynthesis, and the citrate cycle were the top five significant pathways at 12 hpi. Glycolysis/gluconeogenesis, biosynthesis of secondary metabolites, linoleic acid metabolism, pyruvate metabolism, and cyanoamino acid metabolism were the top five at 24 hpi. Some V. dahliae-resistance-related proteins, including phenylpropanoid-pathway-related proteins, stress and defense response proteins, plant-pathogen interaction pathway and pathogenesis-related proteins, cell wall organization and reinforcement-related proteins, phytohormones-signal-pathways-related proteins, and other defense-related proteins were identified. In conclusion, this is the first proteomic analysis of S. sisymbriifolium under V. dahliae stress.

Phylotranscriptomics reveal the spatio-temporal distribution and morphological evolution of Macrozamia, an Australian endemic genus of Cycadales

BACKGROUND AND AIMS

Cycads are regarded as an ancient lineage of living seed plants, and hold important clues to understand the early evolutionary trends of seed plants. The molecular phylogeny and spatio-temporal diversification of one of the species-rich genera of cycads, Macrozamia, have not been well reconstructed.

METHODS

We analysed a transcriptome dataset of 4740 single-copy nuclear genes (SCGs) of 39 Macrozamia species and two outgroup taxa. Based on concatenated (maximum parsimony, maximum likelihood) and multispecies coalescent analyses, we first establish a well-resolved phylogenetic tree of Macrozamia. To identify cyto-nuclear incongruence, the plastid protein coding genes (PCGs) from transcriptome data are extracted using the software HybPiper. Furthermore, we explore the biogeographical history of the genus and shed light on the pattern of floristic exchange between three distinct areas of Australia. Six key diagnostic characters are traced on the phylogenetic framework using two comparative methods, and infra-generic classification is investigated.

KEY RESULTS

The tree topologies of concatenated and multi-species coalescent analyses of SCGs are mostly congruent with a few conflicting nodes, while those from plastid PCGs show poorly supported relationships. The genus contains three major clades that correspond to their distinct distributional areas in Australia. The crown group of Macrozamia is estimated to around 11.80 Ma, with a major expansion in the last 5-6 Myr. Six morphological characters show homoplasy, and the traditional phenetic sectional division of the genus is inconsistent with this current phylogeny.

CONCLUSIONS

This first detailed phylogenetic investigation of Macrozamia demonstrates promising prospects of SCGs in resolving phylogenetic relationships within cycads. Our study suggests that Macrozamia, once widely distributed in Australia, underwent major extinctions because of fluctuating climatic conditions such as cooling and mesic biome disappearance in the past. The current close placement of morphologically distinct species in the phylogenetic tree may be related to neotenic events that occurred in the genus.

Transcriptomic Analysis of Resistant and Susceptible Responses in a New Model Root-Knot Nematode Infection System Using Solanum torvum and Meloidogyne arenaria

Root-knot nematodes (RKNs) are among the most devastating pests in agriculture. Solanum torvum Sw. (Turkey berry) has been used as a rootstock for eggplant (aubergine) cultivation because of its resistance to RKNs, including Meloidogyne incognita and M. arenaria. We previously found that a pathotype of M. arenaria, A2-J, is able to infect and propagate in S. torvum. In vitro infection assays showed that S. torvum induced the accumulation of brown pigments during avirulent pathotype A2-O infection, but not during virulent A2-J infection. This experimental system is advantageous because resistant and susceptible responses can be distinguished within a few days, and because a single plant genome can yield information about both resistant and susceptible responses. Comparative RNA-sequencing analysis of S. torvum inoculated with A2-J and A2-O at early stages of infection was used to parse the specific resistance and susceptible responses. Infection with A2-J did not induce statistically significant changes in gene expression within one day post-inoculation (DPI), but afterward, A2-J specifically induced the expression of chalcone synthase, spermidine synthase, and genes related to cell wall modification and transmembrane transport. Infection with A2-O rapidly induced the expression of genes encoding class III peroxidases, sesquiterpene synthases, and fatty acid desaturases at 1 DPI, followed by genes involved in defense, hormone signaling, and the biosynthesis of lignin at 3 DPI. Both isolates induced the expression of suberin biosynthetic genes, which may be triggered by wounding during nematode infection. Histochemical analysis revealed that A2-O, but not A2-J, induced lignin accumulation at the root tip, suggesting that physical reinforcement of cell walls with lignin is an important defense response against nematodes. The S. torvum-RKN system can provide a molecular basis for understanding plant-nematode interactions.

Transcriptome Analysis of Eggplant Root in Response to Root-Knot Nematode Infection

Eggplant (Solanum melongena L.), which belongs to the Solanaceae family, is an important vegetable crop. However, its production is severely threatened by root-knot nematodes (RKNs) in many countries. Solanum torvum, a wild relative of eggplant, is employed worldwide as rootstock for eggplant cultivation due to its resistance to soil-borne diseases such as RKNs. In this study, to identify the RKN defense mechanisms, the transcriptomic profiles of eggplant and Solanum torvum were compared. A total of 5360 differentially expressed genes (DEGs) were identified for the response to RKN infection. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis showed that these DEGs are mainly involved in the processes of response to stimulus, protein phosphorylation, hormone signal transduction, and plant-pathogen interaction pathways. Many phytohormone-related genes and transcription factors (MYB, WRKY, and NAC) were differentially expressed at the four time points (ck, 7, 14, and 28 days post-infection). The abscisic acid signaling pathway might be involved in plant-nematode interactions. qRT-PCR validated the expression levels of some of the DEGs in eggplant. These findings demonstrate the nematode-induced expression profiles and provide some insights into the nematode resistance mechanism in eggplant.

Genetic engineering of eggplant accumulating β-carotene in fruit

Genetic engineering of eggplant using fruit-specific EEF48 promoter-driven bacterial PSY gene, crtB, confers β-carotene accumulation in fruit. Eggplant (Solanum melongena L.) is globally cultivated especially in Asia and is an important source of nutrients in the diets of low-income consumers in developing countries. Since fruits of eggplant have low provitamin A carotenoid content, it is expected to develop eggplant with high carotenoid content for combatting vitamin A deficiency. To achieve this, the present study implemented a metabolic engineering strategy to modify the carotenoid biosynthetic pathway in eggplant. Expression analysis of carotenogenic genes in eggplant tissues showed that the expression of the endogenous phytoene synthase (PSY) was low in fruit and callus. Orange-colored calluses were generated from ectopic expression of crtB gene, which encodes bacterial PSY, in eggplant cells. The orange calluses accumulated > 20 μg g^-1 FW of β-carotene, which was approximately 150-fold higher than that of the untransformed calluses. These observations suggest that the PSY expression is the rate-limiting step for β-carotene production in callus and fruit. Since the orange calluses did not regenerate plants, we chose eggplant EEF48 gene, which is presumably expressed in fruit. We amplified its promoter region by TAIL-PCR and showed that the EEF48 promoter is indeed active in eggplant fruit. Subsequently, transgenic eggplant lines having EEF48 promoter-driven crtB were produced. Among the transgenic lines produced, one line set fruit containing 1.50 μg g^-1 FW of β-carotene, which was 30-fold higher than that of the untransformed fruits (0.05 μg g^-1 FW). The self-pollinated progenies showed a 3:1 segregation ratio for the presence and absence of the transgene, which was linked to the β-carotene accumulation in fruit. These results provide a strategy for improvement of carotenoid content in eggplant fruit.

Comparison of leaf transcriptome in response to Rhizoctonia solani infection between resistant and susceptible rice cultivars

Background

Sheath blight (SB), caused by Rhizoctonia solani, is a common rice disease worldwide. Currently, rice cultivars with robust resistance to R. solani are still lacking. To provide theoretic basis for molecular breeding of R. solani-resistant rice cultivars, the changes of transcriptome profiles in response to R. solani infection were compared between a moderate resistant cultivar (Yanhui-888, YH) and a susceptible cultivar (Jingang-30, JG).

Results

In the present study, 3085 differentially express genes (DEGs) were detected between the infected leaves and the control in JG, with 2853 DEGs in YH. A total of 4091 unigenes were significantly upregulated in YH than in JG before infection, while 3192 were significantly upregulated after infection. Further analysis revealed that YH and JG showed similar molecular responses to R. solani infection, but the responses were earlier in JG than in YH. Expression levels of trans-cinnamate 4-monooxygenase (C4H), ethylene-insensitive protein 2 (EIN2), transcriptome factor WRKY33 and the KEGG pathway plant-pathogen interaction were significantly affected by R. solani infection. More importantly, these components were all over-represented in YH cultivar than in JG cultivar before and/or after infection.

Conclusions

These genes possibly contribute to the higher resistance of YH to R. solani than JG and were potential target genes to molecularly breed R. solani-resistant rice cultivar.

The Complete Chloroplast Genomes of Two Lespedeza Species: Insights into Codon Usage Bias, RNA Editing Sites, and Phylogenetic Relationships in Desmodieae (Fabaceae: Papilionoideae)

The genus Lespedeza (tribe: Desmodieae) consists of about 40 species that have high medicinal and economic value. However, in this genus, using morphological characters, the species identification is quite complicated, which can be solved by the analysis of the complete chloroplast genomes. As primary organelle genomes, the complete genome sequences of chloroplasts (cp) provide unique molecular information to study the divergence of species, RNA editing, and phylogeny. Therefore, to the best of our knowledge, for the first time, we sequenced the complete cp genomes of two representative Lespedeza species: Lespedeza davurica and Lespedeza cuneata. The cp genomes of both the species were found to be 149,010 bp in length, exhibiting the typical angiosperm chloroplast structure containing four regions. The Lespedeza cp genomes showed similar conserved gene contents, order, and orientations with a total GC content of 35.0%. A total of 128 genes, including 83 protein-coding genes, 37 tRNAs, and eight rRNAs, were identified from each genome. Unique molecular features of the two Lespedeza cp genome sequences were obtained by performing the analysis of repeats, sequence divergence, codon usage, and predicting the RNA editing sites in addition to phylogenetic analysis with other key genera in tribe Desmodieae. Using the two datasets, the phylogenetic relationship of Lespedeza species among Deasmodieae was discovered, suggesting that whole cp genomes provided useful information for phylogenetic studies of these species.

De novo assembly and discovery of genes involved in the response of Solanum sisymbriifolium to Verticillium dahlia

Verticillium wilt, caused by the soil-borne fungus Verticillium dahliae, is a devastating disease of eggplant (Solanum spp.) and causes substantial losses worldwide. Although some genes or biological processes involved in the interaction between eggplant and V. dahliae have been identified in some studies, the underlying molecular mechanism is not yet clear. Here, we monitored the transcriptomic profiles of the roots of resistant S. sisymbriifolium plants challenged with V. dahliae. Based on the measurements of physiological indexes (T-SOD, POD and SSs), three time points were selected and subsequently divided into two stages (S_12 h vs. S_0 h and S_48 h vs. S_12 h). KEGG enrichment analysis of the DEGs revealed several genes putatively involved in regulating plant-V. dahliae interactions, including mitogen-activated protein kinase (MAPK) genes (MEKK1 and MAP2K1), WRKY genes (WRKY22 and WRKY33) and cytochrome P450 (CYP) genes (CYP73A/C4H, CYP98A/C3'H and CYP84A/F5H). In addition, a subset of genes that play an important role in activating V. dahliae defence responses, including Ve genes as well as genes encoding PR proteins and TFs, were screened and are discussed. These results will help to identify key resistance genes and will contribute to a further understanding of molecular mechanisms of the S. sisymbriifolium resistance response to V. dahliae.

Comparative Transcriptome Analysis in Eggplant Reveals Selection Trends during Eggplant Domestication

Eggplant (Solanum melongena L.) is an economically and nutritionally important fruit crop of the Solanaceae family, which was domesticated in India and southern China. However, the genome regions subjected to selective sweeps in eggplant remain unknown. In the present study, we performed comparative transcriptome analysis of cultivated and wild eggplant species with emphasis on the selection pattern during domestication. In total, 44,073 (S. sisymbriifolium) to 58,677 (S. melongena cultivar S58) unigenes were generated for the six eggplant accessions with total lengths of 36.6-46 Mb. The orthologous genes were assessed using the ratio of nonsynonymous (K a) to synonymous (K s) nucleotide substitutions to characterize selective patterns during eggplant domestication. We identified 19 genes under positive selection across the phylogeny that were classified into four groups. The gene (OG12205) under positive selection was possibly associated with fruit-related traits in eggplant, which may have resulted from human manipulation. Eight positive selected genes were potentially involved in stress tolerance or disease resistance, suggesting that environmental changes and biotic stresses were important selective pressures in eggplant domestication. Taken together, our results shed light on the effects of artificial and natural selection on the transcriptomes of eggplant and its wild relatives. Identification of the selected genes will facilitate the understanding of genetic architecture of domesticated-related traits and provide resources for resistant breeding in eggplant.

Whole-Genome Resequencing of Seven Eggplant (Solanum melongena) and One Wild Relative (S. incanum) Accessions Provides New Insights and Breeding Tools for Eggplant Enhancement

Whole-genome resequencing provides information of great relevance for crop genetics, evolution, and breeding. Here, we present the first whole-genome resequencing study using seven eggplant (Solanum melongena) and one wild relative (Solanum incanum) accessions. These eight accessions were selected for displaying a high phenotypic and genetic diversity and for being the founder parents of an eggplant multiparent advanced generation intercrosses population. By resequencing at an average depth of 19.8× and comparing to the high-quality reference genome "67/3" over 10 million high-reliable polymorphisms were discovered, of which over 9 million (84.5%) were single nucleotide polymorphisms and more than 700,000 (6.5%) InDels. However, while for the S. melongena accessions, the variants identified ranged from 0.8 to 1.3 million, over 9 million were detected for the wild S. incanum. This confirms the narrow genetic diversity of the domesticated eggplant and suggests that introgression breeding using wild relatives can efficiently contribute to broadening the genetic basis of this crop. Differences were observed among accessions for the enrichment in genes regulating important biological processes. By analyzing the distribution of the variants, we identified the potential footprints of old introgressions from wild relatives that can help to unravel the unclear domestication and breeding history. The comprehensive annotation of these eight genomes and the information provided in this study represents a landmark in eggplant genomics and allows the development of tools for eggplant genetics and breeding.

Cultivar-Dependent Responses of Eggplant (Solanum melongena L.) to Simultaneous Verticillium dahliae Infection and Drought

Several studies regarding the imposition of stresses simultaneously in plants have shown that plant responses are different under individual and combined stress. Pathogen infection in combination with drought can act both additively and antagonistically, suggesting a tailored-made plant response to these stresses. The aforementioned combination of stresses can be considered as one of the most important factors affecting global crop production. In the present research we studied eggplant responses to simultaneous Verticillium dahliae infection and drought with respect to the application of the individual stresses alone and investigated the extent to which these responses were cultivar dependent. Two eggplant cultivars (Skoutari and EMI) with intermediate resistance to V. dahliae were subjected to combined stress for a 3-week period. Significant differences in plant growth, several physiological and biochemical parameters (photosynthesis rate, leaf gas exchanges, Malondialdehyde, Proline) and gene expression, were found between plants subjected to combined and individual stresses. Furthermore, plant growth and molecular (lipid peroxidation, hydrogen peroxide, gene expression levels) changes highlight a clear discrimination between the two cultivars in response to simultaneous V. dahliae infection and drought. Our results showed that combined stress affects significantly plants responses compared to the application of individual stresses alone and that these responses are cultivar dependent.

Assessment of an Organ-Specific de Novo Transcriptome of the Nematode Trap-Crop, Solanum sisymbriifolium

Solanum sisymbriifolium, also known as “Litchi Tomato” or “Sticky Nightshade,” is an undomesticated and poorly researched plant related to potato and tomato. Unlike the latter species, S. sisymbriifolium induces eggs of the cyst nematode, Globodera pallida, to hatch and migrate into its roots, but then arrests further nematode maturation. In order to provide researchers with a partial blueprint of its genetic make-up so that the mechanism of this response might be identified, we used single molecule real time (SMRT) sequencing to compile a high quality de novo transcriptome of 41,189 unigenes drawn from individually sequenced bud, root, stem, and leaf RNA populations. Functional annotation and BUSCO analysis showed that this transcriptome was surprisingly complete, even though it represented genes expressed at a single time point. By sequencing the 4 organ libraries separately, we found we could get a reliable snapshot of transcript distributions in each organ. A divergent site analysis of the merged transcriptome indicated that this species might have undergone a recent genome duplication and re-diploidization. Further analysis indicated that the plant then retained a disproportionate number of genes associated with photosynthesis and amino acid metabolism in comparison to genes with characteristics of R-proteins or involved in secondary metabolism. The former processes may have given S. sisymbriifolium a bigger competitive advantage than the latter did.

Chloroplast genome analyses and genomic resource development for epilithic sister genera Oresitrophe and Mukdenia (Saxifragaceae), using genome skimming data

Background

Epilithic sister genera Oresitrophe and Mukdenia (Saxifragaceae) have an epilithic habitat (rocky slopes) and a parapatric distribution in East Asia, which makes them an ideal model for a more comprehensive understanding of the demographic and divergence history and the influence of climate changes in East Asia. However, the genetic background and resources for these two genera are scarce.

Results

The complete chloroplast (cp) genomes of two Oresitrophe rupifraga and one Mukdenia rossii individuals were reconstructed and comparative analyses were conducted to examine the evolutionary pattern of chloroplast genomes in Saxifragaceae. The cp genomes ranged from 156,738 bp to 156,960 bp in length and had a typical quadripartite structure with a conserved genome arrangement. Comparative analysis revealed the intron of rpl2 has been lost in Heuchera parviflora, Tiarella polyphylla, M. rossii and O. rupifraga but presents in the reference genome of Penthorum chinense. Seven cp hotspot regions (trnH-psbA, trnR-atpA, atpI-rps2, rps2-rpoC2, petN-psbM, rps4-trnT and rpl33-rps18) were identified between Oresitrophe and Mukdenia, while four hotspots (trnQ-psbK, trnR-atpA, trnS-psbZ and rpl33-rps18) were identified within Oresitrophe. In addition, 24 polymorphic cpSSR loci were found between Oresitrophe and Mukdenia. Most importantly, we successfully developed 126 intergeneric polymorphic gSSR markers between Oresitrophe and Mukdenia, as well as 452 intrageneric ones within Oresitrophe. Twelve randomly selected intergeneric gSSRs have shown that these two genera exhibit a significant genetic structure.

Conclusions

In this study, we conducted genome skimming for Oresitrophe rupifraga and Mukdenia rossii. Using these data, we were able to not only assemble their complete chloroplast genomes, but also develop abundant genetic resources (cp hotspots, cpSSRs, polymorphic gSSRs). The genomic patterns and genetic resources presented here will contribute to further studies on population genetics, phylogeny and conservation biology in Saxifragaceae.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4633-x) contains supplementary material, which is available to authorized users.

Comprehensive Characterization of Simple Sequence Repeats in Eggplant (Solanum melongena L.) Genome and Construction of a Web Resource

We have characterized the simple sequence repeat (SSR) markers of the eggplant (Solanum melongena) using a recent high quality sequence of its whole genome. We found nearly 133,000 perfect SSRs, a density of 125.5 SSRs/Mbp, and also about 178,400 imperfect SSRs. Of the perfect SSRs, 15.6% were complex, with two stretches of repeats separated by an intervening block of <100 nt. Di- and trinucleotide SSRs accounted, respectively, for 43 and 37% of the total. The SSRs were classified according to their number of repeats and overall length, and were assigned to their linkage group. We found 2,449 of the perfect SSRs in 2,086 genes, with an overall density of 18.5 SSRs/Mbp across the gene space; 3,524 imperfect SSRs were present in 2,924 genes at a density of 26.7 SSRs/Mbp. Putative functions were assigned via ontology to genes containing at least one SSR. Using this data we developed an "Eggplant Microsatellite DataBase" (EgMiDB) which permits identification of SSR markers in terms of their location on the genome, type of repeat (perfect vs. imperfect), motif type, sequence, repeat number and genomic/gene context. It also suggests forward and reverse primers. We employed an in silico PCR analysis to validate these SSR markers, using as templates two CDS sets and three assembled transcriptomes obtained from diverse eggplant accessions.

Development and Genetic Characterization of Advanced Backcross Materials and An Introgression Line Population of Solanum incanum in a S. melongena Background

Advanced backcrosses (ABs) and introgression lines (ILs) of eggplant (Solanum melongena) can speed up genetics and genomics studies and breeding in this crop. We have developed the first full set of ABs and ILs in eggplant using Solanum incanum, a wild eggplant that has a relatively high tolerance to drought, as a donor parent. The development of these ABs and IL eggplant populations had a low efficiency in the early stages, because of the lack of molecular markers and genomic tools. However, this dramatically improved after performing genotyping-by-sequencing in the first round of selfing, followed by high-resolution-melting single nucleotide polymorphism genotyping in subsequent selection steps. A set of 73 selected ABs covered 99% of the S. incanum genome, while 25 fixed immortal ILs, each carrying a single introgressed fragment in homozygosis, altogether spanned 61.7% of the S. incanum genome. The introgressed size fragment in the ILs contained between 0.1 and 10.9% of the S. incanum genome, with a mean value of 4.3%. Sixty-eight candidate genes involved in drought tolerance were identified in the set of ILs. This first set of ABs and ILs of eggplant will be extremely useful for the genetic dissection of traits of interest for eggplant, and represents an elite material for introduction into the breeding pipelines for developing new eggplant cultivars adapted to the challenges posed by the climate-change scenario.

Transcriptome analyses provide insights into the difference of alkaloids biosynthesis in the Chinese goldthread (Coptis chinensis Franch.) from different biotopes

Coptis chinensis Franch., the Chinese goldthread (‘Weilian’ in Chinese), one of the most important medicinal plants from the family Ranunculaceae, and its rhizome has been widely used in Traditional Chinese Medicine for centuries. Here, we analyzed the chemical components and the transcriptome of the Chinese goldthread from three biotopes, including Zhenping, Zunyi and Shizhu. We built comprehensive, high-quality de novo transcriptome assemblies of the Chinese goldthread from short-read RNA-Sequencing data, obtaining 155,710 transcripts and 56,071 unigenes. More than 98.39% and 95.97% of core eukaryotic genes were found in the transcripts and unigenes respectively, indicating that this unigene set capture the majority of the coding genes. A total of 520,462, 493,718, and 507,247 heterozygous SNPs were identified in the three accessions from Zhenping, Zunyi, and Shizhu respectively, indicating high polymorphism in coding regions of the Chinese goldthread (∼1%). Chemical analyses of the rhizome identified six major components, including berberine, palmatine, coptisine, epiberberine, columbamine, and jatrorrhizine. Berberine has the highest concentrations, followed by coptisine, palmatine, and epiberberine sequentially for all the three accessions. The drug quality of the accession from Shizhu may be the highest among these accessions. Differential analyses of the transcriptome identified four pivotal candidate enzymes, including aspartate aminotransferaseprotein, polyphenol oxidase, primary-amine oxidase, and tyrosine decarboxylase, were significantly differentially expressed and may be responsible for the difference of alkaloids contents in the accessions from different biotopes.

Comparative Analysis of the Complete Chloroplast Genome of Four Endangered Herbals of Notopterygium

Notopterygium H. de Boissieu (Apiaceae) is an endangered perennial herb endemic to China. A good knowledge of phylogenetic evolution and population genomics is conducive to the establishment of effective management and conservation strategies of the genus Notopterygium. In this study, the complete chloroplast (cp) genomes of four Notopterygium species (N. incisum C. C. Ting ex H. T. Chang, N. oviforme R. H. Shan, N. franchetii H. de Boissieu and N. forrestii H. Wolff) were assembled and characterized using next-generation sequencing. We investigated the gene organization, order, size and repeat sequences of the cp genome and constructed the phylogenetic relationships of Notopterygium species based on the chloroplast DNA and nuclear internal transcribed spacer (ITS) sequences. Comparative analysis of plastid genome showed that the cp DNA are the standard double-stranded molecule, ranging from 157,462 bp (N. oviforme) to 159,607 bp (N. forrestii) in length. The circular DNA each contained a large single-copy (LSC) region, a small single-copy (SSC) region, and a pair of inverted repeats (IRs). The cp DNA of four species contained 85 protein-coding genes, 37 transfer RNA (tRNA) genes and 8 ribosomal RNA (rRNA) genes, respectively. We determined the marked conservation of gene content and sequence evolutionary rate in the cp genome of four Notopterygium species. Three genes (psaI, psbI and rpoA) were possibly under positive selection among the four sampled species. Phylogenetic analysis showed that four Notopterygium species formed a monophyletic clade with high bootstrap support. However, the inconsistent interspecific relationships with the genus Notopterygium were identified between the cp DNA and ITS markers. The incomplete lineage sorting, convergence evolution or hybridization, gene infiltration and different sampling strategies among species may have caused the incongruence between the nuclear and cp DNA relationships. The present results suggested that Notopterygium species may have experienced a complex evolutionary history and speciation process.

Characterization of lignocellulolytic activities from fungi isolated from the deep-sea sponge Stelletta normani

Extreme habitats have usually been regarded as a source of microorganisms that possess robust proteins that help enable them to survive in such harsh conditions. The deep sea can be considered an extreme habitat due to low temperatures (<5°C) and high pressure, however marine sponges survive in these habitats. While bacteria derived from deep-sea marine sponges have been studied, much less information is available on fungal biodiversity associated with these sponges. Following screening of fourteen fungi isolated from the deep-sea sponge Stelletta normani sampled at a depth of 751 metres, three halotolerant strains (TS2, TS11 and TS12) were identified which displayed high CMCase and xylanase activities. Molecular based taxonomic approaches identified these strains as Cadophora sp. TS2, Emericellopsis sp. TS11 and Pseudogymnoascus sp. TS 12. These three fungi displayed psychrotolerance and halotolerant growth on CMC and xylan as sole carbon sources, with optimal growth rates at 20°C. They produced CMCase and xylanase activities, which displayed optimal temperature and pH values of between 50-70°C and pH 5-8 respectively, together with good thermostability and halotolerance. In solid-state fermentations TS2, TS11 and TS12 produced CMCases, xylanases and peroxidase/phenol oxidases when grown on corn stover and wheat straw. This is the first time that CMCase, xylanase and peroxidase/phenol oxidase activities have been reported in these three fungal genera isolated from a marine sponge. Given the biochemical characteristics of these ligninolytic enzymes it is likely that they may prove useful in future biomass conversion strategies involving lignocellulosic materials.

A large set of 26 new reference transcriptomes dedicated to comparative population genomics in crops and wild relatives

We produced a unique large data set of reference transcriptomes to obtain new knowledge about the evolution of plant genomes and crop domestication. For this purpose, we validated a RNA-Seq data assembly protocol to perform comparative population genomics. For the validation, we assessed and compared the quality of de novo Illumina short-read assemblies using data from two crops for which an annotated reference genome was available, namely grapevine and sorghum. We used the same protocol for the release of 26 new transcriptomes of crop plants and wild relatives, including still understudied crops such as yam, pearl millet and fonio. The species list has a wide taxonomic representation with the inclusion of 15 monocots and 11 eudicots. All contigs were annotated using BLAST, prot4EST and Blast2GO. A strong originality of the data set is that each crop is associated with close relative species, which will permit whole-genome comparative evolutionary studies between crops and their wild-related species. This large resource will thus serve research communities working on both crops and model organisms. All the data are available at http://arcad-bioinformatics.southgreen.fr/.

Transcriptomic analysis of heteromorphic stamens in Cassia biscapsularis L

Hermaphroditic flowers have evolved primarily under the selection on male function. Evolutionary modification often leads to stamen differentiation within flowers, or “heteranthery”, a phenomenon intrigued scientists since the 18^th century until recently. However, the genetic basis and molecular regulation mechanism has barely been touched. Here we conducted comparative transcriptome profiling in Cassia biscapsularis L., a heterantherous species with representative patterns of stamen differentiation. Numerous differentially expressed genes (DEGs) were detected between the staminodes (the degenerated stamens) and fertile stamens, while much fewer genes differentially expressed among the three sets of fertile stamens. GO term enrichment and KEGG pathway analysis characterized functional properties of DEGs in different stamen types. Transcripts showing close correlation between expression pattern and stamen types were identified. Transcription factors from the bHLH family were suggested to have taken crucial part in the formation of staminodes. This first global transcriptomic analysis focusing on stamen differentiation opens the door toward a more comprehensive understanding on the molecular regulation of floral organ evolution. Especially, the generated unigene resource would be valuable for developing male sterile lines in agronomy.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2631-4) contains supplementary material, which is available to authorized users.

Comparative Transcriptome and Chloroplast Genome Analyses of Two Related Dipteronia Species

Dipteronia (order Sapindales) is an endangered genus endemic to China and has two living species, D.sinensis and D. dyeriana. The plants are closely related to the genus Acer, which is also classified in the order Sapindales. Evolutionary studies on Dipteronia have been hindered by the paucity of information on their genomes and plastids. Here, we used next generation sequencing to characterize the transcriptomes and complete chloroplast genomes of both Dipteronia species. A comparison of the transcriptomes of both species identified a total of 7814 orthologs. Estimation of selection pressures using Ka/Ks ratios showed that only 30 of 5435 orthologous pairs had a ratio significantly >1, i.e., showing positive selection. However, 4041 orthologs had a Ka/Ks < 0.5 (p < 0.05), suggesting that most genes had likely undergone purifying selection. Based on orthologous unigenes, 314 single copy nuclear genes (SCNGs) were identified. Through a combination of de novo and reference guided assembly, plastid genomes were obtained; that of D. sinensis was 157,080 bp and that of D. dyeriana was 157,071 bp. Both plastid genomes encoded 87 protein coding genes, 40 tRNAs, and 8 rRNAs; no significant differences were detected in the size, gene content, and organization of the two plastomes. We used the whole chloroplast genomes to determine the phylogeny of D. sinensis and D. dyeriana and confirmed that the two species were highly divergent. Overall, our study provides comprehensive transcriptomic and chloroplast genomic resources, which will be valuable for future evolutionary studies of Dipteronia.

A novel workflow correlating RNA-seq data to Phythophthora infestans resistance levels in wild Solanum species and potato clones

Comparative transcriptomics between species can provide valuable understanding of plant-pathogen interactions. Here, we focus on wild Solanum species and potato clones with varying degree of resistance against Phytophthora infestans, which causes the devastating late blight disease in potato. The transcriptomes of three wild Solanum species native to Southern Sweden, Solanum dulcamara, Solanum nigrum, and Solanum physalifolium were compared to three potato clones, Desiree (cv.), SW93-1015 and Sarpo Mira. Desiree and S. physalifolium are susceptible to P. infestans whereas the other four have different degrees of resistance. By building transcript families based on de novo assembled RNA-seq across species and clones and correlating these to resistance phenotypes, we created a novel workflow to identify families with expanded or depleted number of transcripts in relation to the P. infestans resistance level. Analysis was facilitated by inferring functional annotations based on the family structure and semantic clustering. More transcript families were expanded in the resistant clones and species and the enriched functions of these were associated to expected gene ontology (GO) terms for resistance mechanisms such as hypersensitive response, host programmed cell death and endopeptidase activity. However, a number of unexpected functions and transcripts were also identified, for example transmembrane transport and protein acylation expanded in the susceptible group and a cluster of Zinc knuckle family proteins expanded in the resistant group. Over 400 expressed putative resistance (R-)genes were identified and resistant clones Sarpo Mira and SW93-1015 had ca 25% more expressed putative R-genes than susceptible cultivar Desiree. However, no differences in numbers of susceptibility (S-)gene homologs were seen between species and clones. In addition, we identified P. infestans transcripts including effectors in the early stages of P. infestans-Solanum interactions.

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.

Transcriptome characterization of three wild Chinese Vitis uncovers a large number of distinct disease related genes

BACKGROUND

Grape is one of the most valuable fruit crops and can serve for both fresh consumption and wine production. Grape cultivars have been selected and evolved to produce high-quality fruits during their domestication over thousands of years. However, current widely planted grape cultivars suffer extensive loss to many diseases while most wild species show resistance to various pathogens. Therefore, a comprehensive evaluation of wild grapes would contribute to the improvement of disease resistance in grape breeding programs.

RESULTS

We performed deep transcriptome sequencing of three Chinese wild grapes using the Illumina strand-specific RNA-Seq technology. High quality transcriptomes were assembled de novo and more than 93% transcripts were shared with the reference PN40024 genome. Over 1,600 distinct transcripts, which were absent or highly divergent from sequences in the reference PN40024 genome, were identified in each of the three wild grapes, among which more than 1,000 were potential protein-coding genes. Gene Ontology (GO) and pathway annotations of these distinct genes showed those involved in defense responses and plant secondary metabolisms were highly enriched. More than 87,000 single nucleotide polymorphisms (SNPs) and 2,000 small insertions or deletions (indels) were identified between each genotype and PN40024, and approximately 20% of the SNPs caused nonsynonymous mutations. Finally, we discovered 100 to 200 highly confident cis-natural antisense transcript (cis-NAT) pairs in each genotype. These transcripts were significantly enriched with genes involved in secondary metabolisms and plant responses to abiotic stresses.

CONCLUSION

The three de novo assembled transcriptomes provide a comprehensive sequence resource for molecular genetic research in grape. The newly discovered genes from wild Vitis, as well as SNPs and small indels we identified, may facilitate future studies on the molecular mechanisms related to valuable traits possessed by these wild Vitis and contribute to the grape breeding programs. Furthermore, we identified hundreds of cis-NAT pairs which showed their potential regulatory roles in secondary metabolism and abiotic stress responses.