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

De novo assembly of the complete mitochondrial genome of pepino (Solanum muricatum) using PacBio HiFi sequencing: insights into structure, phylogenetic implications, and RNA editing

BACKGROUND

Solanum muricatum is an emerging horticultural fruit crop with rich nutritional and antioxidant properties. Although the chromosome-scale genome of this species has been sequenced, its mitochondrial genome sequence has not been reported to date.

RESULTS

PacBio HiFi sequencing was used to assemble the circular mitogenome of S. muricatum, which was 433,466 bp in length. In total, 38 protein-coding, 19 tRNA, and 3 rRNA genes were annotated. The reticulate mitochondrial conformations with multiple junctions were verified by polymerase chain reaction, and codon usage, sequence repeats, and gene migration from chloroplast to mitochondrial genome were determined. A collinearity analysis of eight Solanum mitogenomes revealed high structural variability. Overall, 585 RNA editing sites in protein coding genes were identified based on RNA-seq data. Among them, mttB was the most frequently edited (52 times), followed by ccmB (46 times). A phylogenetic analysis based on the S. muricatum mitogenome and those of 39 other taxa (including 25 Solanaceae species) revealed the evolutionary and taxonomic status of S. muricatum.

CONCLUSIONS

We provide the first report of the assembled and annotated S. muricatum mitogenome. This information will help to lay the groundwork for future research on the evolutionary biology of Solanaceae species. Furthermore, the results will assist the development of molecular breeding strategies for S. muricatum based on the most beneficial agronomic traits of this species.

Metabolomics analysis of quality components metabolism during the growth process of pepino (Solanum muricatum) fruit

Pepino (Solanum muricatum), a horticultural crop that has experienced significant growth in the highlands of China over the past two decades, is widely embraced by consumers due to its distinctive taste and nutritional advantages. This study focused on the cultivar 'Qingcanxiang' of pepino grown on the Qinghai-Tibetan Plateau was analyzed using UPLC-QTOF-MS and RNA-seq transcriptome sequencing. Fruit samples were collected at three distinct stages of development, and the results of the metabolomics and transcriptomics were compared and correlated. The study's findings indicate that the 'Qingcanxiang' fruit contained a total of 187 metabolites, comprising 12 distinct categories of compounds, including amino acids and their derivatives, organic acids, sugars and alcohols, phenols and phenolic acids. Of these metabolites, 94 were identified as differential. Significant variations in nutrient composition were observed across the three growth stages of the fruit. Specifically, the stage spanning from the growth to the maturation was identified as the critical stages for nutrient accumulation and flavor development. Transcriptome sequencing analysis revealed a set of highly associated genes between aspartate and quinic acid, namely SIR2, IRAK4, RP-L29, and CCNH. These genes are potentially involved in the regulation of both amino acid and phenolic acid synthesis. Through the application of metabolomics and transcriptomics, this investigation elucidates the alterations in metabolites and the underlying molecular regulatory mechanisms of pepino fruits during three growth stages. The findings furnish a theoretical foundation for the evaluation of nutritional quality and the enhancement of breeding strategies for pepino.

Pepino (Solanum muricatum) Metabolic Profiles and Soil Nutrient Association Analysis in Three Growing Sites on the Loess Plateau of Northwestern China

Different soil nutrients affect the accumulation characteristics of plant metabolites. To investigate the differences among the metabolites of pepino grown in greenhouses on the Loess Plateau in northwest China, we investigated the main soil nutrients and their correlation with metabolites. A total of 269 pepino metabolites were identified using UPLC-QTOF-MS to detect metabolites in fruits from three major pepino growing regions and analyze their differential distribution characteristics. A total of 99 of these substances differed among pepino fruits from the three areas, and the main classes of the differential metabolites were, in order of number: amino acids and derivatives, nucleotides and derivatives, organic acids, alkaloids, vitamins, saccharides and alcohols, phenolic acids, lipids and others. An environmental factor analysis identified soil nutrients as the most significant differentiator. Five soil nutrient indicators: TN (total nitrogen), TP (total phosphorus), AP (available phosphorus), AK (available potassium), and OM (organic matter), exhibited significant differences in three growing sites. Metabolite and soil nutrient association analysis using redundancy analysis (RDA) and the Mantel test indicated that TN and OM contributed to the accumulation of amino acids and derivatives, nucleotides and derivatives, and alkaloids while inhibiting organic acids, vitamins coagulation biosynthesis. Moreover, AP and TP were associated with the highest accumulation of saccharides and, alcohols, phenolic acids. Consequently, differences in soil nutrients were reflected in pepino metabolite variability. This study clarified the metabolite variability and the relationship between pepino and soil nutrients in the main planting areas of northwest China. It provides a theoretical basis for the subsequent development of Pepino’s nutritional value and cultivation management.

Revealing the Complete Chloroplast Genome of an Andean Horticultural Crop, Sweet Cucumber (Solanum muricatum), and Its Comparison with Other Solanaceae Species

Sweet cucumber (Solanum muricatum) sect. Basarthrum is a neglected horticultural crop native to the Andean region. It is naturally distributed very close to other two Solanum crops of high importance, potatoes, and tomatoes. To date, molecular tools for this crop remain undetermined. In this study, the complete sweet cucumber chloroplast (cp) genome was obtained and compared with seven Solanaceae species. The cp genome of S. muricatum was 155,681 bp in length and included a large single copy (LSC) region of 86,182 bp and a small single-copy (SSC) region of 18,360 bp, separated by a pair of inverted repeats (IR) regions of 25,568 bp. The cp genome possessed 87 protein-coding genes (CDS), 37 transfer RNA (tRNA) genes, eight ribosomal RNA (rRNA) genes, and one pseudogene. Furthermore, 48 perfect microsatellites were identified. These repeats were mainly located in the noncoding regions. Whole cp genome comparative analysis revealed that the SSC and LSC regions showed more divergence than IR regions. Similar to previous studies, our phylogenetic analysis showed that S. muricatum is a sister species to members of sections Petota + Lycopersicum + Etuberosum. We expect that this first sweet cucumber chloroplast genome will provide potential molecular markers and genomic resources to shed light on the genetic diversity and population studies of S. muricatum, which will allow us to identify varieties and ecotypes. Finally, the features and the structural differentiation will provide us with information about the genes of interest, generating tools for the most precise selection of the best individuals of sweet cucumber, in less time and with fewer resources.

Chromosome-level pepino genome provides insights into genome evolution and anthocyanin biosynthesis in Solanaceae

Pepino (Solanum muricatum, 2n = 2x = 24), a member of the Solanaceae family, is an important globally grown fruit. Herein, we report high-quality, chromosome-level pepino genomes. The 91.67% genome sequence is anchored to 12 chromosomes, with a total length of 1.20 Gb and scaffold N50 of 87.03 Mb. More than half the genome comprises repetitive sequences. In addition to the shared ancient whole-genome triplication (WGT) event in eudicots, an additional new WGT event was present in the pepino. Our findings suggest that pepinos experienced chromosome rearrangements, fusions, and gene loss after a WGT event. The large number of gene removals indicated the instability of Solanaceae genomes, providing opportunities for species divergence and natural selection. The paucity of disease-resistance genes (NBS) in pepino and eggplant has been explained by extensive loss and limited generation of genes after WGT events in Solanaceae. The outbreak of NBS genes was not synchronized in Solanaceae species, which occurred before the Solanaceae WGT event in pepino, tomato, and tobacco, whereas it was almost synchronized with WGT events in the other four Solanaceae species. Transcriptome and comparative genomic analyses revealed several key genes involved in anthocyanin biosynthesis. Although an extra WGT event occurred in Solanaceae, CHS genes related to anthocyanin biosynthesis in grapes were still significantly expanded compared with those in Solanaceae species. Proximal and tandem duplications contributed to the expansion of CHS genes. In conclusion, the pepino genome and annotation facilitate further research into important gene functions and comparative genomic analysis in Solanaceae.

De novo transcriptome assembly and EST-SSR markers development for Zelkova schneideriana Hand.-Mazz. (Ulmaceae)

Zelkova schneideriana Hand.-Mazz. of the Ulmaceae family is a Tertiary relict and economically deciduous tree species endemic to Central and Southern China. In this study, we performed a transcriptome sequencing of Z. schneideriana using high-throughput sequencing approach to detect polymorphic expressed sequence tag-simple sequence repeats (EST-SSR) markers. A total of 3,235 microsatellite loci were detected from 53,517 unigenes. A set of 30 microsatellite markers were randomly selected to validate in 41 individuals from three populations, of which 10 were polymorphic. The number of alleles per locus ranged from 3 to 11. The observed heterozygosity and expected heterozygosity ranged from 0.366 to 0.829 and 0.439 to 0.848, respectively. These polymorphic SSR primers showed good transferability across different Zelkova species, and are valuable for future studies on genetic diversity, conservation, phylogeography, and species delimitation in Z. schneideriana, as well as other Zelkova species.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02968-5.

Transcriptome analysis of two radiated Cycas species and the subsequent species delimitation of the Cycas taiwaniana complex

PREMISE

Cycas is an important gymnosperm genus, and the most diverse of all cycad genera. The C. taiwaniana complex of species are morphologically similar and difficult to distinguish due to a lack of genomic resources.

METHODS

We characterized the transcriptomes of two closely related and endangered Cycas species endemic to Hainan, China: C. hainanensis and C. changjiangensis. Three single-copy nuclear genes in the C. taiwaniana complex were sequenced based on these transcriptomes, enabling us to evaluate the species boundaries using the multispecies coalescent method implemented in the Bayesian Phylogenetics and Phylogeography program.

RESULTS

We obtained 68,184 and 81,561 unigenes for C. changjiangensis and C. hainanensis, respectively. We identified six positively selected genes that are mainly involved in stimulus responses, suggesting that environmental adaptation may have played an important role in the relatively recent divergence of these species. The similar K S distribution peaks at 1.0 observed for the paralogs in the two species indicate a common whole-genome duplication event. Our species delimitation analysis indicated that the C. taiwaniana complex consists of three distinct species, which correspond to the previously reported morphological differences.

DISCUSSION

Our study provides valuable genetic resources for Cycas species and guidance for the taxonomic treatment of the C. taiwaniana complex, as well as new insights into evolution of species within Cycas.

Single Primer Enrichment Technology (SPET) for High-Throughput Genotyping in Tomato and Eggplant Germplasm

Single primer enrichment technology (SPET) is a new, robust, and customizable solution for targeted genotyping. Unlike genotyping by sequencing (GBS), and like DNA chips, SPET is a targeted genotyping technology, relying on the sequencing of a region flanking a primer. Its reliance on single primers, rather than on primer pairs, greatly simplifies panel design, and allows higher levels of multiplexing than PCR-based genotyping. Thanks to the sequencing of the regions surrounding the target SNP, SPET allows the discovery of thousands of closely linked, novel SNPs. In order to assess the potential of SPET for high-throughput genotyping in plants, a panel comprising 5k target SNPs, designed both on coding regions and introns/UTRs, was developed for tomato and eggplant. Genotyping of two panels composed of 400 tomato and 422 eggplant accessions, comprising both domesticated material and wild relatives, generated a total of 12,002 and 30,731 high confidence SNPs, respectively, which comprised both target and novel SNPs in an approximate ratio of 1:1.6, and 1:5.5 in tomato and eggplant, respectively. The vast majority of the markers was transferrable to related species that diverged up to 3.4 million years ago (Solanum pennellii for tomato and S. macrocarpon for eggplant). Maximum Likelihood phylogenetic trees and PCA outputs obtained from the whole dataset highlighted genetic relationships among accessions and species which were congruent with what was previously reported in literature. Better discrimination among domesticated accessions was achieved by using the target SNPs, while better discrimination among wild species was achieved using the whole SNP dataset. Our results reveal that SPET genotyping is a robust, high-throughput technology for genetic fingerprinting, with a high degree of cross-transferability between crops and their cultivated and wild relatives, and allows identification of duplicates and mislabeled accessions in genebanks.

Comparative analysis of two sister Erythrophleum species (Leguminosae) reveal contrasting transcriptome-wide responses to early drought stress

With the ongoing climate change, African rainforests are expected to experience severe drought events in the future. In Africa, the tropical genus Erythrophleum (Fabaceae) includes two forest sister timber tree species displaying contrasting geographical distributions. Erythrophleum ivorense is adapted to wet evergreen Guineo-Congolian forests, whereas E. suaveolens occurs in a wider range of climates, being found in moist dense forests but also in gallery forests under a relatively drier climate. This geographical distribution pattern suggests that the two species might cope differently to drought at the genomic level. Yet, the genetic basis of tolerance response to drought stress in both species is still uncharacterized. To bridge this gap, we performed an RNA-seq approach on seedlings from each species to monitor their transcriptional responses at different levels of drought stress (0, 2 and 6 weeks after stopping watering seedlings). Monitoring of wilting symptoms revealed that E. suaveolens displayed an earlier phenotypic response to drought stress than E. ivorense. At the transcriptomic level, results revealed 2020 (1204 down-regulated/816 up-regulated) and 1495 differentially expressed genes (DEGs) in response to drought stress from a total of 67,432 and 66,605 contigs assembled in E. ivorense and E. suaveolens, respectively. After identifying 30,374 orthologs between species, we found that only 7 of them were DEGs shared between species, while 587 and 458 were differentially expressed only in E. ivorense or E. suaveolens, respectively. GO and KEGG enrichment analysis revealed that the two species differ in terms of significantly regulated pathways as well as the number and expression profile of DEGs (Up/Down) associated with each pathway in the two stress stages. Our results suggested that the two studied species react differently to drought. E. suaveolens seems displaying a prompt response to drought at its early stage strengthened by the down-regulation of many DEGs encoding for signaling and metabolism-related pathways. A considerable up-regulation of these pathways was also found in E. ivorense at the late stage of drought, suggesting this species may be a late responder. Overall, our data may serve as basis for further understanding the genetic control of drought tolerance in tropical trees and favor the selection of crucial genes for genetically enhancing drought resistance.

Comparative transcriptome analysis of the invasive weed Mikania micrantha with its native congeners provides insights into genetic basis underlying successful invasion

Background

Mikania micrantha H.B.K. (Asteraceae) is one of the world’s most invasive weeds which has been rapidly expanding in tropical Asia, including China, while its close relative M. cordata, the only Mikania species native to China, shows no harm to the local ecosystems. These two species are very similar in morphology but differ remarkably in several ecological and physiological traits, representing an ideal system for comparative analysis to investigate the genetic basis underlying invasion success. In this study, we performed RNA-sequencing on the invader M. micrantha and its native congener M. cordata in China, to unravel the genetic basis underlying the strong invasiveness of M. micrantha. For a more robust comparison, another non-invasive congener M. cordifolia was also sequenced and compared.

Results

A total of 52,179, 55,835, and 52,983 unigenes were obtained for M. micrantha, M. cordata, and M. cordifolia, respectively. Phylogenetic analyses and divergence time dating revealed a relatively recent split between M. micrantha and M. cordata, i.e., approximately 4.81 million years ago (MYA), after their divergence with M. cordifolia (8.70 MYA). Gene ontology classifications, pathway assignments and differential expression analysis revealed higher representation or significant up-regulation of genes associated with photosynthesis, energy metabolism, protein modification and stress response in M. micrantha than in M. cordata or M. cordifolia. Analysis of accelerated evolution and positive selection also suggested the importance of these related genes and processes to the adaptability and invasiveness of M. micrantha. Particularly, most (77 out of 112, i.e. 68.75%) positively selected genes found in M. micrantha could be classified into four groups, i.e., energy acquisition and utilization (10 genes), growth and reproduction (13 genes), protection and repair (34 genes), and signal transduction and expression regulation (20 genes), which may have contributed to the high adaptability of M. micrantha to various new environments and the capability to occupy a wider niche, reflected in its high invasiveness.

Conclusions

We characterized the transcriptomes of the invasive species M. micrantha and its non-invasive congeners, M. cordata and M. cordifolia. A comparison of their transcriptomes provided insights into the genetic basis of the high invasiveness of M. micrantha.

Electronic supplementary material

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

Volatile Compound and Gene Expression Analyses Reveal Temporal and Spatial Production of LOX-Derived Volatiles in Pepino (Solanum muricatum Aiton) Fruit and LOX Specificity

Lipoxygenase (LOX) is an important contributor to aroma compounds in most fresh produce; however, little is known about the LOX pathway in pepino (Solanum muricatum Aiton) fruit. We explored the LOX aroma compounds produced by the flesh and the peel and identified eight putative LOX genes expressed in both tissues during fruit growth and development during two consecutive seasons. This study shows that pepino produces C5, C6, and C9 LOX-derived compounds. Odorant C9 volatiles were produced during immature stages with a concomitant decrease when the fruit ripens, whereas C5 and C6 compounds were formed throughout ripening. trans-2-Hexenal and its alcohol were produced in the peel, but not detected in the flesh. The expression of three genes, SmLOXD (putative 13-LOX), SmLOXB, and SmLOX5-like1 (putative 9-LOXs), increased during fruit ripening. These genes may account for aroma volatiles in pepino. Here, we discuss the possible roles of individual LOX genes in pepino.

Whole Transcriptomic Analysis Provides Insights into Molecular Mechanisms for Toxin Biosynthesis in a Toxic Dinoflagellate Alexandrium catenella (ACHK-T)

Paralytic shellfish toxins (PSTs), a group of neurotoxic alkaloids, are the most potent biotoxins for aquatic ecosystems and human health. Marine dinoflagellates and freshwater cyanobacteria are two producers of PSTs. The biosynthesis mechanism of PSTs has been well elucidated in cyanobacteria; however, it remains ambiguous in dinoflagellates. Here, we compared the transcriptome profiles of a toxin-producing dinoflagellate Alexandrium catenella (ACHK-T) at different toxin biosynthesis stages within the cell cycle using RNA-seq. The intracellular toxin content increased gradually in the middle G1 phase and rapidly in the late G1 phase, and then remained relatively stable in other phases. Samples from four toxin biosynthesis stages were selected for sequencing, and finally yielded 110,370 unigenes, of which 66,141 were successfully annotated in the known databases. An analysis of differentially expressed genes revealed that 2866 genes altered significantly and 297 were co-expressed throughout the four stages. These genes participated mainly in protein metabolism, carbohydrate metabolism, and the oxidation-reduction process. A total of 138 homologues of toxin genes were identified, but they altered insignificantly among different stages, indicating that toxin biosynthesis might be regulated translationally or post-translationally. Our results will serve as an important transcriptomic resource to characterize key molecular processes underlying dinoflagellate toxin biosynthesis.