Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gene prediction”

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 127 records · Page 7

Single residue substitutions that change the gating properties of a mechanosensitive channel in Escherichia coli

MscL is a channel that opens a large pore in the Escherichia coli cytoplasmic membrane in response to mechanical stress. Previously, we highly enriched the MscL protein by using patch clamp as a functional assay and cloned the corresponding gene. The predicted protein contains a largely hydrophobic core spanning two-thirds of the molecule and a more hydrophilic carboxyl terminal tail. Because MscL had no homology to characterized proteins, it was impossible to predict functional regions of the protein by simple inspection. Here, by mutagenesis, we have searched for functionally important regions of this molecule. We show that a short deletion from the amino terminus (3 amino acids), and a larger deletion of 27 amino acids from the carboxyl terminus of this protein, had little if any effect in channel properties. We have thus narrowed the search of the core mechanosensitive mechanism to 106 residues of this 136-amino acid protein. In contrast, single residue substitutions of a lysine in the putative first transmembrane domain or a glutamine in the periplasmic loop caused pronounced shifts in the mechano-sensitivity curves and/or large changes in the kinetics of channel gating, suggesting that the conformational structure in these regions is critical for normal mechanosensitive channel gating.

Non-NASA Center↗

Regulators of early maize leaf development inferred from transcriptomes of laser capture microdissection (LCM)-isolated embryonic leaf cells

The superior photosynthetic efficiency of C 4 leaves over C 3 leaves is owing to their unique Kranz anatomy, in which the vein is surrounded by one layer of bundle sheath (BS) cells and one layer of mesophyll (M) cells. Kranz anatomy development starts from three contiguous ground meristem (GM) cells, but its regulators and underlying molecular mechanism are largely unknown. To identify the regulators, we obtained the transcriptomes of 11 maize embryonic leaf cell types from five stages of pre-Kranz cells starting from median GM cells and six stages of pre-M cells starting from undifferentiated cells. Principal component and clustering analyses of transcriptomic data revealed rapid pre-Kranz cell differentiation in the first two stages but slow differentiation in the last three stages, suggesting early Kranz cell fate determination. In contrast, pre-M cells exhibit a more prolonged transcriptional differentiation process. Differential gene expression and coexpression analyses identified gene coexpression modules, one of which included 3 auxin transporter and 18 transcription factor (TF) genes, including known regulators of Kranz anatomy and/or vascular development. In situ hybridization of 11 TF genes validated their expression in early Kranz development. We determined the binding motifs of 15 TFs, predicted TF target gene relationships among the 18 TF and 3 auxin transporter genes, and validated 67 predictions by electrophoresis mobility shift assay. From these data, we constructed a gene regulatory network for Kranz development. Our study sheds light on the regulation of early maize leaf development and provides candidate leaf development regulators for future study.

59 BASIC BIOLOGICAL SCIENCES↗

The relationship between gene traits and transcription in soil microbial communities varies by environmental stimulus

Codon and nucleotide frequencies are known to relate to the rate of gene transcription, yet how these traits shape transcriptional profiles of soil microbial communities remains unclear. Here we test the prediction that functional genes with high codon optimization and energetically lower cost nucleotides (i.e., nucleotides requiring less adenosine triphosphate (ATP) for synthesis) have higher transcriptional expression in a soil microbial community. In laboratory incubations, we subjected an agricultural soil to two separate short-term environmental changes: labile carbon (glucose) addition or a sudden 30-min increase in temperature from 20 °C to 60 °C. Using the total genomic codon frequencies to predict preferred codon usage for each taxon, we then estimated codon optimization for each transcript. On the community level, we found a higher average level of codon optimization after the addition of glucose. Synonymous nucleotide composition in the transcript pool also shifted towards energetically cheaper nucleotides, favoring uracil (U) over adenine (A) and cytosine (C) over guanine (G). Similarly, we found that encoded amino acid usage shifted towards energetically cheaper amino acids in response to labile carbon. In contrast, in communities responding to heat shock, there were no significant differences in the averaged gene traits of expressed transcripts. We used metagenome-assembled-genomes to further examine the ability of gene traits to predict transcriptional responses within and between taxa. We found that traits of individual genes could not reliably predict the level of transcription of a gene within or between taxa—highlighting the limits of this approach. However, we did find that when traits were averaged across several related genes, codon optimization was able to predict levels of transcription in metabolic pathways associated with growth and nutrient uptake in response to glucose. Similar relationships were not observed in response to heat, or for functions associated with stress—such as genes associated with sporulation or heat shock. These results demonstrate that gene traits, such as codon usage, nucleotide selection, and amino acid selection, relate to the transcriptional expression of genes in soil microbial communities and suggests that these relationships may be dependent on both gene function and the specific type of environmental stimuli.

Biological and medical sciences↗

Systems-Level Modeling for CRISPR-Based Metabolic Engineering

The CRISPR-Cas system has enabled the development of sophisticated, multigene metabolic engineering programs through the use of guide RNA-directed activation or repression of target genes. To optimize biosynthetic pathways in microbial systems, we need improved models to inform design and implementation of transcriptional programs. Recent progress has resulted in new modeling approaches for identifying gene targets and predicting the efficacy of guide RNA targeting. Genome-scale and flux balance models have successfully been applied to identify targets for improving biosynthetic production yields using combinatorial CRISPR-interference (CRISPRi) programs. Here, the advent of new approaches for tunable and dynamic CRISPR activation (CRISPRa) promises to further advance these engineering capabilities. Once appropriate targets are identified, guide RNA prediction models can lead to increased efficacy in gene targeting. Developing improved models and incorporating approaches from machine learning may be able to overcome current limitations and greatly expand the capabilities of CRISPR-Cas9 tools for metabolic engineering.

59 BASIC BIOLOGICAL SCIENCES↗

Data augmentation and multimodal learning for predicting drug response in patient-derived xenografts from gene expressions and histology images

Patient-derived xenografts (PDXs) are an appealing platform for preclinical drug studies. A primary challenge in modeling drug response prediction (DRP) with PDXs and neural networks (NNs) is the limited number of drug response samples. We investigate multimodal neural network (MM-Net) and data augmentation for DRP in PDXs. The MM-Net learns to predict response using drug descriptors, gene expressions (GE), and histology whole-slide images (WSIs). We explore whether combining WSIs with GE improves predictions as compared with models that use GE alone. We propose two data augmentation methods which allow us training multimodal and unimodal NNs without changing architectures with a single larger dataset: 1) combine single-drug and drug-pair treatments by homogenizing drug representations, and 2) augment drug-pairs which doubles the sample size of all drug-pair samples. Unimodal NNs which use GE are compared to assess the contribution of data augmentation. The NN that uses the original and the augmented drug-pair treatments as well as single-drug treatments outperforms NNs that ignore either the augmented drug-pairs or the single-drug treatments. In assessing the multimodal learning based on the MCC metric, MM-Net outperforms all the baselines. Our results show that data augmentation and integration of histology images with GE can improve prediction performance of drug response in PDXs.

60 APPLIED LIFE SCIENCES↗

Machine learning reveals genes impacting oxidative stress resistance across yeasts

Reactive oxygen species (ROS) are highly reactive molecules encountered by yeasts during routine metabolism and during interactions with other organisms, including host infection. Here, we characterized the variation in resistance to ROS across the ancient yeast subphylum Saccharomycotina and used machine learning (ML) to identify gene families whose sizes were predictive of ROS resistance.

AI↗

The Gene Fitness Atlas: A Roadmap for Predicting Evolution

We developed a novel, high-throughput microfluidic device design containing “interaction zones” where progeny cell lines compete against each other allowing for accurate analysis of bacterial cell fitness. The goal of the project was to use the device for two applications: 1) gene knockout screening and 2) antibiotic resistance screening. The microfluidic platform was fabricated using photolithography and soft lithography in polydimethylsiloxane (PDMS). E.coli Keio mutants and fluorescent wildtype parent were chosen for the study. Cells were grown overnight and their loading into the devices and seeding in mother machines was optimized. For mutant screening, the least fit mutant and wildtype parent were cultured individually and then added to the microfluidic device. The mother machines which were seeded with mutant and wildtype were imaged through time lapse microscopy and the growth of cells was observed. For antibiotic screening, wildtype E.coli cells which were grown overnight were added to the device and washed with media containing the antibiotic ampicillin. The growth pattern in presence and absence of ampicillin was observed through time lapse microscopy. It was observed that over a period of four hours, both the mutant and the wildtype divided in the mother machine and pushed daughter cells out into the interaction zone. In case of the antibiotic screening experiment, the fluorescent wildtype divided both in the absence and presence of sublethal concentration of ampicillin. This study is a proof of concept demonstration of high- throughput single cell analysis of cells using a novel microfluidics device.

59 BASIC BIOLOGICAL SCIENCES↗

An expanded role for the transcription factor WRINKLED1 in the biosynthesis of triacylglycerols during seed development

The transcription factor WRINKLED1 ( WRI1 ) is known as a master regulator of fatty acid synthesis in developing oilseeds of Arabidopsis thaliana and other species. WRI1 is known to directly stimulate the expression of many fatty acid biosynthetic enzymes and a few targets in the lower part of the glycolytic pathway. However, it remains unclear to what extent and how the conversion of sugars into fatty acid biosynthetic precursors is controlled by WRI 1. To shortlist possible gene targets for future in-planta experimental validation, here we present a strategy that combines phylogenetic foot printing of cis-regulatory elements with additional layers of evidence. Upstream regions of protein-encoding genes in A. thaliana were searched for the previously described DNA-binding consensus for WRI1, the ASML1/WRI1 (AW)-box. For about 900 genes, AW-box sites were found to be conserved across orthologous upstream regions in 11 related species of the crucifer family. For 145 select potential target genes identified this way, affinity of upstream AW-box sequences to WRI1 was assayed by Microscale Thermophoresis. This allowed definition of a refined WRI1 DNA-binding consensus. We find that known WRI1 gene targets are predictable with good confidence when upstream AW-sites are phylogenetically conserved, specifically binding WRI1 in the in vitro assay, positioned in proximity to the transcriptional start site, and if the gene is co-expressed with WRI1 during seed development. When targets predicted in this way are mapped to central metabolism, a conserved regulatory blueprint emerges that infers concerted control of contiguous pathway sections in glycolysis and fatty acid biosynthesis by WRI1. Several of the newly predicted targets are in the upper glycolysis pathway and the pentose phosphate pathway. Of these, plastidic isoforms of fructokinase ( FRK 3) and of phosphoglucose isomerase ( PGI 1) are particularly corroborated by previously reported seed phenotypes of respective null mutations.

59 BASIC BIOLOGICAL SCIENCES↗

Exploring Camelina sativa lipid metabolism regulation by combining gene co-expression and DNA affinity purification analyses

Camelina (Camelina sativa) is an annual oilseed plant that is gaining momentum as a biofuel cover crop. Understanding gene regulatory networks is essential to deciphering plant metabolic pathways, including lipid metabolism. Furthermore, we take advantage of a growing collection of gene expression datasets to predict transcription factors (TFs) associated with the control of Camelina lipid metabolism. We identified approximately 350 TFs highly co-expressed with lipid-related genes (LRGs). These TFs are highly represented in the MYB, AP2/ERF, bZIP, and bHLH families, including a significant number of homologs of well-known Arabidopsis lipid and seed developmental regulators. After prioritizing the top 22 TFs for further validation, we identified DNA-binding sites and predicted target genes for 16 out of the 22 TFs tested using DNA affinity purification followed by sequencing (DAP-seq). Enrichment analyses of targets supported the co-expression prediction for most TF candidates, and the comparison to Arabidopsis revealed some common themes, but also aspects unique to Camelina. Within the top potential lipid regulators, we identified CsaMYB1, CsaABI3AVP1-2, CsaHB1, CsaNAC2, CsaMYB3, and CsaNAC1 as likely involved in the control of seed fatty acid elongation and CsaABI3AVP1-2 and CsabZIP1 as potential regulators of the synthesis and degradation of triacylglycerols (TAGs), respectively. Altogether, the integration of co-expression data and DNA-binding assays permitted us to generate a high-confidence and short list of Camelina TFs involved in the control of lipid metabolism during seed development.

59 BASIC BIOLOGICAL SCIENCES↗

Spatiotemporal Dynamics of the Relative Abundance of Soil Nutrient‐Degrading Enzyme‐Encoding Genes Across Continental US Ecoregions

Understanding the spatiotemporal patterns in the relative abundance of soil extracellular enzyme‐encoding genes is critical for predicting microbial responses to environmental change and their potential role in nutrient cycling. Yet, integrating novel metagenomic observations with spatiotemporal environmental gradients to infer regional patterns and future trajectories has remained unclear. To address this gap, we applied a machine learning (ML) approach, integrating soil metagenomic data with environmental variables—soil properties, topography, vegetation, and climate—to predict the relative abundance of enzyme‐encoding genes for soil carbon (C), nitrogen (N), and phosphorus (P) across surface soils of the continental United States. We assessed potential responses under future emission scenarios (SSP2‐4.5 and SSP5‐8.5) by comparing a baseline (1985–2014) to a future period (2071–2100). The ML model explained 57%–63% of baseline variation. Precipitation was identified as the most influential factor for the relative abundance of C‐ and N‐degrading enzyme‐encoding genes, while slope length, representing horizontal distance that water can travel downslope, was the primary driver for P‐degrading enzyme‐encoding genes abundance. Projections revealed spatially heterogeneous shifts across continental US ecoregions: the relative abundance of C‐ and N‐degrading enzyme‐encoding genes decreased in wetter ecoregions and increased in drier ecoregions under future climate, while P‐degrading enzyme‐encoding genes abundance decreased significantly in semiarid and Mediterranean ecoregions. This study demonstrates the utility of metagenomic data for mapping soil genetic potential and predicting its regional response to environmental change, to inform ecosystem management strategies.

extracellular enzyme-encoding genes↗

Atlas of Transcription Factor Binding Sites from ENCODE DNase Hypersensitivity Data across 27 Tissue Types

Characterizing the tissue-specific binding sites of transcription factors (TFs) is essential to reconstruct gene regulatory networks and predict functions for non-coding genetic variation. DNase-seq footprinting enables the prediction of genome-wide binding sites for hundreds of TFs simultaneously. Despite the public availability of high-quality DNase-seq data from hundreds of samples, a comprehensive, up-to-date resource for the locations of genomic footprints is lacking. Here, we develop a scalable footprinting workflow using two state-of-the-art algorithms: Wellington and HINT. We apply our workflow to detect footprints in 192 ENCODE DNase-seq experiments and predict the genomic occupancy of 1,515 human TFs in 27 human tissues. We validate that these footprints overlap true-positive TF binding sites from ChIP-seq. We demonstrate that the locations, depth, and tissue specificity of footprints predict effects of genetic variants on gene expression and capture a substantial proportion of genetic risk for complex traits.

59 BASIC BIOLOGICAL SCIENCES↗

Disruption of the tagF Orthologue in the epa Locus Variable Region of Enterococcus faecalis Causes Cell Surface Changes and Suppresses an eep -Dependent Lysozyme Resistance Phenotype

The disease-producing capacity of the opportunistic pathogen Enterococcus faecalis is enhanced by the ability of the bacterium to evade killing by antimicrobial agents. Survival of E. faecalis in the presence of the human antimicrobial enzyme lysozyme is mediated in part by the site 2 metalloprotease Eep; however, a complete model of enterococcal lysozyme resistance has not been elucidated. To better understand the molecular basis for lysozyme resistance in E. faecalis, we analyzed Δeep suppressor mutants that acquire resistance to lysozyme through mutation of the gene OG1RF_11713, a predicted teichoic acid biosynthesis-encoding gene located within the variable region of the enterococcal polysaccharide antigen (epa) locus. Sequence comparisons revealed that OG1RF_11713 is most similar to the cytidine-5'-diphosphate (CDP)-glycerol:poly-(glycerolphosphate)glycerophosphotransferase TagF from Staphylococcus epidermidis. Inactivation of OG1RF_11713 in both the wild-type and Δeep genetic backgrounds was sufficient to increase the resistance of E. faecalis OG1RF to lysozyme. Minimal amounts of N-acetylgalactosamine were detectable in cell wall carbohydrate extracts of OG1RF_11713 deletion mutants, and this was associated with a reduction in negative cell surface charge. Targeted disruption of OG1RF_11713 was also associated with increased susceptibility to the antibiotic polymyxin B and membrane-targeting detergents and decreased susceptibility to the lantibiotic nisin. Furthermore, this work implicates OG1RF_11713 as a major determinant of cell envelope integrity and provides further validation that lysozyme resistance is intrinsically linked to the modification of enterococcal cell wall polysaccharides.

59 BASIC BIOLOGICAL SCIENCES↗

Machine learning reveals genes impacting oxidative stress resistance across yeasts

Reactive oxygen species (ROS) are highly reactive molecules encountered by yeasts during routine metabolism and during interactions with other organisms, including host infection. Here, we characterize the variation in resistance to the ROS-inducing compound tert -butyl hydroperoxide across the ancient yeast subphylum Saccharomycotina and use machine learning (ML) to identify gene families whose sizes are predictive of ROS resistance. The most predictive features are enriched in gene families related to cell wall organization and include two reductase gene families. We estimate the quantitative contributions of features to each species’ classification to guide experimental validation and show that overexpression of the old yellow enzyme (OYE) reductase increases ROS resistance in Kluyveromyces lactis , while Saccharomyces cerevisiae mutants lacking multiple mannosyltransferase-encoding genes are hypersensitive to ROS. Altogether, this work provides a framework for how ML can uncover genetic mechanisms underlying trait variation across diverse species and inform trait manipulation for clinical and biotechnological applications.

59 BASIC BIOLOGICAL SCIENCES↗

Identification of carbohydrate gene clusters obtained from in vitro fermentations as predictive biomarkers of prebiotic responses

Prebiotic fibers are non-digestible substrates that modulate the gut microbiome by promoting expansion of microbes having the genetic and physiological potential to utilize those molecules. Although several prebiotic substrates have been consistently shown to provide health benefits in human clinical trials, responder and non-responder phenotypes are often reported. These observations had led to interest in identifying, a priori, prebiotic responders and non-responders as a basis for personalized nutrition. In this study, we conducted in vitro fecal enrichments and applied shotgun metagenomics and machine learning tools to identify microbial gene signatures from adult subjects that could be used to predict prebiotic responders and non-responders. Using short chain fatty acids as a targeted response, we identified genetic features, consisting of carbohydrate active enzymes, transcription factors and sugar transporters, from metagenomic sequencing of in vitro fermentations for three prebiotic substrates: xylooligosacharides, fructooligosacharides, and inulin. A machine learning approach was then used to select substrate-specific gene signatures as predictive features. These features were found to be predictive for XOS responders with respect to SCFA production in an in vivo trial. Our results confirm the bifidogenic effect of commonly used prebiotic substrates along with inter-individual microbial responses towards these substrates. We successfully trained classifiers for the prediction of prebiotic responders towards XOS and inulin with robust accuracy (≥ AUC 0.9) and demonstrated its utility in a human feeding trial. Overall, the findings from this study highlight the practical implementation of pre-intervention targeted profiling of individual microbiomes to stratify responders and non-responders.

59 BASIC BIOLOGICAL SCIENCES↗

Genome Sequence of the Chestnut Blight Fungus Cryphonectria parasitica EP155: A Fundamental Resource for an Archetypical Invasive Plant Pathogen

Cryphonectria parasitica is the causal agent of chestnut blight, a fungal disease that almost entirely eliminated mature American chestnut from North America over a 50-year period. Here, we formally report the genome of C. parasitica EP155 using a Sanger shotgun sequencing approach. After finishing and integration with simple-sequence repeat markers, the assembly was 43.8 Mb in 26 scaffolds (L 50 = 5; N 50 = 4.0Mb). Eight chromosomes are predicted: five scaffolds have two telomeres and six scaffolds have one telomere sequence. In total, 11,609 gene models were predicted, of which 85% show similarities to other proteins. Finally, this genome resource has already increased the utility of a fundamental plant pathogen experimental system through new understanding of the fungal vegetative incompatibility system, with significant implications for enhancing mycovirus-based biological control.

59 BASIC BIOLOGICAL SCIENCES↗

Evaluating the Performance of Random Forest and Iterative Random Forest Based Methods when Applied to Gene Expression Data

Gene-to-gene networks, such as Gene Regulatory Networks (GRN) and Predictive Expression Networks (PEN) capture relationships between genes and are beneficial for use in downstream biological analyses. There exists multiple network inference tools to produce these gene-to-gene networks from matrices of gene expression data. Random Forest-Leave One Out Prediction (RF-LOOP) is a method that has been shown to be efficient at producing these gene-to-gene networks, frequently known as GEne Network Inference with Ensemble of trees (GENIE3). Here we validate that iterative Random Forest-Leave One Out Prediction (iRF-LOOP) produces higher quality networks than GENIE3. We use both synthetic and empirical networks from the Dialogue for Reverse Engineering Assessment and Methods (DREAM) Challenges by Sage Bionetworks, as well as two additional empirical networks created from Arabidopsis thaliana and Populus trichocarpa expression data.

iRF-Loop, expression network, Populus Trichocarpa↗

Evaluating the performance of random forest and iterative random forest based methods when applied to gene expression data

Gene-to-gene networks, such as Gene Regulatory Networks (GRN) and Predictive Expression Networks (PEN) capture relationships between genes and are beneficial for use in downstream biological analyses. There exists multiple network inference tools to produce these gene-to-gene networks from matrices of gene expression data. Random Forest-Leave One Out Prediction (RF-LOOP) is a method that has been shown to be efficient at producing these gene-to-gene networks, frequently known as GEne Network Inference with Ensemble of trees (GENIE3). Random Forest can be replaced in this process by iterative Random Forest (iRF), which performs variable selection and boosting. Here we validate that iterative Random Forest-Leave One Out Prediction (iRF-LOOP) produces higher quality networks than GENIE3 (RF-LOOP). We use both synthetic and empirical networks from the Dialogue for Reverse Engineering Assessment and Methods (DREAM) Challenges by Sage Bionetworks, as well as two additional empirical networks created from Arabidopsis thaliana and Populus trichocarpa expression data.

59 BASIC BIOLOGICAL SCIENCES↗

Transcriptomic and metabolomic changes in postharvest sugarbeet roots reveal widespread metabolic changes in storage and identify genes potentially responsible for respiratory sucrose loss

Endogenous metabolism is primarily responsible for losses in sucrose content and processing quality in postharvest sugarbeet roots. The genes responsible for this metabolism and the transcriptional changes that regulate it, however, are largely unknown. To identify genes and metabolic pathways that participate in postharvest sugarbeet root metabolism and the transcriptional changes that contribute to their regulation, transcriptomic and metabolomic profiles were generated for sugarbeet roots at harvest and after 12, 40 and 120 d storage at 5 and 12°C and gene expression and metabolite concentration changes related to storage duration or temperature were identified. During storage, 8656 genes, or 34% of all expressed genes, and 225 metabolites, equivalent to 59% of detected metabolites, were altered in expression or concentration, indicating extensive transcriptional and metabolic changes in stored roots. These genes and metabolites contributed to a wide range of cellular and molecular functions, with carbohydrate metabolism being the function to which the greatest number of genes and metabolites classified. Because respiration has a central role in postharvest metabolism and is largely responsible for sucrose loss in sugarbeet roots, genes and metabolites involved in and correlated to respiration were identified. Seventy-five genes participating in respiration were differentially expressed during storage, including two bidirectional sugar transporter SWEET17 genes that highly correlated with respiration rate. Weighted gene co-expression network analysis identified 1896 additional genes that positively correlated with respiration rate and predicted a pyruvate kinase gene to be a central regulator or biomarker for respiration rate. Overall, these results reveal the extensive and diverse physiological and metabolic changes that occur in stored sugarbeet roots and identify genes with potential roles as regulators or biomarkers for respiratory sucrose loss.

59 BASIC BIOLOGICAL SCIENCES↗