Bacterial diversity and resistome analysis of drinking water stored in cisterns from two First Nations communities in Manitoba, Canada.

The microbiological content of water is an ongoing concern in First Nations communities in Canada. Many communities lack water treatment plants and continue to be under drinking water advisories. However, lack of access to treatment plants is only a part of the problem as poor water distribution systems also contribute to the failure to provide safe drinking water. Here, we studied the microbial diversity and antibiotic resistome from water stored in cisterns from two First Nations communities in Manitoba, Canada. We found that the cistern water contained a high number of bacteria and showed the presence of diverse antimicrobial resistance genes. Interestingly, the bacterial diversity and antimicrobial resistance genes varied considerably from that of the untreated source water, indicating that the origin of contamination in the cistern water came from within the treatment plant or along the delivery route to the homes. Our study highlights the importance of proper maintenance of the water distribution system in addition to access to water treatment facilities to ensure a supply of safe water to First Nations communities in Canada.IMPORTANCEThe work described addresses a critical issue in First Nations communities in Canada-the microbiological content of water. Many of these communities lack access to water treatment plants and frequently experience drinking water advisories. This study focused on the microbial diversity and antibiotic resistome in water stored in cisterns within two First Nations communities in Manitoba, Canada. These findings reveal that cistern water, a common source of drinking water in these communities, contains a high number of bacteria and a wide range of antimicrobial resistance genes. This highlights a serious health risk as exposure to such water can lead to the spread of drug-resistant infections, posing a threat to the well-being of the residents.

Concentration of Virus Particles from Environmental Water and Wastewater Samples Using Skimmed Milk Flocculation and Ultrafiltration.

Water and wastewater-based epidemiology have emerged as alternative methods to monitor and predict the course of outbreaks in communities. The recovery of microbial fractions, including viruses, bacteria, and microeukaryotes from wastewater and environmental water samples is one of the challenging steps in these approaches. In this study, we focused on the recovery efficiency of sequential ultrafiltration and skimmed milk flocculation (SMF) methods using Armored RNA as a test virus, which is also used as a control by some other studies. Prefiltration with 0.45 µm and 0.2 µm membrane disc filters were applied to eliminate solid particles before ultrafiltration to prevent the clogging of ultrafiltration devices. Test samples, processed with the sequential ultrafiltration method, were centrifuged at two different speeds. An increased speed resulted in lower recovery and positivity rates of Armored RNA. On the other hand, SMF resulted in relatively consistent recovery and positivity rates of Armored RNA. Additional tests conducted with environmental water samples demonstrated the utility of SMF to concentrate other microbial fractions. The partitioning of viruses into solid particles might have an impact on the overall recovery rates, considering the prefiltration step applied before the ultrafiltration of wastewater samples. SMF with prefiltration performed better when applied to environmental water samples due to lower solid concentrations in the samples and thus lower partitioning rates to solids. In the present study, the idea of using a sequential ultrafiltration method arose from the necessity to decrease the final volume of the viral concentrates during the COVID-19 pandemic, when the supply of the commonly used ultrafiltration devices was limited, and there was a need for the development of alternative viral concentration methods.

Quantification of human enteric viruses as alternative indicators of fecal pollution to evaluate wastewater treatment processes.

We investigated the potential use and quantification of human enteric viruses in municipal wastewater samples of Winnipeg (Manitoba, Canada) as alternative indicators of contamination and evaluated the processing stages of the wastewater treatment plant. During the fall 2019 and winter 2020 seasons, samples of raw sewage, activated sludge, effluents, and biosolids (sludge cake) were collected from the North End Sewage Treatment Plant (NESTP), which is the largest wastewater treatment plant in the City of Winnipeg. DNA (Adenovirus and crAssphage) and RNA enteric viruses (Pepper mild mottle virus, Norovirus genogroups GI and GII, Rotavirus Astrovirus, and Sapovirus) as well as the uidA gene found in Escherichia coli were targeted in the samples collected from the NESTP. Total nucleic acids from each wastewater treatment sample were extracted using a commercial spin-column kit. Enteric viruses were quantified in the extracted samples via quantitative PCR using TaqMan assays. Overall, the average gene copies assessed in the raw sewage were not significantly different (p-values ranged between 0.1023 and 0.9921) than the average gene copies assessed in the effluents for DNA and RNA viruses and uidA in terms of both volume and biomass. A significant reduction (p-value ≤ 0.0438) of Adenovirus and Noroviruses genogroups GI and GII was observed in activated sludge samples compared with those for raw sewage per volume. Higher GCNs of enteric viruses were observed in dewatered sludge samples compared to liquid samples in terms of volume (g of sample) and biomass (ng of nucleic acids). Enteric viruses found in gene copy numbers were at least one order of magnitude higher than the E. coli marker uidA, indicating that enteric viruses may survive the wastewater treatment process and viral-like particles are being released into the aquatic environment. Viruses such as Noroviruses genogroups GI and GII, and Rotavirus were detected during colder months. Our results suggest that Adenovirus, crAssphage, and Pepper mild mottle virus can be used confidently as complementary viral indicators of human fecal pollution.

Metagenomic community composition and resistome analysis in a full-scale cold climate wastewater treatment plant.

BACKGROUND: Wastewater treatment plants are an essential part of maintaining the health and safety of the general public. However, they are also an anthropogenic source of antibiotic resistance genes. In this study, we characterized the resistome, the distribution of classes 1-3 integron-integrase genes (intI1, intI2, and intI3) as mobile genetic element biomarkers, and the bacterial and phage community compositions in the North End Sewage Treatment Plant in Winnipeg, Manitoba. Samples were collected from raw sewage, returned activated sludge, final effluent, and dewatered sludge. A total of 28 bacterial and viral metagenomes were sequenced over two seasons, fall and winter. Integron-integrase genes, the 16S rRNA gene, and the coliform beta-glucuronidase gene were also quantified during this time period. RESULTS: Bacterial classes observed above 1% relative abundance in all treatments were Actinobacteria (39.24% ± 0.25%), Beta-proteobacteria (23.99% ± 0.16%), Gamma-proteobacteria (11.06% ± 0.09%), and Alpha-proteobacteria (9.18 ± 0.04%). Families within the Caudovirales order: Siphoviridae (48.69% ± 0.10%), Podoviridae (23.99% ± 0.07%), and Myoviridae (19.94% ± 0.09%) were the dominant phage observed throughout the NESTP. The most abundant bacterial genera (in terms of average percent relative abundance) in influent, returned activated sludge, final effluent, and sludge, respectively, includes Mycobacterium (37.4%, 18.3%, 46.1%, and 7.7%), Acidovorax (8.9%, 10.8%, 5.4%, and 1.3%), and Polaromonas (2.5%, 3.3%, 1.4%, and 0.4%). The most abundant class of antibiotic resistance in bacterial samples was tetracycline resistance (17.86% ± 0.03%) followed by peptide antibiotics (14.24% ± 0.03%), and macrolides (10.63% ± 0.02%). Similarly, the phage samples contained a higher prevalence of macrolide (30.12% ± 0.30%), peptide antibiotic (10.78% ± 0.13%), and tetracycline (8.69% ± 0.11%) resistance. In addition, intI1 was the most abundant integron-integrase gene throughout treatment (1.14 × 104 gene copies/mL) followed by intI3 (4.97 × 103 gene copies/mL) while intI2 abundance remained low (6.4 × 101 gene copies/mL). CONCLUSIONS: Wastewater treatment successfully reduced the abundance of bacteria, DNA phage and antibiotic resistance genes although many antibiotic resistance genes remained in effluent and biosolids. The presence of integron-integrase genes throughout treatment and in effluent suggests that antibiotic resistance genes could be actively disseminating resistance between both environmental and pathogenic bacteria.

Extraction and Detection of Avian Influenza Virus From Wetland Sediment Using Enrichment-Based Targeted Resequencing.

Early virus detection and characterization is key to successful avian influenza virus (AIV) surveillance for the health of humans as well as domestic poultry. We explored a novel sampling approach and molecular strategy using sediment from wetlands and outdoor waterbodies on poultry farms as a population-level proxy of AIV activity in waterfowls. RNA was extracted using the MoBio RNA PowerSoil Total RNA isolation kit with additional chloroform extraction steps to reduce PCR inhibition. AIV matrix protein (MP) gene was detected in 42/345 (12.2%) samples by RT-qPCR; an additional 64 (18.6%) samples showed evidence of amplification below the threshold and were categorized as "suspect positive." Enrichment-based targeted resequencing (TR) identified AIV sequences in 79/345 (22.9%) samples. TR probes were designed for MP, hemagglutinin (HA), and neuraminidase (NA), however PB2 and PA were also identified. Although RT-qPCR and TR only had fair-moderate agreement, RT-qPCR positivity was predictive of TR-positivity both when using only strictly positive RT-qPCR samples (OR = 11.29) and when coding suspect positives as positive (OR = 7.56). This indicates that RT-qPCR could be used as a screening tool to select samples for virus characterization by TR and that future studies should consider RT-qPCR suspect positives to be positive samples for subsequent resequencing when avoiding false negatives is the priority, for instance in a diagnostic test, and to consider suspect positives to be negative samples when cost efficiency over a large number of samples is the priority, for instance in a surveillance program. A total of 13 HA (H1-7, H9-13, H16) and 9 NA (N1-9) subtypes were identified, with a maximum of 8 HA and 8 NA subtypes detected in a single sample. The optimized RNA extraction and targeted resequencing methods provided increased virus detection and subtyping characterization that could be implemented in an AIV surveillance system.

Short-term deleterious effects of standard isolation and cultivation methods on new tropical freshwater microalgae strains.

Algae with potential biotechnological applications in different industries are commonly isolated from the environment in order to obtain pure (axenic) stocks that can be safely stored for long periods of time. To obtain axenic cultures, antibiotics are frequently employed, and cryopreservation is applied to preserve standing stocks. However, many of these now standard methods were developed using strains derived from pristine to near-pristine environments and cold to temperate regions. The potential effect of the said methods on the life cycle and biochemical profile of algae isolates from hyper-eutrophic and constant high-temperature tropical regions is not well understood. These effects could potentially render them unsuitable for their intended biotechnological application. In this study, we conducted a genetic characterization (18S rRNA) and evaluated the effect of purification (the use of the antibiotic chloramphenicol, CAP) and cryopreservation (dimethyl sulfoxide; DMSO-sucrose mix and glycerol) on the growth rate and lipid content of three new tropical freshwater algal isolates: Chorella sp. M2, Chlorella sp. M6, and Scenedesmus sp. R3, obtained from the Ecuadorian coast. The genetic and morphological characterization revealed a clear discrimination between these strains. All strains cultured with CAP exhibited a lower growth rate. Subsequent to cryopreservation, Chorella sp. M2, Chlorella sp. M6, and Scenedesmus sp. R3 presented no significant difference in growth rate between the cryopreservants. Further, a significantly higher lipid content was observed in the biomass cryopreserved with glycerol in relation to the DMSO-sucrose, with Chorella sp. M2 and Chlorella sp. M6 having twice as much as they had in the first treatment. These results highlight the relevance of selecting an appropriate method for storage, as the materials used can affect the biological performance of different tropical species, although it is still to be determined if the effects observed in this study are long lasting in subsequent cultures of these algae.

Human Activity Determines the Presence of Integron-Associated and Antibiotic Resistance Genes in Southwestern British Columbia.

The dissemination of antibiotic resistant bacteria from anthropogenic sources into the environment poses an emerging public health threat. Antibiotic resistance genes (ARGs) and gene-capturing systems such as integron-associated integrase genes (intI) play a key role in alterations of microbial communities and the spread of antibiotic resistant bacteria into the environment. In order to assess the effect of anthropogenic activities on watersheds in southwestern British Columbia, the presence of putative antibiotic resistance and integrase genes was analyzed in the microbiome of agricultural, urban influenced, and protected watersheds. A metagenomics approach and high-throughput quantitative PCR (HT qPCR) were used to screen for elements of resistance including ARGs and intI. Metagenomic sequencing of bacterial genomic DNA was used to characterize the resistome of microbial communities present in watersheds over a 1-year period. There was a low prevalence of ARGs relative to the microbial population (<1%). Analysis of the metagenomic sequences detected a total of 60 elements of resistance including 46 ARGs, intI1, and groEL/intI1 genes and 12 quaternary ammonium compounds (qac) resistance genes across all watershed locations. The relative abundance and richness of ARGs was found to be highest in agriculture impacted watersheds compared to urban and protected watersheds. A downstream transport pattern was observed in the impacted watersheds (urban and agricultural) during dry months. Similar to other reports, this study found a strong association between intI1 and ARGs (e.g., sul1), an association which may be used as a proxy for anthropogenic activities. Chemical analysis of water samples for three major groups of antibiotics was below the detection limit. However, the high richness and gene copy numbers (GCNs) of ARGs in impacted sites suggest that the effects of effluents on microbial communities are occurring even at low concentrations of antimicrobials in the water column. Antibiotic resistance and integrase genes in a year-long metagenomic study showed that ARGs were driven mainly by environmental factors from anthropogenized sites in agriculture and urban watersheds. Environmental factors such as land-use and water quality parameters accounted for 45% of the variability observed in watershed locations.

Beaver Fever: Whole-Genome Characterization of Waterborne Outbreak and Sporadic Isolates To Study the Zoonotic Transmission of Giardiasis.

Giardia causes the diarrheal disease known as giardiasis; transmission through contaminated surface water is common. The protozoan parasite's genetic diversity has major implications for human health and epidemiology. To determine the extent of transmission from wildlife through surface water, we performed whole-genome sequencing (WGS) to characterize 89 Giardia duodenalis isolates from both outbreak and sporadic infections: 29 isolates from raw surface water, 38 from humans, and 22 from veterinary sources. Using single nucleotide variants (SNVs), combined with epidemiological data, relationships contributing to zoonotic transmission were described. Two assemblages, A and B, were identified in surface water, human, and veterinary isolates. Mixes of zoonotic assemblages A and B were seen in all the community waterborne outbreaks in British Columbia (BC), Canada, studied. Assemblage A was further subdivided into assemblages A1 and A2 based on the genetic variation observed. The A1 assemblage was highly clonal; isolates of surface water, human, and veterinary origins from Canada, United States, and New Zealand clustered together with minor variation, consistent with this being a panglobal zoonotic lineage. In contrast, assemblage B isolates were variable and consisted of several clonal lineages relating to waterborne outbreaks and geographic locations. Most human infection isolates in waterborne outbreaks clustered with isolates from surface water and beavers implicated to be outbreak sources by public health. In-depth outbreak analysis demonstrated that beavers can act as amplification hosts for human infections and can act as sources of surface water contamination. It is also known that other wild and domesticated animals, as well as humans, can be sources of waterborne giardiasis. This study demonstrates the utility of WGS in furthering our understanding of Giardia transmission dynamics at the water-human-animal interface.IMPORTANCEGiardia duodenalis causes large numbers of gastrointestinal illness in humans. Its transmission through the contaminated surface water/wildlife intersect is significant, and the water-dwelling rodents beavers have been implicated as one important reservoir. To trace human infections to their source, we used genome techniques to characterize genetic relationships among 89 Giardia isolates from surface water, humans, and animals. Our study showed the presence of two previously described genetic assemblages, A and B, with mixed infections detected from isolates collected during outbreaks. Study findings also showed that while assemblage A could be divided into A1 and A2, A1 showed little genetic variation among animal and human hosts in isolates collected from across the globe. Assemblage B, the most common type found in the study surface water samples, was shown to be highly variable. Our study demonstrates that the beaver is a possible source of human infections from contaminated surface water, while acknowledging that theirs is only one role in the complex cycle of zoonotic spread. Mixes of parasite groups have been detected in waterborne outbreaks. More information on Giardia diversity and its evolution using genomics will further the understanding of the epidemiology of spread of this disease-causing protozoan.

Metagenomic Investigation of Plasma in Individuals with ME/CFS Highlights the Importance of Technical Controls to Elucidate Contamination and Batch Effects.

Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS) is a debilitating disease causing indefinite fatigue. ME/CFS has long been hypothesised to have an infectious cause; however, no specific infectious agent has been identified. We used metagenomics to analyse the RNA from plasma samples from 25 individuals with ME/CFS and compare their microbial content to technical controls as well as three control groups: individuals with alternatively diagnosed chronic Lyme syndrome (N = 13), systemic lupus erythematosus (N = 11), and healthy controls (N = 25). We found that the majority of sequencing reads were removed during host subtraction, thus there was very low microbial RNA content in the plasma. The effects of sample batching and contamination during sample processing proved to outweigh the effects of study group on microbial RNA content, as the few differences in bacterial or viral RNA abundance we did observe between study groups were most likely caused by contamination and batch effects. Our results highlight the importance of including negative controls in all metagenomic analyses, since there was considerable overlap between bacterial content identified in study samples and control samples. For example, Proteobacteria, Firmicutes, Actinobacteria, and Bacteriodes were found in both study samples and plasma-free negative controls. Many of the taxonomic groups we saw in our plasma-free negative control samples have previously been associated with diseases, including ME/CFS, demonstrating how incorrect conclusions may arise if controls are not used and batch effects not accounted for.

A comprehensive method for amplicon-based and metagenomic characterization of viruses, bacteria, and eukaryotes in freshwater samples.

BACKGROUND: Studies of environmental microbiota typically target only specific groups of microorganisms, with most focusing on bacteria through taxonomic classification of 16S rRNA gene sequences. For a more holistic understanding of a microbiome, a strategy to characterize the viral, bacterial, and eukaryotic components is necessary. RESULTS: We developed a method for metagenomic and amplicon-based analysis of freshwater samples involving the concentration and size-based separation of eukaryotic, bacterial, and viral fractions. Next-generation sequencing and culture-independent approaches were used to describe and quantify microbial communities in watersheds with different land use in British Columbia. Deep amplicon sequencing was used to investigate the distribution of certain viruses (g23 and RdRp), bacteria (16S rRNA and cpn60), and eukaryotes (18S rRNA and ITS). Metagenomic sequencing was used to further characterize the gene content of the bacterial and viral fractions at both taxonomic and functional levels. CONCLUSION: This study provides a systematic approach to separate and characterize eukaryotic-, bacterial-, and viral-sized particles. Methodologies described in this research have been applied in temporal and spatial studies to study the impact of land use on watershed microbiomes in British Columbia.

Giardia spp. Are Commonly Found in Mixed Assemblages in Surface Water, as Revealed by Molecular and Whole-Genome Characterization.

Giardia is the most common parasitic cause of gastrointestinal infections worldwide, with transmission through surface water playing an important role in various parts of the world. Giardia duodenalis (synonyms: G. intestinalis and G. lamblia), a multispecies complex, has two zoonotic subtypes, assemblages A and B. When British Columbia (BC), a western Canadian province, experienced several waterborne giardiasis outbreaks due to unfiltered surface drinking water in the late 1980s, collection of isolates from surface water, as well as from humans and beavers (Castor canadensis), throughout the province was carried out. To better understand Giardia in surface water, 71 isolates, including 29 from raw surface water samples, 29 from human giardiasis cases, and 13 from beavers in watersheds from this historical library were characterized by PCR. Study isolates also included isolates from waterborne giardiasis outbreaks. Both assemblages A and B were identified in surface water, human, and beavers samples, including a mixture of both assemblages A and B in waterborne outbreaks. PCR results were confirmed by whole-genome sequencing (WGS) for one waterborne outbreak and supported the clustering of human, water, and beaver isolates within both assemblages. We concluded that contamination of surface water by Giardia is complex, that the majority of our surface water isolates were assemblage B, and that both assemblages A and B may cause waterborne outbreaks. The higher-resolution data provided by WGS warrants further study to better understand the spread of Giardia.

Automated Gel Size Selection to Improve the Quality of Next-generation Sequencing Libraries Prepared from Environmental Water Samples.

Next-generation sequencing of environmental samples can be challenging because of the variable DNA quantity and quality in these samples. High quality DNA libraries are needed for optimal results from next-generation sequencing. Environmental samples such as water may have low quality and quantities of DNA as well as contaminants that co-precipitate with DNA. The mechanical and enzymatic processes involved in extraction and library preparation may further damage the DNA. Gel size selection enables purification and recovery of DNA fragments of a defined size for sequencing applications. Nevertheless, this task is one of the most time-consuming steps in the DNA library preparation workflow. The protocol described here enables complete automation of agarose gel loading, electrophoretic analysis, and recovery of targeted DNA fragments. In this study, we describe a high-throughput approach to prepare high quality DNA libraries from freshwater samples that can be applied also to other environmental samples. We used an indirect approach to concentrate bacterial cells from environmental freshwater samples; DNA was extracted using a commercially available DNA extraction kit, and DNA libraries were prepared using a commercial transposon-based protocol. DNA fragments of 500 to 800 bp were gel size selected using Ranger Technology, an automated electrophoresis workstation. Sequencing of the size-selected DNA libraries demonstrated significant improvements to read length and quality of the sequencing reads.

Year-Long Metagenomic Study of River Microbiomes Across Land Use and Water Quality.

Select bacteria, such as Escherichia coli or coliforms, have been widely used as sentinels of low water quality; however, there are concerns regarding their predictive accuracy for the protection of human and environmental health. To develop improved monitoring systems, a greater understanding of bacterial community structure, function, and variability across time is required in the context of different pollution types, such as agricultural and urban contamination. Here, we present a year-long survey of free-living bacterial DNA collected from seven sites along rivers in three watersheds with varying land use in Southwestern Canada. This is the first study to examine the bacterial metagenome in flowing freshwater (lotic) environments over such a time span, providing an opportunity to describe bacterial community variability as a function of land use and environmental conditions. Characteristics of the metagenomic data, such as sequence composition and average genome size (AGS), vary with sampling site, environmental conditions, and water chemistry. For example, AGS was correlated with hours of daylight in the agricultural watershed and, across the agriculturally and urban-affected sites, k-mer composition clustering corresponded to nutrient concentrations. In addition to indicating a community shift, this change in AGS has implications in terms of the normalization strategies required, and considerations surrounding such strategies in general are discussed. When comparing abundances of gene functional groups between high- and low-quality water samples collected from an agricultural area, the latter had a higher abundance of nutrient metabolism and bacteriophage groups, possibly reflecting an increase in agricultural runoff. This work presents a valuable dataset representing a year of monthly sampling across watersheds and an analysis targeted at establishing a foundational understanding of how bacterial lotic communities vary across time and land use. The results provide important context for future studies, including further analyses of watershed ecosystem health, and the identification and development of biomarkers for improved water quality monitoring systems.