Development of a multiplex droplet digital PCR assay for simultaneous detection and quantification of Escherichia coli, E. marmotae, and E. ruysiae in water samples.

Escherichia coli are widely used by water quality managers as Fecal Indicator Bacteria, but current quantification methods do not differentiate them from benign, environmental Escherichia species such as E. marmotae (formerly named cryptic clade V) or E. ruysiae (cryptic clades III and IV). Reliable and specific techniques for their identification are required to avoid confounding microbial water quality assessments. To address this, a multiplex droplet digital PCR (ddPCR) assay targeting lipB (E. coli and E. ruysiae) and bglC (E. marmotae) was designed. The ddPCR performance was assessed using in silico analysis; genomic DNA from 40 local, international, and reference strains of target and non-target coliforms; and spiked water samples in a range relevant to water quality managers (1 to 1000 cells/100 mL). Results were compared to an analogous quantitative PCR (qPCR) and the Colilert method. Both PCR assays showed excellent sensitivity with a limit of detection of 0.05 pg/µL and 0.005 pg/µl for ddPCR and qPCR respectively, and of quantification of 0.5 pg/µL of genomic DNA. The ddPCR allowed differentiation and quantification of three Escherichia species per run by amplitude multiplexing and showed a high concordance with concentrations measured by Colilert once proportional bias was accounted for. In silico specificity testing underlined the possibility to further detect and distinguish Escherichia cryptic clade VI. Finally, the applicability of the ddPCR was successfully tested on environmental water samples where E. marmotae and E. ruysiae potentially confound E. coli counts based on the Most Probable Number method, highlighting the utility of this novel ddPCR as an efficient and rapid discriminatory test to improve water quality assessments.

Draft genome sequences of Escherichia spp. isolates from New Zealand environmental sources.

Escherichia coli is often used as a fecal indicator bacterium for water quality monitoring. We report the draft genome sequences of 500 Escherichia isolates including newly described Escherichia species, namely Escherichia marmotae, Escherichia ruysiae, and Escherichia whittamii, obtained from diverse environmental sources to assist with improved public health risk assessments.

High-resolution genomic analysis to investigate the impact of the invasive brushtail possum (Trichosurus vulpecula) and other wildlife on microbial water quality assessments.

Escherichia coli are routine indicators of fecal contamination in water quality assessments. Contrary to livestock and human activities, brushtail possums (Trichosurus vulpecula), common invasive marsupials in Aotearoa/New Zealand, have not been thoroughly studied as a source of fecal contamination in freshwater. To investigate their potential role, Escherichia spp. isolates (n = 420) were recovered from possum gut contents and feces and were compared to those from water, soil, sediment, and periphyton samples, and from birds and other introduced mammals collected within the Makirikiri Reserve, Dannevirke. Isolates were characterized using E. coli-specific real-time PCR targeting the uidA gene, Sanger sequencing of a partial gnd PCR product to generate a gnd sequence type (gST), and for 101 isolates, whole genome sequencing. Escherichia populations from 106 animal and environmental sample enrichments were analyzed using gnd metabarcoding. The alpha diversity of Escherichia gSTs was significantly lower in possums and animals compared with aquatic environmental samples, and some gSTs were shared between sample types, e.g., gST535 (in 85% of samples) and gST258 (71%). Forty percent of isolates gnd-typed and 75% of reads obtained by metabarcoding had gSTs shared between possums, other animals, and the environment. Core-genome single nucleotide polymorphism (SNP) analysis showed limited variation between several animal and environmental isolates (<10 SNPs). Our data show at an unprecedented scale that Escherichia clones are shared between possums, other wildlife, water, and the wider environment. These findings support the potential role of possums as contributors to fecal contamination in Aotearoa/New Zealand freshwater. Our study deepens the current knowledge of Escherichia populations in under-sampled wildlife. It presents a successful application of high-resolution genomic methods for fecal source tracking, thereby broadening the analytical toolbox available to water quality managers. Phylogenetic analysis of isolates and profiling of Escherichia populations provided useful information on the source(s) of fecal contamination and suggest that comprehensive invasive species management strategies may assist in restoring not only ecosystem health but also water health where microbial water quality is compromised.

A cross-sectional investigation of Leptospira at the wildlife-livestock interface in New Zealand.

There 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.

Whole genome sequence analysis of ESBL-producing Escherichia coli recovered from New Zealand freshwater sites.

Extended-spectrum beta lactamase (ESBL)-producing Escherichia coli are often isolated from humans with urinary tract infections and may display a multidrug-resistant phenotype. These pathogens represent a target for a One Health surveillance approach to investigate transmission between humans, animals and the environment. This study examines the multidrug-resistant phenotype and whole genome sequence data of four ESBL-producing E. coli isolated from freshwater in New Zealand. All four isolates were obtained from a catchment with a mixed urban and pastoral farming land-use. Three isolates were sequence type (ST) 131 (CTX-M-27-positive) and the other ST69 (CTX-M-15-positive); a phylogenetic comparison with other locally isolated strains demonstrated a close relationship with New Zealand clinical isolates. Genes associated with resistance to antifolates, tetracyclines, aminoglycosides and macrolides were identified in all four isolates, together with fluoroquinolone resistance in two isolates. The ST69 isolate harboured the bla CTX-M-15 gene on a IncHI2A plasmid, and two of the three ST131 isolates harboured the bla CTX-M-27 genes on IncF plasmids. The last ST131 isolate harboured bla CTX-M-27 on the chromosome in a unique site between gspC and gspD. These data highlight a probable human origin of the isolates with subsequent transmission from urban centres through wastewater to the wider environment.

Of Mice, Cattle, and Men: A Review of the Eco-Epidemiology of Leptospira borgpetersenii Serovar Ballum.

In 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.

Disappearance of TBEV Circulation among Rodents in a Natural Focus in Alsace, Eastern France.

Tick-borne encephalitis virus (TBEV) depends mainly on a fragile mode of transmission, the co-feeding between infected nymphs and larvae on rodents, and thus persists under a limited set of biotic and abiotic conditions. If these conditions change, natural TBEV foci might be unstable over time. We conducted a longitudinal study over seven years in a mountain forest in Alsace, Eastern France, located at the western border of known TBEV distribution. The objectives were (i) to monitor the persistence of TBEV circulation between small mammals and ticks and (ii) to discuss the presence of TBEV circulation in relation to the synchronous activity of larvae and nymphs, to the densities of questing nymphs and small mammals, and to potential changes in meteorological conditions and deer densities. Small mammals were trapped five times per year from 2012 to 2018 to collect blood samples and record the presence of feeding ticks, and were then released. Questing nymphs were collected twice a year. Overall, 1344 different small mammals (Myodes glareolus and Apodemus flavicollis) were captured and 2031 serum samples were tested for the presence of antibodies against TBEV using an in-house ELISA. Seropositive rodents (2.1%) were only found from 2012 to 2015, suggesting that the virus disappeared afterwards. In parallel, we observed unusual variations in inter-annual nymph abundance and intra-annual larval activity that could be related to exceptional meteorological conditions. Changes in the densities of questing nymphs and deer associated with the natural stochastic variations in the frequency of contacts between rodents and infected ticks may have contributed to the endemic fadeout of TBEV on the study site. Further studies are needed to assess whether such events occur relatively frequently in the area, which could explain the low human incidence of TBE in Alsace and even in other areas of France.

Tick-Borne Encephalitis Virus: Seasonal and Annual Variation of Epidemiological Parameters Related to Nymph-to-Larva Transmission and Exposure of Small Mammals.

A greater knowledge of the ecology of the natural foci of tick-borne encephalitis virus (TBEV) is essential to better assess the temporal variations of the risk of tick-borne encephalitis for humans. To describe the seasonal and inter-annual variations of the TBEV-cycle and the epidemiological parameters related to TBEV nymph-to-larva transmission, exposure of small mammals to TBEV, and tick aggregation on small mammals, a longitudinal survey in ticks and small mammals was conducted over a 3-year period in a mountain forest in Alsace, eastern France. TBEV prevalence in questing nymphs was lower in 2013 than in 2012 and 2014, probably because small mammals (Myodes glareolus and Apodemus flavicollis) were more abundant in 2012, which reduced tick aggregation and co-feeding transmission between ticks. The prevalence of TBEV in questing nymphs was higher in autumn than spring. Despite these variations in prevalence, the density of infected questing nymphs was constant over time, leading to a constant risk for humans. The seroprevalence of small mammals was also constant over time, although the proportion of rodents infested with ticks varied between years and seasons. Our results draw attention to the importance of considering the complex relationship between small mammal densities, tick aggregation on small mammals, density of infected questing nymphs, and prevalence of infected nymphs in order to forecast the risk of TBEV for humans.

Longitudinal survey of two serotine bat (Eptesicus serotinus) maternity colonies exposed to EBLV-1 (European Bat Lyssavirus type 1): Assessment of survival and serological status variations using capture-recapture models.

This study describes two longitudinal serological surveys of European Bat Lyssavirus type 1 (EBLV-1) antibodies in serotine bat (Eptesicus serotinus) maternity colonies located in the North-East of France. This species is currently considered as the main EBLV-1 reservoir. Multievent capture-recapture models were used to determine the factors influencing bat rabies transmission as this method accounts for imperfect detection and uncertainty in disease states. Considering the period of study, analyses revealed that survival and recapture probabilities were not affected by the serological status of individuals, confirming the capacity of bats to be exposed to lyssaviruses without dying. Five bats have been found with EBLV-1 RNA in the saliva at the start of the study, suggesting they were caught during virus excretion period. Among these bats, one was interestingly recaptured one year later and harbored a seropositive status. Along the survey, some others bats have been observed to both seroconvert (i.e. move from a negative to a positive serological status) and serorevert (i.e. move from a positive to a negative serological status). Peak of seroprevalence reached 34% and 70% in site A and B respectively. On one of the 2 sites, global decrease of seroprevalence was observed all along the study period nuanced by oscillation intervals of approximately 2-3 years supporting the oscillation infection dynamics hypothesized during a previous EBLV-1 study in a Myotis myotis colony. Seroprevalence were affected by significantly higher seroprevalence in summer than in spring. The maximum time observed between successive positive serological statuses of a bat demonstrated the potential persistence of neutralizing antibodies for at least 4 years. At last, EBLV-1 serological status transitions have been shown driven by age category with higher seroreversion frequencies in adults than in juvenile. Juveniles and female adults seemed indeed acting as distinct drivers of the rabies virus dynamics, hypothesis have been addressed but their exact role in the EBLV-1 transmission still need to be specified.

Spatio-temporal dynamics of tularemia in French wildlife: 2002-2013.

Tularemia, caused by Francisella tularensis, is endemic in France. The surveillance of this disease in wildlife is operated by the SAGIR Network and by the National Reference Laboratory for Tularemia. Wild animals found dead or dying collected by the SAGIR network are necropsied and when tularemia is suspected culture and/or PCR are performed to confirm the diagnosis. The aim of this study was to present the results of tularemia surveillance in wildlife and to investigate the spatial and temporal pattern of tularemia observed between the 2002-2003 and 2012-2013 hunting seasons in French wildlife. Fourty-one to 121 cases were collected each hunting season for a total of 693 confirmed cases and 46 additional suspected cases. The main species affected was the European Brown Hare (Lepus europaeus) but 4 rabbits (Oryctolagus cuniculus), 2 roe deer (Capreolus capreolus) and one wild boar (Sus scrofa) were also found positive. The Standard Mortality Ratio was mapped and Kulldorff's retrospective space-time scan statistic was implemented and allowed the detection of several clusters: the nationwide outbreak of 2007-2008; 2 clusters in northern and central-western France in high hare-abundance areas and another in North-eastern France where the abundance of hares is low. Our results confirm the usefulness of brown hare as a sentinel of environmental risk.

Toxoplasmosis in natural populations of ungulates in France: prevalence and spatiotemporal variations.

Toxoplasmosis is characterized by a complex epidemiology. The risk of infection for humans depends on their contact with infective oocysts in a contaminated environment and on the amount of tissue cysts located within consumed meat. Unfortunately, the prevalence of tissue cysts is largely unknown for game species. Although herbivorous game species are a source of infection for humans, the level of infection found in wildlife can also be used to estimate environmental contamination. The aim of this study was to estimate the prevalence of Toxoplasma gondii infection and analyze its temporal dynamics in one population of chamois (Rupicapra rupicapra), one of mouflon (Ovis gmelini musimon), and two of roe deer (Capreolus capreolus) in France, surveyed during a period of 6 to 28 years. Taking into account individual risk factors, we specifically analyzed the relationship between T. gondii prevalence and meteorological conditions that may influence oocyst survival. Serum samples from 101 chamois, 143 mouflons, and 1155 roe deer were tested for antibodies against T. gondii using the modified agglutination test (MAT), an enzyme-linked immunosorbent assay (ELISA) assay, or both. Using MAT with a threshold of 1:6, seroprevalence was 14.7% in mouflon, 16.8% in chamois, and 43.7% in roe deer. In mouflon and roe deer, seroprevalence was positively correlated with age and/or body mass, in accordance with the hypothesis that antibodies have long-term persistence. In roe deer, seropositivity differed between the two populations and changed linearly over time between 1983 and 2010, increasing by a factor 1.75 every 10 years. Moreover, in this species, the highest prevalences were found during dry and cold years or during warm and moist years, depending on the population. Our results suggest that the risk for people to acquire infection through game meat increases over time, but with high variability according to the population of origin and meteorological conditions of the year.

Isolation of Bokeloh bat lyssavirus in Myotis nattereri in France.

Bokeloh bat lyssavirus (BBLV) was found in Myotis nattereri for the first time in northeastern France in July 2012. The complete genome sequence of the virus from the infected Natterer's bat was determined by whole-genome sequencing and compared to that of the first BBLV strain isolated in 2010 in Germany and with those of all currently identified lyssaviruses. The French isolate [KC169985] showed 98.7 % nucleotide sequence identity to the German BBLV strain [JF311903]. Several organs of the infected French bat were examined by classical rabies diagnostic methods: fluorescent antibody test, cell culture inoculation test and RT-qPCR. Antigen, infectious virus and high viral RNA levels were found in both the brain and salivary glands. Traces of genomic RNA were detected in the bladder, kidney and lung tissue. The results of an investigation of the distribution of lyssaviruses with the detection of infectious virus in the salivary glands suggest a possible mode of transmission of the virus.

Active surveillance of bat rabies in France: a 5-year study (2004-2009).

Active surveillance of bats in France started in 2004 with an analysis of 18 of the 45 bat species reported in Europe. Rabies antibodies were detected in six indigenous species, mainly in Eptesicus serotinus and Myotis myotis, suggesting previous contact with the EBLV-1 rabies virus. Nineteen of the 177 tested bats were shown serologically positive in seven sites, particularly in central and south-western France. Neither infectious viral particles nor viral genomes were detected in 173 and 308 tested oral swabs, respectively. The presence of neutralising antibodies in female bats (18.6%) was significantly higher than in males (5.6%).

Leptospirosis in free-ranging endangered European mink (Mustela lutreola) and other small carnivores (Mustelidae, Viverridae) from southwestern France.

To study the possible role of disease in the decline of endangered European mink (Mustela lutreola), we conducted a survey of antibody prevalence and renal carriage of pathogenic leptospira (Leptospira interrogans sensu lato) using serum and kidney samples collected from 1990 to 2007 from several free-ranging small carnivores and farmed American mink (Mustela vison) in southwestern France. An indirect microscopic agglutination test using a panel of 16 serovars belonging to 6 serogroups (Australis, Autumnalis, Icterohæmorrhagiæ, Grippotyphosa, Panama, Sejroe) revealed antibodies in all species, with significant differences in antibody prevalences: 74% in European mink (n=99), 65.4% in European polecats (Mustela putorius, n=133), 86% in American mink (n=74), 89% in stone martens (Martes foina, n=19), 74% in pine martens (Martes martes, n=19), 35% in common genets (Genetta genetta, n=79), and 31% in farmed American mink (n=51). Serogroups Australis and Icterohæmorragiæ were dominant in most free-ranging species; serogroup Grippotyphosa had high prevalences in European mink. Such high antibody prevalences have never been reported. They are probably related to the large number of known reservoirs, rats (Rattus spp.), muskrat (Ondatra zibethicus), and coypu (Myocastor coypu), in the study area. The polymerase chain reaction test specific for pathogenic leptospiral DNA detected renal carriage in 23% of 34 European mink, 22% of 18 polecats, and 15% of 33 free-ranging American mink, with no significant differences. Renal carriage shows that mustelids may shed leptospira for short periods, but their epidemiologic role is probably limited. High antibody prevalences suggest that the disease is unlikely to be highly pathogenic for these species; however, chronic forms of the disease (abortions, renal lesions) could reduce the reproductive success or life span of infected animals. Further studies on the pathogenicity of leptospirosis in these populations are needed to measure its impact on the population dynamics of these rodent predators.