Search NASA⌕ Search

SEARCH · Search NASA

Results for “sequence 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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Nuclear Data Libraries Sensitivity Studies for ORSA Using SCALE

Subcritical assemblies offer valuable training capabilities in nuclear criticality safety (NCS) for individuals handling fissile material. At Oak Ridge National Laboratory(ORNL), the Oak Ridge Subcritical Assembly (ORSA), a new experimental facility, is being established to provide hands-on training for the Nuclear Criticality Safety Program (NCSP).It is essential to accurately determine the neutron multiplication factor (keff) to ensure that ORSA remains subcritical and safe during operations. This study investigated the sensitivity of keff to variations across Evaluated Nuclear Data File (ENDF/B) libraries, consisting of ENDF/B-VII.1, ENDF/B-VIII.0, and ENDF/B-VIII.1. The analysis was conducted using the CSAS6 sequence in the SCALE-6.3 code system. Individual isotopes in the ORSA model were replaced one at a time with ENDF/B-VII.1 as the base library and changing to ENDF/B-VIII.0or ENDF/B-VIII.1. The results demonstrated that the changes in keffof the nuclides associated with the ORSA model were mostly within the uncertainty of the base model (~24 pcm), except for primary nuclides like Uranium-235and H-poly with few other nuclides. The relative delta keff values of Uranium-235 and H-poly, expressed in pcm, were +280 and -212 in ENDF/B-VIII.0 and +332 and -275 in ENDF/B-VIII.1, respectively, which were notable changes in reactivity. These results demonstrated that ORSA was largely insensitive to variations across these nuclear data libraries.

Hong, Evan [North Carolina State University]↗
DOE OSTI2016

Down-regulation of TCF21 by hypermethylation induces cell proliferation, migration and invasion in colorectal cancer

Epigenetic alteration induced loss function of the transcription factor 21 (TCF21) has been associated with different types of human cancers. However, the epigenetic regulation and molecular functions of TCF21 in colorectal cancer (CRC) remain unknown. In this study, TCF21 expression levels and methylation status of its promoter region in CRC cell lines (n = 5) and CRC tissues (n = 151) as well as normal colorectal mucosa (n = 30) were assessed by RTq-PCR and methylation analysis (methylation specific PCR, MSP and bisulfite sequencing PCR, BSP), respectively. The cellular functions of TCF21 on CRC cell proliferation, apoptosis, invasion and migration were investigated in vitro. Our data revealed that TCF21 was frequently silenced by promoter hypermethylation in both tested CRC cell lines and primary CRC, and correlation analysis between methylation status and clinicopathologic parameters found that TCF21 methylation was significantly correlated with lymph node invasion (P = 0.013), while no significant correlation was found in other parameters. In addition, demethylation treatment resulted in re-expression of TCF21 in CRC cell lines, and cellular function experiments revealed that restoration of TCF21 inhibited CRC cell proliferation, promoted apoptosis and suppressed cell invasion and migration, suggesting that TCF21 may function as a tumor suppressor gene, which is downregulated through promoter hypermethylation in CRC development. - Highlights: • TCF21 was frequently silenced by promoter DNA methylation in CRC cells. • TCF21 was frequently methylated in primary CRC and significantly correlated with metastasis. • Restoration of TCF21 promotes cell apoptosis of CRC cells. • Restoration of TCF21 inhibits cell invasion and migration of CRC cells.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Microbial Ecology on Solar Panels in Berkeley, CA, United States

Solar panels can be found practically all over the world and represent a standard surface that can be colonized by microbial communities that are resistant to harsh environmental conditions, including high irradiation, temperature fluctuations and desiccation. These properties make them not only ideal sources of stress-resistant bacteria, but also standard devices to study the microbial communities and their colonization process from different areas of Earth. We report here a comprehensive description of the microbial communities associated with solar panels in Berkeley, CA, United States. Cultivable bacteria were isolated to characterize their adhesive capabilities, and UV- and desiccation-resistance properties. Furthermore, a parallel culture-independent metagenomic and metabolomic approach has allowed us to gain insight on the taxonomic and functional nature of these communities. Metagenomic analysis was performed using the Illumina HiSeq2500 sequencing platform, revealing that the bacterial population of the Berkeley solar panels is composed mainly of Actinobacteria, Bacteroidetes and Proteobacteria, as well as lower amounts of Deinococcus-Thermus and Firmicutes. Furthermore, a clear predominance of Hymenobacter sp. was also observed. A functional analysis revealed that pathways involved in the persistence of microbes on solar panels (i.e., stress response, capsule development, and metabolite repair) and genes assigned to carotenoid biosynthesis were common to all metagenomes. On the other hand, genes involved in photosynthetic pathways and general autotrophic subsystems were rare, suggesting that these pathways are not critical for persistence on solar panels. Metabolomics was performed using a liquid chromatography tandem mass spectrometry (LC-MS/MS) approach. When comparing the metabolome of the solar panels from Berkeley and from Valencia (Spain), a very similar composition in polar metabolites could be observed, although some metabolites appeared to be differentially represented (for example, trigonelline, pantolactone and 5-valerolactone were more abundant in the samples from Valencia than in the ones from Berkeley). Furthermore, triglyceride metabolites were highly abundant in all the solar panel samples, and both locations displayed similar profiles. The comparison of the taxonomic profile of the Californian solar panels with those previously described in Spain revealed striking similarities, highlighting the central role of both selective pressures and the ubiquity of microbial populations in the colonization and establishment of microbial communities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

High-quality Acinetobacter genomes recovered from combat wounds via metagenomic sequencing resemble cultured isolate genomes

The ability to accurately characterize wound pathogens is critical to informing clinical decisions for wound infections with complex treatment requirements. Acinetobacter baumannii is an impactful nosocomial pathogen in combat wounds and civilian hospital-acquired infections. An informed understanding of the phylogenetics and epidemiology of A. baumannii infections in military and civilian environments could guide approaches that improve antibiotic treatment regimens for both military and civilian patients. Whole-genome data for bacterial strains can be difficult to obtain due to challenges in culturing isolates from preserved military specimens. Metagenomic sequencing and assembly create opportunities for genomic analysis of pathogens directly from clinical specimens. The ability to perform comparative analyses between metagenome-derived genomes and culture-derived genomes would support a range of comparative bacterial genomic studies. Wound tissue biopsy and effluent samples from combat injuries were subjected to metagenomic sequencing and assembly. In total, 42 microbial metagenome-assembled genomes (MAGs) were obtained directly from metagenomic sequence data, 36 of which were designated “high” quality. Thirty of these genomes corresponded to Acinetobacter, with 29 mapping specifically to A. baumannii. Other observed genera included Bordetella, Citrobacter, Escherichia, and Pseudomonas. Single-copy and multi-copy orthologs were identified across Acinetobacter MAGs and publicly available isolate genomes derived from military and civilian sources. Both MAG and military isolate genomes were annotated with antimicrobial resistance data, and MAG genomes were statistically comparable to genomes obtained from isolates. Our results highlight the potential of de novo metagenome assembly for enabling high-resolution characterization directly from clinical specimens, thereby improving diagnostic precision, guiding antimicrobial stewardship, and enhancing understanding of pathogen evolution across diverse healthcare and battlefield environments.

Acinetobacter baumannii↗
DOE OSTI2018

Genus-wide assessment of lignocellulose utilization in the extremely thermophilic Caldicellulosiruptor by genomic, pan-genomic and metagenomic analysis

Metagenomic data from Obsidian Pool (Yellowstone National Park, USA) and thirteen genome sequences were used to re-assess genus-wide biodiversity for the extremely thermophilic Caldicellulosiruptor . The updated core-genome contains 1,401 ortholog groups (average genome size for thirteen species = 2,516 genes). The pan-genome, which remains open with a revised total of 3,493 ortholog groups, encodes a variety of multi-domain glycoside hydrolases (GH). These include three cellulases with GH48 domains that are co-located in the Glucan Degradation Locus (GDL) and are specific determinants for microcrystalline cellulose utilization. Three recently sequenced species, Caldicellulosiruptor sp. str. Rt8.B8 (re-named here Caldicellulosiruptor morganii ), Thermoanaerobacter cellulolyticus str. NA10 (re-named here Caldicellulosiruptor naganoensis NA10), and Caldicellulosiruptor sp. str. Wai35.B1 (re-named here Caldicellulosiruptor danielii ) degraded Avicel and lignocellulose (switchgrass). C. morganii was more efficient than C. bescii in this regard and differed from the other twelve species examined here, both based on genome content and organization and in the specific domain features of conserved GHs. Metagenomic analysis of lignocellulose-enriched samples from Obsidian Pool revealed limited new information on genus biodiversity. Enrichments yielded genomic signatures closely related to Caldicellulosiruptor obsidiansis , but there was also evidence for other thermophilic fermentative anaerobes ( Caldanaerobacter , Fervidobacterium , Caloramator , and Clostridium ). One enrichment, containing 89.7% Caldicellulosiruptor and 9.7% Caloramator , had a capacity for switchgrass solubilization comparable to C. bescii . These results refine the known biodiversity of Caldicellulosiruptor and indicate that microcrystalline cellulose degradation at temperatures above 70°C, based on current information, is limited to certain members of this genus that produce GH48 domain-containing enzymes. The genus Caldicellulosiruptor contains the most thermophilic bacteria capable of lignocellulose deconstruction and are promising candidates for consolidated bioprocessing for the production of biofuels and bio-based chemicals. The focus here is on the extant capability of this genus for plant biomass degradation and the extent to which this can be inferred from the core and pan-genomes, based on analysis of thirteen species and metagenomic sequence information from environmental samples. Key to microcrystalline hydrolysis is the content of the Glucan Degradation Locus (GDL), a set of genes encoding glycoside hydrolases (GH), several of which have GH48 and family 3 carbohydrate binding module domains, that function as primary cellulases. Resolving the relationship between the GDL and lignocellulose degradation will inform efforts to identify more prolific members of the genus and to develop metabolic engineering strategies to improve this characteristic

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Circumvention of safety related equipment at nuclear fuel services, Inc. - 14570

On June 17, 2014, at Nuclear Fuel Services, Inc. (NFS), an employee was observed improperly operating two (2) spring return (i.e., spring closed) valves identified as Safety Related Equipment (SRE). These valves were also identified as Items Relied on for Safety (IROFS). The valves were important for both Nuclear Criticality Safety (NCS) and Chemical Safety purposes. The valves were observed to be 'propped' in the open position, thereby circumventing the intended safety function of the valves. The event was reported to the NRC Operations Center on June 18, 2014 (Event no. 50208); and, NFS' investigation concluded that the cause (of the event) was due to a willful act by a Nuclear Process Operator. In addition, when reviewing the event for NRC reportability, there was an organizational focus on NCS accident sequences with limited consideration of other aspects of the Integrated Safety Analysis (ISA). In other words, when determining reportability, the performance requirements for NCS accident sequences were reviewed and determined to be satisfied; however, the chemical safety accident sequences had not been considered during this review. The initial event and the reportability issues were entered into NFS' Corrective Action Program (CAP) and investigated. The investigation also included human performance and safety culture implications. This paper summarizes the circumstances surrounding the event and discusses the lessons learned as a result of the event investigation and review. (authors)

42 ENGINEERING↗
DOE OSTI2021

Predicting drug-metagenome interactions: Variation in the microbial β-glucuronidase level in the human gut metagenomes

Characterizing the gut microbiota in terms of their capacity to interfere with drug metabolism is necessary to achieve drug efficacy and safety. Although examples of drug-microbiome interactions are well-documented, little has been reported about a computational pipeline for systematically identifying and characterizing bacterial enzymes that process particular classes of drugs. The goal of our study is to develop a computational approach that compiles drugs whose metabolism may be influenced by a particular class of microbial enzymes and that quantifies the variability in the collective level of those enzymes among individuals. The present paper describes this approach, with microbial β-glucuronidases as an example, which break down drug-glucuronide conjugates and reactivate the drugs or their metabolites. We identified 100 medications that may be metabolized by β-glucuronidases from the gut microbiome. These medications included morphine, estrogen, ibuprofen, midazolam, and their structural analogues. The analysis of metagenomic data available through the Sequence Read Archive (SRA) showed that the level of β-glucuronidase in the gut metagenomes was higher in males than in females, which provides a potential explanation for the sex-based differences in efficacy and toxicity for several drugs, reported in previous studies. Our analysis also showed that infant gut metagenomes at birth and 12 months of age have higher levels of β-glucuronidase than the metagenomes of their mothers and the implication of this observed variability was discussed in the context of breastfeeding as well as infant hyperbilirubinemia. Overall, despite important limitations discussed in this paper, our analysis provided useful insights on the role of the human gut metagenome in the variability in drug response among individuals. Importantly, this approach exploits drug and metagenome data available in public databases as well as open-source cheminformatics and bioinformatics tools to predict drug-metagenome interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

High quality draft genome sequence of Brachymonas chironomi AIMA4T (DSM 19884T) isolated from a Chironomus sp. egg mass

Brachymonas chironomi strain AIMA4T (Halpern et al., 2009) is a Gram-negative, non-motile, aerobic, chemoorganotroph bacterium. B. chironomi is a member of the Comamonadaceae, a family within the class Betaproteobacteria. This species was isolated from a chironomid (Diptera; Chironomidae) egg mass, sampled from a waste stabilization pond in northern Israel. Phylogenetic analysis based on the 16S rRNA gene sequences placed strain AIMA4T in the genus Brachymonas. Here we describe the features of this organism, together with the complete genome sequence and annotation. We find the DNA GC content is 63.5%. The chromosome length is 2,509,395 bp. It encodes 2,382 proteins and 68 RNA genes. Brachymonas chironomi genome is part of the Genomic Encyclopedia of Type Strains, Phase I: the one thousand microbial genomes (KMG) project.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Ancestral archaea expanded the genetic code with pyrrolysine

The pyrrolysyl-tRNA synthetase (PylRS) facilitates the cotranslational installation of the 22nd amino acid pyrrolysine. Owing to its tolerance for diverse amino acid substrates, and its orthogonality in multiple organisms, PylRS has emerged as a major route to install noncanonical amino acids into proteins in living cells. Recently, a novel class of PylRS enzymes was identified in a subset of methanogenic archaea. Enzymes within this class (ΔPylSn) lack the N-terminal tRNA-binding domain that is widely conserved amongst PylRS enzymes, yet remain active and orthogonal in bacteria and eukaryotes. In this study, we use biochemical and in vivo UAG-readthrough assays to characterize the aminoacylation efficiency and substrate spectrum of a ΔPylSn class PylRS from the archaeon Candidatus Methanomethylophilus alvus. We show that, compared with the full-length enzyme from Methanosarcina mazei, the Ca. M. alvus PylRS displays reduced aminoacylation efficiency but an expanded amino acid substrate spectrum. To gain insight into the evolution of ΔPylSn enzymes, we performed molecular phylogeny using 156 PylRS and 105 pyrrolysine tRNA (tRNA Pyl ) sequences from diverse archaea and bacteria. This analysis suggests that the PylRS•tRNA Pyl pair diverged before the evolution of the three domains of life, placing an early limit on the evolution of the Pyl-decoding trait. Furthermore, our results document the coevolutionary history of PylRS and tRNA Pyl and reveal the emergence of tRNA Pyl sequences with unique A73 and U73 discriminator bases. The orthogonality of these tRNA Pyl species with the more common G73-containing tRNA Pyl will enable future efforts to engineer PylRS systems for further genetic code expansion.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

On the distribution of stacking faults at dissociated medium-angle grain boundaries: Crystallographic geometry and metastability

Grain boundaries in FCC metals with low stacking-fault energy can form in dissociated configurations of stacking faults. Perhaps the most studied example is 9R stacking at boundaries near Σ3{112}, where the distribution of stacking faults is related to the emission of Shockley partial dislocations. Here, we combine atomic-resolution electron microscopy, atomistic simulations, and dislocation theory to demonstrate that boundaries vicinal to Σ33a support the stabilization of dissociated 9R stacking within a narrow range of inclinations. This boundary is interesting since its misorientation (20.05°) lies in the medium-angle regime, just past the upper misorientation limit for low-angle boundaries, motivating questions for how best to describe it in terms of dislocations. Our HAADF-STEM observations of thin film bicrystals, supported by atomistic modeling, reveal that this inclination dependence arises from specific geometric constraints on the arrangement of Shockley partial dislocations at the interface. Quantification of stacking-fault distributions across multiple boundaries indicates that the density and spacing of faults closely follow the ideal 9R motif, with subtle variations reflecting the complex energy landscape of these boundaries. Through energetic analysis, we establish the presence of competing metastable states enabled by variations in stacking sequences, emphasizing the significant role of crystallographic geometry. We generalize our analysis as a function of misorientation, showing how 9R at the Σ33a boundary is related to previous observations and calculations of HCP at a 29.7° boundary. This study provides a crystallographically grounded framework connecting dislocation structures, stacking-fault distributions, and metastability at grain boundaries in FCC metals.

Atomistic modeling↗
DOE OSTI2026

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

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

Genomics↗
DOE OSTI2022

Case Study: NREL Campus Chilled Water Storage Potential: Benchmark Datasets Development and Applications, Task 4 - Use Case Demonstration

The Benchmark Datasets Development and Applications project is a three-year collaboration between the National Renewable Energy Laboratory (NREL), Oak Ridge National Laboratory, Pacific Northwest National Laboratory, and Lawrence Berkeley National Laboratory. The project seeks to collect and curate high-resolution, well-calibrated time series of building operational and indoor/outdoor environmental data, which are crucial to understanding and optimizing building energy efficiency performance and demand flexibility capabilities as well as benchmarking energy algorithms. Project outcomes include approximately twelve high-fidelity building datasets, enhanced data representation tools, and four case studies to illustrate example applications. The goal of these case studies is to define and execute analyses that demonstrate how one or more datasets collected through this project can address a data gap or challenge historically faced by building stakeholders. This technical paper summarizes the findings of one of these case studies, in which we studied the operational efficiencies of the central cooling system at NREL. We looked at three years of data from the three chillers in the Field Test Laboratory Building (FTLB), from 2019 to 2021, to compare equipment operation and demand throughout the time period. Our analysis indicates that all three chillers are operating at or below the optimal loading conditions for most of the operation time, and thus there was no efficiency drop due to loading of the chillers at full capacity. Our recommendation is that no chiller capacity increase is needed; instead, the central plant could benefit from adopting advanced control logics for optimal sequencing of chillers during part load operations. Analysis of adding chilled water thermal storage to the central plant indicated 34% savings in demand cost and 24.5% savings in total cost (energy consumption and demand charge cost). The payback period is estimated to be 11-22 years with an assumed TES cost of $\$$100-$200 per ton. This case study shows how a selected dataset is used to solve a practical building problem - learning the operational status of its components, analyzing the effectiveness of a proposed new technique, and aiding decision-making for the building operations and maintenance team.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Availability of Critical Benchmark Experiments for the Pebble Tanker Transportation Model for Nuclear Criticality Safety Validation of TRISO Pebbles

This study addresses the need for comprehensive investigations into TRi-structural ISOtropic (TRISO) fuel pebble transportation validation. In this work, an exploratory model, the pebble tanker(PT), was developed with the aim of facilitating the validation of nuclear criticality safety calculations in the context of industrial-scale transportation of TRISO fuel. The PT model was designed to investigate the availability and applicability of critical benchmark experiments crucial for assessing the transportation of these pebbles. This work incorporated sensitivity/uncertainty (S/U) similarity studies to quantify the applicability of critical benchmark experiments and to address nuclear data uncertainties in the context of TRISO transportation. Two container models were investigated: one for the Hermes-type pebble and one for the Pebble Bed Modular Reactor (PBMR)–type pebble. The models were simplified, considering fuel, containment, and either water or air, to enable a focus on the underlying physics of applications involving TRISO fuel pebbles using the PT model. A crucial aspect under consideration was the capacity of the transport package to hold pebbles while ensuring subcriticality in the flooded state. An approach in the criticality validation process involves assessing the similarity between systems through an integral index parameter evaluation. This involves calculating a correlation coefficient (referred to as c k ) based on shared nuclear data–induced uncertainty between a benchmark experiment and the application of the PT model. To facilitate this analysis, the SCALE tools, particularly the CSAS6-Shift, TSUNAMI-3D-Shift, and TSUNAMI-IP sequences, were employed for comprehensive studies in neutronics and S/U analysis. Our findings showed that there are sufficient critical experimental benchmarks to perform this validation of the PT model in the most reactive state, i.e. when the tanker is flooded. This paper provides valuable insights into validating a transport package for Generation IV TRISO fuel pebbles.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Synthetic Biology PacBio/JAWS QC Analysis (PBJ) v3.0

This software was designed as a sequence validation tool for the assembly of synthetic constructs. It analyzes FASTQ files against a list of reference sequences, combining the results from eight sequencing libraries to generate a summary, and the files needed to view the results in the Integrative Genomics Viewer (IGV) application for manual verification. This was developed for FASTQ files generated by PacBio sequencing, but could be used on any FASTQ files that do not have paired end reads. It can be used to analyze one - eight libraries at a time, and assumes that each construct sequence in the reference will be in each pool, however, this is not a requirement. This is used to identify which libraries of pooled sequences contains a perfect match, or fixable match to the reference file. This pipeline uses many freely available open source libraries, the value added is that in our application the steps of the pipeline are defined in Workflow Description Language (WDL) and run through the Cromwell workflow engine in Docker containers, for easy distribution and set up, as well as the user friendly html summary that is generated.

Simirenko, Lisa↗
DOE OSTI2024

The auxin efflux carrier PIN1a regulates vascular patterning in cereal roots

Barley (Hordeum vulgare) is an important global cereal crop and a model in genetic studies. Despite advances in characterising barley genomic resources, few mutant studies have identified genes controlling root architecture and anatomy, which plays a critical role in capturing soil resources. Our phenotypic screening of a TILLING mutant collection identified line TM5992 exhibiting a short-root phenotype compared with wild-type (WT) Morex background. Outcrossing TM5992 with barley variety Proctor and subsequent SNP array-based bulk segregant analysis, fine mapped the mutation to a cM scale. Exome sequencing pinpointed a mutation in the candidate gene HvPIN1a, further confirming this by analysing independent mutant alleles. Detailed analysis of root growth and anatomy in Hvpin1a mutant alleles exhibited a slower growth rate, shorter apical meristem and striking vascular patterning defects compared to WT. Expression and mutant analyses of PIN1 members in the closely related cereal brachypodium (Brachypodium distachyon) revealed that BdPIN1a and BdPIN1b were redundantly expressed in root vascular tissues but only Bdpin1a mutant allele displayed root vascular defects similar to Hvpin1a. We conclude that barley PIN1 genes have sub-functionalised in cereals, compared to Arabidopsis (Arabidopsis thaliana), where PIN1a sequences control root vascular patterning.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Distributed Many-to-Many Protein Sequence Alignment using Sparse Matrices

Identifying similar protein sequences is a core step in many computational biology pipelines such as detection of homologous protein sequences, generation of similarity protein graphs for downstream analysis, functional annotation, and gene location. Performance and scalability of protein similarity search have proven to be a bottleneck in many bioinformatics pipelines due to increase in cheap and abundant sequencing data. This work presents a new distributed-memory software PASTIS. PASTIS relies on sparse matrix computations for efficient identification of possibly similar proteins. We use distributed sparse matrices for scalability and show that the sparse matrix infrastructure is a great fit for protein similarity search when coupled with a fully-distributed dictionary of sequences that allow remote sequence requests to be fulfilled. Our algorithm incorporates the unique bias in amino acid sequence substitution in search without altering basic sparse matrix model, and in turn, achieves ideal scaling up to millions of protein sequences.

97 MATHEMATICS AND COMPUTING↗