Real-time PCR assays that detect genes for botulinum neurotoxin A-G subtypes.

The role of Real-Time PCR assays for surveillance and rapid screening for pathogens is garnering more and more attention because of its versatility and ease of adoption. The goal of this study was to design, test, and evaluate Real-Time TaqMan PCR assays for the detection of botulinum neurotoxin (bont/A-G) genes from currently recognized BoNT subtypes. Assays were computationally designed and then laboratory tested for sensitivity and specificity using DNA preparations containing bont genes from 82 target toxin subtypes, including nine bivalent toxin types; 31 strains representing other clostridial species; and an extensive panel that consisted of DNA from a diverse set of prokaryotic (bacterial) and eukaryotic (fungal, protozoan, plant, and animal) species. In addition to laboratory testing, the assays were computationally evaluated using in silico analysis for their ability to detect bont gene sequences from recently identified toxin subtypes. Seventeen specific assays (two for each of the bont/C, bont/D, bont/E, and bont/G subtypes and three for each of the bont/A, bont/B, and bont/F subtypes) were designed and evaluated for their ability to detect bont genes encoding multiple subtypes from all seven serotypes. These assays could provide an additional tool for the detection of botulinum neurotoxins in clinical, environmental and food samples that can complement other existing methods used in clinical diagnostics, regulatory, public health, and research laboratories.

Fast Evaluation of Viral Emerging Risks (FEVER): A computational tool for biosurveillance, diagnostics, and mutation typing of emerging viral pathogens.

Viral pathogens can rapidly evolve, adapt to novel hosts, and evade human immunity. The early detection of emerging viral pathogens through biosurveillance coupled with rapid and accurate diagnostics are required to mitigate global pandemics. However, RNA viruses can mutate rapidly, hampering biosurveillance and diagnostic efforts. Here, we present a novel computational approach called FEVER (Fast Evaluation of Viral Emerging Risks) to design assays that simultaneously accomplish: 1) broad-coverage biosurveillance of an entire group of viruses, 2) accurate diagnosis of an outbreak strain, and 3) mutation typing to detect variants of public health importance. We demonstrate the application of FEVER to generate assays to simultaneously 1) detect sarbecoviruses for biosurveillance; 2) diagnose infections specifically caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2); and 3) perform rapid mutation typing of the D614G SARS-CoV-2 spike variant associated with increased pathogen transmissibility. These FEVER assays had a high in silico recall (predicted positive) up to 99.7% of 525,708 SARS-CoV-2 sequences analyzed and displayed sensitivities and specificities as high as 92.4% and 100% respectively when validated in 100 clinical samples. The D614G SARS-CoV-2 spike mutation PCR test was able to identify the single nucleotide identity at position 23,403 in the viral genome of 96.6% SARS-CoV-2 positive samples without the need for sequencing. This study demonstrates the utility of FEVER to design assays for biosurveillance, diagnostics, and mutation typing to rapidly detect, track, and mitigate future outbreaks and pandemics caused by emerging viruses.

A cross-study analysis of drug response prediction in cancer cell lines.

To enable personalized cancer treatment, machine learning models have been developed to predict drug response as a function of tumor and drug features. However, most algorithm development efforts have relied on cross-validation within a single study to assess model accuracy. While an essential first step, cross-validation within a biological data set typically provides an overly optimistic estimate of the prediction performance on independent test sets. To provide a more rigorous assessment of model generalizability between different studies, we use machine learning to analyze five publicly available cell line-based data sets: National Cancer Institute 60, ancer Therapeutics Response Portal (CTRP), Genomics of Drug Sensitivity in Cancer, Cancer Cell Line Encyclopedia and Genentech Cell Line Screening Initiative (gCSI). Based on observed experimental variability across studies, we explore estimates of prediction upper bounds. We report performance results of a variety of machine learning models, with a multitasking deep neural network achieving the best cross-study generalizability. By multiple measures, models trained on CTRP yield the most accurate predictions on the remaining testing data, and gCSI is the most predictable among the cell line data sets included in this study. With these experiments and further simulations on partial data, two lessons emerge: (1) differences in viability assays can limit model generalizability across studies and (2) drug diversity, more than tumor diversity, is crucial for raising model generalizability in preclinical screening.

Optical Biosensor Platforms Display Varying Sensitivity for the Direct Detection of Influenza RNA.

Detection methods that do not require nucleic acid amplification are advantageous for viral diagnostics due to their rapid results. These platforms could provide information for both accurate diagnoses and pandemic surveillance. Influenza virus is prone to pandemic-inducing genetic mutations, so there is a need to apply these detection platforms to influenza diagnostics. Here, we analyzed the Fast Evaluation of Viral Emerging Risks (FEVER) pipeline on ultrasensitive detection platforms, including a waveguide-based optical biosensor and a flow cytometry bead-based assay. The pipeline was also evaluated in silico for sequence coverage in comparison to the U.S. Centers for Disease Control and Prevention's (CDC) influenza A and B diagnostic assays. The influenza FEVER probe design had a higher tolerance for mismatched bases than the CDC's probes, and the FEVER probes altogether had a higher detection rate for influenza isolate sequences from GenBank. When formatted for use as molecular beacons, the FEVER probes detected influenza RNA as low as 50 nM on the waveguide-based optical biosensor and 1 nM on the flow cytometer. In addition to molecular beacons, which have an inherently high background signal we also developed an exonuclease selection method that could detect 500 pM of RNA. The combination of high-coverage probes developed using the FEVER pipeline coupled with ultrasensitive optical biosensors is a promising approach for future influenza diagnostic and biosurveillance applications.

A public website for the automated assessment and validation of SARS-CoV-2 diagnostic PCR assays.

SUMMARY: Polymerase chain reaction-based assays are the current gold standard for detecting and diagnosing SARS-CoV-2. However, as SARS-CoV-2 mutates, we need to constantly assess whether existing PCR-based assays will continue to detect all known viral strains. To enable the continuous monitoring of SARS-CoV-2 assays, we have developed a web-based assay validation algorithm that checks existing PCR-based assays against the ever-expanding genome databases for SARS-CoV-2 using both thermodynamic and edit-distance metrics. The assay-screening results are displayed as a heatmap, showing the number of mismatches between each detection and each SARS-CoV-2 genome sequence. Using a mismatch threshold to define detection failure, assay performance is summarized with the true-positive rate (recall) to simplify assay comparisons. AVAILABILITY AND IMPLEMENTATION: The assay evaluation website and supporting software are Open Source and freely available at https://covid19.edgebioinformatics.org/#/assayValidation, https://github.com/jgans/thermonucleotide BLAST and https://github.com/LANL-Bioinformatics/assay_validation. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Effects of initial microbial biomass abundance on respiration during pine litter decomposition.

IMPORTANCE: Microbial biomass is one of the most common microbial parameters used in land carbon (C) cycle models, however, it is notoriously difficult to measure accurately. To understand the consequences of mismeasurement, as well as the broader importance of microbial biomass abundance as a direct driver of ecological phenomena, greater quantitative understanding of the role of microbial biomass abundance in environmental processes is needed. Using microcosms, we manipulated the initial biomass of numerous microbial communities across a 100-fold range and measured effects on CO2 production during plant litter decomposition. We found that the effects of initial biomass abundance on CO2 production was largely attenuated within a week, while the effects of community type remained significant over the course of the experiment. Overall, our results suggest that initial microbial biomass abundance in litter decomposition within an ecosystem is a weak driver of long-term C cycling dynamics.

Defending Our Public Biological Databases as a Global Critical Infrastructure.

Progress in modern biology is being driven, in part, by the large amounts of freely available data in public resources such as the International Nucleotide Sequence Database Collaboration (INSDC), the world's primary database of biological sequence (and related) information. INSDC and similar databases have dramatically increased the pace of fundamental biological discovery and enabled a host of innovative therapeutic, diagnostic, and forensic applications. However, as high-value, openly shared resources with a high degree of assumed trust, these repositories share compelling similarities to the early days of the Internet. Consequently, as public biological databases continue to increase in size and importance, we expect that they will face the same threats as undefended cyberspace. There is a unique opportunity, before a significant breach and loss of trust occurs, to ensure they evolve with quality and security as a design philosophy rather than costly "retrofitted" mitigations. This Perspective surveys some potential quality assurance and security weaknesses in existing open genomic and proteomic repositories, describes methods to mitigate the likelihood of both intentional and unintentional errors, and offers recommendations for risk mitigation based on lessons learned from cybersecurity.

Evaluation of DNA extraction methods to detect bacterial targets in aerosol samples.

DNA-based monitoring of pathogens in aerosol samples requires extraction methods that provide high recovery of DNA. To identify a suitable method, we evaluated six DNA extraction methods for recovery of target-specific DNA from samples with four bacterial agents at low abundance (<10,000 genome copies per detection assay). These methods differed in rigor of cell disruption, approach for DNA capture, and extent of DNA purification. The six methods varied 1000-fold in the recovery of DNA from spores or cells of surrogates of Bacillus anthracis, Yersinia pestis, Burkholderia pseudomallei, and Francisella tularensis, each at about 105 CFU per sample. A custom method using paramagnetic Dynabeads for DNA capture greatly outperformed the other five methods. The cDynabead method provided about 80% recovery of target-specific DNA. The cDynabead method and a filtration method were further evaluated for DNA recovery from bacterial agents spiked on filters (c.a. 105 CFU of each agent per filter quadrant) that were subsequently used to collect background outdoor air particulates for 24-h. The filtration method generally failed to recover detectable quantities of target DNA from the spiked filters, suggesting at least a 100-fold loss of target DNA during extraction, whereas the custom cDynabead method consistently yielded DNA sufficient for target detection.

A robust PCR primer design platform applied to the detection of Acidobacteria Group 1 in soil.

Environmental biosurveillance and microbial ecology studies use PCR-based assays to detect and quantify microbial taxa and gene sequences within a complex background of microorganisms. However, the fragmentary nature and growing quantity of DNA-sequence data make group-specific assay design challenging. We solved this problem by developing a software platform that enables PCR-assay design at an unprecedented scale. As a demonstration, we developed quantitative PCR assays for a globally widespread, ecologically important bacterial group in soil, Acidobacteria Group 1. A total of 33,684 Acidobacteria 16S rRNA gene sequences were used for assay design. Following 1 week of computation on a 376-core cluster, 83 assays were obtained. We validated the specificity of the top three assays, collectively predicted to detect 42% of the Acidobacteria Group 1 sequences, by PCR amplification and sequencing of DNA from soil. Based on previous analyses of 16S rRNA gene sequencing, Acidobacteria Group 1 species were expected to decrease in response to elevated atmospheric CO(2). Quantitative PCR results, using the Acidobacteria Group 1-specific PCR assays, confirmed the expected decrease and provided higher statistical confidence than the 16S rRNA gene-sequencing data. These results demonstrate a powerful capacity to address previously intractable assay design challenges.

Simultaneous pathogen detection and antibiotic resistance characterization using SNP-based multiplexed oligonucleotide ligation-PCR (MOL-PCR).

Extensive use of antibiotics in both public health and animal husbandry has resulted in rapid emergence of antibiotic resistance in almost all human pathogens, including biothreat pathogens. Antibiotic resistance has thus become a major concern for both public health and national security. We developed multiplexed assays for rapid, simultaneous pathogen detection and characterization of ciprofloxacin and doxycycline resistance in Bacillus anthracis, Yersinia pestis, and Francisella tularensis. These assays are SNP-based and use Multiplexed Oligonucleotide Ligation-PCR (MOL-PCR). The MOL-PCR assay chemistry and MOLigo probe design process are presented. A web-based tool - MOLigoDesigner (http://MOLigoDesigner.lanl.gov) was developed to facilitate the probe design. All probes were experimentally validated individually and in multiplexed assays, and minimal sets of multiplexed MOLigo probes were identified for simultaneous pathogen detection and antibiotic resistance characterization.

A rapid multiplex assay for nucleic acid-based diagnostics.

We have developed a rapid (under 4 hours), multiplex, nucleic acid assay, adapted to a microsphere array detection platform. We call this assay multiplex oligonucleotide ligation-PCR (MOL-PCR). Unlike other ligation-based assays that require multiple steps, our protocol consists of a single tube reaction, followed by hybridization to a Luminex microsphere array for detection. We demonstrate the ability of this assay to simultaneously detect diverse nucleic acid signatures (e.g., unique sequences, single nucleotide polymorphisms) in a single multiplex reaction. Detection probes consist of modular components that enable target detection, probe amplification, and subsequent capture onto microsphere arrays. To demonstrate the utility of our assay, we applied it to the detection of three biothreat agents, B. anthracis, Y. pestis, and F. tularensis. Combined with the ease and robustness of this assay, the results presented here show a strong potential of our assay for use in diagnostics and surveillance.

Steps toward broad-spectrum therapeutics: discovering virulence-associated genes present in diverse human pathogens.

BACKGROUND: New and improved antimicrobial countermeasures are urgently needed to counteract increased resistance to existing antimicrobial treatments and to combat currently untreatable or new emerging infectious diseases. We demonstrate that computational comparative genomics, together with experimental screening, can identify potential generic (i.e., conserved across multiple pathogen species) and novel virulence-associated genes that may serve as targets for broad-spectrum countermeasures. RESULTS: Using phylogenetic profiles of protein clusters from completed microbial genome sequences, we identified seventeen protein candidates that are common to diverse human pathogens and absent or uncommon in non-pathogens. Mutants of 13 of these candidates were successfully generated in Yersinia pseudotuberculosis and the potential role of the proteins in virulence was assayed in an animal model. Six candidate proteins are suggested to be involved in the virulence of Y. pseudotuberculosis, none of which have previously been implicated in the virulence of Y. pseudotuberculosis and three have no record of involvement in the virulence of any bacteria. CONCLUSION: This work demonstrates a strategy for the identification of potential virulence factors that are conserved across a number of human pathogenic bacterial species, confirming the usefulness of this tool.

Improved assay-dependent searching of nucleic acid sequence databases.

Nucleic acid-based biochemical assays are crucial to modern biology. Key applications, such as detection of bacterial, viral and fungal pathogens, require detailed knowledge of assay sensitivity and specificity to obtain reliable results. Improved methods to predict assay performance are needed for exploiting the exponentially growing amount of DNA sequence data and for reducing the experimental effort required to develop robust detection assays. Toward this goal, we present an algorithm for the calculation of sequence similarity based on DNA thermodynamics. In our approach, search queries consist of one to three oligonucleotide sequences representing either a hybridization probe, a pair of Padlock probes or a pair of PCR primers with an optional TaqMantrade mark probe (i.e. in silico or 'virtual' PCR). Matches are reported if the query and target satisfy both the thermodynamics of the assay (binding at a specified hybridization temperature and/or change in free energy) and the relevant biological constraints (assay sequences binding to the correct target duplex strands in the required orientations). The sensitivity and specificity of our method is evaluated by comparing predicted to known sequence tagged sites in the human genome. Free energy is shown to be a more sensitive and specific match criterion than hybridization temperature.

Genomorama: genome visualization and analysis.

BACKGROUND: The ability to visualize genomic features and design experimental assays that can target specific regions of a genome is essential for modern biology. To assist in these tasks, we present Genomorama, a software program for interactively displaying multiple genomes and identifying potential DNA hybridization sites for assay design. RESULTS: Useful features of Genomorama include genome search by DNA hybridization (probe binding and PCR amplification), efficient multi-scale display and manipulation of multiple genomes, support for many genome file types and the ability to search for and retrieve data from the National Center for Biotechnology Information (NCBI) Entrez server. CONCLUSION: Genomorama provides an efficient computational platform for visualizing and analyzing multiple genomes.

Computational improvements reveal great bacterial diversity and high metal toxicity in soil.

The complexity of soil bacterial communities has thus far confounded effective measurement. However, with improved analytical methods, we show that the abundance distribution and total diversity can be deciphered. Reanalysis of reassociation kinetics for bacterial community DNA from pristine and metal-polluted soils showed that a power law best described the abundance distributions. More than one million distinct genomes occurred in the pristine soil, exceeding previous estimates by two orders of magnitude. Metal pollution reduced diversity more than 99.9%, revealing the highly toxic effect of metal contamination, especially for rare taxa.

APART: automated preprocessing for NMR assignments with reduced tedium.

MOTIVATION: High-throughput NMR structure determination is a goal that will require progress on many fronts, one of which is rapid resonance assignment. An important rate-limiting step in the resonance assignment process is accurate identification of resonance peaks in the NMR spectra. Peak-picking schemes range from incomplete (which lose essential assignment connectivities) to noisy (which obscure true connectivities with many false ones). We introduce an automated preassignment process that removes false peaks from noisy peak lists by requiring consensus between multiple NMR experiments and exploiting a priori information about NMR spectra. This process is designed to accept multiple input formats and generate multiple output formats, in an effort to be compatible with a variety of user preferences. RESULTS: Automated preprocessing with APART rapidly identifies and removes false peaks from initial peak lists, reduces the burden of manual data entry, and documents and standardizes the peak filtering process. Successful preprocessing is demonstrated by the increased number of correct assignments obtained when data are submitted to an automated assignment program. AVAILABILITY: APART is available from http://sir.lanl.gov/NMR/APART.htm CONTACT: npawley@lanl.gov; rmichalczyk@lanl.gov SUPPLEMENTARY INFORMATION: Manual pages with installation instructions, procedures and screen shots can also be found at http://sir.lanl.gov/NMR/APART_Manual1.pdf.

Factors determining the reliable description of global tumbling parameters in solution NMR.

An accurate description of global tumbling of a protein is essential for correct analysis and interpretation of internal dynamics and thermodynamics. The accurate fitting of global tumbling parameters is affected by the number of experimental relaxation data points available for analysis, the distribution of data points over the domain of the function describing the tumbling, the measurement error associated with the data, the error associated with use of an approximate functional form, and errors in the protein structure. We present an analysis of the influence of these factors on the error in global tumbling parameters and the corresponding error in the calculated T(1)/T(2) values. We find that reduction of experimental and approximation error can compensate for a less-than-ideal quantity or distribution of data points, and that accurate parameters can be obtained for proteins with highly anisotropic distributions of bond vectors, as illustrated using the helical bundle protein G-CSF. This indicates that proteins with anisotropic distributions, such as the helical bundle class of proteins, should not summarily be excluded when selecting proteins for dynamic and thermodynamic analyses of (15)N backbone relaxation measurements.