Transcriptomic and proteomic analysis of the virulence inducing effect of ciprofloxacin on enterohemorrhagic Escherichia coli.

Enterohemorrhagic E. coli (EHEC) is considered to be the most dangerous pathotype of E. coli, as it causes severe conditions such as hemorrhagic colitis (HC) and hemolytic uremic syndrome (HUS). Antibiotic treatment of EHEC infections is generally not recommended since it may promote the production of the Shiga toxin (Stx) and lead to worsened symptoms. This study explores how exposure to the fluoroquinolone ciprofloxacin reorganizes the transcriptome and proteome of EHEC O157:H7 strain EDL933, with special emphasis on virulence-associated factors. As expected, exposure to ciprofloxacin caused an extensive upregulation of SOS-response- and Stx-phage proteins, including Stx. A range of other virulence-associated factors were also upregulated, including many genes encoded by the LEE-pathogenicity island, the enterohemolysin gene (ehxA), as well as several genes and proteins involved in LPS production. However, a large proportion of the genes and proteins (17 and 8%, respectively) whose expression was upregulated upon ciprofloxacin exposure (17 and 8%, respectively) are not functionally assigned. This indicates a knowledge gap in our understanding of mechanisms involved in EHECs response to antibiotic-induced stress. Altogether, the results contribute to better understanding of how exposure to ciprofloxacin influences the virulome of EHEC and generates a knowledge base for further studies on how EHEC responds to antibiotic-induced stress. A deeper understanding on how EHEC responds to antibiotics will facilitate development of novel and safer treatments for EHEC infections.

Scoary2: rapid association of phenotypic multi-omics data with microbial pan-genomes.

Unraveling bacterial gene function drives progress in various areas, such as food production, pharmacology, and ecology. While omics technologies capture high-dimensional phenotypic data, linking them to genomic data is challenging, leaving 40-60% of bacterial genes undescribed. To address this bottleneck, we introduce Scoary2, an ultra-fast microbial genome-wide association studies (mGWAS) software. With its data exploration app and improved performance, Scoary2 is the first tool to enable the study of large phenotypic datasets using mGWAS. As proof of concept, we explore the metabolome of yogurts, each produced with a different Propionibacterium reichii strain and discover two genes affecting carnitine metabolism.

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.

Tracing the adaptive evolution of SARS-CoV-2 during vaccine roll-out in Norway.

Vaccination against SARS-CoV-2 has greatly mitigated the impact of the COVID-19 pandemic. However, concerns have been raised about the degree to which vaccination might drive the emergence and selection of immune escape mutations that will hamper the efficacy of the vaccines. In this study, we investigate whether vaccination impacted the micro-scale adaptive evolution of SARS-CoV-2 in the Oslo region of Norway, during the first nine months of 2021, a period in which the population went from near-zero to almost 90 per cent vaccine coverage in the population over 50 years old. Weekly aggregated data stratified by age on vaccine uptake and number of SARS-CoV-2 cases in the area were obtained from the National Immunization Registry and the Norwegian Surveillance System for Communicable Diseases, respectively. A total of 6,438 virus sequences (7.5 per cent of the registered cases) along with metadata were available. We used a causal-driven approach to investigate the relationship between vaccination progress and changes in the frequency of 362 mutations present in at least ten samples, conditioned on the emergence of new lineages, time, and population vaccination coverage. After validating our approach, we identified 21 positive and 12 negative connections between vaccination progress and mutation prevalence, and most of them were outside the Spike protein. We observed a tendency for the mutations that we identified as positively connected with vaccination to decrease as the vaccinated population increased. After modelling the fitness of different competing mutations in a population, we found that our observations could be explained by a clonal interference phenomenon in which high fitness mutations would be outcompeted by the emergence or introduction of other high-fitness mutations.

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.

Invasive meningococcal disease in Norway in the two decades before the COVID-19 pandemic.

OBJECTIVES: Disease caused by the bacterium Neisseria meningitidis remains a worldwide public health challenge, despite the steadily decreasing incidence in Western countries. The objective of this study was to explore the epidemiology of invasive meningococcal disease in Norway over the last two decades. DESIGN: All isolates sent to the National Reference Laboratory from patients with invasive meningococcal disease between the years 2000 and 2019 were analyzed using whole genome sequencing (total: 625). RESULTS: A five-fold decrease in case numbers occurred over this period, and the situation has gone from being dominated by serogroup B to one where serogroups Y and W are more prevalent. Concurrently, the mean age at infection has increased from 18 to 33 years. Among the 350 serogroup B isolates, 87% were an exact match or cross-reactive with one or both the currently available serogroup B vaccines, but the proportion decreased in the past decade. Core genome analyses revealed a high variation in the number of allelic differences accumulated in epidemiologically linked isolates to the point that near-identical isolates were found several years apart. CONCLUSION: Allelic distance is an imprecise metric for the degree of epidemiologic linkage between isolates in N. meningitidis.

Naturally occurring Neisseria gonorrhoeae can have large deletions in housekeeping gene abcZ, making them untypable with multilocus sequence typing.

The abcZ gene is an essential housekeeping gene in all the Neisseria species. It is one of the seven genes used for multilocus sequence typing (MLST) this genus. It encodes the cytosolic component of an ATP-binding cassette (ABC) transporter complex of unknown function. We report here the finding of a strain of Neisseria gonorrhoeae with a 485 base pair deletion in the 5' region of the abcZ gene that truncates the protein product from 636 amino acids to 89 amino acids. A second open reading frame (ORF), encoding the latter 388 amino acids of the abcZ gene, was predicted downstream. The deletion will affect MLST profiling; interrogation of genomic sequences from PubMLST revealed that this isolate is not an anomaly. Deletions in abcZ were identified in 256 Neisseria genomes, roughly 0.6% of isolates. Furthermore, these deletions could leave the abcZ gene in a pseudogenized state. Our strain, isolated from a patient with symptoms of gonorrheal infection, nevertheless behaved normal in terms of growth and in vitro phenotypic properties.

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.

Whole Genome Sequencing of Drug Resistant and Drug Susceptible Mycobacterium tuberculosis Isolates From Tigray Region, Ethiopia.

Background: Tuberculosis, mainly caused by Mycobacterium tuberculosis (Mtb), is an ancient human disease that gravely affects millions of people annually. We wanted to explore the genetic diversity and lineage-specific association of Mtb with drug resistance among pulmonary tuberculosis patients. Methods: Sputum samples were collected from pulmonary tuberculosis patients at six different healthcare institutions in Tigray, Ethiopia, between July 2018 and August 2019. DNA was extracted from 74 Mtb complex isolates for whole-genome sequencing (WGS). All genomes were typed and screened for mutations with known associations with antimicrobial resistance using in silico methods, and results were cross-verified with wet lab methods. Results: Lineage (L) 4 (55.8%) was predominant, followed by L3 (41.2%); L1 (1.5%) and L2 (1.5%) occurred rarely. The most frequently detected sublineage was CAS (38.2%), followed by Ural (29.4%), and Haarlem (11.8%). The recent transmission index (RTI) was relatively low. L4 and Ural strains were more resistant than the other strains to any anti-TB drug (P < 0.05). The most frequent mutations to RIF, INH, EMB, SM, PZA, ETH, FLQs, and 2nd-line injectable drugs occurred at rpoB S450L, katG S315T, embB M306I/V, rpsL K43R, pncA V139A, ethA M1R, gyrA D94G, and rrs A1401G, respectively. Disputed rpoB mutations were also shown in four (16%) of RIF-resistant isolates. Conclusion: Our WGS analysis revealed the presence of diverse Mtb genotypes. The presence of a significant proportion of disputed rpoB mutations highlighted the need to establish a WGS facility at the regional level to monitor drug-resistant mutations. This will help control the transmission of DR-TB and ultimately contribute to the attainment of 100% DST coverage for TB patients as per the End TB strategy.

Comprehensive Contact Tracing, Testing and Sequencing Show Limited Transmission of SARS-CoV-2 between Children in Schools in Norway, August 2020 to May 2021.

The role of children in the spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in schools has been a topic of controversy. In this study among school contacts of SARS-CoV-2 positive children in 43 contact-investigations, we investigated SARS-CoV-2 transmission in Norway, August 2020-May 2021. All participants were tested twice within seven to ten days, using SARS-CoV-2 PCR on home-sampled saliva. Positive samples were whole genome sequenced. Among the 559 child contacts, eight tested positive (1.4%, 95% CI 0.62-2.80), with no significant difference between primary (1.0%, 95% CI 0.27-2.53) and secondary schools (2.6%, 95% CI 0.70-6.39), p = 0.229, nor by viral strain, non-Alpha (1.4%, 95% CI 0.50-2.94) and Alpha variant (B.1.1.7) (1.7%, 95% CI 0.21-5.99), p = 0.665. One adult contact (1/100) tested positive. In 34 index cases, we detected 13 different SARS-CoV-2 Pango lineage variants, with B.1.1.7 being most frequent. In the eight contact-investigations with SARS-CoV-2 positive contacts, four had the same sequence identity as the index, one had no relation, and three were inconclusive. With mitigation measures in place, the spread of SARS-CoV-2 from children in schools is limited. By excluding contact-investigations with adult cases known at the time of enrolment, our data provide a valid estimate on the role of children in the transmission of SARS-CoV-2 in schools.

The impact of global lineage dynamics, border restrictions, and emergence of the B.1.1.7 lineage on the SARS-CoV-2 epidemic in Norway.

As the COVID-19 pandemic swept through an immunologically naïve human population, academics and public health professionals scrambled to establish methods and platforms for genomic surveillance and data sharing. This offered a rare opportunity to study the ecology and evolution of SARS-CoV-2 over the course of the ongoing pandemic. Here, we use population genetic and phylogenetic methodology to characterize the population dynamics of SARS-CoV-2 and reconstruct patterns of virus introductions and local transmission in Norway against this backdrop. The analyses demonstrated that the epidemic in Norway was largely import driven and characterized by the repeated introduction, establishment, and suppression of new transmission lineages. This pattern changed with the arrival of the B.1.1.7 lineage, which was able to establish a stable presence concomitant with the imposition of severe border restrictions.

Correction: Characteristics of pulmonary multidrug-resistant tuberculosis patients in Tigray Region, Ethiopia: A cross-sectional study.

[This corrects the article DOI: 10.1371/journal.pone.0236362.].

Evolution and emergence of multidrug-resistant Mycobacterium tuberculosis in Chisinau, Moldova.

The evolution and emergence of drug-resistant tuberculosis (TB) has been studied extensively in some contexts, but the ecological drivers of these two processes remain poorly understood. This study sought to describe the joint evolutionary and epidemiological histories of a novel multidrug-resistant Mycobacterium tuberculosis strain recently identified in the capital city of the Republic of Moldova (MDR Ural/4.2), where genomic surveillance of drug-resistant M. tuberculosis has been limited thus far. Using whole genome sequence data and Bayesian phylogenomic methods, we reconstruct the stepwise acquisition of drug resistance mutations in the MDR Ural/4.2 strain, estimate its historical bacterial population size over time, and infer the migration history of this strain between Eastern European countries. We infer that MDR Ural/4.2 likely evolved (via acquisition of rpoB S450L, which confers resistance to rifampin) in the early 1990s, during a period of social turmoil following Moldovan independence from the Soviet Union. This strain subsequently underwent substantial population size expansion in the early 2000s, at a time when national guidelines encouraged inpatient treatment of TB patients. We infer exportation of this strain and its isoniazid-resistant ancestral precursor from Moldova to neighbouring countries starting as early as 1985. Our findings suggest temporal and ecological associations between specific public health practices, including inpatient hospitalization of drug-resistant TB cases from the early 2000s until 2013, and the evolution of drug-resistant M. tuberculosis in Moldova. These findings underscore the need for regional coordination in TB control and expanded genomic surveillance efforts across Eastern Europe.

The effect of recombination on the evolution of a population of Neisseria meningitidis.

Neisseria meningitidis (the meningococcus) is a major human pathogen with a history of high invasive disease burden, particularly in sub-Saharan Africa. Our current understanding of the evolution of meningococcal genomes is limited by the rarity of large-scale genomic population studies and lack of in-depth investigation of the genomic events associated with routine pathogen transmission. Here, we fill this knowledge gap by a detailed analysis of 2839 meningococcal genomes obtained through a carriage study of over 50,000 samples collected systematically in Burkina Faso, West Africa, before, during, and after the serogroup A vaccine rollout, 2009-2012. Our findings indicate that the meningococcal genome is highly dynamic, with highly recombinant loci and frequent gene sharing across deeply separated lineages in a structured population. Furthermore, our findings illustrate how population structure can correlate with genome flexibility, as some lineages in Burkina Faso are orders of magnitude more recombinant than others. We also examine the effect of selection on the population, in particular how it is correlated with recombination. We find that recombination principally acts to prevent the accumulation of deleterious mutations, although we do also find an example of recombination acting to speed the adaptation of a gene. In general, we show the importance of recombination in the evolution of a geographically expansive population with deep population structure in a short timescale. This has important consequences for our ability to both foresee the outcomes of vaccination programs and, using surveillance data, predict when lineages of the meningococcus are likely to become a public health concern.

Endospore Appendages: a novel pilus superfamily from the endospores of pathogenic Bacilli.

Bacillus cereus sensu lato is a group of Gram-positive endospore-forming bacteria with high ecological diversity. Their endospores are decorated with micrometer-long appendages of unknown identity and function. Here, we isolate endospore appendages (Enas) from the food poisoning outbreak strain B. cereus NVH 0075-95 and find proteinaceous fibers of two main morphologies: S- and L-Ena. By using cryoEM and 3D helical reconstruction of S-Enas, we show these to represent a novel class of Gram-positive pili. S-Enas consist of single domain subunits with jellyroll topology that are laterally stacked by ß-sheet augmentation. S-Enas are longitudinally stabilized by disulfide bonding through N-terminal connector peptides that bridge the helical turns. Together, this results in flexible pili that are highly resistant to heat, drought, and chemical damage. Phylogenomic analysis reveals a ubiquitous presence of the ena-gene cluster in the B. cereus group, which include species of clinical, environmental, and food importance. We propose Enas to represent a new class of pili specifically adapted to the harsh conditions encountered by bacterial spores.

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.

Frequency and patterns of first- and second-line drug resistance-conferring mutations in Mycobacterium tuberculosis isolated from pulmonary tuberculosis patients in a cross-sectional study in Tigray Region, Ethiopia.

OBJECTIVES: Tuberculosis (TB) is a preventable and treatable infectious disease, but the continuing emergence and spread of multidrug-resistant TB is threatening global TB control efforts. This study aimed to describe the frequency and patterns of drug resistance-conferring mutations of Mycobacterium tuberculosis (MTB) isolates detected from pulmonary TB patients in Tigray Region, Ethiopia. METHODS: A cross-sectional study design was employed to collect sputum samples from pulmonary TB patients between July 2018 to August 2019. Culture and identification tests were done at Tigray Health Research Institute (THRI). Mutations conferring rifampicin (RIF), isoniazid (INH) and fluoroquinolone (FQ) resistance were determined in 227 MTB isolates using GenoType MTBDRplus and GenoType MTBDRsl. RESULTS: Mutations conferring resistance to RIF, INH and FQs were detected in 40/227 (17.6%), 41/227 (18.1%) and 2/38 (5.3%) MTB isolates, respectively. The majority of mutations for RIF, INH and FQs occurred at codons rpoB S531L (70%), katG S315T (78%) and gyrA D94Y/N (100%), respectively. This study revealed a significant number of unknown mutations in the rpoB, katG and inhA genes. CONCLUSION: High rates of mutations conferring resistance to RIF, INH and FQs were observed in this study. A large number of isolates showed unknown mutations, which require further DNA sequencing analysis. Periodic drug resistance surveillance and scaling-up of drug resistance testing facilities are imperative to prevent the transmission of drug-resistant TB in the community.

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.

Characteristics of pulmonary multidrug-resistant tuberculosis patients in Tigray Region, Ethiopia: A cross-sectional study.

BACKGROUND: Tuberculosis (TB) is among the top 10 causes of mortality and the first killer among infectious diseases worldwide. One of the factors fuelling the TB epidemic is the global rise of multidrug resistant TB (MDR-TB). The aim of this study was to determine the magnitude and factors associated with MDR-TB in the Tigray Region, Ethiopia. METHOD: This study employed a facility-based cross-sectional study design, which was conducted between July 2018 and August 2019. The inclusion criteria for the study participants were GeneXpert-positive who were not under treatment for TB, PTB patients' ≥15 years of age and who provided written informed consent. A total of 300 participants were enrolled in the study, with a structured questionnaire used to collect data on clinical, sociodemographic and behavioral factors. Sputum samples were collected and processed for acid-fast bacilli staining, culture and drug susceptibility testing. Drug susceptibility testing was performed using a line probe assay. Logistic regression was used to analyze associations between outcome and predictor variables. RESULTS: The overall proportion of MDR-TB was 16.7% (11.6% and 32.7% for new and previously treated patients, respectively). Of the total MDR-TB isolates, 5.3% were pre-XDR-TB. The proportion of MDR-TB/HIV co-infection was 21.1%. A previous history of TB treatment AOR 3.75; 95% CI (0.7-2.24), cigarette smoking AOR 6.09; CI (1.65-2.50) and patients who had an intermittent fever (AOR = 2.54, 95% CI = 1.21-5.4) were strongly associated with MDR-TB development. CONCLUSIONS: The magnitude of MDR-TB observed among new and previously treated patients is very alarming, which calls for an urgent need for intervention. The high proportion of MDR-TB among newly diagnosed cases indicates ongoing transmission, which suggests the need for enhanced TB control program performance to interrupt transmission. The increased proportion of MDR-TB among previously treated cases indicates a need for better patient management to prevent the evolution of drug resistance. Assessing the TB control program performance gaps and an optimal implementation of the WHO recommended priority actions for the management of drug-resistant TB, is imperative to help reduce the current high MDR-TB burden in the study region.