Integrated metabolome and transcriptome analyses reveal the molecular mechanism underlying dynamic metabolic processes during taproot development of Panax notoginseng.

BACKGROUND: Panax notoginseng (Burk) F. H. Chen is one of the most famous Chinese traditional medicinal plants. The taproot is the main organ producing triterpenoid saponins, and its development is directly linked to the quality and yield of the harvested P. notoginseng. However, the mechanisms underlying the dynamic metabolic changes occurring during taproot development of P. notoginseng are unknown. RESULTS: We carried out metabolomic and transcriptomic analyses to investigate metabolites and gene expression during the development of P. notoginseng taproots. The differentially accumulated metabolites included amino acids and derivatives, nucleotides and derivatives, and lipids in 1-year-old taproots, flavonoids and terpenoids in 2- and 3-year-old taproots, and phenolic acids in 3-year-old taproots. The differentially expressed genes (DEGs) are related to phenylpropanoid biosynthesis, metabolic pathway and biosynthesis of secondary metabolites at all three developmental stages. Integrative analysis revealed that the phenylpropanoid biosynthesis pathway was involved in not only the development of but also metabolic changes in P. notoginseng taproots. Moreover, significant accumulation of triterpenoid saponins in 2- and 3-year-old taproots was highly correlated with the up-regulated expression of cytochrome P450s and uridine diphosphate-dependent glycosyltransferases genes. Additionally, a gene encoding RNase-like major storage protein was identified to play a dual role in the development of P. notoginseng taproots and their triterpenoid saponins synthesis. CONCLUSIONS: These results elucidate the molecular mechanism underlying the accumulation of and change relationship between primary and secondary metabolites in P. notoginseng taproots, and provide a basis for the quality control and genetic improvement of P. notoginseng.

Marsdenia tenacissima genome reveals calcium adaptation and tenacissoside biosynthesis.

Marsdenia tenacissima is a medicinal plant widely distributed in the calcium-rich karst regions of southwest China. However, the lack of a reference genome has hampered the implementation of molecular techniques in its breeding, pharmacology and domestication. We generated the chromosome-level genome assembly in Apocynaceae using combined SMRT sequencing and Hi-C. The genome length was 381.76 Mb, with 98.9% of it found on 11 chromosomes. The genome contained 222.63 Mb of repetitive sequences and 21 899 predicted gene models, with a contig N50 of 6.57 Mb. Phylogenetic analysis revealed that M. tenacissima diverged from Calotropis gigantea at least 13.43 million years ago. Comparative genomics showed that M. tenacissima underwent ancient shared whole-genome duplication. This event, together with tandem duplication, contributed to 70.71% of gene-family expansion. Both pseudogene analysis and selective pressure calculations suggested calcium-related adaptive evolution in the M. tenacissima genome. Calcium-induced differentially expressed genes (DEGs) were mainly enriched in cell-wall-related processes. Domains (e.g. Fasciclin and Amb_all) and cis-elements (e.g. MYB and MYC) frequently occurred in the coding and promoter regions of cell-wall DEGs, respectively, and the expression levels of these genes correlated significantly with those of calcium-signal-related transcription factors. Moreover, calcium addition increased tenacissoside I, G and H contents. The availability of this high-quality genome provides valuable genomic information for genetic breeding and molecular design, and lends insights into the calcium adaptation of M. tenacissima in karst areas.

Biosynthetic pathway of prescription bergenin from Bergenia purpurascens and Ardisia japonica.

Bergenin is a typical carbon glycoside and the primary active ingredient in antitussive drugs widely prescribed for central cough inhibition in China. The bergenin extraction industry relies on the medicinal plant species Bergenia purpurascens and Ardisia japonica as their resources. However, the bergenin biosynthetic pathway in plants remains elusive. In this study, we functionally characterized a shikimate dehydrogenase (SDH), two O-methyltransferases (OMTs), and a C-glycosyltransferase (CGT) involved in bergenin synthesis through bioinformatics analysis, heterologous expression, and enzymatic characterization. We found that BpSDH2 catalyzes the two-step dehydrogenation process of shikimic acid to form gallic acid (GA). BpOMT1 and AjOMT1 facilitate the methylation reaction at the 4-OH position of GA, resulting in the formation of 4-O-methyl gallic acid (4-O-Me-GA). AjCGT1 transfers a glucose moiety to C-2 to generate 2-Glucosyl-4-O-methyl gallic acid (2-Glucosyl-4-O-Me-GA). Bergenin production ultimately occurs in acidic conditions or via dehydration catalyzed by plant dehydratases following a ring-closure reaction. This study for the first time uncovered the biosynthetic pathway of bergenin, paving the way to rational production of bergenin in cell factories via synthetic biology strategies.

[Identification and expression profiling of R2R3-MYB transcription factors in Erigeron breviscapus].

R2 R3-MYB transcription factors are ubiquitous in plants, playing a role in the regulation of plant growth, development, and secondary metabolism. In this paper, the R2 R3-MYB transcription factors were identified by bioinformatics analysis of the genomic data of Erigeron breviscapus, and their gene sequences, structures, physical and chemical properties were analyzed. The functions of R2 R3-MYB transcription factors were predicted by cluster analysis. Meanwhile, the expression patterns of R2 R3-MYB transcription factors in response to hormone treatments were analyzed. A total of 108 R2 R3-MYB transcription factors, named EbMYB1-EbMYB108, were identified from the genome of E. breviscapus. Most of the R2 R3-MYB genes carried 2-4 exons. The phylogenetic tree of MYBs in E. breviscapus and Arabidopsis thaliala was constructed, which classified 234 MYBs into 30 subfamilies. The MYBs in the five MYB subfamilies of A.thaliala were clustered into independent clades, and those in E. breviscapus were clustered into four clades. The transcriptome data showed that MYB genes were differentially expressed in different tissues of E. breviscapus and in response to the treatments with exogenous hormones such as ABA, SA, and GA for different time. The transcription of 13 R2 R3-MYB genes did not change significantly, and the expression patterns of some genes were up-regulated or down-regulated with the extension of hormone treatment time. This study provides a theoretical basis for revealing the mechanisms of R2 R3-MYB transcription factors in regulating the growth and development, stress(hormone) response, and active ingredient accumulation in E. breviscapus.

Identification of two UDP-glycosyltransferases involved in the main oleanane-type ginsenosides in Panax japonicus var. major.

MAIN CONCLUSION: Two UDP-glycosyltransferases from Panax japonicus var. major were identified, and the biosynthetic pathways of three oleanane-type ginsenosides (chikusetsusaponin IVa, ginsenoside Ro, zingibroside R1) were elucidated. Chikusetsusaponin IVa and ginsenoside Ro are primary active components formed by stepwise glycosylation of oleanolic acid in five medicinal plants of the genus Panax. However, the key UDP-glycosyltransferases (UGTs) in the biosynthetic pathway of chikusetsusaponin IVa and ginsenoside Ro are still unclear. In this study, two UGTs (PjmUGT1 and PjmUGT2) from Panax japonicus var. major involved in the biosynthesis of chikusetsusaponin IVa and ginsenoside Ro were identified based on bioinformatics analysis, heterologous expression and enzyme assays. The results show that PjmUGT1 can transfer a glucose moiety to the C-28 carboxyl groups of oleanolic acid 3-O-ß-D-glucuronide and zingibroside R1 to form chikusetsusaponin IVa and ginsenoside Ro, respectively. Meanwhile, PjmUGT2 can transfer a glucose moiety to oleanolic acid 3-O-ß-D-glucuronide and chikusetsusaponin IVa to form zingibroside R1 and ginsenoside Ro. This work uncovered the biosynthetic mechanism of chikusetsusaponin IVa and ginsenoside Ro, providing the rational production of valuable saponins through synthetic biology strategy.

High quality genome of Erigeron breviscapus provides a reference for herbal plants in Asteraceae.

Erigeron breviscapus is an important medicinal plant in Compositae and the first species to realize the whole process from the decoding of the draft genome sequence to scutellarin biosynthesis in yeast. However, the previous low-quality genome assembly has hindered the optimization of candidate genes involved in scutellarin synthesis and the development of molecular-assisted breeding based on the genome. Here, the E. breviscapus genome was updated using PacBio RSII sequencing data and Hi-C data, and increased in size from 1.2 Gb to 1.43 Gb, with a scaffold N50 of 156.82 Mb and contig N50 of 140.95 kb, and a total of 43,514 protein-coding genes were obtained and oriented onto nine pseudo-chromosomes, thus becoming the third plant species assembled to chromosome level after sunflower and lettuce in Compositae. Fourteen genes with evidence for positive selection were identified and found to be related to leaf morphology, flowering and secondary metabolism. The number of genes in some gene families involved in flavonoid biosynthesis in E. breviscapus have been significantly expanded. In particular, additional candidate genes involved in scutellarin biosynthesis, such as flavonoid-7-O-glucuronosyltransferase genes (F7GATs) were identified using updated genome. In addition, three candidate genes encoding indole-3-pyruvate monooxygenase YUCCA2 (YUC2), serine carboxypeptidase-like 18 (SCPL18), and F-box protein (FBP), respectively, were identified to be probably related to leaf development and flowering by resequencing 99 individuals. These results provided a substantial genetic basis for improving agronomic and quality traits of E. breviscapus, and provided a platform for improving other draft genome assemblies to chromosome-level.

[Phosphorus Storage Capacity and Loss Risk in Coastal Reed Wetland Surrounding Bohai Sea].

Coastal wetland, at the intersection of land and sea, is considered as a "sink", "source", and "transformer" of phosphorus (P). Coastal wetland plays an important role in the global P cycle, and its ability to retain excessive P in water receives increasing attention. In this study, the coastal reed wetland sediments surrounding the Bohai Sea were sampled to investigate P adsorption capacity and loss risk by conducting batch experiments. Results show that the maximum P adsorption capacity (Qmax) was 693.7-2117.2 mg·kg-1, with an average of 1468.6 mg·kg-1. The Qmax decreased in the order of Qilihai Wetland > Beidagang Wetland > Nandagang Wetland > Liaohe Delta Wetland > Shouguang Coastal Wetland > Yellow River Delta Wetland. The P adsorption capacity was related to the contents of Ca, Mg, and TOC. The degree of P adsorption saturation (DPS) and loss risk index (ERI) of the coastal wetland were 0.28%-4.50% and 0.53%-10.10%, respectively. The ERI suggested that the P loss risk was relatively low for coastal reed wetland surrounding the Bohai Sea except for the moderate loss risk for Shouguang coastal wetland. In summary, the reed coastal wetland around Bohai Sea demonstrated significantly P storage capacity and served as a P sink for water P. We recommend to fully utilize the coastal wetland to reduce point or non-point source pollution (e.g., P) during the remediation or recovery of the polluted Bohai Sea.

EbARC1, an E3 Ubiquitin Ligase Gene in Erigeron breviscapus, Confers Self-Incompatibility in Transgenic Arabidopsis thaliana.

Erigeron breviscapus (Vant.) Hand.-Mazz. is a famous traditional Chinese medicine that has positive effects on the treatment of cardiovascular and cerebrovascular diseases. With the increase of market demand (RMB 500 million per year) and the sharp decrease of wild resources, it is an urgent task to cultivate high-quality and high-yield varieties of E. breviscapus. However, it is difficult to obtain homozygous lines in breeding due to the self-incompatibility (SI) of E. breviscapus. Here, we first proved that E. breviscapus has sporophyte SI (SSI) characteristics. Characterization of the ARC1 gene in E. breviscapus showed that EbARC1 is a constitutive expression gene located in the nucleus. Overexpression of EbARC1 in Arabidopsis thaliana L. (Col-0) could cause transformation of transgenic lines from self-compatibility (SC) into SI. Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays indicated that EbARC1 and EbExo70A1 interact with each other in the nucleus, and the EbARC1-ubox domain and EbExo70A1-N are the key interaction regions, suggesting that EbARC1 may ubiquitinate EbExo70A to regulate SI response. This study of the SSI mechanism in E. breviscapus has laid the foundation for further understanding SSI in Asteraceae and breeding E. breviscapus varieties.

Comparative transcriptome and metabolome analyses provide new insights into the molecular mechanisms underlying taproot thickening in Panax notoginseng.

BACKGROUND: Taproot thickening is a complex biological process that is dependent on the coordinated expression of genes controlled by both environmental and developmental factors. Panax notoginseng is an important Chinese medicinal herb that is characterized by an enlarged taproot as the main organ of saponin accumulation. However, the molecular mechanisms of taproot enlargement are poorly understood. RESULTS: A total of 29,957 differentially expressed genes (DEGs) were identified during the thickening process in the taproots of P. notoginseng. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment revealed that DEGs associated with "plant hormone signal transduction," "starch and sucrose metabolism," and "phenylpropanoid biosynthesis" were predominantly enriched. Further analysis identified some critical genes (e.g., RNase-like major storage protein, DA1-related protein, and Starch branching enzyme I) and metabolites (e.g., sucrose, glucose, fructose, malate, and arginine) that potentially control taproot thickening. Several aspects including hormone crosstalk, transcriptional regulation, homeostatic regulation between sugar and starch, and cell wall metabolism, were identified as important for the thickening process in the taproot of P. notoginseng. CONCLUSION: The results provide a molecular regulatory network of taproot thickening in P. notoginseng and facilitate the further characterization of the genes responsible for taproot formation in root medicinal plants or crops.

Promoter analysis of the mouse Peg3 gene.

Mouse Peg3 is a paternally expressed gene. Study of methylation of the Peg3 gene in P19 embryonal carcinoma cells suggested that monoallelic methylation of CpG dinucleotides is not only present in the promoter region, but also in the first exon and the first intron. Promoter activity analysis demonstrated that the minimal promoter of the Peg3 gene is located in the region between -827 and +712 and the critical region for promoter activity is between +423 and +712. We further identified the roles of the cis-elements, conserved sequence element (CSE) and YY1-binding sites, in the regulation of Peg3 expression and found that CSE is involved in the inhibition of Peg3 expression, while YY1-binding sites serve as activating cis-elements to antagonize CSE-mediated inhibition.

Cyclophilin a protects Peg3 from hypermethylation and inactive histone modification.

Imprinted genes are expressed from only one of the parental alleles and are marked epigenetically by DNA methylation and histone modifications. Disruption of normal imprinting leads to abnormal embryogenesis, certain inherited diseases, and is associated with various cancers. In the context of screening for the gene(s) responsible for the alteration of phenotype in cyclophilin A knockdown (CypA-KD) P19 cells, we observed a silent paternally expressed gene, Peg3. Treatment of CypA-KD P19 cells with the DNA demethylating agent 5-aza-dC reversed the silencing of Peg3 biallelically. Genomic bisulfite sequencing and methylation-specific PCR revealed DNA hypermethylation in CypA-KD P19 cells, as the normally unmethylated paternal allele acquired methylation that resulted in biallelic methylation of Peg3. Chromatin immunoprecipitation assays indicated a loss of acetylation and a gain of lysine 9 trimethylation in histone 3, as well as enhanced DNA methyltransferase 1 and MBD2 binding on the cytosine-guanine dinucleotide (CpG) islands of Peg3. Our results indicate that DNA hypermethylation on the paternal allele and allele-specific acquisition of histone methylation leads to silencing of Peg3 in CypA-KD P19 cells. This study is the first demonstration of the epigenetic function of CypA in protecting the paternal allele of Peg3 from DNA methylation and inactive histone modifications.