Utility of various genome lengths in diversity and evolution analyses of Hepatitis E virus.

The aim of this study was to investigate to what extent fragments of the HEV genome could be used for accurate diagnostics and inference of viral population-scale processes. For this, we selected all the published whole genome sequences from the NCBI GenBank and trimmed them to various fragment lengths (ORF1,2,3, ORF1, ORF2, ORF3, 493nt in ORF2 and 148 nt in ORF2). Each of the fragment lengths was used to infer the richness and diversity of the viral sequence types, typing accuracy, and potential use in phylodynamics. The results obtained from the different fragments were compared. We observed that, generally, the longer the nucleic acid fragment used in typing, the better the accuracy in predicting the viral subtype. However, the dominant HEV subtypes circulating in Europe were relatively well classified even by the 493nt fragment, with false negative rates as low as 8 in 1000 typed sequences. Most fragments also give comparable results in analyses of population size, albeit with shorter fragments showing a broader 95% highest posterior density interval and less obvious increase of the viral effective population size. The reconstructed phylogenies of a heterochronous subset indicated a good concordance between all the fragments, with the major clades following similar branching patterns. Furthermore, we have used the HEV sequence data from the Netherlands available in the HEVnet database as a case study for reconstruction of population size changes in the past decades. This data showed that molecular and epidemiological results are concordant and point to an increase in the viral effective population size underlying the observed increase in incidence of acute HEV infection cases. In the absence of whole genome sequencing data, the 493bp fragment can be used for analyzing HEV strains currently circulating in Europe, as it is informative for describing short term population-scale processes.

New insights into the epidemiology of Listeria monocytogenes - A cross-sectoral retrospective genomic analysis in the Netherlands (2010-2020).

Introduction: Listeriosis, caused by infection with Listeria monocytogenes (Lm), is a relatively rare but severe disease with one of the highest mortality rates among bacterial foodborne illnesses. A better understanding on the degree of Lm clustering, the temporal distribution of the clusters, and their association with the various food sources is expected to lead to improved source tracing and risk-based sampling. Methods: We investigated the genomic epidemiology of Lm in the Netherlands between 2010 and 2020 by analyzing whole-genome-sequencing (WGS) data of isolates from listerioss patients and food sources from nationwide integrated surveillance and monitoring. WGS data of 756 patient and 770 food/environmental isolates was assessed using core-genome multi-locus sequence typing (cgMLST) with Hamming distance as measure for pairwise distances. Associations of genotype with the epidemiological variables such as patient's age and gender, and systematic use of specific drugs were tested by multinomial logistic regressions. Genetic differentiation of the Lm within and between food categories was calculated based on allele frequencies at the 1701 cgMLST loci in each food category. Results: We confirmed previous results that some clonal complexes (CCs) are overrepresented among clinical isolates but could not identify any epidemiological risk factors. The main findings of this study include the observation of a very weak attribution of Lm types to food categories and a much better attribution to the producer level. In addition, we identified a high degree of temporal persistence of food, patient and mixed clusters, with more than half of the clusters spanning over more than 1 year and up to 10 years. Discussion: Taken together this would indicate that identifying persistent contamination in food production settings, and producers that process a wide variety of raw food produce, could significantly contribute to lowering the Lm disease burden.

Comparing Multiple Locus Variable-Number Tandem Repeat Analyses with Whole-Genome Sequencing as Typing Method for Salmonella Enteritidis Surveillance in The Netherlands, January 2019 to March 2020.

In the Netherlands, whole-genome sequencing (WGS) was implemented as routine typing tool for Salmonella Enteritidis isolates in 2019. Multiple locus variable-number tandem repeat analyses (MLVA) was performed in parallel. The objective was to determine the concordance of MLVA and WGS as typing methods for S. Enteritidis isolates. We included S. Enteritidis isolates from patients that were subtyped using MLVA and WGS-based core-genome Multilocus Sequence Typing (cgMLST) as part of the national laboratory surveillance of Salmonella during January 2019 to March 2020. The concordance of clustering based on MLVA and cgMLST, with a distance of ≤5 alleles, was assessed using the Fowlkes-Mallows (FM) index, and their discriminatory power using Simpson's diversity index. Of 439 isolates in total, 404 (92%) were typed as 32 clusters based on MLVA, with a median size of 4 isolates (range:2 to 141 isolates). Based on cgMLST, 313 (71%) isolates were typed as 48 clusters, with a median size of 3 isolates (range:2 to 39 isolates). The FM index was 0.34 on a scale from 0 to 1, where a higher value indicates greater similarity between the typing methods. The Simpson's diversity index of MLVA and cgMLST was 0.860 and 0.974, respectively. The median cgMLST distance between isolates with the same MLVA type was 27 alleles (interquartile range [IQR]:17 to 34 alleles), and 2 alleles within cgMLST clusters (IQR:1-5 alleles). This study shows the higher discriminatory power of WGS over MLVA and a poor concordance between both typing methods regarding clustering of S. Enteritidis isolates. IMPORTANCE Salmonella is the most frequently reported agent causing foodborne outbreaks and the second most common zoonoses in the European Union. The incidence of the most dominant serotype Enteritidis has increased in recent years. To differentiate between Salmonella isolates, traditional typing methods such as pulsed-field gel electrophoresis (PFGE) and multiple locus variable-number tandem repeat analyses (MLVA) are increasingly replaced with whole-genome sequencing (WGS). This study compared MLVA and WGS-based core-genome Multilocus Sequence Typing (cgMLST) as typing tools for S. Enteritidis isolates that were collected as part of the national Salmonella surveillance in the Netherlands. We found a higher discriminatory power of WGS-based cgMLST over MLVA, as well as a poor concordance between both typing methods regarding clustering of S. Enteritidis isolates. This is especially relevant for cluster delineation in outbreak investigations and confirmation of the outbreak source in trace-back investigations.

Sporadic Occurrence of Enteroaggregative Shiga Toxin-Producing Escherichia coli O104:H4 Similar to 2011 Outbreak Strain.

We describe the recent detection of 3 Shiga toxin-producing enteroaggregative Escherichia coli O104:H4 isolates from patients and 1 from pork in the Netherlands that were genetically highly similar to isolates from the 2011 large-scale outbreak in Europe. Our findings stress the importance of safeguarding food supply production chains to prevent future outbreaks.

Sources and transmission routes of campylobacteriosis: A combined analysis of genome and exposure data.

OBJECTIVES: To determine the contributions of several animal and environmental sources of human campylobacteriosis and identify source-specific risk factors. METHODS: 1417 Campylobacter jejuni/coli isolates from the Netherlands in 2017-2019 were whole-genome sequenced, including isolates from human cases (n = 280), chickens/turkeys (n = 238), laying hens (n = 56), cattle (n = 158), veal calves (n = 49), sheep/goats (n = 111), pigs (n = 110), dogs/cats (n = 100), wild birds (n = 62), and surface water (n = 253). Questionnaire-based exposure data was collected. Source attribution was performed using core-genome multilocus sequence typing. Risk factors were determined on the attribution estimates. RESULTS: Cases were mostly attributed to chickens/turkeys (48.2%), dogs/cats (18.0%), cattle (12.1%), and surface water (8.5%). Of the associations identified, never consuming chicken, as well as frequent chicken consumption, and rarely washing hands after touching raw meat, were risk factors for chicken/turkey-attributable infections. Consuming unpasteurized milk or barbecued beef increased the risk for cattle-attributable infections. Risk factors for infections attributable to environmental sources were open water swimming, contact with dog faeces, and consuming non-chicken/turkey avian meat like game birds. CONCLUSIONS: Poultry and cattle are the main livestock sources of campylobacteriosis, while pets and surface water are important non-livestock sources. Foodborne transmission is only partially consistent with the attributions, as frequency and alternative pathways of exposure are significant.

Genomic epidemiology of emerging ESBL-producing Salmonella Kentucky bla CTX-M-14b in Europe.

Global dissemination of ciprofloxacin-resistant Salmonella Kentucky has been observed over the past decades. In recent years, there have been reports of extended-spectrum ß-lactamase (ESBL) producing S. Kentucky. Routine surveillance at the European Centre for Disease Prevention and Control (ECDC) detected cases with a ciprofloxacin-resistant S. Kentucky with the ESBL-gene bla CTX-M-14b. Ensuing research identified 78 cases in 2013-2018 in eight European countries. Compared to other S. Kentucky and non-typhoidal Salmonella infections, reported to the European Surveillance System, these cases were more likely to be elderly and to present urinary-tract infections. Bayesian time-scaled phylogeny on whole genome sequences of isolates from these cases and supplementary isolates from public sequence databases was used to infer the origin and spread of this clone. We dated the origin of the bla CTX-M-14b clone to approximately 2005 in Northern Africa, most likely in Egypt. The geographic origin predicted by the phylogenetic analysis is consistent with the patients' travel history. Next to multiple introductions of the clone to Europe from Egypt, our analysis suggests that in some parts of Europe the clone might have formed a stable population, from which further spread has occurred. Comparative genomics indicated that the bla CTX-M-14b gene is present on the bacterial chromosome, within the type VI secretion system region. The bla CTX-M-14b gene is integrated downstream of the hcp1 gene, on a 2854 bp plasmid fragment containing also ISEcp1. This is the first report of a chromosomally integrated CTX-M gene in Salmonella spp. in Europe, previous studies having identified similar genes only on plasmids.

Tracing the animal sources of surface water contamination with Campylobacter jejuni and Campylobacter coli.

Campylobacter jejuni and C. coli, the primary agents of human bacterial gastroenteritis worldwide, are widespread in surface water. Several animal sources contribute to surface water contamination with Campylobacter, but their relative contributions thus far remained unclear. Here, the prevalence, genotype diversity, and potential animal sources of C. jejuni and C. coli strains in surface water in the Netherlands were investigated. It was also assessed whether the contribution of the different animal sources varied according to surface water type (i.e. agricultural water, surface water at discharge points of wastewater treatment plants [WWTPs], and official recreational water), season, and local livestock (poultry, pig, ruminant) density. For each surface water type, 30 locations spread over six areas with either high or low density of poultry, ruminants, or pigs, were sampled once every season in 2018-2019. Campylobacter prevalence was highest in agricultural waters (77%), and in autumn and winter (74%), and lowest in recreational waters (46%) and in summer (54%). In total, 76 C. jejuni and 177 C. coli water isolates were whole-genome sequenced. Most C. coli water isolates (78.5%) belonged to hitherto unidentified clones when using the seven-locus sequence type (ST) scheme, while only 11.8% of the C. jejuni isolates had unidentified STs. The origin of these isolates, as defined by core-genome multi-locus sequence typing (cgMLST), was inferred by comparison with Campylobacter strain collections from meat-producing poultry, laying hens, adult cattle, veal calves, small ruminants, pigs, and wild birds. Water isolates were mainly attributed to wild birds (C. jejuni: 60.0%; C. coli: 93.7%) and meat-producing poultry (C. jejuni: 18.9%; C. coli: 5.6%). Wild bird contribution was high among isolates from recreational waters and WWTP discharge points, and in areas with low poultry (C. coli) or high ruminant (C. jejuni) densities. The contribution of meat-producing poultry was high in areas with high density of poultry, springtime, agricultural waters and WWTP discharge points. While wild birds and poultry were the main contributors to Campylobacter contamination in surface water, their contribution differed significantly by water type, season, and local poultry and ruminant densities.

Concordance of SNP- and allele-based typing workflows in the context of a large-scale international Salmonella Enteritidis outbreak investigation.

A large European multi-country Salmonella enterica serovar Enteritidis outbreak associated with Polish eggs was characterized by whole-genome sequencing (WGS)-based analysis, with various European institutes using different analysis workflows to identify isolates potentially related to the outbreak. The objective of our study was to compare the output of six of these different typing workflows (distance matrices of either SNP-based or allele-based workflows) in terms of cluster detection and concordance. To this end, we analysed a set of 180 isolates coming from confirmed and probable outbreak cases, which were representative of the genetic variation within the outbreak, supplemented with 22 unrelated contemporaneous S. enterica serovar Enteritidis isolates. Since the definition of a cluster cut-off based on genetic distance requires prior knowledge on the evolutionary processes that govern the bacterial populations in question, we used a variety of hierarchical clustering methods (single, average and complete) and selected the optimal number of clusters based on the consensus of the silhouette, Dunn2, and McClain-Rao internal validation indices. External validation was done by calculating the concordance with the WGS-based case definition (SNP-address) for this outbreak using the Fowlkes-Mallows index. Our analysis indicates that with complete-linkage hierarchical clustering combined with the optimal number of clusters, as defined by three internal validity indices, the six different allele- and SNP-based typing workflows generate clusters with similar compositions. Furthermore, we show that even in the absence of coordinated typing procedures, but by using an unsupervised machine learning methodology for cluster delineation, the various workflows that are currently in use by six European public-health authorities can identify concordant clusters of genetically related S. enterica serovar Enteritidis isolates; thus, providing public-health researchers with comparable tools for detection of infectious-disease outbreaks.

Comparative genomics reveals a lack of evidence for pigeons as a main source of stx2f-carrying Escherichia coli causing disease in humans and the common existence of hybrid Shiga toxin-producing and enteropathogenic E. coli pathotypes.

BACKGROUND: Wild birds, in particular pigeons are considered a natural reservoir for stx2f-carrying E. coli. An extensive comparison of isolates from pigeons and humans from the same region is lacking, which hampers justifiable conclusions on the epidemiology of these pathogens. Over two hundred human and pigeon stx2f-carrying E. coli isolates predominantly from the Netherlands were analysed by whole genome sequencing and comparative genomic analysis including in silico MLST, serotyping, virulence genes typing and whole genome MLST (wgMLST). RESULTS: Serotypes and sequence types of stx2f-carrying E. coli showed a strong non-random distribution among the human and pigeon isolates with O63:H6/ST583, O113:H6/ST121 and O125:H6/ST583 overrepresented among the human isolates and not found among pigeons. Pigeon isolates were characterized by an overrepresentation of O4:H2/ST20 and O45:H2/ST20. Nearly all isolates harboured the locus of enterocyte effacement (LEE) but different eae and tir subtypes were non-randomly distributed among human and pigeon isolates. Phylogenetic core genome comparison demonstrated that the pigeon isolates and clinical isolates largely occurred in separated clusters. In addition, serotypes/STs exclusively found among humans generally were characterized by high level of clonality, smaller genome sizes and lack of several non-LEE-encoded virulence genes. A bundle-forming pilus operon, including bfpA, indicative for typical enteropathogenic E. coli (tEPEC) was demonstrated in 72.0% of the stx2f-carrying serotypes but with distinct operon types between the main pigeon and human isolate clusters. CONCLUSIONS: Comparative genomics revealed that isolates from mild human disease are dominated by serotypes not encountered in the pigeon reservoir. It is therefore unlikely that zoonotic transmission from this reservoir plays an important role in the contribution to the majority of human disease associated with stx2f-producing E. coli in the Netherlands. Unexpectedly, this study identified the common occurrence of STEC2f/tEPEC hybrid pathotype in various serotypes and STs. Further research should focus on the possible role of human-to-human transmission of Stx2f-producing E. coli.

The genetic diversity of Borrelia afzelii is not maintained by the diversity of the rodent hosts.

BACKGROUND: Small mammals are essential in the enzootic cycle of many tick-borne pathogens (TBP). To understand their contribution to the genetic diversity of Borrelia afzelii, the most prevalent TBP in questing Ixodes ricinus, we compared the genetic variants of B. afzelii at three distinct genetic loci. We chose two plasmid loci, dbpA and ospC, and a chromosomal one, IGS. RESULTS: While the larvae that fed on shrews (Sorex sp.) tested negative for B. afzelii, those fed on bank voles (Myodes glareolus) and wood mice (Apodemus sylvaticus) showed high infection prevalences of 0.13 and 0.27, respectively. Despite the high genetic diversity within B. afzelii, there was no difference between wood mice and bank voles in the number and types of B. afzelii haplotypes they transmit. CONCLUSIONS: The genetic diversity in B. afzelii cannot be explained by separate enzootic cycles in wood mice and bank voles.

Borrelia miyamotoi and Co-Infection with Borrelia afzelii in Ixodes ricinus Ticks and Rodents from Slovakia.

Borrelia miyamotoi causes relapsing fever in humans. The occurrence of this spirochete has been reported in Ixodes ricinus and wildlife, but there are still gaps in the knowledge of its eco-epidemiology and public health impact. In the current study, questing I. ricinus (nymphs and adults) and skin biopsies from rodents captured in Slovakia were screened for the presence of B. miyamotoi and Borrelia burgdorferi s.l. DNA. The prevalence of B. miyamotoi and B. burgdorferi s.l. in questing ticks was 1.7 and 16.9%, respectively. B. miyamotoi was detected in Apodemus flavicollis (9.3%) and Myodes glareolus (4.4%). In contrast, B. burgdorferi s.l. was identified in 11.9% of rodents, with the highest prevalence in Microtus arvalis (68.4%) and a lower prevalence in Apodemus spp. (8.4%) and M. glareolus (12.4%). Borrelia afzelii was the prevailing genospecies infecting questing I. ricinus (37.9%) and rodents (72.2%). Co-infections of B. miyamotoi and B. burgdorferi s.l. were found in 24.1 and 9.3% of the questing ticks and rodents, respectively, whereas the proportion of ticks and rodents co-infected with B. miyamotoi and B. afzelii was 6.9 and 7.0%, respectively. The results suggest that B. miyamotoi and B. afzelii share amplifying hosts. The sequences of the B. miyamotoi glpQ gene fragment from our study showed a high degree of identity with sequences of the gene amplified from ticks and human patients in Europe. The results seem to suggest that humans in Slovakia are at risk of contracting tick-borne relapsing fever, and in some cases together with Lyme borreliosis.

Larvae of Ixodes ricinus transmit Borrelia afzelii and B. miyamotoi to vertebrate hosts.

BACKGROUND: Lyme borreliosis is the most common tick-borne human disease and is caused by Borrelia burgdorferi sensu lato (s.l.). Borrelia miyamotoi, a relapsing fever spirochaete, is transmitted transovarially, whereas this has not been shown for B. burgdorferi (s.l). Therefore, B. burgdorferi (s.l) is considered to cycle from nymphs to larvae through vertebrates. Larvae of Ixodes ricinus are occasionally B. burgdorferi (s.l) infected, but their vector competence has never been studied. METHODS: We challenged 20 laboratory mice with field-collected larvae of I. ricinus. A subset of these larvae was analysed for infections with B. burgdorferi (s.l) and B. miyamotoi. After three to four challenges, mice were sacrificed and skin and spleen samples were analysed for infection by PCR and culture. RESULTS: Field-collected larvae were naturally infected with B. burgdorferi (s.l) (0.62%) and B. miyamotoi (2.0%). Two mice acquired a B. afzelii infection and four mice acquired a B. miyamotoi infection during the larval challenges. CONCLUSION: We showed that larvae of I. ricinus transmit B. afzelii and B. miyamotoi to rodents and calculated that rodents have a considerable chance of acquiring infections from larvae compared to nymphs. As a result, B. afzelii can cycle between larvae through rodents. Our findings further imply that larval bites on humans, which easily go unnoticed, can cause Lyme borreliosis and Borrelia miyamotoi disease.