Browsing by Author "Wilkinson DA"
Now showing 1 - 14 of 14
Results Per Page
Sort Options
- ItemA cross-sectional investigation of Leptospira at the wildlife-livestock interface in New Zealand(PLOS, 2023-09-06) Moinet M; Oosterhof H; Nisa S; Haack N; Wilkinson DA; Aberdein D; Russell JC; Vallée E; Collins-Emerson J; Heuer C; Benschop J; Stevenson BThere has been a recent upsurge in human cases of leptospirosis in New Zealand, with wildlife a suspected emerging source, but up-to-date knowledge on this topic is lacking. We conducted a cross-sectional study in two farm environments to estimate Leptospira seroprevalence in wildlife and sympatric livestock, PCR/culture prevalence in wildlife, and compare seroprevalence and prevalence between species, sex, and age groups. Traps targeting house mice (Mus musculus), black rats (Rattus rattus), hedgehogs (Erinaceus europaeus) and brushtail possums (Trichosurus vulpecula) were set for 10 trap-nights in March-April 2017 on a dairy (A) and a beef and sheep (B) farm. Trapped wild animals and an age-stratified random sample of domestic animals, namely cattle, sheep and working dogs were blood sampled. Sera were tested by microagglutination test for five serogroups and titres compared using a Proportional Similarity Index (PSI). Wildlife kidneys were sampled for culture and qPCR targeting the lipL32 gene. True prevalence in mice was assessed using occupancy modelling by collating different laboratory results. Infection profiles varied by species, age group and farm. At the MAT cut-point of ≥ 48, up to 78% of wildlife species, and 16-99% of domestic animals were seropositive. Five of nine hedgehogs, 23/105 mice and 1/14 black rats reacted to L. borgpetersenii sv Ballum. The sera of 4/18 possums and 4/9 hedgehogs reacted to L. borgpetersenii sv Hardjobovis whilst 1/18 possums and 1/9 hedgehogs reacted to Tarassovi. In ruminants, seroprevalence for Hardjobovis and Pomona ranged 0-90% and 0-71% depending on the species and age group. Titres against Ballum, Tarassovi and Copenhageni were also observed in 4-20%, 0-25% and 0-21% of domestic species, respectively. The PSI indicated rodents and livestock had the most dissimilar serological responses. Three of nine hedgehogs, 31/105 mice and 2/14 rats were carrying leptospires (PCR and/or culture positive). True prevalence estimated by occupancy modelling in mice was 38% [95% Credible Interval 26, 51%] on Farm A and 22% [11, 40%] on Farm B. In the same environment, exposure to serovars found in wildlife species was commonly detected in livestock. Transmission pathways between and within species should be assessed to help in the development of efficient mitigation strategies against Leptospira.
- ItemDensity matters: How population dynamics of house mice (Mus musculus) inform the epidemiology of Leptospira(John Wiley & Sons Ltd on behalf of British Ecological Society, 2024-07-22) Moinet M; Abrahão CR; Gasparotto VPO; Wilkinson DA; Vallée E; Benschop J; Russell JC; Elderd B1. Rodents are maintenance hosts of numerous pathogens, and both their density and the pathogen prevalence determine the risk they pose to other animals or humans. However, density is often overlooked. We investigated a capture-mark-recapture-sampling strategy to study introduced mice (Mus musculus) and Leptospira as a model and demonstrate the advantages of a combined approach. 2. We estimated population density and Leptospira prevalence in mice in a replicated longitudinal survey conducted between 2016 and 2018. Capture-mark-recapture sessions were undertaken at two sites in Spring and Autumn and blood and kidney samples were collected at the end of each session. Mouse density and areas of activity were estimated using spatially explicit capture–recapture (SECR) models and both were compared between Leptospira positive and negative mice. Leptospira exposure and shedding status were estimated using Microscopic Agglutination Test, and a combination of culture and lipL32 PCR on kidneys. 3. Leptospira prevalence was higher in spring (83%–86%) than in autumn (31%–37%) and mouse densities simultaneously varied from 3.6 to 55.9/ha. However, despite these variations in prevalence and density, the density of infected animals remained relatively constant over time (3–8/ha). Shedding or being seropositive was also associated with the activity of mice. Shedding or seropositive mice had a larger activity area, and seropositive mice were trapped on average 1 day earlier than seronegative mice. 4. Synthesis and applications: Our results show how understanding the population dynamics of pathogen-carrying rodents is critical in epidemiology. The wider movement patterns and easier encounters of positive mice highlight the possibility of biases in classical prevalence surveys and have implications for disease transmission within and between species. Importantly, and quite counter-intuitively, Leptospira prevalence was negatively associated with mouse density, resulting in a constant density of shedders that contradicts the conventional view of higher exposure risk at high rodent density. More broadly, such sampling designs can improve animal and disease control policies and better inform modelling studies by providing more parameter estimates than classical prevalence surveys.
- ItemDiverse Epidemiology of Leptospira serovars Notified in New Zealand, 1999-2017(MDPI (Basel, Switzerland), 14/10/2020) Nisa S; Wilkinson DA; Angelin-Bonnet O; Paine S; Cullen K; Wight J; Baker MG; Benschop JLeptospirosis in New Zealand has been strongly associated with animal-contact occupations and with serovars Hardjo and Pomona. However, recent data suggest changes in these patterns, hence, serovar-specific epidemiology of leptospirosis from 1999 to 2017 was investigated. The 19-year average annual incidence is 2.01/100,000. Early (1999-2007) and late (2008-2017) study period comparisons showed a significant increase in notifications with serovar Ballum (IRR: 1.59, 95% CI: 1.22-2.09) in all cases and serovar Tarassovi (IRR: 1.75, 95% CI: 1.13-2.78) in Europeans and a decrease in notifications with serovars Hardjo and Pomona in all cases. Incidences of Ballum peaked in winter, Hardjo peaked in spring and Tarassovi peaked in summer. Incidence was highest in Māori (2.24/100,000) with dominant serovars being Hardjo and Pomona. Stratification by occupation showed meat workers had the highest incidence of Hardjo (57.29/100,000) and Pomona (45.32/100,000), farmers had the highest incidence of Ballum (11.09/100,000) and dairy farmers had the highest incidence of Tarassovi (12.59/100,000). Spatial analysis showed predominance of Hardjo and Pomona in Hawke's Bay, Ballum in West Coast and Northland and Tarassovi in Waikato, Taranaki and Northland. This study highlights the serovar-specific heterogeneity of human leptospirosis in New Zealand that should be considered when developing control and prevention strategies.
- ItemIdentification of pathogenic Leptospira species and serovars in New Zealand using metabarcoding.(2021) Wilkinson DA; Edwards M; Benschop J; Nisa SLeptospirosis is a zoonotic disease of global importance. The breadth of Leptospira diversity associated with both human and animal disease poses major logistical challenges to the use of classical diagnostic techniques, and increasingly molecular diagnostic tools are used for their detection. In New Zealand, this has resulted in an increase in positive cases reported nationally that have not been attributed to the infecting serovar or genomospecies. In this study, we used data from all pathogenic Leptospira genomes to identify a partial region of the glmU gene as a suitable locus for the discrimination of the infecting species and serovars of New Zealand-endemic Leptospira. This method can be used in culture and culture-independent scenarios making it flexible for diagnostics in humans, animals, and environmental samples. We explored the use of this locus as a molecular barcoding tool via the Oxford Nanopore Technology (ONT) sequencing platform MinION. Sequences obtained by this method allowed specific identification of Leptospira species in mixed and enriched environmental cultures, however read error inherent in the MinION sequencing system reduced the accuracy of strain/variant identification. Using this approach to characterise Leptospira in enriched environmental cultures, we detected the likely presence of Leptospira genomospecies that have not been reported in New Zealand to date. This included a strain of L. borgpetersenii that has recently been identified in dairy cattle and sequences similar to those of L. mayottensis. L. tipperaryensis, L. dzianensis and L. alstonii.
- ItemMolecular typing of Leptospira spp. in farmed and wild mammals reveals new host-serovar associations in New Zealand.(Taylor and Francis Group, 2024-01-01) Wilkinson DA; Edwards M; Shum C; Moinet M; Anderson NE; Benschop J; Nisa SAIMS: To apply molecular typing to DNA isolated from historical samples to determine Leptospira spp. infecting farmed and wild mammals in New Zealand. MATERIALS AND METHODS: DNA samples used in this study were extracted from urine, serum or kidney samples (or Leptospira spp. cultures isolated from them) collected between 2007 and 2017 from a range of domestic and wildlife mammalian species as part of different research projects at Massey University. Samples were included in the study if they met one of three criteria: samples that tested positive with a lipL32 PCR for pathogenic Leptospira; samples that tested negative by lipL32 PCR but were recorded as positive to PCR for pathogenic Leptospira in the previous studies; or samples that were PCR-negative in all studies but were from animals with positive agglutination titres against serogroup Tarassovi. DNA samples were typed using PCR that targeted either the glmU or gyrB genetic loci. The resulting amplicons were sequenced and typed relative to reference sequences. RESULTS: We identified several associations between mammalian hosts and Leptospira strains/serovars that had not been previously reported in New Zealand. Leptospira borgpetersenii strain Pacifica was found in farmed red deer (Cervus elaphus) samples, L. borgpetersenii serovars Balcanica and Ballum were found in wild red deer samples, Leptospira interrogans serovar Copenhageni was found in stoats (Mustela erminea) and brushtail possums (Trichosurus vulpecula), and L. borgpetersenii was found in a ferret (Mustela putorius furo). Furthermore, we reconfirmed previously described associations including dairy cattle with L. interrogans serovars Copenhageni and Pomona and L. borgpetersenii serovars Ballum, Hardjo type bovis and strain Pacifica, sheep with L. interrogans serovar Pomona and L. borgpetersenii serovar Hardjo type bovis, brushtail possum with L. borgpetersenii serovar Balcanica, farmed deer with L. borgpetersenii serovar Hardjo type bovis and hedgehogs (Erinaceus europaeus) with L. borgpetersenii serovar Ballum. CONCLUSIONS: This study provides an updated summary of host-Leptospira associations in New Zealand and highlights the importance of molecular typing. Furthermore, strain Pacifica, which was first identified as Tarassovi using serological methods in dairy cattle in 2016, has circulated in animal communities since at least 2007 but remained undetected as serology is unable to distinguish the different genotypes. CLINICAL RELEVANCE: To date, leptospirosis in New Zealand has been diagnosed with serological typing, which is deficient in typing all strains in circulation. Molecular methods are necessary to accurately type strains of Leptospira spp. infecting mammals in New Zealand.
- ItemOf Mice, Cattle, and Men: A Review of the Eco-Epidemiology of Leptospira borgpetersenii Serovar Ballum.(20/10/2021) Moinet M; Wilkinson DA; Aberdein D; Russell JC; Vallée E; Collins-Emerson JM; Heuer C; Benschop JIn New Zealand (NZ), leptospirosis is a mostly occupational zoonosis, with >66% of the recently notified cases being farm or abattoir workers. Livestock species independently maintain Leptospira borgpetersenii serovar Hardjo and L. interrogans serovar Pomona, and both are included in livestock vaccines. The increasing importance in human cases of Ballum, a serovar associated with wildlife, suggests that wildlife may be an overlooked source of infection. Livestock could also act as bridge hosts for humans. Drawing from disease ecology frameworks, we chose five barriers to include in this review based on the hypothesis that cattle act as bridge hosts for Ballum. Using a narrative methodology, we collated published studies pertaining to (a) the distribution and abundance of potential wild maintenance hosts of Ballum, (b) the infection dynamics (prevalence and pathogenesis) in those same hosts, (c) Ballum shedding and survival in the environment, (d) the exposure and competency of cattle as a potential bridge host, and (e) exposure for humans as a target host of Ballum. Mice (Mus musculus), rats (Rattus rattus, R. norvegicus) and hedgehogs (Erinaceus europaeus) were suspected as maintenance hosts of Ballum in NZ in studies conducted in the 1970s-1980s. These introduced species are distributed throughout NZ, and are present on pastures. The role of other wildlife in Ballum (and more broadly Leptospira) transmission remains poorly defined, and has not been thoroughly investigated in NZ. The experimental and natural Ballum infection of cattle suggest a low pathogenicity and the possibility of shedding. The seroprevalence in cattle appears higher in recent serosurveys (3 to 14%) compared with studies from the 1970s (0 to 3%). This review identifies gaps in the knowledge of Ballum, and highlights cattle as a potential spillover host. Further studies are required to ascertain the role that wild and domestic species may play in the eco-epidemiology of Ballum in order to understand its survival in the environment, and to inform control strategies.
- ItemPhylogenomic Characterization of a Novel Corynebacterium Species Associated with Fatal Diphtheritic Stomatitis in Endangered Yellow-Eyed Penguins(American Society for Microbiology, 2021-06-08) Saunderson SC; Nouioui I; Midwinter AC; Wilkinson DA; Young MJ; McInnes KM; Watts J; Sangal V; Segata NYellow-eyed penguins, Megadyptes antipodes, are an endangered species that are endemic to New Zealand. Outbreaks of diphtheritic stomatitis have caused significant mortality for this species, especially among young chicks. In this study, we isolated 16 Corynebacterium sp. isolates from the oral cavities of 2- to 14-day-old chicks at a range of infection stages and sequenced the genomes to understand their virulence mechanisms. Phylogenomic and matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) characterization indicate that these strains belong to a novel Corynebacterium species. A simple multiplex PCR-based diagnostic assay has been developed to identify these strains rapidly and reliably. Similar to other corynebacteria, genomic islands and prophages introduced significant diversity among these strains that has potentially led to minor functional variations between the two lineages. Despite the presence of multiple corynebacterial virulence genes and a spaDEF-type pilus gene cluster among these strains, the survival rate was much higher in Galleria mellonella larvae than in those inoculated with Corynebacterium ulcerans NZRM 818 and Corynebacterium pseudotuberculosis NZRM 3004. Therefore, these strains are opportunistic pathogens causing high mortality among young penguin chicks due to a less-developed immune system. IMPORTANCE Yellow-eyed penguins, Megadyptes antipodes, are endangered species with a sharp decline in the numbers of breeding pairs over the last 2 decades. Diphtheritic stomatitis, characterized by a thick fibrinopurulent exudate in the oral cavities and symptoms, including inanition and significant weight loss, is responsible for significant mortality among the young chicks. These chicks are treated with antibiotics, amoxicillin-clavulanic acid or enrofloxacin, but do not always recover from the infection. The pathogen causing these infections and the mechanism of pathogenesis are unclear. This study has identified a novel Corynebacterium species to be associated with diphtheritic stomatitis in yellow-eyed penguins with potential virulence genes that are likely involved in pathogenesis. Importantly, a gene encoding an exotoxin, phospholipase D, is present among these strains. The inactivated form of this enzyme could potentially be used as an effective vaccine to protect these penguins from infection.
- ItemPopulation Structure and Antimicrobial Resistance in Campylobacter jejuni and C. coli Isolated from Humans with Diarrhea and from Poultry, East Africa.(Centers for Disease Control and Prevention, 2024-10) French NP; Thomas KM; Amani NB; Benschop J; Bigogo GM; Cleaveland S; Fayaz A; Hugho EA; Karimuribo ED; Kasagama E; Maganga R; Melubo ML; Midwinter AC; Mmbaga BT; Mosha VV; Mshana FI; Munyua P; Ochieng JB; Rogers L; Sindiyo E; Swai ES; Verani JR; Widdowson M-A; Wilkinson DA; Kazwala RR; Crump JA; Zadoks RNCampylobacteriosis and antimicrobial resistance (AMR) are global public health concerns. Africa is estimated to have the world's highest incidence of campylobacteriosis and a relatively high prevalence of AMR in Campylobacter spp. from humans and animals. Few studies have compared Campylobacter spp. isolated from humans and poultry in Africa using whole-genome sequencing and antimicrobial susceptibility testing. We explored the population structure and AMR of 178 Campylobacter isolates from East Africa, 81 from patients with diarrhea in Kenya and 97 from 56 poultry samples in Tanzania, collected during 2006-2017. Sequence type diversity was high in both poultry and human isolates, with some sequence types in common. The estimated prevalence of multidrug resistance, defined as resistance to >3 antimicrobial classes, was higher in poultry isolates (40.9%, 95% credible interval 23.6%-59.4%) than in human isolates (2.5%, 95% credible interval 0.3%-6.8%), underlining the importance of antimicrobial stewardship in livestock systems.
- ItemThe Genetic Relatedness and Antimicrobial Resistance Patterns of Mastitis-Causing Staphylococcus aureus Strains Isolated from New Zealand Dairy Cattle(MDPI (Basel, Switzerland), 2021-11-22) Greening SS; Zhang J; Midwinter AC; Wilkinson DA; McDougall S; Gates MC; French NP; Butaye PStaphylococcus aureus is one of the leading causes of bovine mastitis worldwide and is a common indication for use of antimicrobials on dairy farms. This study aims to investigate the association between on-farm antimicrobial usage and the antimicrobial resistance (AMR) profiles of mastitis-causing S. aureus. Whole-genome sequencing was performed on 57 S. aureus isolates derived from cows with either clinical or subclinical mastitis from 17 dairy herds in New Zealand. The genetic relatedness between isolates was examined using the core single nucleotide polymorphism alignment whilst AMR and virulence genes were identified in-silico. The association between gene presence-absence and sequence type (ST), antimicrobial susceptibility and dry cow therapy treatment was investigated using Scoary. Altogether, eight STs were identified with 61.4% (35/57) belonging to ST-1. Furthermore, 14 AMR-associated genes and 76 virulence-associated genes were identified, with little genetic diversity between isolates belonging to the same ST. Several genes including merR1 which is thought to play a role in ciprofloxacin-resistance were found to be significantly overrepresented in isolates sampled from herds using ampicillin/cloxacillin dry cow therapy. Overall, the presence of resistance genes remains low and current antimicrobial usage patterns do not appear to be driving AMR in S. aureus associated with bovine mastitis.
- ItemThe role of seabirds of the Iles Eparses as reservoirs and disseminators of parasites and pathogens(Elsevier, 2015-04) McCoy KD; Dietrich M; Jaeger A; Wilkinson DA; Bastien M; Lagadec E; Boulinier T; Pascalis H; Tortosa P; Le Corre M; othersThe role of birds as reservoirs and disseminators of parasites and pathogens has received much attention over the past several years due to their high vagility. Seabirds are particularly interesting hosts in this respect. In addition to incredible long-distance movements during migration, foraging and prospecting, these birds are long-lived, site faithful and breed in dense aggregations in specific colony locations. These different characteristics can favor both the local maintenance and large-scale dissemination of parasites and pathogens. The Iles Eparses provide breeding and feeding grounds for more than 3 million breeding pairs of seabirds including at least 13 species. Breeding colonies on these islands are relatively undisturbed by human activities and represent natural metapopulations in which seabird population dynamics, movement and dispersal can be studied in relation to that of circulating parasites and pathogens. In this review, we summarize previous knowledge and recently-acquired data on the parasites and pathogens found in association with seabirds of the Iles Eparses. These studies have revealed the presence of a rich diversity of infectious agents (viruses, bacteria and parasites) carried by the birds and/or their local ectoparasites (ticks and louse flies). Many of these agents are widespread and found in other ecosystems confirming a role for seabirds in their large scale dissemination and maintenance. The heterogeneous distribution of parasites and infectious agents among islands and seabird species suggests that relatively independent metacommunities of interacting species may exist within the western Indian Ocean. In this context, we discuss how the patterns and determinants of seabird movements may alter parasite and pathogen circulation. We conclude by outlining key aspects for future research given the baseline data now available and current concerns in eco-epidemiology and biodiversity conservation.
- ItemTransmission dynamics of an antimicrobial resistant Campylobacter jejuni lineage in New Zealand’s commercial poultry network(Elsevier B.V, 2021-12) Greening SS; Zhang J; Midwinter AC; Wilkinson DA; Fayaz A; Williamson DA; Anderson MJ; Gates MC; French NPUnderstanding the relative contribution of different between-farm transmission pathways is essential in guiding recommendations for mitigating disease spread. This study investigated the association between contact pathways linking poultry farms in New Zealand and the genetic relatedness of antimicrobial resistant Campylobacter jejuni Sequence Type 6964 (ST-6964), with the aim of identifying the most likely contact pathways that contributed to its rapid spread across the industry. Whole-genome sequencing was performed on 167C. jejuni ST-6964 isolates sampled from across 30 New Zealand commercial poultry enterprises. The genetic relatedness between isolates was determined using whole genome multilocus sequence typing (wgMLST). Permutational multivariate analysis of variance and distance-based linear models were used to explore the strength of the relationship between pairwise genetic associations among the C. jejuni isolates and each of several pairwise distance matrices, indicating either the geographical distance between farms or the network distance of transportation vehicles. Overall, a significant association was found between the pairwise genetic relatedness of the C. jejuni isolates and the parent company, the road distance and the network distance of transporting feed vehicles. This result suggests that the transportation of feed within the commercial poultry industry as well as other local contacts between flocks, such as the movements of personnel, may have played a significant role in the spread of C. jejuni. However, further information on the historical contact patterns between farms is needed to fully characterise the risk of these pathways and to understand how they could be targeted to reduce the spread of C. jejuni.
- ItemUpdating the genomic taxonomy and epidemiology of Campylobacter hyointestinalis(Springer Nature Limited, 2018-02-05) Wilkinson DA; O'Donnell AJ; Akhter RN; Fayaz A; Mack HJ; Rogers LE; Biggs PJ; French NP; Midwinter ACCampylobacter hyointestinalis is a member of an emerging group of zoonotic Campylobacter spp. that are increasingly identified in both gastric and non-gastric disease in humans. Here, we discovered C. hyointestinalis in three separate classes of New Zealand ruminant livestock; cattle, sheep and deer. To investigate the relevance of these findings we performed a systematic literature review on global C. hyointestinalis epidemiology and used comparative genomics to better understand and classify members of the species. We found that C. hyointestinalis subspecies hyointestinalis has an open pangenome, with accessory gene contents involved in many essential processes such as metabolism, virulence and defence. We observed that horizontal gene transfer is likely to have played an overwhelming role in species diversification, favouring a public-goods-like mechanism of gene ‘acquisition and resampling’ over a tree-of-life-like vertical inheritance model of evolution. As a result, simplistic gene-based inferences of taxonomy by similarity are likely to be misleading. Such genomic plasticity will also mean that local evolutionary histories likely influence key species characteristics, such as host-association and virulence. This may help explain geographical differences in reported C. hyointestinalis epidemiology and limits what characteristics may be generalised, requiring further genomic studies of C. hyointestinalis in areas where it causes disease.
- ItemUsing drivers and transmission pathways to identify SARS-like coronavirus spillover risk hotspots.(Springer Nature Limited, 2023-10-27) Muylaert RL; Wilkinson DA; Kingston T; D'Odorico P; Rulli MC; Galli N; John RS; Alviola P; Hayman DTSThe emergence of SARS-like coronaviruses is a multi-stage process from wildlife reservoirs to people. Here we characterize multiple drivers-landscape change, host distribution, and human exposure-associated with the risk of spillover of zoonotic SARS-like coronaviruses to help inform surveillance and mitigation activities. We consider direct and indirect transmission pathways by modeling four scenarios with livestock and mammalian wildlife as potential and known reservoirs before examining how access to healthcare varies within clusters and scenarios. We found 19 clusters with differing risk factor contributions within a single country (N = 9) or transboundary (N = 10). High-risk areas were mainly closer (11-20%) rather than far ( < 1%) from healthcare. Areas far from healthcare reveal healthcare access inequalities, especially Scenario 3, which includes wild mammals and not livestock as secondary hosts. China (N = 2) and Indonesia (N = 1) had clusters with the highest risk. Our findings can help stakeholders in land use planning, integrating healthcare implementation and One Health actions.
- ItemWhole-Genome Sequencing and Virulome Analysis of Escherichia coli Isolated from New Zealand Environments of Contrasting Observed Land Use(American Society for Microbiology, 2022-05-10) Cookson AL; Marshall JC; Biggs PJ; Rogers LE; Collis RM; Devane M; Stott R; Wilkinson DA; Kamke J; Brightwell G; Elkins CAGeneric Escherichia coli is commonly used as an indicator of fecal contamination to assess water quality and human health risk. Where measured E. coli exceedances occur, the presence of other pathogenic microorganisms, such as Shiga toxin-producing E. coli (STEC), is assumed, but confirmatory data are lacking. Putative E. coli isolates (n = 709) were isolated from water, sediment, soil, periphyton, and feces samples (n = 189) from five sites representing native forest and agricultural environments. Ten E. coli isolates (1.41%) were stx2 positive, 19 (2.7%) were eae positive, and stx1-positive isolates were absent. At the sample level, stx2-positive E. coli (5 of 189, 2.6%) and eae-positive isolates (16 of 189, 8.5%) were rare. Using real-time PCR, these STEC-associated virulence factors were determined to be more prevalent in sample enrichments (stx1, 23.9%; stx2, 31.4%; eae, 53.7%) and positively correlated with generic E. coli isolate numbers (P < 0.05) determined using culture-based methods. Whole-genome sequencing (WGS) was undertaken on a subset of 238 isolates with assemblies representing seven E. coli phylogroups (A, B1, B2, C, D, E, and F), 22 Escherichia marmotae isolates, and 1 Escherichia ruysiae isolate. Virulence factors, including those from extraintestinal pathogenic E. coli, were extremely diverse in isolates from the different locations and were more common in phylogroup B2. Analysis of the virulome from WGS data permitted the identification of gene repertoires that may be involved in environmental fitness and broadly align with phylogroup. Although recovery of STEC isolates was low, our molecular data indicate that they are likely to be widely present in environmental samples containing diverse E. coli phylogroups. IMPORTANCE This study takes a systematic sampling approach to assess the public health risk of Escherichia coli recovered from freshwater sites within forest and farmland. The New Zealand landscape is dominated by livestock farming, and previous work has demonstrated that "recreational exposure to water" is a risk factor for human infection by Shiga toxin-producing Escherichia coli (STEC). Though STEC isolates were rarely isolated from water samples, STEC-associated virulence factors were identified more commonly from water sample culture enrichments and were associated with increased generic E. coli concentrations. Whole-genome sequencing data from both E. coli and newly described Escherichia spp. demonstrated the presence of virulence factors from E. coli pathotypes, including extraintestinal pathogenic E. coli. This has significance for understanding and interpreting the potential health risk from E. coli where water quality is poor and suggests a role of virulence factors in survival and persistence of E. coli and Escherichia spp.