Integrated widely targeted metabolomics and network pharmacology revealed quality disparities between Guizhou and conventional producing areas of Codonopsis Radix.

Introduction: With the internationalization of traditional Chinese medicine, the demand for medicinal and edible Codonopsis Radix (CR) has increased, and its medicinal resources have attracted attention. CR is a well-known traditional Chinese medicine with a long pharmaceutical and edible history. The Guizhou province in China has abundant CR resources, but in the absence of systematic studies on species identification and chemical compositions, the capacity of the capacity of the province to CR resource has not been fully utilized. Methodology: We used plant morphology and DNA barcoding techniques to identify Luodang (LD) and Weidang (WD) species. To investigate the differences in metabolites between LD and WD, as well as three Chinese Pharmacopeia CRs, and to predict pharmacological mechanisms of action for the dominant differential metabolites, we utilized widely targeted metabolomics and network pharmacology. The results also revealed the material basis for the excellent food properties of both LD and WD. Results: The plant traits and DNA barcoding molecular identification results indicated that Luodang and Weidang from Guizhou were Codonopsis tangshen and Codonopsis pilosula, respectively. Widely targeted metabolomics analysis revealed that a total of 1,116 metabolites from 14 categories, including phenolic acids, lipids, flavonoids, were found in five CRs and shared 1,054 (94.4%) metabolites. LD and WD each contained 3 and 10 dominant differential metabolites, respectively, which were primarily flavonoids and amino acids. Amino acids, phenolic acids, and organic acids play important roles in their excellent food attributes. In CR, eight dominant differential metabolites were discovered for the first time, including isoorientin-7-O-(6″-feruloyl) glucoside, N-formyl-L-methionine, and cyclo (Phe-Glu), among others. Network pharmacology analyses showed that, in LD, dominant differential metabolites were closely related to anti-tumor, cardiovascular disease improvement, nervous system protection, and metabolic disease treatment, whereas in WD, they were closely related to nervous system protection and cardiovascular disease improvement. Conclusion: The species of LD and WD were included in the Chinese Pharmacopeia, and their metabolite profiles were remarkably similar to CR from traditional producing areas. Therefore, LD and WD can be used and promoted medicinally as CR, and they have potential value for new drug development. This study enriched the database of CR compounds and provided a reference for quality control, resource development, and new drug development of CR.

Regulatory Mechanisms of the Wnt/ß-Catenin Pathway in Diabetic Cutaneous Ulcers.

Skin ulcers are a serious complication of diabetes. Diabetic patients suffer from vascular lesions and complications such as peripheral neuritis, peripheral vascular lesions, and collagen abnormalities, which result in skin wounds that are refractory and often develop into chronic ulcers. The healing of skin ulcers requires an inflammatory reaction, wound proliferation, remodeling regulation, and control of stem cells. Studies investigating diabetic cutaneous ulcers have focused on cellular and molecular levels. Diabetes can cause nerve and blood vessel damage, and persistent high blood sugar levels can cause systemic multisite nerve damage based on peripheral neuropathy. The long-term hyperglycemia state enables the polyol glucose metabolism pathway to be activated, increasing the accumulation of toxic substances in the vascular injured nerve tissue cells. Sustained hyperglycemia leads to dysfunction of epithelial cells, leading to a decrease in pro-angiogenic signaling and nitric oxide production. In addition, due to impaired leukocyte function in hyperglycemia, immune function is impaired and the immune response at relevant sites is insufficient, making diabetic foot more difficult to heal. The Wnt/ß-catenin pathway is a highly conserved signal transduction pathway involved in a variety of biological processes, such as cell proliferation, apoptosis, and differentiation. It is considered an important pathway involved in the healing of skin wounds. This article summarizes the mechanism of action of the Wnt/ß-catenin pathway involved in the inflammatory responses to diabetic ulcers, wound proliferation, wound remodeling, and stem cells. The interactions between the Wnt signal pathway and other metabolic pathways are also discussed.

Metabolomics and Transcriptomics Analyses Reveal Nitrogen Influences on the Accumulation of Flavonoids and Amino Acids in Young Shoots of Tea Plant ( Camellia sinensis L.) Associated with Tea Flavor.

Tea-specialized metabolites contribute to rich flavors and healthy function of tea. Their accumulation patterns and underlying regulatory mechanism are significantly different under different nitrogen (N) conditions during adaptation stage. Here, we find that flavonoids associated with tea flavor are dominated by different metabolic and transcriptional responses among the four N conditions (N-deficiency, nitrate, ammonia, and nitric oxide). Nitrogen-deficiency tea plants accumulate diverse flavonoids, corresponding with higher expression of hub genes including F3H, FNS, UFGT, bHLH35, and bHLH36. Compared with N-deficiency, N-supply tea plants significantly increase proline, glutamine, and theanine, which are also associated with tea flavor, especially under NH4+-supply. As NH4+-tolerant species, tea plant exploits the adaptive strategy by substantial accumulation of amino acids including theanine to adapt excess NH4+, which attributes to, at least in part, efficient N transport and assimilation, and active protein degradation. A distinct divergence of N reallocation in young shoots of tea plant under different N sources contributes to diverse tea flavor.

De novo transcriptome sequencing of Camellia sasanqua and the analysis of major candidate genes related to floral traits.

Camellia sasanqua is one of the most famous horticultural plants in Camellia (Theaceae) due to its aesthetic appeal as landscape plant. Knowledge regarding the genetic basis of flowering time, floral aroma and color in C. sasanqua is limited, but is essential to breed new varieties with desired floral traits. Here, we described the de novo transcriptome of young leaves, flower buds and flowers of C. sasanqua. A total of 60,127 unigenes were functionally annotated based on the sequence similarity. After analysis, we found that two floral integrator genes, SOC1 and AP1, in flowering time pathway showed evidence of gene family expansion. Compared with 1-deoxy-D-xylulose-5-phosphate pathway, some genes in the mevalonate pathway were most highly expressed, suggesting that this might represent the major pathway for terpenoid biosynthesis related to floral aroma in C. sasanqua. In flavonoid biosynthesis pathway, PAL, CHI, DFR and ANS showing significantly higher expression levels in flowers and flower buds might have important role in regulation of floral color. The top five most transcription factors (TFs) families in C. sasanqua transcriptome were MYB, MIKC, C3H, FAR1 and HD-ZIP, many of which have a direct relationship with floral traits. In addition, we also identified 33,540 simple sequence repeats (SSRs) in the C. sasanqua transcriptome. Collectively, the C. sasanqua transcriptome dataset generated from this study along with the SSR markers provide a new resource for the identification of novel regulatory transcripts and will accelerate the genetic improvement of C. sasanqua breeding programs.

Transcriptome Analysis Identifies Candidate Genes Related to Triacylglycerol and Pigment Biosynthesis and Photoperiodic Flowering in the Ornamental and Oil-Producing Plant, Camellia reticulata (Theaceae).

Camellia reticulata, which is native to Southwest China, is famous for its ornamental flowers and high-quality seed oil. However, the lack of genomic information for this species has largely hampered our understanding of its key pathways related to oil production, photoperiodic flowering process and pigment biosynthesis. Here, we first sequenced and characterized the transcriptome of a diploid C. reticulata in an attempt to identify genes potentially involved in triacylglycerol biosynthesis (TAGBS), photoperiodic flowering, flavonoid biosynthesis (FlaBS), carotenoid biosynthesis (CrtBS) pathways. De novo assembly of the transcriptome provided a catalog of 141,460 unigenes with a total length of ~96.1 million nucleotides (Mnt) and an N50 of 1080 nt. Of them, 22,229 unigenes were defined as differentially expressed genes (DEGs) across five sequenced tissues. A large number of annotated genes in C. reticulata were found to have been duplicated, and differential expression patterns of these duplicated genes were commonly observed across tissues, such as the differential expression of SOC1_a, SOC1_b, and SOC1_c in the photoperiodic flowering pathway. Up-regulation of SAD_a and FATA genes and down-regulation of FAD2_a gene in the TAGBS pathway in seeds may be relevant to the ratio of monounsaturated fatty acid (MUFAs) to polyunsaturated fatty acid (PUFAs) in seed oil. MYBF1, a transcription regulator gene of the FlaBS pathway, was found with great sequence variation and alteration of expression patterns, probably resulting in functionally evolutionary differentiation in C. reticulata. MYBA1_a and some anthocyanin-specific biosynthetic genes in the FlaBS pathway were highly expressed in both flower buds and flowers, suggesting important roles of anthocyanin biosynthesis in flower development. Besides, a total of 40,823 expressed sequence tag simple sequence repeats (EST-SSRs) were identified in the C. reticulata transcriptome, providing valuable marker resources for further basic and applied researches on this economically important Camellia plant.

Autotetraploid rice methylome analysis reveals methylation variation of transposable elements and their effects on gene expression.

Polyploidy, or whole-genome duplication (WGD), serves as a key innovation in plant evolution and is an important genomic feature for all eukaryotes. Neopolyploids have to overcome difficulties in meiosis, genomic alterations, changes of gene expression, and epigenomic reorganization. However, the underlying mechanisms for these processes are poorly understood. One of the most interesting aspects is that genome doubling events increase the dosage of all genes. Unlike allopolyploids entangled by both hybridization and polyploidization, autopolyploids, especially artificial lines, in relatively uniform genetic background offer a model system to understand mechanisms of genome-dosage effects. To investigate DNA methylation effects in response to WGD rather than hybridization, we produced autotetraploid rice with its diploid donor, Oryza sativa ssp. indica cv. Aijiaonante, both of which were independently self-pollinated over 48 generations, and generated and compared their comprehensive transcriptomes, base pair-resolution methylomes, and siRNAomes. DNA methylation variation of transposable elements (TEs) was observed as widespread in autotetraploid rice, in which hypermethylation of class II DNA transposons was predominantly noted in CHG and CHH contexts. This was accompanied by changes of 24-nt siRNA abundance, indicating the role of the RNA-directed DNA methylation pathway. Our results showed that the increased methylation state of class II TEs may suppress the expression of neighboring genes in autotetraploid rice that has obtained double alleles, leading to no significant differences in transcriptome alterations for most genes from its diploid donor. Collectively, our findings suggest that chromosome doubling induces methylation variation in TEs that affect gene expression and may become a "genome shock" response factor to help neoautopolyploids adapt to genome-dosage effects.

De Novo Assembly and Characterization of the Transcriptome of the Chinese Medicinal Herb, Gentiana rigescens.

Gentiana rigescens is an important medicinal herb in China. The main validated medicinal component gentiopicroside is synthesized in shoots, but is mainly found in the plant's roots. The gentiopicroside biosynthetic pathway and its regulatory control remain to be elucidated. Genome resources of gentian are limited. Next-generation sequencing (NGS) technologies can aid in supplying global gene expression profiles. In this study we present sequence and transcript abundance data for the root and leaf transcriptome of G. rigescens, obtained using the Illumina Hiseq2000. Over fifty million clean reads were obtained from leaf and root libraries. This yields 76,717 unigenes with an average length of 753 bp. Among these, 33,855 unigenes were identified as putative homologs of annotated sequences in public protein and nucleotide databases. Digital abundance analysis identified 3306 unigenes differentially enriched between leaf and root. Unigenes found in both tissues were categorized according to their putative functional categories. Of the differentially expressed genes, over 130 were annotated as related to terpenoid biosynthesis. This work is the first study of global transcriptome analyses in gentian. These sequences and putative functional data comprise a resource for future investigation of terpenoid biosynthesis in Gentianaceae species and annotation of the gentiopicroside biosynthetic pathway and its regulatory mechanisms.

CandiSSR: An Efficient Pipeline used for Identifying Candidate Polymorphic SSRs Based on Multiple Assembled Sequences.

Simple sequence repeats (SSRs), also known as microsatellites, are ubiquitous short tandem duplications commonly found in genomes and/or transcriptomes of diverse organisms. They represent one of the most powerful molecular markers for genetic analysis and breeding programs because of their high mutation rate and neutral evolution. However, traditionally experimental screening of the SSR polymorphic status and their subsequent applicability to genetic studies are extremely labor-intensive and time-consuming. Thankfully, the recently decreased costs of next generation sequencing and increasing availability of large genome and/or transcriptome sequences have provided an excellent opportunity and sources for large-scale mining this type of molecular markers. However, current tools are limited. Thus we here developed a new pipeline, CandiSSR, to identify candidate polymorphic SSRs (PolySSRs) based on the multiple assembled sequences. The pipeline allows users to identify putative PolySSRs not only from the transcriptome datasets but also from multiple assembled genome sequences. In addition, two confidence metrics including standard deviation and missing rate of the SSR repetitions are provided to systematically assess the feasibility of the detected PolySSRs for subsequent application to genetic characterization. Meanwhile, primer pairs for each identified PolySSR are also automatically designed and further evaluated by the global sequence similarities of the primer-binding region, ensuring the successful rate of the marker development. Screening rice genomes with CandiSSR and subsequent experimental validation showed an accuracy rate of over 90%. Besides, the application of CandiSSR has successfully identified a large number of PolySSRs in the Arabidopsis genomes and Camellia transcriptomes. CandiSSR and the PolySSR marker sources are publicly available at: http://www.plantkingdomgdb.com/CandiSSR/index.html.

Genome-wide identification and comparative expression analysis reveal a rapid expansion and functional divergence of duplicated genes in the WRKY gene family of cabbage, Brassica oleracea var. capitata.

WRKY transcription factors (TFs), one of the ten largest TF families in higher plants, play important roles in regulating plant development and resistance. To date, little is known about the WRKY TF family in Brassica oleracea. Recently, the completed genome sequence of cabbage (B. oleracea var. capitata) allows us to systematically analyze WRKY genes in this species. A total of 148 WRKY genes were characterized and classified into seven subgroups that belong to three major groups. Phylogenetic and synteny analyses revealed that the repertoire of cabbage WRKY genes was derived from a common ancestor shared with Arabidopsis thaliana. The B. oleracea WRKY genes were found to be preferentially retained after the whole-genome triplication (WGT) event in its recent ancestor, suggesting that the WGT event had largely contributed to a rapid expansion of the WRKY gene family in B. oleracea. The analysis of RNA-Seq data from various tissues (i.e., roots, stems, leaves, buds, flowers and siliques) revealed that most of the identified WRKY genes were positively expressed in cabbage, and a large portion of them exhibited patterns of differential and tissue-specific expression, demonstrating that these gene members might play essential roles in plant developmental processes. Comparative analysis of the expression level among duplicated genes showed that gene expression divergence was evidently presented among cabbage WRKY paralogs, indicating functional divergence of these duplicated WRKY genes.

Rapid diversification of five Oryza AA genomes associated with rice adaptation.

Comparative genomic analyses among closely related species can greatly enhance our understanding of plant gene and genome evolution. We report de novo-assembled AA-genome sequences for Oryza nivara, Oryza glaberrima, Oryza barthii, Oryza glumaepatula, and Oryza meridionalis. Our analyses reveal massive levels of genomic structural variation, including segmental duplication and rapid gene family turnover, with particularly high instability in defense-related genes. We show, on a genomic scale, how lineage-specific expansion or contraction of gene families has led to their morphological and reproductive diversification, thus enlightening the evolutionary process of speciation and adaptation. Despite strong purifying selective pressures on most Oryza genes, we documented a large number of positively selected genes, especially those genes involved in flower development, reproduction, and resistance-related processes. These diversifying genes are expected to have played key roles in adaptations to their ecological niches in Asia, South America, Africa and Australia. Extensive variation in noncoding RNA gene numbers, function enrichment, and rates of sequence divergence might also help account for the different genetic adaptations of these rice species. Collectively, these resources provide new opportunities for evolutionary genomics, numerous insights into recent speciation, a valuable database of functional variation for crop improvement, and tools for efficient conservation of wild rice germplasm.