Search NASA⌕ Search

SEARCH · Search NASA

Results for “gene function”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Controlling CRISPR-Cas9 with ligand-activated and ligand-deactivated sgRNAs

The CRISPR-Cas9 system provides the ability to edit, repress, activate, or mark any gene (or DNA element) by pairing of a programmable single guide RNA (sgRNA) with a complementary sequence on the DNA target. Here we present a new method for small-molecule control of CRISPR-Cas9 function through insertion of RNA aptamers into the sgRNA. We show that CRISPR-Cas9-based gene repression (CRISPRi) can be either activated or deactivated in a dose-dependent fashion over a >10-fold dynamic range in response to two different small-molecule ligands. Since our system acts directly on each target-specific sgRNA, it enables new applications that require differential and opposing temporal control of multiple genes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A riboswitch separated from its ribosome-binding site still regulates translation

Abstract Riboswitches regulate downstream gene expression by binding cellular metabolites. Regulation of translation initiation by riboswitches is posited to occur by metabolite-mediated sequestration of the Shine-Dalgarno sequence (SDS), causing bypass by the ribosome. Recently, we solved a co-crystal structure of a prequeuosine1-sensing riboswitch from Carnobacterium antarcticum that binds two metabolites in a single pocket. The structure revealed that the second nucleotide within the gene-regulatory SDS, G34, engages in a crystal contact, obscuring the molecular basis of gene regulation. Here, we report a co-crystal structure wherein C10 pairs with G34. However, molecular dynamics simulations reveal quick dissolution of the pair, which fails to reform. Functional and chemical probing assays inside live bacterial cells corroborate the dispensability of the C10–G34 pair in gene regulation, leading to the hypothesis that the compact pseudoknot fold is sufficient for translation attenuation. Remarkably, the C. antarcticum aptamer retained significant gene-regulatory activity when uncoupled from the SDS using unstructured spacers up to 10 nucleotides away from the riboswitch—akin to steric-blocking employed by sRNAs. Accordingly, our work reveals that the RNA fold regulates translation without SDS sequestration, expanding known riboswitch-mediated gene-regulatory mechanisms. The results infer that riboswitches exist wherein the SDS is not embedded inside a stable fold.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Maize similar to RCD1 gene induced by salt enhances Arabidopsis thaliana abiotic stress resistance

Highlights: • We cloned a novel SRO class gene, ZmSRO1b, induced by salt stress from maize. • Promoter activities of are ZmSRO1b induced by salt, cadmium and oxidative stress in transgenic Arabidopsis. • Overexpression of ZmSRO1b transgenic Arabidopsis plants enhanced tolerance of abiotic stress. Plant SRO (SIMILAR TO RCD-ONE) proteins play important roles in regulating oxidation and metal ion metabolism. Numbers of SRO proteins have been functional identified in Arabidopsis and rice, but little is known in maize. In this study, we identified a salt induced SRO gene, ZmSRO1b, from maize and analyzed its characteristics. ZmSRO1b expressed mainly in leaf tissues. The ZmSRO1b is encoded by 595 amino acid residues and shared conserved protein models with AtRCD1 and AtSRO1 from Arabidopsis. Promoter-elements analysis showed ZmSRO1b promoter harbored salt and metal stress responsive elements, DRE, GT-like and MRE. Further promoter inductive analysis by GUS staining and quantification confirmed that ZmSRO1b promoter was induced by salt and cadmium (Cd). Methylviologen (MV) simulated oxidative stress showed ZmSRO1b promoter was also induced by MV. Overexpression of ZmSRO1b in Arabidopsis plants showed increased resistance to salt, Cd and oxidative stress. Our results for the first time experimentally validate the function of ZmSRO1b and contribute to the better understanding of SRO genes across different plant species.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

EMT-related long non-coding RNA in hepatocellular carcinoma: A study with TCGA database

Highlights: • Integrated analysis identifies EMT-associated lncRNA in HCC. • Candidate lncRNA-associated genes are enriched in EMT-linked pathways. • LncRNA WDFY3-AS3, MIAT, and MEG3 contribute to the EMT phenotype of HCC cells. Accumulating evidence suggest that dysregulated expression of long non-coding RNA (lncRNA) plays a critical role in human tumorigenesis. However, little is known about the lncRNA implicated in the epithelial-to-mesenchymal transition (EMT) process. In this study, we performed data mining in The Cancer Genome Atlas (TCGA) hepatocellular carcinoma (HCC) data set and identified the a spectrum of differentially expressed lncRNAs implicated the EMT process of HCC, and functionally validated their roles in LM3 cells. Especially, lncRNA WDFY3-AS2-, LINC00472-, MIAT-, and MEG3-associated genes were significantly enriched in EMT-linked pathways. Loss-of-function study showed that genetic silencing of WDFY3-AS3, MIAT, and MEG3, but not LINC00472, resulted in reduced N-cadherin expression, cell migration, and cell invasion. Collectively, our results identify several lncRNAs that regulate the EMT process of HCC, which provides critical information for HCC tumorigenesis and potential therapeutic targets.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Identification of xylan arabinosyl 2- O -xylosyltransferases catalyzing the addition of 2- O -xylosyl residue onto arabinosyl side chains of xylan in grass species

Grass xylan, the major hemicellulose in both primary and secondary cell walls, is heavily decorated with α-1,3-linked arabinofuranosyl (Araf) residues that may be further substituted at O-2 with xylosyl (Xyl) or Araf residues. Although xylan 3-O-arabinosyltransferases (XATs) catalyzing 3-O-Araf addition onto xylan have been characterized, glycosyltransferases responsible for the transfer of 2-O-Xyl or 2-O-Araf onto 3-O-Araf residues of xylan to produce the Xyl-Araf and Araf-Araf disaccharide side chains remain to be identified. In this report, we showed that a rice GT61 member, named OsXAXT1 (xylan arabinosyl 2-O-xylosyltransferase 1) herein, was able to mediate the addition of Xyl-Araf disaccharide side chains onto xylan when heterologously co-expressed with OsXAT2 in the Arabidopsis gux1/2/3 (glucuronic acid substitution of xylan 1/2/3) triple mutant that lacks any glycosyl substitutions. Recombinant OsXAXT1 protein expressed in human embryonic kidney 293 cells exhibited a xylosyltransferase activity catalyzing the addition of Xyl from UDP-Xyl onto arabinosylated xylooligomers. Further, consistent with its function as a xylan arabinosyl 2-O-xylosyltransferase, CRISPR-Cas9-mediated mutations of the OsXAXT1 gene in transgenic rice plants resulted in a reduction in the level of Xyl-Araf disaccharide side chains in xylan. Furthermore, we revealed that XAXT1 close homologs from several other grass species, including switchgrass, maize, and Brachypodium, possessed the same functions as OsXAXT1, indicating functional conservation of XAXTs in grass species. Together, our findings establish that grass XAXTs are xylosyltransferases catalyzing Xyl transfer onto O-2 of Araf residues of xylan to form the Xyl-Araf disaccharide side chains, which furthers our understanding of genes involved in xylan biosynthesis.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2001

Differential expression of decorin and biglycan genes during mouse tooth development

Small leucine-rich proteoglycans (SLRPs) have a number of biological functions and some of them are thought to regulate collagen mineralizaton in bone and tooth. We have previously identified and immunolocalized two members of the SLRPs family, decorin and biglycan, in bovine tooth/periodontium. To investigate their potential roles in tooth development, we examined the mRNA expression patterns of decorin, biglycan and type I collagen in newborn (day 19) mice tooth germs by in situ hybridization. At this developmental stage, the first maxillary and mandibular molars include stages before and after secretion of the predentin matrix, respectively. The expression of decorin mRNA coincided with that of type I collagen mRNA and was mostly observed in secretory odontoblasts, while the biglycan mRNA was expressed throughout the tooth germ, including pre-secretory odontoblasts/ameloblasts, dental papilla and stellate reticulum. However, its signal in secretory odontoblasts was not as evident as that of decorin. In mandibular incisors, where a significant amount of predentin matrix and a small amount of enamel matrix were already secreted, a similar differential expression pattern was observed. In secretory ameloblasts the biglycan mRNA expression was apparent, while that of decorin was not. These differential expression patterns suggest the distinct roles of biglycan and decorin in the process of tooth development.

Non-NASA Center↗
DOE OSTI2024

The influence of microbial colonization on inflammatory versus pro-healing trajectories in combat extremity wounds

A combination of improved body armor, medical transportation, and treatment has led to the increased survival of warfighters from combat extremity injuries predominantly caused by blasts in modern conflicts. Despite advances, a high rate of complications such as wound infections, wound failure, amputations, and a decreased quality of life exist. To study the molecular underpinnings of wound failure, wound tissue biopsies from combat extremity injuries had RNA extracted and sequenced. Wounds were classified by colonization (colonized vs. non-colonized) and outcome (healed vs. failed) status. Differences in gene expression were investigated between timepoints at a gene level, and longitudinally by multi-gene networks, inferred proportions of immune cells, and expression of healing-related functions. Differences between wound outcomes in colonized wounds were more apparent than in non-colonized wounds. Colonized/healed wounds appeared able to mount an adaptive immune response to infection and progress beyond the inflammatory stage of healing, while colonized/failed wounds did not. Although, both colonized and non-colonized failed wounds showed increasing inferred immune and inflammatory programs, non-colonized/failed wounds progressed beyond the inflammatory stage, suggesting different mechanisms of failure dependent on colonization status. Overall, these data reveal gene expression profile differences in healing wounds that may be utilized to improve clinical treatment paradigms.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Structural and Functional Studies of the RBPJ-SHARP Complex Reveal a Conserved Corepressor Binding Site

Notch is a conserved signaling pathway that is essential for metazoan development and homeostasis; dysregulated signaling underlies the pathophysiology of numerous human diseases. Receptor-ligand interactions result in gene expression changes, which are regulated by the transcription factor RBPJ. RBPJ forms a complex with the intracellular domain of the Notch receptor and the coactivator Mastermind to activate transcription, but it can also function as a repressor by interacting with corepressor proteins. Here, we determine the structure of RBPJ bound to the corepressor SHARP and DNA, revealing its mode of binding to RBPJ. We tested structure-based mutants in biophysical and biochemical-cellular assays to characterize the role of RBPJ as a repressor, clearly demonstrating that RBPJ mutants deficient for SHARP binding are incapable of repressing transcription of genes responsive to Notch signaling in cells. Altogether, our structure-function studies provide significant insights into the repressor function of RBPJ.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Codon usage revisited: Lack of correlation between codon usage and the number of tRNA genes in enterobacteria

Highlights: • Average number of tRNA genes correlate with codon usage in highly expressed genes. • Changes in the number of tRNA genes do not correlate to codon usage. • Altered tRNA gene copy numbers do not correlate to changes of codon translation. • Our data indicates that tRNA gene copies adapt to codon usage, but not the contrary. It is widely believed that if a high number of genes are found for any tRNA in a rapidly replicating bacteria, then the cytoplasmic levels of that tRNA will be high and an open reading frame containing a higher frequency of the complementary codon will be translated faster. This idea is based on correlations between the number of tRNA genes, tRNA concentration and the frequency of codon usage observed in a limited number of strains as well as from the fact that artificially changing the number of tRNA genes alters translation efficiency and consequently the amount of properly folded protein synthesized. tRNA gene number may greatly vary in a genome due to duplications, deletions and lateral transfer which in turn would alter the levels and functionality of many proteins. Such changes are potentially deleterious for fitness and as a result it is expected that changes in tRNA gene numbers should be accompanied by a modification of the frequency of codon usage. In contrast to this model, when comparing the number of tRNA genes and the frequency of codon usage of several Salmonella enterica and Escherichia coli strains we found that changes in the number of tRNA genes are not correlated to changes in codon usage. Furthermore, these changes are not correlated with a change in the efficiency of codon translation. These results suggest that once a genome gains or loses tRNA genes, it responds by modulating the concentrations of tRNAs rather than modifying its frequency of codon usage.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

MicroPheno: predicting environments and host phenotypes from 16S rRNA gene sequencing using a k-mer based representation of shallow sub-samples

Microbial communities play important roles in the function and maintenance of various biosystems, ranging from the human body to the environment. A major challenge in microbiome research is the classification of microbial communities of different environments or host phenotypes. The most common and cost-effective approach for such studies to date is 16S rRNA gene sequencing. Recent falls in sequencing costs have increased the demand for simple, efficient and accurate methods for rapid detection or diagnosis with proved applications in medicine, agriculture and forensic science. We describe a reference- and alignment-free approach for predicting environments and host phenotypes from 16S rRNA gene sequencing based on k-mer representations that benefits from a bootstrapping framework for investigating the sufficiency of shallow sub-samples. Deep learning methods as well as classical approaches were explored for predicting environments and host phenotypes. A k-mer distribution of shallow sub-samples outperformed Operational Taxonomic Unit (OTU) features in the tasks of body-site identification and Crohn’s disease prediction. Aside from being more accurate, using k-mer features in shallow sub-samples allows (i) skipping computationally costly sequence alignments required in OTU-picking and (ii) provided a proof of concept for the sufficiency of shallow and short-length 16S rRNA sequencing for phenotype prediction. In addition, k-mer features predicted representative 16S rRNA gene sequences of 18 ecological environments, and 5 organismal environments with high macro-F1 scores of 0.88 and 0.87. For large datasets, deep learning outperformed classical methods such as Random Forest and Support Vector Machine.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Data for A Role for Differential Rubisco Activase Isoform Expression in C4 Bioenergy Grasses at High Temperature

Rubisco activase (Rca) facilitates the release of sugar-phosphate inhibitors at Rubisco catalytic sites during CO2 fixation. Most plant species express two Rca isoforms, the larger Rca-α and the shorter Rca-β, either by alternative splicing from a single gene or expression from separate genes. The mechanism of Rubisco activation by Rca isoforms has been intensively studied in C3 plants. However, the functional role of Rca in C4 plants where Rubisco and Rca are located in a much higher [CO2] compartment is less clear. In this study, we selected four C4 bioenergy grasses and the model C4 grass setaria ( Setaria viridis ) to investigate the role of Rca in C4 photosynthesis. All five C4 grass species contained two Rca genes, one encoding Rca-α and the other Rca-β, which were positioned closely together in the genomes. A variety of abiotic stress-related motifs were identified in the Rca-α promoter of each grass, and while the Rca-β gene was constantly highly expressed at ambient temperature, Rca-α isoforms were expressed only at high temperature but never surpassed 30% of Rca-β content. The pattern of Rca-α induction on transition to high temperature and reduction on return to ambient temperature was the same in all five C4 grasses. In sorghum ( Sorghum bicolor ), sugarcane ( Saccharum officinarum ), and setaria, the induction rate of Rca-α was similar to the recovery rate of photosynthesis and Rubisco activation at high temperature. This association between Rca-α isoform expression and maintenance of Rubisco activation at high temperature suggests that Rca-α has a functional thermo-protective role in carbon fixation in C4 grasses by sustaining Rubisco activation at high temperature.

Genomics↗
DOE OSTI2020

Fatty Acid and Alcohol Metabolism in Pseudomonas putida: Functional Analysis Using Random Barcode Transposon Sequencing.

With its ability to catabolize a wide variety of carbon sources and a growing engineering toolkit, Pseudomonas putida KT2440 is emerging as an important chassis organism for metabolic engineering. Despite advances in our understanding of the organism, many gaps remain in our knowledge of the genetic basis of its metabolic capabilities. The gaps are particularly noticeable in our understanding of both fatty acid and alcohol catabolism, where many paralogs putatively coding for similar enzymes coexist, making biochemical assignment via sequence homology difficult. To rapidly assign function to the enzymes responsible for these metabolisms, we leveraged random barcode transposon sequencing (RB-Tn-Seq). Global fitness analyses of transposon libraries grown on 13 fatty acids and 10 alcohols produced strong phenotypes for hundreds of genes. Fitness data from mutant pools grown on fatty acids of varying chain lengths indicated specific enzyme substrate preferences and enabled us to hypothesize that DUF1302/DUF1329 family proteins potentially function as esterases. From the data, we also postulate catabolic routes for the two biogasoline molecules isoprenol and isopentanol, which are catabolized via leucine metabolism after initial oxidation and activation with coenzyme A (CoA). Because fatty acids and alcohols may serve as both feedstocks and final products of metabolic-engineering efforts, the fitness data presented here will help guide future genomic modifications toward higher titers, rates, and yields. IMPORTANCE To engineer novel metabolic pathways into P. putida , a comprehensive understanding of the genetic basis of its versatile metabolism is essential. Here, we provide functional evidence for the putative roles of hundreds of genes involved in the fatty acid and alcohol metabolism of the bacterium. These data provide a framework facilitating precise genetic changes to prevent product degradation and to channel the flux of specific pathway intermediates as desired.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

A study of mechanistic mapping of novel SNPs to male breast cancer

Alterations in BRCA2, PALB2, CHEK2, and p53 genes have been identified for their association with male breast cancer in various studies. The incidence of male breast cancer in India is consistent with its global rate. The present study was carried out with an aim to evaluate the genetic alterations in male breast cancer patients from Malwa region of Punjab, India. Four male breast cancer patients belonging to different families were recruited from Guru Gobind Singh Medical College and Hospital, Faridkot, India. A total of 51 genes reported with implications in the pathogenesis of breast cancer were screened using next generation sequencing. Germline variations were found in BRCA1, BRCA2, PMS2, p53, and PALB2 genes, previously reported to be associated with MBC as well as FBC. In addition to these, 13 novel missense alterations were detected in eight genes including STK11, FZR1, PALB2, BRCA2, NF2, BAP1, BARD1, and CHEK2. Impact of these missense alterations on structure and function of protein was also analyzed through molecular dynamics simulation. Structural analysis of these single nucleotide polymorphisms (SNPs) revealed significant impact on the encoded protein functioning.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

Zinc chloride affects chondrogenesis via VEGF signaling

Highlights: • ZnCl2 increases expression of various chondrogenic signaling intermediates. • VEGF/VEGFR play a key role in modulating levels of chondrogenesis. • Insulin mimetics are an appealing option for improvement of fracture healing. • Axitinib + ZnCl2 treatment results in decreased transcription of chondrogenic genes. Insulin mimetics, including zinc containing compounds, have previously been shown to influence chondrogenesis as it relates to healing of fractures in various preclinical models. However, the mechanism by which these compounds drive chondrogenic differentiation is yet undefined. Here, via next-generation sequencing (NGS) and in vitro functional validation, we show that Zinc Chloride (ZnCl{sub 2}) induces expression of both chondrogenic genes (Sox9, Runx1, collagen) as well as genes associated with VEGF-mediated signal transduction, including VEGF receptors 1 and 2 and their ligands; VEGF-A and VEGF-B. Noticeably, although insulin was able to also induce expression of these pro-angiogenic and pro-chondrogenic genes, the impact of insulin on expression of VEGF receptor and ligand genes was marginal when compared to that of ZnCl{sub 2.} Furthermore, while the VEGFR antagonist, Axitinib, was able to attenuate the pro-chondrogenic effects of both insulin and ZnCl{sub 2}; a reduction in gene and protein expression was most profoundly observed when the antagonist was applied to cells treated with ZnCl{sub 2.} Taken together, these data suggest an important role for the VEGF-mediated signal transduction pathways in the positive effects observed when applying zinc-based compounds as adjuvants for chondrogenesis-mediated fracture healing. In this regard, further mechanistic evaluation of ZnCl{sub 2} and other zinc-containing insulin mimetics may support rational design of therapies targeted for disease indications associated with impaired fracture healing.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Transporter annotations are holding up progress in metabolic modeling

Mechanistic, constraint-based models of microbial isolates or communities are a staple in the metabolic analysis toolbox, but predictions about microbe-microbe and microbe-environment interactions are only as good as the accuracy of transporter annotations. A number of hurdles stand in the way of comprehensive functional assignments for membrane transporters. These include general or non-specific substrate assignments, ambiguity in the localization, directionality and reversibility of a transporter, and the many-to-many mapping of substrates, transporters and genes. In this perspective, we summarize progress in both experimental and computational approaches used to determine the function of transporters and consider paths forward that integrate both. Investment in accurate, high-throughput functional characterization is needed to train the next-generation of predictive tools toward genome-scale metabolic network reconstructions that better predict phenotypes and interactions. More reliable predictions in this domain will benefit fields ranging from personalized medicine to metabolic engineering to microbial ecology.

Casey, John↗
DOE OSTI2019

cindr, the Drosophila Homolog of the CD2AP Alzheimer's Disease Risk Gene, Is Required for Synaptic Transmission and Proteostasis

The Alzheimer’s disease (AD) susceptibility gene, CD2-associated protein (CD2AP), encodes an actin binding adaptor protein, but its function in the nervous system is largely unknown. Loss of the Drosophila ortholog cindr enhances neurotoxicity of human Tau, which forms neurofibrillary tangle pathology in AD. We show that Cindr is expressed in neurons and present at synaptic terminals. cindr mutants show impairments in synapse maturation and both synaptic vesicle recycling and release. Cindr associates and genetically interacts with 14-3-3?, regulates the ubiquitin-proteasome system, and affects turnover of Synapsin and the plasma membrane calcium ATPase (PMCA). Loss of cindr elevates PMCA levels and reduces cytosolic calcium. Studies of Cd2ap null mice support a conserved role in synaptic proteostasis, and CD2AP protein levels are inversely related to Synapsin abundance in human postmortem brains. Our results reveal CD2AP neuronal requirements with relevance to AD susceptibility, including for proteostasis, calcium handling, and synaptic structure and function.

Alzheimer's disease, risk gene, cindr, Drosophila,↗
DOE OSTI2023

Mycoparasites, Gut Dwellers, and Saprotrophs: Phylogenomic Reconstructions and Comparative Analyses of Kickxellomycotina Fungi

Improved sequencing technologies have profoundly altered global views of fungal diversity and evolution. High-throughput sequencing methods are critical for studying fungi due to the cryptic, symbiotic nature of many species, particularly those that are difficult to culture. However, the low coverage genome sequencing (LCGS) approach to phylogenomic inference has not been widely applied to fungi. Here we analyzed 171 Kickxellomycotina fungi using LCGS methods to obtain hundreds of marker genes for robust phylogenomic reconstruction. Additionally, we mined our LCGS data for a set of nine rDNA and protein coding genes to enable analyses across species for which no LCGS data were obtained. The main goals of this study were to: 1) evaluate the quality and utility of LCGS data for both phylogenetic reconstruction and functional annotation, 2) test relationships among clades of Kickxellomycotina, and 3) perform comparative functional analyses between clades to gain insight into putative trophic modes. In opposition to previous studies, our nine-gene analyses support two clades of arthropod gut dwelling species and suggest a possible single evolutionary event leading to this symbiotic lifestyle. Furthermore, we resolve the mycoparasitic Dimargaritales as the earliest diverging clade in the subphylum and find four major clades of Coemansia species. Finally, functional analyses illustrate clear variation in predicted carbohydrate active enzymes and secondary metabolites (SM) based on ecology, that is biotroph versus saprotroph. Saprotrophic Kickxellales broadly lack many known pectinase families compared with saprotrophic Mucoromycota and are depauperate for SM but have similar numbers of predicted chitinases as mycoparasitic.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Integrated network modeling approach defines key metabolic responses of soil microbiomes to perturbations

The soil environment is constantly changing due to shifts in soil moisture, nutrient availability and other conditions. To contend with these changes, soil microorganisms have evolved a variety of ways to adapt to environmental perturbations, including regulation of gene expression. However, it is challenging to untangle the complex phenotypic response of the soil to environmental change, partly due to the absence of predictive modeling frameworks that can mechanistically link molecular-level changes in soil microorganisms to a community’s functional phenotypes (or metaphenome). Towards filling this gap, we performed a combined analysis of metabolic and gene co-expression networks to explore how the soil microbiome responded to changes in soil moisture and nutrient conditions and to determine which genes were expressed under a given condition. Our integrated modeling approach revealed previously unknown, but critically important aspects of the soil microbiomes’ response to environmental perturbations. Incorporation of metabolomic and transcriptomic data into metabolic reaction networks identified condition-specific signature genes that are uniquely associated with dry, wet, and glycine-amended conditions. A subsequent gene co-expression network analysis revealed that drought-associated genes occupied more central positions in a network model of the soil community, compared to the genes associated with wet, and glycine-amended conditions. These results indicate the occurrence of system-wide metabolic coordination when soil microbiomes cope with moisture or nutrient perturbations. Importantly, the approach that we demonstrate here to analyze large-scale multi-omics data from a natural soil environment is applicable to other microbiome systems for which multi-omics data are available.

54 ENVIRONMENTAL SCIENCES↗