Detection of a historic reservoir of bedaquiline/clofazimine resistance-associated variants in Mycobacterium tuberculosis.

BACKGROUND: Drug resistance in tuberculosis (TB) poses a major ongoing challenge to public health. The recent inclusion of bedaquiline into TB drug regimens has improved treatment outcomes, but this advance is threatened by the emergence of strains of Mycobacterium tuberculosis (Mtb) resistant to bedaquiline. Clinical bedaquiline resistance is most frequently conferred by off-target resistance-associated variants (RAVs) in the mmpR5 gene (Rv0678), the regulator of an efflux pump, which can also confer cross-resistance to clofazimine, another TB drug. METHODS: We compiled a dataset of 3682 Mtb genomes, including 180 carrying variants in mmpR5, and its immediate background (i.e. mmpR5 promoter and adjacent mmpL5 gene), that have been associated to borderline (henceforth intermediate) or confirmed resistance to bedaquiline. We characterised the occurrence of all nonsynonymous mutations in mmpR5 in this dataset and estimated, using time-resolved phylogenetic methods, the age of their emergence. RESULTS: We identified eight cases where RAVs were present in the genomes of strains collected prior to the use of bedaquiline in TB treatment regimes. Phylogenetic reconstruction points to multiple emergence events and circulation of RAVs in mmpR5, some estimated to predate the introduction of bedaquiline. However, epistatic interactions can complicate bedaquiline drug-susceptibility prediction from genetic sequence data. Indeed, in one clade, Ile67fs (a RAV when considered in isolation) was estimated to have emerged prior to the antibiotic era, together with a resistance reverting mmpL5 mutation. CONCLUSIONS: The presence of a pre-existing reservoir of Mtb strains carrying bedaquiline RAVs prior to its clinical use augments the need for rapid drug susceptibility testing and individualised regimen selection to safeguard the use of bedaquiline in TB care and control.

Reference genome assembly and annotation of two Bacillus cereus sensu lato strains from Etosha National Park, Namibia.

Bacillus cereus sensu lato (s.l.) poses health and security issues. Here, we report the reference genome assembly of two Bacillus cereus s.l. strains, isolated from Etosha National Park, Namibia (FFI_BCgr36 and FFI_BCgr46). These unique genomes open for better understanding of environmental diversity and improvements in detection of threatening species.

Evolution of Bordetella pertussis in the acellular vaccine era in Norway, 1996 to 2019.

We described the population structure of Bordetella pertussis (B. pertussis) in Norway from 1996 to 2019 and determined if there were evolutionary shifts and whether these correlated with changes in the childhood immunization program. We selected 180 B. pertussis isolates, 22 from the whole cell vaccine (WCV) era (1996-1997) and 158 from the acellular vaccine (ACV) era (1998-2019). We conducted whole genome sequencing and determined the distribution and frequency of allelic variants and temporal changes of ACV genes. Norwegian B. pertussis isolates were evenly distributed across a phylogenetic tree that included global strains. We identified seven different allelic profiles of ACV genes (A-F), in which profiles A1, A2, and B dominated (89%), all having pertussis toxin (ptxA) allele 1, pertussis toxin promoter (ptxP) allele 3, and pertactin (prn) allele 2 present. Isolates with ptxP1 and prn1 were not detected after 2007, whereas the prn2 allele likely emerged prior to 1972, and ptxP3 before the early 1980s. Allele conversions of ACV genes all occurred prior to the introduction of ACV. Sixteen percent of our isolates showed mutations within the prn gene. ACV and its booster doses (implemented for children in 2007 and adolescents in 2013) might have contributed to evolvement of a more uniform B. pertussis population, with recent circulating strains having ptxA1, ptxP3, and prn2 present, and an increasing number of prn mutations. These strains clearly deviate from ACV strains (ptxA1, ptxP1, prn1), and this could have implications for vaccine efficiency and, therefore, prevention and control of pertussis.

Antibiotic Treatment Regimes as a Driver of the Global Population Dynamics of a Major Gonorrhea Lineage.

The Neisseria gonorrhoeae multilocus sequence type (ST) 1901 is among the lineages most commonly associated with treatment failure. Here, we analyze a global collection of ST-1901 genomes to shed light on the emergence and spread of alleles associated with reduced susceptibility to extended-spectrum cephalosporins (ESCs). The genetic diversity of ST-1901 falls into a minor and a major clade, both of which were inferred to have originated in East Asia. The dispersal of the major clade from Asia happened in two separate waves expanding from ∼1987 and 1996, respectively. Both waves first reached North America, and from there spread to Europe and Oceania, with multiple secondary reintroductions to Asia. The ancestor of the second wave acquired the penA 34.001 allele, which significantly reduces susceptibility to ESCs. Our results suggest that the acquisition of this allele granted the second wave a fitness advantage at a time when ESCs became the key drug class used to treat gonorrhea. Following its establishment globally, the lineage has served as a reservoir for the repeated emergence of clones fully resistant to the ESC ceftriaxone, an essential drug for effective treatment of gonorrhea. We infer that the effective population sizes of both clades went into decline as treatment schemes shifted from fluoroquinolones via ESC monotherapy to dual therapy with ceftriaxone and azithromycin in Europe and the United States. Despite the inferred recent population size decline, the short evolutionary path from the penA 34.001 allele to alleles providing full ceftriaxone resistance is a cause of concern.

Sudden emergence of a Neisseria gonorrhoeae clade with reduced susceptibility to extended-spectrum cephalosporins, Norway.

Neisseria gonorrhoeae multilocus sequence type (ST)-7827 emerged in a dramatic fashion in Norway in the period 2016-2018. Here, we aim to shed light on the provenance and expansion of this ST. ST-7827 was found to be polyphyletic, but the majority of members belonged to a monophyletic clade we termed PopPUNK cluster 7827 (PC-7827). In Norway, both PC-7827 and ST-7827 isolates were almost exclusively isolated from men. Phylogeographical analyses demonstrated an Asian origin of the genogroup, with multiple inferred exports to Europe and the USA. The genogroup was uniformly resistant to fluoroquinolones, and associated with reduced susceptibility to both azithromycin and the extended-spectrum cephalosporins (ESCs) cefixime and ceftriaxone. From a genetic background including the penA allele 13.001, associated with reduced ESC susceptibility, we identified repeated events of acquisition of porB alleles associated with further reduction in ceftriaxone susceptibility. Transmission of the strain was significantly reduced in Norway in 2019, but our results indicate the existence of a recently established global reservoir. The worrisome drug-resistance profile and rapid emergence of PC-7827 calls for close monitoring of the situation.

Neisseria meningitidis: using genomics to understand diversity, evolution and pathogenesis.

Meningococcal disease remains an important cause of morbidity and death worldwide despite the development and increasing implementation of effective vaccines. Elimination of the disease is hampered by the enormous diversity and antigenic variability of the causative agent, Neisseria meningitidis, one of the most variable bacteria in nature. These features are attained mainly through high rates of horizontal gene transfer and alteration of protein expression through phase variation. The recent availability of whole-genome sequencing (WGS) of large-scale collections of N. meningitidis isolates from various origins, databases to facilitate storage and sharing of WGS data and the concomitant development of effective bioinformatics tools have led to a much more thorough understanding of the diversity of the species, its evolution and population structure and how virulent traits may emerge. Implementation of WGS is already contributing to enhanced epidemiological surveillance and is essential to ascertain the impact of vaccination strategies. This Review summarizes the recent advances provided by WGS studies in our understanding of the biology of N. meningitidis and the epidemiology of meningococcal disease.

Rapid phenotypic evolution in multidrug-resistant Klebsiella pneumoniae hospital outbreak strains.

Carbapenem-resistant Klebsiella pneumoniae (CRKP) increasingly cause high-mortality outbreaks in hospital settings globally. Following a patient fatality at a hospital in Beijing due to a blaKPC-2-positive CRKP infection, close monitoring was put in place over the course of 14 months to characterize all blaKPC-2-positive CRKP in circulation in the hospital. Whole genome sequences were generated for 100 isolates from blaKPC-2-positive isolates from infected patients, carriers and the hospital environment. Phylogenetic analyses identified a closely related cluster of 82 sequence type 11 (ST11) isolates circulating in the hospital for at least a year prior to admission of the index patient. The majority of inferred transmissions for these isolates involved patients in intensive care units. Whilst the 82 ST11 isolates collected during the surveillance effort all had closely related chromosomes, we observed extensive diversity in their antimicrobial resistance (AMR) phenotypes. We were able to reconstruct the major genomic changes underpinning this variation in AMR profiles, including multiple gains and losses of entire plasmids and recombination events between plasmids, including transposition of blaKPC-2. We also identified specific cases where variation in plasmid copy number correlated with the level of phenotypic resistance to drugs, suggesting that the number of resistance elements carried by a strain may play a role in determining the level of AMR. Our findings highlight the epidemiological value of whole genome sequencing for investigating multi-drug-resistant hospital infections and illustrate that standard typing schemes cannot capture the extraordinarily fast genome evolution of CRKP isolates.

Global expansion of Mycobacterium tuberculosis lineage 4 shaped by colonial migration and local adaptation.

On the basis of population genomic and phylogeographic analyses of 1669 Mycobacterium tuberculosis lineage 4 (L4) genomes, we find that dispersal of L4 has been completely dominated by historical migrations out of Europe. We demonstrate an intimate temporal relationship between European colonial expansion into Africa and the Americas and the spread of L4 tuberculosis (TB). Markedly, in the age of antibiotics, mutations conferring antimicrobial resistance overwhelmingly emerged locally (at the level of nations), with minimal cross-border transmission of resistance. The latter finding was found to reflect the relatively recent emergence of these mutations, as a similar degree of local restriction was observed for susceptible variants emerging on comparable time scales. The restricted international transmission of drug-resistant TB suggests that containment efforts at the level of individual countries could be successful.

Genome-based transmission modelling separates imported tuberculosis from recent transmission within an immigrant population.

In many countries the incidence of tuberculosis (TB) is low and is largely shaped by immigrant populations from high-burden countries. This is the case in Norway, where more than 80 % of TB cases are found among immigrants from high-incidence countries. A variable latent period, low rates of evolution and structured social networks make separating import from within-border transmission a major conundrum to TB control efforts in many low-incidence countries. Clinical Mycobacterium tuberculosis isolates belonging to an unusually large genotype cluster associated with people born in the Horn of Africa have been identified in Norway over the last two decades. We modelled transmission based on whole-genome sequence data to estimate infection times for individual patients. By contrasting these estimates with time of arrival in Norway, we estimate on a case-by-case basis whether patients were likely to have been infected before or after arrival. Independent import was responsible for the majority of cases, but we estimate that about one-quarter of the patients had contracted TB in Norway. This study illuminates the transmission dynamics within an immigrant community. Our approach is broadly applicable to many settings where TB control programmes can benefit from understanding when and where patients acquired TB.

Re-visiting the evolution, dispersal and epidemiology of Zika virus in Asia.

Based on serological evidence and viral isolation, Zika virus (ZIKV) has circulated for many years relatively benignly in a sylvatic cycle in Africa and an urban cycle in South East Asia (SEA). With the recent availability of limited but novel Indian ZIKV sequences to add to the plethora of SEA sequences, we traced the phylogenetic history and spatio-temporal dispersal pattern of ZIKV in Asia prior to its explosive emergence in the Pacific region and the Americas. These analyses demonstrated that the introduction and dispersal of ZIKV on the Pacific islands were preceded by an extended period of relatively silent transmission in SEA, enabling the virus to expand geographically and evolve adaptively before its unanticipated introduction to immunologically naive populations on the Pacific islands and in the Americas. Our findings reveal new features of the evolution and dispersal of this intriguing virus and may benefit future disease control strategies.

Characterizing the virome of Ixodes ricinus ticks from northern Europe.

RNA viruses are abundant infectious agents and present in all domains of life. Arthropods, including ticks, are well known as vectors of many viruses of concern for human and animal health. Despite their obvious importance, the extent and structure of viral diversity in ticks is still poorly understood, particularly in Europe. Using a bulk RNA-sequencing approach that captures the complete transcriptome, we analysed the virome of the most common tick in Europe - Ixodes ricinus. In total, RNA sequencing was performed on six libraries consisting of 33 I. ricinus nymphs and adults sampled in Norway. Despite the small number of animals surveyed, our virus identification pipeline revealed nine diverse and novel viral species, phylogenetically positioned within four different viral groups - bunyaviruses, luteoviruses, mononegavirales and partitiviruses - and sometimes characterized by extensive genetic diversity including a potentially novel genus of bunyaviruses. This work sheds new light on the virus diversity in I. ricinus, expands our knowledge of potential host/vector-associations and tick-transmitted viruses within several viral groups, and pushes the latitudinal limit where it is likely to find tick-associated viruses. Notably, our phylogenetic analysis revealed the presence of tick-specific virus clades that span multiple continents, highlighting the role of ticks as important virus reservoirs.

Armed conflict and population displacement as drivers of the evolution and dispersal of Mycobacterium tuberculosis.

The "Beijing" Mycobacterium tuberculosis (Mtb) lineage 2 (L2) is spreading globally and has been associated with accelerated disease progression and increased antibiotic resistance. Here we performed a phylodynamic reconstruction of one of the L2 sublineages, the central Asian clade (CAC), which has recently spread to western Europe. We find that recent historical events have contributed to the evolution and dispersal of the CAC. Our timing estimates indicate that the clade was likely introduced to Afghanistan during the 1979-1989 Soviet-Afghan war and spread further after population displacement in the wake of the American invasion in 2001. We also find that drug resistance mutations accumulated on a massive scale in Mtb isolates from former Soviet republics after the fall of the Soviet Union, a pattern that was not observed in CAC isolates from Afghanistan. Our results underscore the detrimental effects of political instability and population displacement on tuberculosis control and demonstrate the power of phylodynamic methods in exploring bacterial evolution in space and time.

An evolutionary analysis of genome expansion and pathogenicity in Escherichia coli.

BACKGROUND: There are several studies describing loss of genes through reductive evolution in microbes, but how selective forces are associated with genome expansion due to horizontal gene transfer (HGT) has not received similar attention. The aim of this study was therefore to examine how selective pressures influence genome expansion in 53 fully sequenced and assembled Escherichia coli strains. We also explored potential connections between genome expansion and the attainment of virulence factors. This was performed using estimations of several genomic parameters such as AT content, genomic drift (measured using relative entropy), genome size and estimated HGT size, which were subsequently compared to analogous parameters computed from the core genome consisting of 1729 genes common to the 53 E. coli strains. Moreover, we analyzed how selective pressures (quantified using relative entropy and dN/dS), acting on the E. coli core genome, influenced lineage and phylogroup formation. RESULTS: Hierarchical clustering of dS and dN estimations from the E. coli core genome resulted in phylogenetic trees with topologies in agreement with known E. coli taxonomy and phylogroups. High values of dS, compared to dN, indicate that the E. coli core genome has been subjected to substantial purifying selection over time; significantly more than the non-core part of the genome (p<0.001). This is further supported by a linear association between strain-wise dS and dN values (ß = 26.94 ± 0.44, R2~0.98, p<0.001). The non-core part of the genome was also significantly more AT-rich (p<0.001) than the core genome and E. coli genome size correlated with estimated HGT size (p<0.001). In addition, genome size (p<0.001), AT content (p<0.001) as well as estimated HGT size (p<0.005) were all associated with the presence of virulence factors, suggesting that pathogenicity traits in E. coli are largely attained through HGT. No associations were found between selective pressures operating on the E. coli core genome, as estimated using relative entropy, and genome size (p~0.98). CONCLUSIONS: On a larger time frame, genome expansion in E. coli, which is significantly associated with the acquisition of virulence factors, appears to be independent of selective forces operating on the core genome.