Search NASA⌕ Search

SEARCH · Search NASA

Results for “genomic analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Evaluating E. coli genome‐scale metabolic model accuracy with high‐throughput mutant fitness data

Abstract The Escherichia coli genome‐scale metabolic model (GEM) is an exemplar systems biology model for the simulation of cellular metabolism. Experimental validation of model predictions is essential to pinpoint uncertainty and ensure continued development of accurate models. Here, we quantified the accuracy of four subsequent E. coli GEMs using published mutant fitness data across thousands of genes and 25 different carbon sources. This evaluation demonstrated the utility of the area under a precision–recall curve relative to alternative accuracy metrics. An analysis of errors in the latest (iML1515) model identified several vitamins/cofactors that are likely available to mutants despite being absent from the experimental growth medium and highlighted isoenzyme gene‐protein‐reaction mapping as a key source of inaccurate predictions. A machine learning approach further identified metabolic fluxes through hydrogen ion exchange and specific central metabolism branch points as important determinants of model accuracy. This work outlines improved practices for the assessment of GEM accuracy with high‐throughput mutant fitness data and highlights promising areas for future model refinement in E. coli and beyond.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

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↗
DOE OSTI2022

The IMG/M data management and analysis system v.7: content updates and new features

The Integrated Microbial Genomes & Microbiomes system at the Department of Energy (DOE) Joint Genome Institute (JGI) continues to provide support for users to perform comparative analysis of isolate and single cell genomes, metagenomes, and metatranscriptomes. In addition to datasets produced by the JGI, IMG v.7 also includes datasets imported from public sources such as NCBI Genbank, SRA, and the DOE National Microbiome Data Collaborative (NMDC), or submitted by external users. In the past couple years, we have continued our effort to help the user community by improving the annotation pipeline, upgrading the contents with new reference database versions, and adding new analysis functionalities such as advanced scaffold search, Average Nucleotide Identity (ANI) for high-quality metagenome bins, new cassette search, improved gene neighborhood display, and improvements to metatranscriptome data display and analysis. Here, we also extended the collaboration and integration efforts with other DOE-funded projects such as NMDC and DOE Biology Knowledgebase (KBase).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Nitrogen fixation in Pampinifervens , a new species-rich genus of Aquificaceae that inhabits a wide pH range in terrestrial hot springs

Members of the Aquificota play important roles in hydrogen, sulfur, and carbon cycles in geothermal systems as ubiquitous and abundant chemolithoautotrophs; however, their roles in the nitrogen cycle are poorly defined. Here, we show that Aquificaceae isolate T-2 T from Tengchong, China, fixes nitrogen based on growth with dinitrogen as the sole nitrogen source and incorporation of 15 N atoms when grown with 15 N 2 . We further show evidence that suggests its nitrogen fixation (nif) genes are transcribed in situ in a nearby spring. We propose the name Pampinifervens diazotrophicum gen. et sp. nov., with strain T-2 T (= JCM 35475 T = DSM 116324 T ) as the nomenclatural type. Its close relative designated Pampinifervens florentissimum sp. nov. T-8 T (= CGMC 1.5214 T = JCM 33569 T ) does not fix nitrogen but respires nitrate. Analysis of 104 metagenome-assembled genomes from 84 metagenomes allowed us to circumscribe the genus Pampinifervens and identify at least 16 sympatric Pampinifervens species that dominate the Aquificaceae populations in diverse hot springs in southwest China, with distinct species globally. Analysis of these Pampinifervens genomes revealed variable presence of nif and denitrification genes and frequent gene loss and horizontal transfer. Genes for pH homeostasis were also present in Pampinifervens genomes, with a proliferation in species from Tengchong, consistent with the wide pH growth range of both Pampinifervens isolates and the prevalence of Pampinifervens in both alkaline and acidic springs. The dynamic evolution of nif and denitrification genes, along with pH homeostasis genes, may facilitate the diversification of Pampinifervens into diverse springs with different nitrogen availability.

Aquificae↗
DOE OSTI2020

Molecular subtyping reveals immune alterations in IDH wild‐type lower‐grade diffuse glioma

Abstract Isocitrate dehydrogenase ( IDH ) wild‐type diffuse lower‐grade glioma (LGG) is usually associated with poor outcome, but there have been disputes over its clinical outcome and classification. We present here a robust gene expression‐based molecular classification of IDH wild‐type diffuse LGG into two subtypes with distinct biological and clinical features. A discovery cohort of 49 IDH wild‐type diffuse LGGs from the Chinese Glioma Genome Atlas (CGGA) was subjected to clustering and function analysis. Seventy‐three tumors from The Cancer Genome Atlas (TCGA) were used to validate our findings. Consensus clustering of transcriptional data uncovered concordant classification of two robust and prognostically significant subtypes of IDH wild‐type LGG. Subtype 1, associated with poorer outcomes, was characterized by significantly higher immune and cytolytic scores, M2 macrophages, and up‐regulation of immune exhaustion markers, while Subtype 2, which had elevated lymphocytes and plasma cells, showed relatively favorable survival. Somatic alteration analysis revealed that Subtype 1 showed more frequently deleted regions, such as the locus of CDKN2A / CDKN2B , DMRTA1 , C9orf53 , and MTAP . Furthermore, we developed and validated a five‐gene signature for better application of this acquired stratification. Our data demonstrate the biological and prognostic heterogeneity within IDH wild‐type diffuse LGGs and deepen our molecular understandi‐g of this tumor entity. © 2020 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Wu, Fan↗
NASA NTRS2018

Variability in Galactic Cosmic Radiation- Induced DNA Damage Response in Inbred Mice Is Modulated by Genetics

In radiation biology, the ability to predict cancer risk associated with exposure to low doses of high-LET (Linear Energy Transfer) ionizing radiation remains a challenge. Epidemiological methods lack the sensitivity and power to provide detailed risk estimates for cancer and ignore individual sensitivity. We have hypothesized that DNA repair capacity is the primary factor differentiating peoples radiation sensitivity. We previously showed in immortalized human cell lines that characterizing the dose and time dependence of p53-binding protein 1 (53BP1) foci formation in the nucleus following X-rays exposure is sufficient to predict DNA repair response to any other LET in the same cell line. We now tested this hypothesis across a population of mice with different genetic background. Fibroblast cells were extracted and cultivated from 76 individual mice from 15 different strains and exposed to HZE (high (H) atomic number (Z) and energy (E) galactic cosmic ray particles) particles and X-rays. Individual radiation sensitivities were investigated by high throughput measurement of DNA repair kinetics that evaluated 53bp1 foci numbers as a surrogate for DNA double-strand breaks at various times post-irradiation. Instead of just counting foci which can be hard to distinguish for high-LET or high doses, we also took into account the track structure of high-LET particles to compute the remaining number of unrepaired tracks as a function of time post-irradiation. As expected, the percentage of unrepaired track over a 48 hours follow-up period increased with LET. In addition, repair rate was modulated by genetics, with animals from the same strain showing small variance while large rate differences were observed between strains. Radiation strain sensitivity ranking was estimated based on repair rates from exposure to each LET evaluated in this work, and ranking for high-LET correlated better with ranking from high dose of X-ray, not low dose. At the in-vivo level, drops in T-cells and B-cells number measured 24 hours after 0.1 Gy (Gray) X-ray exposure, correlated with slower DNA repair kinetic in fibroblast cells of the same strains of mice. At the genomic level, mouse genome wide association (GWA) analysis identified seven significant genetic loci on chromosomes 2, 3, 7, 10, 11, 13 and 19 with different significance depending on the LET. Interestingly, for the two highest LET, a common locus on Chromosome 10 was identified with high enrichment for DNA repair associated genes.Overall, this work suggests that repair kinetics of primary skin fibroblasts is a good surrogate marker for in-vivo radiation sensitivities in other tissues and that this response is modulated by genetics. Our study also confirms that DNA repair kinetics following high doses of X-ray can be used to predict radiation sensitivity to high-LET.

High-LET Radiation↗
NASA NTRS2003

An in silico assessment of gene function and organization of the phenylpropanoid pathway metabolic networks in Arabidopsis thaliana and limitations thereof

The Arabidopsis genome sequencing in 2000 gave to science the first blueprint of a vascular plant. Its successful completion also prompted the US National Science Foundation to launch the Arabidopsis 2010 initiative, the goal of which is to identify the function of each gene by 2010. In this study, an exhaustive analysis of The Institute for Genomic Research (TIGR) and The Arabidopsis Information Resource (TAIR) databases, together with all currently compiled EST sequence data, was carried out in order to determine to what extent the various metabolic networks from phenylalanine ammonia lyase (PAL) to the monolignols were organized and/or could be predicted. In these databases, there are some 65 genes which have been annotated as encoding putative enzymatic steps in monolignol biosynthesis, although many of them have only very low homology to monolignol pathway genes of known function in other plant systems. Our detailed analysis revealed that presently only 13 genes (two PALs, a cinnamate-4-hydroxylase, a p-coumarate-3-hydroxylase, a ferulate-5-hydroxylase, three 4-coumarate-CoA ligases, a cinnamic acid O-methyl transferase, two cinnamoyl-CoA reductases) and two cinnamyl alcohol dehydrogenases can be classified as having a bona fide (definitive) function; the remaining 52 genes currently have undetermined physiological roles. The EST database entries for this particular set of genes also provided little new insight into how the monolignol pathway was organized in the different tissues and organs, this being perhaps a consequence of both limitations in how tissue samples were collected and in the incomplete nature of the EST collections. This analysis thus underscores the fact that even with genomic sequencing, presumed to provide the entire suite of putative genes in the monolignol-forming pathway, a very large effort needs to be conducted to establish actual catalytic roles (including enzyme versatility), as well as the physiological function(s) for each member of the (multi)gene families present and the metabolic networks that are operative. Additionally, one key to identifying physiological functions for many of these (and other) unknown genes, and their corresponding metabolic networks, awaits the development of technologies to comprehensively study molecular processes at the single cell level in particular tissues and organs, in order to establish the actual metabolic context.

NASA Program Fundamental Space Biology↗
DOE OSTI2023

Physiological potential and evolutionary trajectories of syntrophic sulfate-reducing bacterial partners of anaerobic methanotrophic archaea

Sulfate-coupled anaerobic oxidation of methane (AOM) is performed by multicellular consortia of anaerobic methanotrophic archaea (ANME) in obligate syntrophic partnership with sulfate-reducing bacteria (SRB). Diverse ANME and SRB clades co-associate but the physiological basis for their adaptation and diversification is not well understood. In this work, we used comparative metagenomics and phylogenetics to investigate the metabolic adaptation among the 4 main syntrophic SRB clades (HotSeep-1, Seep-SRB2, Seep-SRB1a, and Seep-SRB1g) and identified features associated with their syntrophic lifestyle that distinguish them from their non-syntrophic evolutionary neighbors in the phylum Desulfobacterota. We show that the protein complexes involved in direct interspecies electron transfer (DIET) from ANME to the SRB outer membrane are conserved between the syntrophic lineages. In contrast, the proteins involved in electron transfer within the SRB inner membrane differ between clades, indicative of convergent evolution in the adaptation to a syntrophic lifestyle. Our analysis suggests that in most cases, this adaptation likely occurred after the acquisition of the DIET complexes in an ancestral clade and involve horizontal gene transfers within pathways for electron transfer (CbcBA) and biofilm formation (Pel). We also provide evidence for unique adaptations within syntrophic SRB clades, which vary depending on the archaeal partner. Among the most widespread syntrophic SRB, Seep-SRB1a, subclades that specifically partner ANME-2a are missing the cobalamin synthesis pathway, suggestive of nutritional dependency on its partner, while closely related Seep-SRB1a partners of ANME-2c lack nutritional auxotrophies. Our work provides insight into the features associated with DIET-based syntrophy and the adaptation of SRB towards it.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Genome-Scale Metabolic Modeling Enables In-Depth Understanding of Big Data

Genome-scale metabolic models (GEMs) enable the mathematical simulation of the metabolism of archaea, bacteria, and eukaryotic organisms. GEMs quantitatively define a relationship between genotype and phenotype by contextualizing different types of Big Data (e.g., genomics, metabolomics, and transcriptomics). In this review, we analyze the available Big Data useful for metabolic modeling and compile the available GEM reconstruction tools that integrate Big Data. We also discuss recent applications in industry and research that include predicting phenotypes, elucidating metabolic pathways, producing industry-relevant chemicals, identifying drug targets, and generating knowledge to better understand host-associated diseases. In addition to the up-to-date review of GEMs currently available, we assessed a plethora of tools for developing new GEMs that include macromolecular expression and dynamic resolution. Finally, we provide a perspective in emerging areas, such as annotation, data managing, and machine learning, in which GEMs will play a key role in the further utilization of Big Data.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Modeling Growth Kinetics, Interspecies Cell Fusion, and Metabolism of a Clostridium acetobutylicum / Clostridium ljungdahlii Syntrophic Coculture

ABSTRACT Clostridium acetobutylicum and Clostridium ljungdahlii grown in a syntrophic culture were recently shown to fuse membranes and exchange cytosolic contents, yielding hybrid cells with significant shifts in gene expression and growth phenotypes. Here, we introduce a dynamic genome-scale metabolic modeling framework to explore how cell fusion alters the growth phenotype and panel of metabolites produced by this binary community. Computational results indicate C. ljungdahlii persists in the coculture through proteome exchange during fusing events, which endow C. ljungdahlii cells with expanded substrate utilization, and access to additional reducing equivalents from C. acetobutylicum -evolved H 2 and through acquisition of C. acetobutylicum -native cofactor-reducing enzymes. Simulations predict maximum theoretical ethanol and isopropanol yields that are increased by 0.64 mmol and 0.39 mmol per mmol hexose sugar consumed, respectively, during exponential growth when cell fusion is active. This modeling effort provides a mechanistic explanation for the metabolic outcome of cellular fusion and altered homeostasis achieved in this syntrophic clostridial community. IMPORTANCE Widespread cell fusion and protein exchange between microbial organisms as observed in synthetic C. acetobutylicum / C. ljungdahlii culture is a novel observation that has not been explored in silico . The mechanisms responsible for the observed cell fusion events in this culture are still unknown. In this work, we develop a modeling framework that captures the observed culture composition and metabolic phenotype, use it to offer a mechanistic explanation for how the culture achieves homeostasis, and identify C. ljungdahlii as primary beneficiary of fusion events. The implications for the events described in this study are far reaching, with potential to reshape our understanding of microbial community behavior synthetically and in nature.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Natural transformation as a tool in Acinetobacter baylyi : Streamlined engineering and mutational analysis

Natural transformation and homologous recombination in a soil bacterium, Acinetobacter baylyi ADP1, occur with exceptionally high efficiency. These genetic features can be harnessed to address a wide variety of fundamental and applied scientific topics. Recent advances in synthetic biology and laboratory evolution have led to renewed appreciation for the use of A. baylyi as a model organism. To complement several review articles that highlight new tool sets, this chapter focuses on simple protocols and examples of transformation assays that facilitate genetic analysis and engineering. Whole genome sequencing often reveals extensive genetic variation among closely related isolates that can confound the association of genotypic and phenotypic changes. In A. baylyi, such associations can be deciphered in unique ways by directly transforming cells with linear DNA fragments. The resulting allelic replacement, which occurs at high frequency, rapidly generates desired mutants via targeted chromosomal editing. Diverse screening and selection methods can be used to test hypotheses and streamline experimental strategies to reveal the significance of specific DNA sequences. Moreover, large procedural variations are well tolerated, and techniques can be readily adapted for new purposes. Furthermore, one goal of highlighting natural transformation methodology in A. baylyi is to expand the community of researchers using this versatile bacterial host.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Sequence Analysis of 20,453 Severe Acute Respiratory Syndrome Coronavirus 2 Genomes from the Houston Metropolitan Area Identifies the Emergence and Widespread Distribution of Multiple Isolates of All Major Variants of Concern

Since the beginning of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic, there has been international concern about the emergence of virus variants with mutations that increase transmissibility, enhance escape from the human immune response, or otherwise alter biologically important phenotypes. In late 2020, several variants of concern emerged globally, including the UK variant (B.1.1.7), the South Africa variant (B.1.351), Brazil variants (P.1 and P.2), and two related California variants of interest (B.1.429 and B.1.427). These variants are believed to have enhanced transmissibility. For the South Africa and Brazil variants, there is evidence that mutations in spike protein permit it to escape from some vaccines and therapeutic monoclonal antibodies. On the basis of our extensive genome sequencing program involving 20,453 coronavirus disease 2019 patient samples collected from March 2020 to February 2021, we report identification of all six of these SARS-CoV-2 variants among Houston Methodist Hospital (Houston, TX) patients residing in the greater metropolitan area. Although these variants are currently at relatively low frequency (aggregate of 1.1%) in the population, they are geographically widespread. Houston is the first city in the United States in which active circulation of all six current variants of concern has been documented by genome sequencing. Finally, as vaccine deployment accelerates, increased genomic surveillance of SARS-CoV-2 is essential to understanding the presence, frequency, and medical impact of consequential variants and their patterns and trajectory of dissemination.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

DYNC1H1 regulates NSCLC cell growth and metastasis by IFN-γ-JAK-STAT signaling and is associated with an aberrant immune response

Highlights: • DYNC1H1 is upregulated and its high expression predicted poor outcomes in NSCLC. • Loss of function of DYNC1H1 inhibits cell growth and metastasis of NSCLC cells in vitro. • DYNC1H1 knockdown inhibited metastasis of NSCLC cells in vivo. • DYNC1H1 loss of function suppresses IFN-γ-JAK-STAT signaling pathway in NSCLC. • DYNC1H1 expression is correlated with aberrant immune response in NSCLC indicated by analysis of TCGA immune data. It is urgent to identify new biomarkers and therapeutic targets to ameliorate the clinical prognosis of patients with lung cancer. The functional significance and molecular mechanism of dynein cytoplasmic 1 heavy chain 1 (DYNC1H1) in nonsmall cell lung cancer (NSCLC) progression is still elusive. In our current study, publicly available data and Western blotting experiments confirmed that DYNC1H1 expression was upregulated in lung cancer samples compared with noncancerous samples. Quantitative real-time PCR (qPCR) results indicated that high DYNC1H1 expression in lung cancer tissues was significantly associated with clinical tumor stage and distal metastasis; moreover, its high expression was negatively correlated with prognosis. Functional experiments demonstrated that DYNC1H1 loss of function caused a significant decrease in cell viability and cell proliferative ability, inhibition of the cell cycle, and promotion of both migration potential and invasion potential in vitro. Animal experiments by tail vein injection of lung cancer cells showed that DYNC1H1 knockdown significantly decreased lung cancer metastasis. Mechanistically, the results from a human protein array showed changes in the IFN-γ-JAK-STAT signaling pathway, and analysis of The Cancer Genome Atlas (TCGA) immune data demonstrated that disturbance of the immune microenvironment might be involved in the impaired growth and metastatic ability mediated by DYNC1H1 loss in NSCLC. DYNC1H1 might serve as a promising biological marker of prognosis and a potential clinical therapeutic target for patients with NSCLC.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Metabolic engineering of Synechococcus elongatus 7942 for enhanced sucrose biosynthesis

The capability of cyanobacteria to produce sucrose from CO 2 and light has a remarkable societal and biotechnological impact since sucrose can serve as a carbon and energy source for a variety of heterotrophic organisms and can be converted into value-added products. However, most metabolic engineering efforts have focused on understanding local pathway alterations that drive sucrose biosynthesis and secretion in cyanobacteria rather than analyzing the global flux re-routing that occurs following induction of sucrose production by salt stress. Here, we investigated global metabolic flux alterations in a sucrose-secreting (cscB-overexpressing) strain relative to its wild-type Synechococcus elongatus 7942 parental strain. Here we used targeted metabolomics, 13 C metabolic flux analysis (MFA), and genome-scale modeling (GSM) as complementary approaches to elucidate differences in cellular resource allocation by quantifying metabolic profiles of three cyanobacterial cultures - wild-type S. elongatus 7942 without salt stress (WT), wild-type with salt stress (WT/NaCl), and the cscB-overexpressing strain with salt stress (cscB/NaCl) - all under photoautotrophic conditions. We quantified the substantial rewiring of metabolic fluxes in WT/NaCl and cscB/NaCl cultures relative to WT and identified a metabolic bottleneck limiting carbon fixation and sucrose biosynthesis. This bottleneck was subsequently mitigated through heterologous overexpression of glyceraldehyde-3-phosphate dehydrogenase in an engineered sucrose-secreting strain. Our study also demonstrates that combining 13 C-MFA and GSM is a useful strategy to both extend the coverage of MFA beyond central metabolism and to improve the accuracy of flux predictions provided by GSM.

09 BIOMASS FUELS↗
DOE OSTI2024

RNA language models predict mutations that improve RNA function

Structured RNA lies at the heart of many central biological processes, from gene expression to catalysis. RNA structure prediction is not yet possible due to a lack of high-quality reference data associated with organismal phenotypes that could inform RNA function. We present GARNET (Gtdb Acquired RNa with Environmental Temperatures), a new database for RNA structural and functional analysis anchored to the Genome Taxonomy Database (GTDB). GARNET links RNA sequences to experimental and predicted optimal growth temperatures of GTDB reference organisms. Using GARNET, we develop sequence- and structure-aware RNA generative models, with overlapping triplet tokenization providing optimal encoding for a GPT-like model. Leveraging hyperthermophilic RNAs in GARNET and these RNA generative models, we identify mutations in ribosomal RNA that confer increased thermostability to the Escherichia coli ribosome. The GTDB-derived data and deep learning models presented here provide a foundation for understanding the connections between RNA sequence, structure, and function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Unraveling the functional dark matter through global metagenomics

Metagenomes encode an enormous diversity of proteins, reflecting a multiplicity of functions and activities1,2. Exploration of this vast sequence space has been limited to a comparative analysis against reference microbial genomes and protein families derived from those genomes. Here, to examine the scale of yet untapped functional diversity beyond what is currently possible through the lens of reference genomes, we develop a computational approach to generate reference-free protein families from the sequence space in metagenomes. We analyse 26,931 metagenomes and identify 1.17 billion protein sequences longer than 35 amino acids with no similarity to any sequences from 102,491 reference genomes or the Pfam database3. Using massively parallel graph-based clustering, we group these proteins into 106,198 novel sequence clusters with more than 100 members, doubling the number of protein families obtained from the reference genomes clustered using the same approach. We annotate these families on the basis of their taxonomic, habitat, geographical and gene neighbourhood distributions and, where sufficient sequence diversity is available, predict protein three-dimensional models, revealing novel structures. Overall, our results uncover an enormously diverse functional space, highlighting the importance of further exploring the microbial functional dark matter.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A robust benchmark for detection of germline large deletions and insertions

New technologies and analysis methods are enabling genomic structural variants (SVs) to be detected with ever-increasing accuracy, resolution and comprehensiveness. To help translate these methods to routine research and clinical practice, we developed a sequence-resolved benchmark set for identification of both false-negative and false-positive germline large insertions and deletions. To create this benchmark for a broadly consented son in a Personal Genome Project trio with broadly available cells and DNA, the Genome in a Bottle Consortium integrated 19 sequence-resolved variant calling methods from diverse technologies. The final benchmark set contains 12,745 isolated, sequence-resolved insertion (7,281) and deletion (5,464) calls ≥50 base pairs (bp). The Tier 1 benchmark regions, for which any extra calls are putative false positives, cover 2.51 Gbp and 5,262 insertions and 4,095 deletions supported by ≥1 diploid assembly. We demonstrate that the benchmark set reliably identifies false negatives and false positives in high-quality SV callsets from short-, linked- and long-read sequencing and optical mapping.

59 BASIC BIOLOGICAL SCIENCES↗