RNA-sequencing Reveals Global Transcriptomic Changes in Nicotiana tabacum Responding to Topping and Treatment of Axillary-shoot Control Chemicals

Integrative proteome and metabolome unveil the central role of IAA alteration in axillary bud development following topping in tobacco

Axillary bud is an important aspect of plant morphology, contributing to the final tobacco yield. However, the mechanisms of axillary bud development in tobacco remain largely unknown. To investigate this aspect of tobacco biology, the metabolome and proteome of the axillary buds before and after topping were compared. A total of 569 metabolites were differentially abundant before and 1, 3, and 5 days after topping. KEGG analyses further revealed that the axillary bud was characterized by a striking enrichment of metabolites involved in flavonoid metabolism, suggesting a strong flavonoid biosynthesis activity in the tobacco axillary bud after topping. Additionally, 9035 differentially expressed proteins (DEPs) were identified before and 1, 3, and 5 days after topping. Subsequent GO and KEGG analyses revealed that the DEPs in the axillary bud were enriched in oxidative stress, hormone signal transduction, MAPK signaling pathway, and starch and sucrose metabolism. The integrated proteome and metabolome analysis revealed that the indole-3-acetic acid (IAA) alteration in buds control dormancy release and sustained growth of axillary bud by regulating proteins involved in carbohydrate metabolism, amino acid metabolism, and lipid metabolism. Notably, the proteins related to reactive oxygen species (ROS) scavenging and flavonoid biosynthesis were strongly negatively correlated with IAA content. These findings shed light on a critical role of IAA alteration in regulating axillary bud outgrowth, and implied a potential crosstalk among IAA alteration, ROS homeostasis, and flavonoid biosynthesis in tobacco axillary bud under topping stress, which could improve our understanding of the IAA alteration in axillary bud as an important regulator of axillary bud development.

Transcriptomics of Leaf Development in the Endangered Dioecious Magnolia kwangsiensis: Molecular Basis Underpinning Specialized Metabolism Genes

Magnolia kwangsiensis, a dioecious tree native to China, is recognized not only for its status as an at-risk species but also for its potential in therapeutic applications courtesy of its bioactive compounds. However, the genetic underpinnings of its leaf development and compound biosynthesis are not well documented. Our study aims to bridge this knowledge gap through comparative transcriptomics, analyzing gene expression through different leaf maturation stages. We studied the transcriptome of M. kwangsiensis leaves by applying RNA sequencing at juvenile, tender, and mature phases. We identified differentially expressed genes (DEGs) to explore transcriptional changes accompanying the developmental trajectory. Our analysis delineates the transcriptional landscape of over 20,000 genes with over 6000 DEGs highlighting significant transcriptional shifts throughout leaf maturation. Mature leaves demonstrated upregulation in pathways related to photosynthesis, cell wall formation, and polysaccharide production, affirming their structural integrity and specialized metabolic functions. Our GO and KEGG enrichment analyses underpin these findings. Furthermore, we unveiled coordinated gene activity correlating development with synthesizing therapeutically relevant polysaccharides. We identified four novel glycosyltransferases potentially pivotal in this synergistic mechanism. Our study uncovers the complementary evolutionary forces that concurrently sculpt structural and chemical defenses. These genetic mechanisms calibrate leaf tissue resilience and biochemical efficacy.

Effects of topping and non-topping on growth-regulating hormones of flue-cured tobacco (Nicotiana tabacum L.)-a proteomic analysis

Introduction

Until now, the mechanism underlying the impact of topping on hormone regulation in tobacco plants remains unclear, and most studies investigating the hormone signaling pathways in plants rely on genes or transcriptional pathways.

Methods

This study examines the regulatory mechanisms of hormones in the roots and leaves of tobacco plants with and without topping at the protein level.

Results

The results demonstrate that, compared with non-topped plants, topping leads to a decrease in the levels of IAA (auxin), ABA (abscisic acid), and GA (gibberellin) hormones in the leaves, whereas the content of the JA (jasmonic acid) hormone increases. Furthermore, in the roots, topping results in an increase in the levels of IAA, ABA, and JA hormones, along with a decrease in GA content. In the leaves, a total of 258 significantly different proteins were identified before and after topping, with 128 proteins upregulated and 130 proteins downregulated. In the roots, there were 439 proteins with significantly different quantities before and after topping, consisting of 211 upregulated proteins and 228 downregulated proteins. Notably, these proteins were closely associated with the metabolic and biosynthetic pathways of secondary metabolites, as indicated by functional categorization.

Conclusions

When integrating the hormone changes and the proteomics results, it is evident that topping leads to increased metabolic activity and enhanced hormone synthesis in the root system. This research provides a theoretical foundation for further investigations into the regulation and signaling mechanisms of hormones at the protein level before and after topping in plants.

Upcoming progress of transcriptomics studies on plants: An overview

Transcriptome sequencing or RNA-Sequencing is a high-resolution, sensitive and high-throughput next-generation sequencing (NGS) approach used to study non-model plants and other organisms. In other words, it is an assembly of RNA transcripts from individual or whole samples of functional and developmental stages. RNA-Seq is a significant technique for identifying gene predictions and mining functional analysis that improves gene ontology understanding mechanisms of biological processes, molecular functions, and cellular components, but there is limited information available on this topic. Transcriptomics research on different types of plants can assist researchers to understand functional genes in better ways and regulatory processes to improve breeding selection and cultivation practices. In recent years, several advancements in RNA-Seq technology have been made for the characterization of the transcriptomes of distinct cell types in biological tissues in an efficient manner. RNA-Seq technologies are briefly introduced and examined in terms of their scientific applications. In a nutshell, it introduces all transcriptome sequencing and analysis techniques, as well as their applications in plant biology research. This review will focus on numerous existing and forthcoming strategies for improving transcriptome sequencing technologies for functional gene mining in various plants using RNA- Seq technology, based on the principles, development, and applications.

Metabolic engineering to enhance the accumulation of bioactive flavonoids licochalcone A and echinatin in Glycyrrhiza inflata (Licorice) hairy roots

Echinatin and licochalcone A (LCA) are valuable chalcones preferentially accumulated in roots and rhizomes of licorice (Glycyrrhiza inflata). The licorice chalcones (licochalcones) are valued for their anti-inflammatory, antimicrobial, and antioxidant properties and have been widely used in cosmetic, pharmaceutical, and food industries. However, echinatin and LCA are accumulated in low quantities, and the biosynthesis and regulation of licochalcones have not been fully elucidated. In this study, we explored the potential of a R2R3-MYB transcription factor (TF) AtMYB12, a known regulator of flavonoid biosynthesis in Arabidopsis, for metabolic engineering of the bioactive flavonoids in G. inflata hairy roots. Overexpression of AtMYB12 in the hairy roots greatly enhanced the production of total flavonoids (threefold), echinatin (twofold), and LCA (fivefold). RNA-seq analysis of AtMYB12-overexpressing hairy roots revealed that expression of phenylpropanoid/flavonoid pathway genes, such as phenylalanine ammonia-lyase (PAL), chalcone synthase (CHS), and flavanone 3'-hydroxylase (F3'H), is significantly induced compared to the control. Transient promoter activity assay indicated that AtMYB12 activates the GiCHS1 promoter in plant cells, and mutation to the MYB-binding motif in the GiCHS1 promoter abolished activation. In addition, transcriptomic analysis revealed that AtMYB12 overexpression reprograms carbohydrate metabolism likely to increase carbon flux into flavonoid biosynthesis. Further, AtMYB12 activated the biotic defense pathways possibly by activating the salicylic acid and jasmonic acid signaling, as well as by upregulating WRKY TFs. The transcriptome of AtMYB12-overexpressing hairy roots serves as a valuable source in the identification of potential candidate genes involved in LCA biosynthesis. Taken together, our findings suggest that AtMYB12 is an effective gene for metabolic engineering of valuable bioactive flavonoids in plants.

The chrysanthemum lavandulifolium genome and the molecular mechanism underlying diverse capitulum types

Cultivated chrysanthemum (Chrysanthemum × morifolium Ramat.) is a beloved ornamental crop due to the diverse capitula types among varieties, but the molecular mechanism of capitulum development remains unclear. Here, we report a 2.60 Gb chromosome-scale reference genome of C. lavandulifolium, a wild Chrysanthemum species found in China, Korea and Japan. The evolutionary analysis of the genome revealed that only recent tandem duplications occurred in the C. lavandulifolium genome after the shared whole genome triplication (WGT) in Asteraceae. Based on the transcriptomic profiling of six important developmental stages of the radiate capitulum in C. lavandulifolium, we found genes in the MADS-box, TCP, NAC and LOB gene families that were involved in disc and ray floret primordia differentiation. Notably, NAM and LOB30 homologs were specifically expressed in the radiate capitulum, suggesting their pivotal roles in the genetic network of disc and ray floret primordia differentiation in chrysanthemum. The present study not only provides a high-quality reference genome of chrysanthemum but also provides insight into the molecular mechanism underlying the diverse capitulum types in chrysanthemum.

A spatial-temporal understanding of gene regulatory networks and NtARF-mediated regulation of potassium accumulation in tobacco

MAIN CONCLUSION

After tobacco topping, changes in the auxin content could affect K+ uptake by inhibiting the activity of K+ uptake-related genes through the NtARF genes, thus causing changes in K+ content. Tobacco (Nicotiana tabacum) is a valuable industrial and commercial crop, and the leaf is its primary product. Topping (removing apical buds) is a common agronomic practice that significantly improves the yield of tobacco leaves. Potassium (K+) plays an important physiological role in tobacco growth and leaf traits, including combustibility, aroma, and safety in cigarette products, and its levels are significantly decreased after topping. Here, to present global spatial-temporal gene expression profiles and gene regulatory networks of the core elements of K+ uptake, leaves and roots from topped and untopped plants at short- and long-term time points after topping were sampled for transcriptome analysis. We found that the wounding response was initiated in leaves in the early stages after topping. Then, in the long term, processes related to metabolism and transcription regulation, as well as ion binding and transport, were altered. The expression profiles showed that core elements of K+ uptake and xylem loading were drastically suppressed in roots after topping. Finally, transient expression experiments confirmed that changes in the auxin content could affect K+ uptake by inhibiting the activity of K+ uptake-related genes through the tobacco auxin response factor (NtARF) genes, thus causing changes in the K+ content. These results suggest that some ARFs could be selected as targets to enhance the expressions of K+ uptake transporters, leading to increment of K+ contents and improvement of leaf quality in tobacco breeding.

Protein phosphatase NtPP2C2b and MAP kinase NtMPK4 act in concert to modulate nicotine biosynthesis

Protein phosphatases (PPs) and protein kinases (PKs) regulate numerous developmental, defense, and phytohormone signaling processes in plants. However, the underlying regulatory mechanism governing biosynthesis of specialized metabolites, such as alkaloids, by the combined effects of PPs and PKs, is insufficiently understood. Here, we report the characterization of a group B protein phosphatase type 2C, NtPP2C2b, that likely acts upstream of the NICOTINE2 locus APETALA 2/Ethylene Response Factors (AP2/ERFs), to regulate nicotine biosynthesis in tobacco. Similar to the nicotine pathway genes, NtPP2C2b is highly expressed in roots and induced by jasmonic acid (JA). Overexpression of NtPP2C2b in transgenic hairy roots or stable transgenic tobacco plants repressed nicotine pathway gene expression and reduced nicotine accumulation. Additionally, transient overexpression of NtPP2C2b, together with the NtERF221, repressed transactivation of the quinolinate phosphoribosyltransferase promoter in tobacco cells. We further demonstrate that the JA-responsive tobacco mitogen-activated protein kinase (MAPK) 4 interacts with NtPP2C2b in yeast and plant cells. Conditional overexpression of NtMPK4 in tobacco hairy roots up-regulated nicotine pathway gene expression and increased nicotine accumulation. Our findings suggest that a previously uncharacterized PP-PK module acts to modulate alkaloid biosynthesis, highlighting the importance of post-translational control in the biosynthesis of specialized plant metabolites.

BHLH IRIDOID SYNTHESIS 3 is a member of a bHLH gene cluster regulating terpenoid indole alkaloid biosynthesis in Catharanthus roseus

Basic helix-loop-helix (bHLH) transcription factors (TFs) are key regulators of plant specialized metabolites, including terpenoid indole alkaloids (TIAs) in Catharanthus roseus. Two previously characterized subgroup-IVa bHLH TFs, BIS1 (bHLH Iridoid Synthesis 1) and BIS2 regulate iridoid biosynthesis in the TIA pathway. We reanalyzed the recently updated C. roseus genome sequence and discovered that BIS1 and BIS2 are clustered on the same genomic scaffold with a previously uncharacterized bHLH gene, designated as BIS3. Only a few bHLH gene clusters have been studied to date. Comparative analysis of 49 genome sequences from different plant lineages revealed the presence of analogous bHLH clusters in core angiosperms, including the medicinal plants Calotropis gigantea (giant milkweed) and Gelsemium sempervirens (yellow jessamine), but not in the analyzed basal angiosperm and lower plants. Similar to the iridoid pathway genes, BIS3 is highly expressed in roots and induced by methyl jasmonate. BIS3 activates the promoters of iridoid branch genes, geraniol synthase (GES), geraniol 10-hydroxylase (G10H), 8-hydroxygeraniol oxidoreductase (8HGO), iridoid synthase (IS), 7-deoxyloganetic acid glucosyl transferase (7-DLGT), and 7-deoxyloganic acid hydroxylase (7DLH), but not iridoid oxidase (IO). Transactivation of the promoters was abolished when BIS3 is converted to a dominant repressor by fusing with the ERF-associated amphiphilic repression (EAR) sequence. In addition, BIS3 acts synergistically with BIS1 and BIS2 to activate the G10H promoter in tobacco cells. Mutation of the known bHLH TF binding motif, G-box (CACGTG) in the G10H promoter significantly reduced but did not abolish the transactivation by BIS3. Promoter deletion analysis of G10H suggests that the sequences adjacent to the G-box are also involved in the regulation by BIS3. Overexpression of BIS3 in C. roseus flower petals significantly upregulated the expression of iridoid biosynthetic genes and increased loganic acid accumulation. BIS2 expression was significantly induced by BIS3 although BIS3 did not directly activate the BIS2 promoter. Our results advance our understanding of the regulation of plant specialized metabolites by bHLH TF clusters.

Maleic hydrazide elicits global transcriptomic changes in chemically topped tobacco to influence shoot bud development

MAIN CONCLUSION

Transcriptomic analysis revealed maleic hydrazide suppresses apical and axillary bud development by altering the expression of genes related to meristem development, cell division, DNA replication, DNA damage and recombination, and phytohormone signaling. Topping (removal of apical buds) is a common agricultural practice for some crop plants including cotton, cannabis, and tobacco. Maleic hydrazide (MH) is a systemic suckercide, a chemical that inhibits shoot bud growth, used to control the growth of apical (ApB) and axillary buds (AxB) following topping. However, the influence of MH on gene expression and the underlying molecular mechanism of controlling meristem development are not well studied. Our RNA sequencing analysis showed that MH significantly influences the transcriptomic landscape in ApB and AxB of chemically topped tobacco. Gene ontology (GO) enrichment analysis revealed that upregulated genes in ApB were enriched for phosphorelay signal transduction, and the regulation of transition timing from vegetative to reproductive phase, whereas downregulated genes were largely associated with meristem maintenance, cytokinin metabolism, cell wall synthesis, photosynthesis, and DNA metabolism. In MH-treated AxB, GO terms related to defense response and oxylipin metabolism were overrepresented in upregulated genes. GO terms associated with cell cycle, DNA metabolism, and cytokinin metabolism were enriched in downregulated genes. Expression of KNOX and MADS transcription factor (TF) family genes, known to be involved in meristem development, were affected in ApB and AxB by MH treatment. The promoters of MH-responsive genes are enriched for several known cis-acting elements, suggesting the involvement of a subset of TF families. Our findings suggest that MH affects shoot bud development in chemically topped tobacco by altering the expression of genes related to meristem development, DNA repair and recombination, cell division, and phytohormone signaling.

Plants with genetically encoded autoluminescence

Autoluminescent plants engineered to express a bacterial bioluminescence gene cluster in plastids have not been widely adopted because of low light output. We engineered tobacco plants with a fungal bioluminescence system that converts caffeic acid (present in all plants) into luciferin and report self-sustained luminescence that is visible to the naked eye. Our findings could underpin development of a suite of imaging tools for plants.

Genetic attenuation of alkaloids and nicotine content in tobacco (Nicotiana tabacum)

MAIN CONCLUSION

The role of six alkaloid biosynthesis genes in the process of nicotine accumulation in tobacco was investigated. Downregulation of ornithine decarboxylase, arginine decarboxylase, and aspartate oxidase resulted in viable plants with a significantly lower nicotine content. Attenuation of nicotine accumulation in Nicotiana tabacum was addressed upon the application of RNAi technologies. The approach entailed a downregulation in the expression of six different alkaloid biosynthesis genes encoding upstream enzymes that are thought to function in the pathway of alkaloid and nicotine biosynthesis. Nine different RNAi constructs were designed to lower the expression level of the genes that encode the enzymes arginine decarboxylase, agmatine deiminase, aspartate oxidase, arginase, ornithine decarboxylase, and SAM synthase. Agrobacterium-based transformation of tobacco leaves was applied, and upon kanamycin selection, T0 and subsequently T1 generation seeds were produced. Mature T1 plants in the greenhouse were topped to prevent flowering and leaf nos. 3 and 4 below the topping point were tested for transcript levels and product accumulation. Down-regulation in arginine decarboxylase, aspartate oxidase, and ornithine decarboxylase consistently resulted in lower levels of nicotine in the leaves of the corresponding plants. Transformants with the aspartate oxidase RNAi construct showed the lowest nicotine level in the leaves, which varied from below the limit of quantification (20 μg per g dry leaf weight) to 1.3 mg per g dry leaf weight. The amount of putrescine, the main polyamine related to nicotine biosynthesis, showed a qualitative correlation with the nicotine content in the arginine decarboxylase and ornithine decarboxylase RNAi-expressing transformants. A putative early senescence phenotype and lower viability of the older leaves was observed in some of the transformant lines. The results are discussed in terms of the role of the above-mentioned genes in the alkaloid biosynthetic pathway and may serve to guide efforts to attenuate nicotine content in tobacco leaves.

Genome-wide analysis of glycerol-3-phosphate O-acyltransferase gene family and functional characterization of two cutin group GPATs in Brassica napus

Genome-wide identification, spatio-temporal expression analysis and functional characterization of selected Brassica napus GPATs highlight their roles in cuticular wax biosynthesis and defense against fungal pathogens. Glycerol-3-phosphate 1-O-acyltransferase (GPAT) is a key enzyme in the biosynthesis of glycerolipids, a major component of cellular membranes and extracellular protective layers, such as cuticles in plants. Brassica napus is an economically important crop and cultivated worldwide mostly for its edible oil. The B. napus GPATs (BnGPATs) are insufficiently characterized. Here, we performed genome-wide analysis to identify putative GPATs in B. napus and its diploid progenitors B. rapa and B oleracea. The 32 B. napus BnGPATs are phylogenetically divided into three major groups, cutin, suberin, and diverse ancient groups. Analysis of transcriptomes of different tissues and seeds at different developmental stages revealed the spatial and temporal expression profiles of BnGPATs. The yield and oil quality of B. napus are adversely affected by the necrotrophic fungus, Sclerotinia sclerotiorum. We showed that several BnGPATs, including cutin-related BnGPAT19 and 21, were upregulated in the S. sclerotiorum resistant line. RNAi-mediated suppression of BnGPAT19 and 21 in B. napus resulted in thinner cuticle, leading to rapid water and chlorophyll loss in toluidine blue staining and leaf bleaching assays. In addition, the RNAi plants also developed severe necrotic lesions following fungal inoculation compared to the wild-type plants, indicating that BnGPAT19 and 21 are likely involved in cuticular wax biosynthesis that is critical for initial pathogen defense. Taken together, we provided a comprehensive account of GPATs B. napus and characterized BnGPAT19 and 21 for their potential roles in cuticular wax biosynthesis and defense against fungal pathogens.

Transcriptome sequencing reveals the effect of biochar improvement on the development of tobacco plants before and after topping

The application of biochar is one of the most useful methods for improving soil quality, which is of the utmost significance for the continuous production of crops. As there are no conclusive studies on the specific effects of biochar application on tobacco quality, this study aimed to improve the yield and quality of tobacco as a model crop for economic and genetic research in southern China, by such application. We used transcriptome sequencing to reveal the effects of applied biochar on tobacco development before and after topping. Our results showed that topping affected carbon and nitrogen metabolism, photosynthesis and secondary metabolism in the tobacco plants, while straw biochar-application to the soil resulted in amino acid and lipid synthesis; additionally, it affected secondary metabolism of the tobacco plants through carbon restoration and hormonal action, before and after topping. In addition to the new insights into the impact of biochar on crops, our findings provide a basis for biochar application measures in tobacco and other crops.

Cross-family transcription factor interaction between MYC2 and GBFs modulates terpenoid indole alkaloid biosynthesis

Biosynthesis of medicinally valuable terpenoid indole alkaloids (TIAs) in Catharanthus roseus is regulated by transcriptional activators such as the basic helix-loop-helix factor CrMYC2. However, the transactivation effects are often buffered by repressors, such as the bZIP factors CrGBF1 and CrGBF2, possibly to fine-tune the accumulation of cytotoxic TIAs. Questions remain as to whether and how these factors interact to modulate TIA production. We demonstrated that overexpression of CrMYC2 induces CrGBF expression and results in reduced alkaloid accumulation in C. roseus hairy roots. We found that CrGBF1 and CrGBF2 form homo- and heterodimers to repress the transcriptional activities of key TIA pathway gene promoters. We showed that CrGBFs dimerize with CrMYC2, and CrGBF1 binds to the same cis-elements (T/G-box) as CrMYC2 in the target gene promoters. Our findings suggest that CrGBFs antagonize CrMYC2 transactivation possibly by competitive binding to the T/G-box in the target promoters and/or protein-protein interaction that forms a non-DNA binding complex that prevents CrMYC2 from binding to its target promoters. Homo- and heterodimer formation allows fine-tuning of the amplitude of TIA gene expression. Our findings reveal a previously undescribed regulatory mechanism that governs the TIA pathway genes to balance metabolic flux for TIA production in C. roseus.

Genome-wide identification of hexokinase gene family in Brassica napus: structure, phylogenetic analysis, expression, and functional characterization

Genome-wide identification, expression analysis, and functional characterization of previously uncharacterized hexokinase family of oil crop, Brassica napus, underscore the importance of this gene family in plant growth and development. In plants, the multi-gene family of dual-function hexokinases (HXKs) plays important roles in sugar metabolism and sensing that affect growth and development. Rapeseed (Brassica napus L.) is an important oil crop; however, little is known about the B. napus HXK gene family. We identified 19 putative HXKs in B. napus genome. B. rapa and B. oleracea, the two diploid progenitors of B. napus, contributed almost equally to the BnHXK genes. Phylogenetic analysis divided the 19 BnHXKs into four groups. The exon-intron structures of BnHXKs share high similarity to those of HXKs in Arabidopsis and rice. The group III and IV BnHXKs are highly expressed in roots, whereas group I members preferentially express in leaves. Analysis of seed transcriptomes at different developmental stages showed that most of group I and IV HXKs are highly expressed 2-weeks after pollination (2WAP), compared to 4WAP for group III. BnHKXs are differentially expressed in susceptible and tolerant B. napus cultivars after fungal infection, suggesting the possible involvement in defense response. We generated rapeseed RNAi lines for BnHXK9, a member of relatively less characterized group IV, by pollen-mediated gene transformation. The seedlings of BnHXK9-RNAi lines showed delayed growth compared to the wild type. The RNAi plants were dwarf with curly leaves, suggesting the involvement of BnHXK9 in plant development. Collectively, our findings provides a comprehensive account of BnHXK gene family in an important crop and a starting point for further elucidation of their roles in sugar metabolism and sensing, as well as plant growth and development.

Transcriptomic analysis of topping-induced axillary shoot outgrowth in Nicotiana tabacum

Topping is an important agronomic practice that significantly impacts the yield of various crop plants. Topping and the regulation of axillary shoot outgrowth are common agronomic practices in tobacco. However, the effects of topping on gene expression in tobacco remain unknown. We applied the Illumina HiSeq™ 2000 platform and analyzed differentially expressed genes (DEGs) from untopped and topped plants to study the global changes in gene expression in response to topping. We found that the number of DEGs varied from 7609 to 18,770 based on the reads per kilobase per million mapped reads (RPKM) values. The Gene Ontology (GO) enrichment analysis revealed that the cellular carbohydrate metabolic process and the disaccharide metabolic process, which may contribute to starch accumulation and stress/defense, were overrepresented terms for the DEGs. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that many DEGs were involved in starch and sucrose metabolism, glycolysis/gluconeogenesis, pyruvate metabolism, and plant hormone signal transduction, among other processes. The knowledge gained will improve our understanding of the processes of axillary shoot formation and enlargement at the transcriptional level. This study lays a solid foundation for future studies on molecular mechanisms underlying the growth of axillary shoots.

The miRNAome of Catharanthus roseus: identification, expression analysis, and potential roles of microRNAs in regulation of terpenoid indole alkaloid biosynthesis

MicroRNAs (miRNAs) regulate numerous crucial biological processes in plants. However, information is limited on their involvement in the biosynthesis of specialized metabolites in plants, including Catharanthus roseus that produces a number of pharmaceutically valuable, bioactive terpenoid indole alkaloids (TIAs). Using small RNA-sequencing, we identified 181 conserved and 173 novel miRNAs (cro-miRNAs) in C. roseus seedlings. Genome-wide expression analysis revealed that a set of cro-miRNAs are differentially regulated in response to methyl jasmonate (MeJA). In silico target prediction identified 519 potential cro-miRNA targets that include several auxin response factors (ARFs). The presence of cleaved transcripts of miRNA-targeted ARFs in C. roseus cells was confirmed by Poly(A) Polymerase-Mediated Rapid Amplification of cDNA Ends (PPM-RACE). We showed that auxin (indole acetic acid, IAA) repressed the expression of key TIA pathway genes in C. roseus seedlings. Moreover, we demonstrated that a miRNA-regulated ARF, CrARF16, binds to the promoters of key TIA pathway genes and repress their expression. The C. roseus miRNAome reported here provides a comprehensive account of the cro-miRNA populations, as well as their abundance and expression profiles in response to MeJA. In addition, our findings underscore the importance of miRNAs in posttranscriptional control of the biosynthesis of specialized metabolites.

A differentially regulated AP2/ERF transcription factor gene cluster acts downstream of a MAP kinase cascade to modulate terpenoid indole alkaloid biosynthesis in Catharanthus roseus

Catharanthus roseus produces bioactive terpenoid indole alkaloids (TIAs), including the chemotherapeutics, vincristine and vinblastine. Transcriptional regulation of TIA biosynthesis is not fully understood. The jasmonic acid (JA)-responsive AP2/ERF transcription factor (TF), ORCA3, and its regulator, CrMYC2, play key roles in TIA biosynthesis. ORCA3 forms a physical cluster with two uncharacterized AP2/ERFs, ORCA4 and 5. Here, we report that (1) the ORCA gene cluster is differentially regulated; (2) ORCA4, while overlapping functionally with ORCA3, modulates an additional set of TIA genes. Unlike ORCA3, ORCA4 overexpression resulted in dramatic increase of TIA accumulation in C. roseus hairy roots. In addition, CrMYC2 is capable of activating ORCA3 and co-regulating TIA pathway genes concomitantly with ORCA3. The ORCA gene cluster and CrMYC2 act downstream of a MAP kinase cascade that includes a previously uncharacterized MAP kinase kinase, CrMAPKK1. Overexpression of CrMAPKK1 in C. roseus hairy roots upregulated TIA pathways genes and increased TIA accumulation. This work provides detailed characterization of a TF gene cluster and advances our understanding of the transcriptional and post-translational regulatory mechanisms that govern TIA biosynthesis in C. roseus.

NtWRKY-R1, a Novel Transcription Factor, Integrates IAA and JA Signal Pathway under Topping Damage Stress in Nicotiana tabacum

Topping damage can induce the nicotine synthesis in tobacco roots, which involves the activation of JA and auxin signal transduction. It remains unclear how these hormone signals are integrated to regulate nicotine synthesis. Here we isolated a transcription factor NtWRKY-R1 from the group IIe of WRKY family and it had strong negative correlation with the expression of putrescine N-methyltransferase, the key enzyme of nicotine synthesis pathway. NtWRKY-R1 was specifically and highly expressed in tobacco roots, and it contains two transcriptional activity domains in the N- and C-terminal. The promoter region of NtWRKY-R1 contains two cis-elements which are responding to JA and auxin signals, respectively. Deletion of NtWRKY-R1 promoter showed that JA and auxin signals were subdued by NtWRKY-R1, and the expression of NtWRKY-R1 was more sensitive to auxin than JA. Furthermore, Yeast two-hybrid experiment demonstrated that NtWRKY-R1 can interact with the actin-binding protein. Our data showed that the intensity of JA and auxin signals can be translated into the expression of NtWRKY-R1, which regulates the balance of actin polymerization and depolymerization through binding actin-binding protein, and then regulates the expression of genes related to nicotine synthesis. The results will help us better understand the function of the WRKY-IIe family in the signaling crosstalk of JA and auxin under damage stress.

Comparative Transcriptomic Analysis of Two Brassica napus Near-Isogenic Lines Reveals a Network of Genes That Influences Seed Oil Accumulation

Rapeseed (Brassica napus) is an important oil seed crop, providing more than 13% of the world's supply of edible oils. An in-depth knowledge of the gene network involved in biosynthesis and accumulation of seed oil is critical for the improvement of B. napus. Using available genomic and transcriptomic resources, we identified 1,750 acyl-lipid metabolism (ALM) genes that are distributed over 19 chromosomes in the B. napus genome. B. rapa and B. oleracea, two diploid progenitors of B. napus, contributed almost equally to the ALM genes. Genome collinearity analysis demonstrated that the majority of the ALM genes have arisen due to genome duplication or segmental duplication events. In addition, we profiled the expression patterns of the ALM genes in four different developmental stages. Furthermore, we developed two B. napus near isogenic lines (NILs). The high oil NIL, YC13-559, accumulates significantly higher (∼10%) seed oil compared to the other, YC13-554. Comparative gene expression analysis revealed upregulation of lipid biosynthesis-related regulatory genes in YC13-559, including SHOOTMERISTEMLESS, LEAFY COTYLEDON 1 (LEC1), LEC2, FUSCA3, ABSCISIC ACID INSENSITIVE 3 (ABI3), ABI4, ABI5, and WRINKLED1, as well as structural genes, such as ACETYL-CoA CARBOXYLASE, ACYL-CoA DIACYLGLYCEROL ACYLTRANSFERASE, and LONG-CHAIN ACYL-CoA SYNTHETASES. We observed that several genes related to the phytohormones, gibberellins, jasmonate, and indole acetic acid, were differentially expressed in the NILs. Our findings provide a broad account of the numbers, distribution, and expression profiles of acyl-lipid metabolism genes, as well as gene networks that potentially control oil accumulation in B. napus seeds. The upregulation of key regulatory and structural genes related to lipid biosynthesis likely plays a major role for the increased seed oil in YC13-559.