Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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.

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Author Correction: Expanded encyclopaedias of DNA elements in the human and mouse genomes

In the version of this article initially published, two members of the ENCODE Project Consortium were missing from the author list. Rizi Ai (Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, CA, USA) and Shantao Li (Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA) are now included in the author list. These errors have been corrected in the online version of the article.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Shotgun metagenome data of a defined mock community using Oxford Nanopore, PacBio and Illumina technologies

Metagenomic sequence data from defined mock communities is crucial for the assessment of sequencing platform performance and downstream analyses, including assembly, binning and taxonomic assignment. We report a comparison of shotgun metagenome sequencing and assembly metrics of a defined microbial mock community using the Oxford Nanopore Technologies (ONT) MinION, PacBio and Illumina sequencing platforms. Our synthetic microbial community BMock12 consists of 12 bacterial strains with genome sizes spanning 3.2-7.2 Mbp, 40-73% GC content, and 1.5-7.3% repeats. Size selection of both PacBio and ONT sequencing libraries prior to sequencing was essential to yield comparable relative abundances of organisms among all sequencing technologies. While the Illumina-based metagenome assembly yielded good coverage with few misassemblies, contiguity was greatly improved by both, Illumina + ONT and Illumina + PacBio hybrid assemblies but increased misassemblies, most notably in genomes with high sequence similarity to each other. Our resulting datasets allow evaluation and benchmarking of bioinformatics software on Illumina, PacBio and ONT platforms in parallel.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Improvement of eukaryotic protein predictions from soil metagenomes

During the last decades, metagenomics has highlighted the diversity of microorganisms from environmental or host-associated samples. Most metagenomics public repositories use annotation pipelines tailored for prokaryotes regardless of the taxonomic origin of contigs. Consequently, eukaryotic contigs with intrinsically different gene features, are not optimally annotated. Using a bioinformatics pipeline, we have filtered 7.9 billion contigs from 6,872 soil metagenomes in the JGI’s IMG/M database to identify eukaryotic contigs. We have re-annotated genes using eukaryote-tailored methods, yielding 8 million eukaryotic proteins and over 300,000 orphan proteins lacking homology in public databases. Comparing the gene predictions we made with initial JGI ones on the same contigs, we confirmed our pipeline improves eukaryotic proteins completeness and contiguity in soil metagenomes. The improved quality of eukaryotic proteins combined with a more comprehensive assignment method yielded more reliable taxonomic annotation. This dataset of eukaryotic soil proteins with improved completeness, quality and taxonomic annotation reliability is of interest for any scientist aiming at studying the composition, biological functions and gene flux in soil communities involving eukaryotes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Benchmarking second and third-generation sequencing platforms for microbial metagenomics

Shotgun metagenomic sequencing is a common approach for studying the taxonomic diversity and metabolic potential of complex microbial communities. Current methods primarily use second generation short read sequencing, yet advances in third generation long read technologies provide opportunities to overcome some of the limitations of short read sequencing. Here, we compared seven platforms, encompassing second generation sequencers (Illumina HiSeq 300, MGI DNBSEQ-G400 and DNBSEQ-T7, ThermoFisher Ion GeneStudio S5 and Ion Proton P1) and third generation sequencers (Oxford Nanopore Technologies MinION R9 and Pacific Biosciences Sequel II). We constructed three uneven synthetic microbial communities composed of up to 87 genomic microbial strains DNAs per mock, spanning 29 bacterial and archaeal phyla, and representing the most complex and diverse synthetic communities used for sequencing technology comparisons. Our results demonstrate that third generation sequencing have advantages over second generation platforms in analyzing complex microbial communities, but require careful sequencing library preparation for optimal quantitative metagenomic analysis. Our sequencing data also provides a valuable resource for testing and benchmarking bioinformatics software for metagenomics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Hypersaline sapropels act as hotspots for microbial dark matter

Present-day terrestrial analogue sites are crucial ground truth proxies for studying life in geochemical conditions close to those assumed to be present on early Earth or inferred to exist on other celestial bodies (e.g. Mars, Europa). Although hypersaline sapropels are border-of-life habitats with moderate occurrence, their microbiological and physicochemical characterization lags behind. Here, we study the diversity of life under low water activity by describing the prokaryotic communities from two disparate hypersaline sapropels (Transylvanian Basin, Romania) in relation to geochemical milieu and pore water chemistry, while inferring their role in carbon cycling by matching taxa to known taxon-specific biogeochemical functions. Furthermore, the polyphasic approach combined deep coverage SSU rRNA gene amplicon sequencing and bioinformatics with RT-qPCR and physicochemical investigations. We found that sapropels developed an analogous elemental milieu and harbored prokaryotes affiliated with fifty-nine phyla, among which the most abundant were Proteobacteria , Bacteroidetes and Chloroflexi . Containing thirty-two candidate divisions and possibly undocumented prokaryotic lineages, the hypersaline sapropels were found to accommodate one of the most diverse and novel ecosystems reported to date and may contribute to completing the phylogenetic branching of the tree of life.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Probabilistic variable-length segmentation of protein sequences for discriminative motif discovery (DiMotif) and sequence embedding (ProtVecX)

In this paper, we present peptide-pair encoding (PPE), a general-purpose probabilistic segmentation of protein sequences into commonly occurring variable-length sub-sequences. The idea of PPE segmentation is inspired by the byte-pair encoding (BPE) text compression algorithm, which has recently gained popularity in subword neural machine translation. We modify this algorithm by adding a sampling framework allowing for multiple ways of segmenting a sequence. PPE segmentation steps can be learned over a large set of protein sequences (Swiss-Prot) or even a domain-specific dataset and then applied to a set of unseen sequences. This representation can be widely used as the input to any downstream machine learning tasks in protein bioinformatics. In particular, here, we introduce this representation through protein motif discovery and protein sequence embedding. (i) DiMotif: we present DiMotif as an alignment-free discriminative motif discovery method and evaluate the method for finding protein motifs in three different settings: (1) comparison of DiMotif with two existing approaches on 20 distinct motif discovery problems which are experimentally verified, (2) classification-based approach for the motifs extracted for integrins, integrin-binding proteins, and biofilm formation, and (3) in sequence pattern searching for nuclear localization signal. The DiMotif, in general, obtained high recall scores, while having a comparable F1 score with other methods in the discovery of experimentally verified motifs. Having high recall suggests that the DiMotif can be used for short-list creation for further experimental investigations on motifs. In the classification-based evaluation, the extracted motifs could reliably detect the integrins, integrin-binding, and biofilm formation-related proteins on a reserved set of sequences with high F1 scores. (ii) ProtVecX: we extend k-mer based protein vector (ProtVec) embedding to variablelength protein embedding using PPE sub-sequences. We show that the new method of embedding can marginally outperform ProtVec in enzyme prediction as well as toxin prediction tasks. In addition, we conclude that the embeddings are beneficial in protein classification tasks when they are combined with raw amino acids k-mer features.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Standardized phylogenetic and molecular evolutionary analysis applied to species across the microbial tree of life

There is growing interest in reconstructing phylogenies from the copious amounts of genome sequencing projects that target related viral, bacterial or eukaryotic organisms. To facilitate the construction of standardized and robust phylogenies for disparate types of projects, we have developed a complete bioinformatic workflow, with a web-based component to perform phylogenetic and molecular evolutionary (PhaME) analysis from sequencing reads, draft assemblies or completed genomes of closely related organisms. Furthermore, the ability to incorporate raw data, including some metagenomic samples containing a target organism (e.g. from clinical samples with suspected infectious agents), shows promise for the rapid phylogenetic characterization of organisms within complex samples without the need for prior assembly.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Structural analysis of a new carotenoid-binding protein: the C-terminal domain homolog of the OCP

Abstract The Orange Carotenoid Protein (OCP) is a water-soluble protein that governs photoprotection in many cyanobacteria. The 35 kDa OCP is structurally and functionally modular, consisting of an N-terminal effector domain (NTD) and a C-terminal regulatory domain (CTD); a carotenoid spans the two domains. The CTD is a member of the ubiquitous Nuclear Transport Factor-2 (NTF2) superfamily (pfam02136). With the increasing availability of cyanobacterial genomes, bioinformatic analysis has revealed the existence of a new family of proteins, homologs to the CTD, the C-terminal domain-like carotenoid proteins (CCPs). Here we purify holo-CCP2 directly from cyanobacteria and establish that it natively binds canthaxanthin (CAN). We use small-angle X-ray scattering (SAXS) to characterize the structure of this carotenoprotein in two distinct oligomeric states. A single carotenoid molecule spans the two CCPs in the dimer. Our analysis with X-ray footprinting-mass spectrometry (XFMS) identifies critical residues for carotenoid binding that likely contribute to the extreme red shift (ca. 80 nm) of the absorption maximum of the carotenoid bound by the CCP2 dimer and a further 10 nm shift in the tetramer form. These data provide the first structural description of carotenoid binding by a protein consisting of only an NTF2 domain.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

MULTICOM2 open-source protein structure prediction system powered by deep learning and distance prediction

Protein structure prediction is an important problem in bioinformatics and has been studied for decades. However, there are still few open-source comprehensive protein structure prediction packages publicly available in the field. In this paper, we present our latest open-source protein tertiary structure prediction system—MULTICOM2, an integration of template-based modeling (TBM) and template-free modeling (FM) methods. The template-based modeling uses sequence alignment tools with deep multiple sequence alignments to search for structural templates, which are much faster and more accurate than MULTICOM1. The template-free (ab initio or de novo) modeling uses the inter-residue distances predicted by DeepDist to reconstruct tertiary structure models without using any known structure as template. In the blind CASP14 experiment, the average TM-score of the models predicted by our server predictor based on the MULTICOM2 system is 0.720 for 58 TBM (regular) domains and 0.514 for 38 FM and FM/TBM (hard) domains, indicating that MULTICOM2 is capable of predicting good tertiary structures across the board. It can predict the correct fold for 76 CASP14 domains (95% regular domains and 55% hard domains) if only one prediction is made for a domain. The success rate is increased to 3% for both regular and hard domains if five predictions are made per domain. Moreover, the prediction accuracy of the pure template-free structure modeling method on both TBM and FM targets is very close to the combination of template-based and template-free modeling methods. This demonstrates that the distance-based template-free modeling method powered by deep learning can largely replace the traditional template-based modeling method even on TBM targets that TBM methods used to dominate and therefore provides a uniform structure modeling approach to any protein. Finally, on the 38 CASP14 FM and FM/TBM hard domains, MULTICOM2 server predictors (MULTICOM-HYBRID, MULTICOM-DEEP, MULTICOM-DIST) were ranked among the top 20 automated server predictors in the CASP14 experiment. After combining multiple predictors from the same research group as one entry, MULTICOM-HYBRID was ranked no. 5. The source code of MULTICOM2 is freely available at https://github.com/multicom-toolbox/multicom/tree/multicom_v2.0 .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Transposon sequencing analysis of Bradyrhizobium diazoefficiens 110 spc4

Bradyrhizobium diazoefficiens USDA110 is one of the most effective nitrogen-fixing symbionts of soybeans. Here we carried out a large-scale transposon insertion sequencing (Tn-seq) analysis of strain Bd110spc4, which is derived from USDA110, with the goal of increasing available resources for identifying genes crucial for the survival of this plant symbiont under diverse conditions. We prepared two transposon (Tn) insertion libraries of Bd110spc4 with 155,042 unique Tn insertions when the libraries were combined, which is an average of one insertion every 58.7 bp of the reference USDA110 genome. Application of bioinformatic filtering steps to remove genes too small to be expected to have Tn insertions, resulted in a list of genes that were classified as putatively essential. Comparison of this gene set with genes putatively essential for the growth of the closely related alpha-proteobacterium, Rhodopseudomonas palustris, revealed a small set of five genes that may be collectively essential for closely related members of the family Bradyrhizobiaceae. This group includes bacteria with diverse lifestyles ranging from plant symbionts to animal-associated species to free-living species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Sequence similarity between SARS-CoV-2 nucleocapsid and multiple sclerosis-associated proteins provides insight into viral neuropathogenesis following infection

Abstract The novel coronavirus SARS-CoV-2 continues to cause death and disease throughout the world, underscoring the necessity of understanding the virus and host immune response. From the start of the pandemic, a prominent pattern of central nervous system (CNS) pathologies, including demyelination, has emerged, suggesting an underlying mechanism of viral mimicry to CNS proteins. We hypothesized that immunodominant epitopes of SARS-CoV-2 share homology with proteins associated with multiple sclerosis (MS). Using PEPMatch, a newly developed bioinformatics package which predicts peptide similarity within specific amino acid mismatching parameters consistent with published MHC binding capacity, we discovered that nucleocapsid protein shares significant overlap with 22 MS-associated proteins, including myelin proteolipid protein (PLP). Further computational evaluation demonstrated that this overlap may have critical implications for T cell responses in MS patients and is likely unique to SARS-CoV-2 among the major human coronaviruses. Our findings substantiate the hypothesis of viral molecular mimicry in the pathogenesis of MS and warrant further experimental exploration.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

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

Widespread bacterial diversity within the bacteriome of fungi

Abstract Knowledge of associations between fungal hosts and their bacterial associates has steadily grown in recent years as the number and diversity of examinations have increased, but current knowledge is predominantly limited to a small number of fungal taxa and bacterial partners. Here, we screened for potential bacterial associates in over 700 phylogenetically diverse fungal isolates, representing 366 genera, or a tenfold increase compared with previously examined fungal genera, including isolates from several previously unexplored phyla. Both a 16 S rDNA-based exploration of fungal isolates from four distinct culture collections spanning North America, South America and Europe, and a bioinformatic screen for bacterial-specific sequences within fungal genome sequencing projects, revealed that a surprisingly diverse array of bacterial associates are frequently found in otherwise axenic fungal cultures. We demonstrate that bacterial associations with diverse fungal hosts appear to be the rule, rather than the exception, and deserve increased consideration in microbiome studies and in examinations of microbial interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Structural and mechanistic basis for redox sensing by the cyanobacterial transcription regulator RexT

Abstract Organisms have a myriad of strategies for sensing, responding to, and combating reactive oxygen species, which are unavoidable consequences of aerobic life. In the heterocystous cyanobacterium Nostoc sp. PCC 7120, one such strategy is the use of an ArsR-SmtB transcriptional regulator RexT that senses H 2 O 2 and upregulates expression of thioredoxin to maintain cellular redox homeostasis. Different from many other members of the ArsR-SmtB family which bind metal ions, RexT has been proposed to use disulfide bond formation as a trigger to bind and release DNA. Here, we present high-resolution crystal structures of RexT in the reduced and H 2 O 2 -treated states. These structures reveal that RexT showcases the ArsR-SmtB winged-helix-turn-helix fold and forms a vicinal disulfide bond to orchestrate a response to H 2 O 2 . The importance of the disulfide-forming Cys residues was corroborated using site-directed mutagenesis, mass spectrometry, and H 2 O 2 -consumption assays. Furthermore, an entrance channel for H 2 O 2 was identified and key residues implicated in H 2 O 2 activation were pinpointed. Finally, bioinformatics analysis of the ArsR-SmtB family indicates that the vicinal disulfide “redox switch” is a unique feature of cyanobacteria in the Nostocales order, presenting an interesting case where an ArsR-SmtB protein scaffold has been evolved to showcase peroxidatic activity and facilitate redox-based regulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Characterization of a widespread sugar phosphate-processing bacterial microcompartment

Many prokaryotes form Bacterial Microcompartments (BMCs) that encapsulate segments of specialized metabolic pathways to enhance catalysis. The various functions of metabolosomes, catabolic BMCs, are dictated by the signature enzyme that processes initial substrates of the confined pathway. The components and native functions of several metabolosomes have been experimentally characterized; however one of the most prevalent across all bacteria has yet to be studied. Sugar Phosphate Utilizing (SPU) BMC loci encode enzymes predicted to be involved in sugar phosphate metabolism. The SPU genetic loci are found in organisms occupying habitats ranging from soils to hot springs, highlighting the ubiquity of the SPU BMC. We bioinformatically characterized seven SPU subtypes, all which contain an enzyme unique to SPU BMCs, a deoxyribose 5-phosphate aldolase (DERA). Here, we define the fundamental characteristics of SPU BMCs and have expressed, purified, and characterized a set of SPU core enzymes. These include a protein-protein complex formed between a SPU BMC DERA and a predicted ribose 5-phosphate isomerase. Further, we show that the SPU BMC DERA is catalytically active and propose that it acts as the universal signature enzyme for the SPU BMC, with implications for fundamental understanding and biotechnological applications of SPU BMCs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Sequence Data for Clostridium autoethanogenum using Three Generations of Sequencing Technologies

During the past decade, DNA sequencing output has been mostly dominated by the second generation sequencing platforms which are characterized by low cost, high throughput and shorter read lengths for example, Illumina. The emergence and development of so called third generation sequencing platforms such as PacBio has permitted exceptionally long reads (over 20 kb) to be generated. Due to read length increases, algorithm improvements and hybrid assembly approaches, the concept of one chromosome, one contig and automated finishing of microbial genomes is now a realistic and achievable task for many microbial laboratories. In this paper, we describe high quality sequence datasets which span three generations of sequencing technologies, containing six types of data from four NGS platforms and originating from a single microorganism, Clostridium autoethanogenum. The dataset reported here will be useful for the scientific community to evaluate upcoming NGS platforms, enabling comparison of existing and novel bioinformatics approaches and will encourage interest in the development of innovative experimental and computational methods for NGS data.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Next generation sequencing data of a defined microbial mock community

Generating sequence data of a defined community composed of organisms with complete reference genomes is indispensable for the benchmarking of new genome sequence analysis methods, including assembly and binning tools. Moreover the validation of new sequencing library protocols and platforms to assess critical components such as sequencing errors and biases relies on such datasets. We here report the next generation metagenomic sequence data of a defined mock community (Mock Bacteria ARchaea Community; MBARC-26), composed of 23 bacterial and 3 archaeal strains with finished genomes. These strains span 10 phyla and 14 classes, a range of GC contents, genome sizes, repeat content and encompass a diverse abundance profile. Short read Illumina and long-read PacBio SMRT sequences of this mock community are described. These data represent a valuable resource for the scientific community, enabling extensive benchmarking and comparative evaluation of bioinformatics tools without the need to simulate data. As such, these data can aid in improving our current sequence data analysis toolkit and spur interest in the development of new tools.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Rapid genome resequencing of an atoxigenic strain of Aspergillus carbonarius

In microorganisms, Ion Torrent sequencing technology has been proved to be useful in whole-genome sequencing of bacterial genomes (5 Mbp). In our study, for the first time we used this technology to perform a resequencing approach in a whole fungal genome (36 Mbp), a non-ochratoxin A producing strain of Aspergillus carbonarius . Ochratoxin A (OTA) is a potent nephrotoxin which is found mainly in cereals and their products, but it also occurs in a variety of common foods and beverages. Due to the fact that this strain does not produce OTA, we focused some of the bioinformatics analyses in genes involved in OTA biosynthesis, using a reference genome of an OTA producing strain of the same species. This study revealed that in the atoxigenic strain there is a high accumulation of nonsense and missense mutations in several genes. Importantly, a two fold increase in gene mutation ratio was observed in PKS and NRPS encoding genes which are suggested to be involved in OTA biosynthesis.

59 BASIC BIOLOGICAL SCIENCES↗