Exploring taxonomic and functional microbiome of Hawaiian stream and spring irrigation water systems using Illumina and Oxford Nanopore sequencing platforms

Snapshot of cyanobacterial toxins in Pakistani freshwater bodies

Cyanobacteria are known to produce diverse secondary metabolites that are toxic to aquatic ecosystems and human health. However, data about the cyanotoxins occurrence and cyanobacterial diversity in Pakistan's drinking water reservoirs is scarce. In this study, we first investigated the presence of microcystin, saxitoxin, and anatoxin in 12 water bodies using an enzyme-linked immunosorbent assay (ELISA). The observed cyanotoxin values for the risk quotient (RQ) determined by ELISA indicated a potential risk for aquatic life and human health. Based on this result, we made a more in-depth investigation with a subset of water bodies (served as major public water sources) to analyze the cyanotoxins dynamics and identify potential producers. We therefore quantified the distribution of 17 cyanotoxins, including 12 microcystin congeners using a high-performance liquid chromatography-high-resolution tandem mass spectrometry/mass spectrometry (HPLC-HRMS/MS). Our results revealed for the first time the co-occurrence of multiple cyanotoxins and the presence of cylindrospermopsin in an artificial reservoir (Rawal Lake) and a semi-saline lake (Kallar Kahar). We also quantified several microcystin congeners in a river (Panjnad) with MC-LR and MC-RR being the most prevalent and abundant. To identify potential cyanotoxin producers, the composition of the cyanobacterial community was characterized by shotgun metagenomics sequencing. Despite the noticeable presence of cyanotoxins, Cyanobacteria were not abundant. Synechococcus was the most abundant cyanobacterial genus found followed by a small amount of Anabaena, Cyanobium, Microcystis, and Dolichospermum. Moreover, when we looked at the cyanotoxins genes coverage, we never found a complete microcystin mcy operon. To our knowledge, this is the first snapshot sampling of water bodies in Pakistan. Our results would not only help to understand the geographical spread of cyanotoxin in Pakistan but would also help to improve cyanotoxin risk assessment strategies by screening a variety of cyanobacterial toxins and confirming that cyanotoxin quantification is not necessarily related to producer abundance.

Microbiologic surveys for Baijiu fermentation are affected by experimental design

The traditional Chinese alcoholic beverage Baijiu is produced by spontaneous fermentation of grains under anaerobic conditions. While numerous studies have used metagenomic technology to investigate the microbiome of Baijiu brewing, the microbial succession mechanism of Baijiu brewing has not been fully clarified, and metagenomic strategies for microecology surveys have not been comprehensively evaluated. Using the fermentation process of strong-flavor Baijiu as a model, we compared the data for bacterial communities based on short read 16S rRNA variable regions, V3-V4, and full-length 16S regions, V1-V9, to whole metagenomic shotgun sequencing (WMS) to measure the effect of technology selection on phylogenetic and functional profiles. The results showed differences in bacterial compositions and their relation to volatiles and physicochemical variables between sequencing methods. Furthermore, the percentage of V3-V4 sequences assigned to species level was higher than the percentage of V1-V9 sequences according to SILVA taxonomy, but lower according to NCBI taxonomy (P < 0.05). In both SILVA and NCBI taxonomies, the relative abundances of bacterial communities at both the genus and family levels were more positively correlated with WMS data in the V3-V4 dataset than in the V1-V9 dataset. The WMS identified changes in abundances of multiple metabolic pathways during fermentation (P < 0.05), including "starch and sucrose metabolism," "galactose metabolism," and "fatty acid biosynthesis." Although functional predictions derived from 16S data show similar patterns to WMS, most metabolic pathway changes are uncorrelated (P > 0.05). This study provided new technical and biological insights into Baijiu production that may assist in selection of methodologies for studies of fermentation systems.

Long-read MinION™ sequencing of 16S and 16S-ITS-23S rRNA genes provides species-level resolution of Lactobacillaceae in mixed communities

The Lactobacillaceae are lactic acid bacteria harnessed to deliver important outcomes across numerous industries, and their unambiguous, species-level identification from mixed community environments is an important endeavor. Amplicon-based metataxonomics using short-read sequencing of partial 16S rRNA gene regions is widely used to support this, however, the high genetic similarity among Lactobacillaceae species restricts our ability to confidently describe these communities even at genus level. Long-read sequencing (LRS) of the whole 16S rRNA gene or the near complete rRNA operon (16S-ITS-23S) has the potential to improve this. We explored species ambiguity amongst Lactobacillaceae using in-silico tool RibDif2, which identified allele overlap when various partial and complete 16S rRNA gene and 16S-ITS-23S rRNA regions were amplified. We subsequently implemented LRS by MinION™ to compare the capacity of V3-V4, 16S and 16S-ITS-23S rRNA amplicons to accurately describe the diversity of a 20-species Lactobacillaceae mock community in practice. In-silico analysis identified more instances of allele/species overlap with V3-V4 amplicons (n = 43) compared to the 16S rRNA gene (n = 11) and partial (n = up to 15) or complete (n = 0) 16S-ITS-23S rRNA amplicons. With subsequent LRS of a DNA mock community, 80% of target species were identified using V3-V4 amplicons whilst the 16S rRNA gene and 16S-ITS-23S rRNA region amplicons resulted in 95 and 100% of target species being identified. A considerable reduction in false-positive identifications was also seen with 16S rRNA gene (n = 3) and 16S-ITS-23S rRNA region (n = 9) amplicons compared with V3-V4 amplicons (n = 43). Whilst the target species affected by allele overlap in V3-V4 and 16S rRNA gene sequenced mock communities were predicted by RibDif2, unpredicted species ambiguity was observed in 16S-ITS-23S rRNA sequenced communities. Considering the average nucleotide identity (ANI) between ambiguous species (~97%) and the basecall accuracy of our MinION™ sequencing protocol (96.4%), the misassignment of reads between closely related taxa is to be expected. With basecall accuracy exceeding 99% for recent MinION™ releases, the increased species-level differentiating power promised by longer amplicons like the 16S-ITS-23S rRNA region, may soon be fully realized.