Epistasis, core-genome disharmony, and adaptation in recombining bacteria.

Recombination of short DNA fragments via horizontal gene transfer (HGT) can introduce beneficial alleles, create genomic disharmony through negative epistasis, and create adaptive gene combinations through positive epistasis. For non-core (accessory) genes, the negative epistatic cost is likely to be minimal because the incoming genes have not co-evolved with the recipient genome and are frequently observed as tightly linked cassettes with major effects. By contrast, interspecific recombination in the core genome is expected to be rare because disruptive allelic replacement is likely to introduce negative epistasis. Why then is homologous recombination common in the core of bacterial genomes? To understand this enigma, we take advantage of an exceptional model system, the common enteric pathogens Campylobacter jejuni and C. coli that are known for very high magnitude interspecies gene flow in the core genome. As expected, HGT does indeed disrupt co-adapted allele pairings, indirect evidence of negative epistasis. However, multiple HGT events enable recovery of the genome's co-adaption between introgressing alleles, even in core metabolism genes (e.g., formate dehydrogenase). These findings demonstrate that, even for complex traits, genetic coalitions can be decoupled, transferred, and independently reinstated in a new genetic background-facilitating transition between fitness peaks. In this example, the two-step recombinational process is associated with C. coli that are adapted to the agricultural niche.IMPORTANCEGenetic exchange among bacteria shapes the microbial world. From the acquisition of antimicrobial resistance genes to fundamental questions about the nature of bacterial species, this powerful evolutionary force has preoccupied scientists for decades. However, the mixing of genes between species rests on a paradox: 0n one hand, promoting adaptation by conferring novel functionality; on the other, potentially introducing disharmonious gene combinations (negative epistasis) that will be selected against. Taking an interdisciplinary approach to analyze natural populations of the enteric bacteria Campylobacter, an ideal example of long-range admixture, we demonstrate that genes can independently transfer across species boundaries and rejoin in functional networks in a recipient genome. The positive impact of two-gene interactions appears to be adaptive by expanding metabolic capacity and facilitating niche shifts through interspecific hybridization. This challenges conventional ideas and highlights the possibility of multiple-step evolution of multi-gene traits by interspecific introgression.

Genomic resistant determinants of multidrug-resistant Campylobacter spp. isolates in Peru.

OBJECTIVES: Antimicrobial resistant (AMR) Campylobacter is a global health threat; however, there is limited information on genomic determinants of resistance in low- and middle-income countries. We evaluated genomic determinants of AMR using a collection of whole genome sequenced Campylobacter jejuni and C. coli isolates from Iquitos, Peru. METHODS: Campylobacter isolates from two paediatric cohort studies enriched with isolates that demonstrated resistance to ciprofloxacin and azithromycin were sequenced and mined for AMR determinants. RESULTS: The gyrA mutation leading to the Thr86Ile amino acid change was the only gyrA mutation associated with fluoroquinolone resistance identified. The A2075G mutation in 23S rRNA was present, but three other 23S rRNA mutations previously associated with macrolide resistance were not identified. A resistant-enhancing variant of the cmeABC efflux pump genotype (RE-cmeABC) was identified in 36.1% (35/97) of C. jejuni genomes and 17.9% (12/67) of C. coli genomes. Mutations identified in the CmeR-binding site, an inverted repeat sequence in the cmeABC promoter region that increases expression of the operon, were identified in 24/97 C. jejuni and 14/67 C. coli genomes. The presence of these variants, in addition to RE-cmeABC, was noted in 18 of the 24 C. jejuni and 9 of the 14 C. coli genomes. CONCLUSIONS: Both RE-cmeABC and mutations in the CmeR-binding site were strongly associated with the MDR phenotype in C. jejuni and C. coli. This is the first report of RE-cmeABC in Peru and suggests it is a major driver of resistance to the principal therapies used to treat human campylobacteriosis in this setting.

Recombination-mediated dissemination of Methicillin-resistant S. aureus clonal complex 1 in the Egyptian health care settings.

BACKGROUND: Methicillin-resistant Staphylococcus aureus (MRSA) is a rapidly evolving pathogen that is frequently associated with outbreaks and sustained epidemics. This study investigated the population structure, resistome, virulome, and the correlation between antimicrobial resistance determinants with phenotypic resistance profiles of 36 representative hospital-acquired MRSA isolates recovered from hospital settings in Egypt. RESULTS: The community-acquired MRSA lineage, clonal complex 1 (CC1) was the most frequently detected clone, followed by three other globally disseminated clones, CC121, CC8, and CC22. Most isolates carried SCCmec type V and more than half of isolates demonstrated multi-drug resistant phenotypes. Resistance to linezolid, a last resort antibiotic for treating multidrug resistant MRSA, was observed in 11.11% of the isolates belonging to different genetic backgrounds. Virulome analysis indicated that most isolates harboured a large pool of virulence factors and toxins. Genes encoding aureolysin, gamma hemolysins, and serine proteases were the most frequently detected virulence encoding genes. CC1 was observed to have a high pool of AMR resistance determinants including cfr, qacA, and qacB genes, which are involved in linezolid and quaternary ammonium compounds resistance, as well as high content of virulence-related genes, including both of the PVL toxin genes. Molecular clock analysis revealed that CC1 had the greatest frequency of recombination (compared to mutation) among the four major clones, supporting the role of horizontal gene transfer in modulating AMR and hypervirulence in this clone. CONCLUSIONS: This pilot study provided evidence on the dissemination success of CA-MRSA clone CC1 among Egyptian hospitals. Co-detection of multiple AMR and virulence genes in this lineage pose a broad public health risk, with implications for successful treatment. The results of this study, together with other surveillance studies in Egypt, should be used to develop strategies for controlling MRSA infections in Egyptian health-care settings.

Contrasting genes conferring short- and long-term biofilm adaptation in Listeria.

Listeria monocytogenes is an opportunistic food-borne bacterium that is capable of infecting humans with high rates of hospitalization and mortality. Natural populations are genotypically and phenotypically variable, with some lineages being responsible for most human infections. The success of L. monocytogenes is linked to its capacity to persist on food and in the environment. Biofilms are an important feature that allow these bacteria to persist and infect humans, so understanding the genetic basis of biofilm formation is key to understanding transmission. We sought to investigate the biofilm-forming ability of L. monocytogenes by identifying genetic variation that underlies biofilm formation in natural populations using genome-wide association studies (GWAS). Changes in gene expression of specific strains during biofilm formation were then investigated using RNA sequencing (RNA-seq). Genetic variation associated with enhanced biofilm formation was identified in 273 genes by GWAS and differential expression in 220 genes by RNA-seq. Statistical analyses show that the number of overlapping genes flagged by either type of experiment is less than expected by random sampling. This novel finding is consistent with an evolutionary scenario where rapid adaptation is driven by variation in gene expression of pioneer genes, and this is followed by slower adaptation driven by nucleotide changes within the core genome.

Evidence of potential Campylobacter jejuni zooanthroponosis in captive macaque populations.

Non-human primates share recent common ancestry with humans and exhibit comparable disease symptoms. Here, we explored the transmission potential of enteric bacterial pathogens in monkeys exhibiting symptoms of recurrent diarrhoea in a biomedical research facility in China. The common zoonotic bacterium Campylobacter jejuni was isolated from macaques (Macaca mulatta and Macaca fascicularis) and compared to isolates from humans and agricultural animals in Asia. Among the monkeys sampled, 5 % (44/973) tested positive for C. jejuni, 11 % (5/44) of which displayed diarrhoeal symptoms. Genomic analysis of monkey isolates, and 1254 genomes from various sources in Asia, were used to identify the most likely source of human infection. Monkey and human isolates shared high average nucleotide identity, common MLST clonal complexes and clustered together on a phylogeny. Furthermore, the profiles of putative antimicrobial resistance genes were similar between monkeys and humans. Taken together these findings suggest that housed macaques became infected with C. jejuni either directly from humans or via a common contamination source.

Loss of Pde1 function acts as an evolutionary gateway to penicillin resistance in Streptococcus pneumoniae.

Streptococcus pneumoniae is a major human pathogen and rising resistance to ß-lactam antibiotics, such as penicillin, is a significant threat to global public health. Mutations occurring in the penicillin-binding proteins (PBPs) can confer high-level penicillin resistance but other poorly understood genetic factors are also important. Here, we combined strictly controlled laboratory experiments and population analyses to identify a new penicillin resistance pathway that is independent of PBP modification. Initial laboratory selection experiments identified high-frequency pde1 mutations conferring S. pneumoniae penicillin resistance. The importance of variation at the pde1 locus was confirmed in natural and clinical populations in an analysis of >7,200 S. pneumoniae genomes. The pde1 mutations identified by these approaches reduce the hydrolytic activity of the Pde1 enzyme in bacterial cells and thereby elevate levels of cyclic-di-adenosine monophosphate and penicillin resistance. Our results reveal rapid de novo loss of function mutations in pde1 as an evolutionary gateway conferring low-level penicillin resistance. This relatively simple genomic change allows cells to persist in populations on an adaptive evolutionary pathway to acquire further genetic changes and high-level penicillin resistance.

Oral and gut microbial biomarkers of susceptibility to respiratory tract infection in adults: A feasibility study.

We conducted a feasibility cohort study which aimed to recruit and retain adults from the community to collect saliva (oral) and stool (gut) samples at three time points, at the start of the study (baseline), during a respiratory tract infection (RTI) and post-RTI. Community RTIs place a huge burden on health care services, and a non-invasive microbial diagnostic tool to predict the most vulnerable to respiratory infection would be ideal. To this aim, we analysed oral-gut baseline samples comparing those who reported RTI symptoms to those who remained healthy throughout the study for microbial biomarkers of respiratory susceptibility. Amplicon sequence variants (ASV) were identified by 16S sequence profiling to reveal oral-gut microbes. Reverse transcriptase-polymerase chain reaction (RT-PCR) was applied to target common respiratory microbes. Two general practices were recruited, and the participant recruitment rate was 1.3%. A total of 40 adult participants were retained, of which 19 acquired an RTI whereas 21 remained healthy. In healthy baseline oral and gut samples, ASVs from participants with RTI symptoms compared to those who remained healthy were similar with a high relative abundance of Streptococcus sp., and Blautia sp., respectively. Linear discriminant analysis effect size (LEfSe) revealed baseline oral microbes differed, indicating participants who suffered RTI symptoms had enhanced Streptococcus sobrinus and Megamonas sp., and depletion of Lactobacillus salivarius, Synergistetes, Verrucomicrobia and Dethiosulfovibrio. Furthermore, a random forest model ranked Streptococcus (4.13) as the highest mean decrease in accuracy (MDA) and RT-PCR showed a higher level of carriage of coagulase-negative Staphylococcus. Baseline core gut microbes were similar in both participant groups whereas LEfSe analysis revealed enhanced Veillonella, Rikenellaceae, Enhydobacteria, Eggerthella and Xanthomonsdales and depleted Desulfobulbus and Coprobacillus. Sutterella (4.73) had a high MDA value. Overall, we demonstrated the feasibility of recruiting and retaining adult participants from the community to provide multiple biological samples for microbial profiling. Our analyses identified potential oral-gut microbial biomarkers of respiratory infection susceptibility in otherwise healthy participants.

A high fidelity approach to assembling the complex Borrelia genome.

BACKGROUND: Bacteria of the Borrelia burgdorferi sensu lato (s.l.) complex can cause Lyme borreliosis. Different B. burgdorferi s.l. genospecies vary in their host and vector associations and human pathogenicity but the genetic basis for these adaptations is unresolved and requires completed and reliable genomes for comparative analyses. The de novo assembly of a complete Borrelia genome is challenging due to the high levels of complexity, represented by a high number of circular and linear plasmids that are dynamic, showing mosaic structure and sequence homology. Previous work demonstrated that even advanced approaches, such as a combination of short-read and long-read data, might lead to incomplete plasmid reconstruction. Here, using recently developed high-fidelity (HiFi) PacBio sequencing, we explored strategies to obtain gap-free, complete and high quality Borrelia genome assemblies. Optimizing genome assembly, quality control and refinement steps, we critically appraised existing techniques to create a workflow that lead to improved genome reconstruction. RESULTS: Despite the latest available technologies, stand-alone sequencing and assembly methods are insufficient for the generation of complete and high quality Borrelia genome assemblies. We developed a workflow pipeline for the de novo genome assembly for Borrelia using several types of sequence data and incorporating multiple assemblers to recover the complete genome including both circular and linear plasmid sequences. CONCLUSION: Our study demonstrates that, with HiFi data and an ensemble reconstruction pipeline with refinement steps, chromosomal and plasmid sequences can be fully resolved, even for complex genomes such as Borrelia. The presented pipeline may be of interest for the assembly of further complex microbial genomes.

Repeated out-of-Africa expansions of Helicobacter pylori driven by replacement of deleterious mutations.

Helicobacter pylori lives in the human stomach and has a population structure resembling that of its host. However, H. pylori from Europe and the Middle East trace substantially more ancestry from modern African populations than the humans that carry them. Here, we use a collection of Afro-Eurasian H. pylori genomes to show that this African ancestry is due to at least three distinct admixture events. H. pylori from East Asia, which have undergone little admixture, have accumulated many more non-synonymous mutations than African strains. European and Middle Eastern bacteria have elevated African ancestry at the sites of these mutations, implying selection to remove them during admixture. Simulations show that population fitness can be restored after bottlenecks by migration and subsequent admixture of small numbers of bacteria from non-bottlenecked populations. We conclude that recent spread of African DNA has been driven by deleterious mutations accumulated during the original out-of-Africa bottleneck.

Spatiotemporal Investigation of Antibiotic Resistance in the Urban Water Cycle Influenced by Environmental and Anthropogenic Activity.

With increasing emergence of antimicrobial resistant bacteria (ARB) and the risk this poses to public health, there are growing concerns regarding water pollution contributing to the spread of antimicrobial resistance (AMR) through inadequate amenities and the rapid rate of urbanization. In this study, the impact of different anthropogenic factors on the prevalence of AMR in the urban water cycle in Stellenbosch, South Africa (SA) was examined. Carbapenem, colistin, gentamicin and sulfamethoxazole resistant Gram-negative bacteria were recovered by selectively culturing aqueous, biofilm and sediment samples from sites impacted to varying degrees by informal settlements, residential, industrial, and agricultural activities, as well as a municipal wastewater treatment works (WWTW). A metagenomic approach determined community profiles and dominant AMR genes at various sites, while carbapenem resistant colonies were characterized using whole genome sequencing (WGS). Isolates recovered from agricultural sites exhibited relatively high levels of resistance to carbapenems and colistin, whereas sites impacted by domestic run-off had a higher prevalence of resistance to gentamicin and sulfamethoxazole, corresponding to usage data in SA. Similar microbial taxa were identified in raw sewage, sites downstream of informal settlements, and industrial areas that have limited waste removal infrastructure while WWTW were seen to reduce the prevalence of ARB in treated wastewater when operating efficiently. The results indicate the multiple complex drivers underpinning environmental dissemination of AMR and suggest that WWTW assist in removing AMR from the environment, reinforcing the necessity of adequate waste removal infrastructure and antibiotic stewardship measures to mitigate AMR transmission. IMPORTANCE The results from this study are of importance as they fill a gap in the data available on environmental AMR in South Africa to date. This study was done in parallel with co-investigators focusing on the prevalence of various antimicrobials at the same sites selected in our study, verifying that the sites that are influenced by informal settlements and WWTW influent had higher concentrations of antimicrobials and antimicrobial metabolites. The various locations of the sample sites selected, the frequency of the samples collected over a year, and the different types of samples collected at each site all contribute to informing how AMR in the environment might be affected by anthropogenic activity.

Assessment of community-wide antimicrobials usage in Eastern China using wastewater-based epidemiology.

Wastewater-based epidemiology (WBE) has potential to identify the epidemiological links between people, animals, and the environment, as part of antimicrobial resistance (AMR) surveillance. In this study, we investigated six wastewater treatment plants (WWTPs) serving six communities located in two regions in Eastern China: Site A in Zhejiang and site B in Jiangsu province to assess the public use of antimicrobial agents (AA). Fifty antimicrobials and 24 of their metabolites were quantified using ultraperformance liquid chromatography coupled with triple quadrupole tandem mass spectrometry (UPLC-MS/MS). Spatiotemporal trends were established for measured concentrations, daily loads, and population-normalised daily loads. Daily AA mass loads varied between 1.6 g/day and 324.6 g/day reflecting the WWTP scales, with macrolides and ß-lactams showing the highest overall environmental burden at 223.7 g/day and 173.7 g/day, respectively. Emissions of antibiotic residues from manufacturing have been observed, with the peak daily load 12-fold higher than the overall load from a community serving a population of over 600,000. Community exposure levels of 225.2 ± 156.2 mg/day/1000 inhabitant and 351.9 ± 133.5 mg/day/1000 inhabitant were recorded in site A and B, respectively. Paired parent-metabolites analysis identified a large proportion (64-78%) of un-metabolised metronidazole and clindamycin at site B, indicating improper disposal of unused drugs either in the community or in livestock production. Consumption levels, calculated via WBE, suggested relatively low antimicrobial usage in Eastern China compared to other areas in China. This first application of WBE in Eastern China to assess the community-wide exposure to AAs has potential to inform regional antimicrobial stewardship.

Impacts of Antibiotic Residues in the Environment on Bacterial Resistance and Human Health in Eastern China: An Interdisciplinary Mixed-Methods Study Protocol.

Antibiotic resistance is a global health challenge that threatens human and animal lives, especially among low-income and vulnerable populations in less-developed countries. Its multi-factorial nature requires integrated studies on antibiotics and resistant bacteria in humans, animals, and the environment. To achieve a comprehensive understanding of the situation and management of antibiotic use and environmental transmission, this paper describes a study protocol to document human exposure to antibiotics from major direct and indirect sources, and its potential health outcomes. Our mixed-methods approach addresses both microbiological and pathogen genomics, and epidemiological, geospatial, anthropological, and sociological aspects. Implemented in two rural residential areas in two provinces in Eastern China, linked sub-studies assess antibiotic exposure in population cohorts through household surveys, medicine diaries, and biological sampling; identify the types and frequencies of antibiotic resistance genes in humans and food-stock animals; quantify the presence of antibiotic residues and antibiotic resistance genes in the aquatic environment, including wastewater; investigate the drivers and behaviours associated with human and livestock antibiotic use; and analyse the national and local policy context, to propose strategies and systematic measurements for optimising and monitoring antibiotic use. As a multidisciplinary collaboration between institutions in the UK and China, this study will provide an in-depth understanding of the influencing factors and allow comprehensive awareness of the complexity of AMR and antibiotic use in rural Eastern China.

Signatures of selection in core and accessory genomes indicate different ecological drivers of diversification among Bacillus cereus clades.

Bacterial clades are often ecologically distinct, despite extensive horizontal gene transfer (HGT). How selection works on different parts of bacterial pan-genomes to drive and maintain the emergence of clades is unclear. Focusing on the three largest clades in the diverse and well-studied Bacillus cereus sensu lato group, we identified clade-specific core genes (present in all clade members) and then used clade-specific allelic diversity to identify genes under purifying and diversifying selection. Clade-specific accessory genes (present in a subset of strains within a clade) were characterized as being under selection using presence/absence in specific clades. Gene ontology analyses of genes under selection revealed that different gene functions were enriched in different clades. Furthermore, some gene functions were enriched only amongst clade-specific core or accessory genomes. Genes under purifying selection were often clade-specific, while genes under diversifying selection showed signs of frequent HGT. These patterns are consistent with different selection pressures acting on both the core and the accessory genomes of different clades and can lead to ecological divergence in both cases. Examining variation in allelic diversity allows us to uncover genes under clade-specific selection, allowing ready identification of strains and their ecological niche.

Stereoselective metabolism of chloramphenicol by bacteria isolated from wastewater, and the importance of stereochemistry in environmental risk assessments for antibiotics.

Wastewater treatment plants have been highlighted as a potential hotspot for the development and spread of antibiotic resistance. Although antibiotic resistant bacteria in wastewater present a public health threat, it is also possible that these bacteria play an important role in the bioremediation through the metabolism of antibiotics before they reach the wider environment. Here we address this possibility with a particular emphasis on stereochemistry using a combination of microbiology and analytical chemistry tools including the use of supercritical-fluid chromatography coupled with mass spectrometry for chiral analysis and high-resolution mass spectrometry to investigate metabolites. Due to the complexities around chiral analysis the antibiotic chloramphenicol was used as a proof of concept to demonstrate stereoselective metabolism due to its relatively simple chemical structure and availability over the counter in the U.K. The results presented here demonstrate the chloramphenicol can be stereoselectively transformed by the chloramphenicol acetyltransferase enzyme with the orientation around the first stereocentre being key for this process, meaning that accumulation of two isomers may occur within the environment with potential impacts on ecotoxicity and emergence of bacterial antibiotic resistance within the environment.

Host ecology regulates interspecies recombination in bacteria of the genus Campylobacter.

Horizontal gene transfer (HGT) can allow traits that have evolved in one bacterial species to transfer to another. This has potential to rapidly promote new adaptive trajectories such as zoonotic transfer or antimicrobial resistance. However, for this to occur requires gaps to align in barriers to recombination within a given time frame. Chief among these barriers is the physical separation of species with distinct ecologies in separate niches. Within the genus Campylobacter, there are species with divergent ecologies, from rarely isolated single-host specialists to multihost generalist species that are among the most common global causes of human bacterial gastroenteritis. Here, by characterizing these contrasting ecologies, we can quantify HGT among sympatric and allopatric species in natural populations. Analyzing recipient and donor population ancestry among genomes from 30 Campylobacter species, we show that cohabitation in the same host can lead to a six-fold increase in HGT between species. This accounts for up to 30% of all SNPs within a given species and identifies highly recombinogenic genes with functions including host adaptation and antimicrobial resistance. As described in some animal and plant species, ecological factors are a major evolutionary force for speciation in bacteria and changes to the host landscape can promote partial convergence of distinct species through HGT.

Strain wars and the evolution of opportunistic pathogens.

Bacteria live in complex communities with multiple species and strains competing with each other. Victories and defeats within these microbial wars are largely ignored unless they have a noticeable impact on the environment or the host, for example when a disease causing strain emerges as a winner. Evolutionary theory typically explains pathogen emergence as a trade-off between virulence and transmissibility. However, for opportunistic pathogens the secondary infection niche is often a dead end, as the host is either killed or cured, so a trade-off can not develop. In this context it is difficult to explain the maintenance of virulence genes in the population as they would be costly. Here, current literature is synthesized to address this apparent conundrum. The potential for adaptations to one niche to provide a benefit in another is described for some pathogenic species and this paradigm is extended to include genetic diversity and competition among individual strains. Finally, considering assemblages of strains in fluctuating immune environments with complex micro-niche structure, a scenario is presented in which commensal organisms can be primed for invasive disease should the opportunity arise.