- Research article
- Open Access
High-throughput sequencing of 16S rRNA Gene Reveals Substantial Bacterial Diversity on the Municipal Dumpsite
BMC Microbiology volume 16, Article number: 145 (2016)
Multiple types of solid waste in developing countries is disposed of together in dumpsites where there is interaction between humans, animals and the bacteria in the waste. To study the bacteria at the dumpsite and the associated risks, previous studies have focused on culturable, leaving behind a great number of unculturable bacteria. This study focuses on a more comprehensive approach to study bacteria at the dumpsite. Since the site comprised of unsorted wastes, a qualitative survey was first performed to identify the variety of solid waste as this has influence on the microbial composition. Thus, domestic (Dom), biomedical (Biom), river sludge (Riv), and fecal material of pigs scavenging on the dumpsite (FecD) were sampled. Total DNA was extracted from 78 samples and the v4-16S rRNA amplicons was characterized using an Illumina MiSeq platform.
A total of 8,469,294 sequences passed quality control. Catchall analysis predicted a mean of 8243 species per sample. Diversity was high with an average InvSimpson index of 44.21 ± 1.44. A total of 35 phyla were detected and the predominant were Firmicutes (38 %), Proteobacteria (35 %), Bacteroidetes (13 %) and Actinobacteria (3 %). Overall 76,862 OTUs were detected, however, only 20 % were found more than 10 times. The predominant OTUs were Acinetobacter (12.1 %), Clostridium sensu stricto (4.8 %), Proteinclasticum and Lactobacillus both at (3.4 %), Enterococcus (2.9 %) and Escherichia/Shigella (1.7 %). Indicator analysis (P ≤ 0.05, indicator value ≥ 70) shows that Halomonas, Idiomarina, Tisierella and Proteiniclasticum were associated with Biom; Enterococcus, Bifidobacteria, and Clostridium sensu stricto with FecD and Flavobacteria, Lysobacter and Commamonas to Riv. Acinetobacter and Clostridium sensu stricto were found in 62 % and 49 % of all samples, respectively, at the relative abundance of 1 %. None of OTUs was found across all samples.
This study provides a comprehensive report on the abundance and diversity bacteria in municipal dumpsite. The species richness reported here shows the complexity of this man-made ecosystem and calls for further research to assess for a link between human diseases and the dumpsite. This would provide insight into proper disposal of the waste, as well as, limit the risks to human health associated with the dumpsite.
The amount of solid waste generated has risen due to the increasing urban population in most developing countries [1–3]. Limited resources and inefficient infrastructure prevent proper waste separation leading to waste remaining in their production sites since it does not get transported to the final dumpsite. The lack of proper waste management systems has also created a dumpsite that includes solid waste, such as plastics, organic waste from households, markets, and abattoirs, agricultural waste, industrial waste, and chemical, pharmaceutical, and biomedical waste. Unattended wastes on dumpsites attract insects, birds, and rodents, as well as, domestic and wild animals. Livestock, such as pigs, goats, and cattle scavenge on this waste. Humans, including children also visit the dumpsite to pick waste without protective gear. The unique, man-made ecosystem in the Arusha municipality in Tanzania, representing conditions in other developing cities, has the potential to cause serious impacts to public health.
Studies on the abundance and diversity of bacteria in environments, especially at interfaces like the dumpsite with human and animal activity, would enrich our understanding of the risks from such interfaces. The presence of disease vectors, such as insects, rodents, and several other small wild animals, at dumpsites can potentially spread pathogens from the dumpsite to other habitats. This highlights the importance of a need for a comprehensive study of the abundance and diversity of bacteria in municipal dumpsites.
Previous studies have only focused on the identification of disease vectors and culturable bacteria on the dumpsite. For example, studies by Ahmed et al., and Onyindo et al., [4, 5] report Periplaneta americanus (Cockroach), Musca domestica (House fly), and Ophyra leucostoma (Black garbage fly) and Stomoxys calcitrans (stable fly) the most prevalent disease vectors on the dumpsite. Awisan et al.,  found Staphylococcus aureus, Pseudomonas aerugionosa, Klebsiella pneumonia and Escherichia coli as the aerobic and opportunistic bacteria associated to clinical diseases in Irisan dumpsite, and Emmanuel et al.,  found antibiotic resistant Salmonella spp., Shigella spp., and Vibrio cholerae, Proteus spp., and Pseudomonas spp. from the Utisols dumpsite. However, these methods detect less than 1 % of bacteria found in a particular environment [8, 9]. Other studies have used molecular approaches to identify culturable and unculturable bacteria in a sample, yet, these techniques, such as Sanger sequencing [10, 11], are tedious and inefficient.
In the current study, total genomic DNA was extracted from samples collected from the different types of solid wastes in the dumpsites and the v4 region of the16S rRNA amplicons were sequenced using a high throughput Illumina MiSeq platform. To our knowledge, this is the first report using culture independent approaches and high throughput sequencing to study bacterial abundance and diversity in a municipal dumpsite where interaction between the microbes, animals, disease vectors, and humans is common.
Qualitative survey of the dumpsite
Qualitative survey of the dumpsite revealed the presence of different types of unsorted solid wastes on the dumpsite. Domestic and wild feral animals as well as people were found interacting on the dumpsite. Diversity of animals and solid wastes on the dumpsite is shown in Additional file 1.
Diversities of bacterial communities on dumpsite
A total of 8,469,294 v4 region of 16S rRNA gene sequences of bacteria from 78 solid wastes samples passed all quality control filters. The number of high quality sequences per sample ranged from 329 to 291,482 (mean 108,581, SD 56,553). Catchall analysis of richness predicted an overall mean of 8243 ± 759 species per sample (range 349–21,092, SD 4804). Good’s coverage ranged from 0.9766 to 0.9934 (mean 0.9837, SD 0.0083). Rarefaction curves for different types of solid waste are presented in Additional file 2. The overall diversity of bacteria populations was high with an average inverse Simpson index of 44.21 ± 1.44, Shannon’s evenness average 4.95 ± 0.02, Chao1 richness estimator of 5926 ± 239, and an abundance based coverage estimator (ACE) of 7202 ± 228.
Thirty-five bacterial phyla were detected, however, only seven of these accounted for more than 93 % of all sequences. The predominant phyla were Firmicutes (38 %), Proteobacteria (35 %), Bacteroidetes (13 %), Actinobacteria (3 %), Acidobacteria (2 %), Chloroflexi (2 %) and Spirochaetes (1 %). A large number of phyla were rare and had < 1 % of all sequences (Fig. 1). The unclassified bacteria at phylum level accounted for 2.4 % of all sequences. The heatmap (Additional file 3) shows the abundance and distribution of the predominant bacteria phyla on the dumpsite. When each type of solid waste was separately analysed, Proteobacteria was the most predominant in Biom (37.4 %), Dom (35.7 %) and Riv (42.5 %) solid wastes, while Firmicutes predominated in FecD (59.1 %) only. Majority of the bacterial phyla in all solid waste were rare and contributed < 0.1 % of all sequences. Biom and FecD had only 11 phyla accounted for 98.2 % and 99.7 % of all sequences, respectively. Likewise in Dom and Riv, 12 phyla accounted for 99.6 % and 99.4 % of all sequences, respectively.
A total of 76,862 OTUs were detected, yet, only 20 % (15,272/76,864) were identified more than ten times. At family level; Moraxellaceae (12 %), Clostridiaceae-1 (8 %), Ruminococcaceae (5 %), Lachnospraceae (3.7 %) and Lactobacillaceae (3.4 %) were the most abundant. Unclassified bacteria accounted for about 14 % of all sequences at family level. A total of 1350 genera were detected, of which 16 % (220/1350) were unclassified. Acinetobacter was the most abundant genus accounting for (12.1 %) of all sequences followed by Clostridium sensu stricto (4.8 %), Proteiniclasticum (3.4 %), Lactobacillus (3.4 %), Prevotella (2.6 %), Enterococcus (2.9 %) and Escherichia/Shigella (1.7 %) (Fig. 2).
The core microbiota analyses revealed that none of the OTUs were present across all 78 samples at relative abundance of 1 % or more. Only one OTU, Acinetobacter (assigned sequences 688,875) was found in 62 % (48/78) of samples and Clostridium sensu stricto (assigned sequences 202,494) in 49 % (38/78) of samples. When each type of solid waste was analysed separately, Proteiniclasticum and Acinetobacter were found in 67 % (10/15) of Biom; Clostridium_sensu_stricto in 80 % (16/20) of FecD; Acinetobacter and Proteiniclasticum were found in 73 % (24/33) and 55 % (18/33) of Dom, respectively, and in Riv, 80 % (8/10) of samples had Acinetobacter. Indicator analysis clearly revealed higher affiliation of some bacterial OTUs to specific solid waste on the dumpsite. Thirteen OTUs were significantly associated with Biom (Indicator value ≥ 70 and P ≤ 0.05) amongst which are Halomonas, Alishewanella, and Proteiniclasticum; five were associated to FecD, for example Enterococcus, Bifidobacterium, Clostridium sensu stricto and Cellulosilyticum, and nine were associated with Riv, such as Commamonas, Lysobacter and Flavobacterium (Table 1). None of the OTUs were significantly affiliated with Dom at indicator value ≥ 70 and P ≤ 0.05.
Of the 35 bacterial phyla detected, Metastats revealed 11 phyla significantly different between Dom and FecD. Further, at genus level, out of 1428 differentially abundant genera, 173 were significantly different between Dom and FecD. When Biom and FecD were compared, 8 phyla and 144 genera were significantly different. The Biom and Dom comparison revealed no difference in bacteria community at phylum level (P > 0.05), but 9 bacterial family and 39 genera were significantly different. Of the 1428 genera found in Biom and Dom, 16 % (227/1428) were unclassified. The phylum Lentisphaerae, 66 bacterial family and 180 genera were significantly different between Dom and Riv solid waste (Additional file 4).
Comparison of the bacteria community structure and membership between solid wastes
Phylogenetic tree generated using the Yue & Clayton measures as well as the Jaccard index (Additional file 5) were used in the comparison of the bacterial community structure and membership for different solid wastes. Results of the Parsimony test obtained after the phylogenetic analysis of the Yue and Clayton tree ignoring the branch length revealed a significant difference in bacterial community structure between Dom-FecD (P = 0.011) and Dom-Riv (P = 0.028). There was no difference in community structure between Biom-Dom (P = 0.111), Biom-FecD (P = 0.068), Biom-Riv (P = 0.5240 and FecD-Riv (P = 0.100). When branch length was considered, significantly different structures were found between Dom-Riv (P = 0.001) and FecD-Riv (P = 0.034), while none was detected between Biom-Dom, Biom-FecD, Dom-FecD and Biom-Riv (P > 0.05) using Unweighted UniFrac (Table 2). Further, comparison of the community membership using the phylogenetic tree based on the Jaccard index; the parsimony test revealed that Biom-Dom, Dom-FecD and Biom-Riv (P > 0.05) (Table 3) have the same community membership, while in Biom-FecD (P = 0.016), Dom-Riv (P = 0.002) and FecD-Riv (P = 0.002) were different. When Unweighted UniFrac analysis was performed, only Dom-FecD (P = 0.039) had significantly different community membership. Analysis of molecular variance revealed a significant similarity in bacterial community between Biom-Dom (Yue & Clayton, P = 0.475, Jaccard index, P = 0.012), while the rest of the groups were statistically different (P < 0.008) (Tables 3 and 4).
Further, Fig. 3a and b are the graphic representation of the PCoA plot based on Bray-Curtis distances. The spatial separation between centers of the clouds of the bacteria community structure of different solid waste using Amova have shown statistical difference between Biom-FecD (P < 0.001), Dom-FecD (P = 0.002) and FecD-Riv (P = 0.004), but the same community structure between Biom-Dom (P = 0.012) which is clearly depicted in a PCoA plot constructed from a pool of bacteria community of the same waste type (Fig. 3b). Bacterial OTUs responsible for the difference in clustering of solid waste types were Halomonas, Acinetobacter and Lactobacillus from Biom and Dom solid waste; Enterococcus and Kurthia in FecD, and Lysobacter in Riv.
This study has identified a substantial abundance and diversity of bacteria in a municipal dumpsite in Arusha, Tanzania. The estimated richness (8243 species) per sample and high diversity (InvSimpson index = 44.21) of bacteria on the dumpsite surpasses many of the previous culture based studies. The detection of 76,862 OTUs from 35 bacterial phyla using high throughput sequencing technique gave a more comprehensive look at the bacterial community on the dumpsite. Of the four predominant phyla, Firmicutes was the overall most abundant and includes a variety of gram-positive bacteria. Its particularly dominance in fecal material of pigs is consistent with previous reports by Pajarillo et al.,  as well as in horse fecal microbiota by Costa et al.,  and Shepherd et al., . The predominance of Firmicutes in FecD, especially the genus Clostridium sensu stricto may be due to the feeding habits of pigs on fibre from a variety of unsorted solid waste as reported by Middelbos et al.,  and Yildirim et al.,.
The presence of some bacterial genera exclusively to some types of solid waste, justifies a need to sort and treat solid wastes differently. This would deter a possible genetic material exchange, a process which can result in the emergence and re-emergence of new bacteria of public health importance. The revealed affiliation of Enterococcus to FecD is consistent with previous findings in fecal material of pigs [17–19]. The potential of Enterococcus in nosocomial infections and multidrug resistance [20–22], suggests a need to further examine the antibacterial resistance of isolates from pigs scavenging on dumpsites, and their relation to human and animal pathogens.
The affiliation of Halomonas (1 %) to Biom is consistent with findings in gold mines  and in contaminated heavy metals wastes , as both environments are rich in chemicals. Halomonas metabolize cyanide  and some are reported to carry plasmid pZM3H which confers resistance to different chemicals , leading to their application in soil remediation [25, 26]. The fact that these bacteria are predominant, it would be worthwhile to study their roles in the dumpsite, and examining if they contain plasmids that confer resistance to different chemicals. This could lead into their application in the control of chemical pollutants, especially cyanide at the dumpsite.
Proteobacteria was the predominant second overall phylum, and the most abundant in Biom, Dom and Riv solid wastes. Several studies have linked the predominance of this phylum with human and animals diseases [13, 27, 28]. At the genus level, the overall predominance of Acinetobacter (12 %) on the dumpsite is consistent with report by Saini et al.,  and Hossein et al.,  in solid waste. Previous studies report some species of Acinetobacter associated with human and animal’s diseases. A case example is the existence of Acinetobacter in clinical isolates from intensive care unit , animal with urinary tract infections , and as a causative agent of the nosocomial outbreak in Spain . Its predominance in such extreme environment may be linked to their re-counted capacity to detoxify chemicals as well as their multidrug resistance [33–35] which would support their survival in an environment with diverse chemicals from unsorted solid wastes.
Escherichia and Shigella spp. were among the predominant genera on the dumpsite. These bacteria were also reported in several culture based studies [36–38]. Their capacity to acquire multidrug resistance and their association to human and animal pathogens is well known [39–41]. The abundance of these bacteria within the dumpsite with such extraordinary interaction between animals, humans and microbes may be causing potential health risks. The predominance of these bacteria at the dumpsite underscores a need to study on whether a link exists between bacteria on the dumpsite and the known pathogens. Understanding the potential health risks associated with bacteria from the dumpsite would improve the health of not only the humans interacting with the dumpsite, but with all others that may come in contact with them or other animals.
The significant similarity in bacterial community structure and membership between Biom-Dom, Biom—Riv, and FecD—Riv, implies that similar types of bacteria are found in multiple types of solid wastes. This may be attributed to random disposal of unsorted wastes in the same dumpsite and hence, exchange of bacteria between them. Despite all samples were from the same dumpsite whereby exchange of microbes between solid wastes is possible; distinct affiliation of some bacteria OTUs exclusively to some waste was evident. This phenomenon implies that some bacteria have specific nutritional requirement to survive. It would be interesting to study the bacterial community of sorted solid wastes to examine changes in abundance and diversity.
We report an ever rich and diverse bacterial community in Arusha municipal dumpsite. The species richness reported here shows the complexity of this man-made ecosystem and calls for further research to assess for a link between human diseases and the dumpsite. Understanding the role of the bacteria within the dumpsite and bacteria found within different types of waste will provide insight into proper disposal of the waste, as well as, limit the risks to human health associated with the dumpsite.
Study site and samples
The site for this study was the Arusha municipal dumpsite, where waste from different urban sources is thrown. Sampling was performed in March through June of 2013, whereby prior to sample collection, a qualitative survey was conducted to identify types of most common waste on the dumpsite as these wastes have influence on microbial composition. This comprised waste from households and markets (foods, pampers, clothes, etc.), chemical and biomedical waste (drug containers, used syringes), various plastics and used glassware, waste from abattoirs and brewers, as well as fecal matter from animals scavenging on the dump itself. Samples for this study were the fresh fecal material of pigs scavenging on the dump (FecD, n = 20), domestic solid waste (Dom, n = 33), solid biomedical waste (Biom, n = 15), and run-off water sludge adjoining a nearby river (Riv, n = 8). The core of fresh fecal materials of pigs as well as solid waste and sludge were collected into sterile plastic containers, and, within 1 h, the samples were transported on ice to the laboratory, where total DNA was extracted.
Extraction of total genomic DNA
Total genomic DNA was extracted from about 250 mg of solid waste samples using PowerSoil™ DNA extraction kit (MOBIO Laboratories, Carlsbad, CA) according to the manufacturer’s protocol. Quality and quantity of total DNA was verified with a NanoDrop ND-2000c spectrophotometer (Thermo Scientific) and gel electrophoresis run in 0.8 % agarose and visualized by ultraviolet illumination after staining with gel red™. The DNA was stored at -20oC until further processing.
16S rRNA amplification, Library Construction and Sequencing
The Illumina sequencing preparation guide  was used to prepare a pooled amplicon of the v4 region of 16S rRNA gene for sequencing. Primers (515 F/806R) designed for v4 region of 16S rRNA and protocols were adapted from Caporaso . Duplicate reactions were done in PCR master mix reaction in 20 μl AccuPower® Taq PCR PreMix composed of 0.5 μl of 10pmol/μl each for the forward and reverse primers, 17 μl molecular grade water, and 2 μl DNA template. The PCR program was run on GeneAMP™ PCR system 9700 set at 95 oC for 3 min, 35 cycles of 94 oC for 45 s, 50 oC for 60 s and 72 oC for 90 s and a final extension at 72 oC for 10 min. Amplicon quality was visualized using gel electrophoresis, and then pooled and purified using QIAquick® PCR purification kit (Qiagen, German) following manufacturer’s protocol. Purified PCR products were normalized to 120 ng. DNA was quantified using Qubit® dsDNA assay kit in Qubit fluorometer 2.0 (Invitrogen, Life Technologies) and the quality was assessed using Agilent DNA 1000 Chip in Agilent 2100 Bioanalyzer (Agilent Technologies, Waldbronn, Germany.
Library denaturing, dilution, and PhiX control preparation was done as described in the 16S metagenomic sequencing library preparation guide . Libraries were denatured and primers were used according to the method described in Caporaso . Sequencing of the library was performed with the Illumina MiSeq platform (San Diego, USA) using 2 × 250 paired- end chemistry at the BecA—ILRI Hub genomic platform, Nairobi, Kenya.
Sequence data analysis and statistics
The Mothur package algorithms (v1.34.1) were used for both quality control and sequence data analysis . After paired end reads were assembled, sequences were aligned with the Silva 16S rRNA reference database (www.arb-silva.de) . Sequences that were < 239 bp and > 260 bp in length, contained >2 ambiguous base calls or long runs (>8 bp) of homopolymers, or did not align with the correct region were removed. Chimeras were identified using Uchime  and eliminated. Taxonomy was assigned using the RDP taxonomy database (http://rdp.cme.msu.edu/index.jsp) . Sequences were binned into operational taxonomic units (OTUs) at 97 % sequence similarity cut-off.
Species richness was assessed with Chao1  and abundance based coverage estimator ACE  while evenness and diversity of species were estimated by Shannon , jackknife  and inverse Simpson  indices as well as catchall analysis . All analyses were performed using built-in commands in Mothur v1.34.1. Rarefaction analyses were done at a maximum of 97 % sequence similarity cut-off and was plotted using Phyloseq package  in R version 3.1.2. In order to compare bacterial populations between different solid wastes in the same municipal dumpsite, subsampling of sequences from different wastes was done to normalize them for efficient comparison . This consisted of random selection of a number of sequences from each sample consistent to the lowest abundance in all samples. The community membership was compared using the traditional Jaccard index, while population structure was assessed using the Yue & Clayton measure of dissimilarity. Dendrograms were created using Mothur to compare the similarity of bacterial populations among all sample types using both Jaccard index and Yue & Clayton measure which account for the relative abundances in each sample. Figures were generated by FigTree v1.4.2. .
To check if the bacterial communities differed significantly between solid wastes, the parsimony  and Unweighted UniFrac  tests were done. The MOTHUR commands “parsimony” and “Unifrac.Unweighted”, respectively, were applied to the Jaccard and the Yue & Clayton OTU based phylogenetic tree. The statistical significance of the difference in genetic diversity of bacteria community within each solid waste type from the average genetic diversity of both communities pooled together was also assessed using Analysis of Molecular Variance (AMOVA) .
The core microbiota analysis was performed in Mothur with command “get.coremicrobiome” and it consisted of identification of OTUs (at relative abundance of 1 %) present in all samples when pooled together and also when each type of solid waste was individually analysed. Indicator analysis  was used to test for possible OTUs affiliated to different types of solid waste. Indicator values (IV) ranged from 1 to 100 with higher values for stronger indicators. Though literature considers indicator values > 30 and P-value ≤ 0.05 as good indicators , in this study only OTUs with indicator values ≥ 70 and P ≤ 0.05 were judged as having strong affiliation to particular solid waste. The Metastats program  was used to identify statistically different OTUs among solid wastes. The shared OTUs file, consensus taxonomy file, and metadata file generated in Mothur v.1.34.1 were imported into METAGENassist  where visualization using a heatmap, bar charts and PCA plots was done get more insights on the nature of bacteria present on dumpsite. The BIOM file generated in Mothur was imported into MEGAN5 v5.5.3  where Principal coordinate analysis and relative abundance of different taxa were visualized. A p-value of ≤ 0.05 was considered significant for all comparisons.
ABCF, African Biosciences Challenge Fund; BecA—ILRI, Biosciences eastern and central Africa—International Livestock Research Institute; MEGAN, MetaGenome ANalyser; OTUs, operational taxonomic units
Mungure MJ. Governance and community participation in Municipal Solid Waste management, case of Arusha and Dar es Salaam Tanzania. M.Sc. In Environmental Management. Aalborg, Denmark: Aalborg University; 2008.
Getahun T, Mengistie E, Haddis A, Wasie F, Alemayehu E, Dadi D, Van Gerven T, Van der Bruggen B. Municipal solid waste generation in growing urban areas in Africa: current practices and relation to socioeconomic factors in Jimma, Ethiopia. Environ Monit Assess. 2012;184(10):6337–45.
Guerrero LA, Maas G, Hogland W. Solid waste management challenges for cities in developing countries. Waste Manag. 2013;33(1):220–32.
Ahmed A. Insect Vectors of Pathogens in Selected Undisposed Refuse Dumps in Kaduna Town, Northern Nigeria. Sci World J. 2012;6(4):21–26.
Onyido AE, Azubuike J, Amadi ES, Obiukwu MO, Ozumba NA, Ikpeze OO. A survey of public health disease vectors breeding in Refuse Dumps in Onitsha Metropolis, Anambra State Nigeria. N Y J Sci. 2011;4(9):34–9.
Awisan J, Demicais Jr V, Badiang A, Cariño L, Ferrer DM. Aerobic pathogenic and opportunistic bacteria of public health significance in Irisan dumpsite soil. 2013.
Emmanuel I, Joseph N, Kingsley E-I, Egbebor EM, Lawrence E. Antibiotic susceptibility profiles of enteric bacterial isolates from dumpsite utisols and water sources in a rural community in Cross River State, Southern Nigeria. Nat Sci. 2011;9:46–50.
Vandecandelaere I, Matthijs N, Van Nieuwerburgh F, Deforce D, Vosters P, De Bus L, Nelis HJ, Depuydt P, Coenye T. Assessment of microbial diversity in biofilms recovered from endotracheal tubes using culture dependent and independent approaches. PLoS One. 2012;7(6):e38401.
Simon C, Daniel R. Metagenomic analyses: past and future trends. Appl Environ Microbiol. 2010;77(4):1153–61.
Huang LN, Zhu S, Zhou H, Qu LH. Molecular phylogenetic diversity of bacteria associated with the leachate of a closed municipal solid waste landfill. FEMS Microbiol Lett. 2005;242(2):297–303.
Chen Y, Dumont MG, Cébron A, Murrell JC. Identification of active methanotrophs in a landfill cover soil through detection of expression of 16S rRNA and functional genes. Environ Microbiol. 2007;9(11):2855–69.
Pajarillo EAB, Chae JP, Balolong MP, Kim HB, Seo K-S, Kang D-K. Characterization of the fecal microbial communities of duroc pigs using 16S rRNA gene pyrosequencing. Asian Australas J Anim Sci. 2015;28(4):584–91.
Costa MC, Arroyo LG, Allen-Vercoe E, Stämpfli HR, Kim PT, Sturgeon A, Weese JS. Comparison of the fecal microbiota of healthy horses and horses with colitis by high throughput sequencing of the V3-V5 region of the 16S rRNA gene. PLoS One. 2012;7(7):e41484.
Shepherd ML, Swecker Jr WS, Jensen RV, Ponder MA. Characterization of the fecal bacteria communities of forage-fed horses by pyrosequencing of 16S rRNA V4 gene amplicons. FEMS Microbiol Lett. 2012;326(1):62–8.
Middelbos IS, Boler BMV, Qu A, White BA, Swanson KS, Fahey Jr GC. Phylogenetic characterization of fecal microbial communities of dogs fed diets with or without supplemental dietary fiber using 454 pyrosequencing. PLoS One. 2010;5(3):e9768.
Yildirim S, Yeoman CJ, Sipos M, Torralba M, Wilson BA, Goldberg TL, Stumpf RM, Leigh SR, White BA, Nelson KE. Characterization of the fecal microbiome from non-human wild primates reveals species specific microbial communities. PLoS One. 2010;5(11):e13963.
Dang STT, Van Truong D, Madsen H, Dalsgaard A. Survival of faecal indicator bacteria in treated pig manure stored in clay-covered heaps in Vietnam. Vet Microbiol. 2011;152(3–4):374–8.
Novais C, Freitas AR, Silveira E, Antunes P, Silva R, Coque TM, Peixe L. Spread of multidrug-resistant Enterococcus to animals and humans: an underestimated role for the pig farm environment. J Antimicrob Chemother. 2013;68(12):2746–54. dkt289.
Freitas AR, Coque TM, Novais C, Hammerum AM, Lester CH, Zervos MJ, Donabedian S, Jensen LB, Francia MV, Baquero F. Human and swine hosts share vancomycin-resistant Enterococcus faecium CC17 and CC5 and Enterococcus faecalis CC2 clonal clusters harboring Tn1546 on indistinguishable plasmids. J Clin Microbiol. 2011;49(3):925–31.
Arias CA, Murray BE. The rise of the Enterococcus: beyond vancomycin resistance. Nat Rev Microbiol. 2012;10(4):266–78.
Lebreton F, van Schaik W, McGuire AM, Godfrey P, Griggs A, Mazumdar V, Corander J, Cheng L, Saif S, Young S: Emergence of epidemic multidrug-resistant Enterococcus faecium from animal and commensal strains. mBio. 2013;4(4):e00534-00513.
Marcadé G, Micol J-B, Jacquier H, Raskine L, Donay J-L, Nicolas-Viaud S, Rouveau M, Ribaud P, Dombret H, Leclercq R. Outbreak in a haematology unit involving an unusual strain of glycopeptide-resistant Enterococcus faecium carrying both vanA and vanB genes. J Antimicrob Chemother. 2013;376:1–6.
Lin Y, Fan H, Hao X, Johnstone L, Hu Y, Wei G, Alwathnani HA, Wang G, Rensing C. Draft genome sequence of Halomonas sp. strain HAL1, a moderately halophilic arsenite-oxidizing bacterium isolated from gold-mine soil. J Bacteriol. 2012;194(1):199–200.
Dziewit L, Pyzik A, Matlakowska R, Baj J, Szuplewska M, Bartosik D. Characterization of Halomonas sp. ZM3 isolated from the Zelazny Most post-flotation waste reservoir, with a special focus on its mobile DNA. BMC Microbiol. 2013;13(1):59.
Khamar Z, Makhdoumi A, Mahmudy Gharaie MH. Remediation of cyanide from the gold mine tailing pond by a novel bacterial co-culture. Int Biodeterioration Biodegradation. 2015;99:123–128.
Achal V, Pan X, Zhang D. Bioremediation of strontium (Sr) contaminated aquifer quartz sand based on carbonate precipitation induced by Sr resistant Halomonas sp. Chemosphere. 2012;89(6):764–8.
Suchodolski JS, Markel ME, Garcia-Mazcorro JF, Unterer S, Heilmann RM, Dowd SE, Kachroo P, Ivanov I, Minamoto Y, Dillman EM. The fecal microbiome in dogs with acute diarrhea and idiopathic inflammatory bowel disease. PLoS One. 2012;7(12):e51907.
Suchodolski JS, Xenoulis PG, Paddock CG, Steiner JM, Jergens AE. Molecular analysis of the bacterial microbiota in duodenal biopsies from dogs with idiopathic inflammatory bowel disease. Vet Microbiol. 2010;142(3):394–400.
Saini S, Das BK, Kapil A, Nagarajan SS, Sarma R. The study of bacterial flora of different types in hospital waste: evaluation of waste treatment at AIIMS Hospital, New Delhi. Southeast Asian J Trop Med Public Health. 2004;35(4):986.
Hossain MS, Rahman NNNA, Balakrishnan V, Puvanesuaran VR, Sarker MZI, Kadir MOA. Infectious risk assessment of unsafe handling practices and management of clinical solid waste. Int J Environ Res Public Health. 2013;10(2):556–67.
Omer MI, Gumaa SA, Hassan AA, Idris KH, Ali OA, Osman MM, Saleh MS, Mohamed NA, Khaled MM. Prevalence and resistance profile of acinetobacter baumannii clinical isolates from a private hospital in Khartoum, Sudan. Am J Microbiol Res. 2015;3(2):76–9.
Pomba C, Endimiani A, Rossano A, Saial D, Couto N, Perreten V. First report of OXA-23-mediated carbapenem resistance in sequence type 2 multidrug-resistant acinetobacter baumannii associated with urinary tract infection in a cat. Antimicrob Agents Chemother. 2014;58(2):1267–8.
Mosqueda N, Espinal P, Cosgaya C, Viota S, Plasensia V, Álvarez-Lerma F, Montero M, Gómez J, Horcajada JP, Vila J. Globally expanding carbapenemase finally appears in Spain: nosocomial outbreak of Acinetobacter baumannii producing plasmid-encoded OXA-23 in Barcelona, Spain. Antimicrob Agents Chemother. 2013;57(10):5155–7.
Franzetti A, Gandolfi I, Fracchia L, Van Hamme J, Gkorezis P, Marchant R, Banat IM. Biosurfactant use in heavy metal removal from industrial effluents and contaminated sites. Biosurfactants. 2014;159:361.
Zarrilli R, Pournaras S, Giannouli M, Tsakris A. Global evolution of multidrug-resistant Acinetobacter baumannii clonal lineages. Int J Antimicrob Agents. 2013;41(1):11–9.
Rastogi V, Rastogi P, Bhatia S. Bacteriological profile of biomedical waste: management guidelines. J Indian Acad Forensic Med. 2011;33(2):145–8.
Anitha J, JAYRAAJ IA. Isolation and identification of bacteria in biomedical wastes (BMW). Int J Pharm Pharmaceutical Sci. 2012;4(5):386–8.
Oyeleke S, Istifanus N: The microbiological effects of hospital wastes on the environment. Afr J Biotechnol. 2009;8(7):1253–1257.
Doumith M, Dhanji H, Ellington MJ, Hawkey P, Woodford N. Characterization of plasmids encoding extended-spectrum β-lactamases and their addiction systems circulating among Escherichia coli clinical isolates in the UK. J Antimicrob Chemother. 2012;67(4):878–85.
García A, Fox JG, Besser TE. Zoonotic enterohemorrhagic Escherichia coli: a One Health perspective. Ilar J. 2010;51(3):221–32.
Fritah S, Lhocine N, Golebiowski F, Mounier J, Andrieux A, Jouvion G, Hay RT, Sansonetti P, Dejean A. Sumoylation controls host anti‐bacterial response to the gut invasive pathogen Shigella flexneri. EMBO Rep. 2014;15(9):965–72.
Illumina: 16S Metagenomic Sequencing Library preparation Guide: (http://support.illumina.com/downloads/16s_metagenomic_sequencing_library_preparation.html). 2013. Accessed 6 Mar 2014.
Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, Fierer N, Knight R. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci. 2011;108(Supplement 1):4516–22.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glöckner FO. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41(D1):D590–6.
Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–200.
Cole JR, Wang Q, Fish JA, Chai B, McGarrell DM, Sun Y, Brown CT, Porras-Alfaro A, Kuske CR, Tiedje JM. Ribosomal Database Project: data and tools for high throughput rRNA analysis. Nucleic Acids Res. 2014;42(D1):D633–42.
Chao A. Nonparametric estimation of the number of classes in a community. Scand J Stat. 1984;11:265–70.
Chao A, Lee S-M. Estimating the number of classes via sample coverage. J Am Stat Assoc. 1992;87(417):210–7.
Chao A, Shen T-J. Nonparametric estimation of Shannon’s index of diversity when there are unseen species in sample. Environ Ecol Stat. 2003;10(4):429–43.
Heltshe JF, Forrester NE. Estimating species richness using the jackknife procedure. Biometrics. 1983;39(1):1–11.
Hunter PR, Gaston MA. Numerical index of the discriminatory ability of typing systems: an application of Simpson’s index of diversity. J Clin Microbiol. 1988;26(11):2465–6.
Bunge J. Estimating the number of species with CatchAll. In: Pacific Symposium on Biocomputing: 2011. World Sci. 2011;PBS 2011:121–130.
McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8(4):e61217.
Gihring TM, Green SJ, Schadt CW. Massively parallel rRNA gene sequencing exacerbates the potential for biased community diversity comparisons due to variable library sizes. Environ Microbiol. 2012;14(2):285–90.
Rambaut A: FigTree v1. 3.1: Tree figure drawing tool. Website: http://tree bio ed ac uk/software/figtree. 2009. [Accessed 10 Jul 2015].
Schloss PD, Handelsman J. Introducing TreeClimber, a test to compare microbial community structures. Appl Environ Microbiol. 2006;72(4):2379–84.
Lozupone CA, Hamady M, Kelley ST, Knight R. Quantitative and qualitative β diversity measures lead to different insights into factors that structure microbial communities. Appl Environ Microbiol. 2007;73(5):1576–85.
Michalakis Y, Excoffier L. A generic estimation of population subdivision using distances between alleles with special reference for microsatellite loci. Genetics. 1996;142(3):1061–4.
Dufrêne M, Legendre P. Species assemblages and indicator species: the need for a flexible asymmetrical approach. Ecol Monogr. 1997;67(3):345–66.
White JR, Nagarajan N, Pop M. Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput Biol. 2009;5(4):e1000352.
Arndt D, Xia J, Liu Y, Zhou Y, Guo AC, Cruz JA, Sinelnikov I, Budwill K, Nesbø CL, Wishart DS. METAGENassist: A comprehensive web server for comparative metagenomics. Nucleic Acids Res. 2012;40(Web Server issue):W88–95. gks497.
Mitra S, Stärk M, Huson DH. Analysis of 16S rRNA environmental sequences using MEGAN. BMC Genomics. 2011;12 Suppl 3:S17.
We acknowledge the Government of Tanzania through The Nelson Mandela African Institution of Science and Technology for sponsoring this study; Africa Biosciences Challenge Fund (ABCF) for a fellowship to execute laboratory work at Biosciences eastern and Central Africa—International Livestock Research Institute (BecA—ILRI hub) in Nairobi, Kenya. Dr. Rob Skilton for official facilitation while at BecA-ILRI hub, Dr. Francesca Stomeo for her assistance with sequencing, Mr Alan Orth for linking us to ILRI HPC and the Genome Science Centre at Sokoine University of Agriculture for assistance with microbiology and sampling techniques are highly acknowledged.
This study was funded by the government of Tanzania through the Commission for Sciences and Technology (COSTECH) and the Nelson Mandela African Institute of Science and Technology. Sequencing of DNA and data analyses was partly sponsored by the African Biosciences Challenge Fund (ABCF).
Availability of data and materials
Sequence data discussed in this work have been deposited at the NCBI sequence reads archive (SRA) and assigned with accession numbers. Sequences of bacteria from solid biomedical waste—Biom [GenBank: SRP045926], sequences of bacteria from river sludge near dumpsite—Riv [GenBank: SRP045810], sequences of bacteria from domestic solid waste—Dom [GenBank: SRP046216] and sequences from fecal material of pigs scavenging on dumpsite—FecD [GenBank: SRP046287].
KSM, PSG and VK designed the study; PSG Supervised the study, KSM performed sample collection and DNA extraction; KSM and SM performed sequence library preparation, sequencing, KSM, AS, MS performed data analysis; KSM, SM, AS, MS, VK and PSG were responsible for the draft and final version of the manuscript. All authors read and approved the final draft of the manuscript.
The authors declare no conflict of interest.
Consent for publication
Ethics approval and consent to participate
This study was approved by the research committee of the Nelson Mandela African Institution of Science and Technology. Permits to sample animals and variety of solid wastes on the dumpsite were granted by the Arusha District Veterinary office. Permits to transfer samples between laboratories were given by the Zoosanitary inspectorate services of Tanzania, Arusha (VIC/AR/ZIS/0345) and Veterinary Services under the Ministry of Agriculture Livestock and fisheries of Kenya (RES/POL/VOL.XXIV/506).
Showing different animals and solid wastes on the dumpsite. (DOCX 13 kb)
Rarefaction analysis curves of different solid waste from the dumpsite. (DOCX 129 kb)
Heatmap showing abundance and distribution of predominant bacteria phyla from different solid waste samples. (DOCX 125 kb)
Tables showing significantly different bacteria at genus level between different types of wastes solid. (DOCX 55 kb)
OTU based dendrogram demonstrating the similarity and differences of bacteria communities’ structure and membership from different solid wastes. Results were obtained using the Yue & Clayton measure (Fig a) and the Jaccard index (Fig b). (DOCX 1469 kb)
About this article
Cite this article
Mwaikono, K.S., Maina, S., Sebastian, A. et al. High-throughput sequencing of 16S rRNA Gene Reveals Substantial Bacterial Diversity on the Municipal Dumpsite. BMC Microbiol 16, 145 (2016). https://doi.org/10.1186/s12866-016-0758-8
- Municipal dumpsite
- Solid waste
- Molecular diversity
- 16S rRNA
- Illumina MiSeq