Longitudinal analysis within one hospital in sub-Saharan Africa over 20 years reveals repeated replacements of dominant clones of Klebsiella pneumoniae and stresses the importance to include temporal patterns for vaccine design considerations.

BACKGROUND: Infections caused by multidrug-resistant gram-negative bacteria present a severe threat to global public health. The WHO defines drug-resistant Klebsiella pneumoniae as a priority pathogen for which alternative treatments are needed given the limited treatment options and the rapid acquisition of novel resistance mechanisms by this species. Longitudinal descriptions of genomic epidemiology of Klebsiella pneumoniae can inform management strategies but data from sub-Saharan Africa are lacking. METHODS: We present a longitudinal analysis of all invasive K. pneumoniae isolates from a single hospital in Blantyre, Malawi, southern Africa, from 1998 to 2020, combining clinical data with genome sequence analysis of the isolates. RESULTS: We show that after a dramatic increase in the number of infections from 2016 K. pneumoniae becomes hyperendemic, driven by an increase in neonatal infections. Genomic data show repeated waves of clonal expansion of different, often ward-restricted, lineages, suggestive of hospital-associated transmission. We describe temporal trends in resistance and surface antigens, of relevance for vaccine development. CONCLUSIONS: Our data highlight a clear need for new interventions to prevent rather than treat K. pneumoniae infections in our setting. Whilst one option may be a vaccine, the majority of cases could be avoided by an increased focus on and investment in infection prevention and control measures, which would reduce all healthcare-associated infections and not just one.

Aedes aegypti gut transcriptomes respond differently to microbiome transplants from field-caught or laboratory-reared mosquitoes.

The mosquito microbiome is critical for host development and plays a major role in many aspects of mosquito biology. While the microbiome is commonly dominated by a small number of genera, there is considerable variation in composition among mosquito species, life stages, and geography. How the host controls and is affected by this variation is unclear. Using microbiome transplant experiments, we asked whether there were differences in transcriptional responses when mosquitoes of different species were used as microbiome donors. We used microbiomes from four different donor species spanning the phylogenetic breadth of the Culicidae, collected either from the laboratory or the field. We found that when recipients received a microbiome from a donor reared in the laboratory, the response was remarkably similar regardless of donor species. However, when the donor had been collected from the field, many more genes were differentially expressed. We also found that while the transplant procedure did have some effect on the host transcriptome, this is likely to have had a limited effect on mosquito fitness. Overall, our results highlight the possibility that variation in mosquito microbiome communities is associated with variability in host-microbiome interactions and further demonstrate the utility of the microbiome transplantation technique for investigating host-microbe interactions in mosquitoes.

Mechanisms of transcriptional regulation in Anopheles gambiae revealed by allele specific expression.

Malaria control relies on insecticides targeting the mosquito vector, but this is increasingly compromised by insecticide resistance, which can be achieved by elevated expression of detoxifying enzymes that metabolize the insecticide. In diploid organisms, gene expression is regulated both in cis, by regulatory sequences on the same chromosome, and by trans acting factors, affecting both alleles equally. Differing levels of transcription can be caused by mutations in cis-regulatory modules (CRM), but few of these have been identified in mosquitoes. We crossed bendiocarb resistant and susceptible Anopheles gambiae strains to identify cis-regulated genes that might be responsible for the resistant phenotype using RNAseq, and cis-regulatory module sequences controlling gene expression in insecticide resistance relevant tissues were predicted using machine learning. We found 115 genes showing allele specific expression in hybrids of insecticide susceptible and resistant strains, suggesting cis regulation is an important mechanism of gene expression regulation in Anopheles gambiae. The genes showing allele specific expression included a higher proportion of Anopheles specific genes on average younger than genes those with balanced allelic expression.

Recovery of metagenomic data from the Aedes aegypti microbiome using a reproducible snakemake pipeline: MINUUR.

Background: Ongoing research of the mosquito microbiome aims to uncover novel strategies to reduce pathogen transmission. Sequencing costs, especially for metagenomics, are however still significant. A resource that is increasingly used to gain insights into host-associated microbiomes is the large amount of publicly available genomic data based on whole organisms like mosquitoes, which includes sequencing reads of the host-associated microbes and provides the opportunity to gain additional value from these initially host-focused sequencing projects. Methods: To analyse non-host reads from existing genomic data, we developed a snakemake workflow called MINUUR (Microbial INsights Using Unmapped Reads). Within MINUUR, reads derived from the host-associated microbiome were extracted and characterised using taxonomic classifications and metagenome assembly followed by binning and quality assessment. We applied this pipeline to five publicly available Aedes aegypti genomic datasets, consisting of 62 samples with a broad range of sequencing depths. Results: We demonstrate that MINUUR recovers previously identified phyla and genera and is able to extract bacterial metagenome assembled genomes (MAGs) associated to the microbiome. Of these MAGS, 42 are high-quality representatives with >90% completeness and <5% contamination. These MAGs improve the genomic representation of the mosquito microbiome and can be used to facilitate genomic investigation of key genes of interest. Furthermore, we show that samples with a high number of KRAKEN2 assigned reads produce more MAGs. Conclusions: Our metagenomics workflow, MINUUR, was applied to a range of Aedes aegypti genomic samples to characterise microbiome-associated reads. We confirm the presence of key mosquito-associated symbionts that have previously been identified in other studies and recovered high-quality bacterial MAGs. In addition, MINUUR and its associated documentation are freely available on GitHub and provide researchers with a convenient workflow to investigate microbiome data included in the sequencing data for any applicable host genome of interest.

Aedes aegypti gut transcriptomes respond differently to microbiome transplants from field-caught or laboratory-reared mosquitoes.

The mosquito microbiome is critical for host development and plays a major role in many aspects of mosquito biology. While the microbiome is commonly dominated by a small number of genera, there is considerable variation in composition among mosquito species, life stages, and geography. How the host controls and is affected by this variation is unclear. Using microbiome transplant experiments, we asked whether there were differences in transcriptional responses when mosquitoes of different species were used as microbiome donors. We used microbiomes from four different donor species spanning the phylogenetic breadth of the Culicidae, collected either from the laboratory or field. We found that when recipients received a microbiome from a donor reared in the laboratory, the response was remarkably similar regardless of donor species. However, when the donor had been collected from the field, far more genes were differentially expressed. We also found that while the transplant procedure did have some effect on the host transcriptome, this is likely to have had a limited effect on mosquito fitness. Overall, our results highlight the possibility that variation in mosquito microbiome communities are associated with variability in host-microbiome interactions and further demonstrate the utility of the microbiome transplantation technique.

Colonization dynamics of extended-spectrum beta-lactamase-producing Enterobacterales in the gut of Malawian adults.

Drug-resistant bacteria of the order Enterobacterales which produce extended-spectrum beta-lactamase enzymes (ESBL-Enterobacterales, ESBL-E) are global priority pathogens. Antimicrobial stewardship interventions proposed to curb their spread include shorter courses of antimicrobials to reduce selection pressure but individual-level acquisition and selection dynamics are poorly understood. We sampled stool of 425 adults (aged 16-76 years) in Blantyre, Malawi, over 6 months and used multistate modelling and whole-genome sequencing to understand colonization dynamics of ESBL-E. Models suggest a prolonged effect of antimicrobials such that truncating an antimicrobial course at 2 days has a limited effect in reducing colonization. Genomic analysis shows largely indistinguishable diversity of healthcare-associated and community-acquired isolates, hence some apparent acquisition of ESBL-E during hospitalization may instead represent selection from a patient's microbiota by antimicrobial exposure. Our approach could help guide stewardship protocols; interventions that aim to review and truncate courses of unneeded antimicrobials may be of limited use in preventing ESBL-E colonization.

The genus Serratia revisited by genomics.

The genus Serratia has been studied for over a century and includes clinically-important and diverse environmental members. Despite this, there is a paucity of genomic information across the genus and a robust whole genome-based phylogenetic framework is lacking. Here, we have assembled and analysed a representative set of 664 genomes from across the genus, including 215 historic isolates originally used in defining the genus. Phylogenomic analysis of the genus reveals a clearly-defined population structure which displays deep divisions and aligns with ecological niche, as well as striking congruence between historical biochemical phenotyping data and contemporary genomics data. We highlight the genomic, phenotypic and plasmid diversity of Serratia, and provide evidence of different patterns of gene flow across the genus. Our work provides a framework for understanding the emergence of clinical and other lineages of Serratia.

Vertical and Horizontal Transmission of Cell Fusing Agent Virus in Aedes aegypti.

Cell fusing agent virus (CFAV) is an insect-specific flavivirus (ISF) found in Aedes aegypti mosquitoes. ISFs have demonstrated the ability to modulate the infection or transmission of arboviruses such as dengue, West Nile, and Zika viruses. It is thought that vertical transmission is the main route for ISF maintenance in nature. This has been observed with CFAV, but there is evidence of horizontal and venereal transmission in other ISFs. Understanding the route of transmission can inform strategies to spread ISFs to vector populations as a method of controlling pathogenic arboviruses. We crossed individually reared male and female mosquitoes from both a naturally occurring CFAV-positive Ae. aegypti colony and its negative counterpart to provide information on maternal, paternal, and horizontal transmission. RT-PCR was used to detect CFAV in individual female pupal exuviae and was 89% sensitive, but only 42% in male pupal exuviae. This is a possible way to screen individuals for infection without destroying the adults. Female-to-male horizontal transmission was not observed during this study. However, there was a 31% transmission rate from mating pairs of CFAV-positive males to negative female mosquitoes. Maternal vertical transmission was observed with a filial infection rate of 93%. The rate of paternal transmission was 85% when the female remained negative, 61% when the female acquired CFAV horizontally, and 76% overall. Maternal and paternal transmission of CFAV could allow the introduction of this virus into wild Ae. aegypti populations through male or female mosquito releases, and thus provides a potential strategy for ISF-derived arbovirus control. IMPORTANCE Insect-specific flaviviruses (ISFs), are a group of nonpathogenic flaviviruses that only infect insects. ISFs can have a high prevalence in mosquito populations, but their transmission routes are not well understood. The results of this study confirm maternal transmission of cell fusing agent virus (CFAV) and demonstrate that paternal transmission is also highly efficient. Horizontal transmission of CFAV was also observed, aided by evaluation of the pupal infection status before mating with an infected individual. This technique of detecting infection in discarded pupae exuviae has not been evaluated previously and will be a useful tool for others in the field of studying viral transmission in mosquitoes. Identifying these routes of transmission provides information about how CFAV could be maintained in wild populations of mosquitoes and can aid future studies focusing on interactions of CFAV with their hosts and other viruses that infect mosquitoes.

Multiple Plasmodium falciparum drug resistance polymorphisms identified in a pregnant woman with severe malaria and a concomitant spontaneous abortion in Cross River, Nigeria, West Africa.

BACKGROUND: The development of resistance by Plasmodium falciparum to anti-malarial drugs impedes any benefits of the drug. In addition, absence or delayed availability of current anti-malarial drugs in remote areas has the potential to results to parasite escape and continuous transmission. CASE PRESENTATION: The case of a 29-year old pregnant woman from Biase Local Government Area in Cross River State Nigeria presenting with febrile illness and high body temperature of 38.7 °C was reported. She looked pale and vomited twice on arrival at the health facility. Her blood smear on the first day of hospitalization was positive for P. falciparum by RDT, microscopy (21,960 parasite/µl) and real-time PCR, with a PCV of 18%. She was treated with 600 mg intravenous quinine in 500 ml of 5% Dextrose/0.9% Saline 8-hourly for 24 h. On the second day of hospitalization, she complained of weakness, persistent high-grade fever and vaginal bleeding. A bulging amnion from an extended cervix was observed. Following venous blood collection for laboratory investigations, 600 µg of misoprostol was inserted into the posterior fornix of her vagina as part of her obstetric care. Parenteral quinine was discontinued, and she was given full therapeutic regimen of artemether-lumefantrine 80/480 mg tablets to be taken for 3 days beginning from the second day. Her blood samples on the second and third day of hospitalization remained positive for P. falciparum by all three diagnostic methods. Single nucleotide polymorphism (SNP) assay on all three P. falciparum isolates revealed the presence of variants associated with multiple drug resistant markers. DISCUSSION: Infecting P. falciparum isolates may have been resistant to initial quinine treatment resulting from parasite cross-resistance with other quinoline associated resistant markers such as 86Y and 184 F. CONCLUSIONS: Therefore, the likely transmission of similarly resistant parasites in the study area calls for reinforcement of interventions and adherence to current World Health Organization guidelines in administering only approved drugs to individuals in order to mitigate parasite escape and eventual transmission to other susceptible individuals.

Design-oriented evaluation of the hydrodynamics in a full-scale combined filter-lamella separator for urban stormwater treatment.

The development of compact treatment devices with high removal efficiencies and low space requirements is a key objective of urban stormwater treatment. Thus, many devices utilize a combination of sedimentation and upward flow filtration in a single system. This study, for the first time, evaluates the flow field inside a combined filter-lamella separator via computational fluid dynamics. Herein, three objectives are investigated: (i) the flow field for different structural configurations, (ii) the distribution of particulate matter along the filter bed and (iii) the dynamic clogging in discrete filter zones, which is addressed by a clogging model derived from literature data. The results indicate that a direct combination of a filtration stage with a lamella separator promotes a uniform flow distribution. The distribution of particulate matter along the filter bed varies with configuration and particle size. Clogging, induced by particles in the spectrum <63 µm, creates gradients of hydraulic conductivity along the filter bed. After treating about half of Germany's annual runoff-efficient precipitation at a rainfall intensity of 5 L/(s·ha), the filtration rates increase in the front of the filter bed by +10%. Thus, long-term operating behavior is sensitive to efficient filter utilization in compact treatment devices.

Piperacillin/tazobactam-resistant, cephalosporin-susceptible Escherichia coli bloodstream infections are driven by multiple acquisition of resistance across diverse sequence types.

Resistance to piperacillin/tazobactam (TZP) in Escherichia coli has predominantly been associated with mechanisms that confer resistance to third-generation cephalosporins. Recent reports have identified E. coli strains with phenotypic resistance to piperacillin/tazobactam but susceptibility to third-generation cephalosporins (TZP-R/3GC-S). In this study we sought to determine the genetic diversity of this phenotype in E. coli (n=58) isolated between 2014-2017 at a single tertiary hospital in Liverpool, UK, as well as the associated resistance mechanisms. We compare our findings to a UK-wide collection of invasive E. coli isolates (n=1509) with publicly available phenotypic and genotypic data. These data sets included the TZP-R/3GC-S phenotype (n=68), and piperacillin/tazobactam and third-generation cephalosporin-susceptible (TZP-S/3GC-S, n=1271) phenotypes. The TZP-R/3GC-S phenotype was displayed in a broad range of sequence types, which was mirrored in the same phenotype from the UK-wide collection, and the overall diversity of invasive E. coli isolates. The TZP-R/3GC-S isolates contained a diverse range of plasmids, indicating multiple acquisition events of TZP resistance mechanisms rather than clonal expansion of a particular plasmid or sequence type. The putative resistance mechanisms were equally diverse, including hyperproduction of TEM-1, either via strong promoters or gene amplification, carriage of inhibitor-resistant ß-lactamases, and an S133G blaCTX-M-15 mutation detected for the first time in clinical isolates. Several of these mechanisms were present at a lower abundance in the TZP-S/3GC-S isolates from the UK-wide collection, but without the associated phenotypic resistance to TZP. Eleven (19%) of the isolates had no putative mechanism identified from the genomic data. Our findings highlight the complexity of this cryptic phenotype and the need for continued phenotypic monitoring, as well as further investigation to improve detection and prediction of the TZP-R/3GC-S phenotype from genomic data.

Genomic and antigenic diversity of colonizing Klebsiella pneumoniae isolates mirrors that of invasive isolates in Blantyre, Malawi.

Members of the Klebsiella pneumoniae species complex, particularly K. pneumoniae subsp. pneumoniae are antimicrobial resistance (AMR) associated pathogens of global importance, and polyvalent vaccines targeting Klebsiella O-antigens are in development. Whole-genome sequencing has provided insight into O-antigen distribution in the K. pneumoniae species complex, as well as population structure and virulence determinants, but genomes from sub-Saharan Africa are underrepresented in global sequencing efforts. We therefore carried out a genomic analysis of extended-spectrum beta-lactamase (ESBL)-producing K. pneumoniae species complex isolates colonizing adults in Blantyre, Malawi. We placed these isolates in a global genomic context, and compared colonizing to invasive isolates from the main public hospital in Blantyre. In total, 203 isolates from stool and rectal swabs from adults were whole-genome sequenced and compared to a publicly available multicounty collection and previously sequenced Malawian and Kenyan isolates from blood or sterile sites. We inferred phylogenetic relationships and analysed the diversity of genetic loci linked to AMR, virulence, capsule and LPS O-antigen (O-types). We find that the diversity of Malawian K. pneumoniae subsp. pneumoniae isolates represents the species' population structure, but shows distinct local signatures concerning clonal expansions. Siderophore and hypermucoidy genes were more frequent in invasive versus colonizing isolates (present in 13 % vs 1 %) but still generally lacking in most invasive isolates. O-antigen population structure and distribution was similar in invasive and colonizing isolates, with O4 more common (14%) than in previously published studies (2-5 %). We conclude that host factors, pathogen opportunity or alternate virulence loci not linked to invasive disease elsewhere are likely to be the major determinants of invasive disease in Malawi. Distinct ST and O-type distributions in Malawi highlight the need to sample locations where the burden of invasive Klebsiella disease is greatest to robustly define secular trends in Klebsiella diversity to assist in the development of a useful vaccine. Colonizing and invasive isolates in Blantyre are similar, hence O-typing of colonizing Klebsiella isolates may be a rapid and cost-effective approach to describe global diversity and guide vaccine development.

Spatial and Genomic Data to Characterize Endemic Typhoid Transmission.

BACKGROUND: Diverse environmental exposures and risk factors have been implicated in the transmission of Salmonella Typhi, but the dominant transmission pathways through the environment to susceptible humans remain unknown. Here, we use spatial, bacterial genomic, and hydrological data to refine our view of typhoid transmission in an endemic setting. METHODS: A total of 546 patients presenting to Queen Elizabeth Central Hospital in Blantyre, Malawi, with blood culture-confirmed typhoid fever between April 2015 and January 2017 were recruited to a cohort study. The households of a subset of these patients were geolocated, and 256 S. Typhi isolates were whole-genome sequenced. Pairwise single-nucleotide variant distances were incorporated into a geostatistical modeling framework using multidimensional scaling. RESULTS: Typhoid fever was not evenly distributed across Blantyre, with estimated minimum incidence ranging across the city from <15 to >100 cases per 100 000 population per year. Pairwise single-nucleotide variant distance and physical household distances were significantly correlated (P = .001). We evaluated the ability of river catchment to explain the spatial patterns of genomics observed, finding that it significantly improved the fit of the model (P = .003). We also found spatial correlation at a smaller spatial scale, of households living <192 m apart. CONCLUSIONS: These findings reinforce the emerging view that hydrological systems play a key role in the transmission of typhoid fever. By combining genomic and spatial data, we show how multifaceted data can be used to identify high incidence areas, explain the connections between them, and inform targeted environmental surveillance, all of which will be critical to shape local and regional typhoid control strategies.

Integration of whole genome sequencing and transcriptomics reveals a complex picture of the reestablishment of insecticide resistance in the major malaria vector Anopheles coluzzii.

Insecticide resistance is a major threat to gains in malaria control, which have been stalling and potentially reversing since 2015. Studies into the causal mechanisms of insecticide resistance are painting an increasingly complicated picture, underlining the need to design and implement targeted studies on this phenotype. In this study, we compare three populations of the major malaria vector An. coluzzii: a susceptible and two resistant colonies with the same genetic background. The original colonised resistant population rapidly lost resistance over a 6-month period, a subset of this population was reselected with pyrethroids, and a third population of this colony that did not lose resistance was also available. The original resistant, susceptible and re-selected colonies were subject to RNAseq and whole genome sequencing, which identified a number of changes across the transcriptome and genome linked with resistance. Firstly, an increase in the expression of genes within the oxidative phosphorylation pathway were seen in both resistant populations compared to the susceptible control; this translated phenotypically through an increased respiratory rate, indicating that elevated metabolism is linked directly with resistance. Genome sequencing highlighted several blocks clearly associated with resistance, including the 2Rb inversion. Finally, changes in the microbiome profile were seen, indicating that the microbial composition may play a role in the resistance phenotype. Taken together, this study reveals a highly complicated phenotype in which multiple transcriptomic, genomic and microbiome changes combine to result in insecticide resistance.

Loss of YhcB results in dysregulation of coordinated peptidoglycan, LPS and phospholipid synthesis during Escherichia coli cell growth.

The cell envelope is essential for viability in all domains of life. It retains enzymes and substrates within a confined space while providing a protective barrier to the external environment. Destabilising the envelope of bacterial pathogens is a common strategy employed by antimicrobial treatment. However, even in one of the best studied organisms, Escherichia coli, there remain gaps in our understanding of how the synthesis of the successive layers of the cell envelope are coordinated during growth and cell division. Here, we used a whole-genome phenotypic screen to identify mutants with a defective cell envelope. We report that loss of yhcB, a conserved gene of unknown function, results in loss of envelope stability, increased cell permeability and dysregulated control of cell size. Using whole genome transposon mutagenesis strategies, we report the comprehensive genetic interaction network of yhcB, revealing all genes with a synthetic negative and a synthetic positive relationship. These genes include those previously reported to have a role in cell envelope biogenesis. Surprisingly, we identified genes previously annotated as essential that became non-essential in a ΔyhcB background. Subsequent analyses suggest that YhcB functions at the junction of several envelope biosynthetic pathways coordinating the spatiotemporal growth of the cell, highlighting YhcB as an as yet unexplored antimicrobial target.

Genomic investigation of a suspected Klebsiella pneumoniae outbreak in a neonatal care unit in sub-Saharan Africa.

A special-care neonatal unit from a large public hospital in Malawi was noted as having more frequent, difficult-to-treat infections, and a suspected outbreak of multi-drug-resistant Klebsiella pneumoniae was investigated using genomic characterisation. All K. pneumoniae bloodstream infections (BSIs) from patients in the neonatal ward (n=62), and a subset of K. pneumoniae BSI isolates (n=38) from other paediatric wards in the hospital, collected over a 4 year period were studied. After whole genome sequencing, the strain sequence types (STs), plasmid types, virulence and resistance genes were identified. One ST340 clone, part of clonal complex 258 (CC258) and an ST that drives hospital outbreaks worldwide, harbouring numerous resistance genes and plasmids, was implicated as the likely cause of the outbreak. This study contributes molecular information necessary for tracking and characterizing this important hospital pathogen in sub-Saharan Africa.