Analysis of bronchoalveolar lavage fluid metatranscriptomes among patients with COVID-19 disease.

To better understand the potential relationship between COVID-19 disease and hologenome microbial community dynamics and functional profiles, we conducted a multivariate taxonomic and functional microbiome comparison of publicly available human bronchoalveolar lavage fluid (BALF) metatranscriptome samples amongst COVID-19 (n = 32), community acquired pneumonia (CAP) (n = 25), and uninfected samples (n = 29). We then performed a stratified analysis based on mortality amongst the COVID-19 cohort with known outcomes of deceased (n = 10) versus survived (n = 15). Our overarching hypothesis was that there are detectable and functionally significant relationships between BALF microbial metatranscriptomes and the severity of COVID-19 disease onset and progression. We observed 34 functionally discriminant gene ontology (GO) terms in COVID-19 disease compared to the CAP and uninfected cohorts, and 21 GO terms functionally discriminant to COVID-19 mortality (q < 0.05). GO terms enriched in the COVID-19 disease cohort included hydrolase activity, and significant GO terms under the parental terms of biological regulation, viral process, and interspecies interaction between organisms. Notable GO terms associated with COVID-19 mortality included nucleobase-containing compound biosynthetic process, organonitrogen compound catabolic process, pyrimidine-containing compound biosynthetic process, and DNA recombination, RNA binding, magnesium and zinc ion binding, oxidoreductase activity, and endopeptidase activity. A Dirichlet multinomial mixtures clustering analysis resulted in a best model fit using three distinct clusters that were significantly associated with COVID-19 disease and mortality. We additionally observed discriminant taxonomic differences associated with COVID-19 disease and mortality in the genus Sphingomonas, belonging to the Sphingomonadacae family, Variovorax, belonging to the Comamonadaceae family, and in the class Bacteroidia, belonging to the order Bacteroidales. To our knowledge, this is the first study to evaluate significant differences in taxonomic and functional signatures between BALF metatranscriptomes from COVID-19, CAP, and uninfected cohorts, as well as associating these taxa and microbial gene functions with COVID-19 mortality. Collectively, while this data does not speak to causality nor directionality of the association, it does demonstrate a significant relationship between the human microbiome and COVID-19. The results from this study have rendered testable hypotheses that warrant further investigation to better understand the causality and directionality of host-microbiome-pathogen interactions.

Making Security Viral: Shifting Engineering Biology Culture and Publishing.

The ability to construct, synthesize, and edit genes and genomes at scale and with speed enables, in synergy with other tools of engineering biology, breakthrough applications with far-reaching implications for society. As SARS-CoV-2 spread around the world in early spring of 2020, researchers rapidly mobilized, using these tools in the development of diagnostics, therapeutics, and vaccines for COVID-19. The sharing of knowledge was crucial to making rapid progress. Several publications described the use of reverse genetics for the de novo construction of SARS-CoV-2 in the laboratory, one in the form of a protocol. Given the demonstrable harm caused by the virus, the unequal distribution of mitigating vaccines and therapeutics, their unknown efficacy against variants, and the interest in this research by laboratories unaccustomed to working with highly transmissible pandemic pathogens, there are risks associated with such publications, particularly as protocols. We describe considerations and offer suggestions for enhancing security in the publication of synthetic biology research and techniques. We recommend: (1) that protocol manuscripts for the de novo synthesis of certain pathogenic viruses undergo a mandatory safety and security review; (2) that if published, such papers include descriptions of the discussions or review processes that occurred regarding security considerations in the main text; and (3) the development of a governance framework for the inclusion of basic security screening during the publication process of engineering biology/synthetic biology manuscripts to build and support a safe and secure research enterprise that is able to maximize its positive impacts and minimize any negative outcomes.

Computational analysis of pathogen-borne metallo ß-lactamases reveals discriminating structural features between B1 types.

BACKGROUND: Genes conferring antibiotic resistance to groups of bacterial pathogens are cause for considerable concern, as many once-reliable antibiotics continue to see a reduction in efficacy. The recent discovery of the metallo ß-lactamase blaNDM-1 gene, which appears to grant antibiotic resistance to a variety of Enterobacteriaceae via a mobile plasmid, is one example of this distressing trend. The following work describes a computational analysis of pathogen-borne MBLs that focuses on the structural aspects of characterized proteins. RESULTS: Using both sequence and structural analyses, we examine residues and structural features specific to various pathogen-borne MBL types. This analysis identifies a linker region within MBL-like folds that may act as a discriminating structural feature between these proteins, and specifically resistance-associated acquirable MBLs. Recently released crystal structures of the newly emerged NDM-1 protein were aligned against related MBL structures using a variety of global and local structural alignment methods, and the overall fold conformation is examined for structural conservation. Conservation appears to be present in most areas of the protein, yet is strikingly absent within a linker region, making NDM-1 unique with respect to a linker-based classification scheme. Variability analysis of the NDM-1 crystal structure highlights unique residues in key regions as well as identifying several characteristics shared with other transferable MBLs. CONCLUSIONS: A discriminating linker region identified in MBL proteins is highlighted and examined in the context of NDM-1 and primarily three other MBL types: IMP-1, VIM-2 and ccrA. The presence of an unusual linker region variant and uncommon amino acid composition at specific structurally important sites may help to explain the unusually broad kinetic profile of NDM-1 and may aid in directing research attention to areas of this protein, and possibly other MBLs, that may be targeted for inactivation or attenuation of enzymatic activity.

LAVA: an open-source approach to designing LAMP (loop-mediated isothermal amplification) DNA signatures.

BACKGROUND: We developed an extendable open-source Loop-mediated isothermal AMPlification (LAMP) signature design program called LAVA (LAMP Assay Versatile Analysis). LAVA was created in response to limitations of existing LAMP signature programs. RESULTS: LAVA identifies combinations of six primer regions for basic LAMP signatures, or combinations of eight primer regions for LAMP signatures with loop primers, which can be used as LAMP signatures. The identified primers are conserved among target organism sequences. Primer combinations are optimized based on lengths, melting temperatures, and spacing among primer sites. We compare LAMP signature candidates for Staphylococcus aureus created both by LAVA and by PrimerExplorer. We also include signatures from a sample run targeting all strains of Mycobacterium tuberculosis. CONCLUSIONS: We have designed and demonstrated new software for identifying signature candidates appropriate for LAMP assays. The software is available for download at http://lava-dna.googlecode.com/.

A sample-in-answer-out instrument for the detection of multiple respiratory pathogens in unprepared nasopharyngeal swab samples.

Multiplex RT-PCR suspension array assays provide a powerful tool for identifying the causative agent(s) of respiratory infections. These assays are time consuming and laborious on a time-per-sample basis if only a few samples require processing. To address this shortcoming and provide an automated solution for fast detection and identification of viral pathogens, we developed the first automated multiplex RT-PCR suspension array instrument capable of handling unprepared clinical samples. The instrument requires less than 3 minutes of hands-on time for a result generated in approximately 2.5 hours. In analytical studies, the instrument performed as well as manually performed assays. The performance of the instrument and loaded multiplex viral detection assay was then tested using unprepared nasopharyngeal samples. The instrument-performed assay detected 61 of 71 RSV positive samples, for a sensitivity of 85.9%. Adenovirus (n = 5) and influenza B (n = 3) were less prevalent in the sample set, but detected to similar levels, 80% and 75%, respectively. The same sample set was also tested using FDA approved immuno-assay rapid tests, and the instrument was found to be more sensitive than the rapid tests with the sole exception being influenza A (n = 16), which was poorly detected due to significant sequence mismatches between the influenza A primer/probe set included in the multiplex mixture and the circulating influenza A strains. Overall, these data demonstrate the developed prototype platform performs multiplex array assays as well as hand-performed assays, and that the instrument's sensitivity and specificity are dictated by the quality of the loaded multiplex assay.

Conserved amino acid markers from past influenza pandemic strains.

BACKGROUND: Finding the amino acid mutations that affect the severity of influenza infections remains an open and challenging problem. Of special interest is better understanding how current circulating influenza strains could evolve into a new pandemic strain. Influenza proteomes from distinct viral phenotype classes were searched for class specific amino acid mutations conserved in past pandemics, using reverse engineered linear classifiers. RESULTS: Thirty-four amino acid markers associated with host specificity and high mortality rate were found. Some markers had little impact on distinguishing the functional classes by themselves, however in combination with other mutations they improved class prediction. Pairwise combinations of influenza genomes were checked for reassortment and mutation events needed to acquire the pandemic conserved markers. Evolutionary pathways involving H1N1 human and swine strains mixed with avian strains show the potential to acquire the pandemic markers with a double reassortment and one or two amino acid mutations. CONCLUSION: The small mutation combinations found at multiple protein positions associated with viral phenotype indicate that surveillance tools could monitor genetic variation beyond single point mutations to track influenza strains. Finding that certain strain combinations have the potential to acquire pandemic conserved markers through a limited number of reassortment and mutation events illustrates the potential for reassortment and mutation events to lead to new circulating influenza strains.

Multiplexed molecular assay for rapid exclusion of foot-and-mouth disease.

A nucleic acid-based multiplexed assay was developed that combines detection of foot-and-mouth disease virus (FMDV) with rule-out assays for two other foreign animal diseases and four domestic animal diseases that cause vesicular or ulcerative lesions indistinguishable from FMDV infection in cattle, sheep and swine. The FMDV "look-alike" diagnostic assay panel contains 5 PCR and 12 reverse transcriptase PCR (RT-PCR) signatures for a total of 17 simultaneous PCR amplifications for 7 diseases plus incorporating 4 internal assay controls. It was developed and optimized to amplify both DNA and RNA viruses simultaneously in a single tube and employs Luminex liquid array technology. Assay development including selection of appropriate controls, a comparison of signature performance in single and multiplex testing against target nucleic acids, as well of limits of detection for each of the individual signatures is presented. While this assay is a prototype and by no means a comprehensive test for FMDV "look-alike" viruses, an assay of this type is envisioned to have benefit to a laboratory network in routine surveillance and possibly for post-outbreak proof of freedom from foot-and-mouth disease.

Diagnostic evaluation of multiplexed reverse transcription-PCR microsphere array assay for detection of foot-and-mouth and look-alike disease viruses.

A high-throughput multiplexed assay was developed for the differential laboratory detection of foot-and-mouth disease virus (FMDV) from viruses that cause clinically similar diseases of livestock. This assay simultaneously screens for five RNA and two DNA viruses by using multiplexed reverse transcription-PCR (mRT-PCR) amplification coupled with a microsphere hybridization array and flow-cytometric detection. Two of the 17 primer-probe sets included in this multiplex assay were adopted from previously characterized real-time RT-PCR (rRT-PCR) assays for FMDV. The diagnostic accuracy of the mRT-PCR assay was evaluated using 287 field samples, including 247 samples (213 true-positive samples and 35 true-negative samples) from suspected cases of foot-and-mouth disease collected from 65 countries between 1965 and 2006 and 39 true-negative samples collected from healthy animals. The mRT-PCR assay results were compared to those of two singleplex rRT-PCR assays, using virus isolation with antigen enzyme-linked immunosorbent assays as the reference method. The diagnostic sensitivity of the mRT-PCR assay for FMDV was 93.9% (95% confidence interval [CI], 89.8 to 96.4%), and the sensitivity was 98.1% (95% CI, 95.3 to 99.3%) for the two singleplex rRT-PCR assays used in combination. In addition, the assay could reliably differentiate between FMDV and other vesicular viruses, such as swine vesicular disease virus and vesicular exanthema of swine virus. Interestingly, the mRT-PCR detected parapoxvirus (n = 2) and bovine viral diarrhea virus (n = 2) in clinical samples, demonstrating the screening potential of this mRT-PCR assay to identify viruses in FMDV-negative material not previously recognized by using focused single-target rRT-PCR assays.

On the origin of smallpox: correlating variola phylogenics with historical smallpox records.

Human disease likely attributable to variola virus (VARV), the etiologic agent of smallpox, has been reported in human populations for >2,000 years. VARV is unique among orthopoxviruses in that it is an exclusively human pathogen. Because VARV has a large, slowly evolving DNA genome, we were able to construct a robust phylogeny of VARV by analyzing concatenated single nucleotide polymorphisms (SNPs) from genome sequences of 47 VARV isolates with broad geographic distributions. Our results show two primary VARV clades, which likely diverged from an ancestral African rodent-borne variola-like virus either approximately 16,000 or approximately 68,000 years before present (YBP), depending on which historical records (East Asian or African) are used to calibrate the molecular clock. One primary clade was represented by the Asian VARV major strains, the more clinically severe form of smallpox, which spread from Asia either 400 or 1,600 YBP. Another primary clade included both alastrim minor, a phenotypically mild smallpox described from the American continents, and isolates from West Africa. This clade diverged from an ancestral VARV either 1,400 or 6,300 YBP, and then further diverged into two subclades at least 800 YBP. All of these analyses indicate that the divergence of alastrim and variola major occurred earlier than previously believed.

MannDB - a microbial database of automated protein sequence analyses and evidence integration for protein characterization.

BACKGROUND: MannDB was created to meet a need for rapid, comprehensive automated protein sequence analyses to support selection of proteins suitable as targets for driving the development of reagents for pathogen or protein toxin detection. Because a large number of open-source tools were needed, it was necessary to produce a software system to scale the computations for whole-proteome analysis. Thus, we built a fully automated system for executing software tools and for storage, integration, and display of automated protein sequence analysis and annotation data. DESCRIPTION: MannDB is a relational database that organizes data resulting from fully automated, high-throughput protein-sequence analyses using open-source tools. Types of analyses provided include predictions of cleavage, chemical properties, classification, features, functional assignment, post-translational modifications, motifs, antigenicity, and secondary structure. Proteomes (lists of hypothetical and known proteins) are downloaded and parsed from Genbank and then inserted into MannDB, and annotations from SwissProt are downloaded when identifiers are found in the Genbank entry or when identical sequences are identified. Currently 36 open-source tools are run against MannDB protein sequences either on local systems or by means of batch submission to external servers. In addition, BLAST against protein entries in MvirDB, our database of microbial virulence factors, is performed. A web client browser enables viewing of computational results and downloaded annotations, and a query tool enables structured and free-text search capabilities. When available, links to external databases, including MvirDB, are provided. MannDB contains whole-proteome analyses for at least one representative organism from each category of biological threat organism listed by APHIS, CDC, HHS, NIAID, USDA, USFDA, and WHO. CONCLUSION: MannDB comprises a large number of genomes and comprehensive protein sequence analyses representing organisms listed as high-priority agents on the websites of several governmental organizations concerned with bio-terrorism. MannDB provides the user with a BLAST interface for comparison of native and non-native sequences and a query tool for conveniently selecting proteins of interest. In addition, the user has access to a web-based browser that compiles comprehensive and extensive reports. Access to MannDB is freely available at http://manndb.llnl.gov/.

Cynomolgus macaque as an animal model for severe acute respiratory syndrome.

BACKGROUND: The emergence of severe acute respiratory syndrome (SARS) in 2002 and 2003 affected global health and caused major economic disruption. Adequate animal models are required to study the underlying pathogenesis of SARS-associated coronavirus (SARS-CoV) infection and to develop effective vaccines and therapeutics. We report the first findings of measurable clinical disease in nonhuman primates (NHPs) infected with SARS-CoV. METHODS AND FINDINGS: In order to characterize clinically relevant parameters of SARS-CoV infection in NHPs, we infected cynomolgus macaques with SARS-CoV in three groups: Group I was infected in the nares and bronchus, group II in the nares and conjunctiva, and group III intravenously. Nonhuman primates in groups I and II developed mild to moderate symptomatic illness. All NHPs demonstrated evidence of viral replication and developed neutralizing antibodies. Chest radiographs from several animals in groups I and II revealed unifocal or multifocal pneumonia that peaked between days 8 and 10 postinfection. Clinical laboratory tests were not significantly changed. Overall, inoculation by a mucosal route produced more prominent disease than did intravenous inoculation. Half of the group I animals were infected with a recombinant infectious clone SARS-CoV derived from the SARS-CoV Urbani strain. This infectious clone produced disease indistinguishable from wild-type Urbani strain. CONCLUSIONS: SARS-CoV infection of cynomolgus macaques did not reproduce the severe illness seen in the majority of adult human cases of SARS; however, our results suggest similarities to the milder syndrome of SARS-CoV infection characteristically seen in young children.

Comparative genomics tools applied to bioterrorism defence.

Rapid advances in the genomic sequencing of bacteria and viruses over the past few years have made it possible to consider sequencing the genomes of all pathogens that affect humans and the crops and livestock upon which our lives depend. Recent events make it imperative that full genome sequencing be accomplished as soon as possible for pathogens that could be used as weapons of mass destruction or disruption. This sequence information must be exploited to provide rapid and accurate diagnostics to identify pathogens and distinguish them from harmless near-neighbours and hoaxes. The Chem-Bio Non-Proliferation (CBNP) programme of the US Department of Energy (DOE) began a large-scale effort of pathogen detection in early 2000 when it was announced that the DOE would be providing bio-security at the 2002 Winter Olympic Games in Salt Lake City, Utah. Our team at the Lawrence Livermore National Lab (LLNL) was given the task of developing reliable and validated assays for a number of the most likely bioterrorist agents. The short timeline led us to devise a novel system that utilised whole-genome comparison methods to rapidly focus on parts of the pathogen genomes that had a high probability of being unique. Assays developed with this approach have been validated by the Centers for Disease Control (CDC). They were used at the 2002 Winter Olympics, have entered the public health system, and have been in continual use for non-publicised aspects of homeland defence since autumn 2001. Assays have been developed for all major threat list agents for which adequate genomic sequence is available, as well as for other pathogens requested by various government agencies. Collaborations with comparative genomics algorithm developers have enabled our LLNL team to make major advances in pathogen detection, since many of the existing tools simply did not scale well enough to be of practical use for this application. It is hoped that a discussion of a real-life practical application of comparative genomics algorithms may help spur algorithm developers to tackle some of the many remaining problems that need to be addressed. Solutions to these problems will advance a wide range of biological disciplines, only one of which is pathogen detection. For example, exploration in evolution and phylogenetics, annotating gene coding regions, predicting and understanding gene function and regulation, and untangling gene networks all rely on tools for aligning multiple sequences, detecting gene rearrangements and duplications, and visualising genomic data. Two key problems currently needing improved solutions are: (1) aligning incomplete, fragmentary sequence (eg draft genome contigs or arbitrary genome regions) with both complete genomes and other fragmentary sequences; and (2) ordering, aligning and visualising non-colinear gene rearrangements and inversions in addition to the colinear alignments handled by current tools.

Limitations of TaqMan PCR for detecting divergent viral pathogens illustrated by hepatitis A, B, C, and E viruses and human immunodeficiency virus.

Recent events illustrate the imperative to rapidly and accurately detect and identify pathogens during disease outbreaks, whether they are natural or engineered. Particularly for our primary goal of detecting bioterrorist releases, detection techniques must be both species-wide (capable of detecting all known strains of a given species) and species specific. Due to classification restrictions on the publication of data for species that may pose a bioterror threat, we illustrate the challenges of finding such assays using five nonthreat organisms that are nevertheless of public health concern: human immunodeficiency virus (HIV) and four species of hepatitis viruses. Fluorogenic probe-based PCR assays (TaqMan; Perkin-Elmer Corp., Applied Biosystems, Foster City, Calif.) may be sensitive, fast methods for the identification of species in which the genome is conserved among strains, such as hepatitis A virus. For species such as HIV, however, the strains are highly divergent. We use computational methods to show that nine TaqMan primer and probe sequences, or signatures, are needed to ensure that all strains will be detected, but this is an unfeasible number, considering the cost of TaqMan probes. Strains of hepatitis B, C, and E viruses show intermediate divergence, so that two to three TaqMan signatures are required to detect all strains of each virus. We conclude that for species such as hepatitis A virus with high levels of sequence conservation among strains, signatures can be found computationally for detection by the TaqMan assay, which is a sensitive, rapid, and cost-effective method. However, for species such as HIV with substantial genetic divergence among strains, the TaqMan assay becomes unfeasible and alternative detection methods may be required. We compare the TaqMan assay with some of the alternative nucleic acid-based detection techniques of microarray, chip, and bead technologies in terms of sensitivity, speed, and cost.