Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms

DNA sequencing continues to decrease in cost with the Illumina HiSeq2000 generating up to 600 Gb of paired-end 100 base reads in a ten-day run. Here we present a protocol for community amplicon sequencing on the HiSeq2000 and MiSeq Illumina platforms, and apply that protocol to sequence 24 microbial communities from host-associated and free-living environments. A critical question as more sequencing platforms become available is whether biological conclusions derived on one platform are consistent with what would be derived on a different platform. We show that the protocol developed for these instruments successfully recaptures known biological results, and additionally that biological conclusions are consistent across sequencing platforms (the HiSeq2000 versus the MiSeq) and across the sequenced regions of amplicons.

Rapid whole-genome sequencing for investigation of a neonatal MRSA outbreak

BACKGROUND

Isolates of methicillin-resistant Staphylococcus aureus (MRSA) belonging to a single lineage are often indistinguishable by means of current typing techniques. Whole-genome sequencing may provide improved resolution to define transmission pathways and characterize outbreaks.

METHODS

We investigated a putative MRSA outbreak in a neonatal intensive care unit. By using rapid high-throughput sequencing technology with a clinically relevant turnaround time, we retrospectively sequenced the DNA from seven isolates associated with the outbreak and another seven MRSA isolates associated with carriage of MRSA or bacteremia in the same hospital.

RESULTS

We constructed a phylogenetic tree by comparing single-nucleotide polymorphisms (SNPs) in the core genome to a reference genome (an epidemic MRSA clone, EMRSA-15 [sequence type 22]). This revealed a distinct cluster of outbreak isolates and clear separation between these and the nonoutbreak isolates. A previously missed transmission event was detected between two patients with bacteremia who were not part of the outbreak. We created an artificial "resistome" of antibiotic-resistance genes and demonstrated concordance between it and the results of phenotypic susceptibility testing; we also created a "toxome" consisting of toxin genes. One outbreak isolate had a hypermutator phenotype with a higher number of SNPs than the other outbreak isolates, highlighting the difficulty of imposing a simple threshold for the number of SNPs between isolates to decide whether they are part of a recent transmission chain.

CONCLUSIONS

Whole-genome sequencing can provide clinically relevant data within a time frame that can influence patient care. The need for automated data interpretation and the provision of clinically meaningful reports represent hurdles to clinical implementation. (Funded by the U.K. Clinical Research Collaboration Translational Infection Research Initiative and others.).

Rapid whole-genome sequencing for genetic disease diagnosis in neonatal intensive care units

Monogenic diseases are frequent causes of neonatal morbidity and mortality, and disease presentations are often undifferentiated at birth. More than 3500 monogenic diseases have been characterized, but clinical testing is available for only some of them and many feature clinical and genetic heterogeneity. Hence, an immense unmet need exists for improved molecular diagnosis in infants. Because disease progression is extremely rapid, albeit heterogeneous, in newborns, molecular diagnoses must occur quickly to be relevant for clinical decision-making. We describe 50-hour differential diagnosis of genetic disorders by whole-genome sequencing (WGS) that features automated bioinformatic analysis and is intended to be a prototype for use in neonatal intensive care units. Retrospective 50-hour WGS identified known molecular diagnoses in two children. Prospective WGS disclosed potential molecular diagnosis of a severe GJB2-related skin disease in one neonate; BRAT1-related lethal neonatal rigidity and multifocal seizure syndrome in another infant; identified BCL9L as a novel, recessive visceral heterotaxy gene (HTX6) in a pedigree; and ruled out known candidate genes in one infant. Sequencing of parents or affected siblings expedited the identification of disease genes in prospective cases. Thus, rapid WGS can potentially broaden and foreshorten differential diagnosis, resulting in fewer empirical treatments and faster progression to genetic and prognostic counseling.

A culture-independent sequence-based metagenomics approach to the investigation of an outbreak of Shiga-toxigenic Escherichia coli O104:H4

IMPORTANCE

Identification of the bacterium responsible for an outbreak can aid in disease management. However, traditional culture-based diagnosis can be difficult, particularly if no specific diagnostic test is available for an outbreak strain.

OBJECTIVE

To explore the potential of metagenomics, which is the direct sequencing of DNA extracted from microbiologically complex samples, as an open-ended clinical discovery platform capable of identifying and characterizing bacterial strains from an outbreak without laboratory culture.

DESIGN, SETTING, AND PATIENTS

In a retrospective investigation, 45 samples were selected from fecal specimens obtained from patients with diarrhea during the 2011 outbreak of Shiga-toxigenic Escherichia coli (STEC) O104:H4 in Germany. Samples were subjected to high-throughput sequencing (August-September 2012), followed by a 3-phase analysis (November 2012-February 2013). In phase 1, a de novo assembly approach was developed to obtain a draft genome of the outbreak strain. In phase 2, the depth of coverage of the outbreak strain genome was determined in each sample. In phase 3, sequences from each sample were compared with sequences from known bacteria to identify pathogens other than the outbreak strain.

MAIN OUTCOMES AND MEASURES

The recovery of genome sequence data for the purposes of identification and characterization of the outbreak strain and other pathogens from fecal samples.

RESULTS

During phase 1, a draft genome of the STEC outbreak strain was obtained. During phase 2, the outbreak strain genome was recovered from 10 samples at greater than 10-fold coverage and from 26 samples at greater than 1-fold coverage. Sequences from the Shiga-toxin genes were detected in 27 of 40 STEC-positive samples (67%). In phase 3, sequences from Clostridium difficile, Campylobacter jejuni, Campylobacter concisus, and Salmonella enterica were recovered.

CONCLUSIONS AND RELEVANCE

These results suggest the potential of metagenomics as a culture-independent approach for the identification of bacterial pathogens during an outbreak of diarrheal disease. Challenges include improving diagnostic sensitivity, speeding up and simplifying workflows, and reducing costs.

Biogeographic patterns in below-ground diversity in New York City's Central Park are similar to those observed globally

Soil biota play key roles in the functioning of terrestrial ecosystems, however, compared to our knowledge of above-ground plant and animal diversity, the biodiversity found in soils remains largely uncharacterized. Here, we present an assessment of soil biodiversity and biogeographic patterns across Central Park in New York City that spanned all three domains of life, demonstrating that even an urban, managed system harbours large amounts of undescribed soil biodiversity. Despite high variability across the Park, below-ground diversity patterns were predictable based on soil characteristics, with prokaryotic and eukaryotic communities exhibiting overlapping biogeographic patterns. Further, Central Park soils harboured nearly as many distinct soil microbial phylotypes and types of soil communities as we found in biomes across the globe (including arctic, tropical and desert soils). This integrated cross-domain investigation highlights that the amount and patterning of novel and uncharacterized diversity at a single urban location matches that observed across natural ecosystems spanning multiple biomes and continents.

Protein-inspired antibiotics active against vancomycin- and daptomycin-resistant bacteria

The public health threat posed by a looming ‘post-antibiotic’ era necessitates new approaches to antibiotic discovery. Drug development has typically avoided exploitation of membrane-binding properties, in contrast to nature’s control of biological pathways via modulation of membrane-associated proteins and membrane lipid composition. Here, we describe the rejuvenation of the glycopeptide antibiotic vancomycin via selective targeting of bacterial membranes. Peptide libraries based on positively charged electrostatic effector sequences are ligated to N-terminal lipophilic membrane-insertive elements and then conjugated to vancomycin. These modified lipoglycopeptides, the ‘vancapticins’, possess enhanced membrane affinity and activity against methicillin-resistant Staphylococcus aureus (MRSA) and other Gram-positive bacteria, and retain activity against glycopeptide-resistant strains. Optimised antibiotics show in vivo efficacy in multiple models of bacterial infection. This membrane-targeting strategy has potential to ‘revitalise’ antibiotics that have lost effectiveness against recalcitrant bacteria, or enhance the activity of other intravenous-administered drugs that target membrane-associated receptors.

The antibiotic vancomycin inhibits bacterial cell wall synthesis by binding to a membrane-associated precursor. Here, Blaskovich et al. synthesize vancomycin derivatives containing lipophilic peptide moieties that enhance membrane affinity and in vivo activities against glycopeptide-resistant strains.

Deep Sequencing of B Cell Receptor Repertoires From COVID-19 Patients Reveals Strong Convergent Immune Signatures

Deep sequencing of B cell receptor (BCR) heavy chains from a cohort of 31 COVID-19 patients from the UK reveals a stereotypical naive immune response to SARS-CoV-2 which is consistent across patients. Clonal expansion of the B cell population is also observed and may be the result of memory bystander effects. There was a strong convergent sequence signature across patients, and we identified 1,254 clonotypes convergent between at least four of the COVID-19 patients, but not present in healthy controls or individuals following seasonal influenza vaccination. A subset of the convergent clonotypes were homologous to known SARS and SARS-CoV-2 spike protein neutralizing antibodies. Convergence was also demonstrated across wide geographies by comparison of data sets between patients from UK, USA, and China, further validating the disease association and consistency of the stereotypical immune response even at the sequence level. These convergent clonotypes provide a resource to identify potential therapeutic and prophylactic antibodies and demonstrate the potential of BCR profiling as a tool to help understand patient responses.

Evidence for Community Transmission of Community-Associated but Not Health-Care-Associated Methicillin-Resistant Staphylococcus Aureus Strains Linked to Social and Material Deprivation: Spatial Analysis of Cross-sectional Data

BACKGROUND

Identifying and tackling the social determinants of infectious diseases has become a public health priority following the recognition that individuals with lower socioeconomic status are disproportionately affected by infectious diseases. In many parts of the world, epidemiologically and genotypically defined community-associated (CA) methicillin-resistant Staphylococcus aureus (MRSA) strains have emerged to become frequent causes of hospital infection. The aim of this study was to use spatial models with adjustment for area-level hospital attendance to determine the transmission niche of genotypically defined CA- and health-care-associated (HA)-MRSA strains across a diverse region of South East London and to explore a potential link between MRSA carriage and markers of social and material deprivation.

METHODS AND FINDINGS

This study involved spatial analysis of cross-sectional data linked with all MRSA isolates identified by three National Health Service (NHS) microbiology laboratories between 1 November 2011 and 29 February 2012. The cohort of hospital-based NHS microbiology diagnostic services serves 867,254 usual residents in the Lambeth, Southwark, and Lewisham boroughs in South East London, United Kingdom (UK). Isolates were classified as HA- or CA-MRSA based on whole genome sequencing. All MRSA cases identified over 4 mo within the three-borough catchment area (n = 471) were mapped to small geographies and linked to area-level aggregated socioeconomic and demographic data. Disease mapping and ecological regression models were used to infer the most likely transmission niches for each MRSA genetic classification and to describe the spatial epidemiology of MRSA in relation to social determinants. Specifically, we aimed to identify demographic and socioeconomic population traits that explain cross-area extra variation in HA- and CA-MRSA relative risks following adjustment for hospital attendance data. We explored the potential for associations with the English Indices of Deprivation 2010 (including the Index of Multiple Deprivation and several deprivation domains and subdomains) and the 2011 England and Wales census demographic and socioeconomic indicators (including numbers of households by deprivation dimension) and indicators of population health. Both CA-and HA-MRSA were associated with household deprivation (CA-MRSA relative risk [RR]: 1.72 [1.03-2.94]; HA-MRSA RR: 1.57 [1.06-2.33]), which was correlated with hospital attendance (Pearson correlation coefficient [PCC] = 0.76). HA-MRSA was also associated with poor health (RR: 1.10 [1.01-1.19]) and residence in communal care homes (RR: 1.24 [1.12-1.37]), whereas CA-MRSA was linked with household overcrowding (RR: 1.58 [1.04-2.41]) and wider barriers, which represent a combined score for household overcrowding, low income, and homelessness (RR: 1.76 [1.16-2.70]). CA-MRSA was also associated with recent immigration to the UK (RR: 1.77 [1.19-2.66]). For the area-level variation in RR for CA-MRSA, 28.67% was attributable to the spatial arrangement of target geographies, compared with only 0.09% for HA-MRSA. An advantage to our study is that it provided a representative sample of usual residents receiving care in the catchment areas. A limitation is that relationships apparent in aggregated data analyses cannot be assumed to operate at the individual level.

CONCLUSIONS

There was no evidence of community transmission of HA-MRSA strains, implying that HA-MRSA cases identified in the community originate from the hospital reservoir and are maintained by frequent attendance at health care facilities. In contrast, there was a high risk of CA-MRSA in deprived areas linked with overcrowding, homelessness, low income, and recent immigration to the UK, which was not explainable by health care exposure. Furthermore, areas adjacent to these deprived areas were themselves at greater risk of CA-MRSA, indicating community transmission of CA-MRSA. This ongoing community transmission could lead to CA-MRSA becoming the dominant strain types carried by patients admitted to hospital, particularly if successful hospital-based MRSA infection control programmes are maintained. These results suggest that community infection control programmes targeting transmission of CA-MRSA will be required to control MRSA in both the community and hospital. These epidemiological changes will also have implications for effectiveness of risk-factor-based hospital admission MRSA screening programmes.

Microbial Metabolism Shifts Towards an Adverse Profile with Supplementary Iron in the TIM-2 In vitro Model of the Human Colon

Oral iron administration in African children can increase the risk for infections. However, it remains unclear to what extent supplementary iron affects the intestinal microbiome. We here explored the impact of iron preparations on microbial growth and metabolism in the well-controlled TNO's in vitro model of the large intestine (TIM-2). The model was inoculated with a human microbiota, without supplementary iron, or with 50 or 250 μmol/L ferrous sulfate, 50 or 250 μmol/L ferric citrate, or 50 μmol/L hemin. High resolution responses of the microbiota were examined by 16S rDNA pyrosequencing, microarray analysis, and metagenomic sequencing. The metabolome was assessed by fatty acid quantification, gas chromatography-mass spectrometry (GC-MS), and ¹H-NMR spectroscopy. Cultured intestinal epithelial Caco-2 cells were used to assess fecal water toxicity. Microbiome analysis showed, among others, that supplementary iron induced decreased levels of Bifidobacteriaceae and Lactobacillaceae, while it caused higher levels of Roseburia and Prevotella. Metagenomic analyses showed an enrichment of microbial motility-chemotaxis systems, while the metabolome markedly changed from a saccharolytic to a proteolytic profile in response to iron. Branched chain fatty acids and ammonia levels increased significantly, in particular with ferrous sulfate. Importantly, the metabolite-containing effluent from iron-rich conditions showed increased cytotoxicity to Caco-2 cells. Our explorations indicate that in the absence of host influences, iron induces a more hostile environment characterized by a reduction of microbes that are generally beneficial, and increased levels of bacterial metabolites that can impair the barrier function of a cultured intestinal epithelial monolayer.

Solid-phase oligodeoxynucleotide synthesis: a two-step cycle using peroxy anion deprotection

A novel solid-phase phosphoramidite based oligodeoxynucleotide two-step synthesis method has been developed. Keys to this method are replacement of the 5'-dimethoxytrityl blocking group with an aryloxycarbonyl and the use of N-dimethoxytrityl protection for the exocyclic amines of adenine and cytosine. With these modifications, coupling of each 2'-deoxynucleoside 3'-phosphoramidite to the growing oligodeoxynucleotide on the solid support can be followed by treatment with an aqueous mixture of peroxy anions buffered at pH 9.6. This reagent effectively removes the carbonate protecting group and simultaneously oxidizes the phosphite internucleotide linkage. As a consequence a new two-step synthesis cycle is possible. Oligodeoxynucleotides synthesized using this approach are identical to authentic samples when tested by a variety of analytical techniques.

Frequent Undetected Ward-Based Methicillin-Resistant Staphylococcus aureus Transmission Linked to Patient Sharing Between Hospitals

Background

Recent evidence suggests that hospital transmission of methicillin-resistant Staphylococcus aureus (MRSA) is uncommon in UK centers that have implemented sustained infection control programs. We investigated whether a healthcare-network analysis could shed light on transmission paths currently sustaining MRSA levels in UK hospitals.

Methods

A cross-sectional observational study was performed in 2 National Health Service hospital groups and a general district hospital in Southeast London. All MRSA patients identified at inpatient, outpatient, and community settings between 1 November 2011 and 29 February 2012 were included. We identified genetically defined MRSA transmission clusters in individual hospitals and across the healthcare network, and examined genetic differentiation of sequence type (ST) 22 MRSA isolates within and between hospitals and inpatient or outpatient and community settings, as informed by average and median pairwise single-nucleotide polymorphisms (SNPs) and SNP-based proportions of nearly identical isolates.

Results

Two hundred forty-eight of 610 (40.7%) MRSA patients were linked in 90 transmission clusters, of which 27 spanned multiple hospitals. Analysis of a large 32 patient ST22-MRSA cluster showed that 26 of 32 patients (81.3%) had multiple contacts with one another during ward stays at any hospital. No residential, outpatient, or significant community healthcare contacts were identified. Genetic differentiation between ST22 MRSA inpatient isolates from different hospitals was less than between inpatient isolates from the same hospitals (P ≤ .01).

Conclusions

There is evidence of frequent ward-based transmission of MRSA brought about by frequent patient admissions to multiple hospitals. Limiting in-ward transmission requires sharing of MRSA status data between hospitals.

Clonal variation in high- and low-level phenotypic and genotypic mupirocin resistance of MRSA isolates in south-east London

OBJECTIVES

Both low-level mupirocin resistance (LMR) and high-level mupirocin resistance (HMR) have been identified. The aim of this study was to determine the epidemiology of LMR and HMR in MRSA isolates at five hospitals that have used mupirocin for targeted decolonization as part of successful institutional control programmes.

METHODS

All MRSA identified in three microbiology laboratories serving five central and south-east London hospitals and surrounding communities between November 2011 and February 2012 were included. HMR and LMR were determined by disc diffusion testing. WGS was used to derive multilocus sequence types (MLSTs) and the presence of HMR and LMR resistance determinants.

RESULTS

Prevalence of either HMR or LMR amongst first healthcare episode isolates from 795 identified patients was 9.69% (95% CI 7.72-11.96); LMR was 6.29% (95% CI 4.70-8.21) and HMR was 3.40% (95% CI 2.25-4.90). Mupirocin resistance was not significantly different in isolates identified from inpatients at each microbiology laboratory, but was more common in genotypically defined 'hospital' rather than 'community' isolates (OR 3.17, 95% CI 1.36-9.30, P = 0.002). LMR was associated with inpatient stay, previous history of MRSA and age ≥65 years; HMR was associated with age ≥65 years and residential postcode outside London. LMR and HMR varied by clone, with both being low in the dominant UK MRSA clone ST22 compared with ST8, ST36 and ST239/241 for LMR and with ST8 and ST36 for HMR. V588F mutation and mupA carriage had high specificity (>97%) and area under the curve (>83%) to discriminate phenotypic mupirocin resistance, but uncertainty around the sensitivity point estimate was large (95% CI 52.50%-94.44%). Mutations in or near the mupA gene were found in eight isolates that carried mupA but were not HMR.

CONCLUSIONS

Mupirocin resistance was identified in <10% of patients and varied significantly by clone, implying that changes in clonal epidemiology may have an important role in determining the prevalence of resistance in conjunction with selection due to mupirocin use.

Synthesis and use of a pseudo-cysteine for native chemical ligation

The process of native chemical ligation (NCL) is well described in the literature. An N-terminal cysteine-containing peptide reacts with a C-terminal thioester-containing peptide to yield a native amide bond after transesterification and acyl transfer. An N-terminal cysteine is required as both the N-terminal amino function and the sidechain thiol participate in the ligation reaction. In certain circumstances it is desirable, or even imperative, that the N-terminal region of a peptidic reaction partner remain unmodified, for Instance for the retention of biological activity after ligation. This work discusses the synthesis of a pseudo-N-terminal cysteine building block for incorporation into peptides using standard methods of solid phase synthesis. Upon deprotection, this building block affords a de facto N-terminal cysteine positioned on an amino acid sidechain. which is capable of undergoing native chemical ligation with a thioester. The syntheses of several peptides and structures containing this motif are detailed, their reactions discussed. and further applications of this technology proposed.

Synthesis of DNA using a new two-step cycle

For the first time a new, two-step method is described for synthesizing deoxyribonucleic acid. The approach uses 5'-carbonate protected 2'-deoxynucleoside-3'-phosphoramidites as synthons and a peroxy anion buffer that removes the carbonate protecting group and oxidizes the internucleotide linkage. Following synthesis via this two-step cycle, oligomers are isolated by standard procedures.