Combined physiological and transcriptome analysis revealed the response mechanism of Pogostemon cablin roots to p-hydroxybenzoic acid.

Pogostemon cablin (patchouli) cultivation is challenged by serious soil sickness, of which autotoxins accumulation is a major cause. p-hydroxybenzoic acid (p-HBA) is one of the main autotoxins of patchouli. However, the molecular mechanism underlying the response of patchouli to p-HBA remains unclear. In this study, RNA-sequencing combined with physiological analysis was used to monitor the dynamic transcriptomic and physiological changes in patchouli seedlings 0, 6, 12, 24, 48, and 96 h after p-HBA treatment. p-HBA stress inhibited root biomass accumulation, induced excessive hydrogen peroxide accumulation and lipid peroxidation, and activated most antioxidant enzymes. Compared with that of the control, the osmotic adjustment substance content was elevated with treatment. Subsequently, 15,532, 8,217, 8,946, 2,489, and 5,843 differentially expressed genes (DEGs) at 6, 12, 24, 48, and 96 h after p-HBA treatment, respectively, were identified in patchouli roots. GO functional enrichment analysis showed that the DEGs were enriched mainly in plasma membrane, defense response, response to chitin, DNA-binding transcription factor activity and abscisic acid-activated signaling pathway. The upregulated genes were involved in glycolysis/gluconeogenesis, cysteine and methionine metabolism, starch and sucrose metabolism, biosynthesis of unsaturated fatty acids, and linoleic acid metabolism. Genes associated with MAPK signaling pathway-plant, plant-pathogen interaction, plant hormone signal transduction were downregulated with p-HBA treatment. These pathways are related to root browning and rotting, leading to plant death.

Physiological and Transcriptional Analysis Reveals the Response Mechanism of Camellia vietnamensis Huang to Drought Stress.

Drought stress is considered the main obstacle restricting Camellia vietnamensis Huang (C. vietnamensis) yield. Hainan is the southernmost distribution region of C. vietnamensis in China and experiences a drought period annually. To study the drought-stress-response mechanism of C. vietnamensis, we treated seedlings of drought-tolerant (HD1) and drought-sensitive (WH1) cultivars with PEG-6000 (PEG) to simulate drought stress and compared the physiology and transcriptome of their leaves at 0 d, 3 d and 6 d posttreatment. Under drought stress, the growth of C. vietnamensis was inhibited, the relative water content (RWC) of leaves decreased and the contents of malondialdehyde (MDA), antioxidant enzyme activities, osmotic regulatory substances and secondary metabolites increased. Compared with those of WH1, the leaf RWC, osmotic-regulation substance content (proline, soluble protein and soluble sugar) and antioxidant enzyme activity (superoxide dismutase, peroxidase and catalase) of HD1 were significantly increased, while the relative electrical conductivity and MDA content were significantly decreased. Compared with WH1, 2812, 2070 and 919, differentially expressed genes (DEGs) were detected in HD1 0 d, 3 d and 6 d posttreatment, respectively, and the number of DEGs increased with increasing treatment time. The detected DEGs are involved in the drought stress response of C. vietnamensis mainly through plant-hormone signal transduction and lignin and flavonoid biosynthesis pathways. Drought stress significantly activated the expression of several lignin and flavonoid biosynthesis genes in HD1. Moreover, total flavonoid and total polyphenol contents in HD1 were significantly increased, suggesting that the accumulation of flavonoids may be a key factor in the drought stress response of C. vietnamensis. Additionally, 191 DEGs were associated with coding transcription factors (TFs). This study provides insight into the molecular mechanism of the drought stress response of C. vietnamensis and provides a theoretical basis for the development and cultivation of new drought-resistant cultivars.

Genetic relationships and low diversity among the tea-oil Camellia species in Sect. Oleifera, a bulk woody oil crop in China.

Tea-oil Camellia is one of the four woody oil crops in the world and has high ecological, economic and medicinal values. However, there are great differences in the classification and merging of tea-oil Camellia Sect. Oleifera species, which brings difficulties to the innovative utilization and production of tea-oil Camellia resources. Here, ISSR, SRAP and chloroplast sequence markers were analyzed in 18 populations of tea-oil Camellia Sect. Oleifera species to explore their phylogenetic relationships and genetic diversity. The results showed that their genetic diversity were low, with mean H and π values of 0.16 and 0.00140, respectively. There was high among-population genetic differentiation, with ISSR and SRAP markers showing an Fst of 0.38 and a high Nm of 1.77 and cpDNA markers showing an Fst of 0.65 and a low Nm of 0.27. The C. gauchowensis, C. vietnamensis and Hainan Island populations formed a single group, showing the closest relationships, and supported being the same species for them with the unifying name C. drupifera and classifying the resources on Hainan Island as C. drupifera. The tea-oil Camellia resources of Hainan Island should be classified as a special ecological type or variety of C. drupifera. However, cpDNA marker-based STRUCTURE analysis showed that the genetic components of the C. osmantha population formed an independent, homozygous cluster; hence, C. osmantha should be a new species in Sect. Oleifera. The C. oleifera var. monosperma and C. oleifera populations clustered into two distinct clades, and the C. oleifera var. monosperma populations formed an independent cluster, accounting for more than 99.00% of its genetic composition; however, the C. oleifera populations contained multiple different cluster components, indicating that C. oleifera var. monosperma significantly differs from C. oleifera and should be considered the independent species C. meiocarpa. Haplotype analysis revealed no rapid expansion in the tested populations, and the haplotypes of C. oleifera, C. meiocarpa and C. osmantha evolved from those of C. drupifera. Our results support the phylogenetic classification of Camellia subgenera, which is highly significant for breeding and production in tea-oil Camellia.

Integrated Analysis of Physiological, mRNA Sequencing, and miRNA Sequencing Data Reveals a Specific Mechanism for the Response to Continuous Cropping Obstacles in Pogostemon cablin Roots.

Pogostemon cablin (patchouli) is a commercially important medicinal and industrial crop grown worldwide for its medicinal and aromatic properties. Patchoulol and pogostone, derived from the essential oil of patchouli, are considered valuable components in the cosmetic and pharmaceutical industries. Due to its high application value in the clinic and industry, the demand for patchouli is constantly growing. Unfortunately, patchouli cultivation has suffered due to severe continuous cropping obstacles, resulting in a significant decline in yield and quality. Moreover, the physiological and transcriptional changes in patchouli in response to continuous cropping obstacles remain unclear. This has greatly restricted the development of the patchouli industry. To explore the mechanism underlying the rapid response of patchouli roots to continuous cropping stress, integrated analysis of the transcriptome and miRNA profiles of patchouli roots under continuous and noncontinuous cropping conditions in different growth periods was conducted using RNA sequencing (RNA-seq) and miRNA-seq and complemented with physiological data. The physiological and biochemical results showed that continuous cropping significantly inhibited root growth, decreased root activity, and increased the activity of antioxidant enzymes (superoxide dismutase, peroxidase, and catalase) and the levels of osmoregulators (malondialdehyde, soluble protein, soluble sugar, and proline). Subsequently, we found 4,238, 3,494, and 7,290 upregulated and 4,176, 3,202, and 8,599 downregulated differentially expressed genes (DEGs) in the three growth periods of continuously cropped patchouli, many of which were associated with primary carbon and nitrogen metabolism, defense responses, secondary metabolite biosynthesis, and transcription factors. Based on miRNA-seq, 927 known miRNAs and 130 novel miRNAs were identified, among which 67 differentially expressed miRNAs (DEMIs) belonging to 24 miRNA families were induced or repressed by continuous cropping. By combining transcriptome and miRNA profiling, we obtained 47 miRNA-target gene pairs, consisting of 18 DEMIs and 43 DEGs, that likely play important roles in the continuous cropping response of patchouli. The information provided in this study will contribute to clarifying the intricate mechanism underlying the patchouli response to continuous cropping obstacles. In addition, the candidate miRNAs and genes can provide a new strategy for breeding continuous cropping-tolerant patchouli.

Preparation of Active Peptides from Camellia vietnamensis and Their Metabolic Effects in Alcohol-Induced Liver Injury Cells.

Liver damage seriously affects human health. Over 35% of cases of acute liver damage are caused by alcohol damage. Thus, finding drugs that can inhibit and effectively treat this disease is necessary. This article mainly focuses on the effect of the metabolome physical activity of active peptides in Camellia vietnamensis active peptide (CMAP) and improving liver protection. DEAE Sepharose FF ion-exchange column chromatography was used in separating and purifying crude peptides from Camellia vietnamensis Two components, A1 and A2, were obtained, and the most active A1 was selected. Sephadex G-100 gel column chromatography was used in A1 separation and purification. Three components, Al-1, Al-2, and Al-3, were obtained. Through antioxidant activity in vitro as an index of inspection, the relatively active component A1-2 was removed. Reverse-phase high-performance liquid chromatography showed that the purity of component A1-2 was 93.45%. The extracted CMAPs acted on alcoholic liver injury cells. Metabolomics studies revealed that the up-regulated metabolites were ribothymidine and xanthine; the down-regulated metabolites were hydroxyphenyllactic acid, creatinine, stearoylcarnitine, and inosine. This study provides an effective theoretical support for subsequent research.

Transcriptome analysis of two Pogostemon cablin chemotypes reveals genes related to patchouli alcohol biosynthesis.

Pogostemon cablin, a medicinally and economically important perennial herb, is cultivated around the world due to its medicinal and aromatic properties. Different P. cablin cultivars exhibit different morphological traits and patchouli oil components and contents (especially patchouli alcohol (PA) and pogostone (PO)). According to the signature constituent of the leaf, P. cablin was classified into two different chemotypes, including PA-type and PO-type. To better understand the molecular mechanisms of PA biosynthesis, the transcriptomes of Chinese-cultivated P. cablin cv. PA-type "Nanxiang" (NX) and PO-type "Paixiang" (PX) were analyzed and compared with ribonucleic acid sequencing (RNA-Seq) technology. We obtained a total of 36.83 G clean bases from the two chemotypes, compared them with seven databases and revealed 45,394 annotated unigenes. Thirty-six candidate unigenes participating in the biosynthesis of PA were found in the P. cablin transcriptomes. Overall, 8,390 differentially expressed unigenes were identified between the chemotypes, including 2,467 upregulated and 5,923 downregulated unigenes. Furthermore, six and nine differentially expressed genes (DEGs) were mapped to the terpenoid backbone biosynthetic and sesquiterpenoid and triterpenoid biosynthetic pathways, respectively. One key sesquiterpene synthase gene involved in the sesquiterpenoid and triterpenoid biosynthetic pathways, encoding patchoulol synthase variant 1, was significantly upregulated in NX. Additionally, GC-MS analysis of the two chemotypes in this study showed that the content of PA in NX was significantly higher than that of PX, while the content of PO showed the opposite phenotype. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis showed that the DEG expression tendency was consistent with the transcriptome sequencing results. Overall, 23 AP2/ERF, 13 bHLH, 11 MYB, 11 NAC, three Trihelix, 10 WRKY and three bZIP genes that were differentially expressed may act as regulators of terpenoid biosynthesis. Altogether, 8,314 SSRs were recognized within 6,825 unigenes, with a distribution frequency of 18.32%, among which 1,202 unigenes contained more than one SSR. The transcriptomic characteristics of the two P. cablin chemotypes are comprehensively reported in this study, and these results will contribute to a better understanding of the molecular mechanism of PA biosynthesis. Our transcriptome data also provide a valuable genetic resource for further studies on P. cablin.

Integrative iTRAQ-based proteomic and transcriptomic analysis reveals the accumulation patterns of key metabolites associated with oil quality during seed ripening of Camellia oleifera.

Camellia oleifera (C. oleifera) is one of the four major woody oil-bearing crops in the world and has relatively high ecological, economic, and medicinal value. Its seeds undergo a series of complex physiological and biochemical changes during ripening, which is mainly manifested as the accumulation and transformation of certain metabolites closely related to oil quality, especially flavonoids and fatty acids. To obtain new insights into the underlying molecular mechanisms, a parallel analysis of the transcriptome and proteome profiles of C. oleifera seeds at different maturity levels was conducted using RNA sequencing (RNA-seq) and isobaric tags for relative and absolute quantification (iTRAQ) complemented with gas chromatography-mass spectrometry (GC-MS) data. A total of 16,530 transcripts and 1228 proteins were recognized with significant differential abundances in pairwise comparisons of samples at various developmental stages. Among these, 317 were coexpressed with a poor correlation, and most were involved in metabolic processes, including fatty acid metabolism, α-linolenic acid metabolism, and glutathione metabolism. In addition, the content of total flavonoids decreased gradually with seed maturity, and the levels of fatty acids generally peaked at the fat accumulation stage; these results basically agreed with the regulation patterns of genes or proteins in the corresponding pathways. The expression levels of proteins annotated as upstream candidates of phenylalanine ammonia-lyase (PAL) and chalcone synthase (CHS) as well as their cognate transcripts were positively correlated with the variation in the flavonoid content, while shikimate O-hydroxycinnamoyltransferase (HCT)-encoding genes had the opposite pattern. The increase in the abundance of proteins and mRNAs corresponding to alcohol dehydrogenase (ADH) was associated with a reduction in linoleic acid synthesis. Using weighted gene coexpression network analysis (WGCNA), we further identified six unique modules related to flavonoid, oil, and fatty acid anabolism that contained hub genes or proteins similar to transcription factors (TFs), such as MADS intervening keratin-like and C-terminal (MIKC_MADS), type-B authentic response regulator (ARR-B), and basic helix-loop-helix (bHLH). Finally, based on the known metabolic pathways and WGCNA combined with the correlation analysis, five coexpressed transcripts and proteins composed of cinnamyl-alcohol dehydrogenases (CADs), caffeic acid 3-O-methyltransferase (COMT), flavonol synthase (FLS), and 4-coumarate: CoA ligase (4CL) were screened out. With this exploratory multiomics dataset, our results presented a dynamic picture regarding the maturation process of C. oleifera seeds on Hainan Island, not only revealing the temporal specific expression of key candidate genes and proteins but also providing a scientific basis for the genetic improvement of this tree species.

Complementary transcriptome and proteome profiling in the mature seeds of Camellia oleifera from Hainan Island.

Camellia oleifera Abel. (C. oleifera), as an important woody tree species producing edible oils in China, has attracted enormous attention due to its abundant unsaturated fatty acids and their associated benefits to human health. To reveal novel insights into the characters during the maturation period of this plant as well as the molecular basis of fatty acid biosynthesis and degradation, we conducted a conjoint analysis of the transcriptome and proteome of C. oleifera seeds from Hainan Island. Using RNA sequencing (RNA-seq) technology and shotgun proteomic method, 59,391 transcripts and 40,500 unigenes were obtained by TIGR Gene Indices Clustering Tools (TGICL), while 1691 protein species were identified from Mass Spectrometry (MS). Subsequently, all genes and proteins were employed in euKaryotic Orthologous Groups (KOG) classification, Gene Ontology (GO) annotation, and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis to investigate their essential functions. The results indicated that the most abundant pathways were biological metabolic processes. There were 946 unigenes associated with lipid metabolism at the transcriptome level, with 116 proteins at the proteome level; among these, 38 specific proteins were involved in protein-protein interactions, with the majority being related to fatty acid catabolic process. The expression levels of 21 candidate unigenes encoding target proteins were further detected by quantitative real-time polymerase chain reaction (qRT-PCR). Finally, Gas Chromatography Mass Spectrometry (GC-MS) was carried out to determine the fatty acid composition of C. oleifera oil. These findings not only deepened our understanding about the molecular mechanisms of fatty acid metabolism but also offered new evidence concerning the roles of relevant proteins in oil-bearing crops. Furthermore, the lipid-associated proteins recognized in this research might be helpful in providing a reference for the synthetic regulation of C. oleifera oil quality by genetic engineering techniques, thus resulting in potential application in agriculture.

Analysis of genetic population structure and diversity in Mallotus oblongifolius using ISSR and SRAP markers.

BACKGROUND: Mallotus oblongifolius, an evergreen shrub endemic to Hainan Island, China, is important both medicinally and economically. Due to its special medicinal significance and the continuing rise of market demand, its populations in the wild have been subject to long-term illegal and unrestrained collection. Hence, an evaluation of genetic variability is essential for the conservation and genetic reserve development of this species. METHODS: Sequence-related amplified polymorphism (SRAP) and inter-simple sequence repeat (ISSR) markers were employed to assess the genetic diversity and genetic structure of 20 natural populations of M. oblongifolius growing in different eco-geographical regions of Hainan Island, China. RESULTS: We revealed a considerable genetic diversity (h = 0.336, I = 0.5057, SRAP markers; h = 0.3068, I = 0.4657, ISSR markers) and weak genetic differentiation (Gst = 0.2764 for SRAP, Gst = 0.2709 for ISSR) with the same gene flow (Nm = 1.3092 for SRAP, Nm = 1.346 for ISSR) among the M. oblongifolius populations. The Mantel Test showed that the distribution of genetic variation among populations could not be explained by the pronounced geographical distances (r = 0.01255, p = 0.5538). All results of the Unweighted Pair Group Method with Arithmetic Mean (UPGMA), Neighbor-joining (NJ), Principal Coordinate Analysis (PCoA) and Bayesian analyses supported a habitat-specific genetic clustering model for M. oblongifolius, indicating a local adaptive divergence for the studied populations. DISCUSSION: We suggested that the habitat fragmentation and specificity for M. oblongifolius populations weakened the natural gene flow and promoted an adaptation to special habitats, which was the main reason for local adaptive divergence among M. oblongifolius.