Development and Characterization of High-Throughput EST-Based SSR Markers for Pogostemon cablin Using Transcriptome Sequencing

Population genetic variation and geographic distribution of suitable areas of Coptis species in China

Introduction

The rhizomes of Coptis plants have been used in traditional Chinese medicine over 2000 years. Due to increasing market demand, the overexploitation of wild populations, habitat degradation and indiscriminate artificial cultivation of Coptis species have severely damaged the native germplasms of species in China.

Methods

Genome-wide simple-sequence repeat (SSR) markers were developed using the genomic data of C. chinensis. Population genetic diversity and structure of 345 Coptis accessions collected from 32 different populations were performed based on these SSRs. The distribution of suitable areas for three taxa in China was predicted and the effects of environmental variables on genetic diversity in relation to different population distributions were further analyzed.

Results

22 primer pairs were selected as clear, stable, and polymorphic SSR markers. These had an average of 16.41 alleles and an average polymorphism information content (PIC) value of 0.664. In the neighbor-joining (N-J) clustering analysis, the 345 individuals clustered into three groups, with C. chinensis, C. chinensis var. brevisepala and C. teeta being clearly separated. All C. chinensis accessions were further divided into four subgroups in the population structure analysis. The predicted distributions of suitable areas and the environmental variables shaping these distributions varied considerably among the three species.

Discussion

Overall, the amount of solar radiation, precipitation and altitude were the most important environmental variables influencing the distribution and genetic variation of three species. The findings will provide key information to guide the conservation of genetic resources and construction of a core reserve for species.

Fine-Scale analysis of both wild and cultivated horned galls provides insight into their quality differentiation

BACKGROUND

Galla chinensis is a traditional Chinese medicine (TCM) produced due to the interaction between the Fordinae aphids and the Rhus plant species. Horned galls with high tannin content are the most widely cultivated gall type, and Wufeng county of Hubei province in China is the center of cultivation. However, long-term artificial cultivation and domestication of horned galls to meet the increasing production demand have led to quality degradation. Understanding the reasons underlying quality degradation is urgent for horned gall production and application. The present study used a combination of metabolic, genetic, and ecological analyses to investigate the quality and genetic differentiation of the horned galls under long-term domestication as well as the potential relationships between them.

RESULTS

Analysis of gallic acid content and other three phenotypic traits (fresh weight, gall size, and wall thickness) revealed quality differentiation of horned galls collected from five locations in Wufeng, in which the cultivated samples from Wang Jiaping (WJP) showed the highest degradation. Genetic differentiation between the cultivated and wild Rhus chinensis trees in WJP, and between WJP and the other populations was detected based on SSR molecular markers, however, no significant difference in genetic structure was seen for the aphid populations. Among the various ecological factors examined, temperature was identified as the primary one affecting the quality of horned galls.

CONCLUSIONS

Both genetic and ecological factors caused quality differentiation of horned galls. The collection of diverse germplasm of host trees and aphids will help reduce the quality degradation of horned galls in Wufeng.

Chromosome-level and haplotype-resolved genome provides insight into the tetraploid hybrid origin of patchouli

Patchouli (Pogostemon cablin (Blanco) Benth.), a member of the Lamiaceae family, is an important aromatic plant that has been widely used in medicine and perfumery. Here, we report a 1.94 Gb chromosome-scale assembly of the patchouli genome (contig N50 = 7.97 Mb). The gene annotation reveals that tandem duplication of sesquiterpene biosynthetic genes may be a major contributor to the biosynthesis of patchouli bioactivity components. We further phase the genome into two distinct subgenomes (A and B), and identify a chromosome substitution event that have occurred between them. Further investigations show that a burst of universal LTR-RTs in the A subgenome lead to the divergence between two subgenomes. However, no significant subgenome dominance is detected. Finally, we track the evolutionary scenario of patchouli including whole genome tetraploidization, subgenome divergency, hybridization, and chromosome substitution, which are the key forces to determine the complexity of patchouli genome. Our work sheds light on the evolutionary history of patchouli and offers unprecedented genomic resources for fundamental patchouli research and elite germplasm development.

The ploidy level of patchouli, an aromatic plant in the Lamiaceae family, remain unclear. Here, the authors assemble a chromosome-level and haplotype-resolved genome for patchouli and reveal that it is tetraploid hybrid as well as compensated aneuploidy.

Full-Length SMRT Transcriptome Sequencing and SSR Analysis of Bactrocera dorsalis (Hendel)

Simple Summary

In this study, a full-length transcriptome was analyzed with single-molecule real-time (SMRT) sequencing, which was first used to discover simple sequence repeat (SSR) genetic markers from B. dorsalis. Moreover, SSR markers from isoforms were screened for the identification of species diversity. These results could provide molecular biology methods for further population research.

Abstract

Bactrocera dorsalis (Hendel), as one of the most notorious and destructive invasive agricultural pests in the world, causes damage to over 250 different types of fruits and vegetables throughout tropical and subtropical areas. PacBio single-molecule real-time (SMRT) sequencing was used to generate the full-length transcriptome data of B. dorsalis. A total of 40,319,890 subreads (76.6 Gb, clean reads) were generated, including 535,241 circular consensus sequences (CCSs) and 386,916 full-length non-concatemer reads (FLNCs). Transcript cluster analysis of the FLNC reads revealed 22,780 high-quality reads (HQs). In total, 12,274 transcripts were functionally annotated based on four different databases. A total of 1978 SSR loci were distributed throughout 1714 HQ transcripts, of which 1926 were complete SSRs and 52 were complex SSRs. Among the total SSR loci, 2–3 nucleotide repeats were dominant, occupying 83.62%, of which di- and tri- nucleotide repeats were 39.38% and 44.24%, respectively. We detected 105 repeat motifs, of which AT/AT (50.19%), AC/GT (39.15%), CAA/TTG (32.46%), and ACA/TGT (10.86%) were the most common in di- and tri-nucleotide repeats. The repeat SSR motifs were 12–190 bp in length, and 1638 (88.02%) were shorter than 20 bp. According to the randomly selected microsatellite sequence, 80 pairs of primers were designed, and 174 individuals were randomly amplified by PCR using primers. The number of primers that had amplification products with clear bands and showed good polymorphism came to 41, indicating that this was a feasible way to explore SSR markers from the transcriptomic data of B. dorsalis. These results lay a foundation for developing highly polymorphic microsatellites for researching the functional genomics, population genetic structure, and genetic diversity of B. dorsalis.

Toward understanding of the methoxylated flavonoid biosynthesis pathway in Dracocephalum kotschyi Boiss

Nowadays, with the development and advancement of next-generation sequencing technologies, a new path has been provided for transcriptomic studies. In this study, the transcriptome of Dracocephalum kotschyi Boiss., as an endemic and endangered plant which is contained a large amount of valuable secondary metabolites with antioxidant and anticancer properties, was sequenced. Then functional annotation and gene ontology analysis for 165,597 assembled transcripts were performed, most were associated with the metabolic pathways. This might be because there are various active biochemical pathways in this plant. Furthermore, after comprehensive transcript annotation, the putative genes involved in the main metabolic pathways of D. kotschyi were identified. Then, the biosynthetic pathway of its valuable methoxylated flavones was proposed. Finally, the accumulations of important methoxylated-flavone metabolites in three different tissues were quantified by HPLC. The relative expression of the genes involved in the proposed pathway was investigated by qRT-PCR, which indicated high expression levels in the bud tissue. The present results may lead to the design strategies to preserve the genetic diversity of endangered D. kotschyi plants and apply the new methods for engineering its valuable methoxylated-flavones pathway.

Development of GBTS and KASP Panels for Genetic Diversity, Population Structure, and Fingerprinting of a Large Collection of Broccoli (Brassica oleracea L. var. italica) in China

Broccoli (Brassica oleracea var. italica) is one of the most important and nutritious vegetables widely cultivated in China. In the recent four decades, several improved varieties were bred and developed by Chinese breeders. However, the efforts for improvement of broccoli are hindered by limited information of genetic diversity and genetic relatedness contained within the available germplasms. This study evaluated the genetic diversity, genetic relationship, population structure, and fingerprinting of 372 accessions of broccoli representing most of the variability of broccoli in China. Millions of SNPs were identified by whole-genome sequencing of 23 representative broccoli genotypes. Through several stringent selection criteria, a total of 1,167 SNPs were selected to characterize genetic diversity and population structure. Of these markers, 1,067 SNPs were genotyped by target sequencing (GBTS), and 100 SNPs were genotyped by kompetitive allele specific PCR (KASP) assay. The average polymorphism information content (PIC) and expected heterozygosity (gene diversity) values were 0.33 and 0.42, respectively. Diversity analysis revealed the prevalence of low to moderate genetic diversity in the broccoli accessions indicating a narrow genetic base. Phylogenetic and principal component analyses revealed that the 372 accessions could be clustered into two main groups but with weak groupings. STRUCTURE analysis also suggested the presence of two subpopulations with weak genetic structure. Analysis of molecular variance (AMOVA) identified 13% variance among populations and 87% within populations revealing very low population differentiation, which could be attributed to massive gene flow and the reproductive biology of the crop. Based on high resolving power, a set of 28 KASP markers was chosen for DNA fingerprinting of the broccoli accessions for seed authentication and varietal identification. To the best of our knowledge, this is the first comprehensive study to measure diversity and population structure of a large collection of broccoli in China and also the first application of GBTS and KASP techniques in genetic characterization of broccoli. This work broadens the understanding of diversity, phylogeny, and population structure of a large collection of broccoli, which may enhance future breeding efforts to achieve higher productivity.

Genome survey of Zanthoxylum bungeanum and development of genomic-SSR markers in congeneric species

Zanthoxylum bungeanum, a spice and medicinal plant, is cultivated in many parts of China and some countries in Southeast Asia; however, data on its genome are lacking. In the present study, we performed a whole-genome survey and developed novel genomic-SSR markers of Z. bungeanum. Clean data (∼197.16 Gb) were obtained and assembled into 11185221 scaffolds with an N50 of 183 bp. K-mer analysis revealed that Z. bungeanum has an estimated genome size of 3971.92 Mb, and the GC content, heterozygous rate, and repeat sequence rate are 37.21%, 1.73%, and 86.04%, respectively. These results indicate that the genome of Z. bungeanum is complex. Furthermore, 27153 simple sequence repeat (SSR) loci were identified from 57288 scaffolds with a minimum length > 1 kb. Mononucleotide repeats (19706) were the most abundant type, followed by dinucleotide repeats (5154). The most common motifs were A/T, followed by AT/AT; these SSRs accounted for 71.42% and 11.84% of all repeats, respectively. A total of 21243 non-repeating primer pairs were designed, and 100 were randomly selected and validated by PCR analysis using DNA from 10 Z. bungeanum individuals and 5 Zanthoxylum armatum individuals. Finally, 36 polymorphic SSR markers were developed with polymorphism information content (PIC) values ranging from 0.16 to 0.75. Cluster analysis revealed that Z. bungeanum and Z. armatum could be divided into two major clusters, suggesting that these newly developed SSR markers are useful for genetic diversity and germplasm resource identification in Z. bungeanum and Z. armatum.

Development and Characterization of 15 Novel Genomic SSRs for Viburnum farreri

The Viburnum genus is of particular interest to horticulturalists, phylogeneticists, and biogeographers. Despite its popularity, there are few existing molecular markers to investigate genetic diversity in this large genus, which includes over 160 species. There are also few polymorphic molecular tools that can delineate closely related species within the genus. Viburnum farreri, a member of the Solenotinus subclade and one of the centers of diversity for Viburnum, was selected for DNA sequencing and development of genomic simple sequence repeats (gSSRs). In this study, 15 polymorphic gSSRs were developed and characterized for a collection of 19 V. farreri samples. Number of alleles per locus ranged from two- to- eight and nine loci had four or more alleles. Observed heterozygosity ranged from 0 to 0.84 and expected heterozygosity ranged from 0.10 to 0.80 for the 15 loci. Shannon diversity index values across these loci ranged from 0.21 to 1.62. The markers developed in this study add to the existing molecular toolkit for the genus and will be used in future studies investigating cross-transferability, genetic variation, and species and cultivar delimitation in the Viburnum genus and closely allied genera in the Adoxaceae and Caprifoliaceae.

Development and technical application of SSR-based individual identification system for Chamaecyparis taiwanensis against illegal logging convictions

Chamaecyparis taiwanensis is an endemic plant suffering illegal logging in Taiwan for its high economic value. Lack of direct evidence to correlate stump and timber remains a hurdle for law enforcement. In this report, 23 polymorphic Genomic Simple Sequence Repeat (gSSR) and 12 Expressed Sequence Tag (EST)-SSR markers were developed and their transferability was assessed. The individual identification system built from selected non-linkage 30 SSR markers has a combined probability of identity as 5.596 × 10^–12 equivalents to identifying an individual in a population of up to 18 million C. taiwanensis with 99.99% confidence level. We also applied the system in an actual criminal case by selecting 19 of these markers to correlate illegally felled timbers and victim trees. Our data demonstrate that molecular signals from three timbers hit with three victim trees with confidence level more than 99.99%. This is the first example of successfully applying SSR in C. taiwanensis as a court evidence for law enforcement. The identification system adapted advanced molecular technology and exhibits its great potential for natural resource management on C. taiwanensis.

Constructing a Core Collection of the Medicinal Plant Angelica biserrata Using Genetic and Metabolic Data

Angelica biserrata is an important medicinal plant in Chinese traditional medicine. Its roots, which are known as Duhuo in Chinese, are broadly applied to treat inflammation, arthritis, and headache. With increasing market demand, the wild resources of A. biserrata have been overexploited, and conservation, assessment of genetic resources and breeding for this species is needed. Here, we sequenced the transcriptome of A. biserrata and developed simple sequence repeat (SSR) markers from it to construct a core collection based on 208 samples collected from Changyang-related regions. A total of 132 alleles were obtained for 17 SSR loci used with the polymorphic information content (PIC) ranging from 0.44 to 0.83. Abundant genetic diversity was inferred by Shannon's information index (1.51), observed (0.57) and expected heterozygosity (0.72). The clustering analysis resulted into two sample groups and analysis of molecular variance (AMOVA) showed only 6% genetic variation existed among populations. A further metabolic analysis of these samples revealed the main coumarin contents, such as osthole and columbianadin. According to the genetic and metabolic data, we adopted the least distance stepwise sampling strategy to construct seven preliminary core collections, of which the 20CC collection, which possessed 42 A. biserrata individuals accounting for 90.20% of the genetic diversity of the original germplasm, represented the best core collection. This study will contribute to the conservation and management of A. biserrata wild germplasm resources and provide a material basis for future selection and breeding of this medicinal plant.

De novo assembly, transcriptome characterization, and simple sequence repeat marker development in duckweed Lemna gibba

Lemna gibba is a species of duckweed showing great potential in bioenergy production and wastewater treatment. However, the relevant transcriptomic and genomic resources are very limited for this species, which dramatically hinders its genetic diversity and genome mapping researches. In this work, ~ 233.5 million clean reads were generated from L. gibba by Illumina paired-end sequencing, and subsequently they were de novo assembled into 131,870 unigenes, of which 61,622 were annotated and 43,319 were expressed with Fragments Per Kilobase of transcript per Million fragments mapped (FPKM) > 5. In total, 19,297 simple sequence repeats (SSRs) were identified from 15,261 SSR-containing unigenes. Dinucleotide (78.4%) were the most abundant SSRs, followed by tri- (14.9%), tetra- (4.1%), and penta-nucleotides (1.5%). The top three motifs were AG/CT (69.9%), AC/GT (6.5%), and ATC/ATG (4.9%). Further analysis revealed that the presence of SSR motif was independent of the expression level for a given gene. Based on the sequence of these SSR-containing unigenes, a total of 10,292 SSR markers were developed, of which only 2671 were further retained after removing those derived from unannotated or extra-low expressed (e.g., FPKM ≤ 5) unigenes. Finally, a subset of 70 SSR markers was randomly selected and examined in nine diverse L. gibba genotypes for the PCR amplification and polymorphism, as well as in other duckweed species for the inter-specifically amplifiability. This work is the first report on the transcriptome-based large-scale SSR markers development and analysis in L. gibba. The transcriptome generated and the SSR markers developed in this work will provide a valuable resource for genetic diversity assessment in L. gibba and also for species relationship investigation in Lemnaceae family.

The plastid genome and its implications in barcoding specific-chemotypes of the medicinal herb Pogostemon cablin in China

Pogostemon cablin (Blanco) Benth. (Patchouli) is not only an important essential oil plant, but also a valuable medicinal plant in China. P. cablin in China can be divided into three cultivars (Shipai, Gaoyao, and Hainan) and two chemotypes (pogostone-type and patchoulol-type). The pogostone-type and patchoulol-type are, respectively, used for medicinals and perfumes. In this study, we sequenced and characterized the plastid genomes for all three Chinese cultivars and aimed to develop a chemotype-specific barcode for future quality control. The plastid genomes of P. cablin cultivars ranged from 152,461 to 152,462 bp in length and comprise 114 genes including 80 protein coding genes, 30 tRNA genes, and four rRNA genes. Phylogenetic analyses suggested that P. cablin cultivars clustered with the other two Pogostemon species with strong support. Although extremely conserved in P. cablin plastid genomes, 58 cpSSRs were filtered out among the three cultivars. One single variable locus, cpSSR, was discovered. The cpSSR genotypes successfully matched the chemotypes of Chinese patchouli, which was further supported by PCR-based Sanger sequences in more Chinese patchouli samples. The barcode developed in this study is thought to be a simple and reliable quality control method for Chinese P. cablin on the market.

An RNA Sequencing Transcriptome Analysis and Development of EST-SSR Markers in Chinese Hawthorn through Illumina Sequencing

Chinese hawthorn (Crataegus pinnatifida) is an important ornamental and economic horticultural plant. However, the lack of molecular markers has limited the development and utilization of hawthorn germplasm resources. Simple sequence repeats (SSRs) derived from expressed sequence tags (ESTs) allow precise and effective cultivar characterization and are routinely used for genetic diversity analysis. Thus, we first reported the development of polymorphic EST-SSR markers in C. pinnatifida with perfect repeats using Illumina RNA-Seq technique. In total, we investigated 14,364 unigenes, from which 5091 EST-SSR loci were mined. Di-nucleotides (2012, 39.52%) were the most abundant SSRs, followed by mono- (1989, 39.07%), and tri-nucleotides (1024, 20.11%). On the basis of these EST-SSRs, a total of 300 primer pairs were designed and used for polymorphism analysis in 70 accessions collected from different geographical regions of China. Of 239 (79.67%) pairs of primer-generated amplification products, 163 (54.33%) pairs of primers showed polymorphism. Finally, 33 primers with high polymorphism were selected for genetic diversity analysis and tested on 70 individuals with low-cost fluorescence-labeled M13 primers using capillary electrophoresis genotyping platform. A total of 108 alleles were amplified by 33 SSR markers, with the number of alleles (Na) ranging from 2 to 14 per locus (mean: 4.939), and the effective number of alleles (Ne) ranging from 1.258 to 3.214 (mean: 2.221). The mean values of gene diversity (He), observed heterozygosity (Ho), and polymorphism information content (PIC) were 0.524 (range 0.205–0.689), 0.709 (range 0.132–1.000), and 0.450 (range 0.184–0.642), respectively. Furthermore, the dendrogram constructed based on the EST-SSR separated the cultivars into two main clusters. In sum, our study was the first comprehensive study on the development and analysis of a large set of SSR markers in hawthorn. The results suggested that the use of NGS techniques for SSR development represented a powerful tool for genetic studies. Additionally, fluorescence-labeled M13 markers proved to be a valuable method for genotyping. All of these EST-SSR markers have agronomic potential and constitute a scientific basis for future studies on the identification, classification, and innovation of hawthorn germplasms.

Next Generation Sequencing-Based Molecular Marker Development: A Case Study in Betula Alnoides

Betula alnoides is a fast-growing valuable indigenous tree species with multiple uses in the tropical and warm subtropical regions in South-East Asia and southern China. It has been proved to be tetraploid in most parts of its distribution in China. In the present study, next generation sequencing (NGS) technology was applied to develop numerous SSR markers for B. alnoides, and 64,376 contig sequences of 106,452 clean reads containing 164,357 candidate SSR loci were obtained. Among the derived SSR repeats, mono-nucleotide was the main type (77.05%), followed by di- (10.18%), tetra- (6.12%), tri- (3.56%), penta- (2.14%) and hexa-nucleotide (0.95%). The short nucleotide sequence repeats accounted for 90.79%. Among the 291 repeat motifs, AG/CT (46.33%) and AT/AT (44.15%) were the most common di-nucleotide repeats, while AAT/ATT (48.98%) was the most common tri-nucleotide repeats. A total of 2549 primer sets were designed from the identified putative SSR regions of which 900 were randomly selected for evaluation of amplification successfulness and detection of polymorphism if amplified successfully. Three hundred and ten polymorphic markers were obtained through testing with 24 individuals from B. alnoides natural forest in Jingxi County, Guangxi, China. The number of alleles (NA) of each marker ranged from 2 to 19 with a mean of 5.14. The observed (HO) and expected (HE) heterozygosities varied from 0.04 to 1.00 and 0.04 to 0.92 with their means being 0.64 and 0.57, respectively. Shannon-Wiener diversity index (I) ranged from 0.10 to 2.68 with a mean of 1.12. Cross-species transferability was further examined for 96 pairs of SSR primers randomly selected, and it was found that 48.96⁻84.38% of the primer pairs could successfully amplify each of six related Betula species. The obtained SSR markers can be used to study population genetics and molecular marker assisted breeding, particularly genome-wide association study of these species in the future.