key: cord-0748595-k0nt4k5r
authors: Jani, Kunal; Bandal, Jayshree; Shouche, Yogesh; Shafi, Shuja; Azhar, Esam I.; Zumla, Alimuddin; Sharma, Avinash
title: Extended Ecological Restoration of Bacterial Communities in the Godavari River During the COVID-19 Lockdown Period: a Spatiotemporal Meta-analysis
date: 2021-07-05
journal: Microb Ecol
DOI: 10.1007/s00248-021-01781-0
sha: 658f15d1a2e463c11d2c0c91a78a55aca6033ce9
doc_id: 748595
cord_uid: k0nt4k5r

The unprecedented COVID-19 pandemic has had major impact on human health worldwide. Whilst national and international COVID-19 lockdown and travel restriction measures have had widespread negative impact on economies and mental health, they may have beneficial effect on the environment, reducing air and water pollution. Mass bathing events (MBE) also known as Kumbh Mela are known to cause perturbations of the ecosystem affecting resilient bacterial populations within water of rivers in India. Lockdowns and travel restrictions provide a unique opportunity to evaluate the impact of minimum anthropogenic activity on the river water ecosystem and changes in bacterial populations including antibiotic-resistant strains. We performed a spatiotemporal meta-analysis of bacterial communities of the Godavari River, India. Targeted metagenomics revealed a 0.87-fold increase in the bacterial diversity during the restricted activity of lockdown. A significant increase in the resilient phyla, viz. Proteobacteria (70.6%), Bacteroidetes (22.5%), Verrucomicrobia (1.8%), Actinobacteria (1.2%) and Cyanobacteria (1.1%), was observed. There was minimal incorporation of allochthonous bacterial communities of human origin. Functional profiling using imputed metagenomics showed reduction in infection and drug resistance genes by − 0.71-fold and − 0.64-fold, respectively. These observations may collectively indicate the positive implications of COVID-19 lockdown measures which restrict MBE, allowing restoration of the river ecosystem and minimise the associated public health risk. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s00248-021-01781-0.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a novel beta-CoV group 2B zoonotic infection of humans first identified in December 2019 from a cluster of human coronavirus disease 2019 (COVID- 19) cases in Wuhan, China [1, 2] . The ensuing COVID-19 global pandemic has severely impacted the healthcare systems worldwide, causing ~ 3.39 million deaths as of May 19, 2021 [2-4] . World leaders responded by imposing restriction on the movement of people, lockdowns and implementation of public health measures including social distancing, hand washing and wearing of face masks [5] . India was not an exception and imposed its nationwide lockdown on March 25, 2020. Considering India's population and the population density, an event of community transmission might have led to a havoc [6] . Most countries have implemented lockdown measures with national and international travel restrictions, which have affected all facets of life globally with major negative social and economic consequences [7] . On the other hand, the reduced movement of humans may have played a role in cleansing and restoration of environment, biodiversity and ecosystems that were under the perpetual anthropogenic stress [8] . Hence, it became important to evaluate the implications and/or beneficial impact of COVID-19 pandemic on the environment [9] .

A recent study by Somani et al. [6] , detailing the comprehensive analysis of various pollutants, suggested positive implications of restricted anthropogenic activities during the COVID-19 lockdown, particularly on the quality of water, level of air pollution, noise pollution, waste generation, etc. Restricted movement and lower traffic were also found to aid in dropping the levels of nitrogen dioxide and noise pollution [10] . Similarly, a review by Rume and Islam [9] on environmental implications of COVID-19 underscores the improved quality of air with the reduction in emission of greenhouse gases and decrease in water pollution across the globe. Additionally, the review also highlights the plausibility of ecological restoration especially at the tourist spots. In Indian, the concentrations of pollutant such as PM 2.5 , PM 10 , O 3 , NO 2 , SO 2 and CO were consistently noted to be higher than its permissible limits, but during the course of lockdown, it was found to be under the permissible limits (with some occasional spikes) in the various megacities [5, 11, 12] . Additionally, the stringent lockdown had its contribution in cleansing the rivers and other aquatic ecosystems which are otherwise under the continuous burden of discharges from untreated sewage and industries [13] [14] [15] . Surprisingly, in the matter of weeks of lockdown, the water quality of the longest river of India, i.e. Ganga, was more pronounce than the collective efforts of more than a decade [16] [17] [18] . Singhal and Matto [17] highlighted that the sudden decrease in the number of pilgrims and the 500% decrease in the introduction of domestic and industrial effluents have led the improvement of the water quality. Similarly, in case of rivers Yamuna, Cauvery and Krishna, Central Pollution Control Board [19] reported concentrations of various factors like dissolved oxygen (DO), pH and coliform count to be under the permissible limits during the lockdown. Alternatively, the COVID-19 was associated with the inevitable situation of generation of enormous biomedical waste. Ineffective disposal of the biomedical waste was related to the environmental hazard especially affecting the water bodies [9] . The latter was also supported by a case study from Bangladesh indicating responders suffering effective disposal of COVID-19 preventives, leading to contamination of the water bodies [20] .

Apart from the factors known to influence the water pollution, in India, the recurring pilgrimages that are hosted on the banks of the river are a source of unique perturbations to the aquatic ecosystem [21] . One such pilgrimage is the Kumbh Mela, the world's largest religious mass gathering event hosted along the banks of holy rivers in India, and it attracts millions of pilgrims from India and other countries [22] [23] [24] . The Kumbh Mela involves the ritual practice of bathing in the river which substantially increases the risk associated with the spread of water infectious diseases among the fellow pilgrims [23] and also causes perturbations of the ecosystems. Our spatiotemporal investigation of the river Godavari during the 2015 Kumbh Mela event by using bacterial communities as a proxy revealed the deleterious impact of the mass bathing event (MBE). It was observed to be associated with the loss of diversity, with invasive establishment of allochthonous bacterial communities, enormous introduction of human faecal microbiota and increased virulence and drug resistance genes [25] .

The COVID-19 lockdown measures, travel restrictions and associated scaling down of mass gathering religious events over the past year provide a unique opportunity to evaluate the impact of minimum anthropogenic activity on the river water ecosystem. Recent studies on the river Ganga found a reduction in the biological oxygen demand (BOD), faecal coliform count and nitrate (NO 3 ) concentration and an increase in dissolved oxygen during the lockdown [26, 27] . However, both the studies primarily focus on the quality of the river water and adopted a classical approach of assessment of physiochemical factors and microbial contaminants, and investigation was restricted to the faecal coliform and total coliform counts [26] . However, it fails to provide insight into the bacterial communities thriving in the environment at given time and space. As the bacterial communities are sensitive to subtle variations in their environments, it can serve as a proxy to delineate the extend of ecological disturbances or restoration, if any [28] . In the recent year, the next-generation sequencing has widely been used and has proven its competency in determining such anomalies associated with changing environment or in human health [29] .

Therefore, the current study attempts to elucidate the implication of COVID-19 lockdown on the environment by analysing change in the river bacterial communities. Using targeted metagenomics, we performed a spatiotemporal and meta-analysis of bacterial community structure and composition of the Godavari River, India, during the COVID-19 lockdown period (1) to understand the impact of COVID-19-mediated restrictions on the river ecosystem and (2) to examine the extent of ecosystem restoration, if any.

The current study is an extension of the spatiotemporal investigation of the river bacterial communities under the influence of the mass bathing event at Nashik Kumbh Mela 2015 [25] , and it involved sample collection from same sites to reduce experimental heterogeneity. Five different sampling sites were utilised during the 2015 event [25] representing upstream sites (low human interference) and mass bathing sites (higher human activity). Water samples were collected in triplicate from the mass bathing sites located in Nashik City (Anandvalli, Gharpure Ghat and Tapovan) during the lockdown period (June 2020) implemented as a measure of containment of COVID-19 spread. However, due to travel restrictions imposed due to COVID-19 pandemic, samples were not collected from the upstream sites located in 'Trimbakeshwar' near the 'Anjaneri' forest. The samples collected during the lockdown period were termed as LS3, LS4 and LS5 corresponding to sampling site Anandvalli, Gharpure Ghat and Tapovan, respectively (Fig. 1 ).

Surface water samples were collected in triplicates from the defined bathing areas located at the river bank. The samples were stored in the sterile containers and transported to the laboratory by maintaining the cold chain. Each sample aliquot was subjected to physicochemical analysis at the Department of Microbiology, KTHM College, Nashik. Environmental parameters such as turbidity (OD), total dissolved solid (TDS), total solids (TS), total suspended solids (TSS), DO, BOD, faecal coliforms and most probable number (MPN) of total coliforms were assessed as per the procedures of APHA (1998) whilst pH (Pocket digital pH meter, Milwaukee pH600) and temperature (Vintage USA E8000) were recorded onsite. The remaining water sample was filtered through 0.22-µ Millipore filters (Merck Millipore, USA) and stored at − 20 °C until further processing.

Community DNA extraction was carried out from the stored 0.22-µ Millipore filters using DNeasy Power Water Kit (Qiagen, The Netherlands). The qualitative and quantitative assessment of the DNA was performed in accordance with Sharma et al. [30] . The DNA was processed for the targeted amplicon sequencing by amplification of the V4-hyper variable region of the 16S rRNA gene using primers 515F-806R [31] . Library preparation was performed by following the protocols of Illumina, Inc. (USA). The obtained libraries were sequenced on Illumina MiSeq using 2 × 250 bp v2 chemistry. The sequences are available on the NCBI SRA portal under BioProject Id PRJNA698474.

In order to maintain the uniformity, the amplicon sequences generated during the Kumbh Mela 2015 [25] study were retrieved prior to analysing the current sequencing data (BioProject Id PRJNA383664). The bioinformatics analysis was in accordance with Jani et al. [25] which involved analysis using QIIME (Quantitative Insights Into Microbial Ecology, v1.9) [32] . Assembly of the pairedend reads using FLASH (Fast Length Adjustment of Short reads) tool, followed by adaptor and quality-filtering using Cutadapt (v3.2) and mothur (v1.32), respectively [33] [34] [35] . The operational taxonomic unit (OTU)-picking was performed using a closed reference approach against the Greengenes database with the UCLUST algorithm with sequence similarity of 97% [36, 37] . Taxonomic assignment of the OTUs was performed by adopting the RDP naïve Bayesian classifier against the Greengenes database [38] . Further, downstream analysis precludes sequence rarefaction to the lowest number of reads per sample. Estimates of the alpha diversity indices, viz. OTU richness, Shannon, Chao1, and Good's coverage, were carried out using QIIME (v1.9). Statistical analysis and data visualisation were performed in R using the packages cowplot, corrplot, Hmisc, PerformanceAnalytics, phyloseq, vegan and ggplot2 [39] [40] [41] [42] [43] [44] [45] .

The source and establishment of the human microbiota was assessed by employing an approach of Bayesian mixing model [46] . The model aids in the identification and quantification of potential microbial communities of human origin. Microbiome data (examining various body sites, i.e. skin, oral and stool) of Indian sub-population representing both healthy and disease states was collected and compiled. The human microbiome data was used as a potential source (a reference dataset) of microbial communities whereas the mass bathing data (generated in this study and BioProject ID PRJNA383664) was used as a sink (target dataset) to deduce the origin of microbiota observed during the mass bathing event [47, 48] .

Functional profiling of these bacterial communities was performed by an approach of imputed metagenomics using PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States) [49] . Specifically, the computational argument of "metagenomic contributions" was applied to decipher the functional potential of bacterial communities observed before and during the mass bathing event as well as during the minimal state of perturbation to the river, i.e. COVID-19 lockdown [50] . 1 The map depicts the location of the sampling sites along the bank of the Godavari River, Nashik, India. Samples collected before the event are named as BS3, BS4 and BS5; samples collected during the event are DS3, DS4 and DS5; and samples collected during the COVID lockdown were termed as LS3, LS4 and LS5 that correspond to locations Anandvalli, Gharpure Ghat and Tapovan, respectively 

Assessment of the environmental parameters along mass bathing sites across three points, viz. before MBE, during MBE and lockdown, indicates substantial variability, especially when compared between the groups, i.e. before vs during MBE and during MBE vs lockdown ( Table 1 ). The estimates of the TS, TDS and TSS were found to increase during the course of MBE (by 2.08-, 1.72-and 2.45-fold, respectively) and decrease during the minimal anthropogenic activity, i.e. the period of lockdown (− 0.74-, − 0.65and − 0.81-fold, respectively). The significant increase in the TS (p ≤ 0.01), TDS (p = 0.055) and TSS (p ≤ 0.005) leads to a surge in the turbidity of water. A significant variation in turbidity was noted under the temporal gradient (p ≤ 0.02). The turbidity increased by 2.26-fold during the course of mass bathing event. In contrast, the turbidity was found to degrease by − 0.79-fold during the lockdown period. Increase in the turbidity and solute concentration during MBE corroborated the fact that the event attracts participation by millions of pilgrims at the event. Additionally, the untreated human solid and sewer waste along with the industrial and household waste being routinely discharged into the river lead to deterioration of water quality [22, 23] . Contamination of freshwater ecosystem consequently has a deleterious impact on the aquatic biota, including microorganisms [15, 51] . In contrast, the imposed restriction due the COVID-19 pandemic has been found to improvise the water quality and ecosystem, which was also noted in the case of river Damodar [27] . The observation was further supported by the study on river Ganga, suggesting positive implications of the lockdown towards ecosystem restoration [26] . Further, a study by Cooke et al. [7] highlighted both the benefits and risk due to pandemic on the freshwater fishes. In contrast, pH and temperature of the water were not found to differ significantly (p = 0.65 and p = 0.62, respectively). Further, evaluation of the oxygen relationship depicted a significant variation in the BOD (p ≤ 0.009) under temporal variation across mass bathing sites. The BOD was increased by 1.13-fold during the mass bathing event but was found to decrease by − 0.94-fold during the lockdown. The variations in BOD was the consequence of substantial variation in the concentration of the DO (p = 0.11). The concentration of the DO was decreased by − 0.6-fold but raised by 2.5-fold during the lockdown. Similarly, the microbiological assessment of the river water indicated a 67.7-fold increase in the MPN count (p = 0.34) which was found to reduce to − 0.99-fold during the lockdown. Similarly, the faecal coliform count (p = 0.052) raised by 0.98-fold and lowered to − 0.87 during the lockdown. The increase in the BOD is suggestive of Table 1 Environmental parameters and alpha diversity across the temporal variation of mass bathing event and COVID lockdown Sample ID Event the enormous influx of solutes due to bathing in the river, sediment resuspension, discharge from industries and leaching from unattended solid wastes during the MBE [51] [52] [53] . These changes in the river physiochemistry collectively have deleterious impact on the river microbial communities. Perturbation to the river bacterial communities and augmentation of specific bacterial groups like faecal coliforms pose a risk of spread of gastrointestinal diseases during the event [54, 55] . However, the lockdown was also found to be beneficial concerning the oxygen relationship of the river water which also coincides with observations of Dutta et al. [26] depicting the increased level of DO and reduction in BOD and faecal coliform counts in the river Ganga.

Analysis of high-quality paired-end reads using QIIME yielded 7469 OTUs. The reads were normalised to 105,246 prior to performing further downstream analysis. The estimates of alpha diversity index revealed substantial variability under the spatiotemporal realm; however, it was not be supported by the statistical significance ( Table 1 ). The estimates of the non-parametric Shannon index depicted variance with respect to the time; i.e. it was found to decrease by − 0.53-fold during the mass bathing event whilst it was found to recover during the lockdown by 0.87-fold. This observation suggested potential re-establishment of diversity as minimal divergence of − 0.13-fold was noted between before MBE and lockdown samples. It was further supported by the assessment of Chao1 index and observed OTUs (Table 1 ). In summary, the species richness was found to follow an order of during MBE < lockdown < before MBE samples (Fig. S1 ). Loss of bacterial diversity found during the MBE might be supported by deteriorated quality of the river water and enormous influx of organic and inorganic waste having adverse impact on the microbial communities and river ecosystem [52, 53, 56] . Alternatively, a 0.87-fold increase in the Shannon index implied the events of ecosystem restoration under the improvised water quality during the imposed lockdown period [9, 26] .

The most abundant bacterial phyla found under the spatiotemporal constraint consist of Firmicutes, Proteobacteria, Bacteroidetes, Actinobacteria, Verrucomicrobia, Planctomycetes, Cyanobacteria, [Thermi], Chlorobi and Chloroflexi contributing for 99.9% of bacterial community (Fig. 2) . The river bacterial community composition before the MBE was comprised of varied bacterial taxa having relative abundance ≥ 1%, viz. Firmicutes, Proteobacteria, Bacteroidetes, Actinobacteria, Planctomycetes, Verrucomicrobia and Cyanobacteria. These phyla are known as the major constituents of river autochthonous communities driving key functions of primary production in river ecosystem [57, 58] . Conversely, members of bacterial phylum Actinobacteria are effective players in the degradation of a variety of compounds [59] . However, these phyla were found to be replaced by the selective augmentation of phylum Firmicutes (94.6%) along the mass bathing sites during the course of mass bathing event. It indicated the potential of members of phylum Firmicutes to withstand and flourish in the oxygen-limiting and high-solute conditions [60] . Minimal human activity during the period of lockdown sites featured re-establishment of various bacterial phyla Proteobacteria (70.6%), Bacteroidetes (22.5%), Verrucomicrobia, Actinobacteria and Cyanobacteria having relative abundance of ≥ 1%. It further supported the observed increase in the Shannon index by 0.87-fold during the lockdown period. Similarly, localisation of bacterial communities with respect to the time depicted the predominance of Planococcaceae, followed by Comamonadaceae, Bacillaceae, Flavobacteriaceae, ACK-M1, Moraxellaceae, Verrucomicrobiaceae, Rhodocyclaceae and Oxalobacteraceae before the MBE (relative abundance ≥ 1%). In contrast, enrichment of Planococcaceae (53.1%) (Fig. S2) , Bacillaceae (41.1%) and Sphingomonadaceae (2%) was noted during the course of mass bathing due to its virtue of sustaining the extreme conditions of high-solute concentration and low-oxygen conditions [60, 61] . Furthermore, due to minimal human interference during the lockdown period, augmentation of diverse families, viz. Moraxellaceae, Flavobacteriaceae, [Chromatiaceae], Comamonadaceae, Sphingomonadaceae, Rhodobacteraceae, Pseudomonadaceae, Verrucomicrobiaceae, Oxalobacteraceae and Rhodocyclaceae (relative abundance ≥ 1%), was observed.

At the genus level (relative abundance ≥ 1%), predominance of Bacillus, Lysinibacillus, Flavobacterium, Acinetobacter, Limnohabitans and Planomicrobium was observed as a constituent of autochthonous river microbiota before MBE (Fig. S3) . Conversely, Lysinibacillus (43.8%), Bacillus (41%) and Planomicrobium (4.9%) were found to flourish during the MBE. Moreover, the correlation between the changing environmental factors and bacterial community at the level of genus was carried which revealed a distinct influence of environmental factors on the bacterial communities (Fig. S4) . The variations in the optical density of the water along with the TS, TDS, TSS and BOD were found to favour the bacterial genera Lysinibacillus, Bacillus and Planomicrobium (p ≤ 0.05). However, the dissolved oxygen was found to cast negative impact on the Lysinibacillus, Bacillus and Planomicrobium (p ≤ 0.05). Additionally, the substantial variation in the dissolved oxygen was found to be positively correlated with the Rhodobacter and Novosphingobium (p ≤ 0.05), whereas changes in the estimates of faecal coliforms were found to favour the Lysinibacillus (p ≤ 0.05). Interestingly, the Lysinibacillus, Bacillus and Planomicrobium showed a positive correlation between themselves (p ≤ 0.001); in contrast, they cast negative impact on the bacterial genera, viz. Flavobacterium (p ≤ 0.05), Acinetobacter (p ≤ 0.05), Rheinheimera, Rhodobacter (p ≤ 0.05), Novosphingobium and Limnohabitans. Conversely, Acinetobacter shares a positive correlation between Rhodobacter (p ≤ 0.001) and Novosphingobium (p ≤ 0.01). Similarly, Rhodobacter and Novosphingobium depict a significant positive correlation (p ≤ 0.01). Enumeration of members of Lysinibacillus and Bacillus indicates the potential shedding of human skin-associated microbiota during the MBE. Enrichment of Planomicrobium corroborated the facts of the ability of various Planomicrobium species to thrive in a variety of halophilic environments and their association with recreational waters [28, 60] . Although species of Lysinibacillus represents common members of human skin microbiota, some species of Lysinibacillus have also been associated with the sepsis in immunocompromised patients [61, 62] . Thus, an increase in abundance of Lysinibacillus may pose a serious health risk for the fellow pilgrims bathing in the river. Interestingly, during the lockdown period, re-establishment of few bacterial genera that were observed before the MBE, viz. Acinetobacter, Flavobacterium and Limnohabitans (relative abundance ≥ 1%), was observed.

The relationship between the bacterial communities and the changing physiochemical conditions over the period investigation was derived by applying the linear regression model. . The environmental factors were found to cast negative impact on the bacterial communities during the course of mass bathing [25] . However, the current analysis depicts the positive correlation between above-mentioned factors and the bacterial communities which support composite increase (based on the Shannon index) in river bacterial diversity during the period of lockdown (Fig. 3, Fig. S1 ).

The beta diversity analysis was carried out to determine the divergence in the bacterial communities, under the differential anthropogenic impact across the studied time points. Principal component analysis (PCA) analysis demonstrated three distinct groups based on the temporal variations (i.e. before MGE, during MGE and lockdown) indicating a significant impact of mass bathing on the river bacterial communities. The samples belonging to the before MGE and lockdown period showed lower intra-group variability (Fig. S5) . However, the during MGE samples showed higher intra-group variability which corroborates the fact of samples DS4 and DS5 being sink to other mass bathing sites. Although the Shannon index demonstrates the increase in bacterial diversity during the lockdown period whilst, the distinctive clustering observed during the beta diversity analysis indicates the variation in the bacterial community structure and composition across the studied time points. The latter is supported by the observed changes in bacterial community structure and composition at the level of family and genera (Figs. S2 and S3) . The composite observation of distinct clustering, lower intra-group variability during the lockdown period, significant increase in bacteria diversity (0.87-fold increase in Shannon index) and enumeration of resilient bacterial communities imply a positive impact of lockdown on the river ecosystem and, consequently, on the river bacterial communities [9, 26] .

The attempt to identify the source of diverse bacterial communities using SourceTracker (a Bayesian mixing model) with microbiome profile of healthy Indian subpopulation as a reference dataset revealed a significant increase in skin and stool-associated microbiota in the river during the course of mass bathing event (p ≤ 0.01). mass bathing period (p = 0.004) (Fig. 4a-c) . It is suggestive of the fact that the lockdown had a positive impact, as it indirectly led to cleaning of the river ecosystem. Although health risk associated with such alien communities remains largely unknown however, assessing the source of these microbial communities might aid to devise the mitigation strategies for potential infections related to the eye, skin and bowel [63, 64] . As observed, earlier prevalence of Lysinibacillus during the MBE might have an unfavourable outcome of sepsis [61] . It is in conjunction with the observation from the previous study on sea bathing where a bather itself acts a source of non-point pollution [65] .

Assessment of the functional potential of the river bacterial communities across temporal divergence using an approach of metagenomics imputation depicted substantial variations in the abundance of the gene families involved in the metabolism and genetic information processing. However, gene families related to the infectious diseases and drug resistance were targeted to evaluate the impact of changes in the bacterial communities on the public health. (Fig. 5a) . Likewise, gene families driving the drug resistance in the bacterial communities were also found to rise with the MBE. Specifically, an increase in beta-lactam resistance (1.04-fold), vancomycin resistance (2.49-fold), cationic antimicrobial peptide (CAMP) resistance (2.86-fold) and platinum drug resistance (1.75-fold) was observed during the MBE. The plausible re-establishment of the bacterial communities during the lockdown appeared to have lower prevalence of the gene families associated with drug resistance. A reduction in the levels of beta-lactam resistance (− 0.39-fold), vancomycin resistance (− 0.73-fold), CAMP resistance (− 0.76-fold) and platinum drug resistance (− 0.59-fold) was noted during the lockdown period (Fig. 5b) . Such elevated levels of gene families involved in the infectious diseases and drug resistance during the mass bathing event might pose a serious threat to the attendees which may have catastrophic outcomes like epidemics. Many recent studies on the river ecosystem and/or the religious pilgrimage support the observed increase in the abundance of gene families involved in infectious diseases and drug resistance. Specifically, the study on Ganga targeting the bla NDM-1 gene revealed an increase in antibiotic resistance during the period of pilgrimage [21] . Similarly, in our previous study, we also noted an increase in antibiotic resistance in bacteria isolated across the MBE and described a novel cefotaxime-resistant strain of Corynebacterium godavarianum from the Godavari River, India [66, 67] . The recent study by Yadav et al. [68] involving shotgun metagenomic analysis based on the MinION sequencing supports the observed elevated levels of antibiotic resistance genes during mass bathing.

In summary, a meta-analysis of the river bacterial communities under the extensive anthropogenic impact of mass bathing to minimal impact during the COVID-19 lockdown revealed a significant impact of MBE driving reduction in the bacterial diversity and an increase in potential risk to the human health. Alternatively, the lockdown provides a plausibility of natural restoration of ecosystem as noted by re-establishment of resilient bacterial communities. However, the extent of re-establishment of resilient bacterial communities from the end of previous mass bathing event (year 2015) and before the lockdown period could not be evaluated and could not be ignored.

Our study represents a unique set of baseline data during the lockdown period which will serve as reference for comparative analysis in future studies on the Godavari River.

The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s00248-021-01781-0.

A novel coronavirus from patients with pneumonia in China. NEJM

COVID-19) pandemic. World Health Organization

Immediate psychological responses and associated factors during the initial stage of the 2019 coronavirus disease (COVID-19) epidemic among the general population in China

WHO Coronavirus disease (COVID-19) dashboard

SARS-CoV-2 vaccines: status report

Indirect implications of COVID-19 towards sustainable environment: an investigation in Indian context

A global perspective on the influence of the COVID-19 pandemic on freshwater fish biodiversity

COVID-19 pandemic and associated lockdown as a "global human confinement experiment" to investigate biodiversity conservation

Environmental effects of COVID-19 pandemic and potential strategies of sustainability

Airborne particle levels plummet in Northern India

Effect of restricted emissions during COVID-19 on air quality in India

COVID-19 lockdowns cause global air pollution declines

Evaluation of pollution of Ganga River water using fish as bioindicator

Water quality interaction with alkaline phosphatase in the Ganga River: implications for river health

Anthropogenic activities induce depletion in microbial communities at urban sites of the river Ganges

Impact of lockdown on water quality of river Ganga. Retrieved from: https:// cpcb. nic. in/ upload/ Asses sment-of-Impact-Lockd own-WQ-Major Rivers

COVID-19 lockdown: a ventilator for rivers

Katariy M (2020) 10 things that have happened for the first time in years during the coronavirus lockdown. Scoopwhoop

Central Pollution Control Board (CPCBb) (2020) Report on assessment of impact of lockdown on water quality of river Yamuna-Delhi stretch

Water, sanitation, hygiene and waste disposal practices as COVID-19 response strategy: insights from Bangladesh

Increased waterborne bla NDM-1 resistance gene abundances associated with seasonal human pilgrimages to the Upper Ganges River

Water, sanitation, and hygiene at the world's largest mass gathering

Public health perspectives from the biggest human mass gathering on earth: KumbhMela, India

Mass gatherings medicine: public health issues arising from mass gathering religious and sporting events

World's largest mass bathing event influences the bacterial communities of Godavari, a holy river of India

Cleaning the river Ganga: impact of lockdown on water quality and future implications on river rejuvenation strategies

Cleaning the river Damodar (India): impact of COVID-19 lockdown on water quality and future rejuvenation strategies

The microbiota of recreational freshwaters and the implications for environmental and public health

Contribution of archaea and bacteria in sustaining climate change by oxidizing ammonia and sulfur in an Arctic Fjord

Geochemistry shapes bacterial communities and their metabolic potentials in tertiary coalbed

Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample

QIIME allows analysis of high-throughput community sequencing data

FLASH: fast length adjustment of short reads to improve genome assemblies

Cutadapt removes adapter sequences from highthroughput sequencing reads

Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities

Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB

Search and clustering orders of magnitude faster than BLAST

Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy

phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data

ggplot2: elegant graphics for data analysis

cowplot: streamlined plot theme and plot annotations for 'ggplot2

R package "corrplot": visualization of a correlation matrix (version 0.84)

Hmisc: Harrell Miscellaneous

PerformanceAnalytics: econometric tools for performance and risk analysis

Bayesian community-wide culture-independent microbial source tracking

Salvi S (2020) Gut, oral and skin microbiome of Indian patrilineal families reveal perceptible association with age

The gut microbial diversity of newly diagnosed diabetics but not of prediabetics is significantly different from that of healthy nondiabetics

Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences

Metagenomic biomarker discovery and explanation

A guide to the natural history of freshwater lake bacteria. Microbiol

A dissolved oxygen threshold for shifts in bacterial community structure in a seasonally hypoxic estuary

Tracing particulate matter and associated microorganisms in freshwaters

Attachment of fecal indicator bacteria to particles in the Neuse River Estuary, NC

Population dynamics and ecology of Arcobacter in sewage

Human and environmental impacts on river sediment microbial communities

Actinobacterial 16S rRNA genes from freshwater habitats cluster in four distinct lineages

A unique assemblage of cosmopolitan freshwater bacteria and higher community diversity differentiate an urbanized estuary from oligotrophic Lake Michigan

Characterization of bacterial communities in four freshwater lakes differing in nutrient load and food web structure

Kim BHPlanomicrobium flavidum sp. nov., isolated from a marine solar saltern, and transfer of Planococcus stackebrandtii Mayilraj, et al (2009) 2005 to the genus Planomicrobium as Planomicrobium stackebrandtii comb. nov

Severe sepsis secondary to persistent Lysinibacillus sphaericus, Lysinibacillus fusiformis and Paenibacillus amylolyticus bacteremia

Lysinibacillus timonensis sp. nov., Microbacterium timonense sp. nov., and Erwinia mediterraneensis sp. nov., three new species isolated from the human skin

Guidelines for safe recreational water environments

Elevated bathing-associated disease risks despite certified water quality: a cohort study

Recreational water and infection: a review of recent findings

Corynebacterium godavarianum sp. nov., isolated from the Godavari river

Antimicrobial resistance pattern of microorganisms isolated and identified from Godavari River across the mass gathering event

Comprehensive metagenomic insights into a unique mass gathering and bathing event reveals transient influence on a riverine ecosystem

The authors are grateful to Ms. Vaishali Tile for her help in collecting the samples during the lockdown.Author Contribution K.J., Y.S. and A.S. were involved in the study design. J.B. was involved in the sample collection and analysis of environmental parameters. K.J. and A.S. were involved in the experimental work, bioinformatics and statistical analysis. K.J., S.S., E.I.A., S.A., Y.S. and A.S. were involved in the interpretation of the data. All the authors were involved in the drafting of the manuscript.