K-mer-based High-throughput Analysis of the Adaptive Potential of Campylobacter

Several species of the genus Campylobacter (C.) are zoonotic pathogens, especially C. jejuni and C. coli, the leading causes of foodborne diseases worldwide. Although both species colonize many hosts including poultry, livestock and wild animals, persistence mechanisms enabling the bacteria to adapt towards new ecological niches are not yet fully understood. In this work, novel k-mer-based methods enabling high-throughput analysis of whole-genome sequencing (WGS) data of C. jejuni and C. coli have been developed, extended and applied to investigate the adaptive potential of the distinct species towards different ecological niches and changing environments. A k-mer-based microbial genome-wide association study (GWAS) was set up to identify host-specific genomic signatures of C. jejuni isolated from chicken, cattle, pig and clinical human samples. GWAS revealed a strong association of both, the core and the accessory genome of C. jejuni, with distinct host animal species. Moreover, multiple adaptive trajectories defining the evolution of C. jejuni lifestyle preferences in different ecosystems were identified. In a second approach, WGS data of Campylobacter isolates that showed ambiguous probing results using different polymerase chain reaction (PCR)-based species classification methods during routine-diagnostics were investigated. The Campylobacter genomes were analyzed with respect to their genomic make-up. For this purpose, a k-mer-based method was developed in order to identify recombination events between C. jejuni and C. coli. The identified genes encode proteins that were commonly associated with important pathways involved in chromosome maintenance and DNA repair, membrane transport and stress defense. Overall, the results presented in this work promote molecular surveillance and rapid diagnostics of Campylobacter. In addition, host-specific allelic variants identified among different phylogenetic backgrounds might serve as important maker genes in future source attribution models for fast and precise retrograde outbreak investigation along the food chain.Verschiedene Spezies der Gattung Campylobacter (C.) sind zoonotische Krankheitserreger, die zu den Hauptverursachern von durch Lebensmittel übertragbare Infektionskrankheiten weltweit gehören. Obwohl C. jejuni und C. coli unterschiedliche Wirte wie Geflügel, Vieh und Wildtiere kolonisieren, sind die Mechanismen, die es diesen Bakterien ermöglichen, sich an neue ökologische Nischen anzupassen, nicht vollständig geklärt. In dieser Arbeit wurden neue k-mer-basierte Methoden für Hochdurchsatzanalysen von anzgenomsequenzierungen von C. jejuni und C. coli entwickelt, erweitert und angewendet, um das Anpassungspotenzial der Spezies an unterschiedliche Wirte und Umgebungen zu untersuchen. In der ersten Studie wurde eine auf k-meren basierende mikrobielle genomweite Assoziationsstudie (GWAS) durchgeführt, um wirtsspezifische genomische C. jejuni-Signaturen von Isolaten aus Hühnern, Rindern, Schweinen und humanen klinischen Proben zu identifizieren. Die GWAS zeigte eine starke Assoziation sowohl des Kern- wie auch des akzessorischen C. jejuni Genoms mit verschiedenen Wirtstieren. Durch die in silico Prädiktion von Veränderungen in Peptiden bzw. Proteinen ist es gelungen, mehrere adaptive metabolische Pfade zu identifiziert, welche potentiell die Evolution der Wirtspräferenz von phylogenetisch unterschiedlichen C. jejuni an verschiedene Lebensräume ermöglichen. In einem zweiten Ansatz wurden Ganzgenomsequenzen von Campylobacter Isolaten analysiert, die in der Routinediagnostik mittels Polymerase-Kettenreaktionen (PCR) nicht eindeutig einer genauen Spezies zuzuordnen waren. Die Campylobacter Genome aus diesen Proben wurden hinsichtlich ihres genomischen Aufbaus untersucht. Zu diesem Zweck wurde eine k-mer-basierte Methode entwickelt, um Rekombinationsereignisse zwischen C. jejuni und C. coli zu identifizieren, welche maßgeblich die Ergebnisse der PCR beeinflussten. Die auf diese Weise identifizierten Gene kodieren häufig Proteine mit wichtiger Funktion in der Chromosomenerhaltung bzw. DNA Reparatur, im Membrantransport und Stressabwehr. Die in dieser Arbeit vorgestellten Ergebnisse leisten einen Beitrag zur routinemäßigen Überwachung und schnellen Diagnostik von Campylobacter Ausbrüchen im Sinne einer integrierten molekularen Surveillance. Wirtsspezifische Allele, die in Campylobacter mit unterschiedlichen phylogenetischen Hintergründen identifiziert wurden, können dabei als wichtige Markergene dienen, um die ursprüngliche Quelle des Ausbruchs schnell und präzise retrograd entlang der Lebensmittelkette zu identifizieren

Transmission pathways of campylobacter spp. at broiler farms and their environment in Brandenburg, Germany

Broiler meat is widely known as an important source of foodborne Campylobacter jejuni and Campylobacter coli infections in humans. In this study, we thoroughly investigated transmission pathways that may contribute to possible Campylobacter contamination inside and outside broiler houses. For this purpose we carried out a comprehensive longitudinal sampling approach, using a semi-quantitative cultivation method to identify and quantify transmissions and reservoirs of Campylobacter spp.. Three german broiler farms in Brandenburg and their surrounding areas were intensively sampled, from April 2018 until September 2020. Consecutive fattening cycles and intervening downtimes after cleaning and disinfection were systematically sampled in summer and winter. To display the potential phylogeny of barn and environmental isolates, whole genome sequencing (WGS) and bioinformatic analyses were performed. Results obtained in this study showed very high Campylobacter prevalence in 51/76 pooled feces (67.1%) and 49/76 boot swabs (64.5%). Average counts between 6.4 to 8.36 log10MPN/g were detected in pooled feces. In addition, levels of 4.7 and 4.1 log10MPN/g were detected in boot swabs and litter, respectively. Samples from the barn interior showed mean Campyloacter values in swabs from drinkers 2.6 log10MPN/g, walls 2.0 log10MPN/g, troughs 1.7 log10MPN/g, boards 1.6 log10MPN/g, ventilations 0.9 log10MPN/g and 0.7 log10MPN/g for air samples. However, Campylobacter was detected only in 7/456 (1.5%) of the environmental samples (water bodies, puddles or water-filled wheel tracks; average of 0.6 log10MPN/g). Furthermore, WGS showed recurring Campylobacter genotypes over several consecutive fattening periods, indicating that Campylobacter genotypes persist in the environment during downtime periods. However, after cleaning and disinfection of the barns, we were unable to identify potential sources in the broiler houses. Interestingly, alternating Campylobacter genotypes were observed after each fattening period, also indicating sources of contamination from the wider environment outside the farm. Therefore, the results of this study suggest that a potential risk of Campylobacter transmission may originate from present environmental sources (litter and water reservoirs). However, the sources of Campylobacter transmission may vary depending on the operation and farm environmental conditions

Effects of a Four-Week High-Dosage Zinc Oxide Supplemented Diet on Commensal Escherichia coli of Weaned Pigs

Strategies to reduce economic losses associated with post-weaning diarrhea in pig farming include high-level dietary zinc oxide supplementation. However, excessive usage of zinc oxide in the pig production sector was found to be associated with accumulation of multidrug resistant bacteria in these animals, presenting an environmental burden through contaminated manure. Here we report on zinc tolerance among a random selection of intestinal Escherichia coli comprising of different antibiotic resistance phenotypes and sampling sites isolated during a controlled feeding trial from 16 weaned piglets: In total, 179 isolates from “pigs fed with high zinc concentrations” (high zinc group, [HZG]: n = 99) and a corresponding “control group” ([CG]: n = 80) were investigated with regard to zinc tolerance, antimicrobial- and biocide susceptibilities by determining minimum inhibitory concentrations (MICs). In addition, in silico whole genome screening (WGSc) for antibiotic resistance genes (ARGs) as well as biocide- and heavy metal tolerance genes was performed using an in-house BLAST-based pipeline. Overall, porcine E. coli isolates showed three different ZnCl2 MICs: 128 μg/ml (HZG, 2%; CG, 6%), 256 μg/ml (HZG, 64%; CG, 91%) and 512 μg/ml ZnCl2 (HZG, 34%, CG, 3%), a unimodal distribution most likely reflecting natural differences in zinc tolerance associated with different genetic lineages. However, a selective impact of the zinc-rich supplemented diet seems to be reasonable, since the linear mixed regression model revealed a statistically significant association between “higher” ZnCl2 MICs and isolates representing the HZG as well as “lower ZnCl2 MICs” with isolates of the CG (p = 0.005). None of the zinc chloride MICs was associated with a particular antibiotic-, heavy metal- or biocide- tolerance/resistance phenotype. Isolates expressing the 512 μg/ml MIC were either positive for ARGs conferring resistance to aminoglycosides, tetracycline and sulfamethoxazole-trimethoprim, or harbored no ARGs at all. Moreover, WGSc revealed a ubiquitous presence of zinc homeostasis and – detoxification genes, including zitB, zntA, and pit. In conclusion, we provide evidence that zinc-rich supplementation of pig feed selects for more zinc tolerant E. coli, including isolates harboring ARGs and biocide- and heavy metal tolerance genes – a putative selective advantage considering substances and antibiotics currently used in industrial pork production systems

High-Zinc Supplementation of Weaned Piglets Affects Frequencies of Virulence and Bacteriocin Associated Genes Among Intestinal Escherichia coli Populations

To prevent economic losses due to post-weaning diarrhea (PWD) in industrial pig production, zinc (Zn) feed additives have been widely used, especially since awareness has risen that the regular application of antibiotics promotes buildup of antimicrobial resistance in both commensal and pathogenic bacteria. In a previous study on 179 Escherichia coli collected from piglets sacrificed at the end of a Zn feeding trial, including isolates obtained from animals of a high-zinc fed group (HZG) and a corresponding control group (CG), we found that the isolate collection exhibited three different levels of tolerance toward zinc, i.e., the minimal inhibitory concentration (MIC) detected was 128, followed by 256 and 512 mu g/ml ZnCl2. We further provided evidence that enhanced zinc tolerance in porcine intestinal E. coli populations is clearly linked to excessive zinc feeding. Here we provide insights about the genomic make-up and phylogenetic background of these 179 E. coli genomes. Bayesian analysis of the population structure (BAPS) revealed a lack of association between the actual zinc tolerance level and a particular phylogenetic E. coli cluster or even branch for both, isolates belonging to the HZG and CG. In addition, detection rates for genes and operons associated with virulence (VAG) and bacteriocins (BAG) were lower in isolates originating from the HZG (41 vs. 65% and 22 vs. 35%, p < 0.001 and p = 0.002, resp.). Strikingly, E. coli harboring genes defining distinct pathotypes associated with intestinal disease, i.e., enterotoxigenic, enteropathogenic, and Shiga toxin-producing E. coli (ETEC, EPEC, and STEC) constituted 1% of the isolates belonging to the HZG but 14% of those from the CG. Notably, these pathotypes were positively associated with enhanced zinc tolerance (512 mu g/ml ZnCl2 MIC, p < 0.001). Taken together, zinc excess seems to influence carriage rates of VAGs and BAGs in porcine intestinal E. coli populations, and high-zinc feeding is negatively correlated with enteral pathotype occurrences, which might explain earlier observations concerning the relative increase of Enterobacterales considering the overall intestinal microbiota of piglets during zinc feeding trials while PWD rates have decreased

Frequency, Local Dynamics, and Genomic Characteristics of ESBL-Producing Escherichia coli Isolated From Specimens of Hospitalized Horses

Previous research identified veterinary clinics as hotspots with respect to accumulation and spread of multidrug resistant extended-spectrum β-lactamase (ESBL)-producing Escherichia coli (EC). Therefore, promoting the prudent use of antibiotics to decrease selective pressure in that particular clinical environment is preferable to enhance biosecurity for animal patients and hospital staff. Accordingly, this study comparatively investigated the impact of two distinct perioperative antibiotic prophylaxis (PAP) regimens (short-term versus prolonged) on ESBL-EC carriage of horses subjected to colic surgery. While all horses received a combination of penicillin/gentamicin (P/G) as PAP, they were assigned to either the “single-shot group” (SSG) or the conventional “5-day group” (5DG). Fecal samples collected on arrival (t0), on the 3rd (t1) and on the 10th day after surgery (t2) were screened for ESBL-EC. All isolates were further investigated using whole genome sequences. In total, 81 of 98 horses met the inclusion criteria for this study. ESBL-EC identified in samples available at t0, t1 and t2 were 4.8% (SSG) and 9.7% (5DG), 37% (SSG) and 47.2% (5DG) as well as 55.6% (SSG) and 56.8% (5DG), respectively. Regardless of the P/G PAP regimen, horses were 9.12 times (95% CI 2.79–29.7) more likely to carry ESBL-EC at t1 compared to t0 (p < 0.001) and 15.64 times (95% CI 4.57–53.55) more likely to carry ESBL-EC at t2 compared to t0 (p < 0.001). ESBL-EC belonging to sequence type (ST) 10, ST86, ST641, and ST410 were the most prevalent lineages, with blaCTX–M–1 (60%) being the dominant ESBL gene. A close spatio-temporal relationship between isolates sharing a particular ST was revealed by genome analysis, strongly indicating local spread. Consequently, hospitalization itself has a strong impact on ESBL-EC isolation rates in horses, possibly masking differences between distinct PAP regimens. The results of this study reveal accumulation and spread of multi-drug resistant ESBL-EC among horses subjected to colic surgery with different P/G PAP regimens, challenging the local hygiene management system and work-place safety of veterinary staff. Moreover, the predominance of particular ESBL-EC lineages in clinics providing health care for horses needs further investigation.Peer Reviewe

Genome‑wide insights into population structure and host specifcity of Campylobacter jejuni

The zoonotic pathogen Campylobacter jejuni is among the leading causes of foodborne diseases worldwide. While C. jejuni colonises many wild animals and livestock, persistence mechanisms enabling the bacterium to adapt to host species' guts are not fully understood. In order to identify putative determinants influencing host preferences of distinct lineages, bootstrapping based on stratified random sampling combined with a k-mer-based genome-wide association was conducted on 490 genomes from diverse origins in Germany and Canada. We show a strong association of both the core and the accessory genome characteristics with distinct host animal species, indicating multiple adaptive trajectories defining the evolution of C. jejuni lifestyle preferences in different ecosystems. Here, we demonstrate that adaptation towards a specific host niche ecology is most likely a long evolutionary and multifactorial process, expressed by gene absence or presence and allele variations of core genes. Several host-specific allelic variants from different phylogenetic backgrounds, including dnaE, rpoB, ftsX or pycB play important roles for genome maintenance and metabolic pathways. Thus, variants of genes important for C. jejuni to cope with specific ecological niches or hosts may be useful markers for both surveillance and future pathogen intervention strategies.Peer Reviewe

Silence as a way of niche adaptation: mecC-MRSA with variations in the accessory gene regulator (agr) functionality express kaleidoscopic phenotypes

Functionality of the accessory gene regulator (agr) quorum sensing system is an important factor promoting either acute or chronic infections by the notorious opportunistic human and veterinary pathogen Staphylococcus aureus. Spontaneous alterations of the agr system are known to frequently occur in human healthcare-associated S. aureus lineages. However, data on agr integrity and function are sparse regarding other major clonal lineages. Here we report on the agr system functionality and activity level in mecC-carrying methicillin resistant S. aureus (MRSA) of various animal origins (n = 33) obtained in Europe as well as in closely related human isolates (n = 12). Whole genome analysis assigned all isolates to four clonal complexes (CC) with distinct agr types (CC599 agr I, CC49 agr II, CC130 agr III and CC1943 agr IV). Agr functionality was assessed by a combination of phenotypic assays and proteome analysis. In each CC, isolates with varying agr activity levels were detected, including the presence of completely non-functional variants. Genomic comparison of the agr I–IV encoding regions associated these phenotypic differences with variations in the agrA and agrC genes. The genomic changes were detected independently in divergent lineages, suggesting that agr variation might foster viability and adaptation of emerging MRSA lineages to distinct ecological niches.Peer Reviewe

Masterthesis Appendix A3

Appendix table A3 for the Masterthesis : "Development of a K-mer based Method for High-Throughput Serotyping of S. pneumoniae".<br

simulated Reads

<i>Streptococcus pneumoniae</i> is responsible for nearly 1 million deaths every year worldwide. To target the most dominant serotypes with vaccines it is import to track the distribution and evolution of serotypes within the population. Recent efforts have been made to infer serotypes directly from genomic data. Current software approaches are limited and do not scale well. We introduce a novel method SeroBA which uses a hybrid assembly and mapping approach. We have compared SeroBA against real and simulated data.<br

Comparison of different technologies for the decipherment of the whole genome sequence of Campylobacter jejuni BfR-CA-14430

Background Campylobacter jejuni is a zoonotic pathogen that infects the human gut through the food chain mainly by consumption of undercooked chicken meat, raw chicken cross-contaminated ready-to-eat food or by raw milk. In the last decades, C. jejuni has increasingly become the most common bacterial cause for food-born infections in high income countries, costing public health systems billions of euros each year. Currently, different whole genome sequencing techniques such as short-read bridge amplification and long-read single molecule real-time sequencing techniques are applied for in-depth analysis of bacterial species, in particular, Illumina MiSeq, PacBio and MinION. Results In this study, we analyzed a recently isolated C. jejuni strain from chicken meat by short- and long-read data from Illumina, PacBio and MinION sequencing technologies. For comparability, this strain is used in the German PAC-CAMPY research consortium in several studies, including phenotypic analysis of biofilm formation, natural transformation and in vivo colonization models. The complete assembled genome sequence most likely consists of a chromosome of 1,645,980 bp covering 1665 coding sequences as well as a plasmid sequence with 41,772 bp that encodes for 46 genes. Multilocus sequence typing revealed that the strain belongs to the clonal complex CC-21 (ST-44) which is known to be involved in C. jejuni human infections, including outbreaks. Furthermore, we discovered resistance determinants and a point mutation in the DNA gyrase (gyrA) that render the bacterium resistant against ampicillin, tetracycline and (fluoro-)quinolones. Conclusion The comparison of Illumina MiSeq, PacBio and MinION sequencing and analyses with different assembly tools enabled us to reconstruct a complete chromosome as well as a circular plasmid sequence of the C. jejuni strain BfR-CA-14430. Illumina short-read sequencing in combination with either PacBio or MinION can substantially improve the quality of the complete chromosome and epichromosomal elements on the level of mismatches and insertions/deletions, depending on the assembly program used.Peer Reviewe