Dataset Information

Bacterial Community 16S rRNA Gene Sequencing Characterizes Riverine Microbial Impact on Lake Michigan.

ABSTRACT: Restoration of degraded aquatic habitats is critical to preserve and maintain ecosystem processes and economic viability. Effective restoration requires contaminant sources identification. Microbial communities are increasingly used to characterize fecal contamination sources. The objective was to determine whether nearshore and adjacent beach bacterial contamination originated from the Grand Calumet River, a highly urbanized aquatic ecosystem, and to determine if there were correlations between pathogens/feces associated bacteria in any of the samples to counts of the pathogen indicator species Escherichia coli. Water samples were collected from the river, river mouth, nearshore, and offshore sites along southern Lake Michigan. Comparisons among communities were made using beta diversity distances (weighted and unweighted Unifrac, and Bray Curtis) and Principal Coordinate Analysis of 16S rRNA gene Illumina sequence data that indicated river bacterial communities differed significantly from the river mouth, nearshore lake, and offshore lake samples. These differences were further supported using Source Tracker software that indicated nearshore lake communities differed significantly from river and offshore samples. Among locations, there was separation by sampling date that was associated with environmental factors (e.g., water and air temperature, water turbidity). Although about half the genera (48.1%) were common to all sampling sites, linear discriminant analysis effect size indicated there were several taxa that differed significantly among sites; there were significant positive correlations of feces-associated genera with E. coli most probable numbers. Results collectively highlight that understanding microbial communities, rather than relying solely on select fecal indicators with uncertain origin, are more useful for developing strategies to restore degraded aquatic habitats.

SUBMITTER: Nakatsu CH

PROVIDER: S-EPMC6527805 | biostudies-literature | 2019

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Bacterial Community 16S rRNA Gene Sequencing Characterizes Riverine Microbial Impact on Lake Michigan.

Nakatsu Cindy H CH Byappanahalli Muruleedhara N MN Nevers Meredith B MB

Frontiers in microbiology 20190514

Restoration of degraded aquatic habitats is critical to preserve and maintain ecosystem processes and economic viability. Effective restoration requires contaminant sources identification. Microbial communities are increasingly used to characterize fecal contamination sources. The objective was to determine whether nearshore and adjacent beach bacterial contamination originated from the Grand Calumet River, a highly urbanized aquatic ecosystem, and to determine if there were correlations between ...[more]

PMID: 31139161

Similar Datasets

Project description:Introduction:Microbiome research based on high-throughput sequencing has grown exponentially in recent years, but methodological variations can easily undermine the reproducibility across studies. Objectives:To systematically evaluate the comparability of sequencing results of 16S rRNA gene sequencing (16Ss)- and shotgun metagenomic sequencing (SMs)-based microbial community profiling in laboratories under routine conditions. Methods:We designed a multicenter study across 35 participating laboratories in China using designed mock communities and homogenized fecal samples. Results:A wide range of practices and approaches was reported by the participating laboratories. The observed microbial compositions of the mock communities in 46.2% (12/26) of the 16Ss and 82.6% (19/23) of the SMs laboratories had significant correlations with the expected result (Spearman r>0.59, P <0.05). The results from laboratories with near-identical protocols showed slight interlaboratory deviations. However, a high degree of interlaboratory deviation was found in the observed abundances of specific taxa, such as Bacteroides spp. (range: 0.3%-53.5%), Enterococci spp. (range: 0.8%-43.9%) and Fusobacterium spp. (range: 0.1%-39.8%). SMs performed better than 16Ss in detecting low-abundance bacteria (B. bifidum). The differences in DNA extraction methods, amplified regions and bioinformatics analysis tools (taxonomic classifiers and database) were important factors causing interlaboratory deviations. Addressing laboratory contamination is an urgent task because various sources of unexpected microbes were found in negative control samples. Conclusions:Well-defined control samples, such as the mock communities in this study, should be routinely used in microbiome research for monitoring potential biases. The findings in this study will provide guidance in the choice of more reasonable operating procedures to minimize potential methodological biases in revealing human microbiota composition.

Project description:Bacterial endosymbionts that produce important phenotypic effects on their hosts are common among plant sap-sucking insects. Aphids have become a model system of insect-symbiont interactions. However, endosymbiont research has focused on a few aphid species, making it necessary to make greater efforts to other aphid species through different regions, in order to have a better understanding of the role of endosymbionts in aphids as a group. Aphid endosymbionts have frequently been studied by PCR-based techniques, using species-specific primers, nevertheless this approach may omit other non-target bacteria cohabiting a particular host species. Advances in high-throughput sequencing technologies are complementing our knowledge of microbial communities by allowing us the study of whole microbiome of different organisms. We used a 16S rRNA amplicon sequencing approach to study the microbiome of aphids in order to describe the bacterial community diversity in introduced populations of the cereal aphids, Sitobion avenae and Rhopalosiphum padi in Chile (South America). An absence of secondary endosymbionts and two common secondary endosymbionts of aphids were found in the aphids R. padi and S. avenae, respectively. Of those endosymbionts, Regiella insecticola was the dominant secondary endosymbiont among the aphid samples. In addition, the presence of a previously unidentified bacterial species closely related to a phytopathogenic Pseudomonad species was detected. We discuss these results in relation to the bacterial endosymbiont diversity found in other regions of the native and introduced range of S. avenae and R. padi. A similar endosymbiont diversity has been reported for both aphid species in their native range. However, variation in the secondary endosymbiont infection could be observed among the introduced and native populations of the aphid S. avenae, indicating that aphid-endosymbiont associations can vary across the geographic range of an aphid species. In addition, we discuss the potential role of aphids as vectors and/or alternative hosts of phytopathogenic bacteria.

Project description:BackgroundAccurate identification of bacterial communities is crucial for research applications, diagnostics, and clinical interventions. Although 16S ribosomal RNA (rRNA) gene sequencing is a widely employed technique for bacterial taxonomic classification, it often results in misclassified or unclassified bacterial taxa. This study sought to refine the full-length 16S rRNA gene sequencing protocol using the MinION sequencer, focusing on the V1-V9 regions. Our methodological enquiry examined several factors, including the number of PCR amplification cycles, choice of primers and Taq polymerase, and specific sequence databases and workflows employed. We used a microbial standard comprising eight bacterial strains (five gram-positive and three gram-negative) in known proportions as a validation control.ResultsBased on the MinION protocol, we employed the microbial standard as the DNA template for the 16S rRNA gene amplicon sequencing procedure. Our analysis showed that an elevated number of PCR amplification cycles introduced PCR bias, and the selection of Taq polymerase and primer sets significantly affected the subsequent analysis. Bacterial identification at genus level demonstrated Pearson correlation coefficients ranging from 0.73 to 0.79 when assessed using BugSeq, Kraken-Silva and EPI2ME-16S workflows. Notably, the EPI2ME-16S workflow exhibited the highest Pearson correlation with the microbial standard, minimised misclassification, and increased alignment accuracy. At the species taxonomic level, the BugSeq workflow was superior, with a Pearson correlation coefficient of 0.92.ConclusionsThese findings emphasise the importance of careful selection of PCR settings and a well-structured analytical framework for 16S rRNA full-length gene sequencing. The results showed a robust correlation between the predicted and observed bacterial abundances at both the genus and species taxonomic levels, making these findings applicable across diverse research contexts and with clinical utility for reliable pathogen identification.

Project description:Two consortia (Consortium A and Consortium B) that can use 1,4-dioxane (a groundwater contaminant of emerging concern) as the sole carbon source were enriched from Rice University (Houston, TX, USA) campus soil. Phylogenetic analysis by 16S rRNA sequencing revealed the dominant genus in both of the consortia is Mycobacterium (56% in Consortium A and 49% in Consortium B). The predominance of Mycobacterium spp, in these consortia support the notion that this is an important and commonly encountered genus of dioxane degraders. Among other genera present that make at least 2% of these consortia, only Afipia encompasses a strain (i.e., Afipia sp. D1) that was reported to degrade dioxane as sole carbon and energy source. A nested PCR analysis using two degenerate primers to target the hydroxylase alpha subunit of groups 3 to 6 SDIMOs was performed to gain insights into which enzymes were responsible for dioxane degradation by these consortia. The purified products obtained from the second PCR run were sequenced and compared to genes databases (NCBI) encompassing all of the currently reported SDIMOs. The dominant SDIMO genes in Consortium A corresponded to a group-6 putative propane monooxygenase-like SDIMO (98.8%); while in Consortium B, SDIMO genes from both groups 5 (47.3%) and 6 (51.9%) were observed. In both consortia, the relative abundance of thmA/dxmA gene was negligible (0.03%), which is consistent with the negative amplification of these genes as verified in qPCR. Overall, the high relative abundance of group-6 putative propane monooxygenases in our two consortia suggests the novel finding that group 6-SDIMOs could also play an important role in dioxane degradation. This underscores the need for further research on genes and enzymes involved in dioxane biodegradation to develop novel biomarkers that can be useful for forensic analysis and performance assessment of bioremediation and natural attenuation at dioxane-impacted sites. DNA was extracted from bacteria biomass harvested in exponential growth phase, when half or more of the added dioxane (100 mg/L) was consumed. Total DNA extractions were performed using the UltraClean® Microbial DNA Isolation Kit (MO BIO, Carlsbad, CA, USA) according to the manufacturer’s protocol. The V4 region of the 16S rRNA gene was amplified by PCR using the forward 515F and reverse 806R primers. Sequencing was performed at MR DNA (www.mrdnalab.com, Shallowater, TX, USA) by Illumina MiSeq paired-end sequencing (approximately 2×300 bp as the read length). Sequence data were processed using MR DNA analysis pipeline. Operational taxonomic units (OTUs) were defined by clustering at 3% divergence (97% similarity). Final OTUs were taxonomically classified using BLASTn against the RDPII (http://rdp.cme.msu.edu) and NCBI (www.ncbi.nlm.nih.gov) databases.Previously designed degenerate primers NVC57, NVC58, NVC65 and NVC66 to target conserved regions in the soluble di-iron monooxygenases (SDIMO) alpha subunit gene (groups 3 to 6) were used to examine the presence and diversity of SDIMO genes in these two consortia. A nested PCR strategy was used to increase the PCR product yield. In the first run, the PCR mixture contained 1 µL of NVC65 and NVC58 primer mixture (10 µM), 20 ng of the extracted genomic DNA, 12.5 µL of KAPA HiFi HotStart ReadyMix (2X) (KAPA Biosystems, Wilmington, MA, USA), and nuclease-free water to yield a total volume of 25 µL. PCR was performed in a Bio-Rad Thermal Cycler (Bio-Rad, Hercules, CA, USA) with the following temperature profile: initial denaturation (94°C, 5 min), then 29 amplification cycles (94°C for 30 s, 55°C for 30 s, 72°C for 1 min per kb) and a final extension (72°C for 5 min). The length of the PCR products in the first run was checked by 1% agarose gel and DNA bands of the correct size (1100 bp) were excised and purified. 20 ng of the purified PCR product was used as the DNA template in the second run, with the second set of primers (NVC57 and NVC66). The purified product (420 bp) from the second PCR was sent to MR DNA (www.mrdnalab.com, Shallowater, TX, USA) for Illumina MiSeq paired-end sequencing (approximately 2×300 bp as the read length). Sequence data were processed using MR DNA analysis pipeline. Operational taxonomic units (OTUs) were defined by clustering at 3% divergence (97% similarity). A database including all of the currently reported SDIMO genes on NCBI was created and used to taxonomically classify the final OTUs.

Project description:One of the major methods to identify microbial community composition, to unravel microbial population dynamics, and to explore microbial diversity in environmental samples is high-throughput DNA- or RNA-based 16S rRNA (gene) amplicon sequencing in combination with bioinformatics analyses. However, focusing on environmental samples from contrasting habitats, it was not systematically evaluated (i) which analysis methods provide results that reflect reality most accurately, (ii) how the interpretations of microbial community studies are biased by different analysis methods and (iii) if the most optimal analysis workflow can be implemented in an easy-to-use pipeline. Here, we compared the performance of 16S rRNA (gene) amplicon sequencing analysis tools (i.e., Mothur, QIIME1, QIIME2, and MEGAN) using three mock datasets with known microbial community composition that differed in sequencing quality, species number and abundance distribution (i.e., even or uneven), and phylogenetic diversity (i.e., closely related or well-separated amplicon sequences). Our results showed that QIIME2 outcompeted all other investigated tools in sequence recovery (>10 times fewer false positives), taxonomic assignments (>22% better F-score) and diversity estimates (>5% better assessment), suggesting that this approach is able to reflect the in situ microbial community most accurately. Further analysis of 24 environmental datasets obtained from four contrasting terrestrial and freshwater sites revealed dramatic differences in the resulting microbial community composition for all pipelines at genus level. For instance, at the investigated river water sites Sphaerotilus was only reported when using QIIME1 (8% abundance) and Agitococcus with QIIME1 or QIIME2 (2 or 3% abundance, respectively), but both genera remained undetected when analyzed with Mothur or MEGAN. Since these abundant taxa probably have implications for important biogeochemical cycles (e.g., nitrate and sulfate reduction) at these sites, their detection and semi-quantitative enumeration is crucial for valid interpretations. A high-performance computing conformant workflow was constructed to allow FAIR (Findable, Accessible, Interoperable, and Re-usable) 16S rRNA (gene) amplicon sequence analysis starting from raw sequence files, using the most optimal methods identified in our study. Our presented workflow should be considered for future studies, thereby facilitating the analysis of high-throughput 16S rRNA (gene) sequencing data substantially, while maximizing reliability and confidence in microbial community data analysis.

Dataset Information

Bacterial Community 16S rRNA Gene Sequencing Characterizes Riverine Microbial Impact on Lake Michigan.

Publications

Bacterial Community 16S rRNA Gene Sequencing Characterizes Riverine Microbial Impact on Lake Michigan.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets