Abstract

Background: This study aimed to determine the microbial diversity and Community Structure of Coastal water bodies from some cholera-prone communities in Rivers State. Methods: Water samples were collected from water bodies located in six different cholera prone coastal communities in Rivers state three times each with the determination of coordinates each time using Global Positioning System (GPS). At each collection point, in-situ physicochemical parameters of water sources were determined using H198194, H198195, and H198196 multiparameter meters (Hanna Instrument Inc, USA). The water samples were collected using a sterilized container of 750ml capacity (75cl). DNA was extracted from each water sample collected. The extracted DNA from the water samples were sequenced by Illumina HighSeq by InqabaBiotec, South Africa. The sequencing library was prepared by random fragmentation of the DNA, followed by 5′and 3′ adapter ligation. Adapter-ligated fragments were then PCR-amplified and gel purified. For cluster generation, the library was loaded into a flow cell where fragments were captured on a lawn of surface-bound oligos complementary to the library adapters. Each fragment was then amplified in to distinct, clonal clusters through bridge amplification. When cluster generation was complete, the templates were sequenced using Illumina MiSeq SBS technology. Results: The microbial communities in the coastal water studied are highly diversified. The microorganisms are generally classified into kingdoms containing bacteria, protozoa, viruses, archaea, fungi, and Plantae. For all the coastal waters, kingdom bacteria is common with varying percentage reads which ranged from 45835 reads (46.64%) in Bonny river to 263985 reads (99.71%) in Kaa river. Protozoa were identified in addition to other microbial communities in Andoni and Kaa rivers, both in the South-east senatorial district of Rivers State. The number of viruses detected is more in the Chukwu-Ama river in Ogu/bolo with total reads of 384 (0.25% of microbes present). In addition to other kingdoms present in the Okirika River, fungi were also detected. Alpha diversity of bacterial species across the three regions investigated showed no significant difference while that of individual samples was significantly different at P≤ 0.05 using observed features, Shannon diversity index, Faith’s phylogenetic diversity and Pielou’s evenness tools. Bray-Curtis dissimilarity between the bacterial communities’ structures in the three zones investigated did not reveal distinct differentiation in multivariate space. The observed clustering pattern of samples within the PCoA (Principal coordinate analysis) plot was found to be non-significantly (FDR-adjusted p>0.05) different among the zones in Rivers State (PERMANOVA R2 = 43.3%, p = 0.07). The ASVs detected in all zones of Rivers state spanned 15 phyla (Rivers East = 13; Rivers South-east = 13; Rivers West = 10). Actinobacteriota was the most relatively abundant phyla in both Rivers East (43.9%) and Rivers South-east (35.5%) while Firmicutes (36.2%) was the most abundant phyla in Rivers West. Conclusion: This study has shown that microbial communities present in coastal waters are highly diversified with unknown and unassigned microorganisms detected from all the coastal water examined, showing that these waters are hubs for highly potential microorganisms that could be inimical to human health.

1. Introduction

Cholera is most likely to be found and spread in places with inadequate water treatment, poor sanitation including an improper sewage system, and inadequate hygiene as such places are termed cholera-prone communities. The coastal environment can be immensely biodiverse because the activities of man result in contamination even with the use and expansion along coastlines, notwithstanding, can put weight on these conditions that may restrict the rates of development and availability of living organisms. This can thus lessen the overall biodiversity. Microbial diversity is an essential property of common complexes which decides the working and the origins of a specific indigenous microbial community in a particular environment. Hence, variety analysis turns into the initial phase in portraying the systems. It involves variations in the types of microbes present in an environment and their relative abundance ¹. There have been numerous imaginative, and valuable endeavors being advanced as different strategies to decide the structure and degree of the variety of such environments. While discussing microbial variability, the significant constraint has been most organisms’ non-culturability. Procedures grew at first were reliant upon cultures and thus could deliver just those species which could be developed in the laboratory set-up. The need to catch microorganisms’ historical taxonomy incited the creation of culture-independent procedures like metagenomics ². Notwithstanding, techniques like isolating the DNA, trailed by cloning and sequencing, were extremely bulky albeit valuable in giving extra data about the greater part of the species which couldn't be portrayed before without culture-independent technique. At that point came the innovation which was less work escalated and is getting less expensive step by step—the cutting edge sequencing (Next- generation sequencing). It changed the entire situation and permitted total mechanization with the chance of sequencing entire complexes in a day's time and with costs the same or lesser than the past innovations.

An operational taxonomic unit (OTU) is used to analyze population diversity. It indicates the cluster of similar sequences in a specimen under investigation and defines species. Microbial diversity is measured and expressed as alpha, beta, and gamma diversities at different scales when operational taxonomic units are known. Knowledge of microbial population diversity is critical because studies have shown that disease conditions are linked to microbiome diversity ³.

Phylogenetic distribution-based diversity considers the evolutionary distances and relationships of a variety of microorganisms within the microbial community. It describes the depth of the phylogenetic tree shared by microbes in a given microbial community in a sample. Phylogenetic diversity can be measured using a tool called Unifrac. There are two forms of Unifrac viz Weighed Unifrac and UnweightedUnifrac. While unweighed unifrac utilizes the presence and absence of operational taxonomic units and phylogeny, the weighed Unifrac uses abundance information of OTU and phylogeny ⁴.

The most common organized genetic marker critical in the study of Phylogeny and taxonomy of bacteria is 16SrRNA gene sequencing ⁵. It has been reported that microbial identification at the taxonomic level by 16SrRNA is no longer satisfactory ⁶, thus it becomes important to employ other molecular procedures in order to identify and characterize both culturable and non-culturable microbes at species/ strain levels, such technique is called metagenomics. Metagenomics involves the extraction of microbial genome directly from the specimen without growing the organism(s) in the artificial laboratory medium (culture). It is independent of microbial culture and has helped to resolve the limitations of the genomics study. Metagenome sequencing strives to obtain every genetic information directly from the samples. This technique has been used to characterize microorganisms from both environmental and clinical samples ². All extracted DNA are amplified and sequenced so that the sequence produced are used to determine the taxonomy and the function of all genomes present in a microbial community. By this technique, novel organisms are identified. The reason for metagenomics is to determine microorganisms present and the role they play. Therefore, sequencing alone is not ideal because there are a lot of genes whose functions are not known. The genes help in identifying and studying the organisms more when they are expressed and also aid in enriching the database, hence, providing more opportunities for the development of biotechnology. Metagenomics provides information regarding the functional properties of microbes in microbial communities and their genetic contents. It also reveals microbes not captured by 16SrRNA. The availability of reference genomes in the genebank makes metagenomics studies easier. This study aimed to determine the microbial diversity and Community Structure of Coastal water bodies from some cholera-prone communities in Rivers State.

2. Materials and Methods

Study Design

The research design used in this study is descriptive. The study was done seasonally from February to November 2020. This period covers wet and dry seasons. The reason is that cholera is usually influenced by season, hence seasonal determination of Vibrio species is very crucial. Water samples were collected from coastal communities. The water bodies were situated in coastal communities, two (2) water bodies from each Senatorial District in Rivers State. The water bodies include Marine Base Jetty, Chuku-Ama River, Kaa waterside, Andoni River, Egbolom River, and Bonny River. Water samples were collected in each of these water bodies three times with the determination of coordinates each time using Global Positioning System (GPS). At each collection point, in-situ physicochemical parameters of water sources were determined using H198194, H198195, and H198196 multiparameter meters (Hanna Instrument Inc, USA).The water samples were collected using a sterilized container of 750ml capacity (75cl).

Sample Collection

Sterile plastic containers of 250ml capacity were used to collect water samples from each water source in the six Coastal Communities. Water samples were collected three times from each source. The first collection was done during the dry season at the beginning of the year (February/March). The second collection was done during the middle of the rainy season (June/July) and the third collection was done at the end of the rainy season (October/November 2020). Surface water was collected 300m below the water by tilting and creating a current away from the shore, aseptically to avoid any form of contamination. After the collection, water samples were transported immediately in a cold chain to the laboratory for analysis. The analysis was done seasonally. The water containers were filled up and corked while still under water ⁷.

DNA Extraction from Water Samples

Each 75cl bottle of water sample was duly mixed and dispersed into a corresponding Eppendorf tube of 1.5ml capacity, centrifuged for 5 minutes at 2000rpm. The supernatant was decanted and tubes refilled with corresponding water samples respectively. This process was severally repeated until each water sample became adequately concentrated. Zymo Research bashing bead tubes were arranged and labeled with water sample codes and into each tube was added 250 µl of concentrated water sample and 750µl bashing bead buffer and capped tightly.

Bead beaters filled with a 2ml tube holder were used to secure the tubes and centrifuged at 2000rpm for 5 minutes. Another set of tubes known as zymo-spin III filler in collection tubes were arranged and labeled and 400 µl of supernatant was transferred accordingly into these zymo-spin III-F filter tubes and centrifuge at 8000g for one minute. The filter tubes were discarded leaving the collection tubes. The volume 1200µl of genomiclysis buffer was added to the filtrate on the collection tube.

Another set of tubes called Zymo-spin IC column filled in a collection tubes were arranged and 800µl of the mixture was transferred accordingly into these tubes and centrifuged at 10,000xg for one minute. The flow through from collection was discarded and the step repeated.

200µl of DNA prewashed buffer was added to the Zymo-spin IC column in a new collection tubes and spun for one minute at 10,000xg. The DNA wash buffer was added 500µl into each Zymo-spin IC column and centrifuged again for one minute.

The zymo-spin IC columns were transferred to a clean microcentrifuge tube of 1.5ml capacity and DNA elution buffer was directly added to the column matrix, centrifuged at 10,000xg for 3 minutes to elute the DNA. The elude DNA was transferred to another set of tubes Zymo-spin II-uHRC filter in a clean collection tube, centrifuged at 16,000xg for 3 minutes.

The filtered DNA was quantified and examined for purity using Nanodrop spectrophotometer and then sent to South Africa for metagenomics and biodiversity analysis.

Metagenomics sequencing (Next Generation Sequencing)

The extracted DNA from the water samples were sequenced by Illumina HighSeq by InqabaBiotec, South Africa.

Library Preparation

The sequencing library was prepared by random fragmentation of the DNA, followed by 5′and 3′ adapter ligation. Adapter-ligated fragments were then PCR-amplified and gel purified.

Cluster Generation

For cluster generation, the library was loaded into a flow cell where fragments were captured on a lawn of surface-bound oligos complementary to the library adapters. Each fragment was then amplified in to distinct, clonal clusters through bridge amplification. When cluster generation was complete, the templates were ready for sequencing.

Sequence processing

Method of sequencing was used by applying Illumina MiSeq sequence of processing known as QIIME2 v2018.11 ⁸. In addition, chimeric Sequences, marginal sequence errors and noisy sequences were filtered while picking amplicon sequence variants (ASVs) using DADA2 ⁹. A trained naïve Bayes classifier (SILVA 132) for 16S rRNA V3-V4 hyper variable regions using the q2- feature classifier plugin was carried out using representative sequence taxonomy. Operational Taxonomy Units (OTU) count table was elevated to an even sampling depth of 3845 reads per sample and singleton discarded before determination of alpha-diversity and beta-diversity using the version 2.5-6 Vegan ¹⁰ and R-package (R-Core Team, 2019).

Sequencing was carried out by Illumina SBS technology using a proprietary reversible terminator–based that detected single bases as they were incorporated into DNA-template strands. As all four reversible terminator–bound dNTPs were present during each sequencing cycle, natural competition minimized incorporation bias and greatly reduced raw error rates compared to other technologies.

Statistical Analysis

Statistical analyses for metagenomics study were done using R vs3. 6.1 (R Core Team, 2019). Other parameters such as Alpha diversity differences and physicochemical were determined using the agricolae R package (de Mendiburn, 2020). Permutational multivariate analysis of variance (PERMANOVA) was used to test the significance of the differences in multivariate space; while the comparison of microbial community differences (Beta-diversity) in the three zones of Rivers State was based on Bray-Curtis dissimilarities and Principal Coordinate Analysis (PCoA) using Vegan ¹⁰ R packages. The influence of physicochemical parameters on the community level differentiation was achieved by constrained redundancy analysis (RDA), which was performed using Hellinger transformed bacterial, and environmental data for Vegan R package functions. The significance of the constraining variables was estimated based on a permutation test using the Vegan function.

3. Results

Classification of Microbial Diversity in Coastal Waters

From the metagenomics results obtained, the microbial communities in the coastal water studied are highly diversified. The microorganisms are generally classified into kingdoms (Table 1a) containing bacteria, protozoa, viruses, archaea, fungi, and Plantae. For all the coastal waters, kingdom bacteria is common with varying percentage reads which ranged from 45835 reads (46.64%) in Bonny river to 263985 reads (99.71%) in Kaa river. Protozoa were identified in addition to other microbial communities in Andoni and Kaa rivers, both in the South-east senatorial district of Rivers State. The number of viruses detected is more in the Chukwu-Ama river in Ogu/bolo with total reads of 384 (0.25% of microbes present). In addition to other kingdoms present in the Okirika River, fungi were also detected. There were also microbes of unknown kingdoms detected in each of the coastal water. These also vary but more of such were detected in Chuwu-ama River in Ogu/bolo.

Table 1b showed top Phylum classification of microbial diversity in the coastal waters. The common phyla across the samples are Actinobacteria and Proteobacteria. All the samples showed unknown phyla. Some phyla are also unique to some coastal rivers. Such phyla include Cyanobacteria and Fusobacteria detected in Bonny, Phylum Gemmatimonadetes peculiar to Andoni, and phylum Nitrospira unique to Kaa River. Top-class classification of microbial diversity was also considered in all the water samples (Table 1c). The microbial communities cut across 17 classes with Actinobacteria, Alphaproteobacteria, and Gammaproteobacteriapresent in all the coastal water samples. The samples also demonstrated microorganisms of unknown class. The Abua water contained Verrucomicrobiae and Thermomicrobia unique to it while Bonny River contained microorganisms belonging to classes Clostridia, Bacteroidia,and Fusobacteria. Top ten order was considered and presented in Table 1d with a total of 21 microbial orders identified. The order Actinomycetales is common to all the water samples. The order Bacillales was found in Kaa River, Okirika River, and Bonny River with the highest reads recorded against Kaa River. The order Clostridialeswas seen only in Bonny River. Pseudomonadalesand Flavobacteriales were found only in Abua/Odual River while Pasteurellales, Neisseriales and Lactobacillales were found in Bonny River only. Moreso, Spingobacteriales were seen in Andoni and Okirika while the order Chlamydiales was detected only in the Kaa River. The order Enterobacterialeswas identified in Okirika and Bonny River. At the same time, Rubrobacterales and Caudoviraleswere identified only in theChukwu-ama River.

Table 1e showed the top families of members of microbes from studied coastal waters. A total of thirty-five familieswere identified in addition to the unassigned and unknown families. The top ten families identified include: Nocardioidaceae, Intrasperangiaceae, Acidobacteriaceae,Mycobacteriaceae, Streptomycetaceae, Micrococcaceae, Planctomycetaceae, Methylobacteriaceae, Burkholderiaceae, and Enterobactericeae. Other families of the microbial communities of Medical importance identified are Podoviridae, Spathidiiidae, Pseudomonadaceae, Clostridiaceae, Bacillaceae and Neisseriaceae.

Sequence summary

A total number of both forward and reverse sequence counts was obtained before filtering and clustering into amplicon sequence for each of the water samples with sample CW 25 which has a total number of 271549 sequences each exhibiting the highest number of counts and CW 28 with sequence count of 118664 being the lowest count. After clustering the sequences into different categories based on phylogenetic relationships, the number of amplicon sequence variants (ASV) was obtained for each sample. This ASV is 100% similar for each cluster. The highest number was recorded in sample CW 25 and the lowest was seen in sample CW 28. More so, the frequency of ASV was measured for each sample and sample CW 24 showed the lowest frequency of 3845 ASVs. Rarefaction of the sequences was carried out at the depth of 3845 prior to diversity analysis to remove all bias in the total number of species per sample.

• 34.53% of the ASVs were present in all three regions

• 15.83% were unique to Rivers West

• 19.42% were unique to Rivers South-east

• 13.67% were unique to Rivers East

• 7.19% of the ASVs were present in both Rivers West and Rivers South-east

• 6.47% of ASVs were present in both Rivers South-east and Rivers East

• Only 2.88% of the ASVs were shared between Rivers West and Rivers East

4. Alpha Diversity Analysis

Alpha diversity of bacterial species across the three regions investigated showed no significant difference while that of individual samples was significantly different at P≤ 0.05 using observed features, Shannon diversity index, Faith’s phylogenetic diversity and Pielou’s evenness tools.

The table above showed that pH affected the proportional distribution of species in the coastal waters.

5. Beta-Diversity Analysis

Bray-Curtis dissimilarity between the bacterial communities structures in the three zones investigated did not reveal distinct differentiation in multivariate space (Figure 5). The observed clustering pattern of samples within the PCoA (Principal coordinate analysis) plot was found to be non-significantly (FDR-adjusted p>0.05) different (Table 4a and 4b) among the zones in Rivers State (PERMANOVA R2 = 43.3%, p = 0.07).

The table above showed that the water samples were impacted by the same factor, the non-significant difference on bacterial communities across the regions.

The redundancy analysis (RDA) model for the association between bacterial species and chemical properties (environmental cues) in all three zones of Rivers state was not significant (p>0.05) (Table 4c). Permutational test of the significance of the environmental parameters in the RDA model revealed that only Redox significantly influenced the assembly of bacterial communities (Table 4c). Notably, most individual bacterial phylotypes (genus level) were largely uninfluenced by the measured chemical parameters.

The ASVs detected in all zones of Rivers state spanned 15 phyla (Rivers East = 13; Rivers South-east = 13; Rivers West = 10). Presented in Figure 4.8 are the phyla with at least 1% relative abundance in any of the zones. Actinobacteriota was the most relatively abundant phyla in both Rivers East (43.9%) and Rivers South-east (35.5%) while Firmicutes (36.2%) was the most abundant phyla in Rivers West. Other relatively abundant (>5%) phyla in at least one of the zones includeProteobacteria, Chloroflexi, Planctomycetota, and Acidobacteriota. Notably, Chloroflexi, Acidobacteriota, Patescibacteria, and Gemmatimonadota were not detected in Rivers West unlike in Rivers East and River South-east where the relative abundance of these phyla was relatively high.

At the genus taxonomic rank, Noviherbaspirillum, KD4-96 (Chloroflexi), Nocardioides, Mycobacterium, Rhodococcus, Sphingoaurantiacus, and Sphingomonas were the relatively most abundant (≥3%) genera in Rivers East. In Rivers South-east, KD4-96 (Chloroflexi), Mycobacterium, Sphingomonas, Streptomyces, Rhodococcus, Bacillus, and CandidatusUdaeobacter were the relatively abundant (≥3%) genera whereas, in Rivers West, the most abundant (≥3%) genera were in these order: Mycobacterium > Escherichia-Shigella >Bacillus> Pseudoxanthomonas > Clostridium_sensu_stricto_ 7> Rhodococcus.

Differentially (False Discovery Rate adjusted p<0.05) abundant genera, also referred to as biomarkers for the different zones are presented in Figure 7 and has shown that Acidobacteriota and Chloroflexi were biomarkers in Rivers East and Rivers South-east respectively whereas Firmicutes was biomarker phylum for Rivers West. At the genus taxonomic rank Clostridium, Escherichia-Shigella, and Pseudoxanthomonas were biomarkers in Rivers West while KD4-96 (Chloroflexi) was differentially abundant in Rivers South-east.

Relationship between bacterial communities and environmental parameters

6. Discussion

This study has shown that the microbial communities in coastal waters were highly diversified. The abundance of Actinobacteria and Proteobacteria as the most common phyla in all the coastal water studied is in agreement with the findings of previous works ^{11, 12}. The observations were also made in some coastal water possessing unique phyla such as Cyanobacteria and Fusobacteria in Bonny, Gemmatimonadetes in Andoni, and Nitrospira in Kaa.

Actinobacteria have been reported to play a critical role in the bioremediation of water and soil sites that have been polluted with recalcitrant toxic compounds. It also produces secondary metabolites useful in antibiotic production. Despite these useful functions of Actinobacteria, they pose a serious Public Health challenge as they cause human diseases such as tuberculosis, allergic pneumonia, paratuberculosis, mycetomas, actinomycosis, and various types of abscesses ¹³. The study also agrees with a similar study which reported that Actinobacteria constitute the largest bacteria phyla ubiquitously distributed in both aquatic and soil environments ¹⁴. The highest percentage abundance of Actinbacteria was detected in the Abua River (26.85%). This value is very high compared with that recorded in Bonny which is in the same senatorial district. The difference could be explained by the difference in the Global Positioning System between these two communities. The percentage abundance of Actinobacteria in Ogu/Bolo and Okrika is similar indicating that the two coastal communities have similar activities and environmental factors taking place. Proteobacteriahave three major classes all detected in all the water samples examined in this study. The Alphaproteobacteria which include the Chlamydias and Rickettsia; the Betaproteobacteria which include the Neisseriaspp and the Bordetella pertussis that causes whooping cough; the Gammaproteobacteria which include the Pseudomonasspp, the Pasturellaspp, the Haemophilusducreyi that causes chancroid, the Vibriospecies and intestinal Enterobacteriaceae. All these pathogenic bacteria pose a very serious threat to Public Health ¹⁵.

Fusobacteria causesa wide variety of infections such as Vincent angina, pelvic and postsurgical infection, soft tissue infection, inflammatory bowel disease, Lemierre’s syndrome, and abscess of bone, lungs, and brain. They possess adhesion, slime, and fimbriae as virulence factors. Cyanobacteria are widely present in all aquatic environments because of their ability to source carbon from different environments, possessing neurotoxins, hepatotoxins, and enterotoxins. When contracted through contaminated water cause gastroenteritis. Cyanobacteria also deplete dissolved oxygen when in abundance ¹⁶. This explained why this study reported a low value of dissolved oxygen in Bonny which has the highest abundance of Cyanobacteria.

According to molecular data, Nitrospira are the most diverse nitrifiers that are widely distributed in biological wastewater, and natural ecosystems and yet, are barely studied and uncultured nitrite-oxidizing bacteria ¹⁷. Nitrospira are distributed worldwide in oceans, freshwater habitats, hot spring, and, wastewater treatment plants and is notoriously recalcitrant to isolation via cultural methods ¹⁸. In all the coastal water examined, there was a higher percentage abundance of unknown and unassigned microbial communities. This suggested that the investigated coastal waters are potentially rich hubs for establishing novel microbial strains, and the findings are in tandem with a similar study done in the mangroves region ¹⁹.

In this study, the diversity of microbial communities in the coastal waters located in three senatorial districts in Rivers State was also determined using IllunimaMiSeq sequencing technique. The amplicon sequence variants, the sampling depth, location, and some physicochemical parameters used in shaping microbial communities were also determined. The amplicon sequence variants (ASVs) in this study were 100% similar for each sample, which agrees with Muwawaet al. ¹⁹. The number of ASVs detected in each sample ranged from 61-136. This is far less than the findings of Houet al. (2017) which recorded ASVs range of 659-1835. The difference could be due to different methods used and differences in sampling locations. It could also indicate that the coastal waters studied were disturbed more than that of Hou and co-researchers ²⁰. This study also showed that environmental factors such as redox and pH strongly influence the microbial community structure and distribution of functional groups by shaping the metabolic niches in coastal waters. A similar study recorded temperature as the most environmental contributor to the shaping of microbial community which is in contrast to the findings of this work ²¹.

7. Conclusion

This study used the metagenomics technique to determine the presence of culturable and non-culturable bacteria in surface waters in the coastal area of Rivers State. The microbial communities present in coastal waters are highly diversified. The amplicon sequence variants were obtained following the removal of erroneous sequences generated during PCR amplification and sequencing. The amplicon units obtained in this study shared 100% DNA sequence which was used to classify groups of species, and determine biological and environmental variation as well as ecological patterns.

A lot of unknown and unassigned microorganisms were detected from all the coastal water examined, showing that these waters are hubs for highly potential microorganisms that could be inimical to human health. The relative abundance of Proteobacteria, Actinobacteria, Cyanobacteria, Bacteriodates, Verracomicrobia, and Firmicutes in this study is also a critical pointer showing that the environment is highly polluted with both environmental wastes and human wastes. This study also revealed that the metagenomics method and generated sequences detected the presence of bacteria that are difficult to culture and the existence of highly abundant unknown and unassigned bacteria in water, proving that metagenomics analysis is the better option for detection and identification of microorganisms in a given sample as microbes other than bacteria were also detected.

The microbial diversity also differs from one geographical location to another. This study showed that Acidobacteriota and Chloroflexiwere biomarkers of microbial contamination in Rivers East and Rivers South-east whereas Firmicutes was a biomarker phylum for Rivers West. At the genus taxonomic rank, Clostridium, Escherichia-shigella, and Pseudoxanthomonas were biomarkers of contamination in Rivers West; KD4-96 (Chloroflexi) was differentially abundant in Rivers South-east while Nocardiodes,Mycobacterium, Rhodococcus, Sphingoaurantiacus and Sphingomonas were relatively abundant in Rivers East. It was easy to compare the bacteria structure and diversity in each coastal community and to group them by the use of metagenomic tools.

References