SEARCH · Search NASA

Results for “genomic”

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 541 records · Page 30

All sources NASA NTRS DOE OSTI

DOE OSTI2023

BayFlux: A Bayesian method to quantify metabolic Fluxes and their uncertainty at the genome scale

Metabolic fluxes, the number of metabolites traversing each biochemical reaction in a cell per unit time, are crucial for assessing and understanding cell function. 13 C Metabolic Flux Analysis ( 13 C MFA) is considered to be the gold standard for measuring metabolic fluxes. 13 C MFA typically works by leveraging extracellular exchange fluxes as well as data from 13 C labeling experiments to calculate the flux profile which best fit the data for a small, central carbon, metabolic model. However, the nonlinear nature of the 13 C MFA fitting procedure means that several flux profiles fit the experimental data within the experimental error, and traditional optimization methods offer only a partial or skewed picture, especially in “non-gaussian” situations where multiple very distinct flux regions fit the data equally well. Here, we present a method for flux space sampling through Bayesian inference (BayFlux), that identifies the full distribution of fluxes compatible with experimental data for a comprehensive genome-scale model. This Bayesian approach allows us to accurately quantify uncertainty in calculated fluxes. We also find that, surprisingly, the genome-scale model of metabolism produces narrower flux distributions (reduced uncertainty) than the small core metabolic models traditionally used in 13 C MFA. The different results for some reactions when using genome-scale models vs core metabolic models advise caution in assuming strong inferences from 13 C MFA since the results may depend significantly on the completeness of the model used. Based on BayFlux, we developed and evaluated novel methods (P- 13 C MOMA and P- 13 C ROOM) to predict the biological results of a gene knockout, that improve on the traditional MOMA and ROOM methods by quantifying prediction uncertainty.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Genomes to Structure and Function Workshop Report 2022

The goal of the U.S. Department of Energy (DOE) Biological and Environmental Research (BER) Program is to achieve a predictive understanding of complex biological, earth, and environmental systems with the aim of advancing the nation’s energy and infrastructure security. (https://www.energy.gov/science/ ber/biological-and-environmental-research). To pursue this goal, collaborations among experts in diverse research areas that lead to multidisciplinary projects are indispensable. The roles of DOE’s User Facilities, which offer unique and powerful resources for such research projects, are evolving, and expectations for the facilities are increasing. To respond to Users’ needs, the Joint Genome Institute (JGI) and Environmental Molecular Sciences Laboratory (EMSL) initiated the Facilities Integrating Collaborations for User Science (FICUS) program in 2014. This collaboration has grown into a popular and successful program, advancing more than 100 multidisciplinary projects to date. Similarly, the new interFacility collaborations among the JGI, EMSL, and User resources for BER structural biology and imaging at the Basic Energy Science (BES) Program’s synchrotron and neutron facilities are becoming essential for cutting-edge transdisciplinary science. To further explore the need for the BER research community to combine genomic, functional, and structural approaches to advance their research, an organizing committee was formed to develop and jointly host a 3-part workshop. The committee’s members represented seven DOE National Laboratory User Facilities (Appendix 1 lists the members). The “Genomes to Structure and Function” virtual workshop (see Appendices 2–5) was composed of three sessions. The first session, titled “Molecular Structures” (October 27– 28, 2021), highlighted diverse integrative experimental and computational approaches correlating structural data with sequencing and functional information, as well as predicting protein structures to model complex biological systems. The second session, “Intracellular Organization, and Material Synthesis and Decomposition” (December 15–16, 2021), covered imaging methods for observing, quantifying, and manipulating biosystems. The third session, “Imaging the Rhizosphere and Cellular Organization” (January 26–27, 2022) emphasized advanced and non-invasive imaging techniques applied to plant root-microbe-soil interactions.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Genomes_SET2_178 - Shared

The dataset consists of 178 bacterial genomes in total. Of these, 145 genomes were drawn from published literature linked to Biolog Gen III growth data, while 33 genomes came from soil isolates collected at the Cedar Creek Ecosystem Science Research Reserve (CCESR) in Minnesota. The isolates were sequenced on an Illumina HiSeq platform (~50× coverage), assembled with Megahit, and annotated using RAST at PATRIC, with strict quality filtering (<90% completeness or >5% contamination).

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Leisingera sp. ANG59 & Pseudoalteromonas sp. JC28 Genomes Overview

The Hawaiian bobtail squid, Euprymna scolopes, has a symbiotic bacterial consortium in the accessory nidamental gland (ANG), a female reproductive organ that protects eggs against fouling microorganisms. To test the antibacterial activity of ANG community members, 19 bacterial isolates were screened for their ability to inhibit Gram-negative and Gram-positive bacteria, of which two strains were inhibitory. These two antibacterial isolates, Leisingera sp. ANG59 and Pseudoalteromonas sp. JC28, were subjected to further genomic characterization. Genomic analysis of Leisingera sp. ANG59 revealed a biosynthetic gene cluster encoding the antimicrobial compound indigoidine. The genome of Pseudoalteromonas sp. JC28 had a 14-gene cluster with >95% amino acid identity to a known bromoalterochromide (BAC) cluster. Chemical analysis confirmed production of known BACs, BAC-A/A′ (compounds 1a/1b), as well as two new derivatives, BAC-D/D′ (compounds 2a/2b). Extensive nuclear magnetic resonance (NMR) analyses allowed complete structural elucidation of compounds 2a/2b, and the absolute stereochemistry was unambiguously determined using an optimized Marfey’s method. The BACs were then investigated for in vitro antibacterial, antifungal, and nitric oxide (NO) inhibitory activity. Compounds 1a/1b were active against the marine bacteria Bacillus algicola and Vibrio fischeri, while compounds 2a/2b were active only against B. algicola. Compounds 1a/1b inhibited NO production via lipopolysaccharide (LPS)-induced inflammation in RAW264.7 macrophage cells and also inhibited the pathogenic fungus Fusarium keratoplasticum, which, coupled with their antibacterial activity, suggests that these polyketide-nonribosomal peptides may be used for squid egg defense against potential pathogens and/or fouling microorganisms. These results indicate that BACs may provide Pseudoalteromonas sp. JC28 an ecological niche, facilitating competition against nonsymbiotic microorganisms in the host’s environment.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Whole-Genome Resequencing to Evaluate Life History Variation in Anadromous Migration of Oncorhynchus mykiss

Anadromous fish experience physiological modifications necessary to migrate between vastly different freshwater and marine environments, but some species such as Oncorhynchus mykiss demonstrate variation in life history strategies with some individuals remaining exclusively resident in freshwater, whereas others undergo anadromous migration. Because there is limited understanding of genes involved in this life history variation across populations of this species, we evaluated the genomic difference between known anadromous ( n = 39) and resident ( n = 78) Oncorhynchus mykiss collected from the Klickitat River, WA, USA, with whole-genome resequencing methods. Sequencing of these collections yielded 5.64 million single-nucleotide polymorphisms that were tested for significant differences between resident and anadromous groups along with previously identified candidate gene regions. Although a few regions of the genome were marginally significant, there was one region on chromosome Omy12 that provided the most consistent signal of association with anadromy near two annotated genes in the reference assembly: COP9 signalosome complex subunit 6 (CSN6) and NACHT, LRR, and PYD domain–containing protein 3 (NLRP3). Previously identified candidate genes for anadromy within the inversion region of chromosome Omy05 in coastal steelhead and rainbow trout were not informative for this population as shown in previous studies. Results indicate that the significant region on chromosome Omy12 may represent a minor effect gene for male anadromy and suggests that this life history variation in Oncorhynchus mykiss is more strongly driven by other mechanisms related to environmental rearing such as epigenetic modification, gene expression, and phenotypic plasticity. Further studies into regulatory mechanisms of this trait are needed to understand drivers of anadromy in populations of this protected species.

Collins, Erin E.↗

DOE OSTI2022

Antimicrobial resistance and genomic characterization of Escherichia coli from pigs and chickens in Zhejiang, China

Escherichia coli is considered an opportunistic pathogen and an indicator for antimicrobial resistance (AMR) monitoring. Despite many reports on its AMR monitoring, studies based on genome-based analysis of AMR genes are still insufficient. Here, 181 E. coli strains were isolated from anal swab samples collected from pigs and chickens of animal farms located in Eastern China and sequenced through the Illumina platform. The results showed that 87.85% (159/181) of the E. coli isolates were multidrug-resistant (MDR). Ampicillin (AMP)- spectinomycin (SPT)- tetracycline (TET)- florfenicol (FFC)- sulfisoxazole (SF)- trimethoprim/sulfamethoxazole (SXT) was the predominant AMR pattern. By whole-genome sequencing, we found that ST10 (10.49%, 19/181) and ST48 (7.18%, 13/181) were major sequence types. IncFIB and IncX1 were the most prevalent plasmid replicons. The AMR genes bla NDM-5 (1.10%, 2/181), mcr -1 (1.10%, 2/181), tet (X4) (1.10%, 2/181), and cfr (6.08%, 11/181) were also found in these isolates. In addition, among the 169 virulence genes detected, we identified ast A (37.02%, 67/181), hly A (1.66%, 3/181), hly B (1.66%, 3/181) and hly D (1.66%, 3/181), which were closely related to heat-stable enterotoxin 1 and α-hemolysin. In addition, there were 33 virulence genes associated with the iron uptake system, and 46 were adhesion-related genes. Our study highlighted the need for routine surveillance of AMR with advanced genomic approaches, providing up-to-date data on the prevalence of AMR for the development and execution of antimicrobial stewardship policy.

Zhou, Wei↗

DOE OSTI2024

Metagenome-assembled genomes provide insight into the metabolic potential during early production of Hydraulic Fracturing Test Site 2 in the Delaware Basin

Demand for natural gas continues to climb in the United States, having reached a record monthly high of 104.9 billion cubic feet per day (Bcf/d) in November 2023. Hydraulic fracturing, a technique used to extract natural gas and oil from deep underground reservoirs, involves injecting large volumes of fluid, proppant, and chemical additives into shale units. This is followed by a “shut-in” period, during which the fracture fluid remains pressurized in the well for several weeks. The microbial processes that occur within the reservoir during this shut-in period are not well understood; yet, these reactions may significantly impact the structural integrity and overall recovery of oil and gas from the well. To shed light on this critical phase, we conducted an analysis of both pre-shut-in material alongside production fluid collected throughout the initial production phase at the Hydraulic Fracturing Test Site 2 (HFTS 2) located in the prolific Wolfcamp formation within the Permian Delaware Basin of west Texas, USA. Specifically, we aimed to assess the microbial ecology and functional potential of the microbial community during this crucial time frame. Prior analysis of 16S rRNA sequencing data through the first 35 days of production revealed a strong selection for a Clostridia species corresponding to a significant decrease in microbial diversity. Here, we performed a metagenomic analysis of produced water sampled on Day 33 of production. This analysis yielded three high-quality metagenome-assembled genomes (MAGs), one of which was a Clostridia draft genome closely related to the recently classified Petromonas tenebris. This draft genome likely represents the dominant Clostridia species observed in our 16S rRNA profile. Annotation of the MAGs revealed the presence of genes involved in critical metabolic processes, including thiosulfate reduction, mixed acid fermentation, and biofilm formation. These findings suggest that this microbial community has the potential to contribute to well souring, biocorrosion, and biofouling within the reservoir. Our research provides unique insights into the early stages of production in one of the most prolific unconventional plays in the United States, with important implications for well management and energy recovery.

natural gas↗

DOE OSTI2021

Multiple Wheat Genomes Reveal Novel Gli-2 Sublocus Location and Variation of Celiac Disease Epitopes in Duplicated α-Gliadin Genes

The seed protein α-gliadin is a major component of wheat flour and causes gluten-related diseases. However, due to the complexity of this multigene family with a genome structure composed of dozens of copies derived from tandem and genome duplications, little was known about the variation between accessions, and thus little effort has been made to explicitly target α-gliadin for bread wheat breeding. Here, we analyzed genomic variation in α-gliadins across 11 recently published chromosome-scale assemblies of hexaploid wheat, with validation using long-read data. We unexpectedly found that the Gli-B2 locus is not a single contiguous locus but is composed of two subloci, suggesting the possibility of recombination between the two during breeding. We confirmed that the number of immunogenic epitopes among 11 accessions varied. The D subgenome of a European spelt line also contained epitopes, in agreement with its hybridization history. Evolutionary analysis identified amino acid sites under diversifying selection, suggesting their functional importance. The analysis opens the way for improved grain quality and safety through wheat breeding.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Genomic Analysis of the Natural Variation of Fatty Acid Composition in Seed Oils of Camelina sativa

Camelina sativa is an oilseed crop that has shown strong promise as a biofuel feedstock. The profile of fatty acids greatly influences the oil quality; however, genetic mechanisms that determine the natural variation of fatty acid composition in camelina are not fully understood. A genome wide association study (GWAS) was performed to uncover genetic loci that may contribute to the contents of major fatty acids such as oleic and linolenic acids in camelina seed. Two approaches were taken to improve the GWAS efficiency. First, growing a diversity panel of 212 accessions in four locations and two nitrogen fertilization conditions revealed great variation in fatty acid contents in seeds. Second, using an improved reference genome, abundant markers, including 203,320 single nucleotide polymorphisms (SNPs) and 99,067 insertions/deletions (indels), were developed, which refined the population structure of the diversity panel. GWAS resulted in 118 genetic markers across 31 trait/treatment conditions. Closely linked markers were determined based on linkage decay and by comparing secondarily associated markers when highly associated ones were removed. Candidate genes were examined by comparing the pangenomes of 12 high-quality reference genomes. This study provides new resources to understand seed lipid metabolism and improve camelina oils through molecular breeding.

Life Sciences & Biomedicine - Other Topics↗

DOE OSTI2022

A Comparative Analysis of Weizmannia coagulans Genomes Unravels the Genetic Potential for Biotechnological Applications

The production of biochemicals requires the use of microbial strains with efficient substrate conversion and excellent environmental robustness, such as Weizmannia coagulans species. So far, the genomes of 47 strains have been sequenced. Herein, we report a comparative genomic analysis of nine strains on the full repertoire of Carbohydrate-Active enZymes (CAZymes), secretion systems, and resistance mechanisms to environmental challenges. Moreover, Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) immune system along with CRISPR-associated (Cas) genes, was also analyzed. Overall, this study expands our understanding of the strain’s genomic diversity of W. coagulans to fully exploit its potential in biotechnological applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Genome Sequence and Analysis of the Flavinogenic Yeast Candida membranifaciens IST 626

The ascomycetous yeast Candida membranifaciens has been isolated from diverse habitats, including humans, insects, and environmental sources, exhibiting a remarkable ability to use different carbon sources that include pentoses, melibiose, and inulin. In this study, we isolated four C. membranifaciens strains from soil and investigated their potential to overproduce riboflavin. C. membranifaciens IST 626 was found to produce the highest concentrations of riboflavin. The volumetric production of this vitamin was higher when C. membranifaciens IST 626 cells were cultured in a commercial medium without iron and when xylose was the available carbon source compared to the same basal medium with glucose. Supplementation of the growth medium with 2 g/L glycine favored the metabolization of xylose, leading to biomass increase and consequent enhancement of riboflavin volumetric production that reached 120 mg/L after 216 h of cultivation. To gain new insights into the molecular basis of riboflavin production and carbon source utilization in this species, the first annotated genome sequence of C. membranifaciens is reported in this article, as well as the result of a comparative genomic analysis with other relevant yeast species. A total of 5619 genes were predicted to be present in C. membranifaciens IST 626 genome sequence (11.5 Mbp). Among them are genes involved in riboflavin biosynthesis, iron homeostasis, and sugar uptake and metabolism. This work put forward C. membranifaciens IST 626 as a riboflavin overproducer and provides valuable molecular data for future development of superior producing strains capable of using the wide range of carbon sources, which is a characteristic trait of the species.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Genomic Insights into Moderately Thermophilic Methanotrophs of the Genus Methylocaldum

Considering the increasing interest in understanding the biotic component of methane removal from our atmosphere, it becomes essential to study the physiological characteristics and genomic potential of methanotroph isolates, especially their traits allowing them to adapt to elevated growth temperatures. The genetic signatures of Methylocaldum species have been detected in many terrestrial and aquatic ecosystems. A small set of representatives of this genus has been isolated and maintained in culture. The genus is commonly described as moderately thermophilic, with the growth optimum reaching 50 °C for some strains. Here, we present a comparative analysis of genomes of three Methylocaldum strains—two terrestrial M. szegediense strains (O-12 and Norfolk) and one marine strain, Methylocaldum marinum (S8). The examination of the core genome inventory of this genus uncovers significant redundancy in primary metabolic pathways, including the machinery for methane oxidation (numerous copies of pmo genes) and methanol oxidation (duplications of mxaF, xoxF1-5 genes), three pathways for one-carbon (C1) assimilation, and two methods of carbon storage (glycogen and polyhydroxyalkanoates). We also investigate the genetics of melanin production pathways as a key feature of the genus.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Comparative Genomic Analysis of Three Pseudomonas Species Isolated from the Eastern Oyster (Crassostrea virginica) Tissues, Mantle Fluid, and the Overlying Estuarine Water Column

The eastern oysters serve as important keystone species in the United States, especially in the Gulf of Mexico estuarine waters, and at the same time, provide unparalleled economic, ecological, environmental, and cultural services. One ecosystem service that has garnered recent attention is the ability of oysters to sequester impurities and nutrients, such as nitrogen (N), from the estuarine water that feeds them, via their exceptional filtration mechanism coupled with microbially-mediated denitrification processes. It is the oyster-associated microbiomes that essentially provide these myriads of ecological functions, yet not much is known on these microbiota at the genomic scale, especially from warm temperate and tropical water habitats. Among the suite of bacterial genera that appear to interplay with the oyster host species, pseudomonads deserve further assessment because of their immense metabolic and ecological potential. To obtain a comprehensive understanding on this aspect, we previously reported on the isolation and preliminary genomic characterization of three Pseudomonas species isolated from minced oyster tissue (P. alcaligenes strain OT69); oyster mantle fluid (P. stutzeri strain MF28) and the water collected from top of the oyster reef (P. aeruginosa strain WC55), respectively. In this comparative genomic analysis study conducted on these three targeted pseudomonads, native to the eastern oyster and its surrounding environment, provided further insights into their unique functional traits, conserved gene pools between the selected pseudomonads, as well as genes that render unique characteristics in context to metabolic traits recruited during their evolutionary history via horizontal gene transfer events as well as phage-mediated incorporation of genes. Moreover, the strains also supported extensively developed resistomes, which suggests that environmental microorganisms native to relatively pristine environments, such as Apalachicola Bay, Florida, have also recruited an arsenal of antibiotic resistant gene determinants, thus posing an emerging public health concern.

Pathak, Ashish↗

DOE OSTI2025

Status on Genetic Resistance to Rice Blast Disease in the Post-Genomic Era

Rice blast, caused by Magnaporthe oryzae, is a major threat to global rice production, necessitating the development of resistant cultivars through genetic improvement. Breakthroughs in rice genomics, including the complete genome sequencing of japonica and indica subspecies and the availability of various sequence-based molecular markers, have greatly advanced the genetic analysis of blast resistance. To date, approximately 122 blast-resistance genes have been identified, with 39 of these genes cloned and molecularly characterized. The application of these findings in marker-assisted selection (MAS) has significantly improved rice breeding, allowing for the efficient integration of multiple resistance genes into elite cultivars, enhancing both the durability and spectrum of resistance. Pangenomic studies, along with AI-driven tools like AlphaFold2, RoseTTAFold, and AlphaFold3, have further accelerated the identification and functional characterization of resistance genes, expediting the breeding process. Future rice blast disease management will depend on leveraging these advanced genomic and computational technologies. Emphasis should be placed on enhancing computational tools for the large-scale screening of resistance genes and utilizing gene editing technologies such as CRISPR-Cas9 for functional validation and targeted resistance enhancement and deployment. These approaches will be crucial for advancing rice blast resistance, ensuring food security, and promoting agricultural sustainability.

Pedrozo, Rodrigo↗

DOE OSTI2024

Recommendations for Uniform Variant Calling of SARS-CoV-2 Genome Sequence across Bioinformatic Workflows

Genomic sequencing of clinical samples to identify emerging variants of SARS-CoV-2 has been a key public health tool for curbing the spread of the virus. As a result, an unprecedented number of SARS-CoV-2 genomes were sequenced during the COVID-19 pandemic, which allowed for rapid identification of genetic variants, enabling the timely design and testing of therapies and deployment of new vaccine formulations to combat the new variants. However, despite the technological advances of deep sequencing, the analysis of the raw sequence data generated globally is neither standardized nor consistent, leading to vastly disparate sequences that may impact identification of variants. Here, we show that for both Illumina and Oxford Nanopore sequencing platforms, downstream bioinformatic protocols used by industry, government, and academic groups resulted in different virus sequences from same sample. These bioinformatic workflows produced consensus genomes with differences in single nucleotide polymorphisms, inclusion and exclusion of insertions, and/or deletions, despite using the same raw sequence as input datasets. Here, we compared and characterized such discrepancies and propose a specific suite of parameters and protocols that should be adopted across the field. Consistent results from bioinformatic workflows are fundamental to SARS-CoV-2 and future pathogen surveillance efforts, including pandemic preparation, to allow for a data-driven and timely public health response.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Genomic comparisons of a bacterial lineage that inhabits both marine and terrestrial deep subsurface systems

It is generally accepted that diverse, poorly characterized microorganisms reside deep within Earth’s crust. One such lineage of deep subsurface-dwelling Bacteria is an uncultivated member of the Firmicutes phylum that can dominate molecular surveys from both marine and continental rock fracture fluids, sometimes forming the sole member of a single-species microbiome. Here, we reconstructed a genome from basalt-hosted fluids of the deep subseafloor along the eastern Juan de Fuca Ridge flank and used a phylogenomic analysis to show that, despite vast differences in geographic origin and habitat, it forms a monophyletic clade with the terrestrial deep subsurface genome of “ Candidatus Desulforudis audaxviator” MP104C. While a limited number of differences were observed between the marine genome of “ Candidatus Desulfopertinax cowenii” modA32 and its terrestrial relative that may be of potential adaptive importance, here it is revealed that the two are remarkably similar thermophiles possessing the genetic capacity for motility, sporulation, hydrogenotrophy, chemoorganotrophy, dissimilatory sulfate reduction, and the ability to fix inorganic carbon via the Wood-Ljungdahl pathway for chemoautotrophic growth. Our results provide insights into the genetic repertoire within marine and terrestrial members of a bacterial lineage that is widespread in the global deep subsurface biosphere, and provides a natural means to investigate adaptations specific to these two environments.

Jungbluth, Sean P.↗

DOE OSTI2021

VirION2: a short- and long-read sequencing and informatics workflow to study the genomic diversity of viruses in nature

Microbes play fundamental roles in shaping natural ecosystem properties and functions, but do so under constraints imposed by their viral predators. However, studying viruses in nature can be challenging due to low biomass and the lack of universal gene markers. Though metagenomic short-read sequencing has greatly improved our virus ecology toolkit—and revealed many critical ecosystem roles for viruses—microdiverse populations and fine-scale genomic traits are missed. Some of these microdiverse populations are abundant and the missed regions may be of interest for identifying selection pressures that underpin evolutionary constraints associated with hosts and environments. Though long-read sequencing promises complete virus genomes on single reads, it currently suffers from high DNA requirements and sequencing errors that limit accurate gene prediction. Here we introduce VirION2, an integrated short- and long-read metagenomic wet-lab and informatics pipeline that updates our previous method (VirION) to further enhance the utility of long-read viral metagenomics. Using a viral mock community, we first optimized laboratory protocols (polymerase choice, DNA shearing size, PCR cycling) to enable 76% longer reads (now median length of 6,965 bp) from 100-fold less input DNA (now 1 nanogram). Using a virome from a natural seawater sample, we compared viromes generated with VirION2 against other library preparation options (unamplified, original VirION, and short-read), and optimized downstream informatics for improved long-read error correction and assembly. VirION2 assemblies combined with short-read based data (‘enhanced’ viromes), provided significant improvements over VirION libraries in the recovery of longer and more complete viral genomes, and our optimized error-correction strategy using long- and short-read data achieved 99.97% accuracy. In the seawater virome, VirION2 assemblies captured 5,161 viral populations (including all of the virus populations observed in the other assemblies), 30% of which were uniquely assembled through inclusion of long-reads, and 22% of the top 10% most abundant virus populations derived from assembly of long-reads. Viral populations unique to VirION2 assemblies had significantly higher microdiversity means, which may explain why short-read virome approaches failed to capture them. These findings suggest the VirION2 sample prep and workflow can help researchers better investigate the virosphere, even from challenging low-biomass samples. Our new protocols are available to the research community on protocols.io as a ‘living document’ to facilitate dissemination of updates to keep pace with the rapid evolution of long-read sequencing technology.

Long-reads↗

DOE OSTI2020

Deconvolute individual genomes from metagenome sequences through short read clustering

Metagenome assembly from short next-generation sequencing data is a challenging process due to its large scale and computational complexity. Clustering short reads by species before assembly offers a unique opportunity for parallel downstream assembly of genomes with individualized optimization. However, current read clustering methods suffer either false negative (under-clustering) or false positive (over-clustering) problems. Here we extended our previous read clustering software, SpaRC, by exploiting statistics derived from multiple samples in a dataset to reduce the under-clustering problem. Using synthetic and real-world datasets we demonstrated that this method has the potential to cluster almost all of the short reads from genomes with sufficient sequencing coverage. The improved read clustering in turn leads to improved downstream genome assembly quality.

59 BASIC BIOLOGICAL SCIENCES↗