Host and environmental factors shape upper airway microbiota and respiratory health across the human lifespan.

Our understanding of the normal variation in the upper respiratory tract (URT) microbiota across the human lifespan and how these relate to host, environment, and health is limited. We studied the microbiota of 3,104 saliva (<10 year-olds)/oropharynx (≥10 year-olds) and 2,485 nasopharynx samples of 3,160 Dutch individuals 0-87 years of age, participating in a cross-sectional population-wide study (PIENTER-3) using 16S-rRNA sequencing. The microbiota composition was strongly related to age, especially in the nasopharynx, with maturation occurring throughout childhood and adolescence. Clear niche- and age-specific associations were found between the microbiota composition and host/environmental factors and health outcomes. Among others, social interaction, sex, and season were associated with the nasopharyngeal microbial community. By contrast, the oral microbiota was more related to antibiotics, tobacco, and alcohol use. We present an atlas of the URT microbiota across the lifespan in association with environment and health, establishing a baseline for future research.

Impact of age on pneumococcal colonization of the nasopharynx and oral cavity: an ecological perspective.

Pneumococcal carriage studies have suggested that pneumococcal colonization in adults is largely limited to the oral cavity and oropharynx. In this study, we used total abundance-based ß-diversity (dissimilarity) and ß-diversity components to characterize age-related differences in pneumococcal serotype composition of respiratory samples. quantitative PCR (qPCR) was applied to detect pneumococcal serotypes in nasopharyngeal samples collected from 946 toddlers and 602 adults, saliva samples collected from a subset of 653 toddlers, and saliva and oropharyngeal samples collected from a subset of 318 adults. Bacterial culture rates from nasopharyngeal samples were used to characterize age-related differences in rates of colonizing bacteria. Dissimilarity in pneumococcal serotype composition was low among saliva and nasopharyngeal samples from children. In contrast, respiratory samples from adults exhibited high serotype dissimilarity, which predominantly consisted of abundance gradients and was associated with reduced nasopharyngeal colonization. Age-related serotype dissimilarity was high among nasopharyngeal samples and relatively low for saliva samples. Reduced nasopharyngeal colonization by pneumococcal serotypes coincided with significantly reduced Moraxella catarrhalis and Haemophilus influenzae and increased Staphylococcus aureus nasopharyngeal colonization rates among adults. Findings from this study suggest that within-host environmental conditions, utilized in the upper airways by pneumococcus and other bacteria, undergo age-related changes. It may result in a host-driven ecological succession of bacterial species colonizing the nasopharynx and lead to competitive exclusion of pneumococcus from the nasopharynx but not from the oral habitat. This explains the poor performance of nasopharyngeal samples for pneumococcal carriage among adults and indicates that in adults saliva more accurately represents the epidemiology of pneumococcal carriage than nasopharyngeal samples.

Aspergillus Outbreak in an Intensive Care Unit: Source Analysis with Whole Genome Sequencing and Short Tandem Repeats.

Whole genome sequencing (WGS) is widely used for outbreak analysis of bacteriology and virology but is scarcely used in mycology. Here, we used WGS for genotyping Aspergillus fumigatus isolates from a potential Aspergillus outbreak in an intensive care unit (ICU) during construction work. After detecting the outbreak, fungal cultures were performed on all surveillance and/or patient respiratory samples. Environmental samples were obtained throughout the ICU. WGS was performed on 30 isolates, of which six patient samples and four environmental samples were related to the outbreak, and twenty samples were unrelated, using the Illumina NextSeq 550. A SNP-based phylogenetic tree was created from outbreak samples and unrelated samples. Comparative analysis (WGS and short tandem repeats (STRs), microsatellite loci analysis) showed that none of the strains were related to each other. The lack of genetic similarity suggests the accumulation of Aspergillus spores in the hospital environment, rather than a single source that supported growth and reproduction of Aspergillus fumigatus. This supports the hypothesis that the Aspergillus outbreak was likely caused by release of Aspergillus fumigatus spores during construction work. Indeed, no new Aspergillus cases were observed in the ICU after cessation of construction. This study demonstrates that WGS is a suitable technique for examining inter-strain relatedness of Aspergillus fumigatus in the setting of an outbreak investigation.

Genomic Epidemiology Identifies Azole Resistance Due to TR34/L98H in European Aspergillus fumigatus Causing COVID-19-Associated Pulmonary Aspergillosis.

Aspergillus fumigatus has been found to coinfect patients with severe SARS-CoV-2 virus infection, leading to COVID-19-associated pulmonary aspergillosis (CAPA). The CAPA all-cause mortality rate is approximately 50% and may be complicated by azole resistance. Genomic epidemiology can help shed light on the genetics of A. fumigatus causing CAPA, including the prevalence of resistance-associated alleles. We present a population genomic analysis of 21 CAPA isolates from four European countries with these isolates compared against 240 non-CAPA A. fumigatus isolates from a wider population. Bioinformatic analysis and antifungal susceptibility testing were performed to quantify resistance and identify possible genetically encoded azole-resistant mechanisms. The phylogenetic analysis of the 21 CAPA isolates showed that they were representative of the wider A. fumigatus population with no obvious clustering. The prevalence of phenotypic azole resistance in CAPA was 14.3% (n = 3/21); all three CAPA isolates contained a known resistance-associated cyp51A polymorphism. The relatively high prevalence of azole resistance alleles that we document poses a probable threat to treatment success rates, warranting the enhanced surveillance of A. fumigatus genotypes in these patients. Furthermore, potential changes to antifungal first-line treatment guidelines may be needed to improve patient outcomes when CAPA is suspected.

Molecular characterization of MRSA collected during national surveillance between 2008 and 2019 in the Netherlands.

BACKGROUND: Although the Netherlands is a country with a low endemic level, methicillin-resistant Staphylococcus aureus (MRSA) poses a significant health care problem. Therefore, high coverage national MRSA surveillance has been in place since 1989. To monitor possible changes in the type-distribution and emergence of resistance and virulence, MRSA isolates are molecularly characterized. METHODS: All 43,321 isolates from 36,520 persons, collected 2008-2019, were typed by multiple-locus variable number tandem repeats analysis (MLVA) with simultaneous PCR detection of the mecA, mecC and lukF-PV genes, indicative for PVL. Next-generation sequencing data of 4991 isolates from 4798 persons were used for whole genome multi-locus sequence typing (wgMLST) and identification of resistance and virulence genes. RESULTS: We show temporal change in the molecular characteristics of the MRSA population with the proportion of PVL-positive isolates increasing from 15% in 2008-2010 to 25% in 2017-2019. In livestock-associated MRSA obtained from humans, PVL-positivity increases to 6% in 2017-2019 with isolates predominantly from regions with few pig farms. wgMLST reveals the presence of 35 genogroups with distinct resistance, virulence gene profiles and specimen origin. Typing shows prolonged persistent MRSA carriage with a mean carriage period of 407 days. There is a clear spatial and a weak temporal relationship between isolates that clustered in wgMLST, indicative for regional spread of MRSA strains. CONCLUSIONS: Using molecular characterization, this exceptionally large study shows genomic changes in the MRSA population at the national level. It reveals waxing and waning of types and genogroups and an increasing proportion of PVL-positive MRSA.

A group of bacteria that cause difficult-to-treat infections in humans is methicillin-resistant Staphylococcus aureus (MRSA). The aim of this study was to monitor changes in the spread of MRSA, their disease causing potential and resistance to antibiotics used to treat MRSA infections. MRSA from patients and their contacts in the Netherlands were collected over a period of 12 years and characterized. This revealed new types of MRSA emerged and others disappeared. An increasing number of MRSA produces a protein called PVL toxin, enabling MRSA to cause more severe infections. Also, some people appear to carry MRSA without any disease for more than a year. These findings suggest an increasing disease potential of MRSA and possible unnoticed sources of infection. Consequently, it is important to maintain monitoring of these infections to minimize MRSA spread.

Surveillance of Neisseria meningitidis carriage four years after menACWY vaccine implementation in the Netherlands reveals decline in vaccine-type and rise in genogroup e circulation.

Carriage of Neisseria meningitidisis an accepted endpoint in monitoring meningococcal vaccine effects. We applied molecular methods to assess the impact of menACWY vaccine implementation on meningococcal carriage and genogroup-specific prevalence in young adults in Fall of 2022, four years after the introduction of the tetravalent vaccine in the Netherlands. The overall carriage rate of genogroupable meningococci was not significantly different compared to a pre-menACWY cohort investigated in 2018 (20.8 % or 125 of 601 versus 17.4 % or 52 of 299 individuals, p = 0.25). Of 125 carriers of genogroupable meningococci, 122 (97.6 %) were positive for either vaccine-types menC, menW, menY or genogroups, menB, menE, and menX, which are not targeted by the menACWY vaccine. Compared with a pre-vaccine-implementation cohort, there was 3.8-fold reduction (p < 0.001) in vaccine-type carriage rates and 9.0-fold increase (p < 0.0001) in non-vaccine type menE prevalence. We observe a reduction in menW and menY and an increase in menE, which suggest that implementation of menACWY vaccine affected carriage.

Widespread Detection of Yersiniabactin Gene Cluster and Its Encoding Integrative Conjugative Elements (ICEKp) among Nonoutbreak OXA-48-Producing Klebsiella pneumoniae Clinical Isolates from Spain and the Netherlands.

In this study, we determined the presence of virulence factors in nonoutbreak, high-risk clones and other isolates belonging to less common sequence types associated with the spread of OXA-48-producing Klebsiella pneumoniae clinical isolates from The Netherlands (n = 61) and Spain (n = 53). Most isolates shared a chromosomally encoded core of virulence factors, including the enterobactin gene cluster, fimbrial fim and mrk gene clusters, and urea metabolism genes (ureAD). We observed a high diversity of K-Locus and K/O loci combinations, KL17 and KL24 (both 16%), and the O1/O2v1 locus (51%) being the most prevalent in our study. The most prevalent accessory virulence factor was the yersiniabactin gene cluster (66.7%). We found seven yersiniabactin lineages-ybt 9, ybt 10, ybt 13, ybt 14, ybt 16, ybt 17, and ybt 27-which were chromosomally embedded in seven integrative conjugative elements (ICEKp): ICEKp3, ICEKp4, ICEKp2, ICEKp5, ICEKp12, ICEKp10, and ICEKp22, respectively. Multidrug-resistant lineages-ST11, ST101, and ST405-were associated with ybt 10/ICEKp4, ybt 9/ICEKp3, and ybt 27/ICEKp22, respectively. The fimbrial adhesin kpi operon (kpiABCDEFG) was predominant among ST14, ST15, and ST405 isolates, as well as the ferric uptake system kfuABC, which was also predominant among ST101 isolates. No convergence of hypervirulence and resistance was observed in this collection of OXA-48-producing K. pneumoniae clinical isolates. Nevertheless, two isolates, ST133 and ST792, were positive for the genotoxin colibactin gene cluster (ICEKp10). In this study, the integrative conjugative element, ICEKp, was the major vehicle for yersiniabactin and colibactin gene clusters spreading. IMPORTANCE Convergence of multidrug resistance and hypervirulence in Klebsiella pneumoniae isolates has been reported mostly related to sporadic cases or small outbreaks. Nevertheless, little is known about the real prevalence of carbapenem-resistant hypervirulent K. pneumoniae since these two phenomena are often separately studied. In this study, we gathered information on the virulent content of nonoutbreak, high-risk clones (i.e., ST11, ST15, and ST405) and other less common STs associated with the spread of OXA-48-producing K. pneumoniae clinical isolates. The study of virulence content in nonoutbreak isolates can help us to expand information on the genomic landscape of virulence factors in K. pneumoniae population by identifying virulence markers and their mechanisms of spread. Surveillance should focus not only on antimicrobial resistance but also on virulence characteristics to avoid the spread of multidrug and (hyper)virulent K. pneumoniae that may cause untreatable and more severe infections.

A spitting image: molecular diagnostics applied to saliva enhance detection of Streptococcus pneumoniae and pneumococcal serotype carriage.

Background: Despite strong historical records on the accuracy of saliva testing, oral fluids are considered poorly suited for pneumococcal carriage detection. We evaluated an approach for carriage surveillance and vaccine studies that increases the sensitivity and specificity of pneumococcus and pneumococcal serotype detection in saliva samples. Methods: Quantitative PCR (qPCR)-based methods were applied to detect pneumococcus and pneumococcal serotypes in 971 saliva samples collected from 653 toddlers and 318 adults. Results were compared with culture-based and qPCR-based detection in nasopharyngeal samples collected from children and in nasopharyngeal and oropharyngeal samples collected from adults. Optimal C q cut-offs for positivity in qPCRs were determined via receiver operating characteristic curve analysis and accuracy of different approaches was assessed using a composite reference for pneumococcal and for serotype carriage based on isolation of live pneumococcus from the person or positivity of saliva samples determined with qPCR. To evaluate the inter-laboratory reproducibility of the method, 229 culture-enriched samples were tested independently in the second center. Results: In total, 51.5% of saliva samples from children and 31.8% of saliva samples from adults were positive for pneumococcus. Detection of pneumococcus by qPCR in culture-enriched saliva exhibited enhanced sensitivity and higher agreement with a composite reference compared to diagnostic culture of nasopharyngeal samples in children (Cohen's κ: 0.69-0.79 vs. 0.61-0.73) and in adults (κ: 0.84-0.95 vs. 0.04-0.33) and culture of oropharyngeal samples in adults (κ: 0.84-0.95 vs. -0.12-0.19). Similarly, detection of serotypes with qPCR in culture-enriched saliva exhibited enhanced sensitivity and higher agreement with a composite reference compared to nasopharyngeal culture in children (κ: 0.73-0.82 vs. 0.61-0.73) and adults (κ: 0.90-0.96 vs. 0.00-0.30) and oropharyngeal culture in adults (κ: 0.90-0.96 vs. -0.13 to 0.30). However, results of qPCRs targeting serotype 4, 5, and 17F and serogroups 9, 12, and 35 were excluded due to assays' lack of specificity. We observed excellent quantitative agreement for qPCR-based detection of pneumococcus between laboratories. After exclusion of serotype/serogroup-specific assays with insufficient specificity, moderate agreement (κ 0.68, 95% CI 0.58-0.77) was observed. Conclusion: Molecular testing of culture-enriched saliva samples improves the sensitivity of overall surveillance of pneumococcal carriage in children and adults, but limitations of qPCR-based approaches for pneumococcal serotypes carriage detection should be considered.

Whole-genome sequencing of Shigella for surveillance purposes shows (inter)national relatedness and multidrug resistance in isolates from men who have sex with men.

In the Netherlands, more than half of domestic shigellosis cases are among men who have sex with men (MSM), particularly in the Amsterdam region. However, there is limited insight into which Shigella strains circulate in the Netherlands. Our objective was to assess the added value of whole-genome sequencing (WGS)-based surveillance for Shigella. To this end, we determined the relatedness among Shigella spp. isolates from patients in the Amsterdam region, as well as in an international context, including antimicrobial resistance markers, using WGS. The following criteria were used: it should provide insight into (1) clustering of shigellosis cases and the affected population, (2) the extent of admixture of MSM-associated isolates with those from the broader population and (3) the presence of antimicrobial resistance. It should then lead to more opportunities for targeted control measures. For this study, Shigella isolates from three laboratories in the Amsterdam region obtained between February 2019 and October 2021 were subjected to Illumina WGS at the National Institute for Public Health and the Environment (RIVM). Raw data were quality-checked and assembled, the Shigella serotype was determined with ShigaTyper, and antimicrobial resistance markers were detected using ResFinder and PointFinder. For Shigella sonnei, subclades were determined using Mykrobe. Relatedness of isolates, including 21 international reference genomes, was assessed with core genome multilocus sequence typing. In total, 109 isolates were included, of which 27 were from females (25 %) and 66 were from males (61 %), with which the majority (n=48, 73 %) being from MSM. No information on sex was available for the remaining 16 cases. The WGS data for all isolates, comprising 55 S. sonnei, 52 Shigella flexneri, 1 Shigella boydii and 1 Shigella dysenteriae, met the quality criteria. In total, 14 clusters containing 51 isolates (49 %) were identified, with a median cluster size of 2.5 cases (range: 2-15). Nine out of 14 clusters were MSM-associated, and 8 clusters (57 %) were travel-related. Six of the MSM clusters were related to international reference genomes. The prevalence of antimicrobial resistance markers was higher among isolates from MSM than non-MSM patients, particularly for ciprofloxacin (89 vs 33 %) and azithromycin (58 vs 17 %). In conclusion, about half of Shigella spp. patients were part of a cluster, of which a substantial part were related to international reference genomes, particularly among MSM, and a high prevalence of antimicrobial resistance markers was found. These findings indicate widespread international circulation of Shigella spp., particularly among MSM, with multidrug resistance that hampers treatment of patients. Moreover, the results of this study led to the implementation of a national WGS-based laboratory surveillance programme for Shigella spp. that started in April 2022.

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.

Panton-Valentine Leukocidin-Positive CC398 MRSA in Urban Clinical Settings, the Netherlands.

We report detection of Panton-Valentine leukocidin-positive clonal complex 398 human-origin methicillin-resistant Staphylococcus aureus L2 in the Netherlands. This hypervirulent lineage originated in the Asia-Pacific Region and could become community-acquired in Europe after recurrent travel-related introductions. Genomic surveillance enables early detection to guide control measures and help limit spread of pathogens in urban settings.

Comparative analysis of IMP-4- and OXA-58-containing plasmids of three carbapenemase-producing Acinetobacter ursingii strains in the Netherlands.

OBJECTIVES: A recent occurrence of carbapenemase-producing Acinetobacter ursingii was reported in the Netherlands and comprised three unrelated strains carrying the blaIMP-4 and blaOXA-58 encoding genes. The objective was to investigate a putative common source of the carbapenemase resistance genes and plasmids in these A. ursingii strains. METHODS: Hybrid assembly of short-read and long-read sequencing data was performed using Unicycler and assembled genomes were analysed by ResFinder and PlasmidFinder. RESULTS: Hybrid assemblies of A. ursingii genomes yielded a circular chromosome, a large plasmid harboring blaIMP-4 and blaOXA-58 genes (sizes 259-317kb), and four to five other smaller plasmids. ResFinder analyses revealed 16 other acquired resistance genes on the plasmids carrying the blaIMP-4 and blaOXA-58 genes. These 18 genes encode resistance towards eight antibiotic classes. The smaller plasmids did not carry acquired resistance genes. Comparative analysis showed that the three blaIMP-4/blaOXA-58 plasmids were similar (61%-83%) and shared 13 to 17 of the 18 resistance genes. BLAST analysis showed that the blaIMP-4/blaOXA-58 plasmids were not reported before. However, a close match with a 399 kb plasmid from Acinetobacter johnsonii was found (99% similarity, 80% coverage). This A. johnsonii plasmid contains the blaOXA-58 gene, but lacks blaIMP-4, and it shares eight other resistance genes with those present on the A. ursingii blaIMP-4/blaOXA-58 plasmids. CONCLUSION: Three blaIMP-4/blaOXA-58-carrying plasmids were characterized in three carbapenemase-producing A. ursingii strains. The plasmids were highly similar, suggesting a putative common source or co-selection of resistance genes from A. johnsonii. These results provide initial insights in the dissemination of carbapenem-resistance in A. ursingii in the Netherlands.

Higher off-target amplicon detection rate in MiSeq v3 compared to v2 reagent kits in the context of 16S-rRNA-sequencing.

One of the most widely used techniques in microbiota research is 16S-rRNA-sequencing. Several laboratory processes have been shown to impact sequencing results, especially in low biomass samples. Low biomass samples are prone to off-target amplification, where instead of bacterial DNA, host DNA is erroneously amplified. Knowledge on the laboratory processes influencing off-target amplification and detection is however scarce. We here expand on previous findings by demonstrating that off-target amplification is not limited to invasive biopsy samples, but is also an issue in low bacterial biomass respiratory (mucosal) samples, especially when below 0.3 pg/µL. We show that off-target amplification can partly be mitigated by using gel-based library purification methods. Importantly, we report a higher off-target amplicon detection rate when using MiSeq reagent kit v3 compared to v2 (mean 13.3% vs 0.1% off-target reads/sample, respectively), possibly as a result of differences in reagents or sequencing recipes. However, since after bioinformatic removal of off-target reads, MiSeq reagent kit v3 still results in a twofold higher number of reads when compared to v2, v3 is still preferred over v2. Together, these results add to the growing knowledge base on off-target amplification and detection, allowing researchers to anticipate this problem in 16S-rRNA-based microbiome studies involving low biomass samples.

Changing epidemiology of Salmonella Enteritidis human infections in the Netherlands and Belgium, 2006 to 2019: a registry-based population study.

BackgroundSalmonellosis remains the second most common zoonosis in the European Union despite a long-term decreasing trend. However, this trend has been reported to have stagnated in recent years, particularly for Salmonella enterica serotype Enteritidis (SE).AimTo describe temporal changes in the incidence of SE human infections, and in its associated factors between 2006 and 2019. In addition, we aim to determine which factors influenced the stagnated trend seen in recent years.MethodsData on culture-confirmed SE human infections from national surveillance registries in the Netherlands and Belgium between 2006 and 2019 were analysed using multivariable negative-binomial regression models with restricted cubic splines.ResultsSE incidence was significantly higher in summer and autumn than winter, in persons aged 0-4 years and 5-14 years than in persons ≥ 60 years, and increased with increasing proportions of travel-related and resistant SE infections. SE incidence decreased significantly in both countries until 2015, followed by an increasing trend, which was particularly pronounced in the Netherlands. Potential SE outbreaks in both countries and invasive infections in the Netherlands also increased after 2015.ConclusionThe increase in potential outbreaks and invasive infections since 2015 may partially explain the observed reversal of the decreasing trend. While these results provide insights into the possible causes of this trend reversal, attention should also be given to factors known to influence SE epidemiology at primary (animal) production and pathogen genomic levels.

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.

External validation of WGS-based antimicrobial susceptibility prediction tools, KOVER-AMR and ResFinder 4.1, for Escherichia coli clinical isolates.

OBJECTIVE: To externally validate whole genome sequence-antimicrobial susceptibility testing phenotype prediction tools KOVER-AMR and ResFinder 4.1 for Escherichia coli clinical isolates from Dutch routine care. METHODS: A random sample of 234 E. coli and 283 third generation cephalosporin-resistant E. coli isolates from urine and blood were collected (2014-2017). Culture-antimicrobial susceptibility testing was performed using VITEK 2 and BD Phoenix. Sequences were used as input for KOVER-AMR-SCM, KOVER-AMR-CART, and ResFinder 4.1. The concordance, major error rate (MER), and very major error rate (VMER) were calculated, with subsequent comparison to U.S. Food and Drug Administration (FDA) criteria (MER ≤3% and VMER with a 95% confidence interval ≤1.5-≤7.5%). RESULTS: ResFinder 4.1 performed better than KOVER-AMR-models; however, neither tool achieved overall (V)MERs below FDA criteria. KOVER-AMR-SCM, KOVER-AMR-CART, and ResFinder 4.1, MER (cumulative all antimicrobials) were: 5.1% (4.4-5.9), 4.3% (3.6-5.0), and 5.1% (4.5-5.8), respectively. MERs ≤3% were achieved for 6 (SCM) and 5 (CART) of the 11 tested antimicrobials for KOVER-AMR-models and for 9/13 antimicrobials tested with ResFinder 4.1. KOVER-AMR-SCM, KOVER-AMR-CART, and ResFinder 4.1 cumulative VMERs were: 26% (24-28), 29% (27-31), and 11% (9.2-12). VMERs with a 95% CI ≤ 1.5-≤7.5 were only achieved for 4/13 tested antimicrobials with ResFinder 4.1. DISCUSSION: In this study, whole genome sequence-antimicrobial susceptibility testing phenotype prediction tools KOVER-AMR and ResFinder 4.1 did not meet the FDA criteria needed for clinical diagnostic use in 517 E. coli clinical isolates from Dutch routine care. The tested tools should be further improved before they can be used for clinical decision making.

Global spatial dynamics and vaccine-induced fitness changes of Bordetella pertussis.

As with other pathogens, competitive interactions between Bordetella pertussis strains drive infection risk. Vaccines are thought to perturb strain diversity through shifts in immune pressures; however, this has rarely been measured because of inadequate data and analytical tools. We used 3344 sequences from 23 countries to show that, on average, there are 28.1 transmission chains circulating within a subnational region, with the number of chains strongly associated with host population size. It took 5 to 10 years for B. pertussis to be homogeneously distributed throughout Europe, with the same time frame required for the United States. Increased fitness of pertactin-deficient strains after implementation of acellular vaccines, but reduced fitness otherwise, can explain long-term genotype dynamics. These findings highlight the role of vaccine policy in shifting local diversity of a pathogen that is responsible for 160,000 deaths annually.

Associations and recovery dynamics of the nasopharyngeal microbiota during influenza-like illness in the aging population.

Influenza-like illness (ILI), a disease caused by respiratory pathogens including influenza virus, is a major health concern in older adults. There is little information on changes and recovery dynamics of the nasopharyngeal (NP) microbiota of older adults associated with an ILI. Here, we compared the NP microbiota in older adults reporting (n = 240) or not (n = 157) ILI during the 2014-2015 influenza season at different times of the ILI event. A small but significant effect of the ILI was observed on the microbiota community composition and structure when compared to controls and samples collected at recovery. Corynebacterium was negatively associated with ILI and its abundance increased after recovery. Potential pathobionts such as Haemophilus, Porphyromonas and Gemella had higher abundances during acute-ILI. Stability and changes in the NP microbial community showed individual dynamics. Key core genera, Corynebacterium, Moraxella and Dolosigranulum exhibited higher inter-individual variability in acute-ILI, but showed comparable variability to controls after recovery. Participants in the ILI group with higher core microbiota abundances at the acute phase showed higher microbiota stability after recovery. Our findings demonstrate that acute-ILI is associated with alterations in the phylogenetic structure of the NP microbiota in older adults. The variation in the core microbiota suggests imbalances in the ecosystem, which could potentially play a role in the susceptibility and recovery of the NP microbiota after an ILI event.

Livestock-associated methicillin-resistant Staphylococcus aureus epidemiology, genetic diversity, and clinical characteristics in an urban region.

Objectives: While Livestock-associated methicillin-resistant Staphylococcus aureus (LA-MRSA), defined as CC398, is a well-known pathogen among those working with livestock, there are indications that LA-MRSA prevalence among the general population is increasing. However, the clinical impact in urban areas remains unknown. The aim of this study was to assess the genetic epidemiology and clinical characteristics of LA-MRSA in an urban area with a limited livestock population. Methods: In this retrospective study, we evaluated LA-MRSA strains that were collected between 2014 and 2018 from patients who received clinical care in a single urban area in Netherlands. Patient files were assessed for livestock exposure data, clinical findings, and contact tracing information. Next-generation sequencing (NGS) analysis in combination with wgMLST was conducted to assess genetic diversity and relatedness and to detect virulence and resistance genes. Results: LA-MRSA strains were cultured from 81 patients, comprising 12% of all the MRSA strains found in seven study laboratories between 2014 and 2018. No livestock link was found in 76% of patients (n = 61), and 28% of patients (n = 23) had an infection, mostly of the skin or soft tissue. Contact tracing had been initiated in 14 cases, leading to the identification of two hospital transmissions: a cluster of 9 cases and one of 2 cases. NGS data were available for 91% (n = 75) of the patients. wgMLST confirmed the clusters detected via contact tracing (n = 2) and identified 5 additional clusters without a known epidemiological link. Relevant resistance and virulence findings included the PVL virulence gene (3 isolates) and tetracycline resistance (79 isolates). Conclusion: LA-MRSA may cause a relevant burden of disease in urban areas. Surprisingly, most infections in the present study occurred in the absence of a livestock link, suggesting inter-human transmission. These findings and the presence of PVL and other immune evasive complex virulence genes warrant future surveillance and preventative measures.