The Genetics of Plant Metabolism.

Plant metabolic studies have traditionally focused on the role and regulation of the enzymes catalyzing key reactions within specific pathways. Within the past 20 years, reverse genetic approaches have allowed direct determination of the effects of the deficiency, or surplus, of a given protein on the biochemistry of a plant. In parallel, top-down approaches have also been taken, which rely on screening broad, natural genetic diversity for metabolic diversity. Here, we compare and contrast the various strategies that have been adopted to enhance our understanding of the natural diversity of metabolism. We also detail how these approaches have enhanced our understanding of both specific and global aspects of the genetic regulation of metabolism. Finally, we discuss how such approaches are providing important insights into the evolution of plant secondary metabolism.

Genome-wide association studies identify loci controlling specialized seed metabolites in Arabidopsis.

Plants synthesize specialized metabolites to facilitate environmental and ecological interactions. During evolution, plants diversified in their potential to synthesize these metabolites. Quantitative differences in metabolite levels of natural Arabidopsis (Arabidopsis thaliana) accessions can be employed to unravel the genetic basis for metabolic traits using genome-wide association studies (GWAS). Here, we performed metabolic GWAS on seeds of a panel of 315 A. thaliana natural accessions, including the reference genotypes C24 and Col-0, for polar and semi-polar seed metabolites using untargeted ultra-performance liquid chromatography-mass spectrometry. As a complementary approach, we performed quantitative trait locus (QTL) mapping of near-isogenic introgression lines between C24 and Col-0 for specific seed specialized metabolites. Besides common QTL between seeds and leaves, GWAS revealed seed-specific QTL for specialized metabolites, indicating differences in the genetic architecture of seeds and leaves. In seeds, aliphatic methylsulfinylalkyl and methylthioalkyl glucosinolates associated with the ALKENYL HYDROXYALKYL PRODUCING loci (GS-ALK and GS-OHP) on chromosome 4 containing alkenyl hydroxyalkyl producing 2 (AOP2) and 3 (AOP3) or with the GS-ELONG locus on chromosome 5 containing methylthioalkyl malate synthase (MAM1) and MAM3. We detected two unknown sulfur-containing compounds that were also mapped to these loci. In GWAS, some of the annotated flavonoids (kaempferol 3-O-rhamnoside-7-O-rhamnoside, quercetin 3-O-rhamnoside-7-O-rhamnoside) were mapped to transparent testa 7 (AT5G07990), encoding a cytochrome P450 75B1 monooxygenase. Three additional mass signals corresponding to quercetin-containing flavonols were mapped to UGT78D2 (AT5G17050). The association of the loci and associating metabolic features were functionally verified in knockdown mutant lines. By performing GWAS and QTL mapping, we were able to leverage variation of natural populations and parental lines to study seed specialized metabolism. The GWAS data set generated here is a high-quality resource that can be investigated in further studies.

Decreased metabolic diversity in common beans associated with domestication revealed by untargeted metabolomics, information theory, and molecular networking.

The process of crop domestication leads to a dramatic reduction in the gene expression associated with metabolic diversity. Genes involved in specialized metabolism appear to be particularly affected. Although there is ample evidence of these effects at the genetic level, a reduction in diversity at the metabolite level has been taken for granted despite having never been adequately accessed and quantified. Here we leveraged the high coverage of ultra high performance liquid chromatography-high-resolution mass spectrometry based metabolomics to investigate the metabolic diversity in the common bean (Phaseolus vulgaris). Information theory highlights a shift towards lower metabolic diversity and specialization when comparing wild and domesticated bean accessions. Moreover, molecular networking approaches facilitated a broader metabolite annotation than achieved to date, and its integration with gene expression data uncovers a metabolic shift from specialized metabolism towards central metabolism upon domestication of this crop.

Mass spectrometry-based metabolomics: a guide for annotation, quantification and best reporting practices.

Mass spectrometry-based metabolomics approaches can enable detection and quantification of many thousands of metabolite features simultaneously. However, compound identification and reliable quantification are greatly complicated owing to the chemical complexity and dynamic range of the metabolome. Simultaneous quantification of many metabolites within complex mixtures can additionally be complicated by ion suppression, fragmentation and the presence of isomers. Here we present guidelines covering sample preparation, replication and randomization, quantification, recovery and recombination, ion suppression and peak misidentification, as a means to enable high-quality reporting of liquid chromatography- and gas chromatography-mass spectrometry-based metabolomics-derived data.

Unveiling the regulatory role of DzAGL6-1 in carotenoid biosynthesis during durian (Durio zibethinus) fruit development.

KEY MESSAGE: Approximately 119 MADS-box genes have been identified in durian. Moreover, DzAGL6-1 primarily expressed during fruit development, activates the DzPSY promoter. Transient expression of DzAGL6-1 in tomatoes influences carotenoid production. MADS-box transcription factors play a crucial role in regulating plant biological processes, including fruit ripening and associated events. This study aimed to comprehend the mechanisms involved in durian fruit development and ripening and carotenoid production by conducting a genome-wide analysis of MADS-box proteins in durian (Durio zibethinus L.), an economically important fruit in Southeast Asia. A total of 119 durian MADS-box proteins were identified from the genome of the 'Musang King' cultivar. Based on the phylogenetic analysis, the proteins were classified into types I and II, which exhibited similar conserved motif compositions. Notably, only 16 durian MADS-box genes exhibited fruit-specific expression patterns. Among these genes, DzAGL6-1 was predominantly expressed during fruit development, a stage at which carotenoid biosynthesis is activated. Transient expression of DzAGL6-1 in tomato fruit increased the transcript level of the carotenoid biosynthetic gene phytoene synthase (PSY) and the ß-carotene content. Furthermore, DzAGL6-1 activated the promoter activity of DzPSY, as demonstrated by a dual-luciferase assay. These findings provide insights into the role of MADS-box transcription factors in regulating carotenoid biosynthesis during durian fruit development.

Coordinating the morphogenesis-differentiation balance by tweaking the cytokinin-gibberellin equilibrium.

Morphogenesis and differentiation are important stages in organ development and shape determination. However, how they are balanced and tuned during development is not fully understood. In the compound leaved tomato, an extended morphogenesis phase allows for the initiation of leaflets, resulting in the compound form. Maintaining a prolonged morphogenetic phase in early stages of compound-leaf development in tomato is dependent on delayed activity of several factors that promote differentiation, including the CIN-TCP transcription factor (TF) LA, the MYB TF CLAU and the plant hormone Gibberellin (GA), as well as on the morphogenesis-promoting activity of the plant hormone cytokinin (CK). Here, we investigated the genetic regulation of the morphogenesis-differentiation balance by studying the relationship between LA, CLAU, TKN2, CK and GA. Our genetic and molecular examination suggest that LA is expressed earlier and more broadly than CLAU and determines the developmental context of CLAU activity. Genetic interaction analysis indicates that LA and CLAU likely promote differentiation in parallel genetic pathways. These pathways converge downstream on tuning the balance between CK and GA. Comprehensive transcriptomic analyses support the genetic data and provide insights into the broader molecular basis of differentiation and morphogenesis processes in plants.

AtGH3.10 is another jasmonic acid-amido synthetase in Arabidopsis thaliana.

Jasmonoyl-isoleucine (JA-Ile) is a key signaling molecule that activates jasmonate-regulated flower development and the wound stress response. For years, JASMONATE RESISTANT1 (JAR1) has been the sole jasmonoyl-amino acid synthetase known to conjugate jasmonic acid (JA) to isoleucine, and the source of persisting JA-Ile in jar1 knockout mutants has remained elusive until now. Here we demonstrate through recombinant enzyme assays and loss-of-function mutant analyses that AtGH3.10 functions as a JA-amido synthetase. Recombinant AtGH3.10 could conjugate JA to isoleucine, alanine, leucine, methionine, and valine. The JA-Ile accumulation in the gh3.10-2 jar1-11 double mutant was nearly eliminated in the leaves and flower buds while its catabolism derivative 12OH-JA-Ile was undetected in the flower buds and unwounded leaves. Residual levels of JA-Ile, JA-Ala, and JA-Val were nonetheless detected in gh3.10-2 jar1-11, suggesting the activities of similar promiscuous enzymes. Upon wounding, the accumulation of JA-Ile and 12OH-JA-Ile and the expression of JA-responsive genes OXOPHYTODIENOIC ACID REDUCTASE3 and JASMONATE ZIM-DOMAIN1 observed in WT, gh3.10-1, and jar1-11 leaves were effectively abolished in gh3.10-2 jar1-11. Additionally, an increased proportion of undeveloped siliques associated with retarded stamen development was observed in gh3.10-2 jar1-11. These findings conclusively show that AtGH3.10 contributes to JA-amino acid biosynthesis and functions partially redundantly with AtJAR1 in sustaining flower development and the wound stress response in Arabidopsis.

Cross-Species Metabolomic Analyses in the Brassicaceae Reveals Common Responses to Ultraviolet-B Exposure.

Exposure to UV-B radiation, an intrinsic component of solar light, is detrimental to all living organisms as chromophore units of DNA, RNA and proteins readily absorb high-energy photons. Indirect damage to the same molecules and lipids is mediated by elevated reactive oxygen species (ROS) levels, a side effect of exposure to UV-B stress. To protect themselves from UV-B radiation, plants produce phytochemical sunscreens, among which flavonoids have shown to be particularly effective. The core aglycone of flavonoid molecules is subjected to chemical decoration, such as glycosylation and acylation, further improving sunscreen properties. In particular, acylation, which adds a phenolic ring to flavonoid molecules, enhances the spectral absorption of UV-A and UV-B rays, providing to this class of compounds exceptional shielding power. In this study, we comprehensively analyzed the responses to UV-B radiation in four Brassicaceae species, including Arabidopsis thaliana, Brassica napus, Brassica oleracea, and Brassica rapa. Our study revealed a complete reprogramming of the central metabolic pathway in response to UV-B radiation characterized by increased production of functional precursors of specialized metabolites with UV-B shielding properties, indicating a targeted effort of plant metabolism to provide increased protection. The analysis of specialized metabolites and transcripts revealed the activation of the phenylpropanoid-acetate pathway, leading to the production of specific classes of flavonoids and a cross-species increase in phenylacylated-flavonoid glucosides with synapoyl glycoside decorations. Interestingly, our analysis also revealed that acyltransferase genes of the class of serine carboxypeptidase-like (SCPLs) proteins are costitutively expressed, but downregulated in response to UV-B radiation, possibly independently of the ELONGATED HYPOCOTYL 5 (HY5) signaling pathway.

A comparative transcriptomics and eQTL approach identifies SlWD40 as a tomato fruit ripening regulator.

Although multiple vital genes with strong effects on the tomato (Solanum lycopersicum) ripening process have been identified via the positional cloning of ripening mutants and cloning of ripening-related transcription factors (TFs), recent studies suggest that it is unlikely that we have fully characterized the gene regulatory networks underpinning this process. Here, combining comparative transcriptomics and expression QTLs, we identified 16 candidate genes involved in tomato fruit ripening and validated them through virus-induced gene silencing analysis. To further confirm the accuracy of the approach, one potential ripening regulator, SlWD40 (WD-40 repeats), was chosen for in-depth analysis. Co-expression network analysis indicated that master regulators such as RIN (ripening inhibitor) and NOR (nonripening) as well as vital TFs including FUL1 (FRUITFUL1), SlNAC4 (NAM, ATAF1,2, and CUC2 4), and AP2a (Activating enhancer binding Protein 2 alpha) strongly co-expressed with SlWD40. Furthermore, SlWD40 overexpression and RNAi lines exhibited substantially accelerated and delayed ripening phenotypes compared with the wild type, respectively. Moreover, transcriptome analysis of these transgenics revealed that expression patterns of ethylene biosynthesis genes, phytoene synthase, pectate lyase, and branched chain amino transferase 2, in SlWD40-RNAi lines were similar to those of rin and nor fruits, which further demonstrated that SlWD40 may act as an important ripening regulator in conjunction with RIN and NOR. These results are discussed in the context of current models of ripening and in terms of the use of comparative genomics and transcriptomics as an effective route for isolating causal genes underlying differences in genotypes.

High-energy-level metabolism and transport occur at the transition from closed to open flowers.

During the maturation phase of flower development, the onset of anthesis visibly marks the transition from buds to open flowers, during which petals stretch out, nectar secretion commences, and pollination occurs. Analysis of the metabolic changes occurring during this developmental transition has primarily focused on specific classes of metabolites, such as pigments and scent emission, and far less on the whole network of primary and secondary metabolites. To investigate the metabolic changes occurring at anthesis, we performed multi-platform metabolomics alongside RNA sequencing in individual florets harvested from the main inflorescence of Arabidopsis (Arabidopsis thaliana) ecotype Col-0. To trace metabolic fluxes at the level of the whole inflorescence and individual florets, we further integrated these studies with radiolabeled experiments. These extensive analyses revealed high-energy-level metabolism and transport of carbohydrates and amino acids, supporting intense metabolic rearrangements occurring at the time of this floral transition. These comprehensive data are discussed in the context of our current understanding of the metabolic shifts underlying flower opening. We envision that this analysis will facilitate the introgression of floral metabolic traits promoting pollination in crop species for which a comprehensive knowledge of flower metabolism is still limited.

Expression Atlas of Selaginella moellendorffii Provides Insights into the Evolution of Vasculature, Secondary Metabolism, and Roots.

Selaginella moellendorffii is a representative of the lycophyte lineage that is studied to understand the evolution of land plant traits such as the vasculature, leaves, stems, roots, and secondary metabolism. However, only a few studies have investigated the expression and transcriptional coordination of Selaginella genes, precluding us from understanding the evolution of the transcriptional programs behind these traits. We present a gene expression atlas comprising all major organs, tissue types, and the diurnal gene expression profiles for S. moellendorffii We show that the transcriptional gene module responsible for the biosynthesis of lignocellulose evolved in the ancestor of vascular plants and pinpoint the duplication and subfunctionalization events that generated multiple gene modules involved in the biosynthesis of various cell wall types. We demonstrate how secondary metabolism is transcriptionally coordinated and integrated with other cellular pathways. Finally, we identify root-specific genes and show that the evolution of roots did not coincide with an increased appearance of gene families, suggesting that the development of new organs does not coincide with increased fixation of new gene functions. Our updated database at conekt.plant.tools represents a valuable resource for studying the evolution of genes, gene families, transcriptomes, and functional gene modules in the Archaeplastida kingdom.

A Chimeric TGA Repressor Slows Down Fruit Maturation and Ripening in Tomato.

The bZIP transcription factor (TF) SlTGA2.2 was previously highlighted as a possible hub in a network regulating fruit growth and transition to ripening (maturation phase). It belongs to a clade of TFs well known for their involvement in the regulation of the salicylic acid-dependent systemic acquired resistance. To investigate if this TGA TF plays a role in tomato fruit growth and maturation, we took advantage of the fruit-specific SlPPC2 promoter (PPC2pro) to target the expression of a SlTGA2.2-SRDX chimeric repressor in a developmental window restricted to early fruit growth and maturation. Here, we show that this SlTGA2.2-SRDX repressor alters early fruit development and metabolism, including chloroplast number and structure, considerably extends the time necessary to reach the mature green stage and slows down fruit ripening. RNA sequencing and plant hormone analyses reveal that PPC2pro:SlTGA2.2-SRDX fruits are maintained in an immature stage as long as PPC2pro is active, through early modifications of plant hormonal signaling and down-regulation of MADS-RIN and NAC-NOR ripening regulators. Once PPC2pro becomes inactive and therefore SlTGA2.2-SRDX expression is reduced, ripening can proceed, albeit at a slower pace than normal. Altogether, this work emphasizes the developmental continuum between fruit growth, maturation and ripening and provides a useful tool to alter and study the molecular bases of tomato fruit transition to ripening.

Kingdom-wide analysis of the evolution of the plant type III polyketide synthase superfamily.

The emergence of type III polyketide synthases (PKSs) was a prerequisite for the conquest of land by the green lineage. Within the PKS superfamily, chalcone synthases (CHSs) provide the entry point reaction to the flavonoid pathway, while LESS ADHESIVE POLLEN 5 and 6 (LAP5/6) provide constituents of the outer exine pollen wall. To study the deep evolutionary history of this key family, we conducted phylogenomic synteny network and phylogenetic analyses of whole-genome data from 126 species spanning the green lineage including Arabidopsis thaliana, tomato (Solanum lycopersicum), and maize (Zea mays). This study thereby combined study of genomic location and context with changes in gene sequences. We found that the two major clades, CHS and LAP5/6 homologs, evolved early by a segmental duplication event prior to the divergence of Bryophytes and Tracheophytes. We propose that the macroevolution of the type III PKS superfamily is governed by whole-genome duplications and triplications. The combined phylogenetic and synteny analyses in this study provide insights into changes in the genomic location and context that are retained for a longer time scale with more recent functional divergence captured by gene sequence alterations.

Sulfur deficiency-induced genes affect seed protein accumulation and composition under sulfate deprivation.

Sulfur deficiency-induced proteins SDI1 and SDI2 play a fundamental role in sulfur homeostasis under sulfate-deprived conditions (-S) by downregulating glucosinolates. Here, we identified that besides glucosinolate regulation under -S, SDI1 downregulates another sulfur pool, the S-rich 2S seed storage proteins in Arabidopsis (Arabidopsis thaliana) seeds. We identified that MYB28 directly regulates 2S seed storage proteins by binding to the At2S4 promoter. We also showed that SDI1 downregulates 2S seed storage proteins by forming a ternary protein complex with MYB28 and MYC2, another transcription factor involved in the regulation of seed storage proteins. These findings have significant implications for the understanding of plant responses to sulfur deficiency.

Quantitative trait loci analysis of seed-specialized metabolites reveals seed-specific flavonols and differential regulation of glycoalkaloid content in tomato.

Given the potential health benefits (and adverse effects), of polyphenolic and steroidal glycoalkaloids in the diet there is a growing interest in fully elucidating the genetic control of their levels in foodstuffs. Here we carried out profiling of the specialized metabolites in the seeds of the Solanum pennellii introgression lines identifying 338 putative metabolite quantitative trait loci (mQTL) for flavonoids, steroidal glycoalkaloids and further specialized metabolites. Two putative mQTL for flavonols and one for steroidal glycoalkaloids were cross-validated by evaluation of the metabolite content of recombinants harboring smaller introgression in the corresponding QTL interval or by analysis of lines from an independently derived backcross inbred line population. The steroidal glycoalkaloid mQTL was localized to a chromosomal region spanning 14 genes, including a previously defined steroidal glycoalkaloid gene cluster. The flavonoid mQTL was further validated via the use of transient and stable overexpression of the Solyc12g098600 and Solyc12g096870 genes, which encode seed-specific uridine 5'-diphosphate-glycosyltransferases. The results are discussed in the context of our understanding of the accumulation of polyphenols and steroidal glycoalkaloids, and how this knowledge may be incorporated into breeding strategies aimed at improving nutritional aspects of plants as well as in fortifying them against abiotic stress.

Natural Variation among Arabidopsis Accessions in the Regulation of Flavonoid Metabolism and Stress Gene Expression by Combined UV Radiation and Cold.

Plants are constantly exposed to stressful environmental conditions. Plant stress reactions were mainly investigated for single stress factors. However, under natural conditions plants may be simultaneously exposed to different stresses. Responses to combined stresses cannot be predicted from the reactions to the single stresses. Flavonoids accumulate in Arabidopsis thaliana during exposure to UV-A, UV-B or cold, but the interactions of these factors on flavonoid biosynthesis were unknown. We therefore investigated the interaction of UV radiation and cold in regulating the expression of well-characterized stress-regulated genes, and on transcripts and metabolites of the flavonoid biosynthetic pathway in 52 natural Arabidopsis accessions that differ widely in their freezing tolerance. The data revealed interactions of cold and UV on the regulation of stress-related and flavonoid biosynthesis genes, and on flavonoid composition. In many cases, plant reactions to a combination of cold and UV were unique under combined stress and not predictable from the responses to the single stresses. Strikingly, all correlations between expression levels of flavonoid biosynthesis genes and flavonol levels were abolished by UV-B exposure. Similarly, correlations between transcript levels of flavonoid biosynthesis genes or flavonoid contents, and freezing tolerance were lost in the presence of UV radiation, while correlations with the expression levels of cold-regulated genes largely persisted. This may indicate different molecular cold acclimation responses in the presence or absence of UV radiation.

The Acetate Pathway Supports Flavonoid and Lipid Biosynthesis in Arabidopsis.

The phenylpropanoid pathway of flavonoid biosynthesis has been the subject of considerable research attention. By contrast, the proposed polyketide pathway, also known as the acetate pathway, which provides malonyl-CoA moieties for the C2 elongation reaction catalyzed by chalcone synthase, is less well studied. Here, we identified four genes as candidates for involvement in the supply of cytosolic malonyl-CoA from the catabolism of acyl-CoA, based on coexpression analysis with other flavonoid-related genes. Two of these genes, ACC and KAT5, have been previously characterized with respect to their involvement in lipid metabolism, but no information concerning their relationship to flavonoid biosynthesis is available. To assess the occurrence and importance of the acetate pathway, we characterized the metabolomes of two mutant or transgenic Arabidopsis lines for each of the four enzymes of this putative pathway using a hierarchical approach covering primary and secondary metabolites as well as lipids. Intriguingly, not only flavonoid content but also glucosinolate content was altered in lines deficient in the acetate pathway, as were levels of lipids and most primary metabolites. We discuss these data in the context of our current understanding of flavonoids and lipid metabolism as well as with regard to improving human nutrition.

Multi-tissue integration of transcriptomic and specialized metabolite profiling provides tools for assessing the common bean (Phaseolus vulgaris) metabolome.

Common bean (Phaseolus vulgaris L.) is an important legume species with a rich natural diversity of landraces that originated from the wild forms following multiple independent domestication events. After the publication of its genome, several resources for this relevant crop have been made available. A comprehensive characterization of specialized metabolism in P. vulgaris, however, is still lacking. In this study, we used a metabolomics approach based on liquid chromatography-mass spectrometry to dissect the chemical composition at a tissue-specific level in several accessions of common bean belonging to different gene pools. Using a combination of literature search, mass spectral interpretation, 13 C-labeling, and correlation analyses, we were able to assign chemical classes and/or putative structures for approximately 39% of all measured metabolites. Additionally, we integrated this information with transcriptomics data and phylogenetic inference from multiple legume species to reconstruct the possible metabolic pathways and identify sets of candidate genes involved in the biosynthesis of specialized metabolites. A particular focus was given to flavonoids, triterpenoid saponins and hydroxycinnamates, as they represent metabolites involved in important ecological interactions and they are also associated with several health-promoting benefits when integrated into the human diet. The data are presented here in the form of an accessible resource that we hope will set grounds for further studies on specialized metabolism in legumes.

Enhancement of vitamin B6 levels in rice expressing Arabidopsis vitamin B6 biosynthesis de novo genes.

Vitamin B6 (pyridoxine) is vital for key metabolic reactions and reported to have antioxidant properties in planta. Therefore, enhancement of vitamin B6 content has been hypothesized to be a route to improve resistance to biotic and abiotic stresses. Most of the current studies on vitamin B6 in plants are on eudicot species, with monocots remaining largely unexplored. In this study, we investigated vitamin B6 biosynthesis in rice, with a view to examining the feasibility and impact of enhancing vitamin B6 levels. Constitutive expression in rice of two Arabidopsis thaliana genes from the vitamin B6 biosynthesis de novo pathway, AtPDX1.1 and AtPDX2, resulted in a considerable increase in vitamin B6 in leaves (up to 28.3-fold) and roots (up to 12-fold), with minimal impact on general growth. Rice lines accumulating high levels of vitamin B6 did not display enhanced tolerance to abiotic stress (salt) or biotic stress (resistance to Xanthomonas oryzae infection). While a significant increase in vitamin B6 content could also be achieved in rice seeds (up to 3.1-fold), the increase was largely due to its accumulation in seed coat and embryo tissues, with little enhancement observed in the endosperm. However, seed yield was affected in some vitamin B6 -enhanced lines. Notably, expression of the transgenes did not affect the expression of the endogenous rice PDX genes. Intriguingly, despite transgene expression in leaves and seeds, the corresponding proteins were only detectable in leaves and could not be observed in seeds, possibly pointing to a mode of regulation in this organ.

Genetic Manipulation of Transcriptional Regulators Alters Nicotine Biosynthesis in Tobacco.

The toxic alkaloid nicotine is produced in the roots of Nicotiana species and primarily accumulates in leaves as a specialized metabolite. A series of metabolic and transport genes involved in the nicotine pathway are coordinately upregulated by a pair of jasmonate-responsive AP2/ERF-family transcription factors, NtERF189 and NtERF199, in the roots of Nicotiana tabacum (tobacco). In this study, we explored the potential of manipulating the expression of these transcriptional regulators to alter nicotine biosynthesis in tobacco. The transient overexpression of NtERF189 led to alkaloid production in the leaves of Nicotiana benthamiana and Nicotiana alata. This ectopic production was further enhanced by co-overexpressing a gene encoding a basic helix-loop-helix-family MYC2 transcription factor. Constitutive and leaf-specific overexpression of NtERF189 increased the accumulation of foliar alkaloids in transgenic tobacco plants but negatively affected plant growth. By contrast, in a knockout mutant of NtERF189 and NtERF199 obtained through CRISPR/Cas9-based genome editing, alkaloid levels were drastically reduced without causing major growth defects. Metabolite profiling revealed the impact of manipulating the nicotine pathway on a wide range of nitrogen- and carbon-containing metabolites. Our findings provide insights into the biotechnological applications of engineering metabolic pathways by targeting transcription factors.