skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Integrative analysis of the shikonin metabolic network identifies new gene connections and reveals evolutionary insight into shikonin biosynthesis
Summary Plant specialized 1,4-naphthoquinones present a remarkable case of convergent evolution. Species across multiple discrete orders of vascular plants produce diverse 1,4-naphthoquinones via one of several pathways using different metabolic precursors. Evolution of these pathways was preceded by events of metabolic innovation and many appear to share connections with biosynthesis of photosynthetic or respiratory quinones. Here, we sought to shed light on the metabolic connections linking shikonin biosynthesis with its precursor pathways and on the origins of shikonin metabolic genes. Downregulation of Lithospermum erythrorhizon geranyl diphosphate synthase (LeGPPS), recently shown to have been recruited from a cytoplasmic farnesyl diphosphate synthase (FPPS), resulted in reduced shikonin production and a decrease in expression of mevalonic acid and phenylpropanoid pathway genes. Next, we used LeGPPS and other known shikonin pathway genes to build a coexpression network model for identifying new gene connections to shikonin metabolism. Integrative in silico analyses of network genes revealed candidates for biochemical steps in the shikonin pathway arising from Boraginales-specific gene family expansion. Multiple genes in the shikonin coexpression network were also discovered to have originated from duplication of ubiquinone pathway genes. Taken together, our study provides evidence for transcriptional crosstalk between shikonin biosynthesis and its precursor pathways, identifies several shikonin pathway gene candidates and their evolutionary histories, and establishes additional evolutionary links between shikonin and ubiquinone metabolism. Moreover, we demonstrate that global coexpression analysis using limited transcriptomic data obtained from targeted experiments is effective for identifying gene connections within a defined metabolic network.  more » « less
Award ID(s):
1831493
PAR ID:
10343552
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Horticulture Research
Volume:
9
ISSN:
2052-7276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Goldman, Gustavo H. (Ed.)
    ABSTRACT Fungal secondary metabolites are widely used as therapeutics and are vital components of drug discovery programs. A major challenge hindering discovery of novel secondary metabolites is that the underlying pathways involved in their biosynthesis are transcriptionally silent under typical laboratory growth conditions, making it difficult to identify the transcriptional networks that they are embedded in. Furthermore, while the genes participating in secondary metabolic pathways are typically found in contiguous clusters on the genome, known as biosynthetic gene clusters (BGCs), this is not always the case, especially for global and pathway-specific regulators of pathways’ activities. To address these challenges, we used 283 genome-wide gene expression data sets of the ascomycete cell factory Aspergillus niger generated during growth under 155 different conditions to construct two gene coexpression networks based on Spearman’s correlation coefficients (SCCs) and on mutual rank-transformed Pearson’s correlation coefficients (MR-PCCs). By mining these networks, we predicted six transcription factors, named MjkA to MjkF, to regulate secondary metabolism in A. niger . Overexpression of each transcription factor using the Tet-On cassette modulated the production of multiple secondary metabolites. We found that the SCC and MR-PCC approaches complemented each other, enabling the delineation of putative global (SCC) and pathway-specific (MR-PCC) transcription factors. These results highlight the potential of coexpression network approaches to identify and activate fungal secondary metabolic pathways and their products. More broadly, we argue that drug discovery programs in fungi should move beyond the BGC paradigm and focus on understanding the global regulatory networks in which secondary metabolic pathways are embedded. IMPORTANCE There is an urgent need for novel bioactive molecules in both agriculture and medicine. The genomes of fungi are thought to contain vast numbers of metabolic pathways involved in the biosynthesis of secondary metabolites with diverse bioactivities. Because these metabolites are biosynthesized only under specific conditions, the vast majority of the fungal pharmacopeia awaits discovery. To discover the genetic networks that regulate the activity of secondary metabolites, we examined the genome-wide profiles of gene activity of the cell factory Aspergillus niger across hundreds of conditions. By constructing global networks that link genes with similar activities across conditions, we identified six putative global and pathway-specific regulators of secondary metabolite biosynthesis. Our study shows that elucidating the behavior of the genetic networks of fungi under diverse conditions harbors enormous promise for understanding fungal secondary metabolism, which ultimately may lead to novel drug candidates. 
    more » « less
  2. The evolution of developmental rates may drive morphological change. Caenorhabditis inopinata develops nearly twice as slowly as Caenorhabditis elegans. clk-1 encodes a hydroxylase required for synthesizing ubiquinone, and mutant clk-1 slow growth phenotypes can be rescued by supplying animals with a ubiquinone precursor analogue, 2,4-dihydroxybenzoate. RNA-seq data showing low clk-1 expression raised the possibility that C. inopinata grows slowly because of reduced ubiquinone biosynthesis. C. inopinata did not reveal a clear reduction in the age of maturation when reared on 2,4-dihydroxybenzoate. Further scrutiny of RNA-seq results revealed multiple ubiquinone metabolism genes have low expression in C. inopinata. Divergent clk-1 expression alone may not be a major driver of the evolution of slow development in this species. 
    more » « less
  3. SUMMARY Zea mays(maize) makes phytoalexins such as sesquiterpenoid zealexins, to combat invading pathogens. Zealexins are produced from farnesyl diphosphate in microgram per gram fresh weight quantities. As farnesyl diphosphate is also a precursor for many compounds essential for plant growth, the question arises as to howZ. maysproduces high levels of zealexins without negatively affecting vital plant systems. To examine if specific pools of farnesyl diphosphate are made for zealexin synthesis we made CRISPR/Cas9 knockouts of each of the three farnesyl diphosphate synthases (FPS) inZ. maysand examined the resultant impacts on different farnesyl diphosphate‐derived metabolites. We found that FPS3 (GRMZM2G098569) produced most of the farnesyl diphosphate for zealexins, while FPS1 (GRMZM2G168681) made most of the farnesyl diphosphate for the vital respiratory co‐factor ubiquinone. Indeed,fps1mutants had strong developmental phenotypes such as reduced stature and development of chlorosis. The replication and evolution of thefpsgene family inZ. maysenabled it to produce dedicated FPSs for developmentally related ubiquinone production (FPS1) or defense‐related zealexin production (FPS3). This partitioning of farnesyl diphosphate production between growth and defense could contribute to the ability ofZ. maysto produce high levels of phytoalexins without negatively impacting its growth. 
    more » « less
  4. Extremophilic yeasts have favorable metabolic and tolerance traits for biomanufacturing‐ like lipid biosynthesis, flavinogenesis, and halotolerance – yet the connection between these favorable phenotypes and strain genotype is not well understood. To this end, this study compares the phenotypes and gene expression patterns of biotechnologically relevant yeasts Yarrowia lipolytica, Debaryomyces hansenii, and Debaryomyces subglobosus grown under nitrogen starvation, iron starvation, and salt stress. To analyze the large data set across species and conditions, two approaches were used: a “network‐first” approach where a generalized metabolic network serves as a scaffold for mapping genes and a “cluster‐first” approach where unsupervised machine learning co‐expression analysis clusters genes. Both approaches provide insight into strain behavior. The network‐first approach corroborates that Yarrowia upregulates lipid biosynthesis during nitrogen starvation and provides new evidence that riboflavin overproduction in Debaryomyces yeasts is overflow metabolism that is routed to flavin cofactor production under salt stress. The cluster‐first approach does not rely on annotation; therefore, the coexpression analysis can identify known and novel genes involved in stress responses, mainly transcription factors and transporters. Therefore, this work links the genotype to the phenotype of biotechnologically relevant yeasts and demonstrates the utility of complementary computational approaches to gain insight from transcriptomics data across species and conditions. 
    more » « less
  5. SUMMARY The stilbenoid pathway is responsible for the production of resveratrol in grapevine (Vitis viniferaL.). A few transcription factors (TFs) have been identified as regulators of this pathway but the extent of this control has not been deeply studied. Here we show how DNA affinity purification sequencing (DAP‐Seq) allows for the genome‐wide TF‐binding site interrogation in grape. We obtained 5190 and 4443 binding events assigned to 4041 and 3626 genes for MYB14 and MYB15, respectively (approximately 40% of peaks located within −10 kb of transcription start sites). DAP‐Seq of MYB14/MYB15 was combined with aggregate gene co‐expression networks (GCNs) built from more than 1400 transcriptomic datasets from leaves, fruits, and flowers to narrow down bound genes to a set of high confidence targets. The analysis of MYB14, MYB15, and MYB13, a third uncharacterized member of Subgroup 2 (S2), showed that in addition to the few previously known stilbene synthase (STS) targets, these regulators bind to 30 of 47STSfamily genes. Moreover, all three MYBs bind to severalPAL,C4H, and4CLgenes, in addition to shikimate pathway genes, theWRKY03stilbenoid co‐regulator and resveratrol‐modifying gene candidates among which ROMT2‐3 were validated enzymatically. A high proportion of DAP‐Seq bound genes were induced in the activated transcriptomes of transientMYB15‐overexpressing grapevine leaves, validating our methodological approach for delimiting TF targets. Overall, Subgroup 2 R2R3‐MYBs appear to play a key role in binding and directly regulating several primary and secondary metabolic steps leading to an increased flux towards stilbenoid production. The integration of DAP‐Seq and reciprocal GCNs offers a rapid framework for gene function characterization using genome‐wide approaches in the context of non‐model plant species and stands up as a valid first approach for identifying gene regulatory networks of specialized metabolism. 
    more » « less