Transcriptome Sequencing and Development of Genic SSR Markers of an Endangered Chinese Endemic Genus Dipteronia Oliver (Aceraceae)

Population transcriptome reveals genetic divergence and expression diversity of medicinally effective ingredients-related genes for Rheum palmatum complex derived from the top-geoherb area

MAIN CONCLUSION

The study revealed genetic diversity and moderate differentiation among the R. palmatum complex within the top-geoherb area. RNA-seq-derived SNP datasets hold the potential to trace geographic origins of the core germplasm. The Rheum palmatum complex, the source plant of rhubarb, has been widely used for centuries due to its diverse functions in clinical treatments. However, the wild resources of rhubarb are currently declining and even facing depletion. Therefore, revealing the genetic background of the R. palmatum complex within the top-geoherb area is important for the efficient utilization and conservation of its wild resources. In this study, population transcriptomic analyses were conducted to assess the genetic diversity and gene expression diversity of different populations in the R. palmatum complex within the top-geoherb area. Candidate single nucleotide polymorphisms (SNPs) were identified as specific molecular markers for tracing the origin of the R. palmatum complex from various top-geoherb areas. Based on the reference genome, a total of 30,480 transcripts and 100,966 SNPs were generated across 82 individuals from 17 populations of the R. palmatum complex. Moderate genetic differentiation was detected for the two genetic lineages of the R. palmatum complex derived from the top-geoherb area. Fourteen genes encoding key enzymes were differentially expressed between two genetic lineages. Besides, 26 specific SNPs located on the genes involved in the biosynthesis of the active components were screened out, and these SNPs were highly differentiated between 2 genetic lineages. A large-scale reference-based assembly transcriptome of the R. palmatum complex from the top-geoherb area provided insights into the genetic divergence and expression differentiation of genetic lineages. The results not only help to understand the genetic background of the R. palmatum complex in the top-geoherb area but also contribute to future genetic conservation and directive breeding of rhubarb.

Genetically- and environmentally-dependent processes drive interspecific and intraspecific divergence in the Chinese relict endemic genus Dipteronia

China is a hotspot of relict plant species that were once widespread throughout the Northern Hemisphere. Recent research has demonstrated that the occurrence of long-term stable refugia in the mountainous regions of central and south-western China allowed their persistence through the late Neogene climate fluctuations. One of these relict lineages is Dipteronia, an oligotypic tree genus with a fossil record extending to the Paleocene. Here, we investigated the genetic variability, demographic dynamics and diversification patterns of the two currently recognized Dipteronia species (D ipteronia sinensis and D . dyeriana). Molecular data were obtained from 45 populations of Dipteronia by genotyping three cpDNA regions, two single copy nuclear genes and 15 simple sequence repeat loci. The genetic study was combined with niche comparison analyses on the environmental space, ecological niche modeling, and landscape connectivity analysis. We found that the two Dipteronia species have highly diverged both in genetic and ecological terms. Despite the incipient speciation processes that can be observed in D. sinensis, the occurrence of long-term stable refugia and, particularly, a dispersal corridor along Daba Shan-west Qinling, likely ensured its genetic and ecological integrity to date. Our study will not only help us to understand how populations of Dipteronia species responded to the tectonic and climatic changes of the Cenozoic, but also provide insight into how Arcto-Tertiary relict plants in East Asia survived, evolved, and diversified.

Transcriptome sequencing of medical herb Salvia Rosmarinus (Rosemary) revealed the phenylpropanoid biosynthesis pathway genes and their phylogenetic relationships

BACKGROUND

The Salvia rosmarinus spenn. (rosemary) is considered an economically important ornamental and medicinal plant and is widely utilized in culinary and for treating several diseases. However, the procedure behind synthesizing secondary metabolites-based bioactive compounds at the molecular level in S. rosmarinus is not explored completely.

METHODS AND RESULTS

We performed transcriptomic sequencing of the pooled sample from leaf and stem tissues on the Illumina HiSeqTM X10 platform. The transcriptomics analysis led to the generation of 29,523,608 raw reads, followed by data pre-processing which generated 23,208,592 clean reads, and de novo assembly of S. rosmarinus obtained 166,849 unigenes. Among them, nearly 75.1% of unigenes i.e., 28,757 were interpreted against a non-redundant protein database. The gene ontology-based annotation classified them into 3 main categories and 55 sub-categories, and clusters of orthologous genes annotation categorized them into 23 functional categories. The Kyoto Encyclopedia of Genes and Genomes database-based pathway analysis confirmed the involvement of 13,402 unigenes in 183 biochemical pathways, among these unigenes, 1,186 are involved in the 17 secondary metabolite production pathways. Several key enzymes involved in producing aromatic amino acids and phenylpropanoids were identified from the transcriptome database. Among the identified 48 families of transcription factors from coding unigenes, bHLH, MYB, WRKYs, NAC, C2H2, C3H, and ERF are involved in flavonoids and other secondary metabolites biosynthesis.

CONCLUSION

The phylogenetic analysis revealed the evolutionary relationship between the phenylpropanoid pathway genes of rosemary with other members of Lamiaceae. Our work reveals a new molecular mechanism behind the biosynthesis of phenylpropanoids and their regulation in rosemary plants.

Genetic variation and structure of endemic and endangered wild celery (Kelussia odoratissima Mozaff.) quantified using novel microsatellite markers developed by next-generation sequencing

Kelussia odoratissima Mozaff. (Apiaceae) is a native plant that has been traditionally consumed in Iran's food and pharmaceutical industries. Overharvesting of the taxon, especially at the beginning of the growing season, due to its considerable medicinal and economic value, is believed to be the main reason for the extirpating of this plant. The consequences of the severe anthropogenic impacts on the genetic diversity of populations are poorly known. In order to investigate the level of genetic variation and patterns of the genetic structure of K. odoratissima, we developed novel microsatellite markers using the 454 Roche next-generation sequencing (NGS) platform for the first time. Out of 1,165 microsatellite markers bioinformatically confirmed, twenty-five were tested, of which 23 were used to screen genetic variation across 12 natural populations. Our results showed that the average number of alleles per locus and the polymorphic information content (PIC) were 10.87 (range 7 to 27), and 0.81 (range 0.67 to 0.94), respectively. The mean observed and expected heterozygosities (± SD) across all populations were 0.80 ± 0.31 and 0.72 ± 0.14, respectively. The average pairwise FST among the populations was 0.37 (range 0.04 to 0.81). Bayesian and distance-based clustering, and principal coordinate analyses revealed at least four major genetic clusters. Although high level of structure can be explained by landscape topography and geographic distance, presence of admixed populations can be associated to seed or pollen dispersal. Contrary to expectations, the high level of genetic variation and lack of inbreeding suggest that overexploitation has not yet significantly purged the allelic variability within the natural populations in protected areas.

Transcriptomic divergence of the Rheum palmatum complex derived from top-geoherb and non-geoherb areas provides the insights into geoherbalism properties of rhubarb

Geoherb usually represents high-quality medicinal herbs with better clinical therapeutic effects, and elucidating the geoherbalism is essential for the quality improvement of traditional Chinese Medicine. However, few researches were conducted to clarify the geoherbalism based on a large scale of transcriptomics. In the present study, we compared the transcriptomes of Rheum palmatum complex derived from top-geoherb and non-geoherb areas to show the geoherbalism properties of rhubarb. A total of 412.32 Gb clean reads were obtained with unigene numbers of 100,615 after assembly. Based on the obtained transcriptome datasets, key enzyme-encoding genes involved in the anthraquinones biosynthesis were also obtained. We also found that 21 anthraquinone-related unigenes were differentially expressed between two different groups, and some of these DEGs were correlated to the content accumulation of five free anthraquinones, indicating that the gene expression profiles may promote the geoherbalism formation of rhubarb. In addition, the selective pressure analyses indicated that most paired orthologous genes between these two groups were subject to negative selection, and only a low proportion of orthologs under positive selection were detected. Functional annotation analyses indicated that these positive-selected genes related to the functions such as gene expression, substance transport, stress response and metabolism, indicating that discrepant environment also enhanced the formation of geoherbalism. Our study not only provided insights for the genetic mechanism of geoherbalism of rhubarb, but also laid more genetic cues for the future rhubarb germplasms improvement and utilization.

Comparative transcriptome analyses of three Gentiana species provides signals for the molecular footprints of selection effects and the phylogenetic relationships

Species in Gentiana section Cruciata are important alpine plants with a center of diversity and speciation in Qinghai-Tibet Plateau (QTP), and some of these species are sympatrically distributed in northeastern QTP. Studies on genome features and natural selection signatures of sympatric species in section Crucata have been impeded by a lack of genomic resources. Here, we showed transcript characterizations and molecular footprints of selection effects on G. straminea, G. dahurica and G. officinalis based on the comparative transcriptome. A total of 62.97 Gb clean reads were obtained with unigene numbers per species ranging from 141,819 to 236,408 after assembly. We found that these three species had similar distribution of functional categories in different databases, and key enzyme-encoding genes involved in the iridoids biosynthesis were also obtained. The selective pressure analyses indicated that most paired orthologs between these three species were subject to negative selection, and only a low proportion of the orthologs that underwent positive selection were detected. We found that some positive selected genes were involved in "catalytic activity", "metabolic process", "response to stimulus" and "response to stress". Besides, large numbers of SSR primer pairs with transferabilities were successfully designed based on the available transcriptome datasets of three Gentiana species. The phylogenetic relationships reconstructed based on 352 single-copy nuclear genes provided a rough phylogenetic framework for this genus and confirmed the monophyly of section Cruciata. Our study not only provides insights for the natural selection effects on sympatric Gentiana species, but also enhances future genetic breeding or evolutionary studies on Qinjiao species.

Occurrence of simple sequence repeats in cDNA sequences of safflower (Carthamus tinctorius) reveals the importance of SSR-containing genes for cell biology and dynamic response to environmental cues

Safflower (Carthamus tinctorius) is a diploid crop plant belonging to the family Asteraceae and is well known as one of important oilseed crops due to edible oil containing unsaturated fatty acids. In recent years it is gaining increased attention for food, pharmaceutical and industrial uses, and hence the updating its breeding methods is necessary. Genic simple sequence repeats (SSRs) in addition of being desire molecular markers, are supposed to influence gene function and the respective phenotype. This study aimed to identify SSRs in cDNA sequences and further analysis of the functional features of the SSR-containing genes to elucidate their role in biological and cellular processes. We identified 1,841 SSR regions in 1,667 cDNA sequences. Among all types of repeats, trinucleotide repeats were the most abundant (35.7%), followed by hexanucleotide (29.6%) and dinucleotide repeats (22.0%). Thirty five SSR primer pairs were validated by PCR reaction, detected a high rate of polymorphism (>57%) among safflower accessions, physically mapped on safflower genome and could clearly discriminate the cultivated accessions from wild relatives. The cDNA-derived SSR markers are suitable for evaluation of genetic diversity, linkage and association mapping studies and genome-based breeding programmes. Occurrence of SSR repeats in biologically-important classes of proteins such as kinases, transferases and transcription factors was inferred from functional analyses, which along with variability of their repeat copies, can endow the cell and whole organism the flexibility of facing with continuously changing environment, and indicate a structure-based evolution mechanism of the genome which acts as an up-to-dating tool for the cell and whole origanism, which is realized in GO terms such as involvement of most SSR-containing genes in biological, cellular and metabolic processes, especially in response to stimulus, response to stress, interaction to other organisms and defense responses.

Comparative Transcriptome Analyses of Different Rheum officinale Tissues Reveal Differentially Expressed Genes Associated with Anthraquinone, Catechin, and Gallic Acid Biosynthesis

Rheum officinale Baill. is an important traditional Chinese medicinal herb, its dried roots and rhizomes being widely utilized to cure diverse diseases. However, previous studies mainly focused on the active compounds and their pharmacological effects, and the molecular mechanism underlying the biosynthesis of these ingredients in R. officinale is still elusive. Here, we performed comparative transcriptome analyses to elucidate the differentially expressed genes (DEGs) in the root, stem, and leaf of R. officinale. A total of 236,031 unigenes with N50 of 769 bp was generated, 136,329 (57.76%) of which were annotated. A total of 5884 DEGs was identified after the comparative analyses of different tissues; 175 and 126 key enzyme genes with tissue-specific expression were found in the anthraquinone, catechin/gallic acid biosynthetic pathway, respectively, and some of these key enzyme genes were verified by qRT-PCR. The phylogeny of the PKS III family in Polygonaceae indicated that probably only PL_741 PKSIII1, PL_11549 PKSIII5, and PL_101745 PKSIII6 encoded PKSIII in the polyketide pathway. These results will shed light on the molecular basis of the tissue-specific accumulation and regulation of secondary metabolites in R. officinale, and lay a foundation for the future genetic diversity, molecular assisted breeding, and germplasm resource improvement of this essential medicinal plant.

Transcriptome profiles of yellowish-white and fuchsia colored flowers in the Rheum palmatum complex reveal genes related to color polymorphism

Flower color variation is ubiquitous in many plant species, and several studies have been conducted to elucidate the underlying molecular mechanism. There are two flower color variants (yellowish-white and fuchsia) in the Rheum palmatum complex, however, few studies have investigated this phenomenon. Here, we used transcriptome sequencing of the two color variants to shed light on the molecular and biochemical basis for these color morphs. Comparison of the two transcriptomes identified 9641 differentially expressed unigenes (DEGs), including 6477 up-regulated and 3163 down-regulated genes. Functional analyses indicated that several DEGs were related to the anthocyanin biosynthesis pathway, and the expression profiles of these DEGs were coincident with the qRT-PCR validation results, indicating that expression levels of structural genes have a profound effect on the color variation in the R. palmatum complex. Our results suggested that the interaction of transcription factors (MYB, bHLH and WRKY) also regulated the anthocyanin biosynthesis in the R. palmatum complex. Estimation of selection pressures using the dN/dS ratio showed that 1106 pairs of orthologous genes have undergone positive selection. Of these positively selected genes, 21 were involved in the anthocyanin biosynthetic pathway, indicating that they may encode the proteins for structural alteration and affect flower color in the R. palmatum complex.

Transcriptome Analysis and Identification of a Female-Specific SSR Marker in Pistacia chinensis Based on Illumina Paired-End RNA Sequencing

Pistacia chinensis Bunge (P. chinensis), a dioecious plant species, has been widely found in China. The female P. chinensis plants are more important than male plants in agricultural production, as their seeds can serve as an ideal feedstock for biodiesel. However, the sex of P. chinensis plants is hard to distinguish during the seedling stage due to the scarcity of available transcriptomic and genomic information. In this work, Illumina paired-end RNA sequencing assay was conducted to unravel the transcriptomic profiles of female and male P. chinensis flower buds. In total, 50,925,088 and 51,470,578 clean reads were obtained from the female and male cDNA libraries, respectively. After quality checks and de novo assembly, a total of 83,370 unigenes with a mean length of 1.3 kb were screened. Overall, 64,539 unigenes (77.48%) could be matched in at least one of the NR, NT, Swiss-Prot, COG, KEGG, and GO databases, 71 of which were putatively related to the floral development of P. chinensis. Additionally, 21,662 simple sequence repeat (SSR) motifs were identified in 17,028 unigenes of P. chinensis, and the mononucleotide motif was the most dominant type of repeats (52.59%) in P. chinensis, followed by dinucleotide (22.29%), trinucleotide (20.15%). The most abundant repeats were AG/CT (13.97%), followed by AAC/GTT (6.75%) and AT/TA (6.10%). Based on these SSR, 983 EST-SSR primers were designed, 151 of which were randomly chosen for validation. Of these validated EST-SSR markers, 25 SSR markers were found to be polymorphic between male and female plants. One SSR marker, namelyPCSSR55, displayed excellent specificity in female plants, which could clearly distinguish between male and female P. chinensis. Altogether, our findings not only reveal that the EST-SSR marker is extremely effective in distinguishing between male and female P. chinensis but also provide a solid framework for sex determination of plant seedlings.

Tissue-specific transcriptome for Rheum tanguticum reveals candidate genes related to the anthraquinones biosynthesis

Rheum tanguticum (Maxim. ex Regel) Maxim. ex Balf. is a herbaceous perennial plant indigenous to China, and its root and rhizomes were usually used as an important traditional Chinese medicine. However, the genomic resources are still scarce for R. tanguticum and even for Rheum genus. Transcriptome datasets from different tissues of R. tanguticum were obtained to screen the genes related to anthraquinones biosynthesis, and five free anthraquinones were also determined. Nine cDNA libraries of roots, stems and leaves were generated, and a total of 272 million high-quality reads were assembled into 257,942 unigenes. Based on the functional annotation, A total of 227 candidate enzyme genes involved in the MVA, MEP, shikimate and polyketide pathways were identified, and several differentially expressed genes found functionally associated with anthraquinones biosynthesis showed distinct tissue-specific expression patterns. Especially, we found that the expression levels of PKS III genes might result in the content differences of free anthraquinones in different tissues of R. tanguticum. Besides, 137,400 SSR loci were identified, and 64,081 SSR primer pairs were successfully designed based on these loci. Our results not only provide cues for the genetic mechanism of anthraquinone content differences in different tissues of R. tanguticum, but also lay genomic foundation for the subsequent genetic engineering and breeding for Rheum species.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01099-8.

Transcriptome analysis of two Pogostemon cablin chemotypes reveals genes related to patchouli alcohol biosynthesis

Pogostemon cablin, a medicinally and economically important perennial herb, is cultivated around the world due to its medicinal and aromatic properties. Different P. cablin cultivars exhibit different morphological traits and patchouli oil components and contents (especially patchouli alcohol (PA) and pogostone (PO)). According to the signature constituent of the leaf, P. cablin was classified into two different chemotypes, including PA-type and PO-type. To better understand the molecular mechanisms of PA biosynthesis, the transcriptomes of Chinese-cultivated P. cablin cv. PA-type “Nanxiang” (NX) and PO-type “Paixiang” (PX) were analyzed and compared with ribonucleic acid sequencing (RNA-Seq) technology. We obtained a total of 36.83 G clean bases from the two chemotypes, compared them with seven databases and revealed 45,394 annotated unigenes. Thirty-six candidate unigenes participating in the biosynthesis of PA were found in the P. cablin transcriptomes. Overall, 8,390 differentially expressed unigenes were identified between the chemotypes, including 2,467 upregulated and 5,923 downregulated unigenes. Furthermore, six and nine differentially expressed genes (DEGs) were mapped to the terpenoid backbone biosynthetic and sesquiterpenoid and triterpenoid biosynthetic pathways, respectively. One key sesquiterpene synthase gene involved in the sesquiterpenoid and triterpenoid biosynthetic pathways, encoding patchoulol synthase variant 1, was significantly upregulated in NX. Additionally, GC-MS analysis of the two chemotypes in this study showed that the content of PA in NX was significantly higher than that of PX, while the content of PO showed the opposite phenotype. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis showed that the DEG expression tendency was consistent with the transcriptome sequencing results. Overall, 23 AP2/ERF, 13 bHLH, 11 MYB, 11 NAC, three Trihelix, 10 WRKY and three bZIP genes that were differentially expressed may act as regulators of terpenoid biosynthesis. Altogether, 8,314 SSRs were recognized within 6,825 unigenes, with a distribution frequency of 18.32%, among which 1,202 unigenes contained more than one SSR. The transcriptomic characteristics of the two P. cablin chemotypes are comprehensively reported in this study, and these results will contribute to a better understanding of the molecular mechanism of PA biosynthesis. Our transcriptome data also provide a valuable genetic resource for further studies on P. cablin.

Detection and application of novel SSR markers from transcriptome data for Astronium fraxinifolium Schott, a threatened Brazilian tree species

Astronium fraxinifolium is an endangered tree species from Brazil. Due to its significance in environmental reforestation, as well as the continued exploitation of its wood, it is necessary to develop management programs that support the conservation of the species. Simple sequence repeats (SSR) or microsatellite markers are widely used in population genetic studies across a range of diverse organisms. In this study, we present the first SSR markers developed for A. fraxinifolium as well as their frequency and distribution based on transcriptome data. From transcriptome data, we identified more than 100 thousand sequences presenting microsatellites, with a predominant distribution of trinucleotide repeats. From the initial screening, we selected 20 microsatellite loci which were validated and evaluated for genetic indices in two natural populations. All loci were polymorphic, ranging from four to 11 alleles per locus. The observed and expected heterozygosities ranged from 0 to 1.0 and from 0.533 to 1.0, respectively, while the genetic differentiation (G_ST = 0.363) was greater within than between populations. The developed SSR loci from RNA-Seq data provides a foundation for future studies on genetic diversity and population structure, mating system, and gene flow for A. fraxinifolium populations and related species, aiming at conservation and management.

Discovery of novel genic-SSR markers from transcriptome dataset of an important non-human primate, Macaca fascicularis

Macaca fascicularis, also known as the cynomolgus macaque, is an important non-human primate animal model used in biomedical research. It is an Old-World primate widely distributed in Southeast Asia and is one of the most abundant macaque species in Malaysia. However, the genetic structure of wild cynomolgus macaque populations in Malaysia has not been thoroughly elucidated. In this study, we developed genic-simple sequence repeat (genic-SSR) markers from an in-house transcriptome dataset generated from the Malaysian cynomolgus macaque via RNA sequencing, and applied these markers on 26 cynomolgus macaque individuals. A collection of 14,751 genic-SSRs were identified, where 13,709 were perfect SSRs. Dinucleotide repeats were the most common repeat motifs with a frequency of 65.05%, followed by trinucleotide repeats (20.55%). Subsequently, we designed 300 pairs of primers based on perfect di- and trinucleotide SSRs, in which 105 SSRs were associated with functional genes. A subset of 30 SSR markers were randomly selected and validated, yielding 19 polymorphic markers with an average polymorphism information content value of 0.431. The development of genic-SSR markers in this study is indeed timely to provide useful markers for functional and population genetic studies of the cynomolgus macaque and other related non-human primate species.

De novo assembly of transcriptomes, mining, and development of novel EST-SSR markers in Curcuma alismatifolia (Zingiberaceae family) through Illumina sequencing

Curcuma alismatifolia widely used as an ornamental plant in Thailand and Cambodia. This species of herbaceous perennial from the Zingiberaceae family, includes cultivars with a wide range of colours and long postharvest life, and is used as an ornamental cut flower, as a potted plant, and in exterior landscapes. For further genetic improvement, however, little genomic information and no specific molecular markers are available. The present study used Illumina sequencing and de novo transcriptome assembly of two C. alismatifolia cvs, 'Chiang Mai Pink' and 'UB Snow 701', to develop simple sequence repeat markers for genetic diversity studies. After de novo assembly, 62,105 unigenes were generated and 48,813 (78.60%) showed significant similarities versus six functional protein databases. In addition, 9,351 expressed sequence tag-simple sequence repeats (EST-SSRs) were identified with a distribution frequency of 12.5% total unigenes. Out of 8,955 designed EST-SSR primers, 150 primers were selected for the development of potential molecular markers. Among these markers, 17 EST-SSR markers presented a moderate level of genetic diversity among three C. alismatifolia cultivars, one hybrid, three Curcuma, and two Zingiber species. Three different genetic groups within these species were revealed using EST-SSR markers, indicating that the markers developed in this study can be effectively applied to the population genetic analysis of Curcuma and Zingiber species. This report describes the first analysis of transcriptome data of an important ornamental ginger cultivars, also provides a valuable resource for gene discovery and marker development in the genus Curcuma.

De Novo Transcriptome Assembly and Population Genetic Analyses for an Endangered Chinese Endemic Acer miaotaiense (Aceraceae)

Acer miaotaiense (P. C. Tsoong) is a rare and highly endangered plant in China. Because of the lack of genomic information and the limited number of available molecular markers, there are insufficient tools to determine the genetic diversity of this species. Here, 93,305 unigenes were obtained by multiple assembled contigs with a transcriptome sequencing program. Furthermore, 12,819 expressed sequence tag-derived simple sequence repeat (EST-SSR) markers were generated, 300 were randomly selected and synthesized, 19 primer pairs were identified as highly polymorphic (average number of alleles (N_a) = 8, expected heterozygosity (H_e) = 0.635, polymorphism information content (PIC) = 0.604) and were further used for population genetic analysis. All 261 samples were grouped into two genetic clusters by UPGMA, a principal component analyses and a STRUCTURE analyses. A moderate level of genetic differentiation (genetic differentiation index (F_st) = 0.059–0.116, gene flow = 1.904–3.993) among the populations and the major genetic variance (81.01%) within populations were revealed by the AMOVA. Based on the results, scientific conservation strategies should be established using in situ and ex situ conservation strategies. The study provides useful genetic information for the protection of precious wild resources and for further research on the origin and evolution of this endangered plant and its related species.

Mining and Development of Novel SSR Markers Using Next Generation Sequencing (NGS) Data in Plants

Microsatellites, or simple sequence repeats (SSRs), are one of the most informative and multi-purpose genetic markers exploited in plant functional genomics. However, the discovery of SSRs and development using traditional methods are laborious, time-consuming, and costly. Recently, the availability of high-throughput sequencing technologies has enabled researchers to identify a substantial number of microsatellites at less cost and effort than traditional approaches. Illumina is a noteworthy transcriptome sequencing technology that is currently used in SSR marker development. Although 454 pyrosequencing datasets can be used for SSR development, this type of sequencing is no longer supported. This review aims to present an overview of the next generation sequencing, with a focus on the efficient use of de novo transcriptome sequencing (RNA-Seq) and related tools for mining and development of microsatellites in plants.

Development of microsatellite markers for Myracrodruon urundeuva (F.F. & M.F. Allemão), a highly endangered species from tropical forest based on next-generation sequencing

Myracrodruon urundeuva is a tree species of high economic importance due the strength and durability of its wood. Threatened of extinction in Brazil, it is present only in a few forest remnants, mostly in conservation units. Currently, there is little information on the genetic diversity of natural populations in Brazil and even less information about the genome of this species. Here, new species-specific microsatellite loci were developed based on next-generation sequencing (Illumina). More than 100,000 loci were identified in the run, with di- to hexanucleotides motifs. Of these, 20 loci were selected for validation in 30 individuals, with 15 successfully polymorphic loci detected. The number of alleles ranged among loci from 3 to 16, with an average of 7.73, expected (H _e ) and observed (H _o ) heterozygosity ranged from 0.246 to 0.902 and from 0.103 to 0.867, respectively. These results point out that these new set of markers has a great potential for use in population genetic studies for genetic conservation of the species.

De novo transcriptome assembly of Zanthoxylum bungeanum using Illumina sequencing for evolutionary analysis and simple sequence repeat marker development

Zanthoxylum, an ancient economic crop in Asia, has a satisfying aromatic taste and immense medicinal values. A lack of genomic information and genetic markers has limited the evolutionary analysis and genetic improvement of Zanthoxylum species and their close relatives. To better understand the evolution, domestication, and divergence of Zanthoxylum, we present a de novo transcriptome analysis of an elite cultivar of Z. bungeanum using Illumina sequencing; we then developed simple sequence repeat markers for identification of Zanthoxylum. In total, we predicted 45,057 unigenes and 22,212 protein coding sequences, approximately 90% of which showed significant similarities to known proteins in databases. Phylogenetic analysis indicated that Zanthoxylum is relatively recent and estimated to have diverged from Citrus ca. 36.5–37.7 million years ago. We also detected a whole-genome duplication event in Zanthoxylum that occurred 14 million years ago. We found no protein coding sequences that were significantly under positive selection by Ka/Ks. Simple sequence repeat analysis divided 31 Zanthoxylum cultivars and landraces into three major groups. This Zanthoxylum reference transcriptome provides crucial information for the evolutionary study of the Zanthoxylum genus and the Rutaceae family, and facilitates the establishment of more effective Zanthoxylum breeding programs.

De novo transcriptomic analysis and development of EST-SSRs for Sorbus pohuashanensis (Hance) Hedl

Sorbus pohuashanensis is a native tree species of northern China that is used for a variety of ecological purposes. The species is often grown as an ornamental landscape tree because of its beautiful form, silver flowers in early summer, attractive pinnate leaves in summer, and red leaves and fruits in autumn. However, development and further utilization of the species are hindered by the lack of comprehensive genetic information, which impedes research into its genetics and molecular biology. Recent advances in de novo transcriptome sequencing (RNA-seq) technology have provided an effective means to obtain genomic information from non-model species. Here, we applied RNA-seq for sequencing S. pohuashanensis leaves and obtained a total of 137,506 clean reads. After assembly, 96,213 unigenes with an average length of 770 bp were obtained. We found that 64.5% of the unigenes could be annotated using bioinformatics tools to analyze gene function and alignment with the NCBI database. Overall, 59,089 unigenes were annotated using the Nr database(non-redundant protein database), 35,225 unigenes were annotated using the GO (Gene Ontology categories) database, and 33,168 unigenes were annotated using COG (Cluster of Orthologous Groups). Analysis of the unigenes using the KEGG (Kyoto Encyclopedia of Genes and Genomes) database indicated that 13,953 unigenes were involved in 322 metabolic pathways. Finally, simple sequence repeat (SSR) site detection identified 6,604 unigenes that included EST-SSRs and a total of 7,473 EST-SSRs in the unigene sequences. Fifteen polymorphic SSRs were screened and found to be of use for future genetic research. These unigene sequences will provide important genetic resources for genetic improvement and investigation of biochemical processes in S. pohuashanensis.

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.