UltraSEQ, a Universal Bioinformatic Platform for Information-Based Clinical Metagenomics and Beyond.

Applied metagenomics is a powerful emerging capability enabling the untargeted detection of pathogens, and its application in clinical diagnostics promises to alleviate the limitations of current targeted assays. While metagenomics offers a hypothesis-free approach to identify any pathogen, including unculturable and potentially novel pathogens, its application in clinical diagnostics has so far been limited by workflow-specific requirements, computational constraints, and lengthy expert review requirements. To address these challenges, we developed UltraSEQ, a first-of-its-kind accurate and scalable metagenomic bioinformatic tool for potential clinical diagnostics and biosurveillance utility. Here, we present the results of the evaluation of our novel UltraSEQ pipeline using an in silico-synthesized metagenome, mock microbial community data sets, and publicly available clinical data sets from samples of different infection types, including both short-read and long-read sequencing data. Our results show that UltraSEQ successfully detected all expected species across the tree of life in the in silico sample and detected all 10 bacterial and fungal species in the mock microbial community data set. For clinical data sets, even without requiring data set-specific configuration setting changes, background sample subtraction, or prior sample information, UltraSEQ achieved an overall accuracy of 91%. Furthermore, as an initial demonstration with a limited patient sample set, we show UltraSEQ's ability to provide antibiotic resistance and virulence factor genotypes that are consistent with phenotypic results. Taken together, the above-described results demonstrate that the UltraSEQ platform offers a transformative approach for microbial and metagenomic sample characterization, employing a biologically informed detection logic, deep metadata, and a flexible system architecture for the classification and characterization of taxonomic origin, gene function, and user-defined functions, including disease-causing infections. IMPORTANCE Traditional clinical microbiology-based diagnostic tests rely on targeted methods that can detect only one to a few preselected organisms or slow, culture-based methods. Although widely used today, these methods have several limitations, resulting in rates of cases of an unknown etiology of infection of >50% for several disease types. Massive developments in sequencing technologies have made it possible to apply metagenomic methods to clinical diagnostics, but current offerings are limited to a specific disease type or sequencer workflow and/or require laboratory-specific controls. The limitations associated with current clinical metagenomic offerings result from the fact that the backend bioinformatic pipelines are optimized for the specific parameters described above, resulting in an excess of unmaintained, redundant, and niche tools that lack standardization and explainable outputs. In this paper, we demonstrate that UltraSEQ uses a novel, information-based approach that enables accurate, evidence-based predictions for diagnosis as well as the functional characterization of a sample.

Function-based classification of hazardous biological sequences: Demonstration of a new paradigm for biohazard assessments.

Bioengineering applies analytical and engineering principles to identify functional biological building blocks for biotechnology applications. While these building blocks are leveraged to improve the human condition, the lack of simplistic, machine-readable definition of biohazards at the function level is creating a gap for biosafety practices. More specifically, traditional safety practices focus on the biohazards of known pathogens at the organism-level and may not accurately consider novel biodesigns with engineered functionalities at the genetic component-level. This gap is motivating the need for a paradigm shift from organism-centric procedures to function-centric biohazard identification and classification practices. To address this challenge, we present a novel methodology for classifying biohazards at the individual sequence level, which we then compiled to distinguish the biohazardous property of pathogenicity at the whole genome level. Our methodology is rooted in compilation of hazardous functions, defined as a set of sequences and associated metadata that describe coarse-level functions associated with pathogens (e.g., adherence, immune subversion). We demonstrate that the resulting database can be used to develop hazardous "fingerprints" based on the functional metadata categories. We verified that these hazardous functions are found at higher levels in pathogens compared to non-pathogens, and hierarchical clustering of the fingerprints can distinguish between these two groups. The methodology presented here defines the hazardous functions associated with bioengineering functional building blocks at the sequence level, which provide a foundational framework for classifying biological hazards at the organism level, thus leading to the improvement and standardization of current biosecurity and biosafety practices.

Metagenome-Guided Proteomic Quantification of Reductive Dehalogenases in the Dehalococcoides mccartyi-Containing Consortium SDC-9.

At groundwater sites contaminated with chlorinated ethenes, fermentable substrates are often added to promote reductive dehalogenation by indigenous or augmented microorganisms. Contemporary bioremediation performance monitoring relies on nucleic acid biomarkers of key organohalide-respiring bacteria, such as Dehalococcoides mccartyi (Dhc). Metagenome sequencing of the commercial, Dhc-containing consortium, SDC-9, identified 12 reductive dehalogenase (RDase) genes, including pceA (two copies), vcrA, and tceA, and allowed for specific detection and quantification of RDase peptides using liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS). Shotgun (i.e., untargeted) proteomics applied to the SDC-9 consortium grown with tetrachloroethene (PCE) and lactate identified 143 RDase peptides, and 36 distinct peptides that covered greater than 99% of the protein-coding sequences of the PceA, TceA, and VcrA RDases. Quantification of RDase peptides using multiple reaction monitoring (MRM) assays with 13C-/15N-labeled peptides determined 1.8 × 103 TceA and 1.2 × 102 VcrA RDase molecules per Dhc cell. The MRM mass spectrometry approach allowed for sensitive detection and accurate quantification of relevant Dhc RDases and has potential utility in bioremediation monitoring regimes.

Use of omic tools to assess methyl tert-butyl ether (MTBE) degradation in groundwater.

This study employed innovative technologies to evaluate multiple lines of evidence for natural attenuation (NA) of methyl tertiary-butyl ether (MTBE) in groundwater at the 22 Area of Marine Corps Base (MCB) Camp Pendleton after decommissioning of a biobarrier system. For comparison, data from the 13 Area Gas Station where active treatment of MTBE is occurring was used to evaluate the effectiveness of omic techniques in assessing biodegradation. Overall, the 22 Area Gas Station appeared to be anoxic. MTBE was detected in large portion of the plume. In comparison, concentrations of MTBE at the 13 Area Gas Station were much higher (42,000 µg/L to 2800 µg/L); however, none of the oxygenates were detected. Metagenomic analysis of the indigenous groundwater microbial community revealed the presence of bacterial strains known to aerobically and anaerobically degrade MTBE at both sites. While proteomic analysis at the 22 Area Gas Station showed the presence of proteins of MTBE degrading microorganisms, the MTBE degradative proteins were only found at the 13 Area Gas Station. Taken together, these results provide evidence for previous NA of MTBE in the groundwater at 22 Area Gas Station and demonstrate the effectiveness of innovative-omic technologies to assist monitored NA assessments.

Stability Characterization of a Polysorbate 80-Dimethyl Trisulfide Formulation, a Cyanide Antidote Candidate.

Novel cyanide countermeasures are needed for cases of a mass-exposure cyanide emergency. A lead candidate compound is dimethyl trisulfide (DMTS), which acts as a sulfur donor for rhodanese, thereby assisting the conversion of cyanide into thiocyanate. DMTS is a safe compound for consumption and, in a 15% polysorbate 80 (DMTS-PS80) formulation, has demonstrated good efficacy against cyanide poisoning in several animal models. We performed a stability study that investigated the effect of temperature, location of formulation preparation, and pH under buffered conditions. We found that while the stability of the DMTS component was fairly independent of which laboratory prepared the formulation, the concentration of DMTS in the formulation was reduced 36-58% over the course of 29 weeks when stored at room temperature. This loss typically increased with increasing temperatures, although we did not find statistical differences between the stability at different storage temperatures in all formulations. Further, we found that addition of a light buffer negatively impacted the stability, whereas the pH of that buffer did not impact stability. We investigated the factors behind the reduction of DMTS over time using various techniques, and we suggest that the instability of the formulation is governed at least partially by precipitation and evaporation, although a combination of factors is likely involved.

Next-generation sequencing approach to epigenetic-based tissue source attribution.

The ability to determine the tissue source of biological materials from evidence samples can be highly informative for interpreting forensic data. In this study, a previously published CE-based method to probe locus-specific DNA methylation was modified to accommodate detection using next-generation sequencing (NGS) to perform tissue source attribution. DNA samples (1 ng) from each of four different tissue types were digested with the methylation sensitive restriction endonuclease Hha1 and PCR was used to amplify an optimized subset of ten methylated loci, including positive and negative control loci. The products were prepared as NGS libraries, pooled in a multiplex assay with sample-specific barcodes, sequenced with an Illumina MiSeq, and analyzed using a k-Nearest Neighbor algorithm. With this initial effort a concordance rate of 15/16 was demonstrated from samples of varying types: semen, saliva, skin epidermis, and blood. This method also was designed to be compatible with the workflows published to date for NGS of STRs. Thus, the methylation approach described here is highly accurate and upon further validation and testing may be potentially used in practice as a confirmatory test in conjunction with other NGS protocols used in forensic laboratories.

Non-enzymatically derived minor lipids found in Escherichia coli lipid extracts.

Electrospray ionization mass spectrometry is a powerful technique to analyze lipid extracts especially for the identification of new lipid metabolites. A hurdle to lipid identification is the presence of solvent contaminants that hinder the identification of low abundance species or covalently modify abundant lipid species. We have identified several non-enzymatically derived minor lipid species in lipid extracts of Escherichia coli; phosphatidylmethanol, ethyl and methyl carbamates of PE and N-succinyl PE were identified in lipid extracts of E. coli. Phosphatidylmethanol (PM) was identified by exact mass measurement and collision induced dissociation tandem mass spectrometry (MS/MS). Extraction in the presence of deuterated methanol leads to a 3 atomic mass unit shift in the [M-H](-) ions of PM indicating its formation during extraction. Ethyl and methyl carbamates of PE, also identified by exact mass measurement and MS/MS, are likely to be formed by phosgene, a breakdown product of chloroform. Addition of phosgene to extractions containing synthetic PE significantly increases the levels of PE-MC detected in the lipid extracts by ESI-MS. Extraction in the presence of methylene chloride significantly reduced the levels of these lipid species. N-succinyl PE is formed from reaction of succinyl-CoA with PE during extraction. Interestingly N-succinyl PE can be formed in an aqueous reaction mixture in the absence of added E. coli proteins. This work highlights the reactivity of the amine of PE and emphasizes that careful extraction controls are required to ensure that new minor lipid species identified using mass spectrometry are indeed endogenous lipid metabolites.

Structural differences in lipomannans from pathogenic and nonpathogenic mycobacteria that impact CD1b-restricted T cell responses.

Mannosylated molecules on the Mycobacterium tuberculosis surface are important determinants in the immunopathogenesis of tuberculosis. To date, much attention has been paid to mannose-capped lipoarabinomannan, which mediates phagocytosis and intracellular trafficking of M. tuberculosis by engaging the macrophage mannose receptor and subsequently binds to intracellular CD1b molecules for presentation to T cells. Another important mannosylated lipoglycan on the M. tuberculosis surface is lipomannan (LM). Comparative structural detail of the LMs from virulent and avirulent strains is limited as is knowledge regarding their differential capacity to be recognized by the adaptive immune response. Here, we purified LM from the avirulent M. smegmatis and the virulent M. tuberculosis H(37)R(v), performed a comparative structural biochemical analysis, and addressed their ability to stimulate CD1b-restricted T cell clones. We found that M. tuberculosis H(37)R(v) produces a large neutral LM (TB-LM); in contrast, M. smegmatis produces a smaller linear acidic LM (SmegLM) with a high succinate content. Correspondingly, TB-LM was not as efficiently presented to CD1b-restricted T cells as SmegLM. Thus, here we correlate the structure-function relationships for LMs with CD1b-restricted T cell responses and provide evidence that the structural features of TB-LM contribute to its diminished T cell responsiveness.

Crystal structure and acyl chain selectivity of Escherichia coli LpxD, the N-acyltransferase of lipid A biosynthesis.

LpxD catalyzes the third step of lipid A biosynthesis, the R-3-hydroxyacyl-ACP-dependent N-acylation of UDP-3-O-(acyl)-alpha-D-glucosamine, and is a target for new antibiotic development. Here we report the 2.6 A crystal structure of the Escherichia coli LpxD homotrimer (EcLpxD). As is the case in Chlamydia trachomatis LpxD (CtLxpD), each EcLpxD chain consists of an N-terminal uridine-binding region, a left-handed parallel beta-helix (LbetaH), and a C-terminal alpha-helical domain. The backbones of the LbetaH domains of the two enzymes are similar, as are the positions of key active site residues. The N-terminal nucleotide binding domains are oriented differently relative to the LbetaH regions, but are similar when overlaid on each other. The orientation of the EcLpxD tripeptide (residues 303-305), connecting the distal end of the LbetaH and the proximal end of the C-terminal helical domains, differs from its counterpart in CtLpxD (residues 311-312); this results in a 120 degrees rotation of the C-terminal domain relative to the LbetaH region in EcLpxD versus CtLpxD. M290 of EcLpxD appears to cap the distal end of a hydrophobic cleft that binds the acyl chain of the R-3-hydroxyacyl-ACP donor substrate. Under standard assay conditions, wild-type EcLpxD prefers R,S-3-hydroxymyristoyl-ACP over R,S-3-hydroxypalmitoyl-ACP by a factor of 3, whereas the M290A mutant has the opposite selectivity. Both wild-type and M290A EcLpxD rescue the conditional lethality of E. coli RL25, a temperature-sensitive strain harboring point mutations in lpxD. Complementation with wild-type EcLpxD restores normal lipid A containing only N-linked hydroxymyristate to RL25 at 42 degrees C, as judged by mass spectrometry, whereas the M290A mutant generates multiple lipid A species containing one or two longer hydroxy fatty acids in place of the usual R-3-hydroxymyristate at positions 2 and 2'.

Steady-state kinetics and mechanism of LpxD, the N-acyltransferase of lipid A biosynthesis.

LpxD catalyzes the third step of lipid A biosynthesis, the (R)-3-hydroxymyristoyl-acyl carrier protein ( R-3-OHC14-ACP)-dependent N-acylation of UDP-3-O-[(R)-3-hydroxymyristoyl]-alpha-D-glucosamine [UDP-3-O-(R-3-OHC14)-GlcN]. We have now overexpressed and purified Escherichia coli LpxD to homogeneity. Steady-state kinetics suggest a compulsory ordered mechanism in which R-3-OHC14-ACP binds prior to UDP-3-O-(R-3-OHC14)-GlcN. The product, UDP-2,3-diacylglucosamine, dissociates prior to ACP; the latter is a competitive inhibitor against R-3-OHC14-ACP and a noncompetitive inhibitor against UDP-3-O-(R-3-OHC14)-GlcN. UDP-2-N-[(R)-3-Hydroxymyristoyl]-alpha-D-glucosamine, obtained by mild base hydrolysis of UDP-2,3-diacylglucosamine, is a noncompetitive inhibitor against both substrates. Synthetic (R)-3-hydroxylauroyl-methylphosphopantetheine is an uncompetitive inhibitor against R-3-OHC14-ACP and a competitive inhibitor against UDP-3-O-(R-3-OHC14)-GlcN, but (R)-3-hydroxylauroyl-methylphosphopantetheine is also a very poor substrate. A compulsory ordered mechanism is consistent with the fact that R-3-OHC14-ACP has a high binding affinity for free LpxD whereas UDP-3-O-(R-3-OHC14)-GlcN does not. Divalent cations inhibit R-3-OHC14-ACP-dependent acylation but not (R)-3-hydroxylauroyl-methylphosphopantetheine-dependent acylation, indicating that the acidic recognition helix of R-3-OHC14-ACP contributes to binding. The F41A mutation increases the K(M) for UDP-3-O-(R-3-OHC14)-GlcN 30-fold, consistent with aromatic stacking of the corresponding F43 side chain against the uracil moiety of bound UDP-GlcNAc in the X-ray structure of Chlamydia trachomatis LpxD. Mutagenesis implicates E. coli H239 but excludes H276 as the catalytic base, and neither residue is likely to stabilize the oxyanion intermediate.