ABSTRACT
Humans are radically altering global ecology, and one of the most apparent human-induced effects is urbanization, where high-density human habitats disrupt long-established ecotones. Changes to these transitional areas between organisms, especially enhanced contact among humans and wild animals, provide new opportunities for the spread of zoonotic pathogens. This poses a serious threat to global public health, but little is known about how habitat disruption impacts cross-species pathogen spread. Here, we investigated variation in the zoonotic enteric pathogen Campylobacter jejuni. The ubiquity of C. jejuni in wild bird gut microbiomes makes it an ideal organism for understanding how host behavior and ecology influence pathogen transition and spread. We analyzed 700 C. jejuni isolate genomes from 30 bird species in eight countries using a scalable generalized linear model approach. Comparing multiple behavioral and ecological traits showed that proximity to human habitation promotes lineage diversity and is associated with antimicrobial-resistant (AMR) strains in natural populations. Specifically, wild birds from urban areas harbored up to three times more C. jejuni genotypes and AMR genes. This study provides novel methodology and much-needed quantitative evidence linking urbanization to gene pool spread and zoonoses.
ABSTRACT
Population genomics has revolutionized our ability to study bacterial evolution by enabling data-driven discovery of the genetic architecture of trait variation. Genome-wide association studies (GWAS) have more recently become accompanied by genome-wide epistasis and co-selection (GWES) analysis, which offers a phenotype-free approach to generating hypotheses about selective processes that simultaneously impact multiple loci across the genome. However, existing GWES methods only consider associations between distant pairs of loci within the genome due to the strong impact of linkage-disequilibrium (LD) over short distances. Based on the general functional organisation of genomes it is nevertheless expected that majority of co-selection and epistasis will act within relatively short genomic proximity, on co-variation occurring within genes and their promoter regions, and within operons. Here, we introduce LDWeaver, which enables an exhaustive GWES across both short- and long-range LD, to disentangle likely neutral co-variation from selection. We demonstrate the ability of LDWeaver to efficiently generate hypotheses about co-selection using large genomic surveys of multiple major human bacterial pathogen species and validate several findings using functional annotation and phenotypic measurements. Our approach will facilitate the study of bacterial evolution in the light of rapidly expanding population genomic data.
ABSTRACT
Campylobacter is a leading cause of food-borne gastroenteritis worldwide, linked to the consumption of contaminated poultry meat. Targeting this pathogen at source, vaccines for poultry can provide short-term caecal reductions in Campylobacter numbers in the chicken intestine. However, this approach is unlikely to reduce Campylobacter in the food chain or human incidence. This is likely as vaccines typically target only a subset of the high genomic strain diversity circulating among chicken flocks, and rapid evolution diminishes vaccine efficacy over time. To address this, we used a genomic approach to develop a whole-cell autogenous vaccine targeting isolates harbouring genes linked to survival outside of the host. We hyper-immunised a whole major UK breeder farm to passively target offspring colonisation using maternally-derived antibody. Monitoring progeny, broiler flocks revealed a near-complete shift in the post-vaccination Campylobacter population with an ~50% reduction in isolates harbouring extra-intestinal survival genes and a significant reduction of Campylobacter cells surviving on the surface of meat. Based on these findings, we developed a logistic regression model that predicted that vaccine efficacy could be extended to target 65% of a population of clinically relevant strains. Immuno-manipulation of poultry microbiomes towards less harmful commensal isolates by competitive exclusion, has major potential for reducing pathogens in the food production chain.
ABSTRACT
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.
Subject(s)
Campylobacter coli , Campylobacter jejuni , Epistasis, Genetic , Gene Transfer, Horizontal , Genome, Bacterial , Recombination, Genetic , Campylobacter jejuni/genetics , Campylobacter coli/genetics , Evolution, Molecular , Adaptation, Physiological/genetics , Adaptation, Biological/geneticsABSTRACT
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.
Subject(s)
Anti-Bacterial Agents , Campylobacter , Humans , Child , Anti-Bacterial Agents/pharmacology , Peru , RNA, Ribosomal, 23S/genetics , Drug Resistance, Bacterial/genetics , Macrolides , Campylobacter/genetics , GenomicsABSTRACT
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.
Subject(s)
Methicillin-Resistant Staphylococcus aureus , Staphylococcal Infections , Humans , Methicillin Resistance/genetics , Egypt/epidemiology , Linezolid/pharmacology , Pilot Projects , Staphylococcal Infections/epidemiology , Staphylococcal Infections/drug therapy , Anti-Bacterial Agents/pharmacology , Anti-Bacterial Agents/therapeutic use , Clone Cells , Recombination, Genetic , Delivery of Health Care , Microbial Sensitivity TestsABSTRACT
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.
Subject(s)
Streptococcus pneumoniae , beta-Lactam Resistance , Humans , beta-Lactam Resistance/genetics , Penicillin-Binding Proteins/metabolism , Penicillin Resistance/genetics , Penicillins/pharmacology , Penicillins/metabolism , Bacterial Proteins/metabolism , Anti-Bacterial Agents/pharmacology , Microbial Sensitivity TestsABSTRACT
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.
Subject(s)
Listeria monocytogenes , Listeria , Humans , Listeria/genetics , Genome-Wide Association Study , Biofilms , Listeria monocytogenes/geneticsABSTRACT
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.
Subject(s)
Biomedical Research , Campylobacter jejuni , Animals , Macaca , Campylobacter jejuni/genetics , Multilocus Sequence Typing , Asia , Diarrhea/veterinaryABSTRACT
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.
ABSTRACT
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.
Subject(s)
Borrelia burgdorferi Group , Borrelia burgdorferi , Borrelia , Lyme Disease , Humans , Borrelia/genetics , Genome, Bacterial , Phylogeny , Borrelia burgdorferi/genetics , Lyme Disease/microbiology , Borrelia burgdorferi Group/geneticsABSTRACT
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.
Subject(s)
Helicobacter Infections , Helicobacter pylori , Humans , Helicobacter pylori/genetics , Helicobacter Infections/microbiology , Black People/genetics , Africa , MutationABSTRACT
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.
Subject(s)
Anti-Infective Agents , Drug Resistance, Bacterial , Wastewater , Sewage , Water Cycle , Colistin , Angiotensin Receptor Antagonists , Anthropogenic Effects , Angiotensin-Converting Enzyme Inhibitors , Anti-Bacterial Agents/pharmacology , Carbapenems , Anti-Infective Agents/pharmacology , Gentamicins , SulfamethoxazoleABSTRACT
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.
Subject(s)
Anti-Infective Agents , Water Pollutants, Chemical , Animals , Anti-Bacterial Agents , China , Chromatography, Liquid , Humans , Tandem Mass Spectrometry/methods , Wastewater/chemistry , Wastewater-Based Epidemiological Monitoring , Water Pollutants, Chemical/analysisABSTRACT
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.
Subject(s)
Anti-Bacterial Agents , Wastewater , Animals , Anti-Bacterial Agents/pharmacology , Anti-Bacterial Agents/therapeutic use , Bacteria/genetics , China , Drug Resistance, Bacterial/genetics , Drug Resistance, Microbial/genetics , Humans , Wastewater/microbiologyABSTRACT
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.