Skip to main content
Download PDF

Genome-centric metagenomic insights into the role of Chloroflexi in anammox, activated sludge and methanogenic reactors

BMC Microbiology volume 23, Article number: 45 (2023) Cite this article

Abstract

Background

The phylum Chloroflexi is highly abundant in a wide variety of wastewater treatment bioreactors. It has been suggested that they play relevant roles in these ecosystems, particularly in degrading carbon compounds and on structuring flocs or granules. Nevertheless, their function is not yet well understood as most species have not been isolated in axenic cultures. Here we used a metagenomic approach to investigate Chloroflexi diversity and their metabolic potential in three environmentally different bioreactors: a methanogenic full-scale reactor, a full-scale activated sludge reactor and a lab scale anammox reactor.

Results

Differential coverage binning approach was used to assemble the genomes of 17 new Chloroflexi species, two of which are proposed as new Candidatus genus. In addition, we recovered the first representative genome belonging to the genus ‘Ca. Villigracilis’. Even though samples analyzed were collected from bioreactors operating under different environmental conditions, the assembled genomes share several metabolic features: anaerobic metabolism, fermentative pathways and several genes coding for hydrolytic enzymes. Interestingly, genome analysis from the anammox reactor indicated a putative role of Chloroflexi in nitrogen conversion. Genes related to adhesiveness and exopolysaccharides production were also detected. Complementing sequencing analysis, filamentous morphology was detected by Fluorescent in situ hybridization.

Conclusion

Our results suggest that Chloroflexi participate in organic matter degradation, nitrogen removal and biofilm aggregation, playing different roles according to the environmental conditions.

Peer Review reports

Background

The phylum Chloroflexi comprises an ecologically and physiologically diverse group of bacteria, which have been detected in high abundance in methanogenic [1, 2], anammox [3,4,5] and activated sludge reactors [6, 7]. Despite of being highly abundant in wastewater treatment systems (WWTS), the basic ecophysiology of this group is still largely unknown because only few pure cultures have been investigated to date [8,9,10,11,12,13]. Based on the physiology of isolates, assembled genomes from metagenomes and in situ characterization (e.g., using microautoradiography combined with FISH), it has been proposed that their role is mainly related to the hydrolysis of complex organic matter, fermentation of carbohydrates and proteins, and degradation of debris from lysed bacterial cells [1, 7, 14,15,16,17,18].

In methanogenic reactors it has been hypothesized that the prevalence of this group is based on a syntrophic association with hydrogenotrophic methanogenic Archaea. This idea is supported by the fact that all Chloroflexi species isolated from methanogenic reactors required co-cultivation with methanogenic archaea for efficient growth [8, 10, 19, 20], the co-localization with filamentous Archaea [17], and the positive correlation between Chloroflexi and methanogens in anaerobic reactors [2].

In Anammox reactors, ammonium oxidation and nitrite reduction are coupled to form nitrogen gas under anoxic conditions [21]. In these reactors, the role of Chloroflexi has been associated to the scavenging of organic matter derived from anammox bacterial cell debris [3], or the utilization of soluble microbial products (SMP) and/or extracellular polymeric substances (EPS) produced by autotrophs [22]. This is supported by the high abundance of Chloroflexi in anammox reactors fed with synthetic wastewater containing NH4+ as the sole electron donor without addition of organic carbon compounds [3, 23,24,25]. Recent metagenomic-based works also suggested that members of Chloroflexi could facilitate a nitrite loop with anammox bacteria or support complete denitrification due to the expression of the nitric oxide reductase gene (norZ) and nitrite reductase genes (nirK, nirS). Thus, Chloroflexi would enhance the overall nitrogen removal performance in anammox bioreactors [26, 27]. In addition, some Chloroflexi members encode the function of biosynthesizing sticky macromolecular exopolysaccharide for anammox consortium aggregation [27].

On the other hand, it has been widely speculated that Chloroflexi members might play an important role during sludge granulation or flocculation processes. Evidence for this is their growth as filaments and the fact that some members showed cellular adhesiveness. This would enable Chloroflexi to form a backbone for small sludge particles and granules, in anammox, anaerobic and aerobic reactors [22, 28, 29]. Due to their filamentous growth, several authors have related the Chloroflexi overgrowth with bulking episodes, mainly in activated sludge reactors but also in full scale methanogenic reactors and lab-scale anammox reactors [30,31,32,33,34]. The study of Chloroflexi ecophysiology has direct implications in wastewater treatment plants performance, in particular to prevent and control bulking problems.

Here, we hypothesize that even though methanogenic, activated sludge and anammox reactors harbor different Chloroflexi species, their functions are akin. In order to test this hypothesis, we first performed 16S rRNA gene amplicon sequencing to determine the diversity of Chloroflexi in samples from anammox, activated sludge and methanogenic reactors. Then, we applied phylum-specific Fluorescence in situ hybridization (FISH) to determine the morphology of Chloroflexi in the different samples. Finally, we performed genome-centric metagenomics analysis to reveal their phylogenomic diversity and metabolic potential. The combination of different molecular approaches and the study of samples from different reactors allowed us to obtain information on the putative roles of Chloroflexi in different WWTS. Moreover, 17 new Chloroflexi species were identified and two of them are proposed as new Candidatus genus.

Methods

Sampling and DNA extraction

Samples from the following three reactors were analyzed (Table 1): 1) a full-scale upflow anaerobic sludge blanket (UASB) methanogenic reactor treating effluent from a malt industry (MO) ([35, 36], 2) a full-scale activated sludge reactor treating health-care waste treatment wastewater (Inffluent Chemical oxygen demand (CODinfluent) = 8,000 ppm, CODeffluent = 500 ppm, pH = 7.5, Hydraulic retention time (HRT) = 4.7 days), the wastes were sterilized by autoclave and the wastewater is generated in this process (RH), 3) a lab scale upflow anammox sludge blanket (UAnSB) reactor (S) fed with sewage after a partial nitrification system [37] (the inoculum used on this reactor was also included in the analysis). This inoculum was sludge from an anammox sequencing batch reactor (AM) fed with synthetic wastewater [38]. The full-scale reactors were located in Montevideo (Uruguay) and the lab-scale reactor was operated in the GENOCOV Research Group in Barcelona (Spain).

Table 1 Characteristics of the studied reactors
Full size table

All reactors were operated under mesophilic conditions. In order to apply the differential coverage binning approach to assemble genomes [39] several samples were taken from MO, RH and S (including AM used as inoculum of reactor S). Three samples taken at different time points were collected from the methanogenic (MO1, MO2 and MO3) and activated sludge reactors (RH4, RH5 and RH6). Another three samples were collected at different heights of the sludge bed from the UAnSB anammox reactor (S1, S3 and S5) and one from the sludge used as inoculum (AM).

Biomass samples were stored at -20 °C for sequencing workflows and fixed according to the protocol of Amann et al. [40] for FISH. Detailed procedures of DNA extraction and FISH are described in the Supplementary Material.

Community profiling using 16S rRNA gene amplicon sequencing

The taxonomic composition of the communities was studied by amplicon sequencing of the 16S rRNA gene with primers for V4 region [41, 42]. The raw data analysis was performed using Quantitative Insights Into Microbial Ecology’ pipeline (QIIME2 2020.11 release) [43]. The sequences were classified using MiDAS 3.7 database [44]. Detailed procedures of primers, QIIME2 analyses, data visualization, and phylogenetic analysis of the amplicon sequences are described in the Supplementary Material.

Metagenome sequencing

Metagenomes were shotgun sequenced by Illumina HiSeq 4000 platform (Macrogen, Seoul, Korea) using Library Kit TruSeq Nano DNA Kit (350 bp). The yield was approximately 5 Gb of raw short-read sequences per sample (2 × 100 bp).

Genome assembly, binning and genome annotation

The global quality of the metagenomes reads was checked using FastQC (v0.11.8) [45]. The reads were trimmed to remove adapters and bases below a quality score of 25 using Trimmomatic (v0.39) [46]. The trimmed reads from each reactor (MO and RH, separately) were pooled and assembled using MEGAHIT (v1.1.4–2) [47] with minimum k-mer length 43, maximum k-mer 75, with steps of four. The three samples from the anammox reactor S1, S3 and S5 were pooled and analyzed together with sample AM which was used as inoculum. The contigs obtained shorter than 1000 bp were removed. Quality filtered reads for each metagenome were mapped to the co-assembly contigs (> 1000 bp) using Bowtie2 (v2.3.4.1) [48] with default parameters. Genome bins were recovered using MetaBAT2 (v2.12.1) [49]. The completeness levels and contamination of the bins were assessed using CheckM (v1.0.13) [50]. The bins with an estimated completeness > 50% and contamination < 10% were reclassified using Genome Taxonomy Database GTDB-Tk (version v1.3.0 and GTDB-Tk reference data version r95) [51] to determine if they belonged to the Chloroflexi phylum. Chloroflexi bins were reassembled using SPAdes 3.10.0 [52]. All resulting contigs of > 1,500 bp were clustered using ESOM tools (emergent self-organizing map) [53], on the basis of its tetranucleotide frequency to identify and extract contaminant contigs in each genome. GUNC [54] was used to detect chimerism and to assign the taxonomy of the contigs in MAGs. Contigs assigned to a different phylum of Chloroflexi were removed from the MAGs. The CSS (clade separation score) estimated with gunc was closer to 0 for all genomes. A CSS closer to a value of 0 indicates that a genome is free of contamination and all genes are assigned to the same taxonomy.

Statistics of the reassembled bin were obtained using QUAST v5.0.2 and CheckM lineage wf (v1.0.13). Genes encoding ribosomal RNAs (rRNA) were predicted using the Bacterial rRNA predictor (Barrnap, https://github.com/tseemann/barrnap). The Average Nucleotide Identities (ANI) and the Average Amino Acid Identity (AAI) between MAGs and isolated species were determined using FastANI [55] and EzAAI [56]. The Chloroflexi genomes retrieved from methanogenic, activated sludge and anammox reactors were named as MO, RH and AMX, respectively, with a number suffix.

Genome phylogenetic analyses were performed with the GTDB-Tk tool. Newick format tree file was uploaded onto iTOL, a web-based tool for annotating and editing trees [57].

Metabolic analysis

Protein coding sequences (CDS) were determined with Prokka (1.14.5) [58, 59]. Predicted amino acid sequences were annotated with KOALA (KEGG Orthology And Links Annotation) for K number assignment of KEGG Genes [60] and Metaerg with KEGG, pfam, swiss-prot and tigrfams databases [61]. Pathways and genes completeness were represented as blocks with different colors in Figs. 4 and 5. Carbohydrate-active enzymes (CAZY) were determined using EnrichM v0.5.0 (https://github.com/geronimp/enrichM) with HMMs from dbCAN [62].

Results and discussion

Applying a series of molecular biology tools, the diversity and metabolic potential of the phylum Chloroflexi was studied in three different wastewater treatment systems: a methanogenic reactor (anaerobic, C-removal), an activated sludge reactor (aerobic, C-removal) and an anammox reactor (anaerobic, N-removal).

Diversity of Chloroflexi members in the different wastewater reactors according to 16S rRNA gene amplicon sequencing

The phylum Chloroflexi was detected in all reactor samples. Regarding the relative abundance higher values were detected in the anammox reactor samples compared to the methanogenic and aerobic reactors samples (26.9—33.5% for anammox, 4.9 -22% for activated sludge and 0.8—10.7% for methanogenic reactors) (Fig. 1).

Fig. 1
figure 1

Boxplots showing the relative abundance of the 16 most abundant phyla in the total community according to the 16S rRNA gene amplicon sequence in each reactor type

Full size image

MiDAS database, specifically designed to analyze microbial communities from activated sludge systems and methanogenic reactors [63], was used to classify sequences at different taxonomic levels (non-described microorganisms are labeled with a number according to the MiDAS taxonomy). The results confirmed that Anaerolineae class widely dominated the Chloroflexi community in all reactors (Fig. 2a), which was in accordance with results previously reported in aerobic, anaerobic and anammox reactors [2, 4, 17, 26, 27, 29, 36, 64,65,66,67,68]. Most of the genera detected were specifically associated with one type of reactor.

Fig. 2
figure 2

Heatmaps showing the relative abundance within the phylum Chloroflexi according to the 16S rRNA gene amplicon sequence: A) at class level, B) at genus level (30 most abundant genus) C) at species levels (30 most abundant species). Non-classified microorganisms are labeled with a number according to the MiDAS taxonomy. Black triangles in panels A and B mark genus and species belonging to Dehalococcoidia class, all the remaining microorganisms belonging to the class Anaerolineae. Red triangles in panel A mark ‘Ca. Villigracilis’

Full size image

One exception was ‘Ca. Villigracilis’, which was shared by all reactors with a similar relative abundance (Fig. 2b). This filamentous bacterium has been reported as widely distributed in activated sludge reactors, where it probably contributes to the matrix supporting floc formation [69]. Although ‘Ca. Villigracilis’ has been metabolically characterized by FISH-MAR and reported as the only facultative aerobic member described within Anaerolineae class so far, the genome is not yet available [69].

Several genera classified according to the MiDAS taxonomy, and therefore with no cultured representative, presented high abundance within the Anaerolineaceae and Caldilineaceae families in all reactors. Specifically, genera belonging to the family Caldilineaceae were the most abundant in the activated sludge reactor (Fig. 2b). Meanwhile, genera belonging to the family Anaerolineaceae were the most abundant in anammox and methanogenic reactors.

From all genera detected in our samples, only five harbor cultured representatives: Flexilinea, Bellilinea, Longilinea, Leptolinea and Anaerolinea. All of them have strains isolated from anaerobic reactors which were described as strictly anaerobic chemoorganotroph involved in carbohydrates fermentation [9, 10, 12]. In pure culture, their growth was enhanced in co-cultivation with a hydrogenotrophic methanogens. This could explain why sequences belonging to these genera were almost exclusively detected in the methanogenic reactor.

Moreover, all genera were represented by few dominant species which in most cases were not shared between reactors (Fig. 2c). In the samples studied, all Chloroflexi members presented filamentous morphology except for the anammox reactor samples (these results are described in the Supplementary material, section Morphology determined by FISH).

Several ASVs determined using the 16S rRNA gene sequences, could not be classified at class (e.g., ASV6009, ASV6008, ASV6010), genera (ASV6634) or species level (e.g., ASV6007, e.g., ASV6634, e.g., ASV6400, etc.). Even though MiDAS database represents a great advance in understanding the microbiomes of wastewater treatment systems, there are still undescribed microorganisms not included and more effort is needed to expand this database.

General genomic features of Chloroflexi metagenomic assembled genomes

We used differential coverage binning [39] to obtain metagenomic assembled genomes (MAGs) from the three different reactors. A total of 75, 48 and 56 MAGs with > 50 completeness and < 10% contamination were successfully retrieved from MO, AMX and RH reactors samples, respectively. Seventeen MAGs belonging to the phylum Chloroflexi (9 AMX, 4 MO and 4 RH) were analyzed further. Recovered genomes completeness ranged from 85.59% and 96.79%, and contamination after the reassembly and polishing (manually curated) ranged between 0% and 3.36% (Table 2).

Table 2 Statistics summary for the 17 Chloroflexi MAGs recovered in this study
Full size table

Three of these MAGs (AMX47, AMX55 and AMX56) were > 90% complete, had less than 5% contamination and importantly, included the full-length 16S, 23S, and 5S rRNA genes, and > 18 tRNA genes (Table 2), satisfying the criteria for high-quality (HQ) draft MAGs, according to the minimum information about a MAG (MIMAG) standard [70].

Phylogenomic analysis of Chloroflexi MAGs

To determine the phylogenetic position of the MAGs, a phylogenomic tree based on concatenated alignments of 120 single copy marker genes was constructed using 105 reference Chloroflexi genomes retrieved from NCBI (February 2021) and the 17 Chloroflexi MAGs obtained in our work. Of the 17 MAGs, six had sufficiently high degree of similarity (> 95% ANI to a representative genome [71]) to be classified at species level. The remaining MAGs presented ANI values below 95% with any reported genome and could not be classified at species level. According to the phylogenomic analysis, 7 were classified to genus level and 4 to family level (Table 2). The results showed that most of the MAGs (16 MAGs) were positioned into the Anaerolineae class while only one genome was positioned within the Dehalococcoidia class (Fig. 3a).

Fig. 3
figure 3

Phylogenomic and phylogenetic analysis of Chloroflexi based on the metagenomics results and the 16S rRNA gene sequence analysis: A) Phylogenomic tree constructed with the 17 genomes from this work and reference genomes retrieved from the NCBI, B) Phylogenetic tree constructed with the most abundant ASVs retrieved from the 16S rRNA gene amplicon sequences (with abundance higher than 5% in at least one sample) and the sequences of the 16S rRNA gene from 10 MAGs from this work, and reference sequences of the phylum Chloroflexi retrieved from NCBI. The tree was reconstructed using the ML method and the GTR model. ML bootstrap values greater than or equal to 70% are shown at each node. Bar indicates 0.01 substitutions per site. The ASV and MAGs are colored by reactor: blue for RH, orange for MO and green for AMX. Sequences from the phylum Thermotogota were used as outgroup for rooting trees

Full size image

At the order level, MAGs were distributed within the orders Anaerolineales, SBR1031, B4-G1, Caldilineales, UBA3071, UCB3 and UBA2991 (Table 2).

In addition, a phylogenetic tree was constructed based on the 16S rRNA gene sequences retrieved from 10 MAGs and the ASV sequences retrieved from the amplicon sequencing analysis (Fig. 3b, Fig. S2). Because of the few genomes available in the databases, it was not possible to construct both phylogenetic trees (based on 16S rRNA gene sequences and based on genomes sequences) with the exact same members of Chloroflexi. For example, ‘Ca. Villigracilis’ and ‘Ca. Sarcinatrix’, do not have representative genomes in databases so far, and were therefore not added to the genome-based phylogenetic tree. Moreover, although frequently used for phylogenetic identification, 16S rRNA gene sequences are notoriously difficult to assemble from metagenomes [72]. Thus, 16S rRNA genes were present only in 10 of the 17 MAGs (Table 2). Taking all this into consideration, the phylogenetic placement of 10 MAGs based on the 16S rRNA gene (Fig. 3b) was consistent with the phylogenomic tree based on 120 marker genes genes (Fig. 3a).

Interestingly, five MAGs from the anammox reactor (AMX55, AMX15, AMX39, AMX14 and AMX68), formed a monophyletic clade separated from other genomes in the phylogenomic tree (Fig. 3a). The 16S rRNA gene from AMX55 clustered with the sequence of ‘Ca. Villigracilis’ in the phylogenetic tree, with a 95.68% of identity according to the MiDAS taxonomy (Fig. 3b). Therefore, the retrieved MAG represents the first recovered genome within the ‘Ca. Villigracilis’ genus [73].

Within the ‘Ca. Villigracilis’ cluster, an ASV was classified within the UTCFX1 genus according to MiDAS. UTCFX1 represents a MAG retrieved from a nitritation-anammox sequencing batch reactor. The authors detected the presence nitrate reduction genes in this genome which were actively transcribed according to the metatranscriptomics analysis suggesting an important interaction with anammox microorganisms [26].

New genera and species candidates

In total, three MAGs (AMX47, AMX55 and AMX56) satisfied the criteria for high-quality recovered genomes [70] and met the minimal suggested required standards to be proposed as Candidates of new genera or species [74] (Table 2, Table S1). The full-length 16S rRNA gene sequence from AMX47, AMX55 and AMX56 showed identity values between 90.14% and 96.99% with their closest relatives according to Silva and MiDAS database comparisons (Table S2). The three MAGs presented ANI values < 79.2% with genomes from described species, and the AAI values for AMX47, AMX55 and AMX56 were 63.1, 57.6 and 56.6%, respectively (Table S1). Given the lack of close relatives, AMX47 and AMX56 genomes represent two different novel genera in the phylum, meanwhile, AMX55 represents the first representative genome of the genus ‘Ca. Villigracilis’ (based on the taxonomic sequence identity threshold recommendations of [75].

General analysis of metabolic pathways and genes

Different Chloroflexi genera and species predominate in the different wastewater treatment systems studied but, is their function redundant? For the 17 Chloroflexi MAGs, genes were annotated using a variety of protein databases to infer their metabolic potential. Because most of these MAGs were estimated to be between 85 and 93% complete, genes for additional pathways, might not be identified in this study. To overcome this problem, we also annotated genomes from described species closely related to our assembled genomes according to the phylogenomic tree (Fig. 3a). This allowed us to infer if a metabolic pathway of interest was not present in the genome or if its absence might be due to an incomplete genome recovery.

The metabolic pathway analysis was focused on answering the following questions:

  • 1-Are they potentially capable of performing aerobic and/or anaerobic respiration?

  • 2-Do they have the potential capability to hydrolyze different compounds?

  • 3-Which carbon compound degradation pathways do they have?

  • 4-Do they have N-removal potential?

  • 5-Do they have genes related to biomass adhesion properties?

Aerobic and anaerobic respiration pathways

Reactors were operated in anaerobic (anammox and methanogenic reactors) or aerobic (activated sludge reactor) conditions and therefore, genes related to aerobic and anaerobic respiration, or fermentative metabolisms were searched against the assembled MAG sequences.

Regarding respiration pathways, the genomic analysis showed that most MAGs (10 out of 17) had an incomplete oxidative phosphorylation pathway. Five MAGs (RH43, RH52, AMX55, AMX14 and AMX56) harbor the genes to use O2 and different nitrogenous compounds (N2O, NO3− and/or NO2−) as final electron acceptors, and two of them presented only the genes to use NO3− or NO2− (MO66 and AMX68, respectively) (Fig. 4, Fig. S3, Table S1).

Fig. 4
figure 4

Heatmap showing the completeness and the incompleteness of each metabolic pathway for the 17 MAG. Colors indicate pathways level completeness: dark brown (complete), light brown (1 block missing), beige (2 block missing), light green (more than 2 blocks missing) and white (not present)

Full size image

The complete phosphorylation pathway and the presence of nitrate and nitrite reductases in AMX55 was in accordance with the in situ characterization of ‘Ca. Villigracilis’, as this genus has the ability to take up substrates under anoxic conditions in presence and absence of nitrate/nitrite [69]. As we mentioned, the metatranscriptomics analysis showed that nitrate reduction genes from UTCFX1 (closely related to AMX55 in the phylogenetic tree, Fig. 3b) were expressed in an anammox reactor [26].

The tolerance to aerobic conditions is determined by the presence of oxidative stress protection genes. The presence of these genes was expected in genomes retrieved from the activated sludge reactor as these microorganisms are continuously exposed to oxygen (RH21 and RH38). Nevertheless, these genes were found in all MAGs (Fig. 4, Table S1). This is in concordance with the results obtained for ‘Ca. Brevefilum fermentans’ retrieved from an anaerobic digester [17].

Metabolic pathways involved in polymers hydrolysis, carbon respiration and fermentation

To test the hypothesis that Chloroflexi members are capable of recycling soluble microbial products acting as hydrolytic bacteria, we performed the annotation of genes for glycosyl hydrolases such as cellulases, endohemicellulases, amylases, amino sugar-degrading enzymes, oligosaccharide-degrading enzymes, and also extracellular peptidases. The presence of all these genes indicated that most MAGs have the potential of hydrolyzing cellulose, starch, protein and/or peptides (Fig. 4, Supplementary Data 1). This is in accordance with previous in situ studies, which revealed high level of surface associated hydrolytic enzymes and their involvement in the breakdown of complex organic compounds [14, 76]. This result indicates that the Chloroflexi phylum may play an active role in hydrolyzing complex organic matter in activated sludge, as well as in methanogenic and anammox reactors.

Regarding central carbon metabolism, all MAGs contained multiple transporters for different organic compounds (including sugars, amino acids, proteins and fatty acids) indicating that each species has alternative routes for incorporating and recycling dissolved organic matter (scavenge macronutrients) (Fig. 4, Supplementary Data 1). The Embden-Meyerhof-Parnas (EMP) pathway for glycolysis was complete in all MAGs except in RH21, MO66 and AMX47. RH21 and MO66 are phylogenetically related to Aggregatilinea lenta and Litorilinea aerophila, respectively, which have the complete glycolysis pathway. Thus, incomplete glycolysis pathways in RH21 and MO66 genomes may be due to incomplete genome recoveries. On the other hand, AMX47 was closely related to Dehalococcoides mccartyi, which does not have the complete glycolysis pathway suggesting that AMX47 does not perform glycolysis. The potential of AMX55 to consume glucose was consistent with the results showed for ‘Ca. Villigracilis’ [69]. Genes for pyruvate oxidation to acetyl-CoA were annotated in all MAGs (Fig. 4, Supplementary Data 1). Beta-oxidation was annotated in most of the MAGs and might represent an important metabolic route to obtain carbon and reducing equivalents for all Chloroflexi species. Another interesting finding was that genes for N-acetylglucosamine transportation and metabolism (PTS, nagAB) were annotated in RH43 (Fig. 4, Supplementary Data 1). These results are in accordance with other reports where Chloroflexi members appear to retrieve N-acetylglucosamine from lysed cells revealed by micro-autoradiography or FISH studies [3, 14, 23, 77,78,79,80,81,82]. In addition to the N-acetylglucosamine metabolism, we suggest that scavenging occurs through hydrolysis of complex organic compounds outside the cell, which are then transported to the cytoplasm, and are then metabolized via glycolysis or beta-oxidation (degradation of fatty acids and branched-chain amino acids) pathways. These results support the previous hypothesis that Chloroflexi has an important beneficial role in degradation of lysed bacterial cell debris and EPS.

Fermentation pathways including genes for acetate, ethanol, lactate, acetoin, formate and/or propionate production were present in all MAGs (Fig. 4, Supplementary Data 1). These results were in accordance with previous information of the Anaerolineae class (isolates and MAGS) retrieved from wastewater treatment system which were involved in sugar, amino acid or protein fermentation with different end products [8,9,10, 12].

Metabolic pathways for amino acids degradation and propionate formation (methylmalonyl-CoA pathway and aldehyde:ferredoxin oxidoreductase) was found in five MAGs (RH38, AMX14, AMX47, AMX55 and AMX68). These genes were also annotated in other Chloroflexi members [17, 83,84,85].

Most of the MAGs (13 MAGs) contained a formate dehydrogenase H (fdhF) to convert formic acid decomposition into CO2. Only AMX55 had the ability to convert formic acid into H2 and CO2 by the presence of both fdhF and hydrogenases, under anaerobic conditions in the absence of exogenous electron acceptors as was observed in previous studies [82, 86]. In addition, the formate transporter (focA) was annotated in most of these MAGs. Synergistic interaction of Anaerolineae members with the methanogenic archaea, was previously reported in anammox [87] and methanogenic reactors [2, 17, 88]. This could be a common scenario in these systems where the excess of electrons from organic carbon oxidation by Anaerolineae members could be transferred to methanogenic archaea.

Potential for polymers and lipids storage

The potential to store glycogen, polyhydroxyalkanoates and lipids was more common among MAGs retrieved from activated sludge than from methanogenic or anammox reactors. RH21, RH38, RH43, MO66 and AMX9 had complete glycogen biosynthesis and degradation pathways, suggesting that this polysaccharide may serve as a possible storage compound under unbalanced growth conditions (e.g., when C and/or N is temporally limited) (Fig. 4, Supplementary Data 1) [89, 90]. This result is in accordance with the information reported for ‘Ca. Amarolinea aalborgensis’ and ‘Ca. Promineofilum breve’ both retrieved from activated sludge [7, 16].

In addition to glycogen as storage compound, RH43 and RH52, encodes a potential biosynthesis pathway for polyhydroxyalkanoates (PHA) (Fig. 4, Supplementary Data 1), which are usually formed as carbon and energy storage compounds [91] under conditions of carbon excess and nitrogen or phosphate limitation [92]. Thus, the potential for PHA storage could favor these organisms with intermittent carbon availability present in these systems. Also, the annotation of a putative acyl-CoA:DAG acyltransferases (atfA), which catalyzes the final step in the synthesis of triacylglycerols, indicates the potential for lipid storage [93] in RH21, RH38, AMX56 and AMX57. This capability was also reported for ‘Ca. Promineofilum breve’ [16].

Presence of Nitrogen metabolism pathways

Different genes related to the nitrogen cycle such as the dissimilatory nitrate reduction to ammonia (DNRA) and partial denitrification were mostly annotated in Chloroflexi MAGs from anammox reactors. DNRA is part of the nitrogen cycle, and it has been annotated in several Chloroflexi genomes [7, 16, 26]. In this work, three Chloroflexi species from anammox reactors (AMX14, AMX55 and AMX56) have the potential capability to perform DNRA due to the presence of both dissimilatory nitrate (narGHI) and nitrite reductases (nrfAH) (Fig. 5).

Fig. 5
figure 5

Heatmap showing the nitrogen cycle annotated genes in the 17 MAGs from activated sludge, methanogenic and anammox reactors, Nitrogen cycle, next to each gene, the number of genomes that contain that gene is shown. Blue: activated sludge. Orange: methanogenic. Green: anammox. Dark brown blocks indicate genes that are present and white blocks indicate genes that are not present

Full size image

AMX14 and AMX55 genomes were positioned in the ‘Ca. Villigracilis’ and UTCFX1 cluster. As mentioned above, UTCFFX1 had nitrate reductase genes activity [26]. Thus, probably both species might perform DNRA. The six remaining MAGs have the potential to reduce nitrate to nitrite (narGHI), or nitrite to ammonium (nrfAH). Our results show that all MAGs obtained from the anammox reactor could have a specific role in the nitrogen cycle, supplying nitrite or ammonia to anammox bacteria enhancing the overall nitrogen removal performance.

Two MAGs retrieved from the activated sludge reactor harbor genes to perform the reduction of nitrite to ammonia. These results are in accordance to those reported for ‘Ca. Amarolinea’ (Chloroflexi MAG retrieved from activated sludge) [7]. Despite its potential, the in situ characterization did not confirm the use of nitrate or nitrite [18]. Also, putative nitrite reductase (nrfAH) has been found in Caldilinea aerophila and Anaerolinea thermophile genomes, but their ability to utilize nitrite as an electron acceptor to support anaerobic growth is yet to be assessed [8].

A nitrite reductase and nitric oxide reductase were annotated in RH21, suggesting the potential to denitrify (Fig. 5). This ability has only been proved for one Chloroflexi species, Ardenticatenia maritima [94]. In the anammox reactor, the presence of nirK/nirS and nosZ genes was common among the assembled MAGs. It was consistent with several studies which report that denitrification is more frequent in Chloroflexi species from anammox reactors than in the ones retrieved from aerobic and anaerobic reactors [26, 27, 68]. A recent metagenomic study revealed that most of the heterotrophic organisms in anammox granules encode the ability to respire nitrate via partial denitrification, possibly completing a nitrite loop with anammox and nitrite oxidizing bacteria (NOB) by reducing nitrate back to nitrite [68]. This activity could contribute to the removal of excess nitrate produced from the system during anammox growth or nitrite oxidation by NOB.

Filamentous morphology, adhesion capability and exopolysaccharides production

It has been extensively proposed that Chloroflexi plays an important role in granule and floc formation. In order to deepen our knowledge on this subject, three important capabilities were specifically searched in the MAGs: the filamentous morphology (studied by FISH), the adhesiveness and the production of exopolysaccharides. A complete set of genes for the pilus assembly (pilA, CpaB, CpaE, CpaF, TadB, TadC) which favor adhesiveness, was annotated in MO53, AMX9, AMX47 and AMX56 (Fig. 4, Supplementary Data 1). The other MAGs mainly missed only the pilus assembly protein CpaB. As has been already noted, pili are often involved in facilitating adhesion and colonization in a wide variety of scenarios. In a previous study of Anaerolineae members, the adhesiveness by the expression of the tight adherence protein (Tad) on the active type VI pili indicated its function for cellular attachment, which was further testified to be more likely related to cell aggregation other than cellulose surface adhesion [29].

Regarding the production of exopolysaccharides, it has been suggested that some Chloroflexi-affiliated bacteria encode the function of biosynthesizing sticky macromolecular exopolysaccharides such as UDP-GlcNAcA, GDP-Man, and GDP-Rha from partial nucleotide sugars biosynthesized by anammox bacteria (UDP-ManNAc and CDP-Glc) [27]. We found the same set of genes for sticky macromolecular exopolysaccharide synthesis in all MAGs from RH, 3 MAGs from AMX and 1 MAG from MO, including the enzyme lactate dehydrogenase (prerequisite for exopolysaccharide production) (Fig. 4, Supplementary Data 1). The rest of the MAGs lack some of the necessary genes. Therefore, the results obtained showed that some Chloroflexi have the capability to promote cell aggregation and consequently the formation of cores or carriers, which help to form the initial framework of small sludge particles. Zhao et al. [87] reported that Anaerolineae members affected the nitrogen removal performance through affecting the aggregation because of the positive correlation relationship of this group with the abundance of EPS formation genes. The results obtained showed that key genes involved in B-vitamin biosynthesis were missing in all Chloroflexi MAGs. For example, genes for thiamin (vitamin B1) biosynthesis (thiamine-phosphate synthase and thiamine-monophosphate kinase), biotin (vitamin B7) biosynthesis (adenosylmethionine-8-amino-7-oxononanoate aminotransferase and biotin synthase) and adenosylcobalamin (vitamin B12) biosynthesis (cobalamin synthase and adenosylcobinamide-phosphate synthase) were absent in all Chloroflexi MAGs (Supplementary Data 1) as previously reported by Lawson et al. 2017 [95]. These authors postulated that in anammox reactors Brocadia sp. could supply B-vitamins to Chloroflexi members. This suggests that other microorganisms may support B-vitamin requirements for Chloroflexi community. This hypothesis could be extended to activated sludge and methanogenic reactors.

The role of Chloroflexi species in each reactor

The role of the different MAGs was postulated according to the genome annotation and the reactor´s characteristics (Table 3).

Table 3 Role proposed for the Chloroflexi microorganisms in the different reactors according to the metabolic function detected in the MAGs
Full size table

Although the reactors were fed with very different wastewaters and operated under different conditions (anaerobic, aerobic and anoxic) a common role can be postulated in the degradation of polymers (polysaccharides or proteins) either from the wastewater or from biomass decay or exopolymers secreted by cells. Another important role can be also postulated in the maintenance of cell aggregates as granules or floccules. Yet, Chloroflexi species can grow by using a fermentative way of life (anaerobic), aerobic respiration or nitrate-nitrite reduction, depending on the reactor’s operation.

In conclusion, although the information retrieved from the metagenomics analysis showed shared functions in the different reactors, adaptation to each operation condition was observed, indicating high versatility within this group of microorganisms.

General Chloroflexi features and metabolism

A summary of the metabolic and structural functions from the MAGs retrieved from the three ecosystems is shown in Fig. 6 and Table 3.

Fig. 6
figure 6

Metabolic model of the Chloroflexi MAGs by reactor: Blue: MAGs retrieved from the aerobic reactor, Orange: MAGs retrieved from the methanogenic reactor, Green: MAGs retrieved from the anammox reactor and its inoculum. The identifying number of each MAG appears in the boxes when the gene or metabolic pathway is present

Full size image

The genomic annotation of the 17 MAGs revealed the potential for facultative and strict metabolism. All MAGs had the potential of a versatile metabolism related to the hydrolysis and fermentation of complex and simple organic compounds, suggesting that these hydrolytic and fermentative bacteria occupy a niche in recycling microbial dissolved organic matter. In addition to the filamentous morphology, some species could have a crucial structural role providing the structural matrix around which floc and granules material aggregates due to their adhesiveness and EPS production. Also, these characteristics could represent a selective advantage for the retention of these bacteria in the reactors.

Although the MAGs were retrieved from very different ecosystems, there was functional redundancy in the same reactor and analogous functions between reactors related to carbohydrates metabolism, and a more specific role related to the nitrogen removal in anammox reactors (Fig. 6). The functional redundancy of Chloroflexi members enables functional resistance to reactor operation disturbances due to the presence of multiple species that can perform the same metabolic function.

While MAGs analysis expanded our knowledge on the diversity and potential function of the phylum Chloroflexi, further experiments are necessary to confirm the expression of the identified metabolic functions under certain operational conditions. On the other hand, studying gene expression in bulking episodes would deepen the knowledge about their implication in these events. This knowledge would allow us to understand further which Chloroflexi groups are functionally important in these systems and how their disappearance or overgrowth could affect the processes. And as a final goal manage them and for example reduce and control bulking events. The difficulty of isolating organisms of the phylum Chloroflexi are probably related to their slow growth. New tools have been developed that facilitate the isolation of slow-growing microorganisms such as Droplet-based high-throughput cultivation [96]. Therefore, it is expected in the coming years to have a greater quantity of Chloroflexi members isolated in pure culture, in order to reveal and confirm their role in wastewater treatment systems.

Taxonomic proposal for the new Candidates

‘Ca. Villigracilis nielsenii’ sp. nov.

As a name for the new candidate species (represented by the AMX55 genome) within the genus ‘Ca. Villigracilis’, we suggest ‘Ca. Villigracilis nielsenii’. Niel.se’ni.i. N.L. gen. masc. n. nielsenii, of Nielsen, named in honor of Per Halkjær Nielsen, the Danish environmental microbiologist who made important contributions to research and practices in the field of Microbial Ecology and Water Engineering.

Conclusion

This study provides a first insight into the diversity and metabolic potential of 17 Chloroflexi species in three different wastewater treatment bioreactors, unveiling the ecological role of many species with no representatives in pure culture. The recovered genomes appear in several taxonomic orders, attesting to the broad diversity of Chloroflexi members in WWTS.

Our results suggest that Chloroflexi participate in organic matter degradation, nitrogen removal and biofilm aggregation, playing different roles according to the environmental conditions. Additionally, we proposed two novel genera within the classes Dehalococcoidia and Anaerolineae, and a new species of the genus ‘Ca. Villigracilis’ representing the first representative genome within this genus.

Availability of data and materials

The raw amplicon and metagenome sequences and the 17 MAGs were deposited under the NCBI BioProjects: PRJNA728853 and PRJNA780299.

Scripts are available in GitHub (https://github.com/PatoUru/Bovio-Winkler_etal_2022/blob/main/Bovio-Winkler_etal_2022).

Abbreviations

WWTS:

Wastewater treatment systems

SMP:

Soluble microbial products

EPS:

Extracellular polymeric substances

FISH:

Fluorescence in situ hybridization

UASB:

Upflow anaerobic sludge blanket

UAnSB:

Upflow anammox sludge blanket

MAGs:

Metagenomic assembled genomes

ANI:

Average Nucleotide Identities

AAI:

Average Amino Acid Identity

MIMAG:

Minimum information about a MAG

ASV:

Amplicon sequence variant

References

  1. Petriglieri F, Nierychlo M, Nielsen PH, McIlroy SJ. In situ visualisation of the abundant Chloroflexi populations in full-scale anaerobic digesters and the fate of immigrating species. PLoS ONE. 2018;13(11):1–14.

    Article  Google Scholar 

  2. Bovio-Winkler P, Cabezas A, Etchebehere C, Chandran K H, Lee PK, Halkjaer Nielsen P, et al. Database Mining to Unravel the Ecology of the Phylum Chloroflexi in Methanogenic Full Scale Bioreactors. Front Microbiol. 2021;11(January):603234.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Kindaichi T, Yuri S, Ozaki N, Ohashi A. Ecophysiological role and function of uncultured Chloroflexi in an anammox reactor. Water Sci Technol. 2012;66(12):2556–61.

    Article  CAS  PubMed  Google Scholar 

  4. Cao S, Du R, Li B, Ren N, Peng Y. High-throughput profiling of microbial community structures in an ANAMMOX-UASB reactor treating high-strength wastewater. Appl Microbiol Biotechnol. 2016;100(14):6457–67.

    Article  CAS  PubMed  Google Scholar 

  5. Wang S, Guo J, Lian J, Ngo HH, Guo W, Liu Y, et al. Rapid start-up of the anammox process by denitrifying granular sludge and the mechanism of the anammox electron transport chain. Biochem Eng J. 2016;15(115):101–7.

    Article  Google Scholar 

  6. Speirs LBM, Dyson ZA, Tucci J, Seviour RJ. Eikelboom filamentous morphotypes 0675 and 0041 embrace members of the Chloroflexi: resolving their phylogeny, and design of fluorescence in situ hybridisation probes for their identification. FEMS Microbiol Ecol. 2017;93(10):1–13.

    Article  Google Scholar 

  7. Andersen MH, McIlroy SJ, Nierychlo M, Nielsen PH, Albertsen M. Genomic insights into Candidatus Amarolinea aalborgensis gen. nov., sp. nov., associated with settleability problems in wastewater treatment plants. Syst Appl Microbiol. 2019;42(1):1–8.

    Article  Google Scholar 

  8. Sekiguchi Y, Yamada T, Hanada S, Ohashi A, Harada H, Kamagata Y. Anaerolinea thermophila gen. nov., sp. nov. and Caldilinea aerophila gen. nov., sp. nov., novel filamentous thermophiles that represent a previously uncultured lineage of the domain bacteria at the subphylum level. Int J Syst Evol Microbiol. 2003;53(6):1843–51.

    Article  CAS  PubMed  Google Scholar 

  9. Yamada T, Imachi H, Ohashi A, Harada H, Hanada S, Kamagata Y, et al. Bellilinea caldifistulae gen. nov., sp. nov and Longilinea arvoryzae gen. nov., sp. nov., strictly anaerobic, filamentous bacteria of the phylum Chloroflexi isolated from methanogenic propionate-degrading consortia. Int J Syst Evol Microbiol. 2007;57(10):2299–306.

    Article  CAS  PubMed  Google Scholar 

  10. Yamada T, Sekiguchi Y, Hanada S, Imachi H, Ohashi A, Harada H, et al. Anaerolinea thermolimosa sp. nov., Levilinea saccharolytica gen. nov., sp. nov. and Leptolinea tardivitalis gen. nov., sp. nov., novel filamentous anaerobes, and description of the new classes Anaerolineae classis nov. and Caldilineae classis nov. in the. Int J Syst Evol Microbiol. 2006;56(6):1331–40.

    Article  CAS  PubMed  Google Scholar 

  11. Yoon DN, Park SJ, Kim SJ, Jeon CO, Chae JC, Rhee SK. Isolation, characterization, and abundance of filamentous members of Caldilineae in activated sludge. J Microbiol. 2010;48(3):275–83.

    Article  CAS  PubMed  Google Scholar 

  12. Sun L, Toyonaga M, Ohashi A, Matsuura N, Tourlousse DM, Meng XY, et al. Isolation and characterization of Flexilinea flocculi gen. nov., sp. nov., a filamentous, anaerobic bacterium belonging to the class Anaerolineae in the phylum Chloroflexi. Int J Syst Evol Microbiol. 2016;66:988–96.

    Article  CAS  PubMed  Google Scholar 

  13. Yamada T, Sekiguchi Y. Cultivation of uncultured chloroflexi subphyla: significance and ecophysiology of formerly uncultured chloroflexi “subphylum i” with natural and biotechnological relevance. Microbes and environments / JSME. 2009;24(3):205–16.

    Article  Google Scholar 

  14. Kragelund C, Levantesi C, Borger A, Thelen K, Eikelboom D, Tandoi V, et al. Identity, abundance and ecophysiology of filamentous Chloroflexi species present in activated sludge treatment plants. FEMS Microbiol Ecol. 2007;59(3):671–82.

    Article  CAS  PubMed  Google Scholar 

  15. Nielsen PH, Kragelund C, Seviour RJ, Nielsen JL. Identity and ecophysiology of filamentous bacteria in activated sludge. FEMS Microbiol Rev. 2009;339:69–98.

    Google Scholar 

  16. McIlroy SJ, Karst SM, Nierychlo M, Dueholm MS, Albertsen M, Kirkegaard RH, et al. Genomic and in situ investigations of the novel uncultured Chloroflexi associated with 0092 morphotype filamentous bulking in activated sludge. ISME J. 2016;10(9):2223–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. McIlroy SJ, Kirkegaard RH, Dueholm MS, Fernando E, Karst SM, Albertsen M, et al. Culture-independent analyses reveal novel anaerolineaceae as abundant primary fermenters in anaerobic digesters treating waste activated sludge. Frontiers in microbiology. 2017;8:1134.

  18. Nierychlo M, Miłobȩdzka A, Petriglieri F, McIlroy B, Nielsen PH, McIlroy SJ. The morphology and metabolic potential of the Chloroflexi in full-scale activated sludge wastewater treatment plants. FEMS Microbiol Ecol. 2019;95(2):1–11.

    Article  Google Scholar 

  19. Yamada T, Sekiguchi Y, Imachi H, Kamagata Y, Ohashi A, Harada H. Diversity, localization, and physiological properties of filamentous microbes belonging to Chloroflexi subphylum I in mesophilic and thermophilic methanogenic sludge granules. Appl Environ Microbiol. 2005;71(11):7493–503.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Yamada T, Yamauchi T, Shiraishi K, Hugenholtz P, Ohashi A, Harada H, et al. Characterization of filamentous bacteria, belonging to candidate phylum KSB3, that are associated with bulking in methanogenic granular sludges. ISME J. 2007;1:246–55.

    Article  CAS  PubMed  Google Scholar 

  21. Strous M, Heijnen JJ, Kuenen JG, Jetten MSM. The sequencing batch reactor as a powerful tool for the study of slowly growing anaerobic ammonium-oxidizing microorganisms. Appl Microbiol Biotechnol. 1998;50(5):589–96.

    Article  CAS  Google Scholar 

  22. Chu ZR, Wang K, Li XK, Zhu MT, Yang L, Zhang J. Microbial characterization of aggregates within a one-stage nitritation-anammox system using high-throughput amplicon sequencing. Chemica En J. 2015;262:41–8.

    Article  CAS  Google Scholar 

  23. Kindaichi T, Ito T, Okabe S. Ecophysiological Interaction between Nitrifying Bacteria and Heterotrophic Bacteria in Autotrophic Nitrifying Biofilms as Determined by Microautoradiography-Fluorescence In Situ Hybridization. Appl Environ Microbiol. 2004;70(3):1641–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Strous M, Pelletier E, Mangenot S, Rattei T, Lehner A, Taylor MW, et al. Deciphering the evolution and metabolism of an anammox bacterium from a community genome. Nature. 2006;440(7085):790–4.

    Article  PubMed  Google Scholar 

  25. Cho S, Fujii N, Lee T, Okabe S. Development of a simultaneous partial nitrification and anaerobic ammonia oxidation process in a single reactor. Bioresour Technol. 2011;102(2):652–9.

    Article  CAS  PubMed  Google Scholar 

  26. Lawson CE, Wu S, Bhattacharjee AS, Hamilton JJ, McMahon KD, Goel R, et al. Metabolic network analysis reveals microbial community interactions in anammox granules. Nat Commun. 2017;8(May):1–12.

    Google Scholar 

  27. Zhao Y, Liu SS, Jiang B, Feng Y, Zhu T, Tao H, et al. Genome-Centered Metagenomics Analysis Reveals the Symbiotic Organisms Possessing Ability to Cross-Feed with Anammox Bacteria in Anammox Consortia. Environ Sci Technol. 2018;52(19):11285–96.

    Article  CAS  PubMed  Google Scholar 

  28. Gao D, Liu L, Liang H, Wu WM. Aerobic granular sludge: Characterization, mechanism of granulation and application to wastewater treatment. Critical Rev Biotechnol. 2011;31:137–52.

    Article  CAS  Google Scholar 

  29. Xia Y, Wang Y, Wang Y, Chin FYL, Zhang T. Cellular adhesiveness and cellulolytic capacity in Anaerolineae revealed by omics-based genome interpretation. Biotechnol Biofuels. 2016;9(1):1–13.

    Article  CAS  Google Scholar 

  30. Sekiguchi Y, Takahashi H, Kamagata Y, Ohashi A, Harada H. In Situ Detection, Isolation, and Physiological Properties of a Thin Filamentous Microorganism Abundant in Methanogenic Granular Sludges: A Novel Isolate Affiliated with a Clone Cluster, the Green Non-Sulfur Bacteria. Subdivision I Appl Environ Microbiol. 2001;67(12):5740–9.

    Article  CAS  PubMed  Google Scholar 

  31. Björnsson L, Hugenholtz P, Tyson GW, Blackall LL. Filamentous Chloroflexi (green non-sulfur bacteria) are abundant in wastewater treatment processes with biological nutrient removal. Microbiology (N Y). 2002;148(8):2309–18.

    Google Scholar 

  32. Li J, Hu B, Zheng P, Qaisar M, Mei L. Filamentous granular sludge bulking in a laboratory scale UASB reactor. Bioresour Technol. 2008;99(9):3431–8.

    Article  CAS  PubMed  Google Scholar 

  33. Song YX, Liao Q, Yu C, Xiao R, Tang CJ, Chai LY, et al. Physicochemical and microbial properties of settled and floating anammox granules in upflow reactor. Biochem Eng J. 2017;123:75–85.

    Article  CAS  Google Scholar 

  34. Nierychlo M, McIlroy SJ, Kucheryavskiy S, Jiang C, Ziegler AS, Kondrotaite Z, et al. Candidatus Amarolinea and Candidatus Microthrix Are Mainly Responsible for Filamentous Bulking in Danish Municipal Wastewater Treatment Plants. Front Microbiol. 2020;11(June):1–17.

    Google Scholar 

  35. Borzacconi L, López I, Passeggi M, Etchebehere C, Barcia R. Sludge deterioration in a full scale UASB reactor after a pH drop working under low loading conditions. Water Sci Technol. 2008;57(5):797–802.

    Article  CAS  PubMed  Google Scholar 

  36. Bovio P, Cabezas A, Etchebehere C. Preliminary analysis of Chloroflexi populations in full-scale UASB methanogenic reactors. J Appl Microbiol. 2019;126(2):667–83.

    Article  CAS  PubMed  Google Scholar 

  37. Juan-Díaz X, Pérez J, Carrera J. Effective dampening of temperature effects in an anammox reactor treating real mainstream wastewater. J Water Proc Eng. 2021;1(40):101853.

    Article  Google Scholar 

  38. Isanta E, Bezerra T, Fernández I, Suárez-Ojeda ME, Pérez J, Carrera J. Microbial community shifts on an anammox reactor after a temperature shock using 454-pyrosequencing analysis. Bioresour Technol. 2015;1(181):207–13.

    Article  Google Scholar 

  39. Albertsen M, Hugenholtz P, Skarshewski A, Nielsen KL, Tyson GW, Nielsen PH. Genome sequences of rare, uncultured bacteria obtained by differential coverage binning of multiple metagenomes. Nat Biotechnol. 2013;31(6):533–8.

    Article  CAS  PubMed  Google Scholar 

  40. Amann RI, Ludwig W, Schleifer KH. Phylogenetic identification and in situ detection of individual microbial cells without cultivation. Microbiol Rev. 1995;59(1):143–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Claesson MJ, O’Sullivan O, Wang Q, Nikkilä J, Marchesi JR, Smidt H, et al. Comparative analysis of pyrosequencing and a phylogenetic microarray for exploring microbial community structures in the human distal intestine. PLoS One. 2009;4(8):e6669.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Callejas C, Fernández A, Passeggi M, Wenzel J, Bovio P, Borzacconi L, et al. Microbiota adaptation after an alkaline pH perturbation in a full-scale UASB anaerobic reactor treating dairy wastewater. Bioprocess Biosyst Eng. 2019;42(12):2035–46.

    Article  CAS  PubMed  Google Scholar 

  43. Bolyen E, Rideout JR, Caporaso JG. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):850–2.

    Article  Google Scholar 

  44. Nierychlo M, Andersen KS, Xu Y, Green N, Jiang C, Albertsen M, et al. MiDAS 3: An ecosystem-specific reference database, taxonomy and knowledge platform for activated sludge and anaerobic digesters reveals species-level microbiome composition of activated sludge. Water Res. 2020;1(182):115955. Available from: https://doi.org/10.1016/j.watres.2020.115955.

    Article  CAS  Google Scholar 

  45. Andrews. FastQC: a quality control tool for high throughput sequence data. Babraham Institute. 2010. Available from: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/

  46. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20 Available from: http://www.ncbi.nlm.nih.gov/pubmed/24695404. Cited 16 May 2019.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Li D, Luo R, Liu CM, Leung CM, Ting HF, Sadakane K, et al. MEGAHIT v10: A fast and scalable metagenome assembler driven by advanced methodologies and community practices. Methods. 2016;102:3–11.

    Article  CAS  PubMed  Google Scholar 

  48. Kim D, Langmead B, Salzberg SL. HISAT: A fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Kang DD, Froula J, Egan R, Wang Z. MetaBAT, an efficient tool for accurately reconstructing single genomes from complex microbial communities. PeerJ. 2015;2015(8):1–15.

    Google Scholar 

  50. Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: Assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25(7):1043–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Parks DH, Chuvochina M, Waite DW, Rinke C, Skarshewski A, Chaumeil PA, et al. A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life. Nat Biotechnol. 2018;36(10):996.

    Article  CAS  PubMed  Google Scholar 

  52. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Ultsch A, Fabian M, Mörchen F. ESOM-Maps: tools for clustering, visualization, and classification with Emergent SOM. Univ., 2005. Vol. 46.

  54. Orakov A, Fullam A, Coelho LP, Khedkar S, Szklarczyk D, Mende DR, et al. GUNC: detection of chimerism and contamination in prokaryotic genomes. Genome Biol. 2021;22(1):1–19 Available from: https://genomebiology.biomedcentral.com/articles/10.1186/s13059-021-02393-0. Cited 21 Nov 2022.

    Article  Google Scholar 

  55. Jain C, Rodriguez-R LM, Phillippy AM, Konstantinidis KT, Aluru S. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nature Commun. 2018;9(1):1–8.

    Article  Google Scholar 

  56. Kim D, Park S, Chun J. Introducing EzAAI: a pipeline for high throughput calculations of prokaryotic average amino acid identity. J Microbiol. 2021;59(5):476–80.

    Article  PubMed  Google Scholar 

  57. Letunic I, Bork P. Interactive Tree of Life (iTOL) v4: Recent updates and new developments. Nucleic Acids Res. 2019;47(W1):256–9.

    Article  Google Scholar 

  58. Hyatt D, Chen GL, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: Prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;8:11.

    Google Scholar 

  59. Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinforma. 2014;30(14):2068-9.

  60. Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of Genes and Genomes [Internet]. Nucleic Acids Research. 2000;28:27–30 Oxford University Press. Available from: https://pmc/articles/PMC102409/?report=abstract/. Cited 17 Jan 2021.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Dong X, Strous M. An Integrated Pipeline for Annotation and Visualization of Metagenomic Contigs. Front Genet. 2019;15:10.

    Google Scholar 

  62. Cantarel BL, Coutinho PM, Rancurel C, Bernard T, Lombard V, Henrissat B. The Carbohydrate-Active EnZymes database (CAZy): an expert resource for glycogenomics. Nucleic Acids Res. 2009;37(suppl_1):D233-8.

  63. Nierychlo M, Andersen KS, Xu Y, Green N, Jiang C, Albertsen M, Dueholm MS, Nielsen PH. MiDAS 3: an ecosystem-specific reference database, taxonomy and knowledge platform for activated sludge and anaerobic digesters reveals species-level microbiome composition of activated sludge. Water Res. 2020;182:115955.

  64. Mielczarek AT, Kragelund C, Eriksen PS, Nielsen PH. Population dynamics of filamentous bacteria in Danish wastewater treatment plants with nutrient removal. Water Res. 2012;46(12):3781–95.

    Article  CAS  PubMed  Google Scholar 

  65. Gonzalez-Gil G, Sougrat R, Behzad AR, Lens PNL, Saikaly PE. Microbial Community Composition and Ultrastructure of Granules from a Full-Scale Anammox Reactor. Microb Ecol. 2015;70(1):118–31.

    Article  CAS  PubMed  Google Scholar 

  66. Pereira AD, Cabezas A, Etchebehere C, de Chernicharo CAL, de Araújo JC. Microbial communities in anammox reactors: a review. Environ Technol Rev. 2017;6(1):74–93.

    Article  CAS  Google Scholar 

  67. Campanaro S, Treu L, Rodriguez-R LM, Kovalovszki A, Ziels RM, Maus I, et al. New insights from the biogas microbiome by comprehensive genome-resolved metagenomics of nearly 1600 species originating from multiple anaerobic digesters. Biotechnol Biofuels. 2020;13(1):1–18.

    Article  Google Scholar 

  68. Speth DR, in’t Zandt MH, Guerrero-Cruz S, Dutilh BE, Jetten MS. Genome-based microbial ecology of anammox granules in a full-scale wastewater treatment system. Nat Commun. 2016;7(1):11172.

  69. Nierychlo M, Miłobȩdzk A, Petriglieri F, McIlroy B, Nielsen PH, McIlroy SJ. The morphology and metabolic potential of the Chloroflexi in full-scale activated sludge wastewater treatment plants. FEMS Microbiol Ecol. 2019;95(2):1–11.

    Article  Google Scholar 

  70. Bowers RM, Kyrpides NC, Stepanauskas R, Harmon-Smith M, Doud D, Reddy TBK, et al. Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea [Internet]. Nature Biotechnology. 2017;35:725–31 Nature Publishing Group.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Parks DH, Chuvochina M, Chaumeil PA, Rinke C, Mussig AJ, Hugenholtz P. A complete domain-to-species taxonomy for Bacteria and Archaea. Nature Biotechnol. 2020;38(9):1079–86.

    Article  CAS  Google Scholar 

  72. Blazejak A, Schippers A. High abundance of JS-1- and Chloroflexi-related Bacteria in deeply buried marine sediments revealed by quantitative, real-time PCR. FEMS Microbiol Ecol. 2010;72(2):198–207.

    Article  CAS  PubMed  Google Scholar 

  73. Rosselló-Móra R, Amann R. Past and future species definitions for Bacteria and Archaea. Syst Appl Microbiol. 2015;38(4):209–16.

    Article  PubMed  Google Scholar 

  74. Konstantinidis KT, Rosselló-Móra R, Amann R. Uncultivated microbes in need of their own taxonomy. ISME Journal. 2017;11:2399–406.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Yarza P, Yilmaz P, Pruesse E, Glöckner FO, Ludwig W, Schleifer KH, et al. Uniting the classification of cultured and uncultured bacteria and archaea using 16S rRNA gene sequences. Nat Rev Microbiol. 2014;12(9):635–45.

    Article  CAS  PubMed  Google Scholar 

  76. Xia Y, Kong Y, Nielsen PH. In situ detection of protein-hydrolysing microorganisms in activated sludge. FEMS Microbiol Ecol. 2007;60(1):156–65.

    Article  CAS  PubMed  Google Scholar 

  77. Okabe S, Kindaichi T, Ito T. Fate of 14C-labeled microbial products derived from nitrifying bacteria in autotrophic nitrifying biofilms. Appl Environ Microbiol. 2005;71(7):3987–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Zang K, Kurisu F, Kasuga I, Furumai H, Yagi O. Analysis of the phylogenetic diversity of estrone-degrading bacteria in activated sewage sludge using microautoradiography-fluorescence in situ hybridization. Syst Appl Microbiol. 2008;31(3):206–14.

    Article  CAS  PubMed  Google Scholar 

  79. Miura Y, Okabe S. Quantification of cell specific uptake activity of microbial products by uncultured chloroflexi by microautoradiography combined with fluorescence in situ hybridization. Environ Sci Technol. 2008;42(19):7380–6.

    Article  CAS  PubMed  Google Scholar 

  80. Kragelund C, Thomsen TR, Mielczarek AT, Nielsen PH. Eikelboom’s morphotype 0803 in activated sludge belongs to the genus Caldilinea in the phylum Chloroflexi. FEMS Microbiol Ecol. 2011;76(3):451–62.

    Article  CAS  PubMed  Google Scholar 

  81. Campbell AG, Schwientek P, Vishnivetskaya T, Woyke T, Levy S, Beall CJ, et al. Diversity and genomic insights into the uncultured Chloroflexi from the human microbiota. Environ Microbiol. 2014;16(9):2635–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. McGonigle JM, Lang SQ, Brazelton WJ. Genomic Evidence for Formate Metabolism by Chloroflexi as the Key to Unlocking Deep Carbon in Lost City Microbial Ecosystems. Appl Environ Microbiol. 2019;2020:1–12.

    Google Scholar 

  83. Hug LA, Castelle CJ, Wrighton KC, Thomas BC, Sharon I, Frischkorn KR, et al. Community genomic analyses constrain the distribution of metabolic traits across the Chloroflexi phylum and indicate roles in sediment carbon cycling. Microbiome. 2013;1(1):22.

    Article  PubMed  PubMed Central  Google Scholar 

  84. Wasmund K, Schreiber L, Lloyd KG, Petersen DG, Schramm A, Stepanauskas R, et al. Genome sequencing of a single cell of the widely distributed marine subsurface Dehalococcoidia, phylum Chloroflexi. ISME J. 2014;8(2):383–97.

    Article  CAS  PubMed  Google Scholar 

  85. Heider J, Ma K, Adams MWW. Purification, characterization, and metabolic function of tungsten- containing aldehyde ferredoxin oxidoreductase from the hyperthermophilic and proteolytic archaeon Thermococcus strain ES-1. J Bacteriol. 1995;177(16):4757–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Ali M, Shaw DR, Albertsen M, Saikaly PE. Comparative Genome-Centric Analysis of Freshwater and Marine ANAMMOX Cultures Suggests Functional Redundancy in Nitrogen Removal Processes. Front Microbiol. 2020;11(July):1–16.

    Google Scholar 

  87. Zhao Y, Jiang B, Tang X, Liu S. Metagenomic insights into functional traits variation and coupling effects on the anammox community during reactor start-up. Sci Total Environ. 2019;687(5):50–60.

    Article  CAS  PubMed  Google Scholar 

  88. Narihiro T, Terada T, Ohashi A, Kamagata Y, Nakamura K, Sekiguchi Y. Quantitative detection of previously characterized syntrophic bacteria in anaerobic wastewater treatment systems by sequence-specific rRNA cleavage method. Water Res. 2012;46(7):2167–75.

    Article  CAS  PubMed  Google Scholar 

  89. Singleton C, Petriglieri F, Kristensen J, Kirkegaard R, Michaelsen T, Andersen M, et al. Connecting structure to function with the recovery of over 1000 high-quality activated sludge metagenome-assembled genomes encoding full-length rRNA genes using long-read sequencing. Nat Commun. 2020;2021:1–13.

    Google Scholar 

  90. McIlroy SJ, Onetto CA, McIlroy B, Herbst FA, Dueholm MS, Kirkegaard RH, et al. Genomic and in Situ Analyses reveal the Micropruina spp as Abundant fermentative glycogen accumulating organisms in enhanced biological phosphorus removal systems. Front Microbiol. 2018;9:1004. Available from: pmc/articles/PMC5974061/. Cited 23 May 2021

    Article  PubMed  PubMed Central  Google Scholar 

  91. Lee SY. Bacterial polyhydroxyalkanoates. Biotechnology and bioengineering. 1996;49(1):1-14.

  92. Schlegel HG, Gottschalk G, Von Bartha R. Formation and utilization of poly-β-hydroxybutyric acid by knallgas bacteria (hydrogenomonas). Nature. 1961;191(4787):463–5.

    Article  CAS  PubMed  Google Scholar 

  93. Kalscheuer R, Steinbüchel A. A novel bifunctional wax ester synthase/acyl-CoA: Diacylglycerol acyltransferase mediates wax ester and triacylglycerol biosynthesis in Acinetobacter calcoaceticus ADP1. J Biol Chem. 2003;278(10):8075–82.

    Article  CAS  PubMed  Google Scholar 

  94. Kawaichi S, Yoshida T, Sako Y, Nakamuraa R, Nakamura R. Draft genome sequence of a heterotrophic facultative anaerobic thermophilic bacterium, Ardenticatena maritima strain 110ST. Genome Announc. 2015;3(5):1–2.

    Article  Google Scholar 

  95. Lawson CE, Wu S, Bhattacharjee AS, Hamilton JJ, McMahon KD, Goel R, et al. Metabolic network analysis reveals microbial community interactions in anammox granules. Nat Commun. 2017;8:1–12. https://doi.org/10.1038/ncomms15416.

    Article  CAS  Google Scholar 

  96. Watterson WJ, Tanyeri M, Watson AR, Cham CM, Shan Y, Chang EB, et al. Droplet-based high-throughput cultivation for accurate screening of antibiotic resistant gut microbes. Elife. 2020;9:1–22.

    Article  Google Scholar 

Download references

Acknowledgements

The authors express their gratitude to Uruguayan National Agency for Innovation and Research (ANII, POS_NAC_2015_1_109786) and the Postgraduate Academic Commission (CAP) for providing financial support to PB-W. PB-W acknowledges support from ANII through mobility grant MOV_CA_2017_1_138200. PB-W, CE, and AC were funded by SNI-ANII and PEDECIBA. PO and MESO are grateful for the support received through the Spanish Ministry of Economy and Competitiveness: “HIPATIA: Towards the implementation of the concept of biorefinery and energy self-sustainability in an urban wastewater treatment plant” (CTQ2017-82404-R). PO thanks to the National Agency for Research and Development (ANID) for her PhD Scholarship under the program DOCTORADO BECAS CHILE/2017 – 72180362.

Funding

This work was financially supported by Vaz Ferreira project I/FVF2017/135 and the Spanish Ministry of Economy and Competitiveness.

Author information

Authors and Affiliations

  1. Microbial Ecology Laboratory, Department of Microbial Biochemistry and Genomic, Biological Research Institute “Clemente Estable”, Avenida Italia 3318, CP: 11600, Montevideo, Uruguay

    Patricia Bovio-Winkler & Claudia Etchebehere

  2. Instituto de Investigaciones en Ingeniería Genética Y Biología Molecular “Dr Héctor N. Torres” (INGEBI-CONICET), Buenos Aires, Argentina

    Leandro D. Guerrero & Leonardo Erijman

  3. GENOCOV Research Group, Department of Chemical, Biological and Environmental Engineering, Escola d’Enginyeria, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain

    Pía Oyarzúa & María Eugenia Suárez-Ojeda

  4. Instituto Tecnológico Regional Centro Sur, Universidad Tecnológica, Francisco Antonio Maciel S/N, CP: 97000, Durazno, Uruguay

    Angela Cabezas

Authors
  1. Patricia Bovio-Winkler
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Leandro D. Guerrero
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Leonardo Erijman
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Pía Oyarzúa
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. María Eugenia Suárez-Ojeda
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Angela Cabezas
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Claudia Etchebehere
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

AC and CE conceived and designed the project. PB-W performed the metagenomic and genomic analysis with assistance of LG and LE. PB-W, AC, and CE discussed the results. PB-W wrote the manuscript. All the authors contributed to manuscript revision and approved the final manuscript.

Corresponding author

Correspondence to Claudia Etchebehere.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Bovio-Winkler, P., Guerrero, L.D., Erijman, L. et al. Genome-centric metagenomic insights into the role of Chloroflexi in anammox, activated sludge and methanogenic reactors. BMC Microbiol 23, 45 (2023). https://doi.org/10.1186/s12866-023-02765-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12866-023-02765-5

Keywords

BMC Microbiology

ISSN: 1471-2180

Contact us