Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gene Expression Profiling”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Transcriptome and secretome profiling of sensory neurons reveals sex differences in pathways relevant to insulin sensing and insulin secretion

Sensory neurons in the dorsal root ganglia (DRG) convey somatosensory and metabolic cues to the central nervous system and release substances from stimulated terminal endings in peripheral organs. Sex-biased variations driven by the sex chromosome complement (XX and XY) have been implicated in the sensory–islet crosstalk. However, the molecular underpinnings of these male–female differences are not known. Here, we aim to characterize the molecular repertoire and the secretome profile of the lower thoracic spinal sensory neurons and to identify molecules with sex-biased insulin sensing- and/or insulin secretion-modulating activity that are encoded independently of circulating gonadal sex hormones. We used transcriptomics and proteomics to uncover differentially expressed genes and secreted molecules in lower thoracic T5-12 DRG sensory neurons derived from sexually immature 3-week-old male and female C57BL/6J mice. Comparative transcriptome and proteome analyses revealed differential gene expression and protein secretion in DRG neurons in males and females. Furthermore, the transcriptome analysis identified, among others, higher insulin signaling/sensing capabilities in female DRG neurons; secretome screening uncovered several sex-specific candidate molecules with potential regulatory functions in pancreatic β cells. Together, these data suggest a putative role of sensory interoception of insulin in the DRG–islet crosstalk with implications in sensory feedback loops in the regulation of β-cell activity in a sex-biased manner. Finally, we provide a valuable resource of molecular and secretory targets that can be leveraged for understanding insulin interoception and insulin secretion and inform the development of novel studies/approaches to fathom the role of the sensory–islet axis in the regulation of energy balance in males and females.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Islet β cell: An endocrine cell secreting miRNAs

Highlights: • Islets can secret not only insulin but also miRNAs. • Islets secret miRNAs actively and selectively and may have function on target tissues. • Most of the secret miRNAs by islets were packaged by exosomes. The pancreas is an important endocrine organ. Pancreatic beta (β) cells secrete insulin to regulate the metabolism of glucose and maintain the stability of blood glucose. MicroRNAs (miRNAs), small (20-25 nt) non-coding RNA molecules, function in target gene silencing through post-transcriptional patterns, and multiple miRNAs regulate insulin secretion, insulin signaling pathways in pancreatic β cells and islet activity. In recent years, it has been found that many types of cells can secrete regulatory miRNA, but whether islet cells can secrete miRNA has not been clarified yet. In this study, we detected miRNA expressions in cell culture medium of primary mouse islet cells and islet cell line MIN6. And we found that islet cells can selectively and actively secrete miRNAs, mainly through the way of exosome package, and that miRNA secretion profiling patterns vary according to different insulin secretion stimulating conditions (high Glucose, high Kcl, high Arginine and high Free fat acid (FFA)). Additionally, we found that these miRNAs secreted by islet cells can be transported and have gene-regulating functions on recipient tissue cells, such as liver and muscles. In conclusion, we find another secretory function of the islet β cells: not only insulin, but also miRNAs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Proteomic analysis of hair shafts from monozygotic twins: Expression profiles and genetically variant peptides

Forensic association of hair shaft evidence with individuals is currently assessed by comparing mitochondrial DNA haplotypes of reference and casework samples, primarily for exclusionary purposes. Present work tests and validates more recent proteomic approaches to extract quantitative transcriptional and genetic information from hair samples of monozygotic twin pairs, which would be predicted to partition away from unrelated individuals if the datasets contain identifying information. Protein expression profiles and polymorphic, genetically variant hair peptides were generated from ten pairs of monozygotic twins. Profiling using the protein tryptic digests revealed that samples from identical twins had typically an order of magnitude fewer protein expression differences than unrelated individuals. The data did not indicate that the degree of difference within twin pairs increased with age. In parallel, data from the digests were used to detect genetically variant peptides that result from common nonsynonymous single nucleotide polymorphisms in genes expressed in the hair follicle. Compilation of the variants permitted sorting of the samples by hierarchical clustering, permitting accurate matching of twin pairs. The results demonstrate that genetic differences are detectable by proteomic methods and provide a framework for developing quantitative statistical estimates of personal identification that increase the value of hair shaft evidence.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Collagen esterification enhances the function and survival of pancreatic β cells in 2D and 3D culture systems

Collagen, one of the most important components of the extracellular matrix (ECM), may play a role in the survival of pancreatic islet cells. In addition, chemical modifications that change the collagen charge profile to a net positive charge by esterification have been shown to increase the adhesion and proliferation of various cell types. The purpose of this study was to characterize and compare the effects of native collagen (NC) and esterified collagen (EC) on β cell function and survival. After isolation by the collagenase digestion technique, rat islets were cultured with NC and EC in 2 dimensional (2D) and 3 dimensional (3D) environments for a long-term duration in vitro. The cells were assessed for islet adhesion, morphology, viability, glucose-induced insulin secretion, and mRNA expression of glucose metabolism-related genes, and visualized by scanning electron microscopy (SEM). Islet cells attached tightly in the NC group, but islet cell viability was similar in both the NC and EC groups. Glucose-stimulated insulin secretion was higher in the EC group than in the NC group in both 2D and 3D culture. Furthermore, the mRNA expression levels of glucokinase in the EC group were higher than those in the NC group and were associated with glucose metabolism and insulin secretion. Finally, SEM observation confirmed that islets had more intact component cells on EC sponges than on NC sponges. These results indicate that modification of collagen may offer opportunities to improve function and viability of islet cells. - Highlights: • We changed the collagen charge profile to a net positive charge by esterification. • Islets cultured on esterified collagen improved survival in both 2D and 3D culture. • Islets cultured on esterified collagen enhanced glucose-stimulated insulin release. • High levels of glucokinase mRNA may be associated with increased insulin release.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2013

3D Normal Human Neural Progenitor Tissue-Like Assemblies: A Model of Persistent VZV Infection

Varicella-zoster virus (VZV) is a neurotropic human alphaherpesvirus that causes varicella upon primary infection, establishes latency in multiple ganglionic neurons, and can reactivate to cause zoster. Live attenuated VZV vaccines are available; however, they can also establish latent infections and reactivate. Studies of VZV latency have been limited to the analyses of human ganglia removed at autopsy, as the virus is strictly a human pathogen. Recently, terminally differentiated human neurons have received much attention as a means to study the interaction between VZV and human neurons; however, the short life-span of these cells in culture has limited their application. Herein, we describe the construction of a model of normal human neural progenitor cells (NHNP) in tissue-like assemblies (TLAs), which can be successfully maintained for at least 180 days in three-dimensional (3D) culture, and exhibit an expression profile similar to that of human trigeminal ganglia. Infection of NHNP TLAs with cell-free VZV resulted in a persistent infection that was maintained for three months, during which the virus genome remained stable. Immediate-early, early and late VZV genes were transcribed, and low-levels of infectious VZV were recurrently detected in the culture supernatant. Our data suggest that NHNP TLAs are an effective system to investigate long-term interactions of VZV with complex assemblies of human neuronal cells.

Goodwin, Thomas J.↗
DOE OSTI2024

Profiling expression strategies for a type III polyketide synthase in a lysate-based, cell-free system

Abstract Some of the most metabolically diverse species of bacteria (e.g., Actinobacteria) have higher GC content in their DNA, differ substantially in codon usage, and have distinct protein folding environments compared to tractable expression hosts like Escherichia coli . Consequentially, expressing biosynthetic gene clusters (BGCs) from these bacteria in E. coli often results in a myriad of unpredictable issues with regard to protein expression and folding, delaying the biochemical characterization of new natural products. Current strategies to achieve soluble, active expression of these enzymes in tractable hosts can be a lengthy trial-and-error process. Cell-free expression (CFE) has emerged as a valuable expression platform as a testbed for rapid prototyping expression parameters. Here, we use a type III polyketide synthase from Streptomyces griseus , RppA, which catalyzes the formation of the red pigment flaviolin, as a reporter to investigate BGC refactoring techniques. We applied a library of constructs with different combinations of promoters and rppA coding sequences to investigate the synergies between promoter and codon usage. Subsequently, we assess the utility of cell-free systems for prototyping these refactoring tactics prior to their implementation in cells. Overall, codon harmonization improves natural product synthesis more than traditional codon optimization across cell-free and cellular environments. More importantly, the choice of coding sequences and promoters impact protein expression synergistically, which should be considered for future efforts to use CFE for high-yield protein expression. The promoter strategy when applied to RppA was not completely correlated with that observed with GFP, indicating that different promoter strategies should be applied for different proteins. In vivo experiments suggest that there is correlation, but not complete alignment between expressing in cell free and in vivo. Refactoring promoters and/or coding sequences via CFE can be a valuable strategy to rapidly screen for catalytically functional production of enzymes from BCGs, which advances CFE as a tool for natural product research.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Chlamydomonas mutants lacking chloroplast TRIOSE PHOSPHATE TRANSPORTER3 are metabolically compromised and light sensitive

Abstract Modulation of photoassimilate export from the chloroplast is essential for controlling the distribution of fixed carbon in the cell and maintaining optimum photosynthetic rates. In this study, we identified chloroplast TRIOSE PHOSPHATE/PHOSPHATE TRANSLOCATOR 2 (CreTPT2) and CreTPT3 in the green alga Chlamydomonas (Chlamydomonas reinhardtii), which exhibit similar substrate specificities but whose encoding genes are differentially expressed over the diurnal cycle. We focused mostly on CreTPT3 because of its high level of expression and the severe phenotype exhibited by tpt3 relative to tpt2 mutants. Null mutants for CreTPT3 had a pleiotropic phenotype that affected growth, photosynthetic activities, metabolite profiles, carbon partitioning, and organelle-specific accumulation of H2O2. These analyses demonstrated that CreTPT3 is a dominant conduit on the chloroplast envelope for the transport of photoassimilates. In addition, CreTPT3 can serve as a safety valve that moves excess reductant out of the chloroplast and appears to be essential for preventing cells from experiencing oxidative stress and accumulating reactive oxygen species, even under low/moderate light intensities. Finally, our studies indicate subfunctionalization of the TRIOSE PHOSPHATE/PHOSPHATE TRANSLOCATOR (CreTPT) transporters and suggest that there are differences in managing the export of photoassimilates from the chloroplasts of Chlamydomonas and vascular plants.

Huang, (黄伟超), Weichao (ORCID:0000000189403110)↗
DOE OSTI2017

Circular RNA profiling reveals that circCOL3A1-859267 regulate type I collagen expression in photoaged human dermal fibroblasts

Production of type I collagen declines is a main characteristic during photoaging, but the mechanism is still not fully understood. Circular RNAs (circRNAs) are a class of newly identified non-coding RNAs with regulatory potency by sequestering miRNAs like a sponge. It's more stable than linear RNAs, and would be a useful tool for regulation of gene expression. However, the role of circRNAs in collagen expression during photoaging is still unclear. Here we performed deep sequencing of RNA generated from UVA irradiated and no irradiated human dermal fibroblasts (HDFs) and identified 29 significantly differentially expressed circRNAs (fold change ≥ 1.5, P < 0.05), 12 circRNAs were up-regulated and 17 circRNAs were down-regulated.3 most differentially expressed circRNAs were verified by qRT-PCR and the down-regulated circCOL3A1-859267 exhibited the most significantly altered in photoaged HDFs. Overexpression of circCOL3A1-859267 inhibited UVA-induced decrease of type I collagen expression and silencing of it reduced type I collagen intensity. Via a bioinformatic method, 44 miRNAs were predicted to binding with circCOL3A1-859267, 5 of them have been confirmed or predicted to interact with type I collagen. This study show that circCOL3A1-859267 regulate type I collagen expression in photoaged HDFs, suggesting it may be a novel target for interfering photoaging.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Transcriptome Profiling of Rust Resistance in Switchgrass Using RNA-Seq Analysis

Switchgrass rust caused by Puccinia emaculata is a major limiting factor for switchgrass (Panicum virgatum L.) production, especially in monoculture. Natural populations of switchgrass displayed diverse reactions to P. emaculata when evaluated in an Ardmore, OK, field. In order to identify the differentially expressed genes during the rust infection process and the mechanisms of switchgrass rust resistance, transcriptome analysis using RNA-Seq was conducted in two pseudo-F 1 parents ('PV281' and 'NFGA472'), and three moderately resistant and three susceptible progenies selected from a three-generation, four-founder switchgrass population (K5 x A4) x (AP13 x VS16). On average, 23.5 million reads per sample (leaf tissue was collected at 0, 24, and 60 h post-inoculation (hpi)) were obtained from paired-end (2 x 100 bp) sequencing on the Illumina HiSeq2000 platform. Furthermore, mapping of the RNA-Seq reads to the switchgrass reference genome (AP13 ver. 1.1 assembly) constructed a total of 84,209 transcripts from 98,007 gene loci among all of the samples. Further analysis revealed that host defense- related genes, including the nucleotide binding site-leucinerich repeat domain containing disease resistance gene analogs, play an important role in resistance to rust infection. Rust-induced gene (RIG) transcripts inherited across generations were identified. The rust-resistant gene transcripts can be a valuable resource for developing molecular markers for rust resistance. Finally we identified the rust-resistant genotypes and gene transcripts which can expedite rust-resistant cultivar development in switchgrass.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

TRESK potassium channel in human T lymphoblasts

Highlights: • TRESK (KCNK18) mRNA is present in different T lymphoblastic cell lines. • KCNK18 mRNA was not found in resting peripheral blood lymphocytes. • Clinical samples of T lymphoblastic leukemias and lymphomas were positive for TRESK. • TRESK in T lymphoblasts has dual localization, in plasma membrane and intracellular. -- Abstract: TRESK (TWIK-related spinal cord K{sup +}) channel, encoded by KCNK18 gene, belongs to the double-pore domain K{sup +} channel family and in normal conditions is expressed predominantly in the central nervous system. In our previous patch-clamp study on Jurkat T lymphoblasts we have characterized highly selective K{sup +} channel with pharmacological profile identical to TRESK. In the present work, the presence of KCNK18 mRNA was confirmed in T lymphoblastic cell lines (Jurkat, JCaM, H9) but not in resting peripheral blood lymphocytes of healthy donors. Positive immunostaining for TRESK was demonstrated in lymphoblastic cell lines, in germinal centers of non-tumoral lymph nodes, and in clinical samples of T acute lymphoblastic leukemias/lymphomas. Besides detection in the plasma membrane, intracellular TRESK localization was also revealed. Possible involvement of TRESK channel in lymphocyte proliferation and tumorigenesis is discussed.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Differential DNA methylation profile of key genes in malignant prostate epithelial cells transformed by inorganic arsenic or cadmium

Previous work shows altered methylation patterns in inorganic arsenic (iAs)- or cadmium (Cd)-transformed epithelial cells. Here, the methylation status near the transcriptional start site was assessed in the normal human prostate epithelial cell line (RWPE-1) that was malignantly transformed by 10 μM Cd for 11 weeks (CTPE) or 5 μM iAs for 29 weeks (CAsE-PE), at which time cells showed multiple markers of acquired cancer phenotype. Next generation sequencing of the transcriptome of CAsE-PE cells identified multiple dysregulated genes. Of the most highly dysregulated genes, five genes that can be relevant to the carcinogenic process (S100P, HYAL1, NTM, NES, ALDH1A1) were chosen for an in-depth analysis of the DNA methylation profile. DNA was isolated, bisulfite converted, and combined bisulfite restriction analysis was used to identify differentially methylated CpG sites, which was confirmed with bisulfite sequencing. Four of the five genes showed differential methylation in transformants relative to control cells that was inversely related to altered gene expression. Increased expression of HYAL1 (> 25-fold) and S100P (> 40-fold) in transformants was correlated with hypomethylation near the transcriptional start site. Decreased expression of NES (> 15-fold) and NTM (> 1000-fold) in transformants was correlated with hypermethylation near the transcriptional start site. ALDH1A1 expression was differentially expressed in transformed cells but was not differentially methylated relative to control. In conclusion, altered gene expression observed in Cd and iAs transformed cells may result from altered DNA methylation status. - Highlights: • Cd and iAs are known human carcinogens, yet neither appears directly mutagenic. • Prior data suggest epigenetic modification plays a role in Cd or iAs induced cancer. • Altered methylation of four misregulated genes was found in Cd or iAs transformants. • The resulting altered gene expression may be relevant to cellular transformation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Defining the Boundaries and Characterizing the Landscape of Genome Expression in Vascular Tissues of Populus using Shotgun Proteomics

Current state-of-the-art experimental and computational proteomic approaches were integrated to obtain a comprehensive protein profile of Populus vascular tissue. This featured: 1) a large sample set consisting of two genotypes grown under normal and tension stress conditions, 2) bioinformatics clustering to effectively handle gene duplication, and 3) an informatics approach to track and identify single amino acid polymorphisms (SAAPs). By applying a clustering algorithm to the Populus database, the number of protein entries decreased from 64,689 proteins to a total of 43,069 protein groups, thereby reducing 7,505 identified proteins to a total of 4,226 protein groups, in which 2,016 were singletons. This reduction implies that ~50% of the measured proteins were clustered into groups that shared extensive sequence homology. Using conservative search criteria, we were able to identify 1,354 peptides containing a SAAP and 201 peptides that become tryptic due to a K or R substitution. These newly identified peptides correspond to 502 proteins, including 97 proteins that were not previously identified. In total, the integration of deep proteome measurements on an extensive sample set with protein clustering and peptide sequence variants provided an unprecedented level of proteome characterization for Populus, allowing us to spatially resolve the vascular tissue proteome.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Characterization of growth and development of sorghum genotypes with differential susceptibility to Striga hermonthica

Two sorghum varieties, Shanqui Red (SQR) and SRN39, have distinct levels of susceptibility to the parasitic weed Striga hermonthica, which have been attributed to different strigolactone composition within their root exudates. Root exudates of the Striga-susceptible variety Shanqui Red (SQR) contain primarily 5-deoxystrigol, which has a high efficiency for inducing Striga germination. SRN39 roots primarily exude orobanchol, leading to reduced Striga germination and making this variety resistant to Striga. The structural diversity in exuded strigolactones is determined by a polymorphism in the LOW GERMINATION STIMULANT 1 (LGS1) locus. Yet, the genetic diversity between SQR and SRN39 is broad and has not been addressed in terms of growth and development. Here, we demonstrate additional differences between SQR and SRN39 by phenotypic and molecular characterization. A suite of genes related to metabolism was differentially expressed between SQR and SRN39. Increased levels of gibberellin precursors in SRN39 were accompanied by slower growth rate and developmental delay and we observed an overall increased SRN39 biomass. The slow-down in growth and differences in transcriptome profiles of SRN39 were strongly associated with plant age. Additionally, enhanced lateral root growth was observed in SRN39 and three additional genotypes exuding primarily orobanchol. In summary, we demonstrate that the differences between SQR and SRN39 reach further than the changes in strigolactone profile in the root exudate and translate into alterations in growth and development.

09 BIOMASS FUELS↗
DOE OSTI2021

Knockout of the lignin pathway gene BnF5H decreases the S/G lignin compositional ratio and improves Sclerotinia sclerotiorum resistance in Brassica napus

Ferulate-5-hydroxylase is a key enzyme involved in the conversion of the guaiacyl monolignol to the syringyl monolignol in angiosperms. The monolignol ratio has been proposed to affect biomass recalcitrance and the resistance to plant disease. Stem rot caused by the fungus Sclerotinia sclerotiorum in Brassica napus causes severe losses in its production. To date, there is no information about the effect of the lignin monomer ratio on the resistance to S. sclerotiorum in B. napus. Four dominantly expressed ferulate-5-hydroxylase genes were concertedly knocked out by CRISPR/Cas9 in B. napus, and three mutant lines were generated. The S/G lignin compositional ratio was decreased compared to that of the wild type based on the results of Mӓule staining and 2D-NMR profiling in KO-7. The resistance to S. sclerotiorum in stems and leaves increased for the three f5h mutant lines compared with WT. Furthermore, we found that the stem strength of f5h mutant lines was significantly increased. Overall, we demonstrate for the first time that decreasing the S/G ratio by knocking out of the F5H gene improves S. sclerotiorum resistance in B. napus and increases stem strength.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Double DAP-seq uncovered synergistic DNA binding of interacting bZIP transcription factors

Many eukaryotic transcription factors (TF) form homodimer or heterodimer complexes to regulate gene expression. Dimerization of BASIC LEUCINE ZIPPER (bZIP) TFs are critical for their functions, but the molecular mechanism underlying the DNA binding and functional specificity of homo- versus heterodimers remains elusive. To address this gap, we present the double DNA Affinity Purification-sequencing (dDAP-seq) technique that maps heterodimer binding sites on endogenous genomic DNA. Using dDAP-seq we profile twenty pairs of C/S1 bZIP heterodimers and S1 homodimers in Arabidopsis and show that heterodimerization significantly expands the DNA binding preferences of these TFs. Analysis of dDAP-seq binding sites reveals the function of bZIP9 in abscisic acid response and the role of bZIP53 heterodimer-specific binding in seed maturation. The C/S1 heterodimers show distinct preferences for the ACGT elements recognized by plant bZIPs and motifs resembling the yeast GCN4 cis -elements. This study demonstrates the potential of dDAP-seq in deciphering the DNA binding specificities of interacting TFs that are key for combinatorial gene regulation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Temporal change in chromatin accessibility predicts regulators of nodulation in Medicago truncatula

Symbiotic associations between bacteria and leguminous plants lead to the formation of root nodules that fix nitrogen needed for sustainable agricultural systems. Symbiosis triggers extensive genome and transcriptome remodeling in the plant, yet an integrated understanding of the extent of chromatin changes and transcriptional networks that functionally regulate gene expression associated with symbiosis remains poorly understood. In particular, analyses of early temporal events driving this symbiosis have only captured correlative relationships between regulators and targets at mRNA level. Here, we characterize changes in transcriptome and chromatin accessibility in the model legume Medicago truncatula, in response to rhizobial signals that trigger the formation of root nodules. We profiled the temporal chromatin accessibility (ATAC-seq) and transcriptome (RNA-seq) dynamics of M. truncatula roots treated with bacterial small molecules called lipo-chitooligosaccharides that trigger host symbiotic pathways of nodule development. Using a novel approach, dynamic regulatory module networks, we integrated ATAC-seq and RNA-seq time courses to predict cis-regulatory elements and transcription factors that most significantly contribute to transcriptomic changes associated with symbiosis. Regulators involved in auxin (IAA4-5, SHY2), ethylene (EIN3, ERF1), and abscisic acid (ABI5) hormone response, as well as histone and DNA methylation (IBM1), emerged among those most predictive of transcriptome dynamics. RNAi-based knockdown of EIN3 and ERF1 reduced nodule number in M. truncatula validating the role of these predicted regulators in symbiosis between legumes and rhizobia. Our transcriptomic and chromatin accessibility datasets provide a valuable resource to understand the gene regulatory programs controlling the early stages of the dynamic process of symbiosis. The regulators identified provide potential targets for future experimental validation, and the engineering of nodulation in species is unable to establish that symbiosis naturally.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Improved evidence-based genome-scale metabolic models for maize leaf, embryo, and endosperm

There is a growing demand for genome-scale metabolic reconstructions for plants, fueled by the need to understand the metabolic basis of crop yield and by progress in genome and transcriptome sequencing. Methods are also required to enable the interpretation of plant transcriptome data to study how cellular metabolic activity varies under different growth conditions or even within different organs, tissues, and developmental stages. Such methods depend extensively on the accuracy with which genes have been mapped to the biochemical reactions in the plant metabolic pathways. Errors in these mappings lead to metabolic reconstructions with an inflated number of reactions and possible generation of unreliable metabolic phenotype predictions. Here we introduce a new evidence-based genome-scale metabolic reconstruction of maize, with significant improvements in the quality of the gene-reaction associations included within our model. We also present a new approach for applying our model to predict active metabolic genes based on transcriptome data. This method includes a minimal set of reactions associated with low expression genes to enable activity of a maximum number of reactions associated with high expression genes. We apply this method to construct an organ-specific model for the maize leaf, and tissue specific models for maize embryo and endosperm cells. We validate our models using fluxomics data for the endosperm and embryo, demonstrating an improved capacity of our models to fit the available fluxomics data. All models are publicly available via the DOE Systems Biology Knowledgebase and PlantSEED, and our new method is generally applicable for analysis transcript profiles from any plant, paving the way for further in silico studies with a wide variety of plant genomes.

flux balance analysis↗