Caitlin A. Kuczynski

Overexpression of WRINKLED1 (WRI1), a master regulator of glycolysis and fatty acid biosynthesis, together with DIACYLGLYCEROL ACYLTRANSFERASE1 (DGAT1), which catalyzes the final step of triacylglycerol assembly, is a promising strategy for enhancing seed oil content. However, how these regulators coordinate system-wide metabolic reprogramming at the levels of gene expression, metabolite pools, and fluxes remains poorly understood. To address this, we performed ¹³C-metabolic flux analysis, metabolomics, and transcriptomics on in vitro cultured pennycress (Thlaspi arvense L.) embryos overexpressing the native WRI1 and DGAT1 homologs. In cultured embryos, WRI1/DGAT1 overexpression increased triacylglycerol accumulation by 28% while reducing protein content by 34%, relative to the wild type. Embryos showed ∼20-fold and 50-fold upregulation of WRI1 and DGAT1 along with induction of WRI1 target genes in glycolysis and fatty acid biosynthesis. Genes associated with photosynthesis and Calvin cycle functions were also upregulated, whereas genes encoding ribosomal proteins and seed storage proteins were strongly repressed, consistent with the observed lipid–protein tradeoff. Flux analysis revealed that enhanced triacylglycerol biosynthesis is supported by increased flux through the Rubisco shunt and cytosolic pyruvate kinase, while the oxidative pentose phosphate pathway and malic enzyme contributed little to NADPH or pyruvate supply. Metabolomic profiling revealed extensive perturbations in glycolytic intermediates, tricarboxylic acid cycle metabolites, and amino acids. In plant grown seeds, WRI1/DGAT1 lines also showed a modest but significant increase in total lipid content. Collectively, these findings reveal how WRI1 and DGAT1 reprogram central metabolism to enhance oil accumulation, with relevance to mature seeds.

During triacylglycerol biosynthesis in developing oilseeds of Arabidopsis thaliana , fatty acid production is regulated by the seed-specific transcription factor WRINKLED1 (WRI1). WRI1 is known to directly stimulate the expression of fatty acid biosynthetic enzymes and a few targets in glycolysis. However, it remains unclear to what extent and how the conversion of sugars into fatty acid biosynthetic precursors in seeds is controlled by WRI1. Based on a previously reported DNA binding motif for WRI1, the ASML1/WRI1 (AW)- box, we developed a comparative genomics approach to search for conserved binding motifs in upstream regions of Arabidopsis thaliana protein-encoding genes and orthologous regions of 11 other Brassicaceae species. The AW-box was over-represented across orthologs for 915 Arabidopsis thaliana genes. Among these, 73 genes with functions in the biosynthesis of fatty acids and triacylglycerols and in glycolysis were enriched. For 90 AW-box sequences associated with these target genes, binding affinity to heterologously expressed Arabidopsis thaliana WRI1 protein was determined using Microscale Thermophoresis. Sites with low dissociation constants are preferentially located close to the transcriptional start site and are highly conserved between the 12 Brassicaceae species. Most of the associated genes were found to be co-expressed with WRI1 during seed development. When 46 automatically and manually curated genes containing conserved AW-sites with high binding affinity are mapped to central metabolism, a conserved regulatory blueprint emerges that infers concerted control of contiguous pathway sections in fatty acid biosynthesis and glycolysis. Among unexpectedly identified putative targets of WRI1 are plastidic fructokinase, phosphoglucose isomerase and several transcription factors. One sentence summary: A combined comparative genomics and in-vitro DNA binding assay approach was used to identify conserved binding sites for the WRINKLED1 transcription factor in central metabolism and lipid biosynthesis.

WRINKLED1 (WRI1) is a transcriptional activator that binds to a conserved sequence (designated as AW box) boxes in the promoters of many genes from central metabolism and fatty acid (FA) synthesis, resulting in their transcription. BIOTIN ATTACHMENT DOMAIN-CONTAINING (BADC) proteins lack a biotin-attachment domain and are therefore inactive, but in the presence of excess FA, BADC1 and BADC3 are primarily responsible for the observed long-term irreversible inhibition of ACETYL-COA CARBOXYLASE, and consequently FA synthesis. Here, we tested the interaction of WRI1 with BADC genes in Arabidopsis (Arabidopsis thaliana) and found purified WRI1 bound with high affinity to canonical AW boxes from the promoters of all three BADC genes. Consistent with this observation, both expression of BADC1, BADC2, and BADC3 genes and BADC1 protein levels were reduced in wri1-1 relative to the wild type, and elevated upon WRI1 overexpression. The double mutant badc1 badc2 phenocopied wri1-1 with respect to both reduction in root length and elevation of indole-3-acetic acid-Asp levels relative to the wild type. Overexpression of BADC1 in wri1-1 decreased indole-3-acetic acid-Asp content and partially rescued its short-root phenotype, demonstrating a role for BADCs in seedling establishment. That WRI1 positively regulates genes encoding both FA synthesis and BADC proteins (i.e. conditional inhibitors of FA synthesis), represents a coordinated mechanism to achieve lipid homeostasis in which plants couple the transcription of their FA synthetic capacity with their capacity to biochemically downregulate it.

WRINKLED1 (WRI1) is a transcriptional activator that binds to AW boxes in the promoters of many genes from central metabolism and FA synthesis, resulting in their transcription. BIOTIN ATTACHMENT DOMAIN-CONTAINING (BADC) proteins are homologs of BIOTIN CARBOXYL CARRIER PROTEIN (BCCP) that lack a biotin-attachment domain and are therefore inactive. In the presence of excess FA, BADC1 and BADC3 are primarily responsible for the observed long-term irreversible inhibition of ACETYL-COA CARBOXYLASE (ACCase), and consequently FA synthesis. Purified WRI1 bound with high affinity ( Kd s in the low nanomolar range) to canonical AW-boxes from the promoters of all three BADC genes. Consistent with this observation, the expression of BADC1, BADC2 and BADC3 genes and BADC1 protein levels were reduced in wri1-1 relative to wild type (WT), and BADCs gene expression and BADC1 protein levels also were elevated upon WRI1 overexpression. The double mutant badc1badc2 phenocopied wri1-1 with respect to both reduction in root length, and elevation of indole-3-acetic acid-Asp (IAA-Asp) levels relative to WT. Overexpression of BADC1 in wri1-1 decreased its IAA-Asp and partially rescued its short-root phenotype demonstrating a role for BADCs in seedling establishment. That WRI1 positively regulates genes encoding both FA synthesis and BADCs i.e., conditional inhibitors of FA synthesis, represents a coordinated mechanism to achieve lipid homeostasis in which plants couple the transcription their FA synthetic capacity with their capacity to biochemically downregulate it. One sentence summary: WRI1 regulates genes encoding both fatty acid synthesis and inhibitors of FA synthesis (BADCs), creating a lipid homeostatic mechanism in which the transcription of FA synthetic capacity is coordinated with the capacity to biochemically downregulate FA synthesis.

In this article, we couple the geographic variation in 127 single-nucleotide polymorphism (SNP) frequencies in genes of 46 enzymes of central metabolism with their associated cis-expression variation to predict latitudinal or climatic-driven gene expression changes in the metabolic architecture of Drosophila melanogaster. Forty-two percent of the SNPs in 65% of the genes show statistically significant clines in frequency with latitude across the 20 local population samples collected from southern Florida to Ontario. A number of SNPs in the screened genes are also associated with significant expression variation within the Raleigh population from North Carolina. A principal component analysis of the full variance–covariance matrix of latitudinal changes in SNP-associated standardized gene expression allows us to identify those major genes in the pathway and its associated branches that are likely targets of natural selection. When embedded in a central metabolic context, we show that these apparent targets are concentrated in the genes of the upper glycolytic pathway and pentose shunt, those controlling glycerol shuttle activity, and finally those enzymes associated with the utilization of glutamate and pyruvate. These metabolites possess high connectivity and thus may be the points where flux balance can be best shifted. We also propose that these points are conserved points associated with coupling energy homeostasis and energy sensing in mammals. We speculate that the modulation of gene expression at specific points in central metabolism that are associated with shifting flux balance or possibly energy-state sensing plays a role in adaptation to climatic variation.

Cosmopolitan populations of Drosophila melanogaster have co-opted a form of reproductive diapause to overwinter in northern populations. Polymorphism in the couch potato gene has been implicated in genetic variation for this diapause trait. Using a collection of 20 populations from Florida to Canada and 11 collections from 3 years in a Pennsylvania orchard, we estimated the allele frequencies for 15 single nucleotide polymorphisms (SNPs) in the couch potato gene. These include the specific polymorphism associated with diapause inducability. We find that the SNP polymorphism, 48034(A/T), is correlated with latitude and its frequencies are predicted by the incidence of diapause trait. We find that the clinal patterns for cpo SNPs sampled in 1997 are similar to the same SNPs sampled in 2009-2010. SNPs that show apparent associations with cpo expression are also clinal with the low-expression allele increasing in frequency, as would be predicted from functional knockout studies of cpo. Finally, we see a significant pattern where the frequency of the diapause-causing allele drops in frequency during the summer season, consistent with the drop in the incidence of the diapause trait. The selection required to drive this response is large, roughly 24% to 59% per generation depending on the degree of dominance.

Speciation is often categorized based on geographic modes (allopatric, parapatric or sympatric). Although it is widely accepted that species can arise in allopatry and then later become sympatrically or parapatrically distributed, patterns in the opposite direction are also theoretically possible (e.g. sympatric lineages or ecotypes becoming parapatric), but such patterns have not been shown at a macrogeographic scale. Here, we analyse genetic, climatic, ecological and morphological data and show that two typically sympatric colour morphs of the salamander Plethodon cinereus (redback and leadback) appear to have become parapatrically distributed on Long Island, New York, with pure‐redback populations in the west and pure‐leadback populations in the east (and polymorphic populations in between and on the mainland). In addition, the pure‐leadback populations in eastern Long Island are genetically, ecologically and morphologically divergent from both mainland and other Long Island populations, suggesting the possibility of incipient speciation. This parapatric separation seems to be related to the different ecological preferences of the two morphs, preferences which are present on the mainland and across Long Island. These results potentially support the idea that spatial segregation of sympatric ecotypes may sometimes play an important part in parapatric speciation.

Molecular data offer great potential to resolve the phylogeny of living taxa but can molecular data improve our understanding of relationships of fossil taxa? Simulations suggest that this is possible, but few empirical examples have demonstrated the ability of molecular data to change the placement of fossil taxa. We offer such an example here. We analyze the placement of snakes among squamate reptiles, combining published morphological data (363 characters) and new DNA sequence data (15,794 characters, 22 nuclear loci) for 45 living and 19 fossil taxa. We find several intriguing results. First, some fossil taxa undergo major changes in their phylogenetic position when molecular data are added. Second, most fossil taxa are placed with strong support in the expected clades by the combined data Bayesian analyses, despite each having >98% missing cells and despite recent suggestions that extensive missing data are problematic for Bayesian phylogenetics. Third, morphological data can change the placement of living taxa in combined analyses, even when there is an overwhelming majority of molecular characters. Finally, we find strong but apparently misleading signal in the morphological data, seemingly associated with a burrowing lifestyle in snakes, amphisbaenians, and dibamids. Overall, our results suggest promise for an integrated and comprehensive Tree of Life by combining molecular and morphological data for living and fossil taxa.