Microreact: visualizing and sharing data for genomic epidemiology and phylogeography

Comparative genomics and pangenomics of vancomycin-resistant and susceptible Enterococcus faecium from Irish hospitals

Introduction. Enterococcus faecium has emerged as an important nosocomial pathogen, which is increasingly difficult to treat due to the genetic acquisition of vancomycin resistance. Ireland has a recalcitrant vancomycin-resistant bloodstream infection rate compared to other developed countries.

Hypothesis/Gap statement. Vancomycin resistance rates persist amongst E. faecium isolates from Irish hospitals. The evolutionary genomics governing these trends have not been fully elucidated.

Methodology. A set of 28 vancomycin-resistant isolates was sequenced to construct a dataset alongside 61 other publicly available Irish genomes. This dataset was extensively analysed using in silico methodologies (comparative genomics, pangenomics, phylogenetics, genotypics and comparative functional analyses) to uncover distinct evolutionary, coevolutionary and clinically relevant population trends.

Results. These results suggest that a stable (in terms of genome size, GC% and number of genes), yet genetically diverse population (in terms of gene content) of E. faecium persists in Ireland with acquired resistance arising via plasmid acquisition (vanA) or, to a lesser extent, chromosomal recombination (vanB). Population analysis revealed five clusters with one cluster partitioned into four clades which transcend isolation dates. Pangenomic and recombination analyses revealed an open (whole genome and chromosomal specific) pangenome illustrating a rampant evolutionary pattern. Comparative resistomics and virulomics uncovered distinct chromosomal and mobilomal propensity for multidrug resistance, widespread chromosomal point-mutation-mediated resistance and chromosomally harboured arsenals of virulence factors. Interestingly, a potential difference in biofilm formation strategies was highlighted by coevolutionary analysis, suggesting differential biofilm genotypes between vanA and vanB isolates.

Conclusions. These results highlight the evolutionary history of Irish E. faecium isolates and may provide insight into underlying infection dynamics in a clinical setting. Due to the apparent ease of vancomycin resistance acquisition over time, susceptible E. faecium should be concurrently reduced in Irish hospitals to mitigate potential resistant infections.

Infection and Transmission of SARS-CoV-2 B.1.617.2 Lineage (Delta Variant) among Fully Vaccinated Individuals

The emergence of the SARS-CoV-2 B.1.617.2 lineage (Delta variant) in 2021 was associated with increased case numbers and test positivity rates, including a large number of infections in fully vaccinated individuals. Here, we describe the findings of an investigation conducted in Tompkins County, New York, to evaluate factors underlying a significant uptick in the number of coronavirus disease 2019 (COVID-19) cases observed in the months of July and August 2021. We performed genomic surveillance and genotyping as well as virological assessments to determine infectivity of the virus in a select number of clinical diagnostic samples. Genomic sequence analyses revealed complete replacement of the B.1.1.7 lineage (Alpha variant) with the B.1.617.2 lineage (Delta variant) between July 1 and August 4 2021. We observed a strong association between viral RNA loads detected by real-time reverse transcriptase PCR and infectious virus detected in respiratory secretions by virus titration. A marked increase in positive cases among fully vaccinated individuals was observed. The sequence divergence between two index Delta variant cases in April and May, and the cases after July 1st, revealed independent Delta variant introductions in Tompkins County. Contact tracing information enabled the detection of clusters of connected cases within closely related phylogenetic clusters. We also found evidence of transmission between vaccinated individuals and between vaccinated and unvaccinated individuals. This was confirmed by detection and isolation of infectious virus from a group of individuals within epidemiologically connected transmission clusters, confirming shedding of high viral loads and transmission of the virus by fully vaccinated individuals. IMPORTANCE The SARS-CoV-2 lineage B.1.617.2 (Delta variant) emerged in Asia and rapidly spread to other countries, becoming the dominant circulating lineage. Worldwide infections with B.1.617.2 peaked at a time in which vaccination rates were increasing. In this study, we present data characterizing the emergence of SARS-CoV-2 lineage B.1.617.2 (Delta variant) in Tompkins County, New York, which has one of the highest vaccination rates in the state. We present evidence demonstrating infection, replication, and transmission of SARS-CoV-2 lineage B.1.617.2 (Delta variant) between fully vaccinated individuals. Importantly, infectious virus loads were determined in a subset of samples and demonstrated shedding of high viral titers in respiratory secretions of vaccinated individuals.

Emergence and Molecular Epidemiology of Campylobacter jejuni ST-2993 Associated with a Large Outbreak of Guillain-Barré Syndrome in Peru

Campylobacter jejuni infection is considered the most frequent factor associated with Guillain-Barré syndrome (GBS). In 2019, a large outbreak of GBS was detected in Peru, being associated with C. jejuni detected in stool samples from these patients. The aim of this study was to determine the molecular epidemiology of C. jejuni strains (ST-2993) associated with a large GBS outbreak in Peru. In this study, 26 C. jejuni strains belonging to the ST-2293, obtained from 2019 to 2020, were sequenced using Illumina technology. Five low-quality sequences were removed using bioinformatics, and 21 genomes (17 clinical strains and 4 chicken strains) were considered in the phylogenetic analysis and comparative genomics. Phylogenetic reconstruction, including genomes from international databases, showed a connection between Peruvian and Chinese GBS strains, both of them having lipooligosaccharides (LOS) locus genes related to molecular mimicry with gangliosides in peripheral nerves. Also, ST-2993 was detected in Amazon strains recovered many years before the 2019 outbreak, but with no epidemiological connection with GBS. Besides, a close relationship between human and chicken C. jejuni strains indicated chicken as one of the probable reservoirs. Finally, comparative genomics revealed differences between Chinese and Peruvian strains, including the presence of a prophage inserted into the genome. In conclusion, C. jejuni ST-2993 strains recovered from the GBS outbreak are closely related to Peruvian Amazon strains. Moreover, ST-2993 has been circulated in Peru since 2003 in the Peruvian Amazonia, showing the necessity to reinforce the epidemiological surveillance of C. jejuni to improve the prevention and control of future GBS outbreaks. IMPORTANCE This article describes the molecular epidemiology of C. jejuni strains (ST-2993) associated with a large Guillain-Barré Syndrome (GBS) outbreak in Peru, sequencing several strains recovered from GBS patients and chickens from 2019 to 2020. Phylogenetic analysis showed a connection between Peruvian and Chinese GBS strains, both of them having lipooligosaccharides (LOS) locus genes related to molecular mimicry with gangliosides in peripheral nerves. Also, ST-2993 strains were detected in isolates recovered many years before the 2019 outbreak, but with no epidemiological connection with GBS. Besides, a close relationship between human and chicken strains indicated those animals as a probable reservoir. This information will help to understand the real situation of GBS in Peru and its causal agent, C. jejuni ST-2993, showing the necessity to increase epidemiological tracking of these kinds of pathogens to detect them and avoid GBS outbreaks in the future.

Optimization and Validation of Candida auris Short Tandem Repeat Analysis

Candida auris is an easily transmissible yeast with resistance to different antifungal compounds. Outbreaks of C. auris are mostly observed in intensive care units. To take adequate measures during an outbreak, it is essential to understand the transmission route, which requires isolate genotyping. In 2019, a short tandem repeat (STR) genotyping analysis was developed for C. auris. To determine the discriminatory power of this method, we performed STR analysis of 171 isolates with known whole-genome sequencing (WGS) data using Illumina reads, and we compared their resolutions. We found that STR analysis separated the 171 isolates into four clades (clades I to IV), as was also seen with WGS analysis. Then, to improve the separation of isolates in clade IV, the STR assay was optimized by the addition of 2 STR markers. With this improved STR assay, a total of 32 different genotypes were identified, while all isolates with differences of >50 single-nucleotide polymorphisms (SNPs) were separated by at least 1 STR marker. Altogether, we optimized and validated the C. auris STR panel for clades I to IV and established its discriminatory power, compared to WGS SNP analysis using Illumina reads. IMPORTANCE The emerging fungal pathogen Candida auris poses a threat to public health, mainly causing outbreaks in intensive care units. Genotyping is essential for investigating potential outbreaks and preventing further spread. Previously, we developed a STR genotyping scheme for rapid and high-resolution genotyping, and WGS SNP outcomes for some isolates were compared to STR data. Here, we compared WGS SNP and STR outcomes for a larger sample cohort. Also, we optimized the resolution of this typing scheme with the addition of 2 STR markers. Altogether, we validated and optimized this rapid, reliable, and high-resolution typing scheme for C. auris.

Genomic Diversity of NDM-Producing Klebsiella Species from Brazil, 2013–2022

Background: Since its first report in the country in 2013, NDM-producing Enterobacterales have been identified in all the Brazilian administrative regions. In this study, we characterized by antimicrobial susceptibility testing and by molecular typing a large collection of NDM-producing Klebsiella isolates from different hospitals in Brazil, mainly from the state of Sao Paulo, over the last decade. Methods: Bacterial isolates positive for bla_NDM-genes were identified by MALDI-TOF MS and submitted to antimicrobial susceptibility testing by disk diffusion or broth microdilution (for polymyxin B). All isolates were submitted to pulsed-field gel electrophoresis, and isolates belonging to different clusters were submitted to whole genome sequencing by Illumina technology and downstream analysis. Mating out assays were performed by conjugation, plasmid sizes were determined by S1-PFGE, and plasmid content was investigated by hybrid assembly after MinIon long reads sequencing. Results: A total of 135 NDM-producing Klebsiella were identified, distributed into 107 different pulsotypes; polymyxin B was the only antimicrobial with high activity against 88.9% of the isolates. Fifty-four isolates presenting diversified pulsotypes were distributed in the species K. pneumoniae (70%), K. quasipneumoniae (20%), K. variicola (6%), K. michiganensis (a K. oxytoca Complex species, 2%), and K. aerogenes (2%); bla_NDM-1 was the most frequent allele (43/54, 80%). There was a predominance of Clonal Group 258 (ST11 and ST340) encompassing 35% of K. pneumoniae isolates, but another thirty-one different sequence types (ST) were identified, including three described in this study (ST6244 and ST6245 for K. pneumoniae, and ST418 for K. michiganensis). The bla_NDM-1 and bla_NDM-7 were found to be located into IncF and IncX3 type transferable plasmids, respectively. Conclusions: Both clonal (mainly driven by CG258) and non-clonal expansion of NDM-producing Klebsiella have been occurring in Brazil in different species and clones, associated with different plasmids, since 2013.

Population genomics of emerging Elizabethkingia anophelis pathogens reveals potential outbreak and rapid global dissemination

Elizabethkingia anophelis is an emerging species and have increasingly been reported to cause life-threatening infections and even outbreaks in humans. Nevertheless, there is little data regarding the E. anophelis geographical distribution, phylogenetic structure, and transmission across the globe, especially in Asia. We utilize whole genome sequencing (WGS) data to define a global population framework, phylogenetic structure, geographical distribution, and transmission evaluation of E. anophelis pathogens. The geographical distribution diagram revealed the emerging pathogenic bacteria already distributed in various countries worldwide, especially in the USA and China. Strikingly, phylogenetic analysis showed a part of our China original E. anophelis shared the same ancestor with the USA outbreak strain, which implies the possibility of localized outbreaks and global spread. These closer related strains also contained ICEEaI, which might insert into a disrupted DNA repair mutY gene and made the strain more liable to mutation and outbreak infection. BEAST analysis showed that the most recent common ancestor for ICEEaI E. anophelis was dated twelve years ago, and China might be the most likely recent source of this bacteria. Our study sheds light on the potential possibility of E. anophelis causing the large-scale outbreak and rapid global dissemination. Continued genomic surveillance of the dynamics of E. anophelis populations will generate further knowledge for optimizing future prevent global outbreak infections.

An Unusual Resurgence of Human Metapneumovirus in Western Australia Following the Reduction of Non-Pharmaceutical Interventions to Prevent SARS-CoV-2 Transmission

Non-pharmaceutical interventions (NPIs) to reduce SARS-CoV-2 transmission disrupted respiratory virus seasonality. We examined the unusual return of human metapneumovirus (hMPV) in Western Australia following a period of absence in 2020. We analysed hMPV laboratory testing data from 1 January 2017 to 31 December 2021. Whole-genome sequencing of selected hMPV-positive samples was performed using a tiled-amplicon approach. Following an absence in spring 2020, an unusual hMPV surge was observed during the wet summer season in the tropical Northern region in late 2020. Following a six-month delay, an intense winter season occurred in the subtropical/temperate Southern and Metropolitan regions. Compared to 2017-2019, hMPV incidence in 2021 increased by 3-fold, with a greater than 4-fold increase in children aged 1-4 years. There was a collapse in hMPV diversity in 2020, with the emergence of a single subtype. NPIs contributed to an absent 2020 season and a clonal hMPV resurgence. The summer surge and delayed winter season suggest that prevailing temperature and humidity are keys determinant of hMPV transmission. The increased incidence in 2021 was linked to an expanded cohort of hMPV-naïve 1-4-year-old children and waning population immunity. Further intense and unusual respiratory virus seasons are expected as COVID-19 associated NPIs are removed.

SPREAD 4: online visualisation of pathogen phylogeographic reconstructions

Phylogeographic analyses aim to extract information about pathogen spread from genomic data, and visualising spatio-temporal reconstructions is a key aspect of this process. Here we present SPREAD 4, a feature-rich web-based application that visualises estimates of pathogen dispersal resulting from Bayesian phylogeographic inference using BEAST on a geographic map, offering zoom-and-filter functionality and smooth animation over time. SPREAD 4 takes as input phylogenies with both discrete and continuous location annotation and offers customised visualisation as well as generation of publication-ready figures. SPREAD 4 now features account-based storage and easy sharing of visualisations by means of unique web addresses. SPREAD 4 is intuitive to use and is available online at https://spreadviz.org, with an accompanying web page containing answers to frequently asked questions at https://beast.community/spread4.

Detection of Alpha- and Betacoronaviruses in Frugivorous and Insectivorous Bats in Nigeria

The rise of bat-associated zoonotic viruses necessitates a close monitoring of their natural hosts. Since the detection of severe acute respiratory syndrome coronavirus (SARS-CoV), it is evident that bats are vital reservoirs of coronaviruses (CoVs). In this study, we investigated the presence of CoVs in multiple bat species in Nigeria to identify viruses in bats at high-risk human contact interfaces. Four hundred and nine bats comprising four bat species close to human habitats were individually sampled from five states in Nigeria between 2019 and 2021. Coronavirus detection was done using broadly reactive consensus PCR primers targeting the RNA-dependent RNA polymerase (RdRp) gene of CoVs. Coronavirus RNA was detected in 39 samples (9.5%, CI 95%: [7.0, 12.8]), of which 29 were successfully sequenced. The identified CoVs in Nigerian bats were from the unclassified African alphacoronavirus lineage and betacoronavirus lineage D (Nobecovirus), with one sample from Hipposideros ruber coinfected with alphacoronavirus and betacoronavirus. Different bat species roosting in similar or other places had CoVs from the same genetic lineage. The phylogenetic and evolutionary dynamics data indicated a high CoV diversity in Nigeria, while host switching may have contributed to CoV evolution. Robust sentinel surveillance is recommended to enhance our knowledge of emerging and re-emerging coronaviruses.

Molecular epidemiology and antibiotic resistance profiles of invasive Haemophilus influenzae from Norway 2017-2021

Invasive Haemophilus influenzae (Hi) disease has decreased in countries that included Hi type b (Hib) vaccination in their childhood immunization programs in the 1990s. Non-typeable (NT) and non-b strains are now the leading causes of invasive Hi disease in Europe, with most cases reported in young children and the elderly. Concerningly, no vaccines toward such strains are available and beta-lactam resistance is increasing. We describe the epidemiology of invasive Hi disease reported to the Norwegian Surveillance System for Communicable Diseases (MSIS) (2017-2021, n = 407). Whole-genome sequencing (WGS) was performed on 245 isolates. We investigated the molecular epidemiology (core genome phylogeny) and the presence of antibiotic resistance markers (including chromosomal mutations associated with beta-lactam or quinolone resistance). For isolates characterized with both WGS and phenotypic antibiotic susceptibility testing (AST) (n = 113) we assessed correlation between resistance markers and susceptibility categorization by calculation of sensitivity, specificity, and predictive values. Incidence rates of invasive Hi disease in Norway ranged from 0.7 to 2.3 per 100,000 inhabitants/year (mean 1.5 per 100,000) and declined during the COVID-19 pandemic. The bacterial population consisted of two major phylogenetic groups with subclustering by serotype and multi-locus sequence type (ST). NTHi accounted for 71.8% (176). The distribution of STs was in line with previous European reports. We identified 13 clusters, including four encapsulated and three previously described international NTHi clones with bla TEM-1 (ST103) or altered PBP3 (rPBP3) (ST14/IIA and ST367/IIA). Resistance markers were detected in 25.3% (62/245) of the isolates, with bla TEM-1 (31, 50.0%) and rPBP3 (28, 45.2%) being the most frequent. All isolates categorized as resistant to aminopenicillins, tetracycline or chloramphenicol possessed relevant resistance markers, and the absence of relevant substitutions in PBP3 and GyrA/ParC predicted susceptibility to cefotaxime, ceftriaxone, meropenem and quinolones. Among the 132 WGS-only isolates, one isolate had PBP3 substitutions associated with resistance to third-generation cephalosporins, and one isolate had GyrA/ParC alterations associated with quinolone resistance. The detection of international virulent and resistant NTHi clones underlines the need for a global molecular surveillance system. WGS is a useful supplement to AST and should be performed on all invasive isolates.

Molecular diversity of Paenibacillus larvae strains isolated from Lithuanian apiaries

Paenibacillus larvae bacterium is known to be the causative agent of American foulbrood (AFB), a widespread, highly contagious and fatal disease in honey bees (Apis mellifera). There are four genotypes of Paenibacillus larvae that are named after their enterobacterial repetitive consensus (ERIC), and a fifth ERIC genotype has recently been found. In this study, a total of 108 independent P. larvae isolates from different geographical regions in Lithuania collected between 2011 and 2021 were investigated by molecular methods. The aims of this study were to detect which enterobacterial repetitive intergenic consensus (ERIC) genotype is the most common in Lithuania apiaries, identify and differentiate subtypes of the defined genotype by using multiple-locus variable number of tandem-repeat analysis (MLVA), and review how bacterial molecular diversity has changed over time in different parts of Lithuania. The obtained molecular analysis results showed that 100% of P. larvae bacterial isolates from Lithuania belong to the ERIC I genotype and can be differentiated to nine different subtypes by using the MLVA and capillary electrophoresis methods.

Ancestral lineage of SARS-CoV-2 is more stable in human biological fluids than Alpha, Beta and Omicron variants of concern

SARS-CoV-2 is a zoonotic virus which was first identified in 2019, and has quickly spread worldwide. The virus is primarily transmitted through respiratory droplets from infected persons; however, the virus-laden excretions can contaminate surfaces which can serve as a potential source of infection. Since the beginning of the pandemic, SARS-CoV-2 has continued to evolve and accumulate mutations throughout its genome leading to the emergence of variants of concern (VOCs) which exhibit increased fitness, transmissibility, and/or virulence. However, the stability of SARS-CoV-2 VOCs in biological fluids has not been thoroughly investigated so far. The aim of this study was to determine and compare the stability of different SARS-CoV-2 strains in human biological fluids. Here, we demonstrate that the ancestral strain of Wuhan-like lineage A was more stable than the Alpha VOC B.1.1.7, and the Beta VOC B.1.351 strains in human liquid nasal mucus and sputum. In contrast, there was no difference in stability among the three strains in dried biological fluids. Furthermore, we also show that the Omicron VOC B.1.1.529 strain was less stable than the ancestral Wuhan-like strain in liquid nasal mucus. These studies provide insight into the effect of the molecular evolution of SARS-CoV-2 on environmental virus stability, which is important information for the development of countermeasures against SARS-CoV-2.

Importance

Genetic evolution of SARS-CoV-2 leads to the continuous emergence of novel variants, posing a significant concern to global public health. Five of these variants have been classified so far into variants of concern (VOCs); Alpha, Beta, Gamma, Delta, and Omicron. Previous studies investigated the stability of SARS-CoV-2 under various conditions, but there is a gap of knowledge on the survival of SARS-CoV-2 VOCs in human biological fluids which are clinically relevant. Here, we present evidence that Alpha, Beta, and Omicron VOCs were less stable than the ancestral Wuhan-like strain in human biological fluids. Our findings highlight the potential risk of contaminated human biological fluids in SARS-CoV-2 transmission and contribute to the development of countermeasures against SARS-CoV-2.

Molecular characterisation of Acinetobacter baumannii isolates from bloodstream infections in a tertiary-level hospital in South Africa

Acinetobacter baumannii is an opportunistic pathogen and causes various infections in patients. This study aimed to describe the clinical, epidemiological and molecular characteristics of A. baumannii isolated from BCs in patients at a tertiary-level hospital in South Africa. Ninety-six isolates from bloodstream infections were collected. Clinical characteristics of patients were recorded from patient files. Organism identification and AST was performed using automated systems. PCR screening for the mcr-1 to mcr-5 genes was done. To infer genetic relatedness, a dendrogram was constructed using MALDI-TOF MS. All colistin-resistant isolates (n = 9) were selected for WGS. The patients were divided into three groups, infants (<1 year; n = 54), paediatrics (1–18 years; n = 6) and adults (≥19 years; n = 36) with a median age of 13 days, 1 and 41 years respectively. Of the 96 A. baumannii bacteraemia cases, 96.9% (93/96) were healthcare-associated. The crude mortality rate at 30 days was 52.2% (48/92). The majority of the isolates were multidrug-resistant (MDR). All isolates were PCR-negative for the mcr-1 to mcr-5 genes. The majority of the isolates belonged to cluster 1 (62/96) according to the MALDI-TOF MS dendrogram. Colistin resistance was confirmed in nine A. baumannii isolates (9.4%). The colistin-resistant isolates belonged to sequence type (ST) 1 (5/6) and ST2 (1/6). The majority of ST1 isolates showed low SNP diversity (≤4 SNPs). All the colistin-resistant isolates were resistant to carbapenems, exhibited an XDR phenotype and harboured the bla_OXA–23 gene. The bla_NDM gene was only detected in ST1 colistin-resistant isolates (n = 5). The lpsB gene was detected in all colistin-resistant isolates as well as various efflux pump genes belonging to the RND, the MFS and the SMR families. The lipooligosaccharide OCL1 was detected in all colistin-resistant ST1 and ST2 isolates and the capsular polysaccharide KL3 and KL17 were detected in ST2 and ST1 respectively. This study demonstrated a 9.4% prevalence of colistin-resistant ST1 and ST2 A. baumannii in BC isolates. The detection of the lpsB gene indicates a potential threat and requires close prospective monitoring.

Phylodynamics of HIV in the Mexico City Metropolitan Region

Evolutionary analyses of viral sequences can provide insights into transmission dynamics, which in turn can optimize prevention interventions. Here, we characterized the dynamics of HIV transmission within the Mexico City metropolitan area. HIV pol sequences from persons recently diagnosed at the largest HIV clinic in Mexico City (between 2016 and 2021) were annotated with demographic/geographic metadata. A multistep phylogenetic approach was applied to identify putative transmission clades. A data set of publicly available sequences was used to assess international introductions. Clades were analyzed with a discrete phylogeographic model to evaluate the timing and intensity of HIV introductions and transmission dynamics among municipalities in the region. A total of 6,802 sequences across 96 municipalities (5,192 from Mexico City and 1,610 from the neighboring State of Mexico) were included (93.6% cisgender men, 5.0% cisgender women, and 1.3% transgender women); 3,971 of these sequences formed 1,206 clusters, involving 78 municipalities, including 89 clusters of ≥10 sequences. Discrete phylogeographic analysis revealed (i) 1,032 viral introductions into the region, over one-half of which were from the United States, and (ii) 354 migration events between municipalities with high support (adjusted Bayes factor of ≥3). The most frequent viral migrations occurred between northern municipalities within Mexico City, i.e., Cuauhtémoc to Iztapalapa (5.2% of events), Iztapalapa to Gustavo A. Madero (5.4%), and Gustavo A. Madero to Cuauhtémoc (6.5%). Our analysis illustrates the complexity of HIV transmission within the Mexico City metropolitan area but also identifies a spatially active transmission area involving a few municipalities in the north of the city, where targeted interventions could have a more pronounced effect on the entire regional epidemic. IMPORTANCE Phylogeographic investigation of the Mexico City HIV epidemic illustrates the complexity of HIV transmission in the region. An active transmission area involving a few municipalities in the north of the city, with transmission links throughout the region, is identified and could be a location where targeted interventions could have a more pronounced effect on the entire regional epidemic, compared with those dispersed in other manners.

Comparative Dynamics of Delta and Omicron SARS-CoV-2 Variants across and between California and Mexico

Evolutionary analysis using viral sequence data can elucidate the epidemiology of transmission. Using publicly available SARS-CoV-2 sequence and epidemiological data, we developed discrete phylogeographic models to interrogate the emergence and dispersal of the Delta and Omicron variants in 2021 between and across California and Mexico. External introductions of Delta and Omicron in the region peaked in early July (2021-07-10 [95% CI: 2021-04-20, 2021-11-01]) and mid-December (2021-12-15 [95% CI: 2021-11-14, 2022-01-09]), respectively, 3 months and 2 weeks after first detection. These repeated introductions coincided with domestic migration events with no evidence of a unique transmission hub. The spread of Omicron was most consistent with gravity centric patterns within Mexico. While cross-border events accounted for only 5.1% [95% CI: 4.3-6] of all Delta migration events, they accounted for 20.6% [95% CI: 12.4-29] of Omicron movements, paralleling the increase in international travel observed in late 2021. Our investigations of the Delta and Omicron epidemics in the California/Mexico region illustrate the complex interplay and the multiplicity of viral and structural factors that need to be considered to limit viral spread, even as vaccination is reducing disease burden. Understanding viral transmission patterns may help intra-governmental responses to viral epidemics.

Phylogenomic Analysis of Salmonella enterica Serovar Indiana ST17, an Emerging Multidrug-Resistant Clone in China

Salmonella enterica serovar Indiana (S. Indiana) is an extremely expanded foodborne pathogen in China in recent years. This study aimed to elucidate the national prevalence and phylogenomic characterization of this pathogen in China. Among 5, 287 serotyped Salmonella isolates collected during 2002 to 2018, 466 S. Indiana isolates were found in 15 provinces, and 407 were identified to be ST17, and the rest were ST2040. Among 407 ST17 isolates, 372 (91.4%) were multidrug resistant, and 366 (89.9%) were resistant to ciprofloxacin, 235 (57.7%) were further resistant to ceftriaxone. Phylogenomic analysis revealed that ST17 isolates were classified into four clades (I, II, III and IV), which appeared in international clonal dissemination. ST17 isolates from China fell into Clade IV with part of isolates from the United Kingdom, the United States, South Korea, and Thailand, suggesting their close genetic relationship. Mutations in quinolone resistance-determining regions (QRDR) of GyrA and ParC, and plasmid-mediated quinolone resistance (PMQR) genes aac(6′)-Ib-cr, oqxAB, and qnrS as well as extended spectrum β-lactamases (ESBL) genes bla_CTX-M, bla_OXA, and bla_TEM in isolates from Clade IV were much higher than those from other three clades. Various bla_CTX-M subtypes (bla_CTX-M-65, bla_CTX-M-55, bla_CTX-M-27, bla_CTX-M-14, and bla_CTX-M-123) with ISEcp1, IS903B, ISVsa5, and IS1R were found in ST17 isolates, especially Tn1721 containing ΔISEcp1-bla_CTX-M-27-IS903B in P1-like bacteriophage plasmids. These findings on the prevalent and genomic characterization for the S. Indiana multidrug-resistant ST17 clone in China, which have not been reported yet, provide valuable insights into the potential risk of this high-resistant clone.

IMPORTANCE Fluoroquinolones and cephalosporins are the primary choices for severe salmonellosis treatment. S. Indiana has become one of the most prevalent serovars in breeding poultry and poultry meats in China in recent years. ST17 was recognized as the leading epidemiological importance in S. Indiana because of its high-level resistance to the most of common antibiotics, including ciprofloxacin and ceftriaxone. However, the prevalence and phylogenomic characterization of ST17 isolates are unclear. Here, we did a retrospective screening on a large scale for S. Indiana in China, and performed its phylogenomic analysis. It was found that ST17 isolates had extensive spread in 15 provinces of China and became a multidrug-resistant clone. The international spread of the ST17 isolates was observed among several countries, especially China, the United Kingdom, and the United States. Our study emphasized the importance of surveillance of a high-resistant S. Indiana ST17 clone to combat its threat to public health.

Mating pair stabilization mediates bacterial conjugation species specificity

Bacterial conjugation mediates contact-dependent transfer of DNA from donor to recipient bacteria, thus facilitating the spread of virulence and resistance plasmids. Here we describe how variants of the plasmid-encoded donor outer membrane (OM) protein TraN cooperate with distinct OM receptors in recipients to mediate mating pair stabilization and efficient DNA transfer. We show that TraN from the plasmid pKpQIL (Klebsiella pneumoniae) interacts with OmpK36, plasmids from R100-1 (Shigella flexneri) and pSLT (Salmonella Typhimurium) interact with OmpW, and the prototypical F plasmid (Escherichia coli) interacts with OmpA. Cryo-EM analysis revealed that TraN_pKpQIL interacts with OmpK36 through the insertion of a β-hairpin in the tip of TraN into a monomer of the OmpK36 porin trimer. Combining bioinformatic analysis with AlphaFold structural predictions, we identified a fourth TraN structural variant that mediates mating pair stabilization by binding OmpF. Accordingly, we devised a classification scheme for TraN homologues on the basis of structural similarity and their associated receptors: TraNα (OmpW), TraNβ (OmpK36), TraNγ (OmpA), TraNδ (OmpF). These TraN-OM receptor pairings have real-world implications as they reflect the distribution of resistance plasmids within clinical Enterobacteriaceae isolates, demonstrating the importance of mating pair stabilization in mediating conjugation species specificity. These findings will allow us to predict the distribution of emerging resistance plasmids in high-risk bacterial pathogens.

Combining conjugation and structural analyses, the authors show that TraN-OMP pairings determine bacterial conjugation species specificity, with implications in resistance plasmid distribution within Enterobacteriaceae.

Estimating the genetic structure of Triatoma dimidiata (Hemiptera: Reduviidae) and the transmission dynamics of Trypanosoma cruzi in Boyacá, eastern Colombia

Chagas disease is considered a public health issue in Colombia, where many regions are endemic. Triatoma dimidiata is an important vector after Rhodnius prolixus, and it is gaining importance in Boyacá, eastern Colombia. Following the recent elimination of R. prolixus in the region, it is pivotal to understand the behavior of T. dimidiata and the transmission dynamics of T. cruzi. We used qPCR and Next Generation Sequencing (NGS) to evaluate T. cruzi infection, parasite load, feeding profiles, and T. cruzi genotyping for T. dimidiata specimens collected in nine municipalities in Boyacá and explored T. dimidiata population genetics. We found that T. dimidiata populations are composed by a single population with similar genetic characteristics that present infection rates up to 70%, high parasite loads up to 1.46 × 10⁹ parasite-equivalents/mL, a feeding behavior that comprises at least 17 domestic, synanthropic and sylvatic species, and a wide diversity of TcI genotypes even within a single specimen. These results imply that T. dimidiata behavior is similar to other successful vectors, having a wide variety of blood sources and contributing to the circulation of different genotypes of the parasite, highlighting its importance for T. cruzi transmission and risk for humans. In the light of the elimination of R. prolixus in Boyacá and the results we found, we suggest that T. dimidiata should become a new target for vector control programs. We hope this study provides enough information to enhance surveillance programs and a future effective interruption of T. cruzi vector transmission in endemic regions.

Chagas disease is a complex zoonotic infection caused by the protozoan Trypanosoma cruzi. This pathology is endemic in the Americas and causes a tremendous burden in terms of public health. The feces of triatomine bugs mainly transmit this parasite. A massive diversity of triatomines can be found in the north of South America, where Rhodnius is considered the most epidemiologically relevant genus. However, government efforts have attempted to control the vector transmission of specific regions. That is the case of Boyaca in eastern Colombia, which has several municipalities certified as free of R. prolixus transmission of the parasite. However, other species such as Triatoma dimidiata can occupy the left niche due to R. prolixus elimination. We explored the infection rate, parasite load, feeding preferences, and T. cruzi diversity in T. dimidiata specimens collected in municipalities with no R. prolixus infestation. Our results highlight the preponderant need for increasing serological surveillance and prevention in those communities due to the risk of a plausible reactivation of T. cruzi vector transmission due to T. dimidiata.

HIV-1 Subtype Shift in the Philippines is Associated with High Transmitted Drug Resistance, High Viral Loads and Fast Immunologic Decline

BACKGROUND

The Philippines has one of the fastest growing HIV epidemics. A subtype shift from B to CRF01_AE may have contributed to the rise in cases. We undertook a genotyping and transmitted drug resistance (TDR) study to determine if the dominant subtype has any advantages in resistance and transmission.

METHODS

Treatment-naive Filipinos living with HIV were recruited from two large government treatment hubs March 2016 to August 2018. HIV-1 viral load, CD4 count, genotyping and TDR testing were performed. Demographic and clinical data was collected and compared across subtypes.

RESULTS

Two hundred ninety-eight Filipino PLHIV were recruited. Median CD4 count was 143 cells/µL and HIV viral load was 2,345,431 copies/mL. Sanger-based sequencing showed 230/298 (77.2%) had subtype CRF01_AE, 41 (13.8%) subtype B, and the rest other subtypes or recombinants. Overall TDR was 11.7%. TDR was associated with lower viral loads and no previous HIV testing. CRF01_AE had a higher likelihood of a viral load >100,000 copies/mL and having a baseline CD4 count <50 cells/mm3.

CONCLUSIONS

TDR in the Philippines is high at 11.7%. CRF01_AE was observed to have a higher baseline viral load and lower CD4 counts compared to other co-circulating subtypes. Further research needs to confirm this observation since it suggests that CRF01_AE may have a survival advantage that led to replacement of subtype B as the dominant subtype. Drug-resistance testing is recommended in the Philippines when initiating NNRTI-based anti-retroviral therapy but may not be necessary for INSTI-based regimens.

Outbreak.info genomic reports: scalable and dynamic surveillance of SARS-CoV-2 variants and mutations

The emergence of SARS-CoV-2 variants of concern has prompted the need for near real-time genomic surveillance to inform public health interventions. In response to this need, the global scientific community, through unprecedented effort, has sequenced and shared over 11 million genomes through GISAID, as of May 2022. This extraordinarily high sampling rate provides a unique opportunity to track the evolution of the virus in near real-time. Here, we present outbreak.info, a platform that currently tracks over 40 million combinations of PANGO lineages and individual mutations, across over 7,000 locations, to provide insights for researchers, public health officials, and the general public. We describe the interpretable and opinionated visualizations in the variant and location focussed reports available in our web application, the pipelines that enable the scalable ingestion of heterogeneous sources of SARS-CoV-2 variant data, and the server infrastructure that enables widespread data dissemination via a high performance API that can be accessed using an R package. We present a case study that illustrates how outbreak.info can be used for genomic surveillance and as a hypothesis generation tool to understand the ongoing pandemic at varying geographic and temporal scales. With an emphasis on scalability, interactivity, interpretability, and reusability, outbreak.info provides a template to enable genomic surveillance at a global and localized scale.

Shotgun metagenomic sequencing of the first case of monkeypox virus in Brazil, 2022

Monkeypox virus (MPXV), a zoonotic virus endemic to the African continent, has been reported in 33 non-endemic countries since May 2022. We report an almost complete genome of the first confirmed case of MPXV in Brazil. Shotgun metagenomic sequencing was completed in 18 hours, from DNA extraction to consensus sequence generation.

An outbreak sustained by ST15 Klebsiella pneumoniae carrying 16S rRNA methyltransferases and blaNDM: evaluation of the global dissemination of these resistance determinants

The spread of extremely-drug resistant Klebsiella pneumoniae has become a major health threat worldwide. This is largely mediated by certain lineages, recognized as high-risk clones dispersed in all the world. The analysis of an outbreak of nine ST15, NDM-1 metallo-β-lactamase producing K. pneumoniae was performed. An IncC plasmid carrying the bla_NDM-1 gene also carried the rare rmtC gene, encoding for a 16S rRNA methyltransferases (16RMTases), conferring resistance to all aminoglycosides. We studied the global spread of NDM variants and their association with the 16RMTases among K. pneumoniae complete genomes available in GenBank, producing a complete overview of the association of 16RMTases and NDM in K. pneumoniae genomics. NDM is more and more often associated with16RMTases and both are spreading in K. pneumoniae, conferring resistance to every beta-lactam and aminoglycoside. Our analysis suggest that aminoglycosides have limited future as second line treatment against NDM-producing K. pneumoniae.

Local accessory gene sharing among Egyptian Campylobacter potentially promotes the spread of antimicrobial resistance

Campylobacter is the most common cause of bacterial gastroenteritis worldwide, and diarrhoeal disease is a major cause of child morbidity, growth faltering and mortality in low- and middle-income countries. Despite evidence of high incidence and differences in disease epidemiology, there is limited genomic data from studies in developing countries. In this study, we aimed to quantify the extent of gene sharing in local and global populations. We characterized the genetic diversity and accessory-genome content of a collection of Campylobacter isolates from the Cairo metropolitan area, Egypt. In total, 112 Campylobacter isolates were collected from broiler carcasses (n=31), milk and dairy products (n=24), and patients suffering from gastroenteritis (n=57). Among the most common sequence types (STs), we identified the globally disseminated host generalist ST-21 clonal complex (CC21) and the poultry specialists CC206, CC464 and CC48. Notably, CC45 and the cattle-specialist CC42 were under-represented, with a total absence of CC61. Core- and accessory-genome sharing was compared among isolates from Egypt and a comparable collection from the UK (Oxford). Lineage-specific accessory-genome sharing was significantly higher among isolates from the same country, particularly CC21, which demonstrated greater local geographical clustering. In contrast, no geographical clustering was noted in either the core or accessory genome of CC828, suggesting a highly admixed population. A greater proportion of Campylobacter coli isolates were multidrug resistant compared to Campylobacter jejuni. Our results suggest that there is more horizontal transfer of accessory genes between strains in Egypt. This has strong implications for controlling the spread of antimicrobial resistance among this important pathogen.

The use of core genome multilocus sequence typing to determine the duration of vancomycin-resistant Enterococcus faecium outbreaks

The prevalence of vancomycin-resistant Enterococcus faecium has increased rapidly, and in Denmark, we are facing an endemic outbreak situation in hospitals. The aim of this study was to use whole-genome sequencing (WGS) and core genome multilocus sequencing typing (cgMLST) to determine the duration of VREfm outbreaks and thereby evaluate the effect of our infection control strategies. We included all VREfm isolates from six hospitals in the Capital Region of Denmark that were sequenced between 2012 and 2020. Ward data were collected from our laboratory information system. A ward outbreak was defined as two patient samples from the same ward within a period of 30 days belonging to the same cgMLST cluster. cgMLST complex types were determined using Ridom SeqSphere v7.2.3, where a maximum of 20 allelic differences between isolates defines a cluster. We included 1690 patient isolates between 2012 and 2020. Our collection consisted of 45 unique clusters and 227 ward outbreaks. The median duration of outbreaks was 20 days. We reported a median outbreak duration of VREfm outbreaks based on WGS data to be 20 days, and thus concluded that our infection control precautions are adequate.

Whole genome sequence analysis of Salmonella Typhi provides evidence of phylogenetic linkage between cases of typhoid fever in Santiago, Chile in the 1980s and 2010-2016

Typhoid fever epidemiology was investigated rigorously in Santiago, Chile during the 1980s, when Salmonella enterica serovar Typhi (S. Typhi) caused seasonal, hyperendemic disease. Targeted interventions reduced the annual typhoid incidence rates from 128–220 cases/10⁵ population occurring between 1977–1984 to <8 cases/10⁵ from 1992 onwards. As such, Santiago represents a contemporary example of the epidemiologic transition of an industrialized city from amplified hyperendemic typhoid fever to a period when typhoid is no longer endemic. We used whole genome sequencing (WGS) and phylogenetic analysis to compare the genotypes of S. Typhi cultured from acute cases of typhoid fever occurring in Santiago during the hyperendemic period of the 1980s (n = 74) versus the nonendemic 2010s (n = 80) when typhoid fever was rare. The genotype distribution between “historical” (1980s) isolates and “modern” (2011–2016) isolates was similar, with genotypes 3.5 and 2 comprising the majority of isolations, and 73/80 (91.3%) of modern isolates matching a genotype detected in the 1980s. Additionally, phylogenomically ‘ancient’ genotypes 1.1 and 1.2.1, uncommon in the global collections, were also detected in both eras, with a notable rise amongst the modern isolates. Thus, genotypes of S. Typhi causing acute illness in the modern nonendemic era match the genotypes circulating during the hyperendemic 1980s. The persistence of historical genotypes may be explained by chronic typhoid carriers originally infected during or before the 1980s.

Studies of Salmonella Typhi (the cause of typhoid fever) rarely include isolates collected both before and after the interruption of hyperendemic transmission because this typically occurred decades before modern bacteria preservation methods. After substantial reduction in disease, it was assumed that sporadic cases and infrequent outbreaks were due to either chronic biliary carriers or importations, but this was difficult to characterize with low resolution bacterial typing methods. In Santiago, Chile, typhoid fever persisted at hyperendemic levels through the 1980s until organized control efforts in the 1980s and changes to wastewater policy in 1991 caused annual typhoid incidence to plummet.

In this study, we used whole genome sequencing (WGS) to investigate whether recent sporadic cases occurring in Santiago in the 2010s were genomically similar to S. Typhi circulating in the 1980s, or dissimilar, possibly representing importations of S. Typhi from outside of Chile. We found concordance amongst S. Typhi genotypes between the 1980s and 2010s, and differences from genotypes circulating in Southeast Asia and Africa where typhoid remains hyperendemic. Our findings suggest that a proportion of modern, rare typhoid cases in Santiago are autochthonous, and that chronic carriers or another unknown reservoir likely contribute. Broadly, our findings corroborate the epidemiologic importance of long-term reservoirs of typhoid fever decades after typhoid elimination.

A large outbreak of COVID-19 in a UK prison, October 2020 to April 2021

Prisons are susceptible to outbreaks. Control measures focusing on isolation and cohorting negatively affect wellbeing. We present an outbreak of coronavirus disease 2019 (COVID-19) in a large male prison in Wales, UK, October 2020 to April 2021, and discuss control measures.We gathered case-information, including demographics, staff-residence postcode, resident cell number, work areas/dates, test results, staff interview dates/notes and resident prison-transfer dates. Epidemiological curves were mapped by prison location. Control measures included isolation (exclusion from work or cell-isolation), cohorting (new admissions and work-area groups), asymptomatic testing (case-finding), removal of communal dining and movement restrictions. Facemask use and enhanced hygiene were already in place. Whole-genome sequencing (WGS) and interviews determined the genetic relationship between cases plausibility of transmission.Of 453 cases, 53% (n = 242) were staff, most aged 25-34 years (11.5% females, 27.15% males) and symptomatic (64%). Crude attack-rate was higher in staff (29%, 95% CI 26-64%) than in residents (12%, 95% CI 9-15%).Whole-genome sequencing can help differentiate multiple introductions from person-to-person transmission in prisons. It should be introduced alongside asymptomatic testing as soon as possible to control prison outbreaks. Timely epidemiological investigation, including data visualisation, allowed dynamic risk assessment and proportionate control measures, minimising the reduction in resident welfare.

CABGen: A Web Application for the Bioinformatic Analysis of Bacterial Genomes

Due to recent developments in NGS technologies, genome sequencing is generating large volumes of new data containing a wealth of biological information. Understanding sequenced genomes in a biologically meaningful way and delineating their functional and metabolic landscapes is a first-level challenge. Considering the global antimicrobial resistance (AMR) problem, investments to expand surveillance and improve existing genome analysis technologies are pressing. In addition, the speed at which new genomic data is generated surpasses our capacity to analyze it with available bioinformatics methods, thus creating a need to develop new, user-friendly and comprehensive analytical tools. To this end, we propose a new web application, CABGen,¹ developed with open-source software. CABGen allows storing, organizing, analyzing, and interpreting bioinformatics data in a friendly, scalable, easy-to-use environment and can process data from bacterial isolates of different species and origins. CABGen has three modules: Upload Sequences, Analyze Sequences, and Verify Results. Functionalities include coverage estimation, species identification, de novo genome assembly, and assembly quality, genome annotation, MLST mapping, searches for genes related to AMR, virulence, and plasmids, and detection of point mutations in specific AMR genes. Visualization tools are also available, greatly facilitating the handling of biological data. The reports include those results that are clinically relevant. To illustrate the use of CABGen, whole-genome shotgun data from 181 bacterial isolates of different species collected in 5 Brazilian regions between 2018 and 2020 were uploaded and submitted to the platform’s modules.

STAMINA: Bioinformatics Platform for Monitoring and Mitigating Pandemic Outbreaks

This paper presents the components and integrated outcome of a system that aims to achieve early detection, monitoring and mitigation of pandemic outbreaks. The architecture of the platform aims at providing a number of pandemic-response-related services, on a modular basis, that allows for the easy customization of the platform to address user’s needs per case. This customization is achieved through its ability to deploy only the necessary, loosely coupled services and tools for each case, and by providing a common authentication, data storage and data exchange infrastructure. This way, the platform can provide the necessary services without the burden of additional services that are not of use in the current deployment (e.g., predictive models for pathogens that are not endemic to the deployment area). All the decisions taken for the communication and integration of the tools that compose the platform adhere to this basic principle. The tools presented here as well as their integration is part of the project STAMINA.

Population genomics confirms acquisition of drug-resistant Aspergillus fumigatus infection by humans from the environment

Infections caused by the fungal pathogen Aspergillus fumigatus are increasingly resistant to first-line azole antifungal drugs. However, despite its clinical importance, little is known about how susceptible patients acquire infection from drug-resistant genotypes in the environment. Here, we present a population genomic analysis of 218 A. fumigatus isolates from across the UK and Ireland (comprising 153 clinical isolates from 143 patients and 65 environmental isolates). First, phylogenomic analysis shows strong genetic structuring into two clades (A and B) with little interclade recombination and the majority of environmental azole resistance found within clade A. Second, we show occurrences where azole-resistant isolates of near-identical genotypes were obtained from both environmental and clinical sources, indicating with high confidence the infection of patients with resistant isolates transmitted from the environment. Third, genome-wide scans identified selective sweeps across multiple regions indicating a polygenic basis to the trait in some genetic backgrounds. These signatures of positive selection are seen for loci containing the canonical genes encoding fungicide resistance in the ergosterol biosynthetic pathway, while other regions under selection have no defined function. Lastly, pan-genome analysis identified genes linked to azole resistance and previously unknown resistance mechanisms. Understanding the environmental drivers and genetic basis of evolving fungal drug resistance needs urgent attention, especially in light of increasing numbers of patients with severe viral respiratory tract infections who are susceptible to opportunistic fungal superinfections.

Carbapenemase-Producing Extraintestinal Pathogenic Escherichia coli From Argentina: Clonal Diversity and Predominance of Hyperepidemic Clones CC10 and CC131

Extraintestinal pathogenic Escherichia coli (ExPEC) causes infections outside the intestine. Particular ExPEC clones, such as clonal complex (CC)/sequence type (ST)131, have been known to sequentially accumulate antimicrobial resistance that starts with chromosomal mutations against fluoroquinolones, followed with the acquisition of bla_CTX–M–15 and, more recently, carbapenemases. Here we aimed to investigate the distribution of global epidemic clones of carbapenemase-producing ExPEC from Argentina in representative clinical isolates recovered between July 2008 and March 2017. Carbapenemase-producing ExPEC (n = 160) were referred to the Argentinean reference laboratory. Of these, 71 were selected for genome sequencing. Phenotypic and microbiological studies confirmed the presence of carbapenemases confirmed as KPC-2 (n = 52), NDM-1 (n = 16), IMP-8 (n = 2), and VIM-1 (n = 1) producers. The isolates had been recovered mainly from urine, blood, and abdominal fluids among others, and some were from screening samples. After analyzing the virulence gene content, 76% of the isolates were considered ExPEC, although non-ExPEC isolates were also obtained from extraintestinal sites. Pan-genome phylogeny and clonal analysis showed great clonal diversity, although the first phylogroup in abundance was phylogroup A, harboring CC10 isolates, followed by phylogroup B2 with CC/ST131, mostly H30Rx, the subclone co-producing CTX-M-15. Phylogroups D, B1, C, F, and E were also detected with fewer strains. CC10 and CC/ST131 were found throughout the country. In addition, CC10 nucleated most metalloenzymes, such as NDM-1. Other relevant international clones were identified, such as CC/ST38, CC155, CC14/ST1193, and CC23. Two isolates co-produced KPC-2 and OXA-163 or OXA-439, a point mutation variant of OXA-163, and three isolates co-produced MCR-1 among other resistance genes. To conclude, in this work, we described the molecular epidemiology of carbapenemase-producing ExPEC in Argentina. Further studies are necessary to determine the plasmid families disseminating carbapenemases in ExPEC in this region.

IntegronFinder 2.0: Identification and Analysis of Integrons across Bacteria, with a Focus on Antibiotic Resistance in Klebsiella

Integrons are flexible gene-exchanging platforms that contain multiple cassettes encoding accessory genes whose order is shuffled by a specific integrase. Integrons embedded within mobile genetic elements often contain multiple antibiotic resistance genes that they spread among nosocomial pathogens and contribute to the current antibiotic resistance crisis. However, most integrons are presumably sedentary and encode a much broader diversity of functions. IntegronFinder is a widely used software to identify novel integrons in bacterial genomes, but has aged and lacks some useful functionalities to handle very large datasets of draft genomes or metagenomes. Here, we present IntegronFinder version 2. We have updated the code, improved its efficiency and usability, adapted the output to incomplete genome data, and added a few novel functions. We describe these changes and illustrate the relevance of the program by analyzing the distribution of integrons across more than 20,000 fully sequenced genomes. We also take full advantage of its novel capabilities to analyze close to 4000 Klebsiella pneumoniae genomes for the presence of integrons and antibiotic resistance genes within them. Our data show that K. pneumoniae has a large diversity of integrons and the largest mobile integron in our database of plasmids. The pangenome of these integrons contains a total of 165 different gene families with most of the largest families being related with resistance to numerous types of antibiotics. IntegronFinder is a free and open-source software available on multiple public platforms.

Horizontal Gene Transfer of Fluoroquinolone Resistance-Conferring Genes From Commensal Neisseria to Neisseria gonorrhoeae: A Global Phylogenetic Analysis of 20,047 Isolates

Antimicrobial resistance in Neisseria gonorrhoeae is an important global health concern. The genetically related commensal Neisseria act as a reservoir of resistance genes, and horizontal gene transfer (HGT) has been shown to play an important role in the genesis of resistance to cephalosporins and macrolides in N. gonorrhoeae. In this study, we evaluated if there was evidence of HGT in the genes gyrA/gyrB and parC/parE responsible for fluoroquinolone resistance. Even though the role of gyrB and parE in quinolone resistance is unclear, the subunits gyrB and parE were included as zoliflodacin, a promising new drug to treat N. gonorrhoeae targets the gyrB subunit. We analyzed a collection of 20,047 isolates; 18,800 N. gonorrhoeae, 1,238 commensal Neisseria spp., and nine Neisseria meningitidis. Comparative genomic analyses identified HGT events in genes, gyrA, gyrB, parC, and parE. Recombination events were predicted in N. gonorrhoeae and Neisseria commensals. Neisseria lactamica, Neisseria macacae, and Neisseria mucosa were identified as likely progenitors of the HGT events in gyrA, gyrB, and parE, respectively.

COVID-19 Variant Surveillance and Social Determinants in Central Massachusetts: Development Study (Preprint)

Background

Public health scientists have used spatial tools such as web-based Geographical Information System (GIS) applications to monitor and forecast the progression of the COVID-19 pandemic and track the impact of their interventions. The ability to track SARS-CoV-2 variants and incorporate the social determinants of health with street-level granularity can facilitate the identification of local outbreaks, highlight variant-specific geospatial epidemiology, and inform effective interventions. We developed a novel dashboard, the University of Massachusetts’ Graphical user interface for Geographic Information (MAGGI) variant tracking system that combines GIS, health-associated sociodemographic data, and viral genomic data to visualize the spatiotemporal incidence of SARS-CoV-2 variants with street-level resolution while safeguarding protected health information. The specificity and richness of the dashboard enhance the local understanding of variant introductions and transmissions so that appropriate public health strategies can be devised and evaluated.

Objective

We developed a web-based dashboard that simultaneously visualizes the geographic distribution of SARS-CoV-2 variants in Central Massachusetts, the social determinants of health, and vaccination data to support public health efforts to locally mitigate the impact of the COVID-19 pandemic.

Methods

MAGGI uses a server-client model–based system, enabling users to access data and visualizations via an encrypted web browser, thus securing patient health information. We integrated data from electronic medical records, SARS-CoV-2 genomic analysis, and public health resources. We developed the following functionalities into MAGGI: spatial and temporal selection capability by zip codes of interest, the detection of variant clusters, and a tool to display variant distribution by the social determinants of health. MAGGI was built on the Environmental Systems Research Institute ecosystem and is readily adaptable to monitor other infectious diseases and their variants in real-time.

Results

We created a geo-referenced database and added sociodemographic and viral genomic data to the ArcGIS dashboard that interactively displays Central Massachusetts’ spatiotemporal variants distribution. Genomic epidemiologists and public health officials use MAGGI to show the occurrence of SARS-CoV-2 genomic variants at high geographic resolution and refine the display by selecting a combination of data features such as variant subtype, subject zip codes, or date of COVID-19–positive sample collection. Furthermore, they use it to scale time and space to visualize association patterns between socioeconomics, social vulnerability based on the Centers for Disease Control and Prevention’s social vulnerability index, and vaccination rates. We launched the system at the University of Massachusetts Chan Medical School to support internal research projects starting in March 2021.

Conclusions

We developed a COVID-19 variant surveillance dashboard to advance our geospatial technologies to study SARS-CoV-2 variants transmission dynamics. This real-time, GIS-based tool exemplifies how spatial informatics can support public health officials, genomics epidemiologists, infectious disease specialists, and other researchers to track and study the spread patterns of SARS-CoV-2 variants in our communities.

A curated dataset of peste des petits ruminants virus sequences for molecular epidemiological analyses

Peste des petits ruminants (PPR) is a highly contagious and devastating viral disease infecting predominantly sheep and goats. Tracking outbreaks of disease and analysing the movement of the virus often involves sequencing part or all of the genome and comparing the sequence obtained with sequences from other outbreaks, obtained from the public databases. However, there are a very large number (>1800) of PPRV sequences in the databases, a large majority of them relatively short, and not always well-documented. There is also a strong bias in the composition of the dataset, with countries with good sequencing capabilities (e.g. China, India, Turkey) being overrepresented, and most sequences coming from isolates in the last 20 years. In order to facilitate future analyses, we have prepared sets of PPRV sequences, sets which have been filtered for sequencing errors and unnecessary duplicates, and for which date and location information has been obtained, either from the database entry or from other published sources. These sequence datasets are freely available for download, and include smaller datasets which maximise phylogenetic information from the minimum number of sequences, and which will be useful for simple lineage identification. Their utility is illustrated by uploading the data to the MicroReact platform to allow simultaneous viewing of lineage date and geographic information on all the viruses for which we have information. While preparing these datasets, we identified a significant number of public database entries which contain clear errors, and propose guidelines on checking new sequences and completing metadata before submission.

WHO Critical Priority Escherichia coli as One Health Challenge for a Post-Pandemic Scenario: Genomic Surveillance and Analysis of Current Trends in Brazil

The dissemination of carbapenem-resistant and third generation cephalosporin-resistant pathogens is a critical issue that is no longer restricted to hospital settings. The rapid spread of critical priority pathogens in Brazil is notably worrying, considering its continental dimension, the diversity of international trade, livestock production, and human travel. We conducted a nationwide genomic investigation under a One Health perspective that included Escherichia coli strains isolated from humans and nonhuman sources, over 45 years (1974–2019). One hundred sixty-seven genomes were analyzed extracting clinically relevant information (i.e., resistome, virulome, mobilome, sequence types [STs], and phylogenomic). The endemic status of extended-spectrum β-lactamase (ESBL)-positive strains carrying a wide diversity of bla_CTX-M variants, and the growing number of colistin-resistant isolates carrying mcr-type genes was associated with the successful expansion of international ST10, ST38, ST115, ST131, ST354, ST410, ST648, ST517, and ST711 clones; phylogenetically related and shared between human and nonhuman hosts, and polluted aquatic environments. Otherwise, carbapenem-resistant ST48, ST90, ST155, ST167, ST224, ST349, ST457, ST648, ST707, ST744, ST774, and ST2509 clones from human host harbored bla_KPC-2 and bla_NDM-1 genes. A broad resistome to other clinically relevant antibiotics, hazardous heavy metals, disinfectants, and pesticides was further predicted. Wide virulome associated with invasion/adherence, exotoxin and siderophore production was related to phylogroup B2. The convergence of wide resistome and virulome has contributed to the persistence and rapid spread of international high-risk clones of critical priority E. coli at the human-animal-environmental interface, which must be considered a One Health challenge for a post-pandemic scenario.

IMPORTANCE A One Health approach for antimicrobial resistance must integrate whole-genome sequencing surveillance data of critical priority pathogens from human, animal and environmental sources to track hot spots and routes of transmission and developing effective prevention and control strategies. As part of the Grand Challenges Explorations: New Approaches to Characterize the Global Burden of Antimicrobial Resistance Program, we present genomic data of WHO critical priority carbapenemase-resistant, ESBL-producing, and/or colistin-resistant Escherichia coli strains isolated from humans and nonhuman sources in Brazil, a country with continental proportions and high levels of antimicrobial resistance. The present study provided evidence of epidemiological and clinical interest, highlighting that the convergence of wide virulome and resistome has contributed to the persistence and rapid spread of international high-risk clones of E. coli at the human-animal-environmental interface, which must be considered a One Health threat that requires coordinated actions to reduce its incidence in humans and nonhuman hosts.

Global evolutionary dynamics and resistome analysis of Clostridioides difficile ribotype 017

Clostridioides difficile PCR ribotype (RT) 017 ranks among the most successful strains of C. difficile in the world. In the past three decades, it has caused outbreaks on four continents, more than other ‘epidemic’ strains, but our understanding of the genomic epidemiology underpinning the spread of C. difficile RT 017 is limited. Here, we performed high-resolution phylogenomic and Bayesian evolutionary analyses on an updated and more representative dataset of 282 non-clonal C. difficile RT 017 isolates collected worldwide between 1981 and 2019. These analyses place an estimated time of global dissemination between 1953 and 1983 and identified the acquisition of the ermB-positive transposon Tn6194 as a key factor behind global emergence. This coincided with the introduction of clindamycin, a key inciter of C. difficile infection, into clinical practice in the 1960s. Based on the genomic data alone, the origin of C. difficile RT 017 could not be determined; however, geographical data and records of population movement suggest that C. difficile RT 017 had been moving between Asia and Europe since the Middle Ages and was later transported to North America around 1860 (95 % confidence interval: 1622–1954). A focused epidemiological study of 45 clinical C. difficile RT 017 genomes from a cluster in a tertiary hospital in Thailand revealed that the population consisted of two groups of multidrug-resistant (MDR) C. difficile RT 017 and a group of early, non-MDR C. difficile RT 017. The significant genomic diversity within each MDR group suggests that although they were all isolated from hospitalized patients, there was probably a reservoir of C. difficile RT 017 in the community that contributed to the spread of this pathogen.

Whole genome sequence analysis of Salmonella Typhi in Papua New Guinea reveals an established population of genotype 2.1.7 sensitive to antimicrobials

BACKGROUND

Typhoid fever, a systemic infection caused by Salmonella enterica serovar Typhi, remains a considerable public health threat in impoverished regions within many low- and middle-income settings. However, we still lack a detailed understanding of the emergence, population structure, molecular mechanisms of antimicrobial resistance (AMR), and transmission dynamics of S. Typhi across many settings, particularly throughout the Asia-Pacific islands. Here we present a comprehensive whole genome sequence (WGS) based overview of S. Typhi populations circulating in Papua New Guinea (PNG) over 30 years.

PRINCIPLE FINDINGS

Bioinformatic analysis of 86 S. Typhi isolates collected between 1980-2010 demonstrated that the population structure of PNG is dominated by a single genotype (2.1.7) that appears to have emerged in the Indonesian archipelago in the mid-twentieth century with minimal evidence of inter-country transmission. Genotypic and phenotypic data demonstrated that the PNG S. Typhi population appears to be susceptible to former first line drugs for treating typhoid fever (chloramphenicol, ampicillin and co-trimoxazole), as well as fluoroquinolones, third generation cephalosporins, and macrolides. PNG genotype 2.1.7 was genetically conserved, with very few deletions, and no evidence of plasmid or prophage acquisition. Genetic variation among this population was attributed to either single point mutations, or homologous recombination adjacent to repetitive ribosomal RNA operons.

SIGNIFICANCE

Antimicrobials remain an effective option for the treatment of typhoid fever in PNG, along with other intervention strategies including improvements to water, sanitation and hygiene (WaSH) related infrastructure and potentially the introduction of Vi-conjugate vaccines. However, continued genomic surveillance is warranted to monitor for the emergence of AMR within local populations, or the introduction of AMR associated genotypes of S. Typhi in this setting.

Metagenomics Approaches for Improving Food Safety: A Review

ABSTRACT

Advancements in next-generation sequencing technology have dramatically reduced the cost and increased the ease of microbial whole genome sequencing. This approach is revolutionizing the identification and analysis of foodborne microbial pathogens, facilitating expedited detection and mitigation of foodborne outbreaks, improving public health outcomes, and limiting costly recalls. However, next-generation sequencing is still anchored in the traditional laboratory practice of the selection and culture of a single isolate. Metagenomic-based approaches, including metabarcoding and shotgun and long-read metagenomics, are part of the next disruptive revolution in food safety diagnostics and offer the potential to directly identify entire microbial communities in a single food, ingredient, or environmental sample. In this review, metagenomic-based approaches are introduced and placed within the context of conventional detection and diagnostic techniques, and essential considerations for undertaking metagenomic assays and data analysis are described. Recent applications of the use of metagenomics for food safety are discussed alongside current limitations and knowledge gaps and new opportunities arising from the use of this technology.

HIGHLIGHTS

International links between Streptococcus pneumoniae vaccine serotype 4 sequence type (ST) 801 in Northern European shipyard outbreaks of invasive pneumococcal disease

BACKGROUND

Pneumococcal disease outbreaks of vaccine preventable serotype 4 sequence type (ST)801 in shipyards have been reported in several countries. We aimed to use genomics to establish any international links between them.

METHODS

Sequence data from ST801-related outbreak isolates from Norway (n = 17), Finland (n = 11) and Northern Ireland (n = 2) were combined with invasive pneumococcal disease surveillance from the respective countries, and ST801-related genomes from an international collection (n = 41 of > 40,000), totalling 106 genomes. Raw data were mapped and recombination excluded before phylogenetic dating.

RESULTS

Outbreak isolates were relatively diverse, with up to 100 SNPs (single nucleotide polymorphisms) and a common ancestor estimated around the year 2000. However, 19 Norwegian and Finnish isolates were nearly indistinguishable (0-2 SNPs) with the common ancestor dated around 2017.

CONCLUSION

The total diversity of ST801 within the outbreaks could not be explained by recent transmission alone, suggesting that harsh environmental and associated living conditions reported in the shipyards may facilitate invasion of colonising pneumococci. However, near identical strains in the Norwegian and Finnish outbreaks does suggest that transmission between international shipyards also contributed to those outbreaks. This indicates the need for improved preventative measures in this working population including pneumococcal vaccination.

Genomic epidemiology of SARS-CoV-2 in a UK university identifies dynamics of transmission

Understanding SARS-CoV-2 transmission in higher education settings is important to limit spread between students, and into at-risk populations. In this study, we sequenced 482 SARS-CoV-2 isolates from the University of Cambridge from 5 October to 6 December 2020. We perform a detailed phylogenetic comparison with 972 isolates from the surrounding community, complemented with epidemiological and contact tracing data, to determine transmission dynamics. We observe limited viral introductions into the university; the majority of student cases were linked to a single genetic cluster, likely following social gatherings at a venue outside the university. We identify considerable onward transmission associated with student accommodation and courses; this was effectively contained using local infection control measures and following a national lockdown. Transmission clusters were largely segregated within the university or the community. Our study highlights key determinants of SARS-CoV-2 transmission and effective interventions in a higher education setting that will inform public health policy during pandemics.

Widespread sharing of pneumococcal strains in a rural African setting: proximate villages are more likely to share similar strains that are carried at multiple timepoints

The transmission dynamics of Streptococcus pneumoniae in sub-Saharan Africa are poorly understood due to a lack of adequate epidemiological and genomic data. Here we leverage a longitudinal cohort from 21 neighbouring villages in rural Africa to study how closely related strains of S. pneumoniae are shared among infants. We analysed 1074 pneumococcal genomes isolated from 102 infants from 21 villages. Strains were designated for unique serotype and sequence-type combinations, and we arbitrarily defined strain sharing where the pairwise genetic distance between strains could be accounted for by the mean within host intra-strain diversity. We used non-parametric statistical tests to assess the role of spatial distance and prolonged carriage on strain sharing using a logistic regression model. We recorded 458 carriage episodes including 318 (69.4 %) where the carried strain was shared with at least one other infant. The odds of strain sharing varied significantly across villages (χ²=47.5, df=21, P-value <0.001). Infants in close proximity to each other were more likely to be involved in strain sharing, but we also show a considerable amount of strain sharing across longer distances. Close geographic proximity (<5 km) between shared strains was associated with a significantly lower pairwise SNP distance compared to strains shared over longer distances (P-value <0.005). Sustained carriage of a shared strain among the infants was significantly more likely to occur if they resided in villages within a 5 km radius of each other (P-value <0.005, OR 3.7). Conversely, where both infants were transiently colonized by the shared strain, they were more likely to reside in villages separated by over 15 km (P-value <0.05, OR 1.5). PCV7 serotypes were rare (13.5 %) and were significantly less likely to be shared (P-value <0.001, OR −1.07). Strain sharing was more likely to occur over short geographical distances, especially where accompanied by sustained colonization. Our results show that strain sharing is a useful proxy for studying transmission dynamics in an under-sampled population with limited genomic data. This article contains data hosted by Microreact.

Unveiling the genomic potential of Pseudomonas type strains for discovering new natural products

Microbes host a huge variety of biosynthetic gene clusters that produce an immeasurable array of secondary metabolites with many different biological activities such as antimicrobial, anticarcinogenic and antiviral. Despite the complex task of isolating and characterizing novel natural products, microbial genomic strategies can be useful for carrying out these types of studies. However, although genomic-based research on secondary metabolism is on the increase, there is still a lack of reports focusing specifically on the genus Pseudomonas. In this work, we aimed (i) to unveil the main biosynthetic systems related to secondary metabolism in Pseudomonas type strains, (ii) to study the evolutionary processes that drive the diversification of their coding regions and (iii) to select Pseudomonas strains showing promising results in the search for useful natural products. We performed a comparative genomic study on 194 Pseudomonas species, paying special attention to the evolution and distribution of different classes of biosynthetic gene clusters and the coding features of antimicrobial peptides. Using EvoMining, a bioinformatic approach for studying evolutionary processes related to secondary metabolism, we sought to decipher the protein expansion of enzymes related to the lipid metabolism, which may have evolved toward the biosynthesis of novel secondary metabolites in Pseudomonas. The types of metabolites encoded in Pseudomonas type strains were predominantly non-ribosomal peptide synthetases, bacteriocins, N-acetylglutaminylglutamine amides and ß-lactones. Also, the evolution of genes related to secondary metabolites was found to coincide with Pseudomonas species diversification. Interestingly, only a few Pseudomonas species encode polyketide synthases, which are related to the lipid metabolism broadly distributed among bacteria. Thus, our EvoMining-based search may help to discover new types of secondary metabolite gene clusters in which lipid-related enzymes are involved. This work provides information about uncharacterized metabolites produced by Pseudomonas type strains, whose gene clusters have evolved in a species-specific way. Our results provide novel insight into the secondary metabolism of Pseudomonas and will serve as a basis for the prioritization of the isolated strains. This article contains data hosted by Microreact.

Wild Boars as Reservoir of Highly Virulent Clone of Hybrid Shiga Toxigenic and Enterotoxigenic Escherichia coli Responsible for Edema Disease, France

Edema disease is an often fatal enterotoxemia caused by specific strains of Shiga toxin–producing Escherichia coli (STEC) that affect primarily healthy, rapidly growing nursery pigs. Recently, outbreaks of edema disease have also emerged in France in wild boars. Analysis of STEC strains isolated from wild boars during 2013–2019 showed that they belonged to the serotype O139:H1 and were positive for both Stx2e and F18 fimbriae. However, in contrast to classical STEC O139:H1 strains circulating in pigs, they also possessed enterotoxin genes sta1 and stb, typical of enterotoxigenic E. coli. In addition, the strains contained a unique accessory genome composition and did not harbor antimicrobial-resistance genes, in contrast to domestic pig isolates. These data thus reveal that the emergence of edema disease in wild boars was caused by atypical hybrid of STEC and enterotoxigenic E. coli O139:H1, which so far has been restricted to the wildlife environment.

Dissemination Routes of Carbapenem and Pan-Aminoglycoside Resistance Mechanisms in Hospital and Urban Wastewater Canalizations of Ghana

Wastewater has a major role in antimicrobial resistance (AMR) dynamics and public health. The impact on AMR of wastewater flux at the community-hospital interface in low- and middle-income countries (LMICs) is poorly understood. Therefore, the present study analyzed the epidemiological scenario of resistance genes, mobile genetic elements (MGEs), and bacterial populations in wastewater around the Tamale metropolitan area (Ghana). Wastewater samples were collected from the drainage and canalizations before and after three hospitals and one urban waste treatment plant (UWTP). From all carbapenem/pan-aminoglycoside-resistant bacteria, 36 isolates were selected to determine bacterial species and phenotypical resistance profiles. Nanopore sequencing was used to screen resistance genes and plasmids, whereas, sequence types, resistome and plasmidome contents, pan-genome structures, and resistance gene variants were analyzed with Illumina sequencing. The combination of these sequencing data allowed for the resolution of the resistance gene-carrying platforms. Hospitals and the UWTP collected genetic and bacterial elements from community wastewater and amplified successful resistance gene-bacterium associations, which reached the community canalizations. Uncommon carbapenemase/β-lactamase gene variants, like bla_DIM-1, and novel variants, including bla_VIM-71, bla_CARB-53, and bla_CMY-172, were identified and seem to spread via clonal expansion of environmental Pseudomonas spp. However, bla_NDM-1, bla_CTX-M-15, and armA genes, among others, were associated with MGEs that allowed for their dissemination between environmental and clinical bacterial hosts. In conclusion, untreated hospital wastewater in Ghana is a hot spot for the emergence and spread of genes and gene-plasmid-bacterium associations that accelerate AMR, including to last-resort antibiotics. Urgent actions must be taken in wastewater management in LMICs in order to delay AMR expansion.

IMPORTANCE Antimicrobial resistance (AMR) is one the major threats to public health today, especially resistance to last-resort compounds for the treatment of critical infections, such as carbapenems and aminoglycosides. Innumerable works have focused on the clinical ambit of AMR, but studies addressing the impact of wastewater cycles on the emergence and dissemination of resistant bacteria are still limited. The lack of knowledge is even greater when referring to low- and middle-income countries, where there is an absence of accurate sanitary systems. Furthermore, the combination of short- and long-read sequencing has surpassed former technical limitations, allowing the complete characterization of resistance genes, mobile genetic platforms, plasmids, and bacteria. The present study deciphered the multiple elements and routes involved in AMR dynamics in wastewater canalizations and, therefore, in the local population of Tamale, providing the basis to adopt accurate control measures to preserve and promote public health.

Genomic Analysis and Antimicrobial Resistance of Campylobacter jejuni and Campylobacter coli in Peru

Campylobacter is the leading cause of human bacterial gastroenteritis worldwide and has a major impact on global public health. Whole Genome Sequencing (WGS) is a powerful tool applied in the study of foodborne pathogens. The objective of the present study was to apply WGS to determine the genetic diversity, virulence factors and determinants of antimicrobial resistance of the populations of C. jejuni and C. coli in Peru. A total of 129 Campylobacter strains (108 C. jejuni and 21 C. coli) were sequenced using Illumina Miseq platform. In silico MLST analysis identified a high genetic diversity among those strains with 30 sequence types (STs), several of them within 11 clonal complexes (CC) for C. jejuni, while the strains of C. coli belonged to a single CC with 8 different STs. Phylogeny analysis showed that Peruvian C. jejuni strains were divided into 2 clades with 5 populations, while C. coli formed a single clade with 4 populations. Furthermore, in silico analyses showed the presence of several genes associated with adherence, colonization and invasion among both species: cadF (83.7%), jlpA (81.4%), racR (100%), dnaJ (83.7%), pebA (83.7%), pldA (82.1%), porA (84.5%), ceuE (82.9%), ciaB (78.3%), iamB (86.8%), and flaC (100%). The majority (82.9%) of the Campylobacter strains carried the cdtABC operon which code for cytolethal distending toxin (CDT). Half of them (50.4%) carried genes associated with the presence of T6SS, while the frequency of genes associated with T4SS were relatively low (11.6%). Genetic markers associated with resistance to quinolones, tetracycline (tetO, tetW/N/W), beta-lactamases (bla_oxa–61), macrolides (A2075G in 23S rRNA) were found in 94.5, 21.7, 66.7, 6.2, 69.8, and 18.6% of strains, respectively. The cmeABC multidrug efflux operon was present in 78.3% of strains. This study highlights the importance of using WGS in the surveillance of emerging pathogens associated with foodborne diseases, providing genomic information on genetic diversity, virulence mechanisms and determinants of antimicrobial resistance. The description of several Campylobacter genotypes having many virulence factors and resistance to quinolones and tetracyclines circulating in Peru provides important information which helps in the monitoring, control and prevention strategies of this emerging pathogen in our country.

Genomic Characterization of Carbapenem-Non-susceptible Pseudomonas aeruginosa Clinical Isolates From Saudi Arabia Revealed a Global Dissemination of GES-5-Producing ST235 and VIM-2-Producing ST233 Sub-Lineages

Carbapenem-resistant P. aeruginosa has become a major clinical problem due to limited treatment options. However, studies assessing the trends in the molecular epidemiology and mechanisms of antibiotic resistance in this pathogen are lacking in Saudi Arabia. Here, we reported the genome characterization in a global context of carbapenem non-susceptible clinical isolates from a nationally representative survey. The antibiotic resistance profiles of the isolates (n = 635) collected over 14 months between March 2018 and April 2019 from different geographical regions of Saudi Arabia showed resistance rates to relevant β-lactams, aminoglycosides and quinolones ranging between 6.93 and 27.56%. Overall, 22.52% (143/635) of the isolates exhibited resistance to both imipenem and meropenem that were mainly explained by porin loss and efflux overexpression. However, 18.18% of resistant isolates harbored genes encoding GES (69.23%), VIM (23.07%), NDM (3.85%) or OXA-48-like (3.85%) carbapenemases. Most common GES-positive isolates produced GESs −5, −15 or −1 and all belonged to ST235 whereas the VIM-positive isolates produced mainly VIM-2 and belonged to ST233 or ST257. GES and VIM producers were detected at different sampling periods and in different surveyed regions. Interestingly, a genome-wide comparison revealed that the GES-positive ST235 and VIM-2-positive ST233 genomes sequenced in this study and those available through public databases from various locations worldwide, constituted each a phylogenetically closely related sub-lineage. Profiles of virulence determinants, antimicrobial resistance genes and associated mobile elements confirmed relatedness within each of these two different sub-lineages. Sequence analysis located the bla_GES gene in nearly all studied genomes (95.4%) in the same integrative conjugative element that also harbored the acc(6′)-Ib, aph(3′)-XV, aadA6, sul1, tet(G), and catB resistance genes while bla_VIM–2 in most (98.89%) ST233-positive genomes was co-located with aac(6′)-I1, dfrB-5, and aac(3′)-Id in the same class I integron. The study findings revealed the global spread of GES-5 ST235 and VIM-2 ST233 sub-lineages and highlighted the importance of routine detection of rare β-lactamases.

To catch a hijacker: abundance, evolution and genetic diversity of P4-like bacteriophage satellites

Bacteriophages (phages) are bacterial parasites that can themselves be parasitized by phage satellites. The molecular mechanisms used by satellites to hijack phages are sometimes understood in great detail, but the origins, abundance, distribution and composition of these elements are poorly known. Here, we show that P4-like elements are present in more than 30% of the genomes of Enterobacterales, and in almost half of those of Escherichia coli, sometimes in multiple distinct copies. We identified over 1000 P4-like elements with very conserved genetic organization of the core genome and a few hotspots with highly variable genes. These elements are never found in plasmids and have very little homology to known phages, suggesting an independent evolutionary origin. Instead, they are scattered across chromosomes, possibly because their integrases are often exchanged with other elements. The rooted phylogenies of hijacking functions are correlated and suggest longstanding coevolution. They also reveal broad host ranges in P4-like elements, as almost identical elements can be found in distinct bacterial genera. Our results show that P4-like phage satellites constitute a very distinct, widespread and ancient family of mobile genetic elements. They pave the way for studying the molecular evolution of antagonistic interactions between phages and their satellites. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

First National Genomic Epidemiological Study of Neisseria gonorrhoeae Strains Spreading Across Sweden in 2016

The increasing transmission and antimicrobial resistance (AMR) in Neisseria gonorrhoeae is a global health concern with worrying trends of decreasing susceptibility to also the last-line extended-spectrum cephalosporin (ESC) ceftriaxone. A dramatic increase of reported gonorrhea cases has been observed in Sweden from 2016 and onward. The aim of the present study was to comprehensively investigate the genomic epidemiology of all cultured N. gonorrhoeae isolates in Sweden during 2016, in conjunction with phenotypic AMR and clinical and epidemiological data of patients. In total, 1279 isolates were examined. Etest and whole-genome sequencing (WGS) were performed, and epidemiological data obtained from the Public Health Agency of Sweden. Overall, 51.1%, 1.7%, and 1.3% resistance to ciprofloxacin, cefixime, and azithromycin, respectively, was found. No isolates were resistant to ceftriaxone, however, 9.3% of isolates showed a decreased susceptibility to ceftriaxone and 10.5% to cefixime. In total, 44 penA alleles were found of which six were mosaic (n = 92). Using the typing schemes of MLST, NG-MAST, and NG-STAR; 133, 422, and 280 sequence types, respectively, and 93 NG-STAR clonal complexes were found. The phylogenomic analysis revealed two main lineages (A and B) with lineage A divided into two main sublineages (A1 and A2). Resistance and decreased susceptibility to ESCs and azithromycin and associated AMR determinants, such as mosaic penA and mosaic mtrD, were predominantly found in sublineage A2. Resistance to cefixime and azithromycin was more prevalent among heterosexuals and MSM, respectively, and both were predominantly spread through domestic transmission. Continuous surveillance of the spread and evolution of N. gonorrhoeae, including phenotypic AMR testing and WGS, is essential for enhanced knowledge regarding the dynamic evolution of N. gonorrhoeae and gonorrhea epidemiology.

SARS-CoV-2 Reverse Zoonoses to Pumas and Lions, South Africa

Reverse-zoonotic infections of severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) from humans to wildlife species internationally raise concern over the emergence of new variants in animals. A better understanding of the transmission dynamics and pathogenesis in susceptible species will mitigate the risk to humans and wildlife occurring in Africa. Here we report infection of an exotic puma (July 2020) and three African lions (July 2021) in the same private zoo in Johannesburg, South Africa. One Health genomic surveillance identified transmission of a Delta variant from a zookeeper to the three lions, similar to those circulating in humans in South Africa. One lion developed pneumonia while the other cases had mild infection. Both the puma and lions remained positive for SARS-CoV-2 RNA for up to 7 weeks.