3.1 Estimation of richness and diversity
After quality filtering and removal of chimeric sequences, a total of 203140 high-quality reads were obtained: 36405 for C11, 31108 for C13, 35767 for C16, 32924 for C18, 34251 for LT and 32685 for LT1, with an average length of 419 bp (Online Resource Table 2). The sequences clustered into 2581–5739 OTUs per sample when applying a 97% similarity threshold for inclusion.
The rarefaction analysis of the six leachate samples at an OTU cut-off value of 97% indicated that a full census had not been achieved for any of the samples (appendix Fig. 4), all the curves were far from reaching a plateau, suggesting that the curves largely underestimate the numbers of different types of bacteria in each sample, which did not allow for calculating the total species richness, the description of biological diversity only based on all the identified OTUs.
A series of coverage estimators were computed, including ACE, Chao1, the Good’s coverage, and the Shannon index (Table 2). Good’s estimator indicated nearly complete coverage. The sequence composition analyses revealed that the bacterial diversity estimates for the LFL were much greater than previously published descriptions of LFL diversities. The average number of OTUs for all samples was 3680, which was much greater than that found in previous studies. For example, Song, Wang, Tang, and Lei [5] investigated the bacterial communities of ten LFL samples from five landfill sites in China and found 709–1599 OTUs per sample. The community richness (Chao 1 and ACE) showed the highest for C13 and the lowest for LT1, which correlates with the number of observed OTUs in the samples; the community diversity (Shannon’s H index which considers the evenness of OTU distribution) showed the highest for C13 and the lowest for C16.
Table 2
Coverage estimates for the bacterial LFL communities
Sample | OTU richness | ACE* | Chao1 | Good’s coverage | Shannon |
C11 | 3545 | 13553 | 13226 | 0.983 | 4.225 |
C13 | 5739 | 28221 | 24135 | 0.986 | 4.647 |
C16 | 2960 | 13038 | 11812 | 0.986 | 3.251 |
C18 | 4234 | 16480 | 15116 | 0.985 | 4.019 |
LT | 3019 | 12953 | 11547 | 0.985 | 3.481 |
LT1 | 2581 | 12530 | 10513 | 0.985 | 3.263 |
*ACE abundance-based coverage estimator |
3.3 Microbial community structure
The comparative analysis using the RDP (Ribosomal Database Project) database showed that only 14 to 105 of the sequence tags (0.04 to 0.29%) were affiliated with the domain Archaea. Due to their extremely low abundance, these sequences were excluded from the community structure analysis.
For all six samples (four from landfill cells and one from LFL tank, the most dominant phylum was Proteobacteria (Pseudomonadota) (56–71%, Fig. 2A and Table 3). This was also reflected in the evolutionary trees of the top 100 genera, where the most abundant genera belonged to the phylum Proteobacteria (Pseudomonadota) (appendix Fig. 5). The next most abundant phyla were Bacteroidetes (Bacteroidota) (7 to 21%), Actinobacteria (Actinomycetota) (4 to 9%), Firmicutes (Bacillota) (5 to 8%) and Balneolaeota (Balneolota) (2 to 6%), as shown in Fig. 2A. The remaining phyla that featured in the top 10 abundances for each sample included Spirochaetes (Spirochaetota), Cyanobacteria, Tenericutes, Euryarchaeota, Gemmatimonadetes, Verrucomicrobia (Verrucomicrobiota), Acidobacteria (Acidobacteriota), Synergistetes (Synergistota) and Chloroflexi (Chloroflexota), however, it should be noted that their individual contribution only varied between 0.15 to 2.07% of the bacterial sequence set of each sample (Table 3). The microbial community composition of the top 10 proportions of each sample from superkingdom (L1) to genus (L6) is shown in Online Resource Table 5–8.
Table 3
HiSeq sequencing results (Top 10 Hits Distribution)
Sample | Acido. | Actino. | Bact. | Baln. | Cyano. | Chlo. | Eury. | Firm. | Gemma. | Proteo. | Spir. | Syne. | Tener. | Verr. |
C11 | - | 2689 | 2527 | 2304 | 248 | - | 105 | 1902 | 96 | 25268 | 752 | - | 130 | - |
C13 | 170 | 1927 | 2312 | 1093 | 344 | - | - | 1655 | 396 | 22092 | 335 | - | - | 253 |
C16 | - | 2073 | 7216 | 784 | 54 | - | - | 1855 | 235 | 22745 | 230 | 155 | 70 | - |
C18 | - | 3113 | 2384 | 2111 | 319 | 93 | - | 2436 | 72 | 21554 | 402 | 78 | - | - |
LT | - | 1251 | 8058 | 1005 | 391 | 105 | - | 2636 | 209 | 19319 | 638 | 124 | - | - |
LT1 | - | 2210 | 6984 | 871 | 121 | 66 | - | 1537 | 186 | 19763 | 472 | 100 | - | - |
AcidoAcidobacteria (Acidobacteriota), Actino Actinobacteria (Actinomycetota), Bact Bacteroidetes (Bacteroidota), Baln Balneolaeota (Balneolota), Cyano Cyanobacteria, Chlo Chloroflexi (Chloroflexota), |
Eury Euryarchaeota, Firm Firmicutes (Bacillota), Gemma Gemmatimonadetes (Gemmatimonadota), Proteo Proteobacteria (Pseudomonadota), Spir Spirochaetes (Spirochaetota), Syne Synergistetes (Synergistota), Tener Tenericutes, Verr Verrucomicrobia (Verrucomicrobiota)
The result of OTUs clustering, where both the common and unique information for different samples were analysed is represented in a flower diagram (Fig. 3). This analysis indicated that all samples shared 292 core OTUs. The majority of these OTUs belong to Betaproteobacteria, Gammaproteobacteria, Alphaproteobacteria and Bacteroidetes (Bacteroidota), reflecting the predominance of the Proteobacteria (Pseudomonadota) and Bacteroidetes (Bacteroidota), also observed during whole community analysis (Table 3). In addition, the number of unique OTUs in C11, C13, C16, C18, LT and LT1 were 63, 140, 46, 62, 46 and 32, respectively (Fig. 3).
Each petal in the flower diagram represents a sample. The core number in the centre is for the number of OTUs present in all samples, while the number in the petals for the unique OTUs only showing in each sample.
Active landfill (C16, C18<) and the closed landfill (C11&C13)
Spp. from the phylum Proteobacteria (Pseudomonadota) occupied a significant proportion (up to 71%) of the communities present in both C11 and C13 (landfill cell 11 and landfill cell 13, both closed landfill sites), followed by Bacteroidetes (Bacteroidota) (up to 7%), Actinobacteria (Actinomycetota) (up to 7%), Firmicutes (Bacillota) (up to 5%) and Balneolaeota (Balneolota) (up to 6%) (Fig. 2A and Online Resource Table 4). Samples C16, C18 and LT were sourced from active landfill sites, while C18 was newly filled before sampled, sample LT came from the leachate tank that collected leachate from all active landfill sites. However, the new filled waste caused the microbial composition of C18 to be very different from C16 and LT. Interestingly, the most predominant phyla in C18 were Proteobacteria (Pseudomonadota), Actinobacteria (Actinomycetota), Bacteroidetes (Bacteroidota), Firmicutes (Bacillota) and Balneolaeota (Balneolota), and the most predominant classes were Betaproteobacteria and Gammaproteobacteria, Actinobacteria and Alphaproteobacteria (Table 3), which showed similar bacterial composition with C11 and C13. For samples C16 and LT, the top five phyla were similar to those found in C11 / C13, however, the ratio of each phylum in these samples (C16 and LT) showed a significant difference. The proportion of phylum Proteobacteria (Pseudomonadota) decreased from 69% (C11) / 71% (C13) to 63% (C16) / 56% (LT) mainly caused by the dramatic decreasing of class Betaproteobacteria, Gammaproteobacteria and Alphaproteobacteria, the representation of class Betaproteobacteria dropped from the highest 33.61% (C13) to the lowest 9.4% (C16), as shown in Online Resource Table 8. Meanwhile, it is worth noting that the proportion of class Epsilonproteobacteria in the samples of active landfill sites (C16 and LT) has increased significantly compared with that of closed landfill sites (C11 and C13), class Epsilonproteobacteria has become the most dominant class in C16 and LT. Due to this behaviour, the total Proteobacteria (Pseudomonadota) in samples C16 and LT could be able to remain the most dominant phylum. Contrary to the trend of phylum Proteobacteria (Pseudomonadota), the proportion of phylum Bacteroidetes (Bacteroidota) has increased in active landfill compared to closed landfill, from the lowest 7% (C11) to the highest 23% (LT), members in class Cytophagia contributed main effort for this result, increased from less than 1% in samples C11 and C13 to around 15% in samples C16 and LT (Online Resource Table 5). Genera Pusillimonas and Sulfurimonas, which both are affiliated with Proteobacteria (Pseudomonadota), were predominant at samples in active and closed landfills, respectively; Pusillimonas representing 20% and 21% of the total genera for C11 and C13, Sulfurimonas accounting for 40% and 36% of the total genera for C16 and LT1 (Online Resource Table 8). The second dominant genus for C11 (8%) and C13 (7%) was Pseudomonas, while the genus Spirosoma, which belongs to phylum Bacteroidetes (Bacteroidota) was the second dominant genus for C16 (14%) and LT (16%) (Online Resource Table 8). Interestingly, Spirosoma were not present in the closed LFL samples, Pseudohongiella accounted for 6% and 4% of the genera at C11 and C13 but were not present in active landfill samples (Fig. 2C).
C11 and C13
Both C11 and C13 were sampled from closed aged of landfill cells, with C11 closed two years before C13. Subtle differences between the microbial composition of both samples were observed. Members of the phylum Proteobacteria (Pseudomonadota) occupied a lower proportion of the community present in sample C11 compared to C13, which is caused by the combined effect of an increase of Betaproteobacteria and Alphaproteobacteria and a decrease of Gammaproteobacteria. In contrast, the proportion of Balneolaeota (Balneolota) increased from 3% (C13) to 7% (C11), and the ratio of Spirochaetes (Spirochaetota) increased from 1% (C13) to 2% (C11), other phyla did not show obvious fluctuation between two samples (Online Resource Table 4). At genus level, Pseudomonas (Gammaproteobacteria), Candidatus Cyclonatronum (Balneolaeota (Balneolota)), Pseudohongiella (Gammaproteobacteria) and Leucobacter (Actinobacteria) displayed higher percentage in C11 compared with C13 (Fig. 2C).
LT and LT1
LT and LT1 were collected from the same leachate tank, however, LT1 was sampled 2.5 years before LT. The microbial composition of the two samples presented a significant difference. Proteobacteria (Pseudomonadota) experienced significant growth from 56% (LT) to 60% (LT1), mainly contributed by the increase of Epsilonproteobacteria and Betaproteobacteria. Similarly, Actinobacteria (Actinomycetota) increased from 3% (LT) to 6% (LT1) (Online Resource Tables 4 and 5). In contrast, Bacteroidetes (Bacteroidota) decreased from 23% (LT) to 21% (LT1), which was principally affected by the decrease of Cytophagia and Bacteroidia; with Firmicutes (Bacillota) decreasing from 8% (LT) to 5% (LT1) (Online Resource Tables 4 and 5). At the genus level, Sulfurimonas, which belongs to phylum Bacteroidetes (Bacteroidota), occupied the highest proportion of the bacterial communities present in both samples (36% (LT) and 39% (LT1)). Pusillimonas, from the phylum Proteobacteria (Pseudomonadota), remained constant at 5% (LT) and 6% (LT1). While Spirosoma, from the phylum Bacteroidetes (Bacteroidota), was the second most abundant genus in both samples, decreased with progressing landfill age, from 16% in LT to 14% in LT1 (Online Resource Table 8).