De novo transcriptomic analysis and identification of EST-SSR markers in Stephanandra incisa

Population genetics analysis of Diospyrosmun A.Chev. ex Lecomte (Ebenaceae) based on EST-SSR markers derived from a novel transcriptome

Diospyrosmun A.Chev. ex Lecomte (Ebenaceae), a native evergreen tree in Vietnam, has important economic and ecological values. The absence of effective and reliable molecular markers has hampered the study of D.mun's genetic diversity and population structure, even though it is an endemic and endangered species. Therefore, significant enrichment of genomic resources is urgently needed to uncover and better understand the genetic architecture of D.mun. This study aims to demonstrate an efficient and reliable tool to explore the polymorphism within D.mun germplasm. It provides a valuable platform for the breeding and conservation of this species and other endangered species worldwide. The Illumina HiSeq™ 4000 sequencing technology was applied for the transcriptomic analysis, genetic differentiation and population structure of D.mun in Vietnam. In this study, the transcriptomes of D.mun were analysed using the Illumina HiSeqTM 4000 sequencing system and a total of 5,588,615,700 base pairs were generated. De novo assembly indicated that 91,134 unigenes were generated (average length = 645.55 bp, N50 = 957 bp, Q20 = 98.08% and Q30 = 94.51%). A total of 92,798 and 21,134 unigenes had significant similarities amongst Nr and Swiss-Prot, respectively. In the GO database, 19,929 unigenes were annotated and these genes were divided into three major categories and 50 subcategories. In the KOG analysis, 18,499 unigenes were annotated and divided into 25 gene function categories. In the KEGG analysis, 12,017 unigenes were annotated. According to the related pathways involved, they could be classified into 56 subclasses. In this study, we have identified a total of 9,391 EST-SSR markers. Ten microsatellite loci were employed to assess the genetic diversity and structure of 82 adult D.mun trees across three populations in Vietnam. The results indicated moderate levels of genetic diversity with PIC = 0.77, NA = 3.9, NE = 2.8, Ho = 0.56 and HE = 0.58 and the fixation index value was recorded as positive for three populations (NS, NH and CP). Genetic differentiation among populations was low (FST = 0.045), suggesting limited gene flow (Nm = 5.34). This result indicates gene exchange between the populations of ancient D.mun from different geographical areas and regions. The analysis of molecular variance (AMOVA) showed that high genetic variation existed within individuals (91%) compared to amongst populations (4%). Genetic structure analysis, DAPC and the NJ tree indicated that the three populations were divided into three main clusters. With this study, we provide a molecular resoureces for the breeding and conservation of D.mun.

Genetic diversity, population structure, and DNA fingerprinting of Ailanthus altissima var. erythrocarpa based on EST-SSR markers

Ailanthus altissima var. erythrocarpa is an A. altissima variety with high economic, ecological and ornamental value, but there have been no reports on the development of SSR primers for it. According to the SSR primer information provided by the transcriptome of A. altissima var. erythrocarpa, 120 individuals with different redness levels were used to screen polymorphic primers. Transcriptomic analysis revealed 10,681 SSR loci, of which mononucleotide repeats were dominant (58.3%), followed by dinucleotide and trinucleotide repeats (16.6%, 15.1%) and pentanucleotide repeats (0.2%). Among 140 pairs of randomly selected primers, nineteen pairs of core primers with high polymorphism were obtained. The average number of alleles (Na), average number of effective alleles (Ne), average Shannon's diversity index (I), average observed heterozygosity (Ho), average expected heterozygosity (He), fixation index (F) and polymorphic information content (PIC) were 11.623, 4.098, 1.626, 0.516, 0.696, 0.232 and 0.671, respectively. Nineteen EST-SSR markers were used to study the genetic diversity and population structure of A. altissima var. erythrocarpa. The phylogenetic tree, PCoA, and structure analysis all divided the tested resources into two categories, clearly showing the genetic variation between individuals. The population showed high genetic diversity, mainly derived from intraspecific variation. Among nineteen pairs of primers, 4 pairs (p33, p15, p46, p92) could effectively distinguish and be used for fingerprinting of the tested materials. This study is of great significance for genetic diversity analysis and molecular-assisted breeding of A. altissima var. erythrocarpa.

Development of unigene-derived SSR markers from RNA-seq data of Uraria lagopodioides (Fabaceae) and their application in the genus Uraria Desv. (Fabaceae)

BACKGROUND

Uraria Desv. belongs to the tribe Desmodieae (Fabaceae), a group of legume plants, some of which have medicinal properties. However, due to a lack of genomic information, the interspecific relationships, genetic diversity, population genetics, and identification of functional genes within Uraria species are still unclear.

RESULTS

Using RNA-Seq, a total of 66,026 Uraria lagopodioides unigenes with a total sequence content of 52,171,904 bp were obtained via de novo assembly and annotated using GO, KEGG, and KOG databases. 17,740 SSRs were identified from a set of 66,026 unigenes. Cross-species amplification showed that 54 out of 150 potential unigene-derived SSRs were transferable in Uraria, of which 19 polymorphic SSRs were developed. Cluster analysis based on polymorphisms successfully distinguished seven Uraria species and revealed their interspecific relationships. Seventeen samples of seven Uraria species were clustered into two monophyletic clades, and phylogenetic relationships of Uraria species based on unigene-derived SSRs were consistent with classifications based on morphological characteristics.

CONCLUSIONS

Unigenes annotated in the present study will provide new insights into the functional genomics of Uraria species. Meanwhile, the unigene-derived SSR markers developed here will be invaluable for assessing the genetic diversity and evolutionary history of Uraria and relatives.

Analysis of the genetic diversity and population structure of Monochasma savatieri Franch. ex Maxim using novel EST-SSR markers

Background

Monochasma savatieri Franch. ex Maxim is a medicinally valuable herb. However, the collection and protection of the wild germplasm resources of M. savatieri are still insufficient, and their genetic diversity and population structure have been poorly studied.

Results

We collected and examined 46 M. savatieri individuals from Fujian, Hunan, Jiangxi, and Zhejiang provinces for genetic diversity and population structure, using 33 newly developed expressed sequence tag-simple sequence repeat (EST-SSR) markers. Applying these markers, we detected a total of 208 alleles, with an average of 6.303 alleles per locus. The polymorphic information content varied from 0.138 to 0.884 (average: 0.668), indicating a high level of polymorphism. At the population level, there was a low degree of genetic diversity among populations (I = 0.535, He = 0.342), with Zhejiang individuals showing the highest genetic diversity among the four populations (Fst = 0.497), which indicated little gene flow within the M. savatieri populations (Nm = 0.253). Mantel test analysis revealed a significant positive correlation between geographical and genetic distance among populations (R² = 0.3304, p < 0.05), and structure and principal coordinate analyses supported classification of populations into three clusters, which was consistent with the findings of cluster analysis.

Conclusions

As a rare medicinal plants, the protection of M. savatieri does not look optimistic, and accordingly, protective efforts should be beefed up on the natural wild populations. This study provided novel tools and insights for designing effective collection and conservation strategies for M. savatieri.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08832-x.

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.