Significance Statement Genetic 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. On 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 analyse 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 epistasis in a recipient genome. This challenges conventional ideas and highlights the possibility of single step evolution by saltation. Abstract Recombination of short DNA fragments via horizontal gene transfer (HGT) can both introduce beneficial alleles and create genomic disharmony through negative 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. By contrast, for the core genome, interspecific recombination 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 Campylobacter 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 (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.

Abstract Lumpy skin disease virus (LSDV) causes severe disease in cattle and water buffalo and is transmitted by hematophagous arthropod vectors. Detailed information of the adaptive and innate immune response to LSDV is limited, hampering the development of tools to control the disease. This study provides an in-depth analysis of the immune responses of calves experimentally inoculated with LSDV via either needle-inoculation or arthropod-inoculation using virus-positive Stomoxys calcitrans and Aedes aegypti vectors. Seven out of seventeen needle-inoculated calves (41%) developed clinical disease characterised by multifocal necrotic cutaneous nodules. In comparison 8/10 (80%) of the arthropod-inoculated calves developed clinical disease. A variable LSDV-specific IFN-γ immune response was detected in the needle-inoculated calves from 5 days post inoculation (dpi) onwards, with no difference between clinical calves (developed cutaneous lesions) and nonclinical calves (did not develop cutaneous lesions). In contrast a robust and uniform cell-mediated immune response was detected in all eight clinical arthropod-inoculated calves, with little response detected in the two nonclinical arthropod-inoculated calves. Neutralising antibodies against LSDV were detected in all inoculated cattle from 5-7 dpi. Comparison of the production of anti-LSDV IgM and IgG antibodies revealed no difference between clinical and nonclinical needle-inoculated calves, however a strong IgM response was evident in the nonclinical arthropod-inoculated calves but absent in the clinical arthropod-inoculated calves. This suggests that early IgM production is a correlate of protection in LSD. This study presents the first evidence of differences in the immune response between clinical and nonclinical cattle and highlights the importance of using a relevant transmission model when studying LSD.

ABSTRACT Flagellar motility plays a central role in the bacterial foodborne pathogen Campylobacter jejuni , as flagellar motility is required for reaching the intestinal epithelium and subsequent colonisation or disease. Flagellar proteins also contribute strongly to biofilm formation during transmission. Chemotaxis is the process directing flagellar motility in response to attractant and repellent stimuli, but its role in biofilm formation of C. jejuni is not well understood. Here we show that inactivation of the core chemotaxis genes cheVAWY in C. jejuni strain NCTC 11168 affects both chemotactic motility and biofilm formation. Inactivation of any of the core chemotaxis genes ( cheA , cheY , cheV or cheW ) impaired chemotactic motility but did not affect flagellar assembly or growth. The Δ cheY mutant swam in clockwise loops, while complementation restored normal motility. Inactivation of the core chemotaxis genes interfered with the ability to form a discrete biofilm at the air-media interface, and the Δ cheY mutant displayed reduced dispersal/shedding of bacteria into the planktonic fraction. This suggests that while the chemotaxis system is not required for biofilm formation per se , it is necessary for organized biofilm formation. Hence interference with the Campylobacter chemotaxis system at any level disrupts optimal chemotactic motility and transmission modes such as biofilm formation.

Abstract Ascariasis is the most prevalent helminthic disease affecting both humans and pigs and is caused by the roundworms Ascaris lumbricoides and Ascaris suum . While preventive chemotherapy continues to be the most common control method, recent reports of anthelminthic resistance highlight the need for development of a vaccine against ascariasis. The aim of this study was to use a reverse vaccinology approach to identify potential vaccine candidates for Ascaris . Three Ascaris proteomes predicted from whole-genome sequences were analysed. Candidate proteins were identified using open-access bioinformatic tools (e.g. Vacceed, VaxiJen, Bepipred 2.0) which test for different characteristics such as sub-cellular location, T-cell and B-cell molecular binding, antigenicity, allergenicity and phylogenetic relationship with other nematode proteins. From over 100,000 protein sequences analysed, four transmembrane proteins were predicted to be non-allergen antigens and potential vaccine candidates. The four proteins are a Piezo protein, two voltage-dependent calcium channels and a protocadherin-like protein, are all expressed in either the muscle or ovaries of both Ascaris species, and all contained high affinity epitopes for T-cells and B-cells. The use of a reverse vaccinology approach allowed the prediction of four new potential vaccination targets against ascariasis in humans and pigs. These targets can now be further tested in in vitro and in vivo assays to prove efficacy in both pigs and humans.

Abstract Ascaris species are soil-transmitted helminths that infect humans and livestock mainly in low and middle-income countries. Benzimidazole (BZ) class drugs have predominated for many years in the treatment of Ascaris infections, but persistent use of BZs has already led to widespread resistance in other nematodes, and treatment failure is emerging for Ascaris . Benzimidazoles act by binding to β-tubulin proteins and destabilising microtubules. Three mutations in the β-tubulin protein family are associated with BZ resistance. Seven shared β-tubulin isotypes were identified in Ascaris lumbricoides and A. suum genomes. Benzimidazoles were predicted to bind to all β-tubulin isotypes using in silico docking, demonstrating that the selectivity of BZs to interact with one or two β-tubulin isotypes is likely the result of isotype expression levels affecting the frequency of interaction. Ascaris β-tubulin isotype A clusters with helminth β-tubulins previously shown to interact with BZ. Molecular dynamics simulations using β-tubulin isotype A highlighted the key role of amino acid E198 in BZ-β-tubulin interactions. Simulations indicated that mutations at amino acids E198A and F200Y alter binding of BZ, whereas there was no obvious effect of the F167Y mutation. In conclusion, the key interactions vital for BZ binding with β-tubulins have been identified and show how mutations can lead to resistance in nematodes.

Staphylococcus pseudintermedius is the foremost cause of opportunistic canine skin and mucosal infections worldwide. Multidrug resistant (MDR) and methicillin-resistant S. pseudintermedius (MRSP) lineages have disseminated globally in the last decade and present significant treatment challenges. However, little is known regarding the factors that contribute to the success of MDR lineages. In this study, we compared the genome sequence of 110 UK isolates of S. pseudintermedius to 2,166 genomes of S. pseudintermedius populations from different continents. A novel core genome multi-locus typing scheme was generated to allow large scale, rapid and detailed analysis of S. pseudintermedius phylogenies, and was used to show that the S. pseudintermedius population structure is broadly segregated into an MDR population and a non-MDR population. MRSP lineages are predicted to either encode certain resistance genes chromosomally, or on plasmids, and this is associated with their MLST sequence type. Comparison of lineages most frequently implicated in disease, ST-45 and ST-71, with the phylogenetically related ST-496 lineage that has a comparatively low disease rate, revealed that ST-45 and ST-71 genomes encode distinct combinations of phage-defence systems and concurrently encode a high number of intact prophages. In contrast, ST-496 genomes encode a wider array of phage defence systems and lack intact, complete prophages. Additionally, we show that distinct prophages are widespread in S. pseudintermedius and appear to account for the vast majority of genomic diversity. These findings indicate that MRSP lineages have significant structural genomic differences, and that prophage integration, and differential antiviral systems correlate with the emergence of successful genotypes.

ABSTRACT Campylobacter jejuni and Campylobacter coli are important bacterial causes of human foodborne illness. Despite several years of reduced antibiotics usage in livestock production in the UK and US, high prevalence of antimicrobial resistance (AMR) persists in Campylobacter . Both countries have instigated genome sequencing-based surveillance programs for Campylobacter , and here we have identified AMR genes in 32,256 C. jejuni and 8,776 C. coli publicly available genome sequences to compare the prevalence and trends of AMR in Campylobacter isolated in the UK and US between 2001-2018. AMR markers were detected in 68% of C. coli and 53% of C. jejuni , with 15% of C. coli being multi-drug resistant (MDR) compared to only 2% of C. jejuni . The prevalence of aminoglycoside, macrolide, quinolone and tetracycline resistance remained fairly stable from 2001-2018 in both C. jejuni and C. coli , but statistically significant differences were observed between the UK and US. There was a statistically significant higher prevalence of aminoglycoside and tetracycline resistance for US C. coli and C. jejuni , and macrolide resistance for US C. coli . In contrast, UK C. coli and C. jejuni showed a significantly higher prevalence of quinolone resistance. Specific MLST clonal complexes (e.g. ST-353/464) showed >95% quinolone resistance. This large-scale comparison of AMR prevalence has shown that the prevalence of AMR remains stable for Campylobacter in the UK and the US. This suggests that antimicrobial stewardship and restricted antibiotic usage may help contain further expansion of AMR prevalence in Campylobacter , but are unlikely to reduce it in the short term.

ABSTRACT The Type VIIb protein secretion system (T7SSb) is found in Bacillota (firmicute) bacteria and has been shown to mediate interbacterial competition. EssC is a membrane-bound ATPase that is a critical component of the T7SSb and plays a key role in substrate recognition. Prior analysis of available genome sequences of the foodborne bacterial pathogen Listeria monocytogenes has shown that although the T7SSb was encoded as part of the core genome, EssC could be found as one of seven different sequence variants. While each sequence variant was associated with a specific suite of candidate substrate proteins encoded immediately downstream of essC , many LXG-domain proteins were encoded across multiple essC sequence variants. Here we have extended this analysis using a diverse collection of 37,930 L. monocytogenes genomes. We have identified a rare eighth variant of EssC present in ten L. monocytogenes Lineage III genomes. These genomes also encode a large toxin of the rearrangement hotspot (Rhs) repeat family adjacent to essC8 , along with a probable immunity protein and three small accessory proteins. We have further identified nine novel LXG-domain proteins, and four additional chromosomal hotspots across L. monocytogenes genomes where LXG proteins can be encoded. The eight L. monocytogenes EssC variants were also found in other Listeria species, with additional novel EssC types also identified. Across the genus, species frequently encoded multiple EssC types, indicating that T7SSb diversity is a primary feature of the genus Listeria . DATA SUMMARY All genome sequences used in this study are available via Genbank, and the assembly accession numbers are provided in Table S1. This file also lists relevant metadata (name, source category, country, year and clonal complex). IMPACT STATEMENT Listeria monocytogenes is a soil-borne saprophytic bacterium and a food-borne pathogen of humans. Decomposing plant matter and the human GI tract are rich in diverse microbial species and to colonise these niches L. monocytogenes must be able to compete with other bacteria. The type VII secretion system (T7SS) of Bacillota has been shown to secrete protein toxins that target other bacteria. In this study we have analysed a diverse collection of L. monocytogenes genome sequences to study the diversity of the Listeria T7SS and its putative effector proteins. We show that the EssC component of the L. monocytogenes T7SS is highly diverse, clustering into one of eight sequence variants. Each EssC variant is associated with a specific toxin candidate, and the EssC8 variant T7SS likely secretes a novel rearrangement hotspot (Rhs) repeat toxin. We also identify multiple new LXG-families of T7SS toxins and describe genomic hotspots where they are encoded. We find no link between EssC variants and clinical outcome. In agreement with this, analysis of EssC variability in available genomes of other Listeria species showed that all eight L. monocytogenes EssC variants are present in non-monocytogenes Listeria species.

ABSTRACT Microbial genomes are highly adaptable, with mobile genetic elements (MGEs) such as integrative conjugative elements (ICE) mediating the dissemination of new genetic information throughout bacterial populations. This is countered by defence mechanism such as CRISPR-Cas systems, which limit invading MGEs by sequence-specific targeting. Here we report the distribution of the pVir, pTet and PCC42 plasmids and a new 70-129 kb ICE (CampyICE1) in the foodborne bacterial pathogens Campylobacter jejuni and Campylobacter coli . CampyICE1 contains a degenerated Type II-C CRISPR system consisting of a sole Cas9 protein, which is distinct from the previously described Cas9 proteins from C. jejuni and C. coli . CampyICE1 is conserved in structure and gene order, containing blocks of genes predicted to be involved in recombination, regulation, and conjugation. CampyICE1 was detected in 134/5,829 (2.3%) C. jejuni genomes and 92/1,347 (6.8%) C. coli genomes. Similar ICE were detected in a number of non-jejuni/coli Campylobacter species, although these lacked a CRISPR-Cas system. CampyICE1 carries 3 separate short CRISPR spacer arrays containing a combination of 108 unique spacers and 16 spacer variant families. A total of 69 spacers and 10 spacer variant families (63.7%) were predicted to target Campylobacter plasmids. The presence of a functional CampyICE1 Cas9 protein and matching anti-plasmid spacers was associated with the absence of the pVir, pTet and pCC42 plasmids (188/214 genomes, 87.9%), implicating that the CampyICE1-encoded CRISPR-Cas has contributed to the exclusion of competing plasmids. In conclusion, the characteristics of the CRISPR-Cas9 system on CampyICE1 suggests a history of plasmid warfare in Campylobacter . IMPACT STATEMENT Understanding pathogen evolution is paramount for enhancing food safety and limiting pathogenic disease in humans and animals. Campylobacter species comprise a group of human and animal pathogens with a remarkable success rate, being the most frequent cause of bacterial food-borne disease in high-income countries. A common theme among Campylobacter evolution is genomic plasticity, and a significant proportion of this plasticity is driven by horizontal gene transfer (HGT) that results in acquisition of complex traits in one evolutionary event. Understanding the mechanisms of transfer of MGEs and how MGEs such as integrative conjugative elements (ICE) exclude other MGEs is fundamental to understanding Campylobacter evolution. CRISPR-Cas9 proteins play a role in bacterial immune systems, mediating the defence against bacteriophage, plasmids, and integrative elements. The use of CRISPR-Cas by a mobile element to fight off competing elements, possibly to the advantage or detriment to their host, also increases our understanding of how important selfish genomic islands undergo co-evolution with bacterial pathogens, and generates insight into the complex warfare between MGEs. DATA STATEMENT All genome sequences used in this study are available on the National Center for Biotechnology Information (NCBI) Genome database or in the Campylobacter PubMLST website; the assembly accession numbers (NCBI Genome) or genome ID numbers (Campylobacter PubMLST) are listed in Table S1 (available in the online version of this article). Genome assemblies were quality checked based on N 50 , L 50 , genome size and number of contigs. CRISPR Spacer sequences and predicted targets, Cas9 alignments, presence of mobile elements and plasmids are all included in the Supplementary Information.