De novo transcriptome assembly of Pueraria montana var. lobata and Neustanthus phaseoloides for the development of eSSR and SNP markers: narrowing the US origin(s) of the invasive kudzu

De Novo Assembly and Comparative Analysis of Mitochondrial Genomes of Two Pueraria montana Varieties

Pueraria montana is a species with important medicinal value and a complex genetic background. In this study, we sequenced and assembled the mitochondrial (mt) genomes of two varieties of P. montana. The mt genome lengths of P. montana var. thomsonii and P. montana var. montana were 457,390 bp and 456,731 bp, respectively. Both P. montana mitogenomes showed a multi-branched structure consisting of two circular molecules, with 56 genes annotated, comprising 33 protein-coding genes, 18 tRNA genes (trnC-GCA and trnM-CAU are multi-copy genes), and 3 rRNA genes. Then, 207 pairs of long repeats and 96 simple sequence repeats (SSRs) were detected in the mt genomes of P. montana, and 484 potential RNA-editing sites were found across the 33 mitochondrial protein-coding genes of each variety. Additionally, a syntenic sequence analysis showed a high collinearity between the two mt genomes. This work is the first to analyze the mt genomes of P. montana. It can provide information that can be used to analyze the structure of mt genomes of higher plants and provide a foundation for future comparative genomic studies and evolutionary biology research in related species.

Pueraria montana Population Structure and Genetic Diversity Based on Chloroplast Genome Data

Despite having a generally conserved structure, chloroplast genome data have been helpful for plant population genetics and evolution research. To mine Pueraria montana chloroplast genome variation architecture and phylogeny, we investigated the chloroplast variation architecture of 104 P. montana accessions from across China. P. montana's chloroplast genome showed high diversity levels, with 1674 variations, including 1118 single nucleotide polymorphisms and 556 indels. The intergenic spacers, psbZ-trnS and ccsA-ndhD, are the two mutation hotspot regions in the P. montana chloroplast genome. Phylogenetic analysis based on the chloroplast genome dataset supported four P. montana clades. P. montana variations were conserved among and within clades, which showed high gene flow levels. Most P. montana clades were estimated to have diverged at 3.82-5.17 million years ago. Moreover, the East Asian summer monsoon and South Asian summer monsoon may have accelerated population divergence. Our results show that chloroplast genome sequences were highly variable and can be used as molecular markers to assess genetic variation and relationships in P. montana.

Comparative plastomes of Pueraria montana var. lobata (Leguminosae: Phaseoleae) and closely related taxa: insights into phylogenomic implications and evolutionary divergence

BACKGROUND

Pueraria montana var. lobata (kudzu) is an important food and medicinal crop in Asia. However, the phylogenetic relationships between Pueraria montana var. lobata and the other two varieties (P. montana var. thomsonii and P. montana var. montana) remain debated. Although there is increasing evidence showing that P. montana var. lobata adapts to various environments and is an invasive species in America, few studies have systematically investigated the role of the phylogenetic relationships and evolutionary patterns of plastomes between P. montana var. lobata and its closely related taxa.

RESULTS

26 newly sequenced chloroplast genomes of Pueraria accessions resulted in assembled plastomes with sizes ranging from 153,360 bp to 153,551 bp. Each chloroplast genome contained 130 genes, including eight rRNA genes, 37 tRNA genes, and 85 protein-coding genes. For 24 newly sequenced accessions of these three varieties of P. montana, we detected three genes and ten noncoding regions with higher nucleotide diversity (π). After incorporated publically available chloroplast genomes of Pueraria and other legumes, 47 chloroplast genomes were used to construct phylogenetic trees, including seven P. montana var. lobata, 14 P. montana var. thomsonii and six P. montana var. montana. Phylogenetic analysis revealed that P. montana var. lobata and P. montana var. thomsonii formed a clade, while all sampled P. montana var. montana formed another cluster based on cp genomes, LSC, SSC and protein-coding genes. Twenty-six amino acid residues were identified under positive selection with the site model. We also detected six genes (accD, ndhB, ndhC, rpl2, rpoC2, and rps2) that account for among-site variation in selective constraint under the clade model between accessions of the Pueraria montana var. lobata clade and the Pueraria montana var. montana clade.

CONCLUSION

Our data provide novel comparative plastid genomic insights into conservative gene content and structure of cp genomes pertaining to P. montana var. lobata and the other two varieties, and reveal an important phylogenetic clue and plastid divergence among related taxa of P. montana come from loci that own moderate variation and underwent modest selection.

Chromosomal-level genome and multi-omics dataset of Pueraria lobata var. thomsonii provide new insights into legume family and the isoflavone and puerarin biosynthesis pathways

Pueraria lobata var. thomsonii (hereinafter abbreviated as Podalirius thomsonii), a member of legumes, is one of the important traditional Chinese herbal medicines, and its puerarin extraction is widely used in health and pharmaceutical industry. Here, we assembled a high-quality genome of P. thomsonii using long-read single-molecule sequencing and Hi-C technologies. The genome assembly is approximately 1.37 Gb in size and consists of 5145 contigs with a contig N50 of 593.70 Kb, further clustered into 11 pseudochromosomes. The genome structural annotation resulted in about 869.33 Mb (about 62.70% of the genome) repeat regions and 45 270 protein-coding genes. Genome evolution analysis revealed that P. thomsonii is most closely related to soybean and underwent two ancient whole-genome duplication events, one was in the common ancestor shared by legume species, the other occurred independently at around 7.2 million years ago after its specification. A total of 2373 gene families were found unique in P. thomsonii comparing to five other legume species. Genes and metabolites related to puerarin content in tuberous tissues were characterized. A total of 572 genes upregulated in the puerarin biosynthesis pathway were identified, and 235 candidate genes were further enriched by omics data. Furthermore, we identified 6 8-C-glucosyltransferase (8-C-GT) candidate genes significantly involved in puerarin metabolism. Our study filled in a key genomic gap in legume family, and provided valuable multi-omic resources for the genetic improvement of P. thomsonii.

Identification of Pueraria spp. through DNA barcoding and comparative transcriptomics

BACKGROUND

Kudzu is a term used generically to describe members of the genus Pueraria. Kudzu roots have been used for centuries in traditional Chinese medicine in view of their high levels of beneficial isoflavones including the unique 8-C-glycoside of daidzein, puerarin. In the US, kudzu is seen as a noxious weed causing ecological and economic damage. However, not all kudzu species make puerarin or are equally invasive. Kudzu remains difficult to identify due to its diverse morphology and inconsistent nomenclature.

RESULTS

We have generated sequences for the internal transcribed spacer 2 (ITS2) and maturase K (matK) regions of Pueraria montana lobata, P. montana montana, and P. phaseoloides, and identified two accessions previously used for differential analysis of puerarin biosynthesis as P. lobata and P. phaseoloides. Additionally, we have generated root transcriptomes for the puerarin-producing P. m. lobata and the non-puerarin producing P. phaseoloides. Within the transcriptomes, microsatellites were identified to aid in species identification as well as population diversity.

CONCLUSIONS

The barcode sequences generated will aid in fast and efficient identification of the three kudzu species. Additionally, the microsatellites identified from the transcriptomes will aid in genetic analysis. The root transcriptomes also provide a molecular toolkit for comparative gene expression analysis towards elucidation of the biosynthesis of kudzu phytochemicals.

Identification of Nutritional Ingredients and Medicinal Components of Pueraria lobata and Its Varieties Using UPLC-MS/MS-Based Metabolomics

Pueraria lobata and its variety P. lobata var. thomsonii are both traditional Chinese medicines that have high nutritional and medical value; whereas another variety, P. lobata var. montana has low nutritional and medicinal value and can cause ecological disasters. The material basis of different nutritional and medicinal values, which are caused by metabolite differences among these varieties, remains to be further clarified. Here, we performed ultra performance liquid chromatography-tandem mass spectrometry based widely targeted metabolome analysis on Pueraria lobata, P. lobata var. thomsonii, and P. lobata var. montana. Among them, a total of 614 metabolites were identified, and distinguished from each other using orthogonal partial least squares discriminant analysis. Our results suggest that the nutritional differences between P. lobata and its varieties can be explained by variations in the abundance of amino acids, nucleotides, saccharides, and lipids; differences in flavonoids, isoflavones, phenolic acids, organic acids, and coumarins contents caused the differences in the medicinal quality of P. lobata and its varieties. Additionally, the key metabolites responsible for the classification of the three Pueraria varieties were identified. This study provides new insights into the underlying metabolic causes of nutritional and medicinal variation in P. lobata and its varieties.

The development of SSR markers based on RNA-sequencing and its validation between and within Carex L. species

BACKGROUND

Carex L. is one of the largest genera in the Cyperaceae family and an important vascular plant in the ecosystem. However, the genetic background of Carex is complex and the classification is not clear. In order to investigate the gene function annotation of Carex, RNA-sequencing analysis was performed. Simple sequence repeats (SSRs) were generated based on the Illumina data and then were utilized to investigate the genetic characteristics of the 79 Carex germplasms.

RESULTS

In this study, 36,403 unigenes with a total length of 41,724,615 bp were obtained and annotated based on GO, KOG, KEGG, NR databases. The results provide a theoretical basis for gene function exploration. Out of 8776 SSRs, 96 pairs of primers were randomly selected. One hundred eighty polymorphic bands were amplified with a polymorphism rate of 100% based on 42 pairs of primers with higher polymorphism levels. The average band number was 4.3 per primer, the average distance value was 0.548, and the polymorphic information content was ranged from 0.133 to 0.494. The number of observed alleles (Na), effective alleles (Ne), Nei's (1973) gene diversity (H), and the Shannon information index (I) were 2.000, 1.376, 0.243, and 0.391, respectively. NJ clustering divided into three groups and the accessions from New Zealand showed a similar genetic attribute and clustered into one group. UPGMA and PCoA analysis also revealed the same result. The analysis of molecular variance (AMOVA) revealed a superior genetic diversity within accessions than between accessions based on geographic origin cluster and NJ cluster. What's more, the fingerprints of 79 Carex species are established in this study. Different combinations of primer pairs can be used to identify multiple Carex at one time, which overcomes the difficulties of traditional identification methods.

CONCLUSIONS

The transcriptomic analysis shed new light on the function categories from the annotated genes and will facilitate future gene functional studies. The genetic characteristics analysis indicated that gene flow was extensive among 79 Carex species. These markers can be used to investigate the evolutionary history of Carex and related species, as well as to serve as a guide in future breeding projects.

In-depth transcriptome reveals the potential biotechnological application of Bothrops jararaca venom gland

Background:

Lack of complete genomic data of Bothrops jararaca impedes molecular biology research focusing on biotechnological applications of venom gland components. Identification of full-length coding regions of genes is crucial for the correct molecular cloning design.

Methods:

RNA was extracted from the venom gland of one adult female specimen of Bothrops jararaca. Deep sequencing of the mRNA library was performed using Illumina NextSeq 500 platform. De novo assembly of B. jararaca transcriptome was done using Trinity. Annotation was performed using Blast2GO. All predicted proteins after clustering step were blasted against non-redundant protein database of NCBI using BLASTP. Metabolic pathways present in the transcriptome were annotated using the KAAS-KEGG Automatic Annotation Server. Toxins were identified in the B. jararaca predicted proteome using BLASTP against all protein sequences obtained from Animal Toxin Annotation Project from Uniprot KB/Swiss-Pro database. Figures and data visualization were performed using ggplot2 package in R language environment.

Results:

We described the in-depth transcriptome analysis of B. jararaca venom gland, in which 76,765 de novo assembled isoforms, 96,044 transcribed genes and 41,196 unique proteins were identified. The most abundant transcript was the zinc metalloproteinase-disintegrin-like jararhagin. Moreover, we identified 78 distinct functional classes of proteins, including toxins, inhibitors and tumor suppressors. Other venom proteins identified were the hemolytic lethal factors stonustoxin and verrucotoxin.

Conclusion:

It is believed that the application of deep sequencing to the analysis of snake venom transcriptomes may represent invaluable insight on their biotechnological potential focusing on candidate molecules.