Subgroup 4 R2R3-MYBs in conifer trees: gene family expansion and contribution to the isoprenoid- and flavonoid-oriented responses

PaMYB11 promotes suberin deposition in Norway spruce embryogenic tissue during cryopreservation: A novel resistance mechanism against osmosis

The osmotic resistance mechanism has been extensively studied in whole plants or plant tissues. However, little is known about it in embryogenic tissue (ET) which is widely used in plant-based biotechnological systems. Suberin, a cell wall aliphatic and aromatic heteropolymer, plays a critical role in plant cells against osmosis stress. The suberin regulatory biosynthesis has rarely been studied in gymnosperms. Here, PaMYB11, a subgroup 11 R2R3-MYB transcription factor, plays a key role in the osmotic resistance of Norway spruce (Picea abies) ETs during cryoprotectant pretreatment. Thus, RNA-seq, histological, and analytical chemical analyses are performed on the stable transformations of PaMYB11-OE and PaMYB11-SRDX in Norway spruce ETs. DAP-seq, Y1H, and LUC are further combined to explore the PaMYB11 targets. Activation of PaMYB11 is necessary and sufficient for suberin lamellae deposition on Norway spruce embryogenic cell walls, which plays a decisive role in ET survival under osmotic stress. Transcriptome analysis shows that PaMYB11 enhances suberin lamellae monomer synthesis by promoting very long-chain fatty acid (VLCFA) synthesis. PaPOP, PaADH1, and PaTET8L, the first two (PaADH1 and PaPOP, included) involved in VLCFA synthesis, are proved to be the direct targets of PaMYB11. Our study identified a novel osmotic response directed by PaMYB11 in Norway spruce ET, which provides a new understanding of the resistance mechanism against osmosis in gymnosperms.

Chromosome-scale genome assembly of Docynia delavayi provides new insights into the α-farnesene biosynthesis

Docynia delavayi is an economically significant fruit species with a high market potential due to the special aroma of its fruit. Here, a 653.34 Mb high-quality genome of D. delavayi was first reported, of which 93.8 % of the sequences (612.98 Mb) could be anchored to 17 chromosomes, containing 48,325 protein-coding genes. Ks analysis proved that two whole genome duplication (WGD) events occurred in D. delavayi, resulting in the expansion of genes associated with terpene biosynthesis, which promoted its fruit-specific aroma production. Combined multi-omics analysis, α-farnesene was detected as the most abundant aroma substance emitted by D. delavayi fruit during storage, meanwhile one α-farnesene synthase gene (AFS) and 15 transcription factors (TFs) were identified as the candidate genes potentially involved in α-farnesene biosynthesis. Further studies for the regulation network of α-farnesene biosynthesis revealed that DdebHLH, DdeERF1 and DdeMYB could activate the transcription of DdeAFS. To our knowledge, it is the first report that MYB TF plays a regulatory role in α-farnesene biosynthesis, which will greatly facilitate future breeding programs for D. delavayi.

Conifers Concentrate Large Numbers of NLR Immune Receptor Genes on One Chromosome

Nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes form a major line of defense in plants, acting in both pathogen recognition and resistance machinery activation. NLRs are reported to form large gene clusters in limber pine (Pinus flexilis), but it is unknown how widespread this genomic architecture may be among the extant species of conifers (Pinophyta). We used comparative genomic analyses to assess patterns in the abundance, diversity, and genomic distribution of NLR genes. Chromosome-level whole genome assemblies and high-density linkage maps in the Pinaceae, Cupressaceae, Taxaceae, and other gymnosperms were scanned for NLR genes using existing and customized pipelines. The discovered genes were mapped across chromosomes and linkage groups and analyzed phylogenetically for evolutionary history. Conifer genomes are characterized by dense clusters of NLR genes, highly localized on one chromosome. These clusters are rich in TNL-encoding genes, which seem to have formed through multiple tandem duplication events. In contrast to angiosperms and nonconiferous gymnosperms, genomic clustering of NLR genes is ubiquitous in conifers. NLR-dense genomic regions are likely to influence a large part of the plant's resistance, informing our understanding of adaptation to biotic stress and the development of genetic resources through breeding.

High-density genetic linkage mapping in Sitka spruce advances the integration of genomic resources in conifers

In species with large and complex genomes such as conifers, dense linkage maps are a useful resource for supporting genome assembly and laying the genomic groundwork at the structural, populational, and functional levels. However, most of the 600+ extant conifer species still lack extensive genotyping resources, which hampers the development of high-density linkage maps. In this study, we developed a linkage map relying on 21,570 single nucleotide polymorphism (SNP) markers in Sitka spruce (Picea sitchensis [Bong.] Carr.), a long-lived conifer from western North America that is widely planted for productive forestry in the British Isles. We used a single-step mapping approach to efficiently combine RAD-seq and genotyping array SNP data for 528 individuals from 2 full-sib families. As expected for spruce taxa, the saturated map contained 12 linkages groups with a total length of 2,142 cM. The positioning of 5,414 unique gene coding sequences allowed us to compare our map with that of other Pinaceae species, which provided evidence for high levels of synteny and gene order conservation in this family. We then developed an integrated map for P. sitchensis and Picea glauca based on 27,052 markers and 11,609 gene sequences. Altogether, these 2 linkage maps, the accompanying catalog of 286,159 SNPs and the genotyping chip developed, herein, open new perspectives for a variety of fundamental and more applied research objectives, such as for the improvement of spruce genome assemblies, or for marker-assisted sustainable management of genetic resources in Sitka spruce and related species.

Long-insert sequence capture detects high copy numbers in a defence-related beta-glucosidase gene βglu-1 with large variations in white spruce but not Norway spruce

Conifers are long-lived and slow-evolving, thus requiring effective defences against their fast-evolving insect natural enemies. The copy number variation (CNV) of two key acetophenone biosynthesis genes Ugt5/Ugt5b and βglu-1 may provide a plausible mechanism underlying the constitutively variable defence in white spruce (Picea glauca) against its primary defoliator, spruce budworm. This study develops a long-insert sequence capture probe set (Picea_hung_p1.0) for quantifying copy number of βglu-1-like, Ugt5-like genes and single-copy genes on 38 Norway spruce (Picea abies) and 40 P. glauca individuals from eight and nine provenances across Europe and North America respectively. We developed local assemblies (Piabi_c1.0 and Pigla_c.1.0), full-length transcriptomes (PIAB_v1 and PIGL_v1), and gene models to characterise the diversity of βglu-1 and Ugt5 genes. We observed very large copy numbers of βglu-1, with up to 381 copies in a single P. glauca individual. We observed among-provenance CNV of βglu-1 in P. glauca but not P. abies. Ugt5b was predominantly single-copy in both species. This study generates critical hypotheses for testing the emergence and mechanism of extreme CNV, the dosage effect on phenotype, and the varying copy number of genes with the same pathway. We demonstrate new approaches to overcome experimental challenges in genomic research in conifer defences.

Assembly and annotation of the black spruce genome provide insights on spruce phylogeny and evolution of stress response

Black spruce (Picea mariana [Mill.] B.S.P.) is a dominant conifer species in the North American boreal forest that plays important ecological and economic roles. Here, we present the first genome assembly of P. mariana with a reconstructed genome size of 18.3 Gbp and NG50 scaffold length of 36.0 kbp. A total of 66,332 protein-coding sequences were predicted in silico and annotated based on sequence homology. We analyzed the evolutionary relationships between P. mariana and 5 other spruces for which complete nuclear and organelle genome sequences were available. The phylogenetic tree estimated from mitochondrial genome sequences agrees with biogeography; specifically, P. mariana was strongly supported as a sister lineage to P. glauca and 3 other taxa found in western North America, followed by the European Picea abies. We obtained mixed topologies with weaker statistical support in phylogenetic trees estimated from nuclear and chloroplast genome sequences, indicative of ancient reticulate evolution affecting these 2 genomes. Clustering of protein-coding sequences from the 6 Picea taxa and 2 Pinus species resulted in 34,776 orthogroups, 560 of which appeared to be specific to P. mariana. Analysis of these specific orthogroups and dN/dS analysis of positive selection signatures for 497 single-copy orthogroups identified gene functions mostly related to plant development and stress response. The P. mariana genome assembly and annotation provides a valuable resource for forest genetics research and applications in this broadly distributed species, especially in relation to climate adaptation.

Integrative analysis of non-targeted metabolome and transcriptome reveals the mechanism of volatile formation in pepper fruit

Introduction: Aroma is a key inherent quality attributes of pepper fruit, yet the underlying mechanisms of aroma compound biosynthesis remain unclear. Methods: In this study, the volatile profile of the QH (cultivated Capsicum chinense) and WH (cultivated Capsicum annuum) pepper varieties were putatively identified during fruit development using gas chromatography-mass spectrometry (GC-MS). Results and discussion: The results identified 203 volatiles in pepper, and most of the esters, terpenes, aldehydes and alcohols were significantly down-regulated with fruit ripening. The comparison of volatile components between varieties revealed that aldehydes and alcohols were highly expressed in the WH fruit, while esters and terpenes with fruity or floral aroma were generally highly accumulated in the QH fruit, providing QH with a fruity odor. Transcriptome analysis demonstrated the close relationship between the synthesis of volatiles and the fatty acid and terpene metabolic pathways, and the high expression of the ADH, AAT and TPS genes was key in determining the accumulation of volatiles in pepper fruit. Furthermore, integrative metabolome and transcriptome analysis revealed that 208 differentially expressed genes were highly correlated with 114 volatiles, and the transcription factors of bHLH, MYB, ARF and IAA were identified as fundamental for the regulation of volatile synthesis in pepper fruit. Our results extend the understanding of the synthesis and accumulation of volatiles in pepper fruit.

Wound-Induced Temporal Reprogramming of Gene Expression during Agarwood Formation in Aquilaria sinensis

Agarwood is a resinous heartwood of Aquilaria sinensis that is formed in response to mechanical wounding. However, the transcriptional response of A. sinensis to mechanical wounding during the agarwood formation process is still unclear. Here, three five-year-old A. sinensis trees were mechanically damaged by a chisel, and time-series transcriptomic analysis of xylem tissues in the treated area (TA) was performed at 15 (TA1), 70 (TA2) and 180 days after treatment (TA3). Samples from untreated areas at the corresponding time points (UA1, UA2, UA3, respectively) were collected as controls. A total of 1862 (TA1 vs. UA1), 961 (TA2 vs. UA2), 1370 (TA3 vs. UA3), 3305 (TA2 vs. TA1), 2625 (TA3 vs. TA1), 2899 (TA3 vs. TA2), 782 (UA2 vs. UA1), 4443 (UA3 vs. UA1) and 4031 (UA3 vs. UA2) genes were differentially expressed (DEGs). Functional enrichment analysis showed that DEGs were significantly enriched for secondary metabolic processes, signal transduction and transcriptional regulation processes. Most of the genes involved in lignin biosynthesis were more abundant in the TA groups, which included phenylalanine ammonia-lyase, 4-coumarate CoA ligase, cinnamate 4-hydroxylase, caffeoyl-CoA O-methyltransferase and cinnamoyl-CoA reductase. DEGs involved in sesquiterpene biosynthesis were also identified. Hydroxymethylglutaryl-CoA synthase, 3-hydroxy-3-methylglutaryl-coenzyme A reductase, phosphomevalonate kinase and terpene synthase genes were significantly increased in the TA groups, promoting sesquiterpene biosynthesis in the wounded xylem tissues. The TF-gene transcriptomic networks suggested that MYB DNA-binding, NAM, WRKY, HLH and AP2 TFs co-expressed with genes related to lignin and sesquiterpene synthesis, indicating their critical regulatory roles in the biosynthesis of these compounds. Overall, our study reveals a dynamic transcriptional response of A. sinensis to mechanical wounding, provides a resource for identifying candidate genes for molecular breeding of agarwood quality, and sheds light on the molecular mechanisms of agarwood formation in A. sinensis.

PnMYB4 negatively modulates saponin biosynthesis in Panax notoginseng through interplay with PnMYB1

Saponins are the main triterpenoid ingredients from Panax notoginseng, a well-known Chinese medicine, and are important sources for producing drugs to prevent and treat cerebrovascular and cardiovascular diseases. However, the transcriptional regulatory network of saponin biosynthesis in P. notoginseng is largely unknown. In the present study we demonstrated that one R2R3-MYB transcription factor, designated PnMYB4, acts as a repressor of saponin accumulation. Suppression of PnMYB4 in P. notoginseng calli significantly increased the saponin content and the expression level of saponin biosynthetic genes. PnMYB4 directly bound to the promoters of key saponin biosynthetic genes, including PnSS, PnSE, and PnDS, to repress saponin accumulation. PnMYB4 and the activator PnMYB1 could interacted with PnbHLH, which is a positive regulator of saponin biosynthesis, to modulate the biosynthesis of saponin. PnMYB4 competed with PnMYB1 for binding to PnbHLH, repressing activation of the promoters of saponin structural genes induced by the PnMYB1-PnbHLH complex. Our study reveals that a complex regulatory module of saponin biosynthesis is associated with positive and negative MYB transcriptional regulators and provides a theoretical basis for improving the content of saponins and efficacy of P. notoginseng.

Identification and functional analysis of ZmDLS associated with the response to biotic stress in maize

Terpenes are the main class of secondary metabolites produced in response to pest and germ attacks. In maize (Zea mays L.), they are the essential components of the herbivore-induced plant volatile mixture, which functioned as a direct or indirect defense against pest and germ attacks. In this study, 43 maize terpene synthase gene (ZmTPS) family members were systematically identified and analyzed through the whole genomes of maize. Nine genes, including Zm00001d032230, Zm00001d045054, Zm00001d024486, Zm00001d004279, Zm00001d002351, Zm00001d002350, Zm00001d053916, Zm00001d015053, and Zm00001d015054, were isolated for their differential expression pattern in leaves after corn borer (Ostrinia nubilalis) bite. Additionally, six genes (Zm00001d045054, Zm00001d024486, Zm00001d002351, Zm00001d002350, Zm00001d015053, and Zm00001d015054) were significantly upregulated in response to corn borer bite. Among them, Zm00001d045054 was cloned. Heterologous expression and enzyme activity assays revealed that Zm00001d045054 functioned as d-limonene synthase. It was renamed ZmDLS. Further analysis demonstrated that its expression was upregulated in response to corn borer bites and Fusarium graminearum attacks. The mutant of ZmDLS downregulated the expressions of Zm00001d024486, Zm00001d002351, Zm00001d002350, Zm00001d015053, and Zm00001d015054. It was more attractive to corn borer bites and more susceptible to F. graminearum infection. The yeast one-hybrid assay and dual-luciferase assay showed that ZmMYB76 and ZmMYB101 could upregulate the expression of ZmDLS by binding to the promoter region. This study may provide a theoretical basis for the functional analysis and transcriptional regulation of terpene synthase genes in crops.

Trichome-Specific Analysis and Weighted Gene Co-Expression Correlation Network Analysis (WGCNA) Reveal Potential Regulation Mechanism of Artemisinin Biosynthesis in Artemisia annua

Trichomes are attractive cells for terpenoid biosynthesis and accumulation in Artemisia annua. However, the molecular process underlying the trichome of A. annua is not yet fully elucidated. In this study, an analysis of multi-tissue transcriptome data was performed to examine trichome-specific expression patterns. A total of 6646 genes were screened and highly expressed in trichomes, including artemisinin biosynthetic genes such as amorpha-4,11-diene synthase (ADS) and cytochrome P450 monooxygenase (CYP71AV1). Mapman and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that trichome-specific genes were mainly enriched in lipid metabolism and terpenoid metabolism. These trichome-specific genes were analyzed by a weighted gene co-expression network analysis (WGCNA), and the blue module linked to terpenoid backbone biosynthesis was determined. Hub genes correlated with the artemisinin biosynthetic genes were selected based on TOM value. ORA, Benzoate carboxyl methyltransferase (BAMT), Lysine histidine transporter-like 8 (AATL1), Ubiquitin-like protease 1 (Ulp1) and TUBBY were revealed as key hub genes induced by methyl jasmonate (MeJA) for regulating artemisinin biosynthesis. In summary, the identified trichome-specific genes, modules, pathways and hub genes provide clues and shed light on the potential regulatory mechanisms of artemisinin biosynthesis in trichomes in A. annua.

A long-read and short-read transcriptomics approach provides the first high-quality reference transcriptome and genome annotation for Pseudotsuga menziesii (Douglas-fir)

Douglas-fir (Pseudotsuga menziesii) is native to western North America. It grows in a wide range of environmental conditions and is an important timber tree. Although there are several studies on the gene expression responses of Douglas-fir to abiotic cues, the absence of high-quality transcriptome and genome data is a barrier to further investigation. Like for most conifers, the available transcriptome and genome reference dataset for Douglas-fir remains fragmented and requires refinement. We aimed to generate a highly accurate, and complete reference transcriptome and genome annotation. We deep-sequenced the transcriptome of Douglas-fir needles from seedlings that were grown under nonstress control conditions or a combination of heat and drought stress conditions using long-read (LR) and short-read (SR) sequencing platforms. We used 2 computational approaches, namely de novo and genome-guided LR transcriptome assembly. Using the LR de novo assembly, we identified 1.3X more high-quality transcripts, 1.85X more "complete" genes, and 2.7X more functionally annotated genes compared to the genome-guided assembly approach. We predicted 666 long noncoding RNAs and 12,778 unique protein-coding transcripts including 2,016 putative transcription factors. We leveraged the LR de novo assembled transcriptome with paired-end SR and a published single-end SR transcriptome to generate an improved genome annotation. This was conducted with BRAKER2 and refined based on functional annotation, repetitive content, and transcriptome alignment. This high-quality genome annotation has 51,419 unique gene models derived from 322,631 initial predictions. Overall, our informatics approach provides a new reference Douglas-fir transcriptome assembly and genome annotation with considerably improved completeness and functional annotation.

Lend Me Your EARs: A Systematic Review of the Broad Functions of EAR Motif-Containing Transcriptional Repressors in Plants

The ethylene-responsive element binding factor-associated amphiphilic repression (EAR) motif, defined by the consensus sequence patterns LxLxL or DLNx(x)P, is found in a diverse range of plant species. It is the most predominant form of active transcriptional repression motif identified so far in plants. Despite its small size (5 to 6 amino acids), the EAR motif is primarily involved in the negative regulation of developmental, physiological and metabolic functions in response to abiotic and biotic stresses. Through an extensive literature review, we identified 119 genes belonging to 23 different plant species that contain an EAR motif and function as negative regulators of gene expression in various biological processes, including plant growth and morphology, metabolism and homeostasis, abiotic stress response, biotic stress response, hormonal pathways and signalling, fertility, and ripening. Positive gene regulation and transcriptional activation are studied extensively, but there remains much more to be discovered about negative gene regulation and the role it plays in plant development, health, and reproduction. This review aims to fill the knowledge gap and provide insights into the role that the EAR motif plays in negative gene regulation, and provoke further research on other protein motifs specific to repressors.

Discovery of specific catalytic activity toward IAA/FA by LaSABATHs based on genome-wide phylogenetic and enzymatic analysis of SABATH gene family from Larix kaempferi

The SABATH methyltransferases catalyze methylation of small-molecule metabolites, which participate in plant growth, development and defense response. Given lack of genome-wide studies on gymnosperms SABATH family, the formation and functional differentiation mechanism of the Larix kaempferi SABATH gene family was systematically and exhaustively explored by analyzing gene sequence characteristics, phylogenetic relationship, expression pattern, and enzyme activities. Phylogenetic analysis showed that 247 SABATH genes from 14 land plants were divided into 4 clades, and lineage-specific gene duplication events were important factors that contributed to the evolution of the SABATH gene family in gymnosperms and angiosperms. Substrate specificity analysis of 18 Larix SABATH proteins showed that LaSABATHs could catalyze O-methylation of indole-3-acetic acid (IAA) and farnesic acid (FA), N-methylation of theobromine, and S-methylation of thiobenzoic acid. Furthermore, only LaSABATH2 and LaSABATH29 could catalyze O-methylation of FA, and only LaSABATH30 could catalyze O-methylation of IAA. Homology modeling and molecular docking studies showed the hydrogen bond formed between the His188 of LaSABATH30 and IAA and the noticeable hydrophobic IAA-binding pocket may be helpful for IAA methylation. In this study, identification of proteins with significant specific catalytic activity toward FA and IAA provided high-quality candidate genes for forest genetics and breeding.

Genome-wide characterization and functional identification of MYB genes in Malus sieversii infected by Valsa mali

Among the most important transcription factors in plants, the v-myb avian myeloblastosis viral oncogene homolog (MYB) regulates the expression network of response genes under stresses such as fungal infection. In China, the canker disease Valsa mali threatens the survival of Malus sieversii, an ancestor of cultivated apples. Using the M. sieversii genome, we identified 457 MsMYB and 128 R2R3-MsMYB genes that were randomly distributed across 17 chromosomes. Based on protein sequence and structure, the R2R3-MsMYB genes were phylogenetically divided into 29 categories, and 26 conserved motifs were identified. We further predicted cis-elements in the 2000-kb promoter region of R2R3-MsMYBs based on the genome. Transcriptome analysis of M. sieversii under V. mali infection showed that 27 R2R3-MsMYBs were significantly differentially expressed, indicating their key role in the response to V. mali infection. Using transient transformation, MsMYB14, MsMYB24, MsMYB39, MsMYB78, and MsMYB108, which were strongly induced by V. mali infection, were functionally identified. Among the five MsMYBs, MsMYB14 and MsMYB78 were both important in enhancing resistance to diseases, whereas MsMYB24 inhibited resistance. Based on the results of this study, we gained a better understanding of the MsMYB transcription factor family and laid the foundation for a future research program on disease prevention strategies in M. sieversii.

Plagiodera versicolora feeding induces systemic and sexually differential defense responses in poplars

Dioecious plants have evolved effective defense strategies to deal with various biotic and abiotic stresses. However, little is known regarding sexual differences in their defense against herbivores. In this study, we investigated the mechanism of systemic defense responses in male and female Populus cathayana attacked by Plagiodera versicolora Laicharting. The results revealed that P. cathayana exhibits sexually differential responses to a defoliator. The percentage of damaged leaf area was greater in males than in females. Furthermore, the observed saccharide changes imply that males and females exhibit different response times to defoliators. The contents of flavonoids and anthocyanins were significantly increased in both sexes but were higher in females. Specifically, the jasmonic acid (JA) pathway plays an important role. Expression of pest-related genes further revealed that hormones induce changes in downstream genes and metabolites, and upregulation of JA ZIM-domain (JAZ) and CORONATINE INSENSITIVE 1 (COI1) was more significant in females. In the undamaged adjacent leaves, metabolite and gene changes displayed similar patterns to the damaged local leaves, but levels of JA, JAZ1, and COI1 were higher in females. Therefore, our data confirmed that plants initiate the JA pathway to defend against herbivores, that there is systematic signal transduction, and that this ability is stronger in females than in males. This study provides new insights into the resistance of dioecious plants to herbivory and adds a new theoretical basis for the systemic signal transduction of plants in response to biotic stress.

The SmMYB36-SmERF6/SmERF115 module regulates the biosynthesis of tanshinones and phenolic acids in salvia miltiorrhiza hairy roots

Tanshinone and phenolic acids are the most important active substances of Salvia miltiorrhiza, and the insight into their transcriptional regulatory mechanisms is an essential process to increase their content in vivo. SmMYB36 has been found to have important regulatory functions in the synthesis of tanshinone and phenolic acid; paradoxically, its mechanism of action in S. miltiorrhiza is not clear. Here, we demonstrated that SmMYB36 functions as a promoter of tanshinones accumulation and a suppressor of phenolic acids through the generation of SmMYB36 overexpressed and chimeric SmMYB36-SRDX (EAR repressive domain) repressor hairy roots in combination with transcriptomic-metabolomic analysis. SmMYB36 directly down-regulate the key enzyme gene of primary metabolism, SmGAPC, up-regulate the tanshinones biosynthesis branch genes SmDXS2, SmGGPPS1, SmCPS1 and down-regulate the phenolic acids biosynthesis branch enzyme gene, SmRAS. Meanwhile, SmERF6, a positive regulator of tanshinone synthesis activating SmCPS1, was up-regulated and SmERF115, a positive regulator of phenolic acid biosynthesis activating SmRAS, was down-regulated. Furthermore, the seven acidic amino acids at the C-terminus of SmMYB36 are required for both self-activating domain and activation of target gene expression. As a consequence, this study contributes to reveal the potential relevance of transcription factors synergistically regulating the biosynthesis of tanshinone and phenolic acid.

Structural and interactions analysis of a transcription factor PnMYB2 in Panax notoginseng

The main active ingredients of the traditional Chinese medicinal plant, Panax notoginseng, are the Panax notoginseng saponins (PNS). They can be synthesized via the mevalonate pathway; PnSS and PnSE1 are the key rate-limiting enzymes in this pathway. In this study, an interaction between PnMYB2 and the key enzymes was identified and characterized from the P. notoginseng cDNA library using the Y1H technique. Subsequently, X-α-gal color reaction confirmed the interaction between PnMYB2 and the upstream sequences of PnSS and PnSE1 promoters. Full-length cDNA sequence of PnMYB2 was isolated and characterized. PnMYB2 has an open reading frame of 864 bp, encoding 287 amino acids. 3D structural analysis of PnMYB2 indicated that its structure was similar to that of the template. Phylogenetic analysis revealed that PnMYB2 and PgMYB2 are highly homologous and belong to the R2R3 MYB transcription factor (TF). Subcellular localization analysis showed that PnMYB2 was localized in the nucleus. The recombinant protein PnMYB2 was successfully obtained through prokaryotic expression and was confirmed to be an inclusion body protein. Furthermore, electrophoretic mobility shift assay (EMSA) experiments demonstrated that PnMYB2 specifically binds to MYB core and AC-rich elements. This study provides a theoretical basis for transcriptional regulation of saponin biosynthesis in P. notoginseng.

Identification and characterization of a novel gene involved in glandular trichome development in Nepeta tenuifolia

Nepeta tenuifolia is a medicinal plant rich in terpenoids and flavonoids with antiviral, immunoregulatory, and anti-inflammatory activities. The peltate glandular trichome (PGT) is a multicellular structure considered to be the primary storage organ for monoterpenes; it may serve as an ideal model for studying cell differentiation and the development of glandular trichomes (GTs). The genes that regulate the development of GTs have not yet been well studied. In this study, we identified NtMIXTA1, a GT development-associated gene from the R2R3 MYB SBG9 family. NtMIXTA1 overexpression in tobacco resulted in the production of longer and denser GTs. Virus-induced gene silencing of NtMIXTA1 resulted in lower PGT density, a significant reduction in monoterpene concentration, and the decreased expression of genes related to monoterpene biosynthesis. Comparative transcriptome and widely targeted metabolic analyses revealed that silencing NtMIXTA1 significantly influenced the expression of genes, and the production of metabolites involved in the biosynthesis of terpenoids, flavonoids, and lipids. This study provides a solid foundation describing a mechanism underlying the regulation of GT development. In addition, this study further deepens our understanding of the regulatory networks involved in GT development and GT development-associated metabolite flux, as well as provides valuable reference data for studying plants with a high medicinal value without genetic transformation.

Analysis of Transcriptome and Terpene Constituents of Scots Pine Genotypes Inherently Resistant or Susceptible to Heterobasidion annosum

Root and stem rot caused by Heterobasidion annosum is a severe problem in boreal Scots pine. Dissecting the features of disease resistance is generally an essential step in resistance breeding in plants and forest trees. In this study, we explored inherent resistance factors of Scots pine against H. annosum. A total of 236 families consisting of 85 full-sib (FS), 35 half-sib population mix (HSpm), and 116 half-sib (HS) families of Scots pine seedlings were inoculated with a H. annosum isolate. We sampled needle tissues before inoculation for terpene measurements and RNA sequencing. Based on the lesion area, the extremes of 12 resistant and 12 susceptible families were selected for further analyses. Necrotic lesions resulting from fungal infection were in a weak to moderate relationship with the plant height. Monoterpenes were the principal terpene compounds observed in Scots pine seedlings. Concentrations of 3-carene were significantly higher in pine genotypes inherently resistant compared with susceptible seedlings. By contrast, susceptible genotypes had significantly higher proportions of α-pinene. Gene ontology analysis of differential expressed transcripts (DETs) revealed that response to biotic factors was enriched in resistant seedlings. Functional characterization of individual DETs revealed that higher expression of transcripts involved in response to abiotic stress was common in susceptible genotypes. This observation was supported by the annotation of hub genes in a key module that was significantly correlated with the lesion trait through weighted gene co-expression network analysis (WGCNA) of 16 HS and HSpm samples. These findings contribute to our understanding of constitutive resistance factors of Scots pine against Heterobasidion root and stem rot diseases.

Transcriptional cascades in the regulation of 2-AP biosynthesis under Zn supply in fragrant rice

Transcription factors (TFs) regulate gene expression to control certain genetic programs, such as growth and development, phytohormone regulation, and environmental stresses. 2-acetyl-1-pyrroline (2-AP) is the key element involved in aroma biosynthesis pathway, and the application of micronutrients can increase the 2-AP levels. However, little is known about the micronutrient-induced TFs involved in 2-AP biosynthesis. Here, we identify a number of TF families in two fragrant rice varieties, "Meixiangzhan-2" (M) and "Xiangyaxiangzhan" (X), in response to Zinc (Zn) application through transcriptomic analysis. A total of ~678 TFs were identified and grouped into 26 TF families, each of which was found to be involved in numerous signaling pathways. The WRKY TF family was found to be the most abundant, followed by bHLH and MYB. Furthermore, members of the WRKY, bHLH, MYB, ERF, HSF, MADS-box, NFY, and AP2 TF families were significantly upregulated and may be involved in the transcriptional regulation of aroma biosynthesis. In brief, this study enhances our understanding of the molecular mechanism of 2-AP biosynthesis and highlights the key TFs potentially involved in the production of aroma in fragrant rice.

Molecular and Chemical Screening for Inherent Disease Resistance Factors of Norway Spruce (Picea abies) Clones Against Conifer Stem Rot Pathogen Heterobasidion parviporum

Root and stem rot of conifer trees caused by Heterobasidion annosum species complex leads to huge economic losses in Europe, yet not much is known about the molecular and chemical basis for host resistance. To identify inherent chemical or molecular markers in clones found to be either resistant or susceptible, we sampled needle tissues of all the clones before pathogen inoculation. We conducted a short-term resistance screening by using the pathogen H. parviporum to inoculate 70 Norway spruce clones. Based on lesion size, subsets of highly susceptible and resistant clones were further analyzed. Terpene detection and RNA sequencing were performed to explore inherent variations in genotypes differing in resistance to pathogenic challenge at chemical and transcriptional levels. A negative correlation emerged between resistance and growth. Terpene profiles of resistant clones showed higher content of monoterpenes and sesquiterpenes, with concomitant increased transcript abundance of genes involved in the terpenoid pathway. A set of upregulated genes relevant to flavonoid biosynthesis was observed in resistant genotypes, whereas higher transcripts of lignin biosynthetic genes were prevalent in susceptible clones. Genes involved in flavonoid and lignin biosynthesis as well as terpene content may have a role in facilitating resistance of Norway spruce against H. parviporum. Our results provide strong support on the feasibility of sampling needle tissues before pathogen inoculation, and the approach could be of value for large-scale screening of novel biomarkers for durable resistance. The additional insights could form a basis for further research on resistance screening in this pathosystem.

The Potential of Cold Plasma and Electromagnetic Field as Stimulators of Natural Sweeteners Biosynthesis in Stevia rebaudiana Bertoni

Stevioside (Stev) and rebaudioside A (RebA) are the most abundant steviol glycosides (SGs) responsible for the sweetness of Stevia rabaudiana Bertoni. As compared to Stev, RebA has a higher sweetening potency, better taste and therefore is the most preferred component of the stevia leaf extracts. The aim of this study was to determine the effect of pre-sowing seed treatment with abiotic stressors cold plasma (CP) and electromagnetic field (EMF) on the amount and ratio of RebA and Stev in the leaves of stevia. Additionally, the effect on total phenolic content, flavonoid content and antioxidant activity was investigated. Seeds were treated 5 and 7 min with cold plasma (CP5 and CP7 groups) and 10 min with electromagnetic field (EMF10 group) six days before sowing. The germination tests in vitro demonstrated that all treatments slightly increased germination rate and percentage. HPLC analysis revealed that CP and EMF had strong stimulating effect on SGs accumulation. All treatments increased RebA concentration approximately 1.6-fold; however, the ratio of RebA/Stev decreased from 8.5 in the control to 1.9, 2.5 and 1.1 in CP5, CP7 and EMF10 groups respectively, since the concentration of Stev increased more than RebA, 7.1, 4.6 and 11.0-fold, respectively, compared to control. However, treatments had opposite effect on total phenolic content, flavonoid content, and antioxidant activity. We have demonstrated for the first time that short time pre-sowing treatment of stevia seeds with CP and EMF can be a powerful tool for the enhancement of biosynthesis of RebA and Stev, however it can have negative impact on the content of other secondary metabolites.

Full-length transcriptomic identification of R2R3-MYB family genes related to secondary cell wall development in Cunninghamia lanceolata (Chinese fir)

BACKGROUND

R2R3-MYB is a class of transcription factor crucial in regulating secondary cell wall development during wood formation. The regulation of wood formation in gymnosperm has been understudied due to its large genome size. Using Single-Molecule Real-Time sequencing, we obtained full-length transcriptomic libraries from the developmental stem of Cunninghamia lanceolata, a perennial conifer known as Chinese fir. The R2R3-MYB of C. lanceolata (hereafter named as ClMYB) associated with secondary wall development were identified based on phylogenetic analysis, expression studies and functional study on transgenic line.

RESULTS

The evolutionary relationship of 52 ClMYBs with those from Arabidopsis thaliana, Eucalyptus grandis, Populus trichocarpa, Oryza sativa, two gymnosperm species, Pinus taeda, and Picea glauca were established by neighbour-joining phylogenetic analysis. A large number of ClMYBs resided in the woody-expanded subgroups that predominated with the members from woody dicots. In contrast, the woody-preferential subgroup strictly carrying the members of woody dicots contained only one candidate. The results suggest that the woody-expanded subgroup emerges before the gymnosperm/angiosperm split, while most of the woody-preferential subgroups are likely lineage-specific to woody dicots. Nine candidates shared the same subgroups with the A. thaliana orthologs, with known function in regulating secondary wall development. Gene expression analysis inferred that ClMYB1/2/3/4/5/26/27/49/51 might participate in secondary wall development, among which ClMYB1/2/5/26/27/49 were significantly upregulated in the highly lignified compression wood region, reinforcing their regulatory role associated with secondary wall development. ClMYB1 was experimentally proven a transcriptional activator that localised in the nucleus. The overexpression of ClMYB1 in Nicotiana benthamiana resulted in an increased lignin deposition in the stems. The members of subgroup S4, ClMYB3/4/5 shared the ERF-associated amphiphilic repression motif with AtMYB4, which is known to repress the metabolism of phenylpropanoid derived compounds. They also carried a core motif specific to gymnosperm lineage, suggesting divergence of the regulatory process compared to the angiosperms.

CONCLUSIONS

This work will enrich the collection of full-length gymnosperm-specific R2R3-MYBs related to stem development and contribute to understanding their evolutionary relationship with angiosperm species.

Characterization of transcription factor MYB59 and expression profiling in response to low K+ and NO3- in indica rice (Oryza sativa L.)

Background

Nitrogen and potassium are crucial supplements for plant development and growth. Plants can detect potassium and nitrate ions in soils and in like way, they modify root-to-shoot transport of these ions to adjust the conveyance among roots and shoots. Transcription factor MYB59 plays essential roles in numerous physiological processes inclusive of hormone response, abiotic stress tolerance, plant development, and metabolic regulation. In this study, we retrieved 56 MYB59 proteins from different plant species. Multiple sequence alignment, phylogenetic tree, conserved motif, chromosomal localization, and cis-regulatory elements of the retrieved sequences were analyzed. Gene structure, protein 3D structure, and DNA binding of OsMYB59 indica were also predicted. Finally, we characterized OsMYB59 and its function under low K⁺/NO₃⁻ conditions in Oryza sativa subsp. indica.

Results

Data analysis showed that MYB59s from various groups separated in terms of conserved functional domains and gene structure, where members of genus Oryza clustered together. Plants showed reduced height and yellowish appearance when grown on K⁺ and NO₃⁻ deficient medium. Quantitative real-time PCR uncovered that the OsMYB59 reacted to abiotic stresses where its expression was increased in BRRI dhan56 but decreased in other varieties on K⁺ deficient medium. In addition, OsMYB59 transcript level increased on NO₃⁻ deficient medium.

Conclusions

Our results can help to explain the biological functions of indica rice MYB59 protein and gave a theoretical premise to additionally describe its biological roles in response to abiotic stresses particularly drought.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-021-00248-6.

Integrated Four Comparative-Omics Reveals the Mechanism of the Terpenoid Biosynthesis in Two Different Overwintering Cryptomeria fortunei Phenotypes

Chinese cedar (Cryptomeria fortunei) is a tree species with important ornamental, medicinal, and economic value. Terpenoids extracted from the essential oil of C. fortunei needles have been considered valuable ingredients in the pharmaceutical and cosmetic industries. However, the possible gene regulation mechanisms that limit terpenoid biosynthesis in this genus are poorly understood. Here, we adopted integrated metabolome analysis, transcriptome, small-RNA (sRNA), and degradome sequencing to analyze the differences in terpenoid regulatory mechanisms in two different overwintering C. fortunei phenotypes (wild-type and an evergreen mutant). A total of 1447/6219 differentially synthesized metabolites (DSMs)/unigenes (DEGs) were detected through metabolome/transcriptome analyses, and these DSMs/DEGs were significantly enriched in flavonoid and diterpenoid biosynthesis pathways. In C. fortunei needles, 587 microRNAs (miRNAs), including 67 differentially expressed miRNAs (DERs), were detected. Among them, 8346 targets of 571 miRNAs were predicted using degradome data, and a 72-miRNA-target regulatory network involved in the metabolism of terpenoids and polyketides was constructed. Forty-one targets were further confirmed to be involved in terpenoid backbone and diterpenoid biosynthesis, and target analyses revealed that two miRNAs (i.e., aly-miR168a-5p and aof-miR396a) may be related to the different phenotypes and to differential regulation of diterpenoid biosynthesis. Overall, these results reveal that C. fortunei plants with the evergreen mutation maintain high terpenoid levels in winter through miRNA-target regulation, which provides a valuable resource for essential oil-related bioengineering research.

Functional Characterization of Hedychium coronarium J. Koenig MYB132 Confers the Potential Role in Floral Aroma Synthesis

The R2R3-MYB transcription factors (TFs) play several key roles in numerous plant biological processes. Hedychium coronarium is an important ornamental plant well-known for its elegant flower shape and abundant aroma type. The floral aroma of H. coronarium is due to the presence of a large amount of terpenes and benzenoids. However, less is known about the role of R2R3-MYB TFs in the regulatory mechanism of floral aroma production in this breed. Herein, we isolate and functionally characterize the R2R3-MYB TF HcMYB132, which is potentially involved in regulating floral aroma synthesis. Sequence alignment analysis revealed that it includes a nuclear localization signal NLS(s) and a 2R, 3R motif signature in the sequences. A subcellular localization assay revealed that HcMYB132 protein localizes to the nucleus. Real-time qPCR assays showed that HcMYB132 is specifically expressed in flowers and its expression pattern correlates with the emission of floral volatile compounds. In HcMYB132-silenced flowers, the levels of floral volatile compounds were significantly reduced, and the expression of key structural volatile synthesis genes was downregulated compared to control. Collectively, these results suggest that HcMYB132 might play a significant role in the regulation of terpenoid biosynthesis in H. coronarium.

Melatonin alters the secondary metabolite profile of grape berry skin by promoting VvMYB14-mediated ethylene biosynthesis

The interplay between melatonin and ethylene in the regulation of fruit metabolism and the underlying molecular mechanism of this interplay remain largely unclear. Here, widely targeted metabolomics analysis revealed a total of 464 metabolites present in berry skin. Among them, 27 significantly differentially accumulated metabolites (DAMs) were produced in response to melatonin treatment in the presence or absence of 1-MCP. Most of the DAMs were secondary metabolites, including flavonoids, phenolic acids, stilbenes, and flavonols. Additionally, the accumulation of 25 DAMs was regulated by melatonin via ethylene. RNA-seq analysis indicated that melatonin primarily regulated the pathways of plant hormone signal transduction and secondary metabolite biosynthesis via ethylene. Gene-metabolite association analysis showed that melatonin regulated the expression of the VvSTS1, VvF3H, VvLAR2, and VvDFR genes, suggesting that these genes may play key roles in regulating secondary metabolites in the skin; additionally, VvMYB14 and VvACS1 were suggested to be involved in the regulation of secondary metabolites. Further experiments revealed that melatonin induced the expression of VvMYB14 and that VvMYB14 increased ethylene production by transcriptionally activating VvACS1, thereby affecting the accumulation of secondary metabolites. Collectively, melatonin promotes ethylene biosynthesis and alters secondary metabolite accumulation through the regulation of VvACS1 by VvMYB14.

Genome-Wide Analysis Reveals the Potential Role of MYB Transcription Factors in Floral Scent Formation in Hedychium coronarium

The MYB gene family is one of the largest groups of transcription factors (TFs) playing diverse roles in several biological processes. Hedychium coronarium (white ginger lily) is a renowned ornamental plant both in tropical and subtropical regions due to its flower shape and strong floral scent mainly composed of terpenes and benzenoids. However, there is no information available regarding the role of the MYB gene family in H. coronarium. In the current study, the MYB gene family was identified and extensively analyzed. The identified 253 HcMYB genes were unevenly mapped on 17 chromosomes at a different density. Promoter sequence analysis showed numerous phytohormones related to cis-regulatory elements. The majority of HcMYB genes contain two to three introns and motif composition analysis showed their functional conservation. Phylogenetic analysis revealed that HcMYBs could be classified into 15 distinct clades, and the segmental duplication events played an essential role in the expansion of the HcMYB gene family. Tissue-specific expression patterns of HcMYB genes displayed spatial and temporal expression. Furthermore, seven HcMYB (HcMYB7/8/75/79/145/238/248) were selected for further investigation. Through RT-qPCR, the response of candidates HcMYB genes toward jasmonic acid methyl ester (MeJA), abscisic acid (ABA), ethylene, and auxin was examined. Yeast one-hybrid (Y1H) assays revealed that candidate genes directly bind to the promoter of bottom structural volatile synthesis genes (HcTPS1, HcTPS3, HcTPS10, and HcBSMT2). Moreover, yeast two-hybrid (Y2H) assay showed that HcMYB7/8/75/145/248 interact with HcJAZ1 protein. In HcMYB7/8/79/145/248-silenced flowers, the floral volatile contents were decreased and downregulated the expression of key structural genes, suggesting that these genes might play crucial roles in floral scent formation in H. coronarium by regulating the expression of floral scent biosynthesis genes. Collectively, these findings indicate that HcMYB genes might be involved in the regulatory mechanism of terpenoids and benzenoid biosynthesis in H. coronarium.

Auxin-Responsive R2R3-MYB Transcription Factors HcMYB1 and HcMYB2 Activate Volatile Biosynthesis in Hedychium coronarium Flowers

Auxin, an important plant hormone, induces the biosynthesis of various secondary metabolites by modulating the expression of auxin-responsive genes. In the ornamental plant Hedychium coronarium, linalool and methyl benzoate are biosynthesized by the terpene synthase (TPS) HcTPS5 and the benzoic/salicylic acid methyltransferase (BSMT) HcBSMT2, respectively. However, the transcriptional regulation of this process remains unclear. Here, we identified and functionally characterized the R2R3-MYB transcription factors HcMYB1 and HcMYB2 in regulating the biosynthesis of these floral aroma compounds. HcMYB1 and HcMYB2 are specifically expressed in flowers, their expression is correlated with the emission of volatile compounds in flowers, and is induced by auxin. Moreover, HcMYB1 and HcMYB2 interact with the HcBSMT2 promoter region. HcMYB2 activates the expression of the linalool synthase gene HcTPS5. In flowers with HcMYB1 or HcMYB2 silenced, the levels of floral scent compounds were significantly reduced, and HcBSMT2 and HcTPS5 were downregulated compared with the wild type. Moreover, HcMYB1 form protein-protein interaction with key scent-related HcIAA4 protein to regulate floral aroma production. Taken together, these results indicate that HcMYB1 and HcMYB2 play crucial roles in regulating the formation of scent compounds in Hedychium coronarium (H. coronarium) flowers in response to auxin signaling.

Identification of Transcription Factor Genes and Functional Characterization of PlMYB1 From Pueraria lobata

Kudzu, Pueraria lobata, is a traditional Chinese food and medicinal herb that has been commonly used since ancient times. Kudzu roots are rich sources of isoflavonoids, e.g., puerarin, with beneficial effects on human health. To gain global information on the isoflavonoid biosynthetic regulation network in kudzu, de novo transcriptome sequencings were performed using two genotypes of kudzu with and without puerarin accumulation in roots. RNAseq data showed that the genes of the isoflavonoid biosynthetic pathway were significantly represented in the upregulated genes in the kudzu with puerarin. To discover regulatory genes, 105, 112, and 143 genes encoding MYB, bHLH, and WD40 transcription regulators were identified and classified, respectively. Among them, three MYB, four bHLHs, and one WD40 gene were found to be highly identical to their orthologs involved in flavonoid biosynthesis in other plants. Notably, the expression profiles of PlMYB1, PlHLH3-4, and PlWD40-1 genes were closely correlated with isoflavonoid accumulation profiles in different tissues and cell cultures of kudzu. Over-expression of PlMYB1 in Arabidopsis thaliana significantly increased the accumulation of anthocyanins in leaves and proanthocyanidins in seeds, by activating AtDFR, AtANR, and AtANS genes. Our study provided valuable comparative transcriptome information for further identification of regulatory or structural genes involved in the isoflavonoid pathway in P. lobata, as well as for bioengineering of bioactive isoflavonoid compounds.

Combining QTL Mapping and Transcriptomics to Decipher the Genetic Architecture of Phenolic Compounds Metabolism in the Conifer White Spruce

Conifer forests worldwide are becoming increasingly vulnerable to the effects of climate change. Although the production of phenolic compounds (PCs) has been shown to be modulated by biotic and abiotic stresses, the genetic basis underlying the variation in their constitutive production level remains poorly documented in conifers. We used QTL mapping and RNA-Seq to explore the complex polygenic network underlying the constitutive production of PCs in a white spruce (Picea glauca) full-sib family for 2 years. QTL detection was performed for nine PCs and differentially expressed genes (DEGs) were identified between individuals with high and low PC contents for five PCs exhibiting stable QTLs across time. A total of 17 QTLs were detected for eight metabolites, including one major QTL explaining up to 91.3% of the neolignan-2 variance. The RNA-Seq analysis highlighted 50 DEGs associated with phenylpropanoid biosynthesis, several key transcription factors, and a subset of 137 genes showing opposite expression patterns in individuals with high levels of the flavonoids gallocatechin and taxifolin glucoside. A total of 19 DEGs co-localized with QTLs. Our findings represent a significant step toward resolving the genomic architecture of PC production in spruce and facilitate the functional characterization of genes and transcriptional networks responsible for differences in constitutive production of PCs in conifers.

Metabolite and Transcriptome Profiling Analysis Revealed That Melatonin Positively Regulates Floral Scent Production in Hedychium coronarium

Melatonin is a pleiotropic molecule that regulates a variety of developmental processes. Floral volatiles are important features of flowers that facilitate flower-visitor interactions by attracting pollinators, structure flower-visitor communities, and play defensive roles against plant and flower antagonists. Aside from their role in plants, floral volatiles are an essential ingredient in cosmetics, perfumes, pharmaceuticals, and flavorings. Herein, integrated metabolomic and transcriptomic approaches were carried out to analyze the changes triggered by melatonin exposure during the Hedychium coronarium flower development stages. Quantitative analysis of the volatiles of H. coronarium flowers revealed that volatile organic compound emission was significantly enhanced after melatonin exposure during the half bloom (HS), full bloom (FB) and fade stage (FS). Under the melatonin treatment, the emission of volatile contents was highest during the full bloom stage of the flower. Variable importance in projection (VIP) analysis and partial least-squares discriminant analysis (PLS-DA) identified 15 volatile compounds with VIP > 1 that were prominently altered by the melatonin treatments. According to the transcriptome sequencing data of the HS, FB, and FS of the flowers, 1,372, 1,510, and 1,488 differentially expressed genes were identified between CK-HS and 100MT-HS, CK-FB and 100MT-FB, and CK-FS and 100MT-FS, respectively. Among the significant differentially expressed genes (DEGs), 76 were significantly upregulated and directly involved in the floral scent biosynthesis process. In addition, certain volatile organic compounds were substantially linked with various DEGs after combining the metabolome and transcriptome datasets. Moreover, some transcription factors, such as MYB and bHLH, were also significantly upregulated in the comparison, which might be related to the floral aroma mechanism. Our results suggested that melatonin increased floral aroma production in H. coronarium flowers by modifying the expression level of genes involved in the floral scent biosynthesis pathway. These findings serve as a foundation for future research into the molecular mechanisms underlying the dynamic changes in volatile contents induced by melatonin treatment in H. coronarium.

The spatio-temporal biosynthesis of floral flavonols is controlled by differential phylogenetic MYB regulators in Freesia hybrida

Floral flavonols play specific pivotal roles in pollinator attraction, pollen germination and fertility, in addition to other functions in vegetative organs. For many plants, the process of flavonol biosynthesis in late flower development stages and in mature flower tissues is poorly understood, in contrast to early flower development stages. It is thought that this process may be regulated independently of subgroup 7 R2R3 MYB (SG7 MYB) transcription factors. In this study, two FLS genes were shown to be expressed synchronously with the flower development-specific and tissue-specific biosynthesis of flavonols in Freesia hybrida. FhFLS1 contributed to flavonol biosynthesis in early flower buds, toruses and calyxes, and was regulated by four well-known SG7 MYB proteins, designated as FhMYBFs, with at least partial regulatory redundancy. FhFLS2 accounted for flavonols in late developed flowers and in the petals, stamens and pistils, and was targeted directly by non SG7 MYB protein FhMYB21L2. In parallel, AtMYB21 and AtMYB24 also activated AtFLS1, a gene highly expressed in Arabidopsis anthers and pollen, indicating the conserved regulatory roles of MYB21 against FLS genes in these two evolutionarily divergent angiosperm plants. Our results reveal a novel regulatory and synthetic mechanism underlying flavonol biosynthesis in floral organs and tissues which may be exploited to investigate supplementary roles of flavonols in flowers.

Plant Volatile Organic Compounds Evolution: Transcriptional Regulation, Epigenetics and Polyploidy

Volatile organic compounds (VOCs) are emitted by plants as a consequence of their interaction with biotic and abiotic factors, and have a very important role in plant evolution. Floral VOCs are often involved in defense and pollinator attraction. These interactions often change rapidly over time, so a quick response to those changes is required. Epigenetic factors, such as DNA methylation and histone modification, which regulate both genes and transcription factors, might trigger adaptive responses to these evolutionary pressures as well as regulating the rhythmic emission of VOCs through circadian clock regulation. In addition, transgenerational epigenetic effects and whole genome polyploidy could modify the generation of VOCs’ profiles of offspring, contributing to long-term evolutionary shifts. In this article, we review the available knowledge about the mechanisms that may act as epigenetic regulators of the main VOC biosynthetic pathways, and their importance in plant evolution.

Boosting leaf contents of rutin and solanesol in bio-waste of Solanum lycopersicum

In tomato production, the accruing green biomass shows promising potential as source of health-promoting compounds, such as rutin and solanesol, that are of high interest due to their medicinal properties. Naturally, they accumulate in plants growing in suboptimal growing conditions, e.g. influenced by biotic and abiotic stressors. With the aim to evaluate the potential use of tomato residues as source, we analyzed both leaf metabolites during a complete cultivation cycle, while applying single and combined stresses practically realized in greenhouse production. In the late season, contents of both metabolites were significantly enhanced by nutrient deficit in combination with 2 °C colder nights for 4 weeks and prolonged for in total 9 weeks. Particularly, higher solanesol contents were achieved by salt stress and elevated temperature after one week, even stronger when combined with drought. At harvest, stressed plants consist of less green biomass reducing the overall economic potential. However, practicable abiotic stresses should be considered as potential tool to induce the accumulation of beneficial compounds. Extracting profitable metabolites from the green biomass of the model crop tomato supports the overall goal to promote sustainable approaches in horticultural production.

Ectopic Expression of Production of Anthocyanin Pigment 1 (PAP1) Improves the Antioxidant and Anti-Melanogenic Properties of Ginseng (Panax ginseng C.A. Meyer) Hairy Roots

The development of genetically engineered cell cultures has been suggested as a potential approach for the production of target compounds from medicinal plants. In this study, we generated PAP1 (production of anthocyanin pigment 1)-overexpressing ginseng (Panax ginseng C.A. Meyer) hairy roots to improve the production of anthocyanins, as well as the bioactivity (e.g., antioxidant and whitening activities) of ginseng. Based on differentially expressed gene analysis, we found that ectopic expression of PAP1 induced the expression of genes involved in the ‘phenylpropanoid biosynthesis’ (24 genes), and ‘flavonoid biosynthesis’ (17 genes) pathways, resulting in 191- to 341-fold increases in anthocyanin production compared to transgenic control (TC) hairy roots. Additionally, PAP1-overexpressing ginseng hairy roots exhibited an approximately seven-fold higher DPPH-free radical scavenging activity and 10-fold higher ORAC value compared to the TC. In α-melanocyte-stimulating hormone-stimulated B16F10 cells, PAP1-overexpressing ginseng hairy roots strongly inhibited the accumulation of melanin by 50 to 59% compared to mock-control. Furthermore, results obtained by quantitative real-time PCR, western blot, and tyrosinase inhibition assay suggested that the anti-melanogenic activity of PAP1-overexpressing ginseng hairy roots is mediated by tyrosinase activity inhibition. Taken together, our results suggested that the ectopic expression of PAP1 is an effective strategy for the enhancement of anthocyanin production, which improves the biological activities of ginseng root cultures.

MYB21 interacts with MYC2 to control the expression of terpene synthase genes in flowers of Freesia hybrida and Arabidopsis thaliana

Previously, linalool was found to be the most abundant component among the cocktail of volatiles released from flowers of Freesia hybrida. Linalool formation is catalysed by monoterpene synthase TPS1. However, the regulatory network developmentally modulating the expression of the TPS1 gene in Freesia hybrida remains unexplored. In this study, three regulatory genes, FhMYB21L1, FhMYB21L2, and FhMYC2, were screened from 52 candidates. Two MYB transcription factor genes were synchronously expressed with FhTPS1 and could activate its expression significantly when overexpressed, and the binding of FhMYB21L2 to the MYBCORE sites in the FhTPS1 promoter was further confirmed, indicating a direct role in activation. FhMYC2 showed an inverse expression pattern compared with FhTPS1; its expression led to a decreased binding of FhMYB21 to the FhTPS1 promoter to reduce its activation capacity when co-expressed, suggesting a role for an MYB-bHLH complex in the regulation of the FhTPS1 gene. In Arabidopsis, both MYB21 and MYC2 regulators were shown to activate the expression of sesquiterpene synthase genes, and the regulatory roles of AtMYB21 and AtMYC2 in the expression of the linalool synthase gene were also confirmed, implying conserved functions of the MYB-bHLH complex in these two evolutionarily divergent plants. Moreover, the expression ratio between MYB21 and MYC2 orthologues might be a determinant factor in floral linalool emission.

Function analysis and stress-mediated cis-element identification in the promoter region of VqMYB15

The transcription factor MYB15 plays an important role in grape basal immunity, and its promoter can be used as a potential target in resistance breeding. However, the regulatory mechanisms of cis-elements in its promoter region under a variety of stresses remain unclear. In this study, we identified some putative cis-regulatory elements present upstream of MYB15 in Vitis quinquangularis Shanyang (pVqMYB15_SY) and subsequently characterized the function of these elements using reporter assays. Our results showed that TCA-elements 1 and 2, ABRE, MYC and 3-AF1 binding site 1 are key cis-regulatory elements in pVqMYB15_SY and play important roles in plant bio/abiotic stress resistance.

Sites of biosynthesis and storage of Taxol in Taxus media (Rehder) plants: Mechanism of accumulation

Taxol is a cytotoxic agent against various types of cancers. The cytotoxic activities of Taxol can be extended to its synthesizing plant. Here, Taxol is shown to have special synthesis, storage and transport mechanisms that avoid the toxic effects on its source plant. The sites of Taxol biosynthesis, transport and storage were revealed by quantification of plant Taxol, its intermediate baccatin III, the polyphenol side chain precursor , gene expression analysis of the major Taxol biosynthetic genes and in situ immuno-labeling. Although the biosynthesis of Taxol was limited by the expression of its biosynthetic genes and the presence of baccatin III, its presence did not correlate to baccatin III accumulation, nor to the expression of biosynthetic genes. However, Taxol presence positively correlated to polyphenol accumulation (late stage in Taxol assembly) and the resin-like hydrophobic bodies (HB, storage organelles). These results indicate that the presence of Taxol requires two complementary steps, biosynthesis followed by storage. Each step is limited by the availability of different precursors, which differ in their localization within the plant. Thus, the sites of biosynthesis, transport and storage of Taxol are different. Taxus media (Rehder) plant wood showed high concentrations of baccatin III and the expression of biosynthetic genes. However, the concentrations of Taxol, polyphenol and HB were very high in the plant outer layers including phloem and dead bark (rhytidome). Furthermore, in situ immuno-labeling showed that taxadiene synthase (the rate-limiting enzyme in Taxol biosynthesis) was mainly found in the wood, while Taxol primarily localized to the outer tissues. Conclusively, wood can be considered as the site of Taxol biosynthesis. Our data also propose that Taxol then accumulates into HB in order to permit its transport within the living plant tissues without causing toxic effects. This is followed by Taxol storage in the outer tissues including phloem and dead bark.

The MYB transcription factor CiMYB42 regulates limonoids biosynthesis in citrus

BACKGROUND

Limonoids are major bioactive compounds that are produced by the triterpenoid metabolic pathway. The detailed biochemical process of limonoid biosynthesis and the mechanism of its molecular regulation remain elusive. The identification of transcription factors that regulate limonoid biosynthetic pathways is very important for understanding the underlying regulatory mechanisms. This information could also provide tools for manipulating biosynthesis genes to modulate limonoid production.

RESULTS

In this study, the CiMYB42 transcription factor was isolated to identify its role in limonoid biosynthesis. Multiple alignment analysis and phylogenetic analysis demonstrated that CiMYB42 is a typical R2R3MYB transcription factor that shares high similarity of its amino acid sequence with AtMYB42. Limonoids contents were higher in Citrus sinensis and Citrus grandis than in other species. Limonoid accumulation during leaf development also showed diverse trends in different genotypes. The expression of CiMYB42 was significantly related to the limonoid content and the expression of CiOSC in some citrus accessions. The overexpression of CiMYB42 in sweet orange resulted in significant accumulation of limonin, whereas the downregulation of CiMYB42 by RNAi resulted in a dwarf phenotype and less nomilin accumulation. Furthermore, the results of a yeast one-hybrid assay and EMSA indicated that CiMYB42 binds exclusively to the TTGTTG sequence (type II MYB core) in the promoter of CiOSC. Together, these results suggest that CiMYB42 positively regulates limonoid biosynthesis by regulating the expression of CiOSC by binding to the TTGTTG sequence (type II MYB core) of its promoter.

CONCLUSIONS

CiMYB42 is an important transcription activator involved in limonoid biosynthesis that regulates the expression of CiOSC by binding to the TTGTTG sequence (type II MYB core).

Parallel Transcriptional Regulation of Artemisinin and Flavonoid Biosynthesis

Plants regulate the synthesis of specialized compounds through the actions of individual transcription factors (TFs) or sets of TFs. One such compound, artemisinin from Artemisia annua, is widely used as a pharmacological product in the first-line treatment of malaria. However, the emergence of resistance to artemisinin in Plasmodium species, as well as its low production rates, have required innovative treatments such as exploiting the synergistic effects of flavonoids with artemisinin. We overview current knowledge about flavonoid and artemisinin transcriptional regulation in A. annua, and review the dual action of TFs and structural genes that can regulate both pathways simultaneously. Understanding the concerted action of these TFs and their associated structural genes can guide the development of strategies to further improve flavonoid and artemisinin production.

VviWRKY40, a WRKY Transcription Factor, Regulates Glycosylated Monoterpenoid Production by VviGT14 in Grape Berry

Glycosylated volatile precursors are important, particularly in wine grape berries, as they contribute to the final aroma in wines by releasing volatile aglycones during yeast fermentation and wine storage. Previous study demonstrated that VviGT14 was functioned as a critical monoterpene glucosyltransferase in grape berry, while the transcriptional regulation mechanism of VviGT14 was still unknown. Here we identified VviWRKY40 as a binding factor of VviGT14 promoter by both DNA pull-down and yeast one-hybrid screening, followed by a series of in vitro verification. VviWRKY40 expression pattern negatively correlated with that of VviGT14 in grape berries. And the suppressor role of VviWRKY40 was further confirmed by using the dual luciferase assay with Arabidopsis protoplast and grape cell suspension system. Furthermore, the grape suspension cell ABA treatment study showed that ABA downregulated VviWRKY40 transcript level but promoted that of VviGT14, indicating that VviWRKY40 was at the downstream of ABA signal transduction network to regulate monoterpenoid glycosylation. These data extend our knowledge of transcriptional regulation of VviGT14, and provide new targets for grape breeding to alter monoterpenoid composition.

Phenylpropanoid Pathway Engineering: An Emerging Approach towards Plant Defense

Pathogens hitting the plant cell wall is the first impetus that triggers the phenylpropanoid pathway for plant defense. The phenylpropanoid pathway bifurcates into the production of an enormous array of compounds based on the few intermediates of the shikimate pathway in response to cell wall breaches by pathogens. The whole metabolomic pathway is a complex network regulated by multiple gene families and it exhibits refined regulatory mechanisms at the transcriptional, post-transcriptional, and post-translational levels. The pathway genes are involved in the production of anti-microbial compounds as well as signaling molecules. The engineering in the metabolic pathway has led to a new plant defense system of which various mechanisms have been proposed including salicylic acid and antimicrobial mediated compounds. In recent years, some key players like phenylalanine ammonia lyases (PALs) from the phenylpropanoid pathway are proposed to have broad spectrum disease resistance (BSR) without yield penalties. Now we have more evidence than ever, yet little understanding about the pathway-based genes that orchestrate rapid, coordinated induction of phenylpropanoid defenses in response to microbial attack. It is not astonishing that mutants of pathway regulator genes can show conflicting results. Therefore, precise engineering of the pathway is an interesting strategy to aim at profitably tailored plants. Here, this review portrays the current progress and challenges for phenylpropanoid pathway-based resistance from the current prospective to provide a deeper understanding.

Genome-Wide Identification, Classification, and Expression Profiling Reveals R2R3-MYB Transcription Factors Related to Monoterpenoid Biosynthesis in Osmanthus fragrans

Osmanthus fragrans is widely grown for the purpose of urban greening and the pleasant aroma emitted from its flowers. The floral scent is determined by several monoterpenoid volatiles, such as linalool and its oxides, which are a few of the most common volatiles and the main components of the essential oils in most sweet osmanthus cultivars. In addition, the relative contents of cis- and trans-linalool oxide (furan) may affect the aromas and quality of the essential oils. MYB proteins represent the largest family of transcription factors in plants and participate in regulating secondary metabolites. Several cis-elements, especially AC-rich regions, are known to be bound by 2R-MYBs and could be found in the promoter of the enzyme genes in the terpenoid metabolic pathway. However, there has to date been no investigation into the 2R-MYB family genes involved in regulating terpenoid biosynthesis in O. fragrans. Here, 243 non-redundant 2R-MYB proteins were grouped into 33 clusters based on the phylogeny and exon-intron distribution. These genes were unevenly distributed on 23 chromosomes. Ka/Ks analysis showed that the major mode of 2R-MYB gene evolution was purifying selection. Expression analysis indicated that 2R-MYB genes in O. fragrans exhibited varied expression patterns. A total of 35 OfMYBs representing the highest per kilobase per million mapped reads in the flower were selected for quantitative real-time PCR analysis. The correlation analysis between the expression level and the contents of fragrant compounds at different flowering stages suggested that OfMYB19/20 exhibited remarkably positive correlation with the accumulation of cis-linalool oxides. OfMYB51/65/88/121/137/144 showed significantly negative correlations with one or more linalool oxides. Characterization of these proteins revealed that OfMYB19 and OfMYB137 were localized in the nuclei, but did not show transcriptional activation in the yeast system, which suggested that they may be bound to other transcription factors to exert regulatory functions. These findings provide useful information for further functional investigation of the 2R-MYBs and offer a foundation for clarifying the 2R-MYB transcription factors involved in the molecular mechanism of the regulation of monoterpenoid biosynthesis in Osmanthus fragrans.

Structural insights into target DNA recognition by R2R3-MYB transcription factors

As the largest group of MYB family transcription factors, R2R3-MYB proteins play essential roles during plant growth and development. However, the structural basis underlying how R2R3-MYBs recognize the target DNA remains elusive. Here, we report the crystal structure of Arabidopsis WEREWOLF (WER), an R2R3-MYB protein, in complex with its target DNA. Structural analysis showed that the third α-helices in both the R2 and R3 repeats of WER fit in the major groove of the DNA, specifically recognizing the DNA motif 5'-AACNGC-3'. In combination with mutagenesis, in vitro binding and in vivo luciferase assays, we showed that K55, N106, K109 and N110 are critical for the function of WER. Although L59 of WER is not involved in DNA binding in the structure, ITC analysis suggested that L59 plays an important role in sensing DNA methylation at the fifth position of cytosine (5mC). Like 5mC, methylation at the sixth position of adenine (6mA) in the AAC element also inhibits the interaction between WER and its target DNA. Our study not only unravels the molecular basis of how WER recognizes its target DNA, but also suggests that 5mC and 6mA modifications may block the interaction between R2R3-MYB transcription factors and their target genes.

Functional and morphological evolution in gymnosperms: A portrait of implicated gene families

Gymnosperms diverged from their sister plant clade of flowering plants 300 Mya. Morphological and functional divergence between the two major seed plant clades involved significant changes in their reproductive biology, water-conducting systems, secondary metabolism, stress defense mechanisms, and small RNA-mediated epigenetic silencing. The relatively recent sequencing of several gymnosperm genomes and the development of new genomic resources have enabled whole-genome comparisons within gymnosperms, and between angiosperms and gymnosperms. In this paper, we aim to understand how genes and gene families have contributed to the major functional and morphological differences in gymnosperms, and how this information can be used for applied breeding and biotechnology. In addition, we have analyzed the angiosperm versus gymnosperm evolution of the pleiotropic drug resistance (PDR) gene family with a wide range of functionalities in plants' interaction with their environment including defense mechanisms. Some of the genes reviewed here are newly studied members of gene families that hold potential for biotechnological applications related to commercial and pharmacological value. Some members of conifer gene families can also be exploited for their potential in phytoremediation applications.

Transcriptome and Metabolite Profiling Reveal Novel Insights into Volatile Heterosis in the Tea Plant (Camellia Sinensis)

Tea aroma is a key indicator for evaluating tea quality. Although notable success in tea aroma improvement has been achieved with heterosis breeding technology, the molecular basis underlying heterosis remains largely unexplored. Thus, the present report studies the tea plant volatile heterosis using a high-throughput next-generation RNA-seq strategy and gas chromatography-mass spectrometry. Phenotypically, we found higher terpenoid volatile and green leaf volatile contents by gas chromatography-mass spectrometry in the F1 hybrids than in their parental lines. Volatile heterosis was obvious in both F1 hybrids. At the molecular level, the comparative transcriptomics analysis revealed that approximately 41% (9027 of 21,995) of the genes showed non-additive expression, whereas only 7.83% (1723 of 21,995) showed additive expression. Among the non-additive genes, 42.1% showed high parental dominance and 17.6% showed over-dominance. Among different expression genes with high parental dominance and over-dominance expression patterns, KEGG and GO analyses found that plant hormone signal transduction, tea plant physiological process related pathways and most pathways associated with tea tree volatiles were enriched. In addition, we identified multiple genes (CsDXS, CsAATC2, CsSPLA2, etc.) and transcription factors (CsMYB1, CsbHLH79, CsWRKY40, etc.) that played important roles in tea volatile heterosis. Based on transcriptome and metabolite profiling, we conclude that non-additive action plays a major role in tea volatile heterosis. Genes and transcription factors involved in tea volatiles showing over-dominance expression patterns can be considered candidate genes and provide novel clues for breeding high-volatile tea varieties.

Functional Analysis of the PgCesA3 White Spruce Cellulose Synthase Gene Promoter in Secondary Xylem

Cellulose is an essential structural component of the plant cell wall. Its biosynthesis involves genes encoding cellulose synthase enzymes and a complex transcriptional regulatory network. Three cellulose synthases have been identified in conifers as being potentially involved in secondary cell wall biosynthesis because of their preferential expression in xylem tissues; however, no direct functional association has been made to date. In the present work, we characterized the white spruce [Picea glauca (Moench) Voss] cellulose synthase PgCesA3 gene and 5' regulatory elements. Phylogenetic analysis showed that PgCesA1-3 genes grouped with secondary cell wall-associated Arabidopsis cellulose synthase genes, such as AtCesA8, AtCesA4, and AtCesA7. We produced transgenic spruce expressing the GUS reporter gene driven by the PgCesA3 promoter. We observed blue staining in differentiating xylem cells from stem and roots, and in foliar guard cells indicating that PgCesA3 is clearly involved in secondary cell wall biosynthesis. The promoter region sequence of PgCesA3 contained several putative MYB cis-regulatory elements including AC-I like motifs and secondary wall MYB-responsive element (SMRE); however, it lacked SMRE4, 7 and 8 that correspond to the sequences of AC-I, II, and III. Based on these findings and results of previous transient trans-activation assays that identified interactions between the PgCesA3 promoter and different MYB transcription factors, we performed electrophoretic mobility shift assays with MYB recombinant proteins and cis-regulatory elements present in the PgCesA3 promoter. We found that PgMYB12 bound to a canonical AC-I element identified in the Pinus taeda PAL promoter and two AC-I like elements. We hypothesized that the PgMYB12 could regulate PgCesA3 in roots based on previous expression results. This functional study of PgCesA3 sequences and promoter opens the door for future studies on the interaction between PgMYBs and the PgCesA3 regulatory elements.