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

High-integrity Pueraria montana var. lobata genome and population analysis revealed the genetic diversity of Pueraria genus

Pueraria montana var. lobata (P. lobata) is a traditional medicinal plant belonging to the Pueraria genus of Fabaceae family. Pueraria montana var. thomsonii (P. thomsonii) and Pueraria montana var. montana (P. montana) are its related species. However, evolutionary history of the Pueraria genus is still largely unknown. Here, a high-integrity, chromosome-level genome of P. lobata and an improved genome of P. thomsonii were reported. It found evidence for an ancient whole-genome triplication and a recent whole-genome duplication shared with Fabaceae in three Pueraria species. Population genomics of 121 Pueraria accessions demonstrated that P. lobata populations had substantially higher genetic diversity, and P. thomsonii was probably derived from P. lobata by domestication as a subspecies. Selection sweep analysis identified candidate genes in P. thomsonii populations associated with the synthesis of auxin and gibberellin, which potentially play a role in the expansion and starch accumulation of tubers in P. thomsonii. Overall, the findings provide new insights into the evolutionary and domestication history of the Pueraria genome and offer a valuable genomic resource for the genetic improvement of these species.

Integrated Metabolomic and Transcriptomic Analysis of Puerarin Biosynthesis in Pueraria montana var. thomsonii at Different Growth Stages

Puerarin, a class of isoflavonoid compounds concentrated in the roots of Puerarias, has antipyretic, sedative, and coronary blood-flow-increasing properties. Although the biosynthetic pathways of puerarin have been investigated by previous researchers, studies focusing on the influence of different growth stages on the accumulation of metabolites in the puerarin pathway are not detailed, and it is still controversial at the last step of the 8-C-glycosylation reaction. In this study, we conducted a comprehensive analysis of the metabolomic and transcriptomic changes in Pueraria montana var. thomsonii during two growing years, focusing on the vigorous growth and dormant stages, to elucidate the underlying mechanisms governing the changes in metabolite and gene expression within the puerarin biosynthesis pathway. In a comparison of the two growth stages in the two groups, puerarin and daidzin, the main downstream metabolites in the puerarin biosynthesis pathway, were found to accumulate mainly during the vigorous growth stage. We also identified 67 common differentially expressed genes in this pathway based on gene expression differences at different growth stages. Furthermore, we identified four candidate 8-C-GT genes that potentially contribute to the conversion of daidzein into puerarin and eight candidate 7-O-GT genes that may be involved in the conversion of daidzein into daidzin. A co-expression network analysis of important UGTs and HIDs along with daidzein and puerarin was conducted. Overall, our study contributes to the knowledge of puerarin biosynthesis and offers information about the stage at which the 8-C-glycosylation reaction occurs in biosynthesis. These findings provide valuable insights into the cultivation and quality enhancement of Pueraria montana var. thomsonii.

A telomere-to-telomere reference genome provides genetic insight into the pentacyclic triterpenoid biosynthesis in Chaenomeles speciosa

Chaenomeles speciosa (2n = 34), a medicinal and edible plant in the Rosaceae, is commonly used in traditional Chinese medicine. To date, the lack of genomic sequence and genetic studies has impeded efforts to improve its medicinal value. Herein, we report the use of an integrative approach involving PacBio HiFi (third-generation) sequencing and Hi-C scaffolding to assemble a high-quality telomere-to-telomere genome of C. speciosa. The genome comprised 650.4 Mb with a contig N50 of 35.5 Mb. Of these, 632.3 Mb were anchored to 17 pseudo-chromosomes, in which 12, 4, and 1 pseudo-chromosomes were represented by a single contig, two contigs, and four contigs, respectively. Eleven pseudo-chromosomes had telomere repeats at both ends, and four had telomere repeats at a single end. Repetitive sequences accounted for 49.5% of the genome, while a total of 45 515 protein-coding genes have been annotated. The genome size of C. speciosa was relatively similar to that of Malus domestica. Expanded or contracted gene families were identified and investigated for their association with different plant metabolisms or biological processes. In particular, functional annotation characterized gene families that were associated with the biosynthetic pathway of oleanolic and ursolic acids, two abundant pentacyclic triterpenoids in the fruits of C. speciosa. Taken together, this telomere-to-telomere and chromosome-level genome of C. speciosa not only provides a valuable resource to enhance understanding of the biosynthesis of medicinal compounds in tissues, but also promotes understanding of the evolution of the Rosaceae.

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.

Genome-wide identification and expression pattern analysis of R2R3-MYB transcription factor gene family involved in puerarin biosynthesis and response to hormone in Pueraria lobata var. thomsonii

BACKGROUND

R2R3-MYB transcription factors regulate secondary metabolism, stress responses and development in various plants. Puerarin is a bioactive ingredient and most abundant secondary metabolite isolated from Pueraria lobata. The biosynthesis of puerarin proceeds via the phenylpropanoid pathway and isoflavonoids pathway, in which 9 key enzymes are involved. The expression of these structural genes is under control of specific PtR2R3-MYB genes in different plant tissues. However, how PtR2R3-MYB genes regulates structural genes in puerarin biosynthesis remains elusive. This study mined the PtR2R3-MYB genes involved in puerarin biosynthesis and response to hormone in Pueraria lobata var. thomsonii.

RESULTS

A total of 209 PtR2R3-MYB proteins were identified, in which classified into 34 subgroups based on the phylogenetic topology and the classification of the R2R3-MYB superfamily in Arabidopsis thaliana. Furtherly physical and chemical characteristics, gene structure, and conserved motif analysis were also used to further analyze PtR2R3-MYBs. Combining puerarin content and RNA-seq data, speculated on the regulated puerarin biosynthesis of PtR2R3-MYB genes and structural genes, thus 21 PtR2R3-MYB genes and 25 structural genes were selected for validation gene expression and further explore its response to MeJA and GSH treatment by using qRT-PCR analysis technique. Correlation analysis and cis-acting element analysis revealed that 6 PtR2R3-MYB genes (PtMYB039, PtMYB057, PtMYB080, PtMYB109, PtMYB115 and PtMYB138) and 7 structural genes (PtHID2, PtHID9, PtIFS3, PtUGT069, PtUGT188, PtUGT286 and PtUGT297) were directly or indirectly regulation of puerarin biosynthesis in ZG11. It is worth noting that after MeJA and GSH treatment for 12-24 h, the expression changes of most candidate genes were consistent with the correlation of puerarin biosynthesis, which also shows that MeJA and GSH have the potential to mediate puerarin biosynthesis by regulating gene expression in ZG11.

CONCLUSIONS

Overall, this study provides a comprehensive understanding of the PtR2R3-MYB and will paves the way to reveal the transcriptional regulation of puerarin biosynthesis and response to phytohormone of PtR2R3-MYB genes in Pueraria lobata var. thomsonii.

Comparative Transcriptome Analysis of Pueraria lobata Provides Candidate Genes Involved in Puerarin Biosynthesis and Its Regulation

Pueraria lobata is a traditional Chinese herb in which an isoflavone C-glucoside, namely puerarin, has received the utmost interest due to its medicinal properties. To date, the biogenesis of puerarin, especially its C-glucosyl reaction in the pathway, remains poorly understood. Moreover, the transcription factors (TFs) that regulate puerarin biosynthesis in P. lobata have not been reported. Here, we performed phytochemical studies on the different developmental stages of the root, stem, and leaf tissues of two P. lobata cultivars, which suggested that both the roots and stems of P. lobata were the sites of puerarin biosynthesis. RNA-sequencing was conducted with the root and stem tissues of the two cultivars under different stages, and the clean reads were mapped to the recently published genome of P. lobata var. thomsonii, yielding the transcriptome dataset. A detailed analysis of the gene expression data, gene coexpression network, and phylogeny proposed several C-GTs that likely participate in puerarin biosynthesis. The first genome-wide analysis of the whole MYB superfamily in P. lobata presented here identified a total of 123 nonredundant PlMYB genes that were significantly expressed in the analyzed tissues. The phylogenetic analysis of PlMYBs with other plant MYB proteins revealed strong PlMYB candidates that may regulate the biosynthesis of isoflavones, such as puerarin.

Chromosome-level and graphic genomes provide insights into metabolism of bioactive metabolites and cold-adaption of Pueraria lobata var. montana

Pueraria lobata var. montana (P. montana) belongs to the genus Pueraria and originated in Asia. Compared with its sister P. thomsonii, P. montana has stronger growth vigour and cold-adaption but contains less bioactive metabolites such as puerarin. To promote the investigation of metabolic regulation and genetic improvement of Pueraria, the present study reports a chromosome-level genome of P. montana with length of 978.59 Mb and scaffold N50 of 80.18 Mb. Comparative genomics analysis showed that P. montana possesses smaller genome size than that of P. thomsonii owing to less repeat sequences and duplicated genes. A total of 6,548 and 4,675 variety-specific gene families were identified in P. montana and P. thomsonii, respectively. The identified variety-specific and expanded/contracted gene families related to biosynthesis of bioactive metabolites and microtubules are likely the causes for the different characteristics of metabolism and cold-adaption of P. montana and P. thomsonii. Moreover, a graphic genome was constructed based on 11 P. montana accessions. Total 92 structural variants were identified and most of which are related to stimulus-response. In conclusion, the chromosome-level and graphic genomes of P. montana will not only facilitate the studies of evolution and metabolic regulation, but also promote the breeding of Pueraria.