Exploring the pediatric nasopharyngeal bacterial microbiota with culture-based MALDI-TOF mass spectrometry and targeted metagenomic sequencing

The nasopharynx is an important reservoir of disease-associated and antimicrobial-resistant bacterial species. This proof-of-concept study assessed the utility of a combined culture, matrix-assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOF MS), and targeted metagenomic sequencing workflow for the study of the pediatric nasopharyngeal bacterial microbiota. Nasopharyngeal swabs and clinical metadata were collected from Cambodian children during a hospital outpatient visit and then biweekly for 12 weeks. Swabs were cultured on chocolate and blood-gentamicin agar, and all colony morphotypes were identified by MALDI-TOF MS. Metagenomic sequencing was done on a scrape of all colonies from a chocolate agar culture and processed using the mSWEEP pipeline. One hundred one children were enrolled, yielding 620 swabs. MALDI-TOF MS identified 106 bacterial species/40 genera: 20 species accounted for 88.5% (2,190/2,474) of isolates. Colonization by Moraxella catarrhalis (92.1% of children on ≥1 swab), Haemophilus influenzae (87.1%), and Streptococcus pneumoniae (83.2%) was particularly common. In S. pneumoniae-colonized children, a median of two serotypes [inter-quartile range (IQR) 1-2, range 1-4] was detected. For the 21 bacterial species included in the mSWEEP database and identifiable by MALDI-TOF, detection by culture + MALDI-TOF MS and culture + mSWEEP was highly concordant with a median species-level agreement of 96.9% (IQR 86.8%-98.8%). mSWEEP revealed highly dynamic lineage-level colonization patterns for S. pneumoniae which were quite different to those for S. aureus. A combined culture, MALDI-TOF MS, targeted metagenomic sequencing approach for the exploration of the young child nasopharyngeal microbiome was technically feasible, and each component yielded complementary data.

IMPORTANCE

The human upper respiratory tract is an important source of disease-causing and antibiotic-resistant bacteria. However, understanding the interactions and stability of these bacterial populations is technically challenging. We used a combination of approaches to determine colonization patterns over a 3-month period in 101 Cambodian children. The combined approach was feasible to implement, and each component gave complementary data to enable a better understanding of the complex patterns of bacterial colonization.

Altered oral microbiome, but normal human papilloma virus prevalence in cartilage-hair hypoplasia patients

BACKGROUND

Cartilage-hair hypoplasia (CHH) is a rare syndromic immunodeficiency with metaphyseal chondrodysplasia and increased risk of malignancy. In this cross-sectional observational study, we examined HPV status and oral microbiome in individuals with CHH. Oral brush samples were collected from 20 individuals with CHH (aged 5-59 years) and 41 controls (1-69 years). Alpha HPVs (43 types) were tested by nested PCR followed by bead-based probe hybridization. Separately, beta-, gamma-, mu- and nu- HPV types were investigated, and a genome-based bacterial microbiome sequencing was performed.

RESULTS

We found a similar alpha HPV prevalence in individuals with CHH (45%) and controls (36%). The HPV types of individuals with CHH were HPV-16 (25%), 27, 28, and 78, and of controls HPV-3, 16 (21%), 27, and 61. Beta HPV positivity and combined beta/gamma/mu/nu prevalence was detected in 11% and 11% of individuals with CHH and in 5% and 3% of the controls, respectively. Individuals with CHH differed from the controls in bacterial microbiota diversity, richness, and in microbial composition. Individuals with CHH had lower abundance of species Mitsuokella sp000469545, Parascardovia denticolens, Propionibacterium acidifaciens, UMGS1907 sp004151455, Salinicola halophilus, Haemophilus_A paraphrohaemolyticus, Fusobacterium massiliense, and Veillonella parvula, and higher abundance of Slackia exigua.

CONCLUSIONS

Individuals with CHH exhibit similar prevalence of HPV DNA but different bacterial microbiota on their oral mucosa compared to healthy controls. This may partly explain the previously observed high prevalence of oral diseases in CHH, and regular oral examination is warranted.

Differences in Gut Microbiota Profiles and Microbiota Steroid Hormone Biosynthesis in Men with and Without Prostate Cancer

Background

Although prostate cancer (PCa) is the most common cancer in men in Western countries, there is significant variability in geographical incidence. This might result from genetic factors, discrepancies in screening policies, or differences in lifestyle. Gut microbiota has recently been associated with cancer progression, but its role in PCa is unclear.

Objective

Characterization of the gut microbiota and its functions associated with PCa.

Design setting and participants

In a prospective multicenter clinical trial (NCT02241122), the gut microbiota profiles of 181 men with a clinical suspicion of PCa were assessed utilizing 16S rRNA sequencing.

Outcome measurements and statistical analysis

Sequences were assigned to operational taxonomic units, differential abundance analysis, and α- and β-diversities, and predictive functional analyses were performed. Plasma steroid hormone levels corresponding to the predicted microbiota steroid hormone biosynthesis profiles were investigated.

Results and limitations

Of 364 patients, 181 were analyzed, 60% of whom were diagnosed with PCa. Microbiota composition and diversity were significantly different in PCa, partially affected by Prevotella 9, the most abundant genus of the cohort, and significantly higher in PCa patients. Predictive functional analyses revealed higher 5-α-reductase, copper absorption, and retinol metabolism in the PCa-associated microbiome. Plasma testosterone was associated negatively with the predicted microbial 5-α-reductase level.

Conclusions

Gut microbiota of the PCa patients differed significantly compared with benign individuals. Microbial 5-α-reductase, copper absorption, and retinol metabolism are potential mechanisms of action. These findings support the observed association of lifestyle, geography, and PCa incidence.

Patient summary

In this report, we found that several microbes and potential functions of the gut microbiota are altered in prostate cancer compared with benign cases. These findings suggest that gut microbiota could be the link between environmental factors and prostate cancer.

Modulation of multidrug-resistant clone success in Escherichia coli populations: a longitudinal, multi-country, genomic and antibiotic usage cohort study

BACKGROUND

The effect of antibiotic usage on the success of multidrug-resistant (MDR) clones in a population remains unclear. With this genomics-based molecular epidemiology study, we aimed to investigate the contribution of antibiotic use to Escherichia coli clone success, relative to intra-strain competition for colonisation and infection.

METHODS

We sequenced all the available E coli bloodstream infection isolates provided by the British Society for Antimicrobial Chemotherapy (BSAC) from 2012 to 2017 (n=718) and combined these with published data from the UK (2001-11; n=1090) and Norway (2002-17; n=3254). Defined daily dose (DDD) data from the European Centre for Disease Prevention and Control (retrieved on Sept 21, 2021) for major antibiotic classes (β-lactam, tetracycline, macrolide, sulfonamide, quinolone, and non-penicillin β-lactam) were used together with sequence typing, resistance profiling, regression analysis, and non-neutral Wright-Fisher simulation-based modelling to enable systematic comparison of resistance levels, clone success, and antibiotic usage between the UK and Norway.

FINDINGS

Sequence type (ST)73, ST131, ST95, and ST69 accounted for 892 (49·3%) of 1808 isolates in the BSAC collection. In the UK, the proportion of ST69 increased between 2001-10 and 2011-17 (p=0·0004), whereas the proportions of ST73 and ST95 did not vary between periods. ST131 expanded quickly after its emergence in 2003 and its prevalence remained consistent throughout the study period (apart from a brief decrease in 2009-10). The extended-spectrum β-lactamase (ESBL)-carrying, globally disseminated MDR clone ST131-C2 showed overall greater success in the UK (154 [56·8%] of 271 isolates in 2003-17) compared with Norway (51 [18·3%] of 278 isolates in 2002-17; p<0·0001). DDD data indicated higher total use of antimicrobials in the UK, driven mainly by the class of non-penicillin β-lactams, which were used between 2·7-times and 5·1-times more in the UK per annum (ratio mean 3·7 [SD 0·8]). This difference was associated with the higher success of the MDR clone ST131-C2 (pseudo-R2 69·1%). A non-neutral Wright-Fisher model replicated the observed expansion of non-MDR and MDR sequence types under higher DDD regimes.

INTERPRETATION

Our study indicates that resistance profiles of contemporaneously successful clones can vary substantially, warranting caution in the interpretation of correlations between aggregate measures of resistance and antibiotic usage. Our study further suggests that in countries with low-to-moderate use of antibiotics, such as the UK and Norway, the extent of non-penicillin β-lactam use modulates rather than determines the success of widely disseminated MDR ESBL-carrying E coli clones. Detailed understanding of underlying causal drivers of success is important for improved control of resistant pathogens.

FUNDING

Trond Mohn Foundation, Marie Skłodowska-Curie Actions, European Research Council, Royal Society, and Wellcome Trust.

Optimal sampling and analysis methods for clinical diagnostics of vaginal microbiome

Next-generation sequencing-based microbiological analysis is a complex way to profile vaginal microbiome samples since each step affects the results gained. Methodologies for sample collection lack golden standards. We compared Puritan DNA/RNA swab (PS) and Copan FLOQ swab (CS) and provided consistent and reliable microbiome profiles analyzed by 16S rRNA gene sequencing. We collected two consecutive vaginal samples utilizing PS with room temperature storing and CS with instant freezing from 26 women. Variable region 4 of bacterial 16S rRNA gene was amplified with single PCR by custom-designed dual-indexed primers and sequenced with Illumina MiSeq system. Read quality control, operational taxonomic unit tables, and alpha and beta diversities analysis were performed, and community richness, diversity, and evenness were evaluated and compared between the two samplings and tests. Nineteen sample pairs produced detectable, intact DNA during the extraction protocol and/or further microbial profiles. Alpha bacterial diversity indices were independent on the collection protocol. No significant statistical differences were found in the measured beta diversity metrics between the collection methods. Of the women, 43% had Lactobacillus-dominated vaginal microbiome profile despite of collection method. Previously reported important vaginal microbiome phyla Actinobacteria, Bacteroidetes, Firmicutes, Fusobacteria, and Proteobacteria were present in the sample set although their relative abundances varied among individuals. PS and CS enable constant vaginal microbiota sampling. The PS method with no need for instant freezing is suitable for on-site collections at clinics. Furthermore, it seems to be possible to take two samples instead of one with constant microbiological results.

Appendiceal microbiome in uncomplicated and complicated acute appendicitis: A prospective cohort study

Background

Uncomplicated and complicated acute appendicitis seem to be two different forms of this common abdominal emergency. The contribution of appendiceal microbiota to appendicitis pathogenesis has been suggested, but differences between uncomplicated and complicated appendicitis are largely unknown. We compared the appendiceal microbiota in uncomplicated and complicated acute appendicitis.

Methods

This prospective single-center clinical cohort study was conducted as part of larger multicenter MAPPAC trial enrolling adult patients with computed tomography or clinically confirmed uncomplicated or complicated acute appendicitis. The microbial composition of the appendiceal lumen was determined using 16S rRNA gene amplicon sequencing.

Results

Between April 11, 2017, and March 29, 2019, 118 samples (41 uncomplicated and 77 complicated appendicitis) were available. After adjusting for age, sex, and BMI, alpha diversity in complicated appendicitis was higher (Shannon p = 0.011, Chao1 p = 0.006) compared to uncomplicated appendicitis. Microbial compositions were different between uncomplicated and complicated appendicitis (Bray-Curtis distance, P = 0.002). Species poor appendiceal microbiota composition with specific predominant bacteria was present in some patients regardless of appendicitis severity.

Conclusion

Uncomplicated and complicated acute appendicitis have different appendiceal microbiome profiles further supporting the disconnection between these two different forms of acute appendicitis.

Study registration

ClinicalTrials.gov NCT03257423.

Gut Microbiota and Serum Metabolome in Elite Cross-Country Skiers: A Controlled Study

Exercise has been shown to affect gut the microbiome and metabolic health, with athletes typically displaying a higher microbial diversity. However, research on the gut microbiota and systemic metabolism in elite athletes remains scarce. In this study, we compared the gut microbiota profiles and serum metabolome of national team cross-country skiers at the end of an exhausting training and competitive season to those of normally physically-active controls. The gut microbiota were analyzed using 16S rRNA amplicon sequencing. Serum metabolites were analyzed using nuclear magnetic resonance. Phylogenetic diversity and the abundance of several mucin-degrading gut microbial taxa, including Akkermansia, were lower in the athletes. The athletes had a healthier serum lipid profile than the controls, which was only partly explained by body mass index. Butyricicoccus associated positively with HDL cholesterol, HDL2 cholesterol and HDL particle size. The Ruminococcus torques group was less abundant in the athlete group and positively associated with total cholesterol and VLDL and LDL particles. We found the healthier lipid profile of elite athletes to co-occur with known health-beneficial gut microbes. Further studies should elucidate these links and whether athletes are prone to mucin depletion related microbial changes during the competitive season.

External Quality Assessment of SARS-CoV-2 Sequencing: an ESGMD-SSM Pilot Trial across 15 European Laboratories

This first pilot trial on external quality assessment (EQA) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) whole-genome sequencing, initiated by the European Society of Clinical Microbiology and Infectious Diseases (ESCMID) Study Group for Genomic and Molecular Diagnostics (ESGMD) and the Swiss Society for Microbiology (SSM), aims to build a framework between laboratories in order to improve pathogen surveillance sequencing. Ten samples with various viral loads were sent out to 15 clinical laboratories that had free choice of sequencing methods and bioinformatic analyses. The key aspects on which the individual centers were compared were the identification of (i) single nucleotide polymorphisms (SNPs) and indels, (ii) Pango lineages, and (iii) clusters between samples. The participating laboratories used a wide array of methods and analysis pipelines. Most were able to generate whole genomes for all samples. Genomes were sequenced to various depths (up to a 100-fold difference across centers). There was a very good consensus regarding the majority of reporting criteria, but there were a few discrepancies in lineage and cluster assignments. Additionally, there were inconsistencies in variant calling. The main reasons for discrepancies were missing data, bioinformatic choices, and interpretation of data. The pilot EQA was overall a success. It was able to show the high quality of participating laboratories and provide valuable feedback in cases where problems occurred, thereby improving the sequencing setup of laboratories. A larger follow-up EQA should, however, improve on defining the variables and format of the report. Additionally, contamination and/or minority variants should be a further aspect of assessment.

A large-scale genomic snapshot of Klebsiella spp. isolates in Northern Italy reveals limited transmission between clinical and non-clinical settings

The Klebsiella group, found in humans, livestock, plants, soil, water and wild animals, is genetically and ecologically diverse. Many species are opportunistic pathogens and can harbour diverse classes of antimicrobial resistance genes. Healthcare-associated Klebsiella pneumoniae clones that are non-susceptible to carbapenems can spread rapidly, representing a high public health burden. Here we report an analysis of 3,482 genome sequences representing 15 Klebsiella species sampled over a 17-month period from a wide range of clinical, community, animal and environmental settings in and around the Italian city of Pavia. Northern Italy is a hotspot for hospital-acquired carbapenem non-susceptible Klebsiella and thus a pertinent setting to examine the overlap between isolates in clinical and non-clinical settings. We found no genotypic or phenotypic evidence for non-susceptibility to carbapenems outside the clinical environment. Although we noted occasional transmission between clinical and non-clinical settings, our data point to a limited role of animal and environmental reservoirs in the human acquisition of Klebsiella spp. We also provide a detailed genus-wide view of genomic diversity and population structure, including the identification of new groups.

Bacterial genomic epidemiology with mixed samples

Genomic epidemiology is a tool for tracing transmission of pathogens based on whole-genome sequencing. We introduce the mGEMS pipeline for genomic epidemiology with plate sweeps representing mixed samples of a target pathogen, opening the possibility to sequence all colonies on selective plates with a single DNA extraction and sequencing step. The pipeline includes the novel mGEMS read binner for probabilistic assignments of sequencing reads, and the scalable pseudoaligner Themisto. We demonstrate the effectiveness of our approach using closely related samples in a nosocomial setting, obtaining results that are comparable to those based on single-colony picks. Our results lend firm support to more widespread consideration of genomic epidemiology with mixed infection samples.

High-resolution sweep metagenomics using fast probabilistic inference

Determining the composition of bacterial communities beyond the level of a genus or species is challenging because of the considerable overlap between genomes representing close relatives. Here, we present the mSWEEP pipeline for identifying and estimating the relative sequence abundances of bacterial lineages from plate sweeps of enrichment cultures. mSWEEP leverages biologically grouped sequence assembly databases, applying probabilistic modelling, and provides controls for false positive results. Using sequencing data from major pathogens, we demonstrate significant improvements in lineage quantification and detection accuracy. Our pipeline facilitates investigating cultures comprising mixtures of bacteria, and opens up a new field of plate sweep metagenomics.

Dynamics of intestinal multidrug-resistant bacteria colonisation contracted by visitors to a high-endemic setting: a prospective, daily, real-time sampling study

BACKGROUND

Antimicrobial resistance is highly prevalent in low-income and middle-income countries. International travel contributes substantially to the global spread of intestinal multidrug-resistant Gram-negative bacteria. Hundreds of millions of annual visitors to low-income and middle-income countries are all exposed to intestinal multidrug-resistant Gram-negative bacteria resulting in 30-70% of them being colonised at their return. The colonisation process in high-exposure environments is poorly documented because data have only been derived from before travel and after travel sampling. We characterised colonisation dynamics by exploring daily stool samples while visiting a low-income and middle-income countries.

METHODS

In this prospective, daily, real-time sampling study 20 European visitors to Laos volunteered to provide daily stool samples and completed daily questionnaires for 22 days. Samples were initially assessed at Mahosot Hospital, Vientiane, Laos, for acquisition of extended-spectrum β-lactamase-producing (ESBL) Gram-negative bacteria followed by whole-genome sequencing of isolates at MicrobesNG, University of Birmingham, Birmingham, UK. The primary outcome of the study was to obtain data on the dynamics of intestinal multidrug-resistant bacteria acquisition.

FINDINGS

Between Sept 18 and Sept 20, 2015, 23 volunteers were recruited, of whom 20 (87%) European volunteers were included in the final study population. Although colonisation rates were 70% at the end of the study, daily sampling revealed that all participants had acquired ESBL-producing Gram-negative bacteria at some point during the study period; the colonisation status varied day by day. Whole-genome sequencing analysis ascribed the transient pattern of colonisation to sequential acquisition of new strains, resulting in a loss of detectable colonisation by the initial multidrug-resistant Gram-negative strains. 19 (95%) participants acquired two to seven strains. Of the 83 unique strains identified (53 Escherichia coli, 10 Klebsiella spp, and 20 other ESBL-producing Gram-negative bacteria), some were shared by as many as four (20%) participants.

INTERPRETATION

To our knowledge, this is the first study to characterise in real-time the dynamics of acquiring multidrug-resistant Gram-negative bacterial colonisation during travel. Our data show multiple transient colonisation events indicative of constant microbial competition and suggest that travellers are exposed to a greater burden of multidrug-resistant bacteria than previously thought. The data emphasise the need for preventing travellers' diarrhoea and limiting antibiotic use, addressing the two major factors predisposing colonisation.

FUNDING

The Finnish Governmental Subsidy for Health Science Research, The Scandinavian Society for Antimicrobial Chemotherapy, the Sigrid Jusélius Foundation, Biotechnology and Biological Sciences Research Council; Wellcome Trust, Medical Research Council; The Royal Society; Joint Programming Initiative on Antimicrobial Resistance, and European Research Council.

Genome of Superficieibacter maynardsmithii, a novel, antibiotic susceptible representative of Enterobacteriaceae

During a citywide microbiological screening project in Pavia (Italy) a bacterial strain isolated from the surface of an Automated Teller Machine was classified as a Klebsiella sp. by MALDI-TOF spectrometry, and shown to be susceptible to the most antimicrobial classes by phenotypic testing. After Illumina genome sequencing and subsequent assembly, a high-quality draft genome was obtained (size = 5,051,593 bp, N50 = 615,571 bp, largest contig = 1,328,029 bp, N_contig = 17, GC content = 51.58%, coverage = 141.42), absence of antimicrobial resistance genes was confirmed, but the strain resulted to be highly divergent from all Klebsiella, and more related to other Enterobacteriaceae. The higher values of 16S rRNA identity were with members of the genera Citrobacter, Salmonella, and “Superficieibacter.” An ortholog-based phylogenomic analysis indicated a sister group relationship with “Superficieibacter electus,” in a distinct clade from other members of the Enterobacteriaceae family. In order to evaluate whether the novel genome represents a new species of “Superficiebacter,” average nucleotide identity (ANI) and Hadamard analysis were performed on a dataset of 78 Enterobacteriaceae. The novel genome showed an ANI of 87.51% with S. electus, which compared on identity values between other members of the family, clearly indicates that the genome represents a new species within the genus “Superficieibacter.” We propose for the new species the name “Superficieibacter maynardsmithii.”

The binding mechanism of the virulence factor Streptococcus suis adhesin P subtype to globotetraosylceramide is associated with systemic disease

Streptococcus suis is part of the pig commensal microbiome but strains can also be pathogenic, causing pneumonia and meningitis in pigs as well as zoonotic meningitis. According to genomic analysis, S. suis is divided into asymptomatic carriage, respiratory and systemic strains with distinct genomic signatures. Because the strategies to target pathogenic S. suis are limited, new therapeutic approaches are needed. The virulence factor S. suis adhesin P (SadP) recognizes the galabiose Galα1-4Gal-oligosaccharide. Based on its oligosaccharide fine specificity, SadP can be divided into subtypes P_N and P_O We show here that subtype P_N is distributed in the systemic strains causing meningitis, whereas type P_O is found in asymptomatic carriage and respiratory strains. Both types of SadP are shown to predominantly bind to pig lung globotriaosylceramide (Gb3). However, SadP adhesin from systemic subtype P_N strains also binds to globotetraosylceramide (Gb4). Mutagenesis studies of the galabiose-binding domain of type P_N SadP adhesin showed that the amino acid asparagine 285, which is replaced by an aspartate residue in type P_O SadP, was required for binding to Gb4 and, strikingly, was also required for interaction with the glycomimetic inhibitor phenylurea-galabiose. Molecular dynamics simulations provided insight into the role of Asn-285 for Gb4 and phenylurea-galabiose binding, suggesting additional hydrogen bonding to terminal GalNAc of Gb4 and the urea group. Thus, the Asn-285-mediated molecular mechanism of type P_N SadP binding to Gb4 could be used to selectively target S. suis in systemic disease without interfering with commensal strains, opening up new avenues for interventional strategies against this pathogen.

High-resolution sweep metagenomics using fast probabilistic inference

Determining the composition of bacterial communities beyond the level of a genus or species is challenging because of the considerable overlap between genomes representing close relatives. Here, we present the mSWEEP pipeline for identifying and estimating the relative sequence abundances of bacterial lineages from plate sweeps of enrichment cultures. mSWEEP leverages biologically grouped sequence assembly databases, applying probabilistic modelling, and provides controls for false positive results. Using sequencing data from major pathogens, we demonstrate significant improvements in lineage quantification and detection accuracy. Our pipeline facilitates investigating cultures comprising mixtures of bacteria, and opens up a new field of plate sweep metagenomics.

Genomic analysis of Klebsiella pneumoniae isolates from Malawi reveals acquisition of multiple ESBL determinants across diverse lineages

OBJECTIVES

ESBL-producing Klebsiella pneumoniae (KPN) pose a major threat to human health globally. We carried out a WGS study to understand the genetic background of ESBL-producing KPN in Malawi and place them in the context of other global isolates.

METHODS

We sequenced genomes of 72 invasive and carriage KPN isolates collected from patients admitted to Queen Elizabeth Central Hospital, Blantyre, Malawi. We performed phylogenetic and population structure analyses on these and previously published genomes from Kenya (n = 66) and from outside sub-Saharan Africa (n = 67). We screened for presence of antimicrobial resistance (AMR) genetic determinants and carried out association analyses by genomic sequence cluster, AMR phenotype and time.

RESULTS

Malawian isolates fit within the global population structure of KPN, clustering into the major lineages of KpI, KpII and KpIII. KpI isolates from Malawi were more related to those from Kenya, with both collections exhibiting more clonality than isolates from the rest of the world. We identified multiple ESBL genes, including blaCTX-M-15, several blaSHV, blaTEM-63 and blaOXA-10, and other AMR genes, across diverse lineages of the KPN isolates from Malawi. No carbapenem resistance genes were detected; however, we detected IncFII and IncFIB plasmids that were similar to the carbapenem resistance-associated plasmid pNDM-mar.

CONCLUSIONS

There are multiple ESBL genes across diverse KPN lineages in Malawi and plasmids in circulation that are capable of carrying carbapenem resistance. Unless appropriate interventions are rapidly put in place, these may lead to a high burden of locally untreatable infection in vulnerable populations.

Emergence of ceftazidime-avibactam-resistant Klebsiella pneumoniae during treatment, Finland, December 2018

In December 2018, a ceftazidime-avibactam (CAZ-AVI)-resistant KPC-2-producing Klebsiella pneumoniae strain was isolated in Finland. CAZ-AVI resistance was observed 34 days after CAZ-AVI treatment in a trauma patient transferred from a hospital in Greece who had been colonised with blaKPC-2-producing K. pneumoniae ST39, and later developed a bloodstream infection. The CAZ-AVI-resistant strain contained a novel 15 amino acid insertion in the KPC-2 protein causing structural changes proximal to the KPC-2 active site.

Diversification of Colonization Factors in a Multidrug-Resistant Escherichia coli Lineage Evolving under Negative Frequency-Dependent Selection

Escherichia coli is a major cause of bloodstream and urinary tract infections globally. The wide dissemination of multidrug-resistant (MDR) strains of extraintestinal pathogenic E. coli (ExPEC) poses a rapidly increasing public health burden due to narrowed treatment options and increased risk of failure to clear an infection. Here, we present a detailed population genomic analysis of the ExPEC ST131 clone, in which we seek explanations for its success as an emerging pathogenic strain beyond the acquisition of antimicrobial resistance (AMR) genes. We show evidence for evolution toward separate ecological niches for the main clades of ST131 and differential evolution of anaerobic metabolism, key colonization, and virulence factors. We further demonstrate that negative frequency-dependent selection acting across accessory loci is a major mechanism that has shaped the population evolution of this pathogen.IMPORTANCE Infections with multidrug-resistant (MDR) strains of Escherichia coli are a significant global public health concern. To combat these pathogens, we need a deeper understanding of how they evolved from their background populations. By understanding the processes that underpin their emergence, we can design new strategies to limit evolution of new clones and combat existing clones. By combining population genomics with modelling approaches, we show that dominant MDR clones of E. coli are under the influence of negative frequency-dependent selection, preventing them from rising to fixation in a population. Furthermore, we show that this selection acts on genes involved in anaerobic metabolism, suggesting that this key trait, and the ability to colonize human intestinal tracts, is a key step in the evolution of MDR clones of E. coli.

The first human report of mobile colistin resistance gene, mcr-1, in Finland

Colistin resistance mediated by mobile mcr-1 gene has raised concern during the last years. After steep increase in mcr-1 reports, other mcr-gene variants (mcr-2 to mcr-5) have been revealed as well. In 2016, a clinical study was conducted on asymptomatic stool carriage of extended spectrum beta-lactamase (ESBL) producing Escherichia coli and Klebsiella pneumoniae among Finnish adults. All suspected ESBL producing bacterial isolates were first tested by phenotypic ESBL-confirmation methods, and then further analyzed with whole genome sequencing to identify the resistance genes. We found one study subject carrying a colistin resistant E. coli with a transferrable mcr-1 gene. This multi-drug resistant isolate, although initially suspected to be an ESBL producer, did not carry any ESBL genes, but was proven to carry several other resistance genes by using whole genome sequencing. Sequence type was ST93. The mcr-1 gene was connected to IncX4 plasmid which suggests that the colistin resistance gene locates in the respective plasmid. Here, we report the finding of a mcr-1 harboring human E. coli isolate from Finland. Clinical antimicrobial resistance (AMR) rates are low in Finland, and mobile colistin resistance has not been reported previously. This highlights the importance of AMR surveillance also in populations with low levels of resistance.

Genomic landscape of extended-spectrum β-lactamase resistance in Escherichia coli from an urban African setting

Objectives

Efforts to treat Escherichia coli infections are increasingly being compromised by the rapid, global spread of antimicrobial resistance (AMR). Whilst AMR in E. coli has been extensively investigated in resource-rich settings, in sub-Saharan Africa molecular patterns of AMR are not well described. In this study, we have begun to explore the population structure and molecular determinants of AMR amongst E. coli isolates from Malawi.

Methods

Ninety-four E. coli isolates from patients admitted to Queen's Hospital, Malawi, were whole-genome sequenced. The isolates were selected on the basis of diversity of phenotypic resistance profiles and clinical source of isolation (blood, CSF and rectal swab). Sequence data were analysed using comparative genomics and phylogenetics.

Results

Our results revealed the presence of five clades, which were strongly associated with E. coli phylogroups A, B1, B2, D and F. We identified 43 multilocus STs, of which ST131 (14.9%) and ST12 (9.6%) were the most common. We identified 25 AMR genes. The most common ESBL gene was bla CTX-M-15 and it was present in all five phylogroups and 11 STs, and most commonly detected in ST391 (4/4 isolates), ST648 (3/3 isolates) and ST131 [3/14 (21.4%) isolates].

Conclusions

This study has revealed a high diversity of lineages associated with AMR, including ESBL and fluoroquinolone resistance, in Malawi. The data highlight the value of longitudinal bacteraemia surveillance coupled with detailed molecular epidemiology in all settings, including low-income settings, in describing the global epidemiology of ESBL resistance.

Whole genome sequencing of ESBL-producing Escherichia coli isolated from patients, farm waste and canals in Thailand

BACKGROUND

Tackling multidrug-resistant Escherichia coli requires evidence from One Health studies that capture numerous potential reservoirs in circumscribed geographic areas.

METHODS

We conducted a survey of extended β-lactamase (ESBL)-producing E. coli isolated from patients, canals and livestock wastewater in eastern Thailand between 2014 and 2015, and analyzed isolates using whole genome sequencing.

RESULTS

The bacterial collection of 149 isolates consisted of 84 isolates from a single hospital and 65 from the hospital sewer, canals and farm wastewater within a 20 km radius. E. coli ST131 predominated the clinical collection (28.6%), but was uncommon in the environment. Genome-based comparison of E. coli from infected patients and their immediate environment indicated low genetic similarity overall between the two, although three clinical-environmental isolate pairs differed by ≤ 5 single nucleotide polymorphisms. Thai E. coli isolates were dispersed throughout a phylogenetic tree containing a global E. coli collection. All Thai ESBL-positive E. coli isolates were multidrug resistant, including high rates of resistance to tobramycin (77.2%), gentamicin (77.2%), ciprofloxacin (67.8%) and trimethoprim (68.5%). ESBL was encoded by six different CTX-M elements and SHV-12. Three isolates from clinical samples (n = 2) or a hospital sewer (n = 1) were resistant to the carbapenem drugs (encoded by NDM-1, NDM-5 or GES-5), and three isolates (clinical (n = 1) and canal water (n = 2)) were resistant to colistin (encoded by mcr-1); no isolates were resistant to both carbapenems and colistin.

CONCLUSIONS

Tackling ESBL-producing E. coli in this setting will be challenging based on widespread distribution, but the low prevalence of resistance to carbapenems and colistin suggests that efforts are now required to prevent these from becoming ubiquitous.

Longitudinal genomic surveillance of multidrug-resistant Escherichia coli carriage in a long-term care facility in the United Kingdom

BACKGROUND

Residents of long-term care facilities (LTCF) may have high carriage rates of multidrug-resistant pathogens, but are not currently included in surveillance programmes for antimicrobial resistance or healthcare-associated infections. Here, we describe the value derived from a longitudinal epidemiological and genomic surveillance study of drug-resistant Escherichia coli in a LTCF in the United Kingdom (UK).

METHODS

Forty-five of 90 (50%) residents were recruited and followed for six months in 2014. Participants were screened weekly for carriage of extended-spectrum beta-lactamase (ESBL) producing E. coli. Participants positive for ESBL E. coli were also screened for ESBL-negative E. coli. Phenotypic antibiotic susceptibility of E. coli was determined using the Vitek2 instrument and isolates were sequenced on an Illumina HiSeq2000 instrument. Information was collected on episodes of clinical infection and antibiotic consumption.

RESULTS

Seventeen of 45 participants (38%) carried ESBL E. coli. Twenty-three of the 45 participants (51%) had 63 documented episodes of clinical infection treated with antibiotics. Treatment with antibiotics was associated with higher risk of carrying ESBL E. coli. ESBL E. coli was mainly sequence type (ST)131 (16/17, 94%). Non-ESBL E. coli from these 17 cases was more genetically diverse, but ST131 was found in eight (47%) cases. Whole-genome analysis of 297 ST131 E. coli from the 17 cases demonstrated highly related strains from six participants, indicating acquisition from a common source or person-to-person transmission. Five participants carried highly related strains of both ESBL-positive and ESBL-negative ST131. Genome-based comparison of ST131 isolates from the LTCF study participants with ST131 associated with bloodstream infection at a nearby acute hospital and in hospitals across England revealed sharing of highly related lineages between the LTCF and a local hospital.

CONCLUSIONS

This study demonstrates the power of genomic surveillance to detect multidrug-resistant pathogens and confirm their connectivity within a healthcare network.

Systematic longitudinal survey of invasive Escherichia coli in England demonstrates a stable population structure only transiently disturbed by the emergence of ST131

Escherichia coli associated with urinary tract infections and bacteremia has been intensively investigated, including recent work focusing on the virulent, globally disseminated, multidrug-resistant lineage ST131. To contextualize ST131 within the broader E. coli population associated with disease, we used genomics to analyze a systematic 11-yr hospital-based survey of E. coli associated with bacteremia using isolates collected from across England by the British Society for Antimicrobial Chemotherapy and from the Cambridge University Hospitals NHS Foundation Trust. Population dynamics analysis of the most successful lineages identified the emergence of ST131 and ST69 and their establishment as two of the five most common lineages along with ST73, ST95, and ST12. The most frequently identified lineage was ST73. Compared to ST131, ST73 was susceptible to most antibiotics, indicating that multidrug resistance was not the dominant reason for prevalence of E. coli lineages in this population. Temporal phylogenetic analysis of the emergence of ST69 and ST131 identified differences in the dynamics of emergence and showed that expansion of ST131 in this population was not driven by sequential emergence of increasingly resistant subclades. We showed that over time, the E. coli population was only transiently disturbed by the introduction of new lineages before a new equilibrium was rapidly achieved. Together, these findings suggest that the frequency of E. coli lineages in invasive disease is driven by negative frequency-dependent selection occurring outside of the hospital, most probably in the commensal niche, and that drug resistance is not a primary determinant of success in this niche.

Whole-genome sequencing reveals the effect of vaccination on the evolution of Bordetella pertussis

Herd immunity can potentially induce a change of circulating viruses. However, it remains largely unknown that how bacterial pathogens adapt to vaccination. In this study, Bordetella pertussis, the causative agent of whooping cough, was selected as an example to explore possible effect of vaccination on the bacterial pathogen. We sequenced and analysed the complete genomes of 40 B. pertussis strains from Finland and China, as well as 11 previously sequenced strains from the Netherlands, where different vaccination strategies have been used over the past 50 years. The results showed that the molecular clock moved at different rates in these countries and in distinct periods, which suggested that evolution of the B. pertussis population was closely associated with the country vaccination coverage. Comparative whole-genome analyses indicated that evolution in this human-restricted pathogen was mainly characterised by ongoing genetic shift and gene loss. Furthermore, 116 SNPs were specifically detected in currently circulating ptxP3-containing strains. The finding might explain the successful emergence of this lineage and its spread worldwide. Collectively, our results suggest that the immune pressure of vaccination is one major driving force for the evolution of B. pertussis, which facilitates further exploration of the pathogenicity of B. pertussis.

Proviral integration site for Moloney murine leukemia virus (PIM) kinases promote human T helper 1 cell differentiation

The differentiation of human primary T helper 1 (Th1) cells from naïve precursor cells is regulated by a complex, interrelated signaling network. The identification of factors regulating the early steps of Th1 cell polarization can provide important insight in the development of therapeutics for many inflammatory and autoimmune diseases. The serine/threonine-specific proviral integration site for Moloney murine leukemia virus (PIM) kinases PIM1 and PIM2 have been implicated in the cytokine-dependent proliferation and survival of lymphocytes. We have established that the third member of this family, PIM3, is also expressed in human primary Th cells and identified a new function for the entire PIM kinase family in T lymphocytes. Although PIM kinases are expressed more in Th1 than Th2 cells, we demonstrate here that these kinases positively influence Th1 cell differentiation. Our RNA interference results from human primary Th cells also suggest that PIM kinases promote the production of IFNγ, the hallmark cytokine produced by Th1 cells. Consistent with this, they also seem to be important for the up-regulation of the critical Th1-driving factor, T box expressed in T cells (T-BET), and the IL-12/STAT4 signaling pathway during the early Th1 differentiation process. In summary, we have identified PIM kinases as new regulators of human primary Th1 cell differentiation, thus providing new insights into the mechanisms controlling the selective development of human Th cell subsets.

SNP-based typing: a useful tool to study Bordetella pertussis populations

To monitor changes in Bordetella pertussis populations, mainly two typing methods are used; Pulsed-Field Gel Electrophoresis (PFGE) and Multiple-Locus Variable-Number Tandem Repeat Analysis (MLVA). In this study, a single nucleotide polymorphism (SNP) typing method, based on 87 SNPs, was developed and compared with PFGE and MLVA. The discriminatory indices of SNP typing, PFGE and MLVA were found to be 0.85, 0.95 and 0.83, respectively. Phylogenetic analysis, using SNP typing as Gold Standard, revealed false homoplasies in the PFGE and MLVA trees. Further, in contrast to the SNP-based tree, the PFGE- and MLVA-based trees did not reveal a positive correlation between root-to-tip distance and the isolation year of strains. Thus PFGE and MLVA do not allow an estimation of the relative age of the selected strains. In conclusion, SNP typing was found to be phylogenetically more informative than PFGE and more discriminative than MLVA. Further, in contrast to PFGE, it is readily standardized allowing interlaboratory comparisons. We applied SNP typing to study strains with a novel allele for the pertussis toxin promoter, ptxP3, which have a worldwide distribution and which have replaced the resident ptxP1 strains in the last 20 years. Previously, we showed that ptxP3 strains showed increased pertussis toxin expression and that their emergence was associated with increased notification in the Netherlands. SNP typing showed that the ptxP3 strains isolated in the Americas, Asia, Australia and Europe formed a monophyletic branch which recently diverged from ptxP1 strains. Two predominant ptxP3 SNP types were identified which spread worldwide. The widespread use of SNP typing will enhance our understanding of the evolution and global epidemiology of B. pertussis.

Effect of vaccination on Bordetella pertussis strains, China

Strains in China may differ from those in countries that have long histories of high vaccination coverage.

Whole-cell pertussis vaccine was introduced in China in the early 1960s. We used standard typing methods to compare 96 Bordetella pertussis isolates collected before and after introduction of vaccination, during 1953–2005. The following vaccine-type alleles of the pertussis toxin (ptx) gene were characteristic for all prevaccination strains: ptxA2, ptxA3, and ptxA4. The shift to ptxA1 occurred since 1963. All isolates collected since 1983 contained ptxA1. Pertactin (prn) allele 1, prn1, was predominant, although prn2 and prn3 have been detected since 2000. Serotypes fimbriae (Fim) 2 and Fim2,3 were found in all isolates collected before 1986. During 1997–2005, Fim3 became prevalent. Although changes in electrophoresis profiles over time were observed, the predominant profiles during 1997–2005 resembled those during the prevaccine era and those found in Europe before the 1990s. B. pertussis strains in China may differ from those in countries that have a long history of high vaccine coverage.

Bordetella pertussis vaccine strains and circulating isolates in Serbia

In Serbia, whole cell pertussis vaccine was introduced in 1957. Current composition of the vaccine has been used since 1985 and contains four autochthonous strains of Bordetella pertussis isolated from 1957 to 1984. To monitor changes in bacterial population, 70 isolates collected from 1953 to 2000 were studied together with the vaccine strains. The methods included serotyping of fimbriae (Fim), genotyping of pertactin (prn) and pertussis toxin S1 subunit (ptxA), and pulsed-field gel electrophoresis analysis. Shift from ptxA2 to ptxA1 has been observed in isolates since the late of 1960s. All isolates from 1980 to 1984 harbored ptxA1. Re-appearance of the ptxA2 allele followed an addition of the two strains harboring ptxA1 in the vaccine in 1985. The allele prn1 was predominant among the Serbian isolates, though prn3 and prn11 have been detected since 1981 and 1984. The allele prn2 was found only in two strains isolated in 2000. Serotype Fim2.3 disappeared before 1980 and serotype Fim2 became predominant since then. The Serbian vaccine strains showed differences in ptxA and prn. The results of this present study indicate that the B. pertussis population in Serbia is different from other vaccinated populations and that this difference may be related to the vaccine used.

Comparative genomics of Bordetella pertussis reveals progressive gene loss in Finnish strains

Background

Bordetella pertussis is a Gram-negative bacterium that infects the human respiratory tract and causes pertussis or whooping cough. The disease has resurged in many countries including Finland where the whole-cell pertussis vaccine has been used for more than 50 years. Antigenic divergence has been observed between vaccine strains and clinical isolates in Finland. To better understand genome evolution in B. pertussis circulating in the immunized population, we developed an oligonucleotide-based microarray for comparative genomic analysis of Finnish strains isolated during the period of 50 years.

Methodology/Principal Findings

The microarray consisted of 3,582 oligonucleotides (70-mer) and covered 94% of 3,816 ORFs of Tohama I, the strain of which the genome has been sequenced [1]. Twenty isolates from 1953 to 2004 were studied together with two Finnish vaccine strains and two international reference strains. The isolates were selected according to their characteristics, e.g. the year and place of isolation and pulsed-field gel electrophoresis profiles. Genomic DNA of the tested strains, along with reference DNA of Tohama I strain, was labelled and hybridized. The absence of genes as established with microarrays, was confirmed by PCR. Compared with the Tohama I strain, Finnish isolates lost 7 (8.6 kb) to 49 (55.3 kb) genes, clustered in one to four distinct loci. The number of lost genes increased with time, and one third of lost genes had functions related to inorganic ion transport and metabolism, or energy production and conversion. All four loci of lost genes were flanked by the insertion sequence element IS481.

Conclusion/Significance

Our results showed that the progressive gene loss occurred in Finnish B. pertussis strains isolated during a period of 50 years and confirmed that B. pertussis is dynamic and is continuously evolving, suggesting that the bacterium may use gene loss as one strategy to adapt to highly immunized populations.

Enrichment of nucleofected primary human CD4+ T cells: a novel and efficient method for studying gene function and role in human primary T helper cell differentiation

Identification of key factors mediating the differentiation of naïve CD4(+) T helper cells into Th1 and Th2 subsets is important for understanding the molecular mechanisms of the development of autoimmune diseases as well as asthma and allergy. Functional importance of a given gene in the initiation of human T helper cell differentiation has been hard to study due to the difficulty in transfecting primary resting human T lymphocytes. In this study we have successfully transfected human primary CD4(+) T helper cells using Amaxa's Nucleofection technology. To overcome the background caused by untransfected cells, we have developed a system for enriching nucleofected unstimulated human primary T helper cells that express the gene of interest. This is achieved by introducing a plasmid construct containing a bicistronic unit coding for a truncated mouse MHC class l H-2K(k) cell surface marker followed by selection of H-2K(k) positive cells using antibody coated beads. We demonstrate that the nucleofected and enriched H-2K(k) positive T helper cells differentiate into Th1 and Th2 cells as well as the non-transfected control cells. We also show that by using this novel method, introduction of an shRNA targeting Stat6, a key molecule driving the Th2 cell development, results in impaired Th2 cell differentiation, as expected. The method described here, enables fast and feasible preparation of highly pure transfected primary CD4(+) T cell cultures ideal for studying the influence of overexpression or knockdown of a given gene on T helper cell differentiation and other primary human T cell functions.