Microbial diversity in water and animal faeces: a metagenomic analysis to assess public health risk

Health impacts of poor water quality on an endangered shorebird breeding programme in Aotearoa New Zealand

CASE HISTORY

Two clusters of mortality among endangered tūturuatu/tchūriwat'/shore plover (Thinornis novaeseelandiae) have occurred at captive breeding facilities around New Zealand in recent years. In the first, four chicks died at Pūkaha National Wildlife Centre (Mount Bruce, NZ) in February 2016, and in the second five adult birds at the Cape Sanctuary (Cape Kidnappers, NZ) died in 2022.

CLINICAL FINDINGS

In 2016, four chicks were noted to become weak, have increased vocalisations and closed eyes prior to death. The remaining chicks were treated for 5 days with amoxycillin/clavulanate orally twice daily. Water containers and brooders were cleaned and disinfected with chlorhexidine. No further mortality was seen.In the 2022 cluster, three adult breeding birds died acutely and five others showed inappetence, weight loss and diarrhoea approximately 10 days after heavy rains flooded the local river. The five birds were treated with amoxycillin/clavulanate orally twice daily and oral fluids for 5 days. Two birds died and three survived. No breeding occurred in the aviaries in the following season.

PATHOLOGICAL FINDINGS

In 2016, the chicks showed pulmonary changes ranging from congestion and oedema to heterophilic inflammation consistent with septicaemia.In 2022, the adult birds showed proliferation of bacteria in the distal small intestine associated with mucosal ulceration and heterophilic infiltration. Acid-fast staining of the caecal contents in one bird showed organisms consistent with Cryptosporidium spp.

LABORATORY FINDINGS

Aerobic bacterial cultures of the lung and liver of two affected chicks carried out in 2016 showed heavy growth of Plesiomonas shigelloides. The same organism was cultured from water trays and holding tanks containing water boatmen (Sigara arguta) on which the chicks were fed.In 2022, cultures from the livers of three dead birds each showed a mixed bacterial growth with differing dominant organisms (Aeromonas sobria, Hafnia alvei, Citrobacter freundii and an Enterococcus sp.). PCR and sequencing confirmed Cryptosporidium parvum in the caecum of one bird. Fresh faeces from 24 breeding birds from the captive breeding facilities were negative by PCR for Cryptosporidium spp.The captive breeding facilities obtain water for the aviaries and aquatic invertebrates to feed to the chicks from local freshwater sources. Water quality testing at the Cape Sanctuary revealed concentrations of faecal indicator bacteria in excess of safe drinking water guidelines, with peaks following heavy rainfall.

CLINICAL RELEVANCE

Fluctuations in water quality associated with mammalian faecal bacteria can adversely affect bird health and impact on captive rearing of endangered wildlife.

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 Mākirikiri 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.

Omics-based ecosurveillance for the assessment of ecosystem function, health, and resilience

Current environmental monitoring efforts often focus on known, regulated contaminants ignoring the potential effects of unmeasured compounds and/or environmental factors. These specific, targeted approaches lack broader environmental information and understanding, hindering effective environmental management and policy. Switching to comprehensive, untargeted monitoring of contaminants, organism health, and environmental factors, such as nutrients, temperature, and pH, would provide more effective monitoring with a likely concomitant increase in environmental health. However, even this method would not capture subtle biochemical changes in organisms induced by chronic toxicant exposure. Ecosurveillance is the systematic collection, analysis, and interpretation of ecosystem health-related data that can address this knowledge gap and provide much-needed additional lines of evidence to environmental monitoring programs. Its use would therefore be of great benefit to environmental management and assessment. Unfortunately, the science of ‘ecosurveillance’, especially omics-based ecosurveillance is not well known. Here, we give an overview of this emerging area and show how it has been beneficially applied in a range of systems. We anticipate this review to be a starting point for further efforts to improve environmental monitoring via the integration of comprehensive chemical assessments and molecular biology-based approaches. Bringing multiple levels of omics technology-based assessment together into a systems-wide ecosurveillance approach will bring a greater understanding of the environment, particularly the microbial communities upon which we ultimately rely to remediate perturbed ecosystems.

Human Land-Use Impacts Viral Diversity and Abundance in a New Zealand River

Although water-borne viruses have important implications for the health of humans and other animals, little is known about the impact of human land use on viral diversity and evolution in water systems such as rivers. We used metatranscriptomic sequencing to compare the diversity and abundance of viruses at sampling sites along a single river in New Zealand that differed in human land-use impacts, ranging from pristine to urban. From this, we identified 504 putative virus species, of which 97 per cent were novel. Many of the novel viruses were highly divergent and likely included a new subfamily within the Parvoviridae. We identified at least sixty-three virus species that may infect vertebrates—most likely fish and water birds—from the Astroviridae, Birnaviridae, Parvoviridae, and Picornaviridae. No putative human viruses were detected. Importantly, we observed differences in the composition of viral communities at sites impacted by human land use (farming and urban) compared to native forest sites (pristine). At the viral species level, the urban sites had higher diversity (327 virus species) than the farming (n = 150) and pristine sites (n = 119), and more viruses were shared between the urban and farming sites (n = 76) than between the pristine and farming or urban sites (n = 24). The two farming sites had a lower viral abundance across all host types, while the pristine sites had a higher abundance of viruses associated with animals, plants, and fungi. We also identified viruses linked to agriculture and human impact at the river sampling sites in farming and urban areas that were not present at the native forest sites. Although based on a small sample size, our study suggests that human land use can impact viral communities in rivers, such that further work is needed to reduce the impact of intensive farming and urbanisation on water systems.

Investigation of the microbial community structure and diversity in the environment surrounding a veterinary antibiotic production factory

The ecological toxicity caused by antibiotic residues and resistance genes in the environment affects the community structures and activities of environmental microorganisms; the ecological toxicity effects of a long-term exposure to low doses antibiotic residues on microorganisms have not however been well-studied. In this work, sequence analysis and species annotation of the full-length 16S rRNA gene were carried out on the extracted whole genome by a 3-generation sequencing method to analyze the diversity of the microbial populations and the population differences among different sampling sites in the environment surrounding a veterinary antibiotic production factory. A total of 1720 OTUs (Operational Taxonomic Unit, OTU) were found, of which 1055 OTUs were in the river samples and 993 OTUs were in the soil samples. 643 and 438 bacterial strains were identified in the river water and soil samples respectively. The bacterial populations are classified into 29 phylum, 612 genus, and 849 species. The dominant phylum of bacteria was Proteobacteria, which was also the absolutely dominant phylum. Shannon diversity index of the bacteria showed that the bacterial abundance in downstream river was significantly higher than that in an upstream non-polluted area (P < 0.001), but the difference of bacterial abundance between soil samples was not significant. There were 61 biomarkers in the river water samples from different sampling points, and 14 biomarkers in soil samples. It was found by the difference statistics of the microbial community that there are multiple biomarkers between this veterinary drug production site and the upstream non-polluted area. Significant differences between multiple functional genes were also found in metabolic pathways of the microorganisms. A similar trend was found for the distribution of antibiotic resistance genes (ARGs). It is concluded that the population composition of microorganisms and diversity are likely related to antibiotic residues and to the distributions of ARGs in the environment surrounding the antibiotic production factory.

First report of novel assemblages and mixed infections of Giardia duodenalis in human isolates from New Zealand

Giardia duodenalis (syn. G. intestinalis and G. lamblia) is a protozoan parasite that cause disease (giardiasis) in humans and other animals. The pathogen is classified into eight assemblages, further divided into sub-assemblages, based on genetic divergence and host specificities. There are two zoonotic subtypes known as assemblages A and B, whilst assemblages from C to H are mainly found in domesticated animals, rodents and marine mammals. Here, we report for the first time the presence of assemblage E and sub-assemblage AIII in human isolates from the South Island in New Zealand. We identified a > 99% nucleotide similarity of assemblage E and sub-assemblage AIII with sequences of the gdh gene available in GenBank from individual human samples collected in Dunedin and Christchurch, respectively. We also performed a deep sequencing approach to assess intra-host assemblage variation. The sample from Dunedin showed evidence of mixed assemblage E and zoonotic sub-assemblage BIV. The report of two novel assemblages and mixed infections provides insights into the genetic diversity, epidemiology and transmission dynamics of Giardia duodenalis in New Zealand.