Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gene regulatory network”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Adaptive laboratory evolution and metabolic engineering of Cupriavidus necator for improved catabolism of volatile fatty acids

Bioconversion of high-volume waste streams into value-added products will be an integral component of the growing bioeconomy. Volatile fatty acids (VFAs) (e.g., butyrate, valerate, and hexanoate) are an emerging and promising waste-derived feedstock for microbial carbon upcycling. Cupriavidus necator H16 is a favorable host for conversion of VFAs into various bioproducts due to its diverse carbon metabolism, ease of metabolic engineering, and use at industrial scales. Here, in this study, we report that a common strategy to improve product titers in C. necator, deletion of the polyhydroxybutyrate (PHB) biosynthetic operon, results in a significant growth defect on VFA substrates. Using adaptive laboratory evolution, we identify mutations to the regulator gene phaR, the two-component response regulator-histidine kinase pair encoded by H16_A1372/H16_A1373, and the tripartite transporter assembly encoded by H16_A2296-A2298 as causative for improved growth on VFA substrates. Deletion of phaR and H16_A1373 led to significantly reduced NADH abundance accompanied by large changes to expression of genes involved in carbon metabolism, balance of electron carriers, and oxidative stress tolerance that may be responsible for improved growth of these engineered strains. These results provide insight into the role of PHB biosynthesis in carbon and energy metabolism and highlight a key role for the regulator PhaR in global regulatory networks. By combining mutations, we generated platform strains with significant growth improvements on VFAs, which can enable improved conversion of waste-derived VFA substrates to target bioproducts.

09 BIOMASS FUELS↗
DOE OSTI2024

maranasgroup/yeastsMFA

The parameterization of kinetic models requires measurement of fluxes and/or metabolite levels for a base strain and a few genetic perturbations thereof. Unlike stoichiometric models that are mostly invariant to the specific strain, it remains unclear whether kinetic models constructed for different strains of the same species have similar or significantly different kinetic parameters. This important question underpins the applicability range and prediction limits of kinetic reconstructions. To this end, herein we parameterize two separate large-scale kinetic models using K-FIT with genome-wide coverage corresponding to two distinct strains of Saccharomyces cerevisiae: CEN.PK 113-7D strain (model k-sacce306-CENPK), and growth-deficient BY4741 (isogenic to S288c; model k-sacce306-BY4741). The metabolic network for each model contains 306 reactions, 230 metabolites, and 119 substrate-level regulatory interactions. The two models (for CEN.PK and BY4741) recapitulate, within one standard deviation, 77% and 75% of the fitted dataset fluxes, respectively, determined by 13C metabolic flux analysis for wild-type and eight single-gene knockout mutants of each strain. Strain-specific kinetic parameterization results indicate that key enzymes in the TCA cycle, glycolysis, and arginine and proline metabolism drive the metabolic differences between these two strains of S. cerevisiae. Our results suggest that although kinetic models cannot be readily used across strains as stoichiometric models, they can capture species-specific information through the kinetic parameterization process.

Dinh, Hoang↗
DOE OSTI2023

maranasgroup/k-sacce306

The parameterization of kinetic models requires measurement of fluxes and/or metabolite levels for a base strain and a few genetic perturbations thereof. Unlike stoichiometric models that are mostly invariant to the specific strain, it remains unclear whether kinetic models constructed for different strains of the same species have similar or significantly different kinetic parameters. This important question underpins the applicability range and prediction limits of kinetic reconstructions. To this end, herein we parameterize two separate large-scale kinetic models using K-FIT with genome-wide coverage corresponding to two distinct strains of Saccharomyces cerevisiae: CEN.PK 113-7D strain (model k-sacce306-CENPK), and growth-deficient BY4741 (isogenic to S288c; model k-sacce306-BY4741). The metabolic network for each model contains 306 reactions, 230 metabolites, and 119 substrate-level regulatory interactions. The two models (for CEN.PK and BY4741) recapitulate, within one standard deviation, 77% and 75% of the fitted dataset fluxes, respectively, determined by 13C metabolic flux analysis for wild-type and eight single-gene knockout mutants of each strain. Strain-specific kinetic parameterization results indicate that key enzymes in the TCA cycle, glycolysis, and arginine and proline metabolism drive the metabolic differences between these two strains of S. cerevisiae. Our results suggest that although kinetic models cannot be readily used across strains as stoichiometric models, they can capture species-specific information through the kinetic parameterization process.

Dinh, Hoang↗
DOE OSTI2023

Comparative study of two Saccharomyces cerevisiae strains with kinetic models at genome-scale

The parameterization of kinetic models requires measurement of fluxes and/or metabolite levels for a base strain and a few genetic perturbations thereof. Unlike stoichiometric models that are mostly invariant to the specific strain, it remains unclear whether kinetic models constructed for different strains of the same species have similar or significantly different kinetic parameters. This important question underpins the applicability range and prediction limits of kinetic reconstructions. To this end, herein we parameterize two separate large-scale kinetic models using K-FIT with genome-wide coverage corresponding to two distinct strains of Saccharomyces cerevisiae: CEN.PK 113-7D strain (model k-sacce306-CENPK), and growth-deficient BY4741 (isogenic to S288c; model k-sacce306-BY4741). The metabolic network for each model contains 306 reactions, 230 metabolites, and 119 substrate-level regulatory interactions. The two models (for CEN.PK and BY4741) recapitulate, within one standard deviation, 77% and 75% of the fitted dataset fluxes, respectively, determined by 13 C metabolic flux analysis for wild-type and eight single-gene knockout mutants of each strain. Strain-specific kinetic parameterization results indicate that key enzymes in the TCA cycle, glycolysis, and arginine and proline metabolism drive the metabolic differences between these two strains of S. cerevisiae. Our results suggest that although kinetic models cannot be readily used across strains as stoichiometric models, they can capture species-specific information through the kinetic parameterization process.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Spurious regulatory connections dictate the expression-fitness landscape of translation factors

During steady-state cell growth, individual enzymatic fluxes can be directly inferred from growth rate by mass conservation, but the inverse problem remains unsolved. Perturbing the flux and expression of a single enzyme could have pleiotropic effects that may or may not dominate the impact on cell fitness. Here, we quantitatively dissect the molecular and global responses to varied expression of translation termination factors (peptide release factors, RFs) in the bacterium Bacillus subtilis. While endogenous RF expression maximizes proliferation, deviations in expression lead to unexpected distal regulatory responses that dictate fitness reduction. Molecularly, RF depletion causes expression imbalance at specific operons, which activates master regulators and detrimentally overrides the transcriptome. Through these spurious connections, RF abundances are thus entrenched by focal points within the regulatory network, in one case located at a single stop codon. Such regulatory entrenchment suggests that predictive bottom-up models of expression-fitness landscapes will require near-exhaustive characterization of parts.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Relevance of the Unfolded Protein Response to Spaceflight-Induced Transcriptional Reprogramming in Arabidopsis

Plants are primary producers of food and oxygen on Earth and will likewise be indispensable to the establishment of large-scale sustainable ecosystems and human survival in space. To contribute to the understanding of how plants respond to spaceflight stress, we examined the significance of the unfolded protein response (UPR), a conserved signaling cascade that responds to a number of unfavorable environmental stresses, in the model plant Arabidopsis thaliana. To do so, we performed a large-scale comparative transcriptome profiling in wild type and various UPR-defective mutants during the SpaceX-CRS12 mission to the International Space Station. We established that orbital culture substantially alters the expression of hundreds of stress-related genes compared with ground control conditions. Although expression of those genes varied in the UPR mutants on the ground, it was largely similar across the genotypes in the spaceflight condition. Furthermore, our results have yielded new information on how plants respond to growth in orbit and support the hypothesis that spaceflight induces the activation of signaling pathways that compensate for the loss of UPR regulators in the control of downstream transcriptional regulatory networks.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2000

Diversity and regulation of plant Ca2+ pumps: insights from expression in yeast

The spatial and temporal regulation of calcium concentration in plant cells depends on the coordinate activities of channels and active transporters located on different organelles and membranes. Several Ca2+ pumps have been identified and characterized by functional expression of plant genes in a yeast mutant (K616). This expression system has opened the way to a genetic and biochemical characterization of the regulatory and catalytic features of diverse Ca2+ pumps. Plant Ca(2+)-ATPases fall into two major types: AtECA1 represents one of four or more members of the type IIA (ER-type) Ca(2+)-ATPases in Arabidopsis, and AtACA2 is one of seven or more members of the type IIB (PM-type) Ca(2+)-ATPases that are regulated by a novel amino terminal domain. Type IIB pumps are widely distributed on membranes, including the PM (plasma membrane), vacuole, and ER (endoplasmic reticulum). The regulatory domain serves multiple functions, including autoinhibition, calmodulin binding, and sites for modification by phosphorylation. This domain, however, is considerably diverse among several type IIB ATPases, suggesting that the pumps are differentially regulated. Understanding of Ca2+ transporters at the molecular level is providing insights into their roles in signaling networks and in regulating fundamental processes of cell biology.

Review, Academic↗
DOE OSTI2026

Lignin’s Indispensable role in orchestrating seed stone formation: Insights from jujuba, peach and pear with future prospective on pitless fruits

A hard endocarp (i.e., stone) inside fruit is a characteristic of drupe fruits such as jujube, peach, mango, etc. Hard stone significantly affects the quality and downstream processing of fruits. The complex aromatic polymer lignin deposition in the secondary cell wall determines stone hardness. Lignin comprises phenylpropanoid units formed by hydroxycinnamoyl alcohol, which includes coniferyl, sinapyl, and p-coumaroyl alcohols. Lignin biosynthesis pathway involves a series of complex enzymatic reactions initiated from phenylalanine ammonia lyase and ends up polymerizing lignin monomers by laccase and peroxidase enzymes. Phytohormones, particularly auxin, gibberellins, and Ca²⁺ signaling, further modulate endocarp lignification by regulating transcriptional networks and lignin biosynthetic genes, thereby fine-tuning secondary cell wall thickening and stone hardness in drupe fruits. Lignin biosynthesis is controlled by both structural genes and transcriptional regulators. The structural genes encoding lignin biosynthetic enzymes include LAC12–1, PAL2, C4H, C3H, CSE, CCoAOMT, F5H, CAD, and PRX1. In addition, several transcription factors regulating secondary cell wall and lignin deposition, such as MYB24, bZIP48, and bZIP33 play key regulatory roles. Conversely, delignification or suppression of stone formation is associated with transcription factors (Pistillata, MYB32, FUL, and REPLUMLESS) and post-transcriptional regulators, including miR397a, miR31-3p, and miR8-5p. Accurate alteration in the expression of these genes will result in the attainment of stoneless fruits for cheap and hazel-free downstream processing.

Fruit endocarp↗
DOE OSTI2021

Heterospecific Neighbor Plants Impact Root Microbiome Diversity and Molecular Function of Root Fungi

Within the forest community, competition and facilitation between adjacent-growing conspecific and heterospecific plants are mediated by interactions involving common mycorrhizal networks. The ability of plants to alter their neighbor’s microbiome is well documented, but the molecular biology of plant-fungal interactions during competition and facilitation has not been previously examined. We used a common soil-plant bioassay experiment to study molecular plant-microbial interactions among rhizosphere communities associated with Pinus taeda (native host) and Populus trichocarpa (non-native host). Gene expression of interacting fungal and bacterial rhizosphere communities was compared among three plant-pairs: Populus growing with Populus , Populus with Pinus , and Pinus with Pinus . Our results demonstrate that heterospecific plant partners affect the assembly of root microbiomes, including the changes in the structure of host specific community. Comparative metatranscriptomics reveals that several species of ectomycorrhizal fungi (EMF) and saprotrophic fungi exhibit different patterns of functional and regulatory gene expression with these two plant hosts. Heterospecific plants affect the transcriptional expression pattern of EMF host-specialists (e.g., Pinus -associated Suillus spp.) on both plant species, mainly including the genes involved in the transportation of amino acids, carbohydrates, and inorganic ions. Alteration of root microbiome by neighboring plants may help regulate basic plant physiological processes via modulation of molecular functions in the root microbiome.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1999

Early animal evolution: emerging views from comparative biology and geology

The Cambrian appearance of fossils representing diverse phyla has long inspired hypotheses about possible genetic or environmental catalysts of early animal evolution. Only recently, however, have data begun to emerge that can resolve the sequence of genetic and morphological innovations, environmental events, and ecological interactions that collectively shaped Cambrian evolution. Assembly of the modern genetic tool kit for development and the initial divergence of major animal clades occurred during the Proterozoic Eon. Crown group morphologies diversified in the Cambrian through changes in the genetic regulatory networks that organize animal ontogeny. Cambrian radiation may have been triggered by environmental perturbation near the Proterozoic-Cambrian boundary and subsequently amplified by ecological interactions within reorganized ecosystems.

Non-NASA Center↗
DOE OSTI2023

Regulatory response to a hybrid ancestral nitrogenase in Azotobacter vinelandii

Biological nitrogen fixation, the microbial reduction of atmospheric nitrogen to bioavailable ammonia, represents both a major limitation on biological productivity and a highly desirable engineering target for synthetic biology. However, the engineering of nitrogen fixation requires an integrated understanding of how the gene regulatory dynamics of host diazotrophs respond across sequence-function space of its central catalytic metalloenzyme, nitrogenase. Here, we interrogate this relationship by analyzing the transcriptome of Azotobacter vinelandii engineered with a phylogenetically inferred ancestral nitrogenase protein variant. The engineered strain exhibits reduced cellular nitrogenase activity but recovers wild-type growth rates following an extended lag period. We find that expression of genes within the immediate nitrogen fixation network is resilient to the introduced nitrogenase sequence-level perturbations. Rather the sustained physiological compatibility with the ancestral nitrogenase variant is accompanied by reduced expression of genes that support trace metal and electron resource allocation to nitrogenase. Our results spotlight gene expression changes in cellular processes adjacent to nitrogen fixation as productive engineering considerations to improve compatibility between remodeled nitrogenase proteins and engineered host diazotrophs.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2023

Regulatory response to a hybrid ancestral nitrogenase in Azotobacter vinelandii

Biological nitrogen fixation, the microbial reduction of atmospheric nitrogen to bioavailable ammonia, represents both a major limitation on biological productivity and a highly desirable engineering target for synthetic biology. However, the engineering of nitrogen fixation requires an integrated understanding of how the gene regulatory dynamics of host diazotrophs respond across sequence-function space of its central catalytic metalloenzyme, nitrogenase. Here, we interrogate this relationship by analyzing the transcriptome of Azotobacter vinelandii engineered with a phylogenetically inferred ancestral nitrogenase protein variant. The engineered strain exhibits reduced cellular nitrogenase activity but recovers wild-type growth rates following an extended lag period. We find that expression of genes within the immediate nitrogen fixation network is resilient to the introduced nitrogenase sequence-level perturbations. Rather the sustained physiological compatibility with the ancestral nitrogenase variant is accompanied by reduced expression of genes that support trace metal and electron resource allocation to nitrogenase. Our results spotlight gene expression changes in cellular processes adjacent to nitrogen fixation as productive engineering considerations to improve compatibility between remodeled nitrogenase proteins and engineered host diazotrophs.

nitrogen fixation↗
DOE OSTI2025

Physiological roles of an Acinetobacter -specific σ factor

ABSTRACT The Gram-negative pathogen Acinetobacter baumannii is considered an “urgent threat” to human health due to its propensity to become antibiotic resistant. Understanding the distinct regulatory paradigms used by A. baumannii to mitigate cellular stresses may uncover new therapeutic targets. Many γ-proteobacteria use the extracytoplasmic function (ECF) σ factor, RpoE, to invoke envelope homeostasis networks in response to stress. Acinetobacter species contain the poorly characterized ECF “SigAb”; however, it is unclear if SigAb has the same physiological role as RpoE. Here, we show that SigAb is a metal stress-responsive ECF that appears unique to Acinetobacter species and distinct from RpoE-like ECFs. We combine promoter mutagenesis, motif scanning, and chromatin immunoprecipitation-sequencing (ChIP-seq) to define the direct SigAb regulon, which consists of genes encoding SigAb itself, the stringent response mediator, RelA, and the uncharacterized small RNA, “SabS.” However, RNA-seq of strains overexpressing SigAb revealed a large, indirect regulon containing hundreds of genes. Metal resistance genes are key elements of the indirect regulon, as CRISPRi knockdown of sigAb or sabS resulted in increased copper sensitivity and excess copper-induced SigAb-dependent transcription. Furthermore, we found that two uncharacterized genes in the sigAb operon, “ aabA ” and “ aabB ,” have anti-SigAb activity. Finally, employing a targeted Tn-seq approach that uses CRISPR-associated transposons, we show that sigAb , aabA , and aabB are important for fitness even during optimal growth conditions. Our work reveals new physiological roles for SigAb and SabS, provides a novel approach for assessing gene fitness, and highlights the distinct regulatory architecture of A. baumannii . IMPORTANCE Acinetobacter baumannii is a hospital-acquired pathogen, and many strains are resistant to multiple antibiotics. Understanding how A. baumannii senses and responds to stress may uncover novel routes to treat infections. Here, we examine how the Acinetobacter -specific transcription factor, SigAb, mitigates stress. We find that SigAb directly regulates only a small number of genes, but indirectly controls hundreds of genes that have substantial impacts on cell physiology. We show that SigAb is required for maximal growth, even during optimal conditions, and is acutely required during growth in the presence of elevated copper. Given that copper toxicity plays roles in pathogenesis and on copper-containing surfaces in hospitals, we speculate that SigAb function may be important in clinically relevant contexts.

Bacon, Emily E. (ORCID:0000000180907689)↗
DOE OSTI2023

TroR is the primary regulator of the iron homeostasis transcription network in the halophilic archaeon Haloferax volcanii

Abstract Maintaining the intracellular iron concentration within the homeostatic range is vital to meet cellular metabolic needs and reduce oxidative stress. Previous research revealed that the haloarchaeon Halobacterium salinarum encodes four diphtheria toxin repressor (DtxR) family transcription factors (TFs) that together regulate the iron response through an interconnected transcriptional regulatory network (TRN). However, the conservation of the TRN and the metal specificity of DtxR TFs remained poorly understood. Here we identified and characterized the TRN of Haloferax volcanii for comparison. Genetic analysis demonstrated that Hfx. volcanii relies on three DtxR transcriptional regulators (Idr, SirR, and TroR), with TroR as the primary regulator of iron homeostasis. Bioinformatics and molecular approaches revealed that TroR binds a conserved cis-regulatory motif located ∼100 nt upstream of the start codon of iron-related target genes. Transcriptomics analysis demonstrated that, under conditions of iron sufficiency, TroR repressed iron uptake and induced iron storage mechanisms. TroR repressed the expression of one other DtxR TF, Idr. This reduced DtxR TRN complexity relative to that of Hbt. salinarum appeared correlated with natural variations in iron availability. Based on these data, we hypothesize that variable environmental conditions such as iron availability appear to select for increasing TRN complexity.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2017

Systemic Microgravity Response: Utilizing GeneLab to Develop Hypotheses for Spaceflight Risks

Biological risks associated with microgravity is a major concern for space travel. Although determination of risk has been a focus for NASA research, data examining systemic (i.e., multi- or pan-tissue) responses to space flight are sparse. The overall goal of our work is to identify potential master regulators responsible for such responses to microgravity conditions. To do this we utilized the NASA GeneLab database which contains a wide array of omics experiments, including data from: 1) different flight conditions (space shuttle (STS) missions vs. International Space Station (ISS); 2) different tissues; and 3) different types of assays that measure epigenetic, transcriptional, and protein expression changes. We have performed meta-analysis identifying potential master regulators involved with systemic responses to microgravity. The analysis used 7 different murine and rat data sets, examining the following tissues: liver, kidney, adrenal gland, thymus, mammary gland, skin, and skeletal muscle (soleus, extensor digitorum longus, tibialis anterior, quadriceps, and gastrocnemius). Using a systems biology approach, we were able to determine that p53 and immune related pathways appear central to pan-tissue microgravity responses. Evidence for a universal response in the form of consistency of change across tissues in regulatory pathways was observed in both STS and ISS experiments with varying durations; while degree of change in expression of these master regulators varied across species and strain, some change in these master regulators was universally observed. Interestingly, certain skeletal muscle (gastrocnemius and soleus) show an overall down-regulation in these genes, while in other types (extensor digitorum longus, tibialis anterior and quadriceps) they are up-regulated, suggesting certain muscle tissues may be compensating for atrophy responses caused by microgravity. Studying these organtissue-specific perturbations in molecular signaling networks, we demonstrate the value of GeneLab in characterizing potential master regulators associated with biological risks for spaceflight.

Microgravity↗
DOE OSTI2022

Machine Learning of All Mycobacterium tuberculosis H37Rv RNA-seq Data Reveals a Structured Interplay between Metabolism, Stress Response, and Infection

Mycobacterium tuberculosis is one of the most consequential human bacterial pathogens, posing a serious challenge to 21st century medicine. A key feature of its pathogenicity is its ability to adapt its transcriptional response to environmental stresses through its transcriptional regulatory network (TRN). While many studies have sought to characterize specific portions of the M. tuberculosis TRN, and some studies have performed system-level analysis, few have been able to provide a network-based model of the TRN that also provides the relative shifts in transcriptional regulator activity triggered by changing environments. Here, we compiled a compendium of nearly 650 publicly available, high quality M. tuberculosis RNA-sequencing data sets and applied an unsupervised machine learning method to obtain a quantitative, top-down TRN. It consists of 80 independently modulated gene sets known as “iModulons,” 41 of which correspond to known regulons. These iModulons explain 61% of the variance in the organism’s transcriptional response. We show that iModulons (i) reveal the function of poorly characterized regulons, (ii) describe the transcriptional shifts that occur during environmental changes such as shifting carbon sources, oxidative stress, and infection events, and (iii) identify intrinsic clusters of regulons that link several important metabolic systems, including lipid, cholesterol, and sulfur metabolism. This transcriptome-wide analysis of the M. tuberculosis TRN informs future research on effective ways to study and manipulate its transcriptional regulation and presents a knowledge-enhanced database of all published high-quality RNA-seq data for this organism to date.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Genomic analyses implicate noncoding de novo variants in congenital heart disease

A genetic etiology is identified for one-third of patients with congenital heart disease (CHD), with 8% of cases attributable to coding de novo variants (DNVs). To assess the contribution of noncoding DNVs to CHD, we compared genome sequences from 749 CHD probands and their parents with those from 1,611 unaffected trios. Neural network prediction of noncoding DNV transcriptional impact identified a burden of DNVs in individuals with CHD (n = 2,238 DNVs) compared to controls (n = 4,177; P = 8.7 × 10 -4 ). Independent analyses of enhancers showed an excess of DNVs in associated genes (27 genes versus 3.7 expected, P = 1 × 10 -5 ). We observed significant overlap between these transcription-based approaches (odds ratio (OR) = 2.5, 95% confidence interval (CI) 1.1-5.0, P = 5.4 × 10 -3 ). CHD DNVs altered transcription levels in 5 of 31 enhancers assayed. Finally, we observed a DNV burden in RNA-binding-protein regulatory sites (OR = 1.13, 95% CI 1.1-1.2, P = 8.8 × 10 -5 ). Overall, our findings demonstrate an enrichment of potentially disruptive regulatory noncoding DNVs in a fraction of CHD at least as high as that observed for damaging coding DNVs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

A systems genetics approach to deciphering the effect of dosage variation on leaf morphology in Populus

Abstract Gene copy number variation is frequent in plant genomes of various species, but the impact of such gene dosage variation on morphological traits is poorly understood. We used a large population of Populus carrying genomically characterized insertions and deletions across the genome to systematically assay the effect of gene dosage variation on a suite of leaf morphology traits. A systems genetics approach was used to integrate insertion and deletion locations, leaf morphology phenotypes, gene expression, and transcriptional network data, to provide an overview of how gene dosage influences morphology. Dosage-sensitive genomic regions were identified that influenced individual or pleiotropic morphological traits. We also identified cis-expression quantitative trait loci (QTL) within these dosage QTL regions, a subset of which modulated trans-expression QTL as well. Integration of data types within a gene co-expression framework identified co-expressed gene modules that are dosage sensitive, enriched for dosage expression QTL, and associated with morphological traits. Functional description of these modules linked dosage-sensitive morphological variation to specific cellular processes, as well as candidate regulatory genes. Together, these results show that gene dosage variation can influence morphological variation through complex changes in gene expression, and suggest that frequently occurring gene dosage variation has the potential to likewise influence quantitative traits in nature.

Biochemistry & Molecular Biology↗