Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nucleotide”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Characterization of a 2016-2017 Human Seasonal H3 Influenza A Virus Spillover Now Endemic to U.S. Swine

In 2017, the Iowa State University Veterinary Diagnostic Laboratory detected a reverse-zoonotic transmission of a human seasonal H3 influenza A virus into swine (IAV-S) in Oklahoma. Pairwise comparison between the recently characterized human seasonal H3 IAV-S (H3.2010.2) hemagglutinin (HA) sequences detected in swine and the most similar 2016–2017 human seasonal H3 revealed 99.9% nucleotide identity. To elucidate the origin of H3.2010.2 IAV-S, 45 HA and 27 neuraminidase (NA) sequences from 2017 to 2020 as well as 11 whole-genome sequences (WGS) were genetically characterized. Time to most recent common human ancestor was estimated between August and September 2016. The N2 NA was of human origin in all but one strain from diagnostic submissions with NA sequences, and the internal gene segments from WGS consisted of matrix genes originating from the 2009 pandemic H1N1 and another 5 internal genes of triple reassortant swine origin (TTTTPT). Pigs experimentally infected with H3.2010.2 demonstrated efficient nasal shedding and replication in the lungs, mild pneumonia, and minimal microscopic lung lesions and transmitted the virus to indirect contact swine. Antigenically, H3.2010.2 viruses were closer to a human seasonal vaccine strain, A/Hong Kong/4801/2014, than to the H3.2010.1 human seasonal H3 viruses detected in swine in 2012. This was the second sustained transmission of a human seasonal IAV into swine from the 2010 decade after H3.2010.1. Monitoring the spillover and detection of novel IAV from humans to swine may help vaccine antigen selection and could impact pandemic preparedness.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Suppressor Mutations in Type II Secretion Mutants of Vibrio cholerae : Inactivation of the VesC Protease

The type II secretion system (T2SS) is a conserved transport pathway responsible for the secretion of a range of virulence factors by many pathogens, including Vibrio cholerae. Disruption of the T2SS genes in V. cholerae results in loss of secretion, changes in cell envelope function, and growth defects. While T2SS mutants are viable, high-throughput genomic analyses have listed these genes among essential genes. To investigate whether secondary mutations arise as a consequence of T2SS inactivation, we sequenced the genomes of six V. cholerae T2SS mutants with deletions or insertions in either the epsG, epsL, or epsM genes and identified secondary mutations in all mutants. Two of the six T2SS mutants contain distinct mutations in the gene encoding the T2SS-secreted protease VesC. Other mutations were found in genes coding for V. cholerae cell envelope proteins. Subsequent sequence analysis of the vesC gene in 92 additional T2SS mutant isolates identified another 19 unique mutations including insertions or deletions, sequence duplications, and single-nucleotide changes resulting in amino acid substitutions in the VesC protein. Analysis of VesC variants and the X-ray crystallographic structure of wild-type VesC suggested that all mutations lead to loss of VesC production and/or function. One possible mechanism by which V. cholerae T2SS mutagenesis can be tolerated is through selection of vesC-inactivating mutations, which may, in part, suppress cell envelope damage, establishing permissive conditions for the disruption of the T2SS. Other mutations may have been acquired in genes encoding essential cell envelope proteins to prevent proteolysis by VesC.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Predicting variable gene content in Escherichia coli using conserved genes

Having the ability to predict the protein-encoding gene content of an incomplete genome or metagenome-assembled genome is important for a variety of bioinformatic tasks. In this study, as a proof of concept, we built machine learning classifiers for predicting variable gene content in Escherichia coli genomes using only the nucleotide k-mers from a set of 100 conserved genes as features. Protein families were used to define orthologs, and a single classifier was built for predicting the presence or absence of each protein family occurring in 10%–90% of all E. coli genomes. The resulting set of 3,259 extreme gradient boosting classifiers had a per-genome average macro F1 score of 0.944 [0.943–0.945, 95% CI]. We show that the F1 scores are stable across multi-locus sequence types and that the trend can be recapitulated by sampling a smaller number of core genes or diverse input genomes. Surprisingly, the presence or absence of poorly annotated proteins, including “hypothetical proteins” was accurately predicted (F1 = 0.902 [0.898–0.906, 95% CI]). Models for proteins with horizontal gene transfer-related functions had slightly lower F1 scores but were still accurate (F1s = 0.895, 0.872, 0.824, and 0.841 for transposon, phage, plasmid, and antimicrobial resistance-related functions, respectively). Finally, using a holdout set of 419 diverse E. coli genomes that were isolated from freshwater environmental sources, we observed an average per-genome F1 score of 0.880 [0.876–0.883, 95% CI], demonstrating the extensibility of the models. Overall, this study provides a framework for predicting variable gene content using a limited amount of input sequence data.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Identification of hidden N4-like viruses and their interactions with hosts

The N4-like viruses, which were recently assigned to the novel viral family Schitoviridae in 2021, belong to a podoviral-like viral lineage and possess conserved genomic characteristics and a unique replication mechanism. Despite their significance, our understanding of N4-like viruses is primarily based on viral isolates. To address this knowledge gap, this study has established a comprehensive N4-like viral data sets comprising 342 high-quality N4-like viruses/proviruses (144 viral isolates, 158 uncultured viruses, and 40 integrated N4-like proviruses). These viruses were classified into 97 subfamilies (89 of which are newly identified), 148 genera (100 of which are newly identified), and 253 species (177 of which are newly identified). The study reveals that N4-like viruses inhibit the polar region, oligotrophic open oceans, and the human gut, where they infect various bacterial lineages, such as Alpha/Beta/Gamma/Epsilon-proteobacteria in the Proteobacteria phylum. Although N4-like viral endogenization appears to be prevalent in Proteobacteria, it has also been observed in Firmicutes. Additionally, the phylogenetic analysis has identified evolutionary divergence within the hallmark genes of N4-like viruses, indicating a complex origin of the different conserved parts of viral genomes. Moreover, 1,101 putative auxiliary metabolic genes (AMGs) were identified in the N4-like viral pan-proteome, which mainly participate in nucleotide and cofactor/vitamin metabolisms. Of these AMGs, 27 were found to be associated with virulence, suggesting their potential involvement in the spread of bacterial pathogenicity. The findings of this study are significant, as N4-like viruses represent a unique viral lineage with a distinct replication mechanism and a conserved core genome. This work has resulted in a comprehensive global map of the entire N4-like viral lineage, including information on their distribution in different biomes, evolutionary divergence, genomic diversity, and the potential for viral-mediated host metabolic reprogramming. As such, this work significantly contributes to our understanding of the ecological function and viral-host interactions of bacteriophages.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Comprehensive characterization of extracellular vesicles produced by environmental (Neff) and clinical (T4) strains of Acanthamoeba castellanii

We conducted a comprehensive comparative analysis of extracellular vesicles (EVs) from two Acanthamoeba castellanii strains, Neff (environmental) and T4 (clinical). Morphological analysis via transmission electron microscopy revealed slightly larger Neff EVs (average = 194.5 nm) compared to more polydisperse T4 EVs (average = 168.4 nm). Nanoparticle tracking analysis (NTA) and dynamic light scattering validated these differences. Proteomic analysis of the EVs identified 1,352 proteins, with 1,107 common, 161 exclusive in Neff, and 84 exclusively in T4 EVs. Gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) mapping revealed distinct molecular functions and biological processes and notably, the T4 EVs enrichment in serine proteases, aligned with its pathogenicity. Lipidomic analysis revealed a prevalence of unsaturated lipid species in Neff EVs, particularly triacylglycerols, phosphatidylethanolamines (PEs), and phosphatidylserine, while T4 EVs were enriched in diacylglycerols and diacylglyceryl trimethylhomoserine, phosphatidylcholine and less unsaturated PEs, suggesting differences in lipid metabolism and membrane permeability. Metabolomic analysis indicated Neff EVs enrichment in glycerolipid metabolism, glycolysis, and nucleotide synthesis, while T4 EVs, methionine metabolism. Furthermore, RNA-seq of EVs revealed differential transcript between the strains, with Neff EVs enriched in transcripts related to gluconeogenesis and translation, suggesting gene regulation and metabolic shift, while in the T4 EVs transcripts were associated with signal transduction and protein kinase activity, indicating rapid responses to environmental changes. In this novel study, data integration highlighted the differences in enzyme profiles, metabolic processes, and potential origins of EVs in the two strains shedding light on the diversity and complexity of A. castellanii EVs and having implications for understanding host-pathogen interactions and developing targeted interventions for Acanthamoeba-related diseases.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Distinguishing Leptothrix and Sphaerotilus genera by an integrated genomic-phenotypic analysis supported by new Leptothrix genomes

The Sphaerotilus-Leptothrix group of bacteria includes one of the first described microorganisms, Leptothrix ochracea, an uncultured type strain, plus isolates of Leptothrix and Sphaerotilus. This group is unified by the ability to form sheaths and oxidize metals, although L. ochracea exhibits obvious ecological, morphological, and functional differences from the rest of Sphaerotilus-Leptothrix. Recently, there have been calls to combine the group into one genus, Sphaerotilus; however, these studies lacked adequate genomic representation of L. ochracea. Here, we present a comprehensive comparative genomic analysis of the Sphaerotilus-Leptothrix group, including expanded representation of L. ochracea, a closely related novel species, Leptothrix toolikensis, and two new isolates (Leptothrix mechoopdaensis). Analysis of 38 genomes resolves three phylogenetic and functional groups: the ochracea-type Leptothrix (Group 1), the mobilis-type Leptothrix (Group 2), and Sphaerotilus (Group 3). Group 1 genomes form a separate genus based on average nucleotide identity and alignment fraction. The genomes clearly diverge from the rest of Sphaerotilus-Leptothrix in phylogeny, size, and metabolic potential. Group 1 genomes are much smaller (2.59–3.04 Mb) than those of Groups 2 (4.55–6.06 Mb) and 3 (3.94–5.07 Mb), while encoding more metal oxidases and fewer carbohydrate-active enzymes. Group 2 clusters with Group 3 phylogenetically and is similar in organic carbon metabolisms but maintains more metal oxidation genes. Group 2 members lack homogeneity in phenotype and genotype, suggesting that additional isolates and genomes are needed for confident classification. However, Group 1 genomes (L. ochracea and L. toolikensis) show clear divergence, precluding their inclusion in Sphaerotilus and supporting the retention of the genus Leptothrix.

Leptothrix↗
DOE OSTI2023

Application of quasimetagenomics methods to define microbial diversity and subtype Listeria monocytogenes in dairy and seafood production facilities

Microorganisms frequently colonize surfaces within food production facilities. Detection of Listeria monocytogenes in this setting relies on culture-dependent methods, but the complex dynamics of bacterial interactions within these environments and their impact on pathogen detection remain largely unexplored. To address this challenge, we applied both 16S rRNA and shotgun quasimetagenomic (enriched microbiome) sequencing of swab culture enrichments from five seafood and seven dairy production environments. Utilizing 16S rRNA amplicon sequencing, we observed variability between 355 samples taken from these 12 production facilities and a distinctive microbiome for each environment. With shotgun quasimetagenomic sequencing, we were able to assemble L. monocytogenes metagenome-assembled genomes (MAGs) from 28 of the 32 culture-positive samples. We compared these MAGs to their corresponding whole-genome sequencing assemblies, which resulted in two polyphyletic clades consisting of L. monocytogenes lineages I and II with 13,195 and 25,556 single-nucleotide polymorphism sites, respectively. The remaining four MAGs did not produce sufficient genome coverage. To understand and establish limits for pathogen detection and subtyping using shotgun quasimetagenomics, these same data sets were downsampled in slilico to produce a titration series of abundances of L. monocytogenes and analyzed. Pathogen detection was achieved for all downsampled data sets, even those with only 3× genome coverage. This study contributes to the understanding of microbial diversity within food production environments and presents insights into the level of genome coverage needed in a metagenome sequencing data set to detect, subtype, and source track a foodborne pathogen.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

SourceFinder: a Machine-Learning-Based Tool for Identification of Chromosomal, Plasmid, and Bacteriophage Sequences from Assemblies

High-throughput genome sequencing technologies enable the investigation of complex genetic interactions, including the horizontal gene transfer of plasmids and bacteriophages. However, identifying these elements from assembled reads remains challenging due to genome sequence plasticity and the difficulty in assembling complete sequences. In this study, we developed a classifier, using random forest, to identify whether sequences originated from bacterial chromosomes, plasmids, or bacteriophages. The classifier was trained on a diverse collection of 23,211 chromosomal, plasmid, and bacteriophage sequences from hundreds of bacterial species. In order to adapt the classifier to incomplete sequences, each complete sequence was subsampled into 5,000 nucleotide fragments and further subdivided into k-mers. This three-class classifier succeeded in identifying chromosomes, plasmids, and bacteriophages using k-mer distributions of complete and partial genome sequences, including simulated metagenomic scaffolds with minimum performance of 0.939 area under the receiver operating characteristic curve (AUC). This classifier, implemented as SourceFinder, has been made available as an online web service to help the community with predicting the chromosomal, plasmid, and bacteriophage sources of assembled bacterial sequence data (https://cge.food.dtu.dk/services/SourceFinder/).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Cyanobacteriochromes: A Rainbow of Photoreceptors

Widespread phytochrome photoreceptors use photoisomerization of linear tetrapyrrole (bilin) chromophores to measure the ratio of red to far-red light. Cyanobacteria also contain distantly related cyanobacteriochrome (CBCR) proteins that share the bilin-binding GAF domain of phytochromes but sense other colors of light. CBCR photocycles are extremely diverse, ranging from the near-UV to the near-IR. Photoisomerization of the bilin triggers photoconversion of the CBCR input, thereby modulating the biochemical signaling state of output domains such as histidine kinase bidomains that can interface with cellular signal transduction pathways. CBCRs thus can regulate several aspects of cyanobacterial photobiology, including phototaxis, metabolism of cyclic nucleotide second messengers, and optimization of the cyanobacterial light-harvesting apparatus. This review examines spectral tuning, photoconversion, and photobiology of CBCRs and recent developments in understanding their evolution and in applying them in synthetic biology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Exploring the Genetic Basis of Wild Boar ( Sus scrofa ) and Its Connection to Classical Swine Fever Spread

Classical swine fever (CSF) is the one of the most devastating contagious diseases in domestic swine and wild boar/pigs (Sus scrofa). Population genetics is often used to estimate animal dispersal and can also help evaluate host population connectivity, which is crucial for understanding pathogen dispersal. We surveyed genetic population structure of boars using MIG-seq analysis to clarify the geographic barriers that influence boar dispersal in north-central Japan and to demonstrate the relationship between the spread of CSF infection among boars and their population structure. We obtained 382 single-nucleotide polymorphisms from 348 wild boar samples, and the results of STRUCTURE analysis indicated that the highest ΔK value was at K = 2, followed by K = 4. Based on these results, it is evident that the Abukuma river, a major river within north-central Japan, does not act as a barrier to the gene flow of boars, but rather that human infrastructure hinders their dispersal. Further, according to the time series change in the capture site of CSF-infected wild boar and the sum of the probability of belonging to each of the four clades in individual CSF-infected wild boar, our results indicated that the genetic structure of boar populations was correlated with the outbreak pathway of CSF across our study region. Our study suggests that predictions of disease spread, especially for widely distributed host species, is challenging because of the risk of cryptic breaks and changes in wide range connectivity; however, understanding the genetic population structure of wild boar can be a useful tool for predicting the spread of CSF. We concluded that genetic analysis of host population structure may have the possibility to improve predictions of the future dynamics of disease spread.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

FutureTense

Protective vaccines and reliable diagnostics are essential tools for controlling viral diseases. However, the efficacy of these tools can be diminished by mutations in viral genomes. The delay between the emergence of new viral strains and the redesign of vaccines and diagnostics allows for continued viral transmission. Is it possible to address this challenge by computationally predicting viral genome sequence evolution? Can we “future-proof” vaccines and diagnostics by targeting both current and anticipated future sequence variants? While predicting viral evolution is still an unsolved, “grand challenge” problem in biology, the large, and rapidly growing, number of SARS-CoV-2 genome sequences provide an opportunity to quantify the ability of machine learning to predict viral genome sequence evolution. Towards this end, we have developed a simple computational model for predicting viral evolution at the level of individual nucleotides. The key metric for quantifying the per-base, prediction accuracy for viral evolution is the Mann-Whitney U statistic (or, equivalently, the area under the receiver operator curve). Since the Mann-Whitney U statistic is not a differentiable function, existing deep leaning packages (like Pytorch and Keras/TensorFlow) are not useful, as they require that the accuracy metric/objective function be analytically differentiable with respect to the model parameters. To overcome this challenge, we have implemented custom software, “FutureTense”, that can train a machine learning model by maximizing the non-differentiable Mann-Whitney U statistic. This software trains a machine learning model by exploring along the direction of the discrete gradient of the Mann-Whitney U statistic in the model parameter space. Parallel computing and genome sequence-specific optimizations are used to accelerate model training. The resulting machine learning model learns the observed high C->U mutation rates in the SARS-CoV-2 genome (which are potentially induced by host defenses) and provides prediction accuracies that are significantly better than one would expect from random chance. While predicting viral evolution is still quite far from a solved problem, the surprising performance of this simple model gives hope that the accuracy of predicting viral genome evolution can be further increased by more sophisticated approaches.

Gans, Jason↗
DOE OSTI2021

jialiu232/MetaFunPrimer_paper_info

Genes belonging to the same functional group may include numerous and variable gene sequences, making characterizing and quantifying difficult. Therefore, high-throughput design tools are needed to simultaneously create primers for improved quantification of target genes. We developed MetaFunPrimer, a bioinformatic pipeline, to design primers for numerous genes of interest. This tool also enables gene target prioritization based on ranking the presence of genes in user-defined references, such as environment-specific metagenomes. Given inputs of protein and nucleotide sequences for gene targets of interest and an accompanying set of reference metagenomes or genomes, MetaFunPrimer generates primers for ranked genes of interest. To demonstrate the usage and benefits of MetaFunPrimer, a total of 78 primer pairs were designed to target observed ammonia monooxygenase subunit A (amoA) genes of ammonia-oxidizing bacteria (AOB) in 1,550 publicly available soil metagenomes. We demonstrate computationally that these amoA-AOB primers can cover 94% of the amoA-AOB genes observed in the 1,550 soil metagenomes compared with a 49% estimated coverage by previously published primers. Finally, we verified the utility of these primer sets in incubation experiments that used long-term nitrogen fertilized or unfertilized soils. High-throughput quantitative PCR (qPCR) results and statistical analyses showed significant differences in relative quantification patterns between the two soils, and subsequent absolute quantifications also confirmed that target genes enumerated by six selected primer pairs were significantly more abundant in the nitrogen-fertilized soils. This new tool gives microbial ecologists a new approach to assess functional gene abundance and related microbial community dynamics quickly and affordably.

Liu, Jia↗
DOE OSTI2024

NuclPred v1

This tool takes a genome assembly as input and predicts per-site nucleosome occupancy as output. Trained on physical maps of nucleosome binding preferences across the fungal kingdom, NuclPred can be applied broadly across fungi (and other eukaryotes). This breadth, combined with its accuracy, means it could have both basic and applied biological implications, for example in understanding eukaryotic gene regulation and genetic engineering. Almost universally across eukaryotes, nucleosomes - each wrapping ~150 base pairs of DNA - serve to package DNA inside the nucleus, with major consequences on DNA access, gene activity and DNA integration. NuclPred was generated using a supervised deep learning approach combining convolutional and recurrent neural networks to take DNA features (nucleotides, GC content and structural information) as input, then use that information to predict the physical attractiveness DNA sequences might have for forming nucleosomes. With this information at hand, researchers can design more efficient CRISPR constructs, explore the interplay between DNA signatures and other regulators impact nucleosome locations, predict expression patterns, etc. This tool will be published as part of a manuscript currently under revision at iScience (draft attached).

Mondo, Stephen↗
DOE OSTI2026

Genomic Language model for Annotation of Repetitive Elements (GLARE) v1.0

GLARE (Genomic Language model for Annotation of Repetitive Elements) is a tool that classifies transposable elements (TEs)—the mobile, repetitive DNA sequences that make up large fractions of eukaryotic genomes. GLARE fine-tunes the NTv3-650M genomic language model on a harmonized collection of curated TE sequences from the PanTEon and Repbase reference databases, assigning each input sequence to one of 11 orders and 32 superfamilies in a Wicker-compatible taxonomy. Features. From nucleotide FASTA input, GLARE outputs per-sequence predictions, class summaries, composition figures, and an annotated FASTA. It provides calibrated confidence scores with optional abstention and runs on CPU or GPU. Uses. GLARE serves as a classification component in genome-annotation pipelines, downstream of TE discovery, supporting genome annotation and comparative and evolutionary genomics. Advantages. GLARE is the first repeat-element classifier to leverage a pretrained genomic language model. Combined with multi-database training, this approach outperformed all nine classifiers in the PanTEon benchmark, generalized better to unseen taxonomic clades, and remained robust to sequence orientation—a common failure mode of existing tools.

Bruna, Tomas [Lawrence Berkeley National Laborator↗
DOE OSTI2024

An FDA-approved drug structurally and phenotypically corrects the K210del mutation in genetic cardiomyopathy models

Dilated cardiomyopathy (DCM) due to genetic disorders results in decreased myocardial contractility, leading to high morbidity and mortality rates. There are several therapeutic challenges in treating DCM, including poor understanding of the underlying mechanism of impaired myocardial contractility and the difficulty of developing targeted therapies to reverse mutation-specific pathologies. In this report, we focused on K210del, a DCM-causing mutation, due to 3-nucleotide deletion of sarcomeric troponin T (TnnT), resulting in loss of Lysine210. We resolved the crystal structure of the troponin complex carrying the K210del mutation. K210del induced an allosteric shift in the troponin complex resulting in distortion of activation Ca 2+ -binding domain of troponin C (TnnC) at S69, resulting in calcium discoordination. Next, we adopted a structure-based drug repurposing approach to identify bisphosphonate risedronate as a potential structural corrector for the mutant troponin complex. Cocrystallization of risedronate with the mutant troponin complex restored the normal configuration of S69 and calcium coordination. Risedronate normalized force generation in K210del patient-induced pluripotent stem cell–derived (iPSC-derived) cardiomyocytes and improved calcium sensitivity in skinned papillary muscles isolated from K210del mice. Systemic administration of risedronate to K210del mice normalized left ventricular ejection fraction. Collectively, these results identify the structural basis for decreased calcium sensitivity in K210del and highlight structural and phenotypic correction as a potential therapeutic strategy in genetic cardiomyopathies.

Research & Experimental Medicine↗
DOE OSTI2020

ADEPT: a domain independent sequence alignment strategy for gpu architectures

Bioinformatic workflows frequently make use of automated genome assembly and protein clustering tools. At the core of most of these tools, a significant portion of execution time is spent in determining optimal local alignment between two sequences. This task is performed with the Smith-Waterman algorithm, which is a dynamic programming based method. With the advent of modern sequencing technologies and increasing size of both genome and protein databases, a need for faster Smith-Waterman implementations has emerged. Multiple SIMD strategies for the Smith-Waterman algorithm are available for CPUs. However, with the move of HPC facilities towards accelerator based architectures, a need for an efficient GPU accelerated strategy has emerged. Existing GPU based strategies have either been optimized for a specific type of characters (Nucleotides or Amino Acids) or for only a handful of application use-cases. In this paper, we present ADEPT, a new sequence alignment strategy for GPU architectures that is domain independent, supporting alignment of sequences from both genomes and proteins. Our proposed strategy uses GPU specific optimizations that do not rely on the nature of sequence. We demonstrate the feasibility of this strategy by implementing the Smith-Waterman algorithm and comparing it to similar CPU strategies as well as the fastest known GPU methods for each domain. ADEPT’s driver enables it to scale across multiple GPUs and allows easy integration into software pipelines which utilize large scale computational systems. We have shown that the ADEPT based Smith-Waterman algorithm demonstrates a peak performance of 360 GCUPS and 497 GCUPs for protein based and DNA based datasets respectively on a single GPU node (8 GPUs) of the Cori Supercomputer. Overall ADEPT shows 10x faster performance in a node-to-node comparison against a corresponding SIMD CPU implementation. ADEPT demonstrates a performance that is either comparable or better than existing GPU strategies. We demonstrated the efficacy of ADEPT in supporting existing bionformatics software pipelines by integrating ADEPT in MetaHipMer a high-performance denovo metagenome assembler and PASTIS a high-performance protein similarity graph construction pipeline. Our results show 10% and 30% boost of performance in MetaHipMer and PASTIS respectively.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

G4Boost: a machine learning-based tool for quadruplex identification and stability prediction

Background: G-quadruplexes (G4s), formed within guanine-rich nucleic acids, are secondary structures involved in important biological processes. Although every G4 motif has the potential to form a stable G4 structure, not every G4 motif would, and accurate energy-based methods are needed to assess their structural stability. Here, we present a decision tree-based prediction tool, G4Boost, to identify G4 motifs and predict their secondary structure folding probability and thermodynamic stability based on their sequences, nucleotide compositions, and estimated structural topologies. Results: G4Boost predicted the quadruplex folding state with an accuracy greater then 93% and an F1-score of 0.96, and the folding energy with an RMSE of 4.28 and R 2 of 0.95 only by the means of sequence intrinsic feature. G4Boost was successfully applied and validated to predict the stability of experimentally-determined G4 structures, including for plants and humans. Conclusion: G4Boost outperformed the three machine-learning based prediction tools, DeepG4, Quadron, and G4RNA Screener, in terms of both accuracy and F1-score, and can be highly useful for G4 prediction to understand gene regulation across species including plants and humans.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Profiling the BLAST bioinformatics application for load balancing on high-performance computing clusters

Abstract Background The Basic Local Alignment Search Tool (BLAST) is a suite of commonly used algorithms for identifying matches between biological sequences. The user supplies a database file and query file of sequences for BLAST to find identical sequences between the two. The typical millions of database and query sequences make BLAST computationally challenging but also well suited for parallelization on high-performance computing clusters. The efficacy of parallelization depends on the data partitioning, where the optimal data partitioning relies on an accurate performance model. In previous studies, a BLAST job was sped up by 27 times by partitioning the database and query among thousands of processor nodes. However, the optimality of the partitioning method was not studied. Unlike BLAST performance models proposed in the literature that usually have problem size and hardware configuration as the only variables, the execution time of a BLAST job is a function of database size, query size, and hardware capability. In this work, the nucleotide BLAST application BLASTN was profiled using three methods: shell-level profiling with the Unix “time” command, code-level profiling with the built-in “profiler” module, and system-level profiling with the Unix “gprof” program. The runtimes were measured for six node types, using six different database files and 15 query files, on a heterogeneous HPC cluster with 500+ nodes. The empirical measurement data were fitted with quadratic functions to develop performance models that were used to guide the data parallelization for BLASTN jobs. Results Profiling results showed that BLASTN contains more than 34,500 different functions, but a single function, RunMTBySplitDB, takes 99.12% of the total runtime. Among its 53 child functions, five core functions were identified to make up 92.12% of the overall BLASTN runtime. Based on the performance models, static load balancing algorithms can be applied to the BLASTN input data to minimize the runtime of the longest job on an HPC cluster. Four test cases being run on homogeneous and heterogeneous clusters were tested. Experiment results showed that the runtime can be reduced by 81% on a homogeneous cluster and by 20% on a heterogeneous cluster by re-distributing the workload. Discussion Optimal data partitioning can improve BLASTN’s overall runtime 5.4-fold in comparison with dividing the database and query into the same number of fragments. The proposed methodology can be used in the other applications in the BLAST+ suite or any other application as long as source code is available.

59 BASIC BIOLOGICAL SCIENCES↗