Development of SSR markers for a Tibetan medicinal plant, Lancea tibetica (Phrymaceae), based on RAD sequencing

SSR Marker Acquisition and Application from Transcriptome of Captive Chinese Forest Musk Deer (Moschus berezovskii)

Forest musk deer (Moschus berezovskii) is one of the most endangered medicinally important wild animals in the world. Forest musk deer farming is the main way of production of musk. However, the single provenance and lack of genetic information lead to reduced genetic diversity of forest musk deer. Therefore, more SSR markers need to be developed to identify forest musk deer germplasm. In this study, bone marrow derived mesenchymal cells were used to construct cDNA library for transcriptome sequencing. The datasets were de novo assembled and annotated. 9 polymorphic simple sequence repeat (SSR) markers were finally identified and used to detect population genetic diversity. 6.07 Gb clean data were generated using Illumina sequencing technology, and de novo assembled into 138,591 transcripts and 81,553 unigenes. 5,777 simple sequence repeats (SSRs) were identified, in which there were 578 repeating motif types, with mono-nucleotide and tri-nucleotides comprising 55.88% and 25.60%, respectively. 100 primer pairs were designed to validate amplification and polymorphism using DNA from fecal samples. 9 polymorphic SSRs were developed and used to detect population genetic diversity of 122 forest musk deer in 2 farms. The average number of alleles per locus varied from 4 to 15 (average = 8.3). The observed heterozygosity (HO) per locus ranged from 0.102 to 0.941, while the expected heterozygosity (HE) per locus was from 0.111 to 0.651. All loci deviated significantly from the Hardy-Weinberg equilibrium (p < 0.001). The polymorphism information content (PIC) of these loci varied from 0.108 to 0.619. 9 polymorphic SSR markers were developed in this research. These sites can be used for breeding planning and conservation of germplasm resources.

The genetic legacy of fragmentation and overexploitation in the threatened medicinal African pepper-bark tree, Warburgia salutaris

The pepper-bark tree (Warburgia salutaris) is one of the most highly valued medicinal plant species worldwide. Native to southern Africa, this species has been extensively harvested for the bark, which is widely used in traditional health practices. Illegal harvesting coupled with habitat degradation has contributed to fragmentation of populations and a severe decline in its distribution. Even though the species is included in the IUCN Red List as Endangered, genetic data that would help conservation efforts and future re-introductions are absent. We therefore developed new molecular markers to understand patterns of genetic diversity, structure, and gene flow of W. salutaris in one of its most important areas of occurrence (Mozambique). In this study, we have shown that, despite fragmentation and overexploitation, this species maintains a relatively high level of genetic diversity supporting the existence of random mating. Two genetic groups were found corresponding to the northern and southern locations. Our study suggests that, if local extinctions occurred in Mozambique, the pepper-bark tree persisted in sufficient numbers to retain a large proportion of genetic diversity. Management plans should concentrate on maintaining this high level of genetic variability through both in and ex-situ conservation actions.

Development of novel microsatellite markers using RAD sequencing technology for diversity assessment of rambutan (Nephelium lappaceum L.) germplasm

The trend of microsatellite marker discovery and development revolved as a result of the advancement of next generation sequencing (NGS) technology as it has developed numerous microsatellites within a short period of time at a low cost. This study generated microsatellite markers using RAD sequencing technologies for the understudied Nephelium lappaceum. A total of 1403 microsatellite markers were successfully designed, which consisted of 853 di-, 525 tri-, 17 tetra-, 5 penta-, and 3 hexanucleotide microsatellite markers. Subsequently, selection of 39 microsatellites was made for the evaluation of genetic diversity of the selected 22 rambutan varieties. Twelve microsatellites, which exhibited high call rates across the samples, were used to assess the diversity of the aforementioned rambutan varieties. The analysis of 12 microsatellites revealed the presence of 72 alleles and six alleles per locus in average. Furthermore, the polymorphic information content (PIC) value ranged from 0.326 (NlaSSR20) to 0.832 (NlaSSR32), which included an average of 0.629 per locus, while the generated Neighbour Joining dendrogram showed two major clusters. The pairwise genetic distance of shared alleles exhibited a range of values from 0.046 (R134↔R170) to 0.818 (R5↔R170), which suggested highest dissimilarity detected between R5 and R170. Notably, these research findings would useful for varietal identification, proper management and conservation of the genetic resources, and exploitation and utilization in future breeding programs.

Genetic Structure and Eco-Geographical Differentiation of Lancea tibetica in the Qinghai-Tibetan Plateau

The uplift of the Qinghai-Tibetan Plateau (QTP) had a profound impact on the plant speciation rate and genetic diversity. High genetic diversity ensures that species can survive and adapt in the face of geographical and environmental changes. The Tanggula Mountains, located in the central of the QTP, have unique geographical significance. The aim of this study was to investigate the effect of the Tanggula Mountains as a geographical barrier on plant genetic diversity and structure by using Lancea tibetica. A total of 456 individuals from 31 populations were analyzed using eight pairs of microsatellite makers. The total number of alleles was 55 and the number per locus ranged from 3 to 11 with an average of 6.875. The polymorphism information content (PIC) values ranged from 0.2693 to 0.7761 with an average of 0.4378 indicating that the eight microsatellite makers were efficient for distinguishing genotypes. Furthermore, the observed heterozygosity (H_o), the expected heterozygosity (H_e), and the Shannon information index (I) were 0.5277, 0.4949, and 0.9394, respectively, which indicated a high level of genetic diversity. We detected high genetic differentiation among all sampling sites and restricted gene flow among populations. Bayesian-based cluster analysis (STRUCTURE), principal coordinates analysis (PCoA), and Neighbor-Joining (NJ) cluster analysis based on microsatellite markers grouped the populations into two clusters: the southern branch and the northern branch. The analysis also detected genetic barriers and restricted gene flow between the two groups separated by the Tanggula Mountains. This study indicates that the geographical isolation of the Tanggula Mountains restricted the genetic connection and the distinct niches on the two sides of the mountains increased the intraspecific divergence of the plants.

Traditional Tibetan medicinal plants: a highlighted resource for novel therapeutic compounds

Around 70-80% of drugs used in traditional Tibetan medicine (TTM) come from Qinghai Tibet Plateau, the majority of which are plants. The biological and medicinal culture diversity on Qinghai Tibet Plateau are amazing and constitute a less tapped resource for innovative drug research and development. Meanwhile, the problem of the exhausting Tibetan medicine resources is worrying. Here, the latest awareness, as well as the gaps of the traditional Tibetan medicinal plant issues in drug development and clinical usage of TTM compounds, was systematically reviewed and highlighted. The TTM resource studies should be enhanced within the context of deeper and more extensive investigations of molecular biology and genomics of TTM plants, phytometabolites and metabolomics and ethnopharmacology-based bioactivity, thus enabling the sustainable conservation and exploitation of Tibetan medicinal resource.

Development of microsatellite markers and analysis of genetic diversity of Barbatia virescens in the southern coasts of China

BACKGROUND

The blood clam Barbatia virescens is an ecologically and economically important species in the southern coast of China. Understanding of the genetic structure of B. virescens populations is vital to breeding strategies and conservation programs.

OBJECTIVE

To develop and characterize a set of microsatellites loci primers for B. virescens, and provide helpful information for reasonable utilization and protection of B. virescens natural resources.

METHODS

The microsatellites of B. virescens were detected using a RAD-seq approach based on an Illumina sequencing platform. For the test of microsatellite development, we calculated the number of alleles (Na), observed heterozygosities (Ho), expected heterozygosities (He) and exact tests for deviations from Hardy-Weinberg equilibrium (HWE). Twelve polymorphic loci were used to access the genetic diversity and population structure of four B. virescens populations.

RESULTS

In this study, 50,729 microsatellites of B. virescens were detected. Twenty-two polymorphic microsatellite loci were developed for B. virescens. The number of alleles per locus ranged from 6 to 15, and expected heterozygosities varied from 0. 567 to 0.911. All the PIC values of the 22 loci were greater than 0.5, indicating that these markers were highly informative for further genetic analysis. Twelve loci were selected to analyze genetic diversity and population structure of four B. virescens populations collected from different geographical regions along the southern coast of China. The results showed moderate to high levels of genetic diversity in the four populations (mean Ar = 7.756-8.133; mean Ho = 0.575-0.639; mean He = 0.754-0.775). Pairwise F_ST estimates indicated that there was significant divergence among the four populations.

CONCLUSION

This study not only provides a large scale of sequence information of microsatellites which are valuable for future genetic mapping, trait association and kinship among B. virescens, but also offers useful information for the sustainable management of natural stocks and the development of breeding industry of B. virescens.

Deep Intraspecific Divergence in the Endemic Herb Lancea tibetica (Mazaceae) Distributed Over the Qinghai-Tibetan Plateau

Qinghai-Tibetan Plateau (QTP) is an important biodiversity hub, which is very sensitive to climate change. Here in this study, we investigated genetic diversity and past population dynamics of Lancea tibetica (Mazaceae), an endemic herb to QTP and adjacent highlands. We sequenced chloroplast and nuclear ribosomal DNA fragments for 429 individuals, collected from 29 localities, covering their major distribution range at the QTP. A total of 19 chloroplast haplotypes and 13 nuclear genotypes in two well-differentiated lineages, corresponding to populations into two groups isolated by Tanggula and Bayangela Mountains. Meanwhile, significant phylogeographical structure was detected among sampling range of L. tibetica, and 61.50% of genetic variations was partitioned between groups. Gene flow across the whole region appears to be restricted by high mountains, suggesting a significant role of geography in the genetic differences between the two groups. Divergence time between the two lineages dated to 8.63 million years ago, which corresponded to the uplifting of QTP during the late Miocene and Pliocene. Ecological differences were found between both the lineages represent species-specific characteristics, sufficient to keep the lineages separated to a high degree. The simulated distribution from the last interglacial period to the current period showed that the distribution of L. tibetica experienced shrinkage and expansion. Climate changes during the Pleistocene glacial-interglacial cycles had a dramatic effect on L. tibetica distribution ranges. Multiple refugia of L. tibetica might have remained during the species history, to south of the Tanggula and north of Bayangela Mountains, both appeared as topological barrier and contributed to restricting gene flow between the two lineages. Together, geographic isolation and climatic factors have played a fundamental role in promoting diversification and evolution of L. tibetica.

Trends in herbgenomics

From Shen Nong's Herbal Classic (Shennong Bencao Jing) to the Compendium of Materia Medica (Bencao Gangmu) and the first scientific Nobel Prize for the mainland of China, each milestone in the historical process of the development of traditional Chinese medicine (TCM) involves screening, testing and integrating. After thousands of years of inheritance and development, herbgenomics (bencaogenomics) has bridged the gap between TCM and international advanced omics studies, promoting the application of frontier technologies in TCM. It is a discipline that uncovers the genetic information and regulatory networks of herbs to clarify their molecular mechanism in the prevention and treatment of human diseases. The main theoretical system includes genomics, functional genomics, proteomics, transcriptomics, metabolomics, epigenomics, metagenomics, synthetic biology, pharmacogenomics of TCM, and bioinformatics, among other fields. Herbgenomics is mainly applicable to the study of medicinal model plants, genomic-assisted breeding, herbal synthetic biology, protection and utilization of gene resources, TCM quality evaluation and control, and TCM drug development. Such studies will accelerate the application of cutting-edge technologies, revitalize herbal research, and strongly promote the development and modernization of TCM.

Genetic Evaluation of Natural Populations of the Endangered Conifer Thuja koraiensis Using Microsatellite Markers by Restriction-Associated DNA Sequencing

Thuja koraiensis Nakai is an endangered conifer of high economic and ecological value in Jilin Province, China. However, studies on its population structure and conservation genetics have been limited by the lack of genomic data. Here, 37,761 microsatellites (simple sequence repeat, SSR) were detected based on 875,792 de novo-assembled contigs using a restriction-associated DNA (RAD) approach. Among these SSRs, 300 were randomly selected to test for polymorphisms and 96 obtained loci were able to amplify a fragment of expected size. Twelve polymorphic SSR markers were developed to analyze the genetic diversity and population structure of three natural populations. High genetic diversity (mean NA = 5.481, HE = 0.548) and moderate population differentiation (pairwise Fst = 0.048–0.078, Nm = 2.940–4.958) were found in this species. Molecular variance analysis suggested that most of the variation (83%) existed within populations. Combining the results of STRUCTURE, principal coordinate, and neighbor-joining analysis, the 232 individuals were divided into three genetic clusters that generally correlated with their geographical distributions. Finally, appropriate conservation strategies were proposed to protect this species. This study provides genetic information for the natural resource conservation and utilization of T. koraiensis and will facilitate further studies of the evolution and phylogeography of the species.