Whole genome sequencing of human Borrelia burgdorferi isolates reveals linked blocks of accessory genome elements located on plasmids and associated with human dissemination.

Lyme disease is the most common vector-borne disease in North America and Europe. The clinical manifestations of Lyme disease vary based on the genospecies of the infecting Borrelia burgdorferi spirochete, but the microbial genetic elements underlying these associations are not known. Here, we report the whole genome sequence (WGS) and analysis of 299 B. burgdorferi (Bb) isolates derived from patients in the Eastern and Midwestern US and Central Europe. We develop a WGS-based classification of Bb isolates, confirm and extend the findings of previous single- and multi-locus typing systems, define the plasmid profiles of human-infectious Bb isolates, annotate the core and strain-variable surface lipoproteome, and identify loci associated with disseminated infection. A core genome consisting of ~900 open reading frames and a core set of plasmids consisting of lp17, lp25, lp36, lp28-3, lp28-4, lp54, and cp26 are found in nearly all isolates. Strain-variable (accessory) plasmids and genes correlate strongly with phylogeny. Using genetic association study methods, we identify an accessory genome signature associated with dissemination in humans and define the individual plasmids and genes that make up this signature. Strains within the RST1/WGS A subgroup, particularly a subset marked by the OspC type A genotype, have increased rates of dissemination in humans. OspC type A strains possess a unique set of strongly linked genetic elements including the presence of lp56 and lp28-1 plasmids and a cluster of genes that may contribute to their enhanced virulence compared to other genotypes. These features of OspC type A strains reflect a broader paradigm across Bb isolates, in which near-clonal genotypes are defined by strain-specific clusters of linked genetic elements, particularly those encoding surface-exposed lipoproteins. These clusters of genes are maintained by strain-specific patterns of plasmid occupancy and are associated with the probability of invasive infection.

Whole genome sequencing of Borrelia burgdorferi isolates reveals linked clusters of plasmid-borne accessory genome elements associated with virulence.

Lyme disease is the most common vector-borne disease in North America and Europe. The clinical manifestations of Lyme disease vary based on the genospecies of the infecting Borrelia burgdorferi spirochete, but the microbial genetic elements underlying these associations are not known. Here, we report the whole genome sequence (WGS) and analysis of 299 patient-derived B. burgdorferi sensu stricto ( Bbss ) isolates from patients in the Eastern and Midwestern US and Central Europe. We develop a WGS-based classification of Bbss isolates, confirm and extend the findings of previous single- and multi-locus typing systems, define the plasmid profiles of human-infectious Bbss isolates, annotate the core and strain-variable surface lipoproteome, and identify loci associated with disseminated infection. A core genome consisting of ∻800 open reading frames and a core set of plasmids consisting of lp17, lp25, lp36, lp28-3, lp28-4, lp54, and cp26 are found in nearly all isolates. Strain-variable (accessory) plasmids and genes correlate strongly with phylogeny. Using genetic association study methods, we identify an accessory genome signature associated with dissemination and define the individual plasmids and genes that make up this signature. Strains within the RST1/WGS A subgroup, particularly a subset marked by the OspC type A genotype, are associated with increased rates of dissemination. OspC type A strains possess a unique constellation of strongly linked genetic changes including the presence of lp56 and lp28-1 plasmids and a cluster of genes that may contribute to their enhanced virulence compared to other genotypes. The patterns of OspC type A strains typify a broader paradigm across Bbss isolates, in which genetic structure is defined by correlated groups of strain-variable genes located predominantly on plasmids, particularly for expression of surface-exposed lipoproteins. These clusters of genes are inherited in blocks through strain-specific patterns of plasmid occupancy and are associated with the probability of invasive infection.

Enhanced Virus Detection and Metagenomic Sequencing in Patients with Meningitis and Encephalitis.

Meningitis and encephalitis are leading causes of central nervous system (CNS) disease and often result in severe neurological compromise or death. Traditional diagnostic workflows largely rely on pathogen-specific tests, sometimes over days to weeks, whereas metagenomic next-generation sequencing (mNGS) profiles all nucleic acid in a sample. In this single-center, prospective study, 68 hospitalized patients with known (n = 44) or suspected (n = 24) CNS infections underwent mNGS from RNA and DNA to identify potential pathogens and also targeted sequencing of viruses using hybrid capture. Using a computational metagenomic classification pipeline based on KrakenUniq and BLAST, we detected pathogen nucleic acid in cerebrospinal fluid (CSF) from 22 subjects, 3 of whom had no clinical diagnosis by routine workup. Among subjects diagnosed with infection by serology and/or peripheral samples, we demonstrated the utility of mNGS to detect pathogen nucleic acid in CSF, importantly for the Ixodes scapularis tick-borne pathogens Powassan virus, Borrelia burgdorferi, and Anaplasma phagocytophilum. We also evaluated two methods to enhance the detection of viral nucleic acid, hybrid capture and methylated DNA depletion. Hybrid capture nearly universally increased viral read recovery. Although results for methylated DNA depletion were mixed, it allowed the detection of varicella-zoster virus DNA in two samples that were negative by standard mNGS. Overall, mNGS is a promising approach that can test for multiple pathogens simultaneously, with efficacy similar to that of pathogen-specific tests, and can uncover geographically relevant infectious CNS disease, such as tick-borne infections in New England. With further laboratory and computational enhancements, mNGS may become a mainstay of workup for encephalitis and meningitis. IMPORTANCE Meningitis and encephalitis are leading global causes of central nervous system (CNS) disability and mortality. Current diagnostic workflows remain inefficient, requiring costly pathogen-specific assays and sometimes invasive surgical procedures. Despite intensive diagnostic efforts, 40 to 60% of people with meningitis or encephalitis have no clear cause of CNS disease identified. As diagnostic uncertainty often leads to costly inappropriate therapies, the need for novel pathogen detection methods is paramount. Metagenomic next-generation sequencing (mNGS) offers the unique opportunity to circumvent these challenges using unbiased laboratory and computational methods. Here, we performed comprehensive mNGS from 68 prospectively enrolled patients with known (n = 44) or suspected (n = 24) CNS viral infection from a single center in New England and evaluated enhanced methods to improve the detection of CNS pathogens, including those not traditionally identified in the CNS by nucleic acid detection. Overall, our work helps elucidate how mNGS can become integrated into the diagnostic toolkit for CNS infections.

Negative Selection by Spiral Inertial Microfluidics Improves Viral Recovery and Sequencing from Blood.

In blood samples from patients with viral infection, it is often important to separate viral particles from human cells, for example, to minimize background in performing viral whole genome sequencing. Here, we present a microfluidic device that uses spiral inertial microfluidics with continuous circulation to separate host cells from viral particles and free nucleic acid. We demonstrate that this device effectively reduces white blood cells, red blood cells, and platelets from both whole blood and plasma samples with excellent recovery of viral nucleic acid. Furthermore, microfluidic separation leads to greater viral genome coverage and depth, highlighting an important application of this device in processing clinical samples for viral genome sequencing.

Rapid Detection of Powassan Virus in a Patient With Encephalitis by Metagenomic Sequencing.

We describe a patient with severe and progressive encephalitis of unknown etiology. We performed rapid metagenomic sequencing from cerebrospinal fluid and identified Powassan virus, an emerging tick-borne flavivirus that has been increasingly detected in the United States.

A global map of genetic diversity in Babesia microti reveals strong population structure and identifies variants associated with clinical relapse.

Human babesiosis caused by Babesia microti is an emerging tick-borne zoonosis of increasing importance due to its rising incidence and expanding geographic range(1). Infection with this organism, an intraerythrocytic parasite of the phylum Apicomplexa, causes a febrile syndrome similar to malaria(2). Relapsing disease is common among immunocompromised and asplenic individuals(3,4) and drug resistance has recently been reported(5). To investigate the origin and genetic diversity of this parasite, we sequenced the complete genomes of 42 B. microti samples from around the world, including deep coverage of clinical infections at endemic sites in the continental USA. Samples from the continental USA segregate into a Northeast lineage and a Midwest lineage, with subsequent divergence of subpopulations along geographic lines. We identify parasite variants that associate with relapsing disease, including amino acid substitutions in the atovaquone-binding regions of cytochrome b (cytb) and the azithromycin-binding region of ribosomal protein subunit L4 (rpl4). Our results shed light on the origin, diversity and evolution of B. microti, suggest possible mechanisms for clinical relapse, and create the foundation for further research on this emerging pathogen.

Biomolecular characterization and protein sequences of the Campanian hadrosaur B. canadensis.

Molecular preservation in non-avian dinosaurs is controversial. We present multiple lines of evidence that endogenous proteinaceous material is preserved in bone fragments and soft tissues from an 80-million-year-old Campanian hadrosaur, Brachylophosaurus canadensis [Museum of the Rockies (MOR) 2598]. Microstructural and immunological data are consistent with preservation of multiple bone matrix and vessel proteins, and phylogenetic analyses of Brachylophosaurus collagen sequenced by mass spectrometry robustly support the bird-dinosaur clade, consistent with an endogenous source for these collagen peptides. These data complement earlier results from Tyrannosaurus rex (MOR 1125) and confirm that molecular preservation in Cretaceous dinosaurs is not a unique event.

Molecular phylogenetics of mastodon and Tyrannosaurus rex.

We report a molecular phylogeny for a nonavian dinosaur, extending our knowledge of trait evolution within nonavian dinosaurs into the macromolecular level of biological organization. Fragments of collagen alpha1(I) and alpha2(I) proteins extracted from fossil bones of Tyrannosaurus rex and Mammut americanum (mastodon) were analyzed with a variety of phylogenetic methods. Despite missing sequence data, the mastodon groups with elephant and the T. rex groups with birds, consistent with predictions based on genetic and morphological data for mastodon and on morphological data for T. rex. Our findings suggest that molecular data from long-extinct organisms may have the potential for resolving relationships at critical areas in the vertebrate evolutionary tree that have, so far, been phylogenetically intractable.

A label-free quantification method by MS/MS TIC compared to SILAC and spectral counting in a proteomics screen.

In order to assess the biological function of proteins and their modifications for understanding signaling mechanisms within cells as well as specific biomarkers to disease, it is important that quantitative information be obtained under different experimental conditions. Stable isotope labeling is a powerful method for accurately determining changes in the levels of proteins and PTMs; however, isotope labeling experiments suffer from limited dynamic range resulting in signal change ratios of less than approximately 20:1 using most commercial mass spectrometers. Label-free approaches to relative quantification in proteomics such as spectral counting have gained popularity since no additional chemistries are needed. Here, we show a label-free method for relative quantification based on the TIC from peptide MS/MS spectra collected from data-dependent runs can be used effectively as a quantitative measure and expands the dynamic range over isotope labeling experiments allowing for abundance differences up to approximately 60:1 in a screen for proteins that bind to phosphotyrosine residues.

Protein sequences from mastodon and Tyrannosaurus rex revealed by mass spectrometry.

Fossilized bones from extinct taxa harbor the potential for obtaining protein or DNA sequences that could reveal evolutionary links to extant species. We used mass spectrometry to obtain protein sequences from bones of a 160,000- to 600,000-year-old extinct mastodon (Mammut americanum) and a 68-million-year-old dinosaur (Tyrannosaurus rex). The presence of T. rex sequences indicates that their peptide bonds were remarkably stable. Mass spectrometry can thus be used to determine unique sequences from ancient organisms from peptide fragmentation patterns, a valuable tool to study the evolution and adaptation of ancient taxa from which genomic sequences are unlikely to be obtained.