SEARCH · Search NASA

Results for “GENES”

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.

SourceResearch connectionFrom yearThrough year

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

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Gene expression of functionally-related genes coevolves across fungal species: detecting coevolution of gene expression using phylogenetic comparative methods

Researchers often measure changes in gene expression across conditions to better understand the shared functional roles and regulatory mechanisms of different genes. Analogous to this is comparing gene expression across species, which can improve our understanding of the evolutionary processes shaping the evolution of both individual genes and functional pathways. One area of interest is determining genes showing signals of coevolution, which can also indicate potential functional similarity, analogous to co-expression analysis often performed across conditions for a single species. However, as with any trait, comparing gene expression across species can be confounded by the non-independence of species due to shared ancestry, making standard hypothesis testing inappropriate. We compared RNA-Seq data across 18 fungal species using a multivariate Brownian Motion phylogenetic comparative method (PCM), which allowed us to quantify coevolution between protein pairs while directly accounting for the shared ancestry of the species. Our work indicates proteins which physically-interact show stronger signals of coevolution than randomly-generated pairs. Interactions with stronger empirical and computational evidence also showing stronger signals of coevolution. We examined the effects of number of protein interactions and gene expression levels on coevolution, finding both factors are overall poor predictors of the strength of coevolution between a protein pair. Simulations further demonstrate the potential issues of analyzing gene expression coevolution without accounting for shared ancestry in a standard hypothesis testing framework. Furthermore, our simulations indicate the use of a randomly-generated null distribution as a means of determining statistical significance for detecting coevolving genes with phylogenetically-uncorrected correlations, as has previously been done, is less accurate than PCMs, although is a significant improvement over standard hypothesis testing. These methods are further improved by using a phylogenetically-corrected correlation metric. Our work highlights potential benefits of using PCMs to detect gene expression coevolution from high-throughput omics scale data. This framework can be built upon to investigate other evolutionary hypotheses, such as changes in transcription regulatory mechanisms across species.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2018

Improved fatty acid profiles in seeds of Camelina sativa by artificial microRNA mediated FATB gene suppression

The fatty acid profile of plant oils determines their quality and uses. Saturated fatty acids are often not desirable from the standpoints of nutrition and some industrial applications. Camelina sativa is a re-emerged oilseed crop, however its oil needs to be improved to meet different application requirements. In this study, saturated fatty acids were greatly reduced by down-regulating genes encoding the fatty acyl-ACP thioesterases (FATB). An artificial microRNA (amiFATB) was created by replacing a microRNA sequence in the camelina Csa-miR159a gene with a FATB gene specific sequence. Seed-specific expression of amiFATB caused a 45% reduction of palmitic acid (16:0) and a 38% reduction of stearic acid (18:0) compared to wildtype seeds. The total saturated fatty acid content was decreased by 35% from 14.6% to 9.4% of total fatty acids. When amiFATB was expressed in a high-oleic acid transgenic line, it caused further increased oleic acid content. This work demonstrates that the FATB genes in camelina can be effectively knocked down by an artificial microRNA targeting gene-specific sequences, thus provides an additional tool to improve seed oils for desired properties.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2019

Changes in the Metastatic Properties of MDA-MB-231 Cells after IGFBP6 Gene Knockdown Is Associated with Increased Expression of miRNA Genes Controlling INSR, IGF1R, and CCND1 Genes

Metastatic cascade is associated with the process of epithelial-mesenchymal transition accompanied by changes in cell proliferation, migration, adhesion, and invasiveness mediated by the insulin-like growth factor (IGF) signal pathway. IGFBP6 protein binds IGF and prevents its interaction with receptors. IGFBP6 gene knockdown through RNA-interference inhibits cell migration and increased the rate of proliferation of breast cancer MDA-MB-231 cells. IGFBP6 knockdown cells are characterized by increased expression of MIR100 and MIRLET7A2 genes encoding hsa-miR-100-3p, hsa-miR-100-5p, hsa-let-7a-5p, and hsa-let-7a-2-3p miRNA. The target genes of these microRNAs are IGF2, IGF1R, INSR, and CCND1 associated with IGF signaling pathway and proliferative and migratory activity during the metastatic cascade. A significant decrease in the expression of INSR and CCND1 genes was demonstrated by PCR and microarray analysis.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2019

Deletion of either the regulatory gene ara1 or metabolic gene xki1 in Trichoderma reesei leads to increased CAZyme gene expression on crude plant biomass

Background:Trichoderma reesei is one of the major producers of enzymes for the conversion of plant biomass to sustainable fuels and chemicals. Crude plant biomass can induce the production of CAZymes in T. reesei, but there is limited understanding of how the transcriptional response to crude plant biomass is regulated. In addition, it is unknown whether induction on untreated recalcitrant crude plant biomass (with a large diversity of inducers) can be sustained for longer. We investigated the transcriptomic response of T. reesei to the two industrial feedstocks, corn stover (CS) and soybean hulls (SBH), over time (4 h, 24 h and 48 h), and its regulatory basis using transcription factor deletion mutants (Δxyr1 and Δara1). We also investigated whether deletion of a xylulokinase gene (Δxki1) from the pentose catabolic pathway that converts potential inducers could lead to increased CAZyme gene expression. Results:By analyzing the transcriptomic responses using clustering as well as differential and cumulative expression of plant biomass degrading CAZymes, we found that corn stover induced a broader range and higher expression of CAZymes in T. reesei, while SBH induced more pectinolytic and mannanolytic transcripts. XYR1 was the major TF regulating CS utilization, likely due to the significant amount of d-xylose in this substrate. In contrast, ARA1 had a stronger effect on SBH utilization, which correlates with a higher abundance of l-arabinose in SBH that activates ARA1. Blocking pentose catabolism by deletion of xki1 led to higher expression of CAZyme encoding genes on both substrates at later time points. Surprisingly, this was also observed for Δara1 at later time points. Many of these genes were XYR1 regulated, suggesting that inducers for this regulator accumulated over time on both substrates. Conclusion:Our data demonstrates the complexity of the regulatory system related to plant biomass degradation in T. reesei and the effect the feedstock composition has on this. Furthermore, this dataset provides leads to improve the efficiency of a T. reesei enzyme cocktail, such as by the choice of substrate or by deleting xki1 to obtain higher production of plant biomass degrading CAZymes.

09 BIOMASS FUELS↗

DOE OSTI2023

JGI Plant Gene Atlas: an updateable transcriptome resource to improve functional gene descriptions across the plant kingdom

Abstract Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Thousands of small, novel genes predicted in global phage genomes

Small genes (<150nucleotides) have been systematically overlooked in phage genomes. We employ a large scale comparative genomics approach to predict >40,000 small-gene families in 2.3 million phage genome contigs. We find that small genes in phage genomes are approximately 3-fold more prevalent than in host prokaryotic genomes. Our approach enriches for small genes that are translated in microbiomes, suggesting the small genes identified are coding. More than 9,000 families encode potentially secreted or transmembrane proteins, more than 5,000families encode predicted anti-CRISPR proteins, and more than500families encode predicted antimicrobial proteins. By combining homology and genomic-neighborhood analyses, we reveal substantial novelty and diversity within phage biology, including small phage genes found in multiple host phyla, small genes encoding proteins that play essential roles in host infection, and small genes that share genomic neighborhoods and whose encoded proteins may share related functions.

Fremin, Brayon↗

DOE OSTI2021

The Distinctive Evolution of orfX Clostridium parabotulinum Strains and Their Botulinum Neurotoxin Type A and F Gene Clusters Is Influenced by Environmental Factors and Gene Interactions via Mobile Genetic Elements

Of the seven currently known botulinum neurotoxin-producing species of Clostridium, C. parabotulinum, or C. botulinum Group I, is the species associated with the majority of human botulism cases worldwide. Phylogenetic analysis of these bacteria reveals a diverse species with multiple genomic clades. The neurotoxins they produce are also diverse, with over 20 subtypes currently represented. The existence of different bont genes within very similar genomes and of the same bont genes/gene clusters within different bacterial variants/species indicates that they have evolved independently. The neurotoxin genes are associated with one of two toxin gene cluster types containing either hemagglutinin (ha) genes or orfX genes. These genes may be located within the chromosome or extrachromosomal elements such as large plasmids. Although BoNT-producing C parabotulinum bacteria are distributed globally, they are more ubiquitous in certain specific geographic regions. Notably, northern hemisphere strains primarily contain ha gene clusters while southern hemisphere strains have a preponderance of orfX gene clusters. OrfX C. parabotulinum strains constitute a subset of this species that contain highly conserved bont gene clusters having a diverse range of bont genes. While much has been written about strains with ha gene clusters, less attention has been devoted to those with orfX gene clusters. The recent sequencing of 28 orfX C. parabotulinum strains and the availability of an additional 91 strains for analysis provides an opportunity to compare genomic relationships and identify unique toxin gene cluster characteristics and locations within this species subset in depth. The mechanisms behind the independent processes of bacteria evolution and generation of toxin diversity are explored through the examination of bacterial relationships relating to source locations and evidence of horizontal transfer of genetic material among different bacterial variants, particularly concerning bont gene clusters. Analysis of the content and locations of the bont gene clusters offers insights into common mechanisms of genetic transfer, chromosomal integration, and development of diversity among these genes.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Concentration-response gene expression analysis in zebrafish reveals phenotypically-anchored transcriptional responses to retene

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous environmental contaminants and are associated with human disease. Canonically, many PAHs induce toxicity via activation of the aryl hydrocarbon receptor (AHR) pathway. While the interaction between PAHs and the AHR is well-established, understanding which AHR-regulated transcriptional effects directly result in observable phenotypes and which are adaptive or benign is important to better understand PAH toxicity. Retene is a frequently detected PAH in environmental sampling and has been associated with AHR2-dependent developmental toxicity in zebrafish, though its mechanism of toxicity has not been fully elucidated. To interrogate transcriptional changes causally associated with retene toxicity, we conducted whole-animal RNA sequencing at 48 hours post-fertilization after exposure to eight retene concentrations. The concentrations were selected to produce effects ranging from no phenotype to mortality and malformations in 100% of animals at 5 days post-fertilization. We identified a concentration-response relationship between retene teratogenicity and differential gene expression in both number of DEGs and magnitude of expression change. Elevated expression of cyp1a at retene concentrations below the threshold for teratogenicity suggested that while cyp1a expression is a sensitive biomarker of AHR activation, it may be too sensitive to serve as a biomarker of AHR-dependent teratogenicity. Genes differentially expressed at only non-teratogenic concentrations were enriched for transforming growth factor-ß (TGF-ß) signaling pathway disruption while DEGs identified at only teratogenic concentrations were significantly enriched for response to xenobiotic stimulus and reduction-oxidation reaction activity. DEGs which spanned both non-teratogenic and teratogenic concentrations showed similar disrupted biological processes to those unique to teratogenic concentrations, indicating these processes were disrupted at low exposure concentrations. Gene co-expression network analysis identified several gene modules, including those associated with PAHs and AHR2 activation. One, Module 7, was strongly enriched for AHR2-associated genes and contained the strongest responses to retene. Benchmark concentration (BMC) of Module 7 genes identified a median BMC of 7.5 µM, nearly the highest retene concentration with no associated teratogenicity, supporting the hypothesis that Module 7 genes are largely responsible for retene toxicity.

Toxin, Zebrafish, Retene, Transcriptomics, network↗

DOE OSTI2026

Data for FUN-PROSE: A Deep Learning Approach to Predict Condition-Specific Gene Expression in Fungi

mRNA levels of all genes in a genome is a critical piece of information defining the overall state of the cell in a given environmental condition. Being able to reconstruct such condition-specific expression in fungal genomes is particularly important to metabolically engineer these organisms to produce desired chemicals in industrially scalable conditions. Most previous deep learning approaches focused on predicting the average expression levels of a gene based on its promoter sequence, ignoring its variation across different conditions. Here we present FUN-PROSE—a deep learning model trained to predict differential expression of individual genes across various conditions using their promoter sequences and expression levels of all transcription factors. We train and test our model on three fungal species and get the correlation between predicted and observed condition-specific gene expression as high as 0.85. We then interpret our model to extract promoter sequence motifs responsible for variable expression of individual genes. We also carried out input feature importance analysis to connect individual transcription factors to their gene targets. A sizeable fraction of both sequence motifs and TF-gene interactions learned by our model agree with previously known biological information, while the rest corresponds to either novel biological facts or indirect correlations.

Genomics↗

DOE OSTI2018

Annotation of the Corymbia terpene synthase gene family shows broad conservation but dynamic evolution of physical clusters relative to Eucalyptus

Terpenes are economically and ecologically important phytochemicals. Their synthesis is controlled by the terpene synthase (TPS) gene family, which is highly diversified throughout the plant kingdom. The plant family Myrtaceae are characterised by especially high terpene concentrations, and considerable variation in terpene profiles. Many Myrtaceae are grown commercially for terpene products including the eucalypts Corymbia and Eucalyptus. Eucalyptus grandis has the largest TPS gene family of plants currently sequenced, which is largely conserved in the closely related E. globulus. However, the TPS gene family has been well studied only in these two eucalypt species. The recent assembly of two Corymbia citriodora subsp. variegata genomes presents an opportunity to examine the conservation of this important gene family across more divergent eucalypt lineages. Manual annotation of the TPS gene family in C. citriodora subsp. variegata revealed a similar overall number, and relative subfamily representation, to that previously reported in E. grandis and E. globulus. Many of the TPS genes were in physical clusters that varied considerably between Eucalyptus and Corymbia, with several instances of translocation, expansion/contraction and loss. Notably, there was greater conservation in the subfamilies involved in primary metabolism than those involved in secondary metabolism, likely reflecting different selective constraints. The variation in cluster size within subfamilies and the broad conservation between the eucalypts in the face of this variation are discussed, highlighting the potential contribution of selection, concerted evolution and stochastic processes. Finally, these findings provide the foundation to better understand terpene evolution within the ecologically and economically important Myrtaceae.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Identification and selection of optimal reference genes for qPCR-based gene expression analysis in Fucus distichus under various abiotic stresses

Quantitative gene expression analysis is an important tool in the scientist’s belt. The identification of evenly expressed reference genes is necessary for accurate quantitative gene expression analysis, whether by traditional RT-PCR (reverse-transcription polymerase chain reaction) or by qRT-PCR (quantitative real-time PCR; qPCR). In the Stramenopiles (the major line of eukaryotes that includes brown algae) there is a noted lack of known reference genes for such studies, largely due to the absence of available molecular tools. Here we present a set of nine reference genes (Elongation Factor 1 alpha (EF1A), Elongation Factor 2 alpha (EF2A), Elongation Factor 1 beta (EF1B), 14-3-3 Protein, Ubiquitin Conjugating Enzyme (UBCE2), Glyceraldehyde-3-phosphate Dehydrogenase (GAPDH), Actin Related Protein Complex (ARP2/3), Ribosomal Protein (40s; S23), and Actin) for the brown alga Fucus distichus . These reference genes were tested on adult sporophytes across six abiotic stress conditions (desiccation, light and temperature modification, hormone addition, pollutant exposure, nutrient addition, and wounding). Suitability of these genes as reference genes was quantitatively evaluated across conditions using standard methods and the majority of the tested genes were evaluated favorably. However, we show that normalization genes should be chosen on a condition-by-condition basis. We provide a recommendation that at least two reference genes be used per experiment, a list of recommended pairs for the conditions tested here, and a procedure for identifying a suitable set for an experimenter’s unique design. With the recent expansion of interest in brown algal biology and accompanied molecular tools development, the variety of experimental conditions tested here makes this study a valuable resource for future work in basic biology and understanding stress responses in the brown algal lineage.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2019

Determining the Reliability of Measuring Mercury Cycling Gene Abundance with Correlations with Mercury and Methylmercury Concentrations

Methylmercury (MeHg) is a bioaccumulative toxic contaminant in many ecosystems but factors governing its production are poorly understood. Recent work has shown that the anaerobic microbial conversion of mercury (Hg) to MeHg requires the Hg-methylation genes hgcAB and that these genes can be used as biomarkers in PCR-based estimators of Hg-methylator abundance. In an effort to determine reliable methods for assessing hgcA abundance and diversity and linking it to MeHg concentrations, multiple approaches were compared including metagenomic shotgun sequencing, 16S rRNA gene pyrosequencing and cloning/sequencing hgcAB gene products. Hg-methylator abundance was also determined by quantitative hgcA qPCR amplification and metaproteomics for comparison to the above measurements. Samples from eight Hg-contaminated sites were examined over a range of total Hg (HgT; 0.03–14 mg kg-1 wet wt. soil) and MeHg (0.05–27 µg kg-1 wet wt. soil) concentrations. In the metagenome and amplicon sequencing of hgcAB diversity, the Deltaproteobacteria were the dominant Hg-methylators, while Firmicutes and methanogenic Archaea were typically ~50% less abundant. This was consistent with metaproteomics estimates where the Deltaproteobacteria were consistently higher. The 16S rRNA gene pyrosequencing did not have sufficient resolution to identify hgcAB+ species. Metagenomic and hgcAB results were similar for Hg-methylator diversity and showed a tractable relationship between Hg-methylating bacteria and soil Hg concentrations among similar environmental samples. However, clade-specific qPCR-based approaches for hgcA are only appropriate when comparing the abundance of a particular clade across various samples.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Gene-by-environment interactions influence the fitness cost of gene copy-number variation in yeast

Abstract Variation in gene copy number can alter gene expression and influence downstream phenotypes; thus copy-number variation provides a route for rapid evolution if the benefits outweigh the cost. We recently showed that genetic background significantly influences how yeast cells respond to gene overexpression, revealing that the fitness costs of copy-number variation can vary substantially with genetic background in a common-garden environment. But the interplay between copy-number variation tolerance and environment remains unexplored on a genomic scale. Here, we measured the tolerance to gene overexpression in four genetically distinct Saccharomyces cerevisiae strains grown under sodium chloride stress. Overexpressed genes that are commonly deleterious during sodium chloride stress recapitulated those commonly deleterious under standard conditions. However, sodium chloride stress uncovered novel differences in strain responses to gene overexpression. West African strain NCYC3290 and North American oak isolate YPS128 are more sensitive to sodium chloride stress than vineyard BC187 and laboratory strain BY4743. Consistently, NCYC3290 and YPS128 showed the greatest sensitivities to overexpression of specific genes. Although most genes were deleterious, hundreds were beneficial when overexpressed—remarkably, most of these effects were strain specific. Few beneficial genes were shared between the sodium chloride-sensitive isolates, implicating mechanistic differences behind their sodium chloride sensitivity. Transcriptomic analysis suggested underlying vulnerabilities and tolerances across strains, and pointed to natural copy-number variation of a sodium export pump that likely contributes to strain-specific responses to overexpression of other genes. Our results reveal extensive strain-by-environment interactions in the response to gene copy-number variation, raising important implications for the accessibility of copy-number variation-dependent evolutionary routes under times of stress.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

Use of virus‐induced gene silencing to characterize genes involved in modulating hypersensitive cell death in maize

Abstract Plant disease resistance proteins (R‐proteins) detect specific pathogen‐derived molecules, triggering a defence response often including a rapid localized cell death at the point of pathogen penetration called the hypersensitive response (HR). The maize Rp1‐D21 gene encodes a protein that triggers a spontaneous HR causing spots on leaves in the absence of any pathogen. Previously, we used fine mapping and functional analysis in a Nicotiana benthamiana transient expression system to identify and characterize a number of genes associated with variation in Rp1‐D21 ‐induced HR. Here we describe a system for characterizing genes mediating HR, using virus‐induced gene silencing (VIGS) in a maize line carrying Rp1‐D21 . We assess the roles of 12 candidate genes. Three of these genes, SGT1 , RAR1 , and HSP90 , are required for HR induced by a number of R‐proteins across several plant–pathogen systems. We confirmed that maize HSP90 was required for full Rp1‐D21 ‐induced HR. However, suppression of SGT1 expression unexpectedly increased the severity of Rp1‐D21 ‐induced HR while suppression of RAR1 expression had no measurable effect. We confirmed the effects on HR of two genes we had previously validated in the N. benthamiana system, hydroxycinnamoyltransferase and caffeoyl CoA O‐methyltransferase . We further showed the suppression the expression of two previously uncharacterized, candidate genes, IQ calmodulin binding protein ( IQM3 ) and vacuolar protein sorting protein 37 , suppressed Rp1‐D21 ‐induced HR. This approach is an efficient way to characterize the roles of genes modulating the hypersensitive defence response and other dominant lesion phenotypes in maize.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Population‐level gene expression can repeatedly link genes to functions in maize

SUMMARY Transcriptome‐wide association studies (TWAS) can provide single gene resolution for candidate genes in plants, complementing genome‐wide association studies (GWAS) but efforts in plants have been met with, at best, mixed success. We generated expression data from 693 maize genotypes, measured in a common field experiment, sampled over a 2‐h period to minimize diurnal and environmental effects, using full‐length RNA‐seq to maximize the accurate estimation of transcript abundance. TWAS could identify roughly 10 times as many genes likely to play a role in flowering time regulation as GWAS conducted data from the same experiment. TWAS using mature leaf tissue identified known true‐positive flowering time genes known to act in the shoot apical meristem, and trait data from a new environment enabled the identification of additional flowering time genes without the need for new expression data. eQTL analysis of TWAS‐tagged genes identified at least one additional known maize flowering time gene through trans ‐eQTL interactions. Collectively these results suggest the gene expression resource described here can link genes to functions across different plant phenotypes expressed in a range of tissues and scored in different experiments.

Torres‐Rodríguez, J. Vladimir↗

DOE OSTI2021

Hypergraph Models of Biological Networks to Identify Genes Critical to Pathogenic Viral Response

Motivation: Representing biological networks as graphs is a powerful approach to reveal underlying patterns, signatures, and critical components from high-throughput biomolecular data. However, graphs do not natively capture the multi-way relationships present among genes and proteins in biological systems such as protein complexes, metabolic reactions, and signal transduction pathways. Hypergraphs are generalizations of graphs that naturally model multi-way interactions in data, and we therefore seek to understand how they can more faithfully identify, and potentially predict, complex relationships in genomic expression data sets. Results: We compiled a novel data set of transcriptional host response to pathogenic viral infections and formulated relationships between genes as a hypergraph where hyperedges are differentially expressed genes and vertices represent conditions. We find that hypergraph betweenness centrality is a superior method for identification of genes important to viral response when compared with graph centrality. Our results demonstrate the utility of using hypergraphs to represent complex biological systems, and highlight potentially interesting biological results about host response to highly pathogenic viruses.

systems biology, hypergraph, viral infection, biol↗

DOE OSTI2020

FeGenie: a comprehensive tool for the identification of iron genes and iron gene neighborhoods in genomes and metagenome assemblies

Iron is a micronutrient for nearly all life on Earth. It can be used as an electron donor and electron acceptor by iron-oxidizing and iron-reducing microorganisms and is used in a variety of biological processes, including photosynthesis and respiration. While it is the fourth most abundant metal in the Earth’s crust, iron is often limiting for growth in oxic environments because it is readily oxidized and precipitated. Much of our understanding of how microorganisms compete for and utilize iron is based on laboratory experiments. However, the advent of next-generation sequencing and surge in publicly available sequence data has made it possible to probe the structure and function of microbial communities in the environment. To bridge the gap between our understanding of iron acquisition, iron redox cycling, iron storage, and magnetosome formation in model microorganisms and the plethora of sequence data available from environmental studies, we have created a comprehensive database of hidden Markov models (HMMs) based on genes related to iron acquisition, storage, and reduction/oxidation in Bacteria and Archaea. Along with this database, we present FeGenie, a bioinformatics tool that accepts genome and metagenome assemblies as input and uses our comprehensive HMM database to annotate provided datasets with respect to iron-related genes and gene neighborhood. An important contribution of this tool is the efficient identification of genes involved in iron oxidation and dissimilatory iron reduction, which have been largely overlooked by standard annotation pipelines. We validated FeGenie against a selected set of 28 isolate genomes and showcase its utility in exploring iron genes present in 27 metagenomes, 4 isolate genomes from human oral biofilms, and 17 genomes from candidate organisms, including members of the candidate phyla radiation. We show that FeGenie accurately identifies iron genes in isolates. Furthermore, analysis of metagenomes using FeGenie demonstrates that the iron gene repertoire and abundance of each environment is correlated with iron richness. While this tool will not replace the reliability of culture-dependent analyses of microbial physiology, it provides reliable predictions derived from the most up-to-date genetic markers. FeGenie’s database will be maintained and continually updated as new genes are discovered.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2018

Weighted gene co-expression network analysis reveals potential genes involved in early metamorphosis process in sea cucumber Apostichopus japonicus

Highlights: • Gene expression profiles of sea cucumber were systematically examined. • Gene co-expression modules of sea cucumber were analyzed for the first time. • The module MEdarkmagenta is highly related to the early metamorphosis process. • Genes associated with the early metamorphosis process were identified. Sea cucumbers, one main class of Echinoderms, have a very fast and drastic metamorphosis process during their development. However, the molecular basis under this process remains largely unknown. Here we systematically examined the gene expression profiles of Japanese common sea cucumber (Apostichopus japonicus) for the first time by RNA sequencing across 16 developmental time points from fertilized egg to juvenile stage. Based on the weighted gene co-expression network analysis (WGCNA), we identified 21 modules. Among them, MEdarkmagenta was highly expressed and correlated with the early metamorphosis process from late auricularia to doliolaria larva. Furthermore, gene enrichment and differentially expressed gene analysis identified several genes in the module that may play key roles in the metamorphosis process. Our results not only provide a molecular basis for experimentally studying the development and morphological complexity of sea cucumber, but also lay a foundation for improving its emergence rate.

60 APPLIED LIFE SCIENCES↗