EDGE Bioinformatics
This report gives an overview of Empowering the Development of Genomics Expertise (EDGE) Bioinformatics software.
59 BASIC BIOLOGICAL SCIENCES↗
SEARCH · Search NASA
Results for “EDGE Bioinformatics”
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.
Modularization of EDGE Workflows Using Nextflow: Improving the Efficiency and Maintainability of Bioinformatics Software
EDGE is a bioinformatics platform developed in 2016 by researchers at Los Alamos National Laboratory (LANL) to facilitate the analysis of next-generation sequencing data by researchers with varying levels of experience in bioinformatics (Li et al., 2017). Users with single-end, paired-end or long-read sequencing data can provide their reads as input to EDGE and select the combination of workflows to run that are most useful for their research (e.g., quality control of reads, genome assembly, or the taxonomic classification of input reads). Table 1 summarizes the modules available in EDGE. EDGE is available as a web platform at https://edgebioinformatics.org, as installable source code maintained on GitHub under a GPLv3 license, and as a publicly hosted Docker image.
59 BASIC BIOLOGICAL SCIENCES↗
Transboundary determinants of avian zoonotic infectious diseases: challenges for strengthening research capacity and connecting surveillance networks
As the climate changes, global systems have become increasingly unstable and unpredictable. This is particularly true for many disease systems, including subtypes of highly pathogenic avian influenzas (HPAIs) that are circulating the world. Ecological patterns once thought stable are changing, bringing new populations and organisms into contact with one another. Wild birds continue to be hosts and reservoirs for numerous zoonotic pathogens, and strains of HPAI and other pathogens have been introduced into new regions via migrating birds and transboundary trade of wild birds. With these expanding environmental changes, it is even more crucial that regions or counties that previously did not have surveillance programs develop the appropriate skills to sample wild birds and add to the understanding of pathogens in migratory and breeding birds through research. For example, little is known about wild bird infectious diseases and migration along the Mediterranean and Black Sea Flyway (MBSF), which connects Europe, Asia, and Africa. Focusing on avian influenza and the microbiome in migratory wild birds along the MBSF, this project seeks to understand the determinants of transboundary disease propagation and coinfection in regions that are connected by this flyway. Through the creation of a threat reduction network for avian diseases (Avian Zoonotic Disease Network, AZDN) in three countries along the MBSF (Georgia, Ukraine, and Jordan), this project is strengthening capacities for disease diagnostics; microbiomes; ecoimmunology; field biosafety; proper wildlife capture and handling; experimental design; statistical analysis; and vector sampling and biology. Here, we cover what is required to build a wild bird infectious disease research and surveillance program, which includes learning skills in proper bird capture and handling; biosafety and biosecurity; permits; next generation sequencing; leading-edge bioinformatics and statistical analyses; and vector and environmental sampling. Creating connected networks for avian influenzas and other pathogen surveillance will increase coordination and strengthen biosurveillance globally in wild birds.
54 ENVIRONMENTAL SCIENCES↗
Does Collection Time Bias the Ecology of Cleanroom Air Samples?
Microbial monitoring of astromaterials collections has taken on increased importance with the return of biologically sensitive samples from the asteroids Ryugu and Bennu and the initiation of the Mars Sample Return Program. Terrestrial bacteria and fungi can alter the mineralogy and organic composition of our collections causing irreversible contamination of pristine samples and increasing the risk of false positives for life detection measurements. NASA has conducted routine microbial monitoring of its existing collections since 20181. Initial monitoring focused on surface samples collected with foam swabs. Although, airborne microbiology is often decoupled from surface microbiology in the built environment2 culture-based air sampling techniques like impactors were not compliant with existing contamination control requirements. Bringing organic rich media, gelatin or liquids into curation cleanrooms presents an unacceptable risk to pristine samples. In 2022 NASA purchased a materials complaint air sampler and began collecting air samples from the cleanrooms in addition to surface samples3. The new instrument uses an electret filter to collect samples that are suitable for cultivating organisms or for direct DNA sequencing. Preliminary DNA sequencing results appeared to indicate that longer sampling times biased the microbial community in favor of hearty, spore-forming bacteria3. We present the results of a study comparing overnight sampling (17 hours) to short (1 hour) sampling of unoccupied curation cleanrooms. The results will help us optimize our monitoring protocols and develop a more detailed inventory of the ecology of astromaterials curation cleanrooms. Methods: We analyzed 72 paired air samples from six different cleanrooms including the meteorite processing lab (ISO 7 equivalent, 16 samples), the lunar lab (ISO 6 equivalent, 10 samples), the stardust lab (ISO 5 equivalent 14 samples), the OSIRIS-REx lab (ISO 5 equivalent, 12 samples), the Hayabusa2 lab (ISO 5 equivalent, 14 samples), and the Genesis lab (ISO 4 equivalent, 6 samples). All the samples were collected with an InnovaPrep Bobcat air sampler operating at a sampling rate of 200 L/min. The sampler operates for 5 minutes out of every 20 minute period. Half of the samples were collected by filtering 3,000L (15 min. of active sampling) of air across an electret filter for one hour. The rest of the samples were collected by filtering approximately 51,000 L air across the filter overnight (~17 hours, 255 min. of active sampling). Cells were eluted from the filter using 6-7 ml of pressurized 0.15% tween 20 in PBS (phosphate buffered saline). This liquid was used to cultivate bacteria according to previously published methods1,4,5 and for DNA extraction and next generation sequencing. DNA was extracted with a Qiagen MagAttract PowerMicrobiome kit6. To identify bacteria and archaea, the 16S rRNA gene was amplified using Earth Microbiome primers for the V4 region 7. The amplified DNA was sequenced on an Illumina MiSeq using a V3 reagent kit. The resulting sequences were processed using DADA2 and QIIME2 as implemented on the EDGE bioinformatics platform8–10. Results: Only two of the 72 samples had no amplifiable DNA. Amplified DNA concentrations ranged from 2.67 – 0.272 ng/µl. The median concentration of amplified DNA for the 1 hour samples was 0.770 ± 0.368 ng/µl. The median concentration of amplified DNA for the overnight samples was 0.877 ± 0.434 ng/µl. On average the overnight samples had slightly more sequences (58,960 vs. 59,456) and ASV’s (amplicon sequence variants) (60 vs 64.5) than the one hour samples, but these differences are not statistically significant. The most abundant ASV in every sample mapped to the genus Cupravidus. ASV’s mapping to the genuses Bacillus, Schlegelella, Thermus, and Staphylococcus were also common. Discussion and Future Work: Alpha diversity statistics like Shannon Entropy and Faith Phylogenetic Diversity are used to describe the diversity of organisms in a single sample. If a longer sampling time was biasing the data, we would expect to see a change in these diversity statistics vs. sample time. However, we did not observe this in our data. The median Shannon entropy was slightly higher for the overnight samples (3.773 vs 3.611) as was the Faith Phylogenetic Diversity (4.042 vs 3.596), but both values were within a standard deviation of each other for the two sampling times (Fig. 1). It is unlikely, that the longer sampling time is introducing bias into our data. We do observe a significant decrease in diversity when comparing the air samples by lab. The Genesis lab (ISO 4 equivalent) has a lower median number of ASV’s (45.5) than the other labs (62). Median values for Shannon Entropy (3.717 vs. 3.430) and Faith Phylogenetic Diversity (3.796 vs. 3.548) are also lower for Genesis, but those values are with one standard deviation of each other for the different sampling times. This is consistent with previous culture-based results suggesting that the environment in cleanrooms tends to select for a core group of organisms capable of surviving under dry, low nutrient, conditions. The presence of the ASV’s mapping to Cupravidus and Thermus in our sequencing blanks and controls suggests that several of the most common organisms in our samples represent contaminants from the reagents used to perform the DNA extractions and sequencing. Further work is needed to identify these contaminants, remove them from our data and recalculate the diversity statistics. This is a systematic error. Therefore, we do not expect removing the sequencing contaminants to change our conclusions. Longer air sample collection times appear to result in slightly higher diversity and do not bias the results towards “hardy” bacteria like spore-formers. Based on these preliminary results we conclude that sampling at least 3,000 liters of air is sufficient to capture the microbial diversity of cleanrooms, and that air samples can also be collected overnight without negatively impacting diversity. These results allow us to be flexible when designing microbial monitoring plans so that they do not interfere with routine lab activity. References: 1. Regberg, A. B. et al. 49th Lunar and Planetary Science Conference (2018). 2. The United States Pharmacopeial Convention. USP General Chapter <1116> (2013). 3. Regberg, A. B., et al. 54th Lunar and Planetary Science Conference (2023). 4. Regberg, A. B. et al. 53rd Lunar and Planetary Science Conference ( 2022). 5. Davis, R. E.,et al. 50th Lunar and Planetary Science Conference (2019). 6. Qiagen. MagAttract® PowerMicrobiome® DNA/RNA EP Kit Handbook. (2018). 7. Walters, W. et al. mSystems 1, (2015). 8. Callahan, B. J. et al. Nat. Methods 13, 581–583 (2016). 9. Hall, M. & Beiko, R. G. Microbiome Analysis: Methods and Protocols113–129 (Springer, 2018). 10. Philipson, C. et al. Bio-Protoc. 7, e2622 (2017).
A. B. Regberg↗
Does the International Space Station Leak DNA? Preliminary Results from the ISS External Microorganisms Payload
Existing crewed spacecraft like the ISS (International Space Station) leak by design. The ISS routinely releases gas to maintain life support systems and when astronauts exit the station to perform space walks. The chemical component of this leakage is well characterized, but the biological components are not. The ISS is not subject to planetary protection requirements, but planned missions to Mars will use similar systems and will be subject to planetary protection requirements. If detectable microorganisms are escaping through vents and or airlocks we may need to redesign our crewed habitats to minimize this type of contamination. To test the hypothesis that microorganisms from inside ISS are detectable on exterior surfaces an astronaut used the ISS External Microorganisms sampling kit (Rucker et al. 2018) to sample exterior surfaces of the ISS during an EVA (Extra Vehicular Activity) in January of 2025. These samples were returned to Earth for DNA extraction and sequencing. We successfully, extracted and sequenced bacterial, fungal and viral DNA from these samples that was not present in the negative controls. These results should help NASA refine the planetary protection requirements for crewed missions. Methods: The samples were collected using sterile, DNA free, buccal swabs (23 mm. diameter) housed in custom canisters. Each canister uses a 0.2 μm Teflon filter to maintain sterility as the caddy, holding 8 swabs moves in and out of vacuum. The astronaut sampled the: 1) airlock vestibule, 2) airlock thermal cover, 3) a gap in the micrometeorite shielding near the airlock, 4) a handrail near the airlock, 5) the Carbon Dioxide Removal Assembly vent, and 6) the Vacuum Exhaust System vent. The seventh swab was exposed to vacuum during the EVA without touching it to a surface. The eighth swab, a negative control, was not opened until the caddy returned to Earth. DNA was extracted from the swabs using a QIamp UCP Pathogen kit and prepared for sequencing on an Aviti (Element Biosciences) sequencer (Arslan et al. 2024). The resulting sequences were analyzed using the EDGE Bioinformatics platform (Li et al. 2017). The sequences were analyzed individually using tools like BLAST, GOTTCHA2, Kraken2, and PanGIA. The data were also assembled into metagenome assembled genomes) using tools like CONCOCT, MaxBin2 and MetaBAT2. Results: We successfully extracted and sequenced bacterial, archaeal, fungal and viral DNA from all seven samples. The handrail swab had the lowest number of reads (768,651) and the airlock thermal cover had the highest number of reads (8,819,230). These samples contain DNA from human associated bacteria (e.g. Crynebacterium riegelii ), fungi (.e.g. Penicillium rubens ), and viruses (e.g Alphapapillomavirus ). Conclusion: Preliminary interpretation suggest that the airlock and the space suits themselves are the largest sources of contaminant DNA. Most if not all of the DNA is from organisms known to be present inside the ISS. Vents attached to life support systems may be a lesser source of biological contamination. Further analysis should help NASA address planetary protection knowledge gaps for crewed missions.
Aaron B Regberg↗
Standardized and accessible multi-omics bioinformatics workflows through the NMDC EDGE resource
Not Available
59 BASIC BIOLOGICAL SCIENCES↗
Enabling Open and Interoperable Science: Multi-Omics Data Processing Platform with NASA GeneLab Standardized Bioinformatics Workflows for Space and Earth Research
Multi-omics biological data continues to be generated at an astounding pace. Genomics, transcriptomics, metabolomics, and proteomics, or collectively known as multi-omics data, are used to assess biological functions, and provide invaluable insights into human, animal, plant, and environmental health both on Earth and in Space. Despite the abundance of these valuable data, the need for bioinformatics expertise, particularly as it relates to the niche filed of space biology, and a lack of accessible resources for processing these data limit their usefulness in deriving biological insights. The NASA Open Science Data Repository (OSDR) provides access to omics data from various spaceflight and analog studies. To enhance the accessibility and reusability of these data, GeneLab (part of OSDR) designs and implements standardized, community-driven, open-source bioinformatics workflows to transform raw omics data into standardized processed data. Currently, GeneLab-processed data from hundreds of space studies have been reused for meta-analyses. This has led to new insights and scientific publications that extend beyond the initial research, thereby enriching our understanding of molecular-scale biological responses to the space environment. To make these bioinformatics workflows open and accessible, GeneLab teamed up with DOE-funded initiatives, including the National Microbiome Data Collaborative (NMDC), to create the NASA EDGE [Empowering the Development of Genomics Expertise] Bioinformatics web-based platform. NASA EDGE utilizes shared compute resources to run the GeneLab standardized bioinformatics workflows, which eliminates the need for researchers to have their own high performance computing cluster. The web-based platform makes complicated biological analyses incredibly easy to perform, thus expanding the reach of these analyses to bioinformatics novices, students, and even citizen scientists enabling them to contribute to scientific discoveries and progress. The authors will demonstrate how the NASA EDGE platform can be used to process microbial omics data hosted on OSDR as well as user-generated omics datasets using GeneLab’s standard workflows.
Amanda M. Saravia-Butler↗
Sequencing and analysis of 131 SARS-CoV-2 isolates in previously sampled and unsampled regions of Jordan from 2020 to 2023
The Hashemite Kingdom of Jordan remains an understudied country for next generation sequencing analysis of SARS-CoV-2 genomes collected during the 2019 pandemic. Here we provide 131 additional reference genomes collected between 2020–2023 from SARS-CoV-2-positive patients across Jordan. Phylogenetic analysis supports existing pandemic narratives of changing clade dominance over time and adds genomes in novel Jordanian locations and timepoints to make Jordan SARS-CoV-2 databases more comprehensive. Samples from the less-sequenced cities of Ajloun, Jaresh, Karak, and Madaba identified previously unreported lineages while Amman, Irbid, and Zarqa have existing sequencing efforts bolstered. Despite many incomplete patient records and a relatively small sample size, we observe interesting symptom patterns that support existing global and Jordanian pandemic narratives. We note how in-country COVID-19 pandemic genomic studies showcase Jordan’s efforts to expand next generation sequencing capabilities, especially through the leveraging of EDGE COVID-19, a bioinformatics platform for performing rapid, batched analysis of SARS-CoV-2 sequencing that streamlines sample processing prepared from a network of hospital locations.
60 APPLIED LIFE SCIENCES↗
HDBind: encoding of molecular structure with hyperdimensional binary representations
Traditional methods for identifying “hit” molecules from a large collection of potential drug-like candidates rely on biophysical theory to compute approximations to the Gibbs free energy of the binding interaction between the drug and its protein target. These approaches have a significant limitation in that they require exceptional computing capabilities for even relatively small collections of molecules. Increasingly large and complex state-of-the-art deep learning approaches have gained popularity with the promise to improve the productivity of drug design, notorious for its numerous failures. However, as deep learning models increase in their size and complexity, their acceleration at the hardware level becomes more challenging. Hyperdimensional Computing (HDC) has recently gained attention in the computer hardware community due to its algorithmic simplicity relative to deep learning approaches. The HDC learning paradigm, which represents data with high-dimension binary vectors, allows the use of low-precision binary vector arithmetic to create models of the data that can be learned without the need for the gradient-based optimization required in many conventional machine learning and deep learning methods. This algorithmic simplicity allows for acceleration in hardware that has been previously demonstrated in a range of application areas (computer vision, bioinformatics, mass spectrometery, remote sensing, edge devices, etc.). To the best of our knowledge, our work is the first to consider HDC for the task of fast and efficient screening of modern drug-like compound libraries. We also propose the first HDC graph-based encoding methods for molecular data, demonstrating consistent and substantial improvement over previous work. We compare our approaches to alternative approaches on the well-studied MoleculeNet dataset and the recently proposed LIT-PCBA dataset derived from high quality PubChem assays. We demonstrate our methods on multiple target hardware platforms, including Graphics Processing Units (GPUs) and Field Programmable Gate Arrays (FPGAs), showing at least an order of magnitude improvement in energy efficiency versus even our smallest neural network baseline model with a single hidden layer. Our work thus motivates further investigation into molecular representation learning to develop ultra-efficient pre-screening tools. We make our code publicly available at https://github.com/LLNL/hdbind.
59 BASIC BIOLOGICAL SCIENCES↗
Using Machine Learning for Quantum Annealing Accuracy Prediction
Quantum annealers, such as the device built by D-Wave Systems, Inc., offer a way to compute solutions of NP-hard problems that can be expressed in Ising or quadratic unconstrained binary optimization (QUBO) form. Although such solutions are typically of very high quality, problem instances are usually not solved to optimality due to imperfections of the current generations quantum annealers. In this contribution, we aim to understand some of the factors contributing to the hardness of a problem instance, and to use machine learning models to predict the accuracy of the D-Wave 2000Q annealer for solving specific problems. We focus on the maximum clique problem, a classic NP-hard problem with important applications in network analysis, bioinformatics, and computational chemistry. By training a machine learning classification model on basic problem characteristics such as the number of edges in the graph, or annealing parameters, such as the D-Wave’s chain strength, we are able to rank certain features in the order of their contribution to the solution hardness, and present a simple decision tree which allows to predict whether a problem will be solvable to optimality with the D-Wave 2000Q. We extend these results by training a machine learning regression model that predicts the clique size found by D-Wave.
97 MATHEMATICS AND COMPUTING↗
cuAlign: Scalable Network Alignment on GPU Accelerators
Given two graphs, the objective of network alignment is to find the best one-to-one mapping of vertices in one graph (??) to vertices in the other (??), such that the number of overlaps is maximized. We say that edges(??, ??) ???and(??', ??') ??? are overlapped if ?? is mapped to ??' and ?? is mapped to??'. Network alignment is an important optimization problem with several applications in bioinformatics, computer vision and ontology matching. Since it is an NP-hard problem, efficient heuristics and scalable implementations are necessary. In this work, we introduce a new framework that combines the concepts of intra-network proximity using vertex embedding,Belief Propagation (BP) and approximate weighted matching, and provides qualitative improvements up to22%over state-of-the-art approaches. We also provide scalable implementations on GPU accelerators, demonstrating up to19×speedup for Belief Propagation and 3× speedup for approximate weighted matching relative to previous multithreaded implementation. A combination of combinatorial and algebraic kernels within the network alignment algorithm poses significant hurdles for parallelization. Load imbalance and irregular DRAM traffic limit achievable performance on GPUs. Our novel approach identifies and exploits unique structural proper-ties of the BP-based algorithm and employs code fusion to reduce data movement between different steps of the algorithm. Using a diverse set of inputs, we demonstrate qualitative improvements of our algorithms, and performance gains of our GPU-accelerated implementation. We believe that our work will enable algorithmic improvements and practical applications of network alignment.
Xiang, Lizhi↗
Scalable deep text comprehension for Cancer surveillance on high-performance computing
Background: Deep Learning (DL) has advanced the state-of-the-art capabilities in bioinformatics applications which has resulted in trends of increasingly sophisticated and computationally demanding models trained by larger and larger data sets. This vastly increased computational demand challenges the feasibility of conducting cutting-edge research. One solution is to distribute the vast computational workload across multiple computing cluster nodes with data parallelism algorithms. In this study, we used a High-Performance Computing environment and implemented the Downpour Stochastic Gradient Descent algorithm for data parallelism to train a Convolutional Neural Network (CNN) for the natural language processing task of information extraction from a massive dataset of cancer pathology reports. We evaluated the scalability improvements using data parallelism training and the Titan supercomputer at Oak Ridge Leadership Computing Facility. To evaluate scalability, we used different numbers of worker nodes and performed a set of experiments comparing the effects of different training batch sizes and optimizer functions.Results: We found that Adadelta would consistently converge at a lower validation loss, though requiring over twice as many training epochs as the fastest converging optimizer, RMSProp. The Adam optimizer consistently achieved a close 2nd place minimum validation loss significantly faster; using a batch size of 16 and 32 allowed the network to converge in only 4.5 training epochs.Conclusions: We demonstrated that the networked training process is scalable across multiple compute nodes communicating with message passing interface while achieving higher classification accuracy compared to a traditional machine learning algorithm.
Qiu, John X.↗
Temporal network alignment via GoT-WAVE
Abstract Motivation Network alignment (NA) finds conserved regions between two networks. NA methods optimize node conservation (NC) and edge conservation. Dynamic graphlet degree vectors are a state-of-the-art dynamic NC measure, used within the fastest and most accurate NA method for temporal networks: DynaWAVE. Here, we use graphlet-orbit transitions (GoTs), a different graphlet-based measure of temporal node similarity, as a new dynamic NC measure within DynaWAVE, resulting in GoT-WAVE. Results On synthetic networks, GoT-WAVE improves DynaWAVE’s accuracy by 30% and speed by 64%. On real networks, when optimizing only dynamic NC, the methods are complementary. Furthermore, only GoT-WAVE supports directed edges. Hence, GoT-WAVE is a promising new temporal NA algorithm, which efficiently optimizes dynamic NC. We provide a user-friendly user interface and source code for GoT-WAVE. Availability and implementation http://www.dcc.fc.up.pt/got-wave/ Supplementary information Supplementary data are available at Bioinformatics online.
Aparício, David↗
DOE BSSD Performance Management Metrics Report Q1
Microbes play key roles in our biosphere, from driving global nutrient cycling to impacting plant, animal and human health and disease. Complex data from microbial genomes, proteins, and metabolites provide a window into these tiny engines that drive life on our planet. Yet these data are dispersed among researchers’ laboratories and various repositories, making it difficult to access. This calls for new ways of managing data, improving data interoperability, advancing community standards, and creating an infrastructure where data are shared efficiently. We have built the National Microbiome Data Collaborative (NMDC) to advance how scientists create, use, and reuse data to redefine the way we understand and harness the power of microbes. The vision of the National Microbiome Data Collaborative (NMDC) is to drive a microbiome data sharing network connecting data, people, and ideas to advance microbiome innovation and discovery. The NMDC was launched in 2019 and brought together DOE National Laboratories to collaborate across resources, capabilities, and expertise. The NMDC team was strategically assembled to include software developers, microbial researchers, metadata experts, and multi-omics specialists. The diversity of the NMDC team reflects the inherently interdisciplinary nature of microbiome science, and we leverage the strengths of the DOE National Laboratory system. Towards BER’s goal of advancing an iterative systems biology approach to the understanding of microbial genomes, the NMDC serves as a foundation for infrastructure, data standards, and community building. Together with the flagship DOE User Facilities, the Joint Genome Institute (JGI) and the Environmental Molecular Sciences Laboratory (EMSL), we are developing core capabilities in metadata standards for environmental descriptors and sample handling and processing; standardized bioinformatic workflows; an interface for data search and access; and robust community engagement activities. The NMDC production platform supports long-term data infrastructure and community building for BER’s bioenergy and environmental research goals. Our approach leverages lessons learned and an ambitious framework for collaborative, interdisciplinary data infrastructure to support microbiome research. The NMDC supports data, information, and knowledge access through three defined software tools – the Submission Portal, NMDC EDGE, and the Data Portal – driven by community needs. Herein, we describe the value proposition for the microbiome research community, our overarching strategy, and challenges and opportunities for developing the NMDC as both an infrastructure and community engagement program.
59 BASIC BIOLOGICAL SCIENCES↗
EDGE COVID-19: a web platform to generate submission-ready genomes from SARS-CoV-2 sequencing efforts
Abstract Summary Genomics has become an essential technology for surveilling emerging infectious disease outbreaks. A range of technologies and strategies for pathogen genome enrichment and sequencing are being used by laboratories worldwide, together with different and sometimes ad hoc, analytical procedures for generating genome sequences. A fully integrated analytical process for raw sequence to consensus genome determination, suited to outbreaks such as the ongoing COVID-19 pandemic, is critical to provide a solid genomic basis for epidemiological analyses and well-informed decision making. We have developed a web-based platform and integrated bioinformatic workflows that help to provide consistent high-quality analysis of SARS-CoV-2 sequencing data generated with either the Illumina or Oxford Nanopore Technologies (ONT). Using an intuitive web-based interface, this workflow automates data quality control, SARS-CoV-2 reference-based genome variant and consensus calling, lineage determination and provides the ability to submit the consensus sequence and necessary metadata to GenBank, GISAID and INSDC raw data repositories. We tested workflow usability using real world data and validated the accuracy of variant and lineage analysis using several test datasets, and further performed detailed comparisons with results from the COVID-19 Galaxy Project workflow. Our analyses indicate that EC-19 workflows generate high-quality SARS-CoV-2 genomes. Finally, we share a perspective on patterns and impact observed with Illumina versus ONT technologies on workflow congruence and differences. Availability and implementation https://edge-covid19.edgebioinformatics.org, and https://github.com/LANL-Bioinformatics/EDGE/tree/SARS-CoV2. Supplementary information Supplementary data are available at Bioinformatics online.
59 BASIC BIOLOGICAL SCIENCES↗
Pathways to Carbon Liberation: A Systems Approach to Understanding Carbon Transformations and Losses from Thawing Permafrost
Our objective in this project was to discover how microbial communities mediate the fate of carbon in thawing permafrost under climate change. We proposed a systems approach integrating (a) molecular microbial and viral ecology, (b) molecular organic chemistry and stable and radiocarbon isotopes, and (c) state-of-the-art modeling, along a chronosequence of permafrost thaw in subarctic Sweden. The fate of carbon (C) in thawing permafrost is an outstanding challenge of modern biogeochemistry and climate change. Permafrost C pools are large (~1700 PgC), and C dynamics of thawing permafrost complex: old C decomposes as it is liberated from thawing permafrost as CO 2 or CH 4 , with a significant fraction cycling through lake sediments, even as new C accumulates due to thaw-initiated ecological succession. Our work is allowing better prediction of the net effect of these processes. Microbes mediate C dynamics in thawing permafrost, but a mechanistic understanding of how to scale microbial population dynamics, genomic potential, and expression to ecosystem-scale processes has been missing. A key question was: What is the interplay of microbial communities and organic matter chemical structure in the decomposition/preservation of organic C across a thaw gradient? And intriguingly, what (if any) is the role of phage (viruses that infect prokaryotic cells) in mediating these processes? Viruses appear to play a large role in driving oceanic function, but these phenomena are virtually unstudied in terrestrial systems. This endeavor linking microbial and viral dynamics, organic geochemistry and trace gas production will improve models of C cycling in permafrost systems, and clarify the fate of C under future climates. Technical Approach – We conducted a study at Stordalen Mire, Sweden, along a permafrost thaw chronosequence encompassing permafrost palsas and their initial collapse stage, thawing bog sites (dominated by Sphagnum moss), fully-thawed and inundated fen sites (dominated by Eriophorum spp or Carex spp.), and lakes. Our project used cutting edge technologies in both biogeochemistry and molecular microbial ecology to advance systems biology research on microbial carbon cycling through: (a) systems-level mapping of chemical states and ages of organic matter (via FT-ICR-MS and 14 C analysis) along thaw gradients to associated microbial communities, biochemical potentials, and activities (via meta-genomics, -transcriptomics, –proteomics, and viral genomics), and to CO 2 and CH 4 fluxes; (b) experimental incubations using Quantum Dot Probing, to test key hypotheses, arising from (a), about particular microbially-driven biochemical degradation pathways; (c) bioinformatics designed to simultaneously enhance the DOE Knowledgebase, and (d) integrated ecosystem C-cycle modeling testable by soil organic chemical and microbial data. Results from this project spanned 4 categories: (i) integrative across disciplines, and disciplinarily-focused in (ii) microbiology, (iii) biogeochemistry, and (iv) modeling. In integrated understanding of the system, we answered the questions posed above, mapping resident microbes to carbon transformations to reveal a carbon-processing shift with permafrost thaw, and discovering abundant, diverse and habitat-specific viruses that appear to be impacting carbon cycling, while also relating microbial communities to modeled processes. In microbiology-specific papers, we dug into viromes, microbial community networks, and informatic methods, while also introducing a new permafrost-associated phylum and its ecogenomic context. Relevance to DOE FOA - With a focus on discovering how microbial communities mediate the fate of old and new carbon in permafrost systems with climate change, our work directly responds to the call for “–omics driven basic research on the contribution of... microbial communities to C cycling processes in terrestrial ecosystems.” We bridge all three of the focal areas of the call, including (i) “Systems biology studies”, and development of –omics approaches (ii) “to investigate microbial community functional processes”, and (iii) “for imaging and analysis of microbially-mediated carbon cycling.”
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
From Structure to Functions: Metagenomics-Enabled Predictive Understanding of Soil Microbial Feedbacks to Climate Change
Understanding the responses, adaptations and feedback mechanisms of biological communities to climate change is critical to project future states of earth and climate systems. Although a significant amount of knowledge is available on the feedbacks of aboveground communities to climate change, the responses of belowground microbial communities are still poorly understood. This is due to challenges in analyzing soil microbial communities in terms of their diversity, structure and function. Thus, the ultimate goal of this study is to advance system-level predictive understanding of the feedbacks of belowground microbial communities to multiple climate change factors and their impacts on soil carbon (C) cycling processes. For this, cutting-edge integrated metagenomic technologies are employed to establish mechanistic linkages among aboveground plant communities, belowground microbial communities, ecosystem processes and functions and climate change. Towards this goal, the following four specific objectives are pursued: (i) To determine the responses of microbial community structure, functions and activities to climate warming, altered precipitation, soil moisture regime and/or clipping in the tundra and temperate grassland ecosystems; (ii) To determine the temperature sensitivity and the effect of substrate priming on recalcitrant C decomposition; (iii) To determine microbiological bases underlying temperature sensitivity of recalcitrant C decomposition; and (iv) To develop integrated bioinformatics and modeling approaches to scale information across different organizational levels towards predictive understanding of ecosystem responses to multiple climate change factors in collaboration with the Knowledge Base (KBase) and integrated into their database. Warming experiments in an Alaskan tundra permafrost soil (AK) as well as an Oklahoma temperate grassland soil (OK) provide a more complete picture of the effects of warming on belowground microbial communities of different (key) soil ecosystems. Numerous research activities have been carried out by the institutions under the support of this program. For brevity, a few key results and their significance are highlighted in this report.
54 ENVIRONMENTAL SCIENCES↗
Autotrophic biofilms sustained by deeply sourced groundwater host diverse bacteria implicated in sulfur and hydrogen metabolism
Abstract Background Biofilms in sulfide-rich springs present intricate microbial communities that play pivotal roles in biogeochemical cycling. We studied chemoautotrophically based biofilms that host diverse CPR bacteria and grow in sulfide-rich springs to investigate microbial controls on biogeochemical cycling. Results Sulfide springs biofilms were investigated using bulk geochemical analysis, genome-resolved metagenomics, and scanning transmission X-ray microscopy (STXM) at room temperature and 87 K. Chemolithotrophic sulfur-oxidizing bacteria, including Thiothrix and Beggiatoa , dominate the biofilms, which also contain CPR Gracilibacteria, Absconditabacteria, Saccharibacteria, Peregrinibacteria, Berkelbacteria, Microgenomates, and Parcubacteria. STXM imaging revealed ultra-small cells near the surfaces of filamentous bacteria that may be CPR bacterial episymbionts. STXM and NEXAFS spectroscopy at carbon K and sulfur L 2,3 edges show that filamentous bacteria contain protein-encapsulated spherical elemental sulfur granules, indicating that they are sulfur oxidizers, likely Thiothrix . Berkelbacteria and Moranbacteria in the same biofilm sample are predicted to have a novel electron bifurcating group 3b [NiFe]-hydrogenase, putatively a sulfhydrogenase, potentially linked to sulfur metabolism via redox cofactors. This complex could potentially contribute to symbioses, for example, with sulfur-oxidizing bacteria such as Thiothrix that is based on cryptic sulfur cycling. One Doudnabacteria genome encodes adjacent sulfur dioxygenase and rhodanese genes that may convert thiosulfate to sulfite. We find similar conserved genomic architecture associated with CPR bacteria from other sulfur-rich subsurface ecosystems. Conclusions Our combined metagenomic, geochemical, spectromicroscopic, and structural bioinformatics analyses of biofilms growing in sulfide-rich springs revealed consortia that contain CPR bacteria and sulfur-oxidizing Proteobacteria, including Thiothrix , and bacteria from a new family within Beggiatoales. We infer roles for CPR bacteria in sulfur and hydrogen cycling.
59 BASIC BIOLOGICAL SCIENCES↗