Abstract
Advances in null-model approaches have resulted in a deeper understanding of community assembly mechanisms for a variety of complex microbiomes. One under-explored application is assembly of communities from the built-environment, especially during process disturbances. Anaerobic digestion for biological wastewater treatment is often underpinned by retaining millions of active granular biofilm aggregates. Flotation of granules is a major problem, resulting in process failure. Anaerobic aggregates were sampled from three identical bioreactors treating dairy wastewater. Microbiome structure was analysed using qPCR and 16S rRNA gene amplicon sequencing from DNA and cDNA. A comprehensive null-model approach quantified assembly mechanisms of floating and settled communities. Significant differences in diversity were observed between floating and settled granules, in particular, we highlight the changing abundances of Methanosaeta and Lactococcus. Both stochastic and deterministic processes were important for community assembly. Homogeneous selection was the primary mechanism for all categories, but dispersal processes also contributed. The lottery model was used to identify clade-level competition driving community assembly. Lottery “winners” were identified with different winners between floating and settled groups. Some groups changed their winner status when flotation occurred. Spirochaetaceae, for example, was only a winner in settled biomass (cDNA-level) and lost its winner status during flotation. Alternatively, Arcobacter butzerli gained winner status during flotation. This analysis provides a deeper understanding of changes that occur during process instabilities and identified groups which may be washed out—an important consideration for process control.
Introduction
Even Earth’s most extreme environments are teeming with complex microbial communities, performing vital eco-system functions. Microorganisms persist in nearly every environment, both natural and built, regulating global biogeochemical cycling of critical nutrients (; ). Understanding and measuring biodiversity and the way these complex communities assemble and continue to develop has intrigued ecologists for decades (). The use of null-models for identifying and quantifying community assembly has continued to gain traction with several advances in both methodology and interpretation (Presley et al., 2010; ; Stegen et al., 2012; Tucker et al., 2016; Zhou and Ning, 2017; Verster and Borenstein, 2018; ; Vass et al., 2020). Recently the assembly in rock pool communities was assessed using a combination of null-models, finding that dispersal limitation was a previously underestimated process driving assembly in these communities (Vass et al., 2020). Alternatively, a lottery-based approach was implemented to explore clade-based assembly in the human gut microbiome by identifying lottery “winners” which were out-competing other closely related taxa (Verster and Borenstein, 2018).
Using a combination of several approaches and ecological frameworks, we can begin to piece together the dynamics of assembly processes over space and time, and under varying environmental conditions. These processes include both stochastic and deterministic mechanisms. Stochastic processes include ecological drift driven by random birth-death events and random colonisation. Conversely, deterministic processes are driven by both abiotic (environmental filtering: changes in pH, temperature, salinity, etc.) and biotic factors (competition, facilitation, mutualism, predation, etc.) (). Moreover, we understand that several distinct mechanisms can act on the microbiome simultaneously—especially in complex communities—and that the relative proportion of these processes can change (Stegen et al., 2013; Vass et al., 2020). For example, several studies have identified shifts from stochasticity to determinism following an environmental disturbance (; Zhang et al., 2016). To date, there are conflicting reports as to the importance of these types of processes in engineered biological systems—complicated by the fact that each study uses a different ecological framework and methodology (Ofiţeru et al., 2010; Zhou et al., 2013; Vanwonterghem et al., 2014; ; ). However, as a whole, the understanding of how communities assemble in ecosystems of the built environment remains under-studied. Yet, they are wide-spread applications of biotechnology, as well as a source of highly replicated microbial communities (; Trego et al., 2020).
Anaerobic granules, for example, form through self-immobilisation of bacteria and archaea () and are crucial to the success of high-rate, upflow anaerobic systems such as the upflow anaerobic sludge bed (UASB) bioreactor (). Each individual granule contains the entire microbial community necessary to convert organic pollutants to a methane-based biogas via the anaerobic digestion (AD) process. Granule flotation, however, is a frequently reported problem associated with UASB bioreactors (; ; ; Wang et al., 2018) which can lead to washout of active granules and can dramatically reduce system capacity (Yoda and Nishimura, 1997). Furthermore, persistent flotation and washout can lead to severe biomass loss and, eventually, process failure ().
The causes of flotation are several-fold and complex, but a number of studies have reported that the main issue is attachment of biogas bubbles to the biofilm (; Yoda and Nishimura, 1997; ). This occurs when biogas within the granule is produced faster than it can be released (Rinzema et al., 1994; ; ; ; ; Wang et al., 2018; Singh et al., 2019). Interestingly, the archaeal and bacterial composition, and the location of these organisms within the biofilm seems to be an important factor, although again there are conflicting ideas about their role. Some studies report that aggregates with hydrophobic surfaces have a higher affinity for biogas adsorption, thus granules with methanogenic archaea located at the exterior, may be more susceptible to flotation (; ). Conversely, other studies have demonstrated that granules with methanogens located in the interior were more likely to float due to entrapped biogas within the biofilm (Saiki et al., 2002). Finally, others have demonstrated that excessive growth of Methanosaeta can lead to filamentous bulking, which entraps biogas and results in flotation (). Missing from these microbiome-based explanations is a deeper analysis of the differences between communities of floating and settled granules and the identification and quantification of the exact ecological mechanisms driving community assembly.
Hence, the aim of this study was twofold: (i) to examine several physiological characteristics of both floating and settled aggregates; and (ii) to resolve patterns in the microbial community structure and assembly that may facilitate our understanding of granule flotation. Not only did we uncover significant differences in microbiome diversity and structure between floating and settled aggregates, but we also found that the assembly patterns differed, resulting in a distinct “floating microbiome.”
Materials and Methods
Bioreactor Operation and Biomass Sampling
Anaerobic granules were sourced from a full-scale (1,500 m3) internal circulation reactor at Carbery Milk Products (CMP; Ballineen, Co., Cork, Ireland). At CMP the granules were used to treat ethanol production wastewater under mesophilic (37°C) conditions. The volatile solids (VS) concentration of the biomass was 91 g VS L–1. The granules were 0.8–3.2 mm in diameter, spherical, and black in colour. For this study they were used to inoculate triplicate, laboratory-scale (3.5-L working volume; Supplementary Figure 1) upflow anaerobic bioreactors, which were operated under identical conditions (). Briefly, the bioreactors were inoculated with 20 g VS L–1 of granular biomass. Synthetic dairy wastewater was continuously supplied, consisting of skimmed milk powder (2.5 g COD L–1) supplemented with both macro- and micronutrients (Shelton and Tiedje, 1984). The pH was buffered using 1.2 g L–1 NaHCO3. The operating temperature, controlled using external water jackets, was maintained at 37°C until day 115, when it was decreased for low-temperature operation to 15°C for the following 75 days. This decrease in temperature was to examine the feasibility of low-temperature anaerobic treatment of dairy wastewater, on which we have previously reported (). Separately, granular biomass was sampled from the settled and floating biomass layers from each of the triplicate bioreactors after 75 days of operation at 15°C, when sludge flotation was frequently observed. Samples for DNA/RNA extractions were immediately flash-frozen in liquid nitrogen and stored at −80°C.
Physico-Chemical Characterisation
Granule diameter was measured for individual granules using digital callipers (RS Components Ltd, Northants, United Kingdom; accurate to 0.03 mm). The total solids (TS) and VS concentrations of granules (n = 3) from both floating and settled biomass were determined using the standard loss-on-ignition technique (). VS was calculated as a percentage of the biomass wet weight. Settling velocity of granules (n = 341) was determined by measuring the time required for a single granule to travel 0.3 m down a clear, acrylic tube filled with deionised water. The settling velocity was the distance divided by the settling time. Density was calculated using Stokes’ law. Long chain fatty acid (LCFA) content from granular biomass collected from the floating and settled layers was determined on a Varian Saturn 2000 GC/MS system (Varian Inc., Walnut Creek, CA) using a method adapted from with details available in Supplementary Material. Since non-parametric tests are distribution-free tests, to find significant differences between floating and settled granules, we used the Kruskal-Wallis test.
DNA/RNA Co-extraction and cDNA Synthesis
DNA and RNA were co-extracted from granular biomass collected from floating (n = 6 DNA and n = 6 RNA) and settled (n = 6 DNA and n = 6 RNA) layers of the bioreactors (total of n = 24 samples). Specifically, two samples were used from each of the three replicated bioreactors to constitute the six for each category. For each sample, nucleic acids from 3 g of wet biomass, consisting of multiple granules, were extracted on ice following the phenol-chloroform based procedure previously described (; ). Concentrations were determined using a Qubit fluorometer (Invitrogen, Carlsbad, CA, United States). An aliquot of DNA was stored at −80°C. cDNA was then synthesised from the RNA. DNA was removed using the Turbo DNA-free kit (Ambion—Invitrogen, Carlsbad, CA, United States). PCR, with universal bacterial and archaeal primers 515F and 806R (), confirmed the samples to be DNA-free. cDNA was then synthesised using the SuperScript III Reverse Transcriptase Kit (Thermo Fisher Scientific, Waltham, MA, United States). Successful cDNA generation was confirmed by PCR amplification using primers 515F and 806R and cDNA was stored at −80°C.
Quantitative Real-Time Polymerase Chain Reaction (qPCR)
Bacterial and archaeal domains were separately targeted and quantified using qPCR. Reactions were performed on both DNA and cDNA samples from floating (n = 6) and settled biomass (n = 6) following the complete description by using bacterial primer pair 1369F and 1492R and Taqman probe TM1389F (Suzuki et al., 2000), and archaeal primer pair 787F and 1059R and the Taqman probe TM915F (Yu et al., 2005). Quantitative standard curves used Escherichia coli as a representative bacterial isolate and Methanosarcina barkeri as the representative archaeal isolate. Final gene copy numbers were determined per gram wet biomass. We used the non-parametric Mann-Whitney test to compare qPCR values between different categories.
High-Throughput Sequencing
The V4 region of the 16S rRNA gene was amplified using universal bacterial and archaeal primer set 515F and 806R () with indexed barcodes on the forward primer. Nucleic acids were normalised to 20 μg mL–1. Normalised samples were combined and run in triplicate on a 2% agarose gel. The ∼300 bp bands were excised and purified using the Wizard SV gel and PCR clean-up kit (Promega, Madison, Wisconsin, United States). Purified PCR products were normalised to 7.1 ng μL–1. Sequencing was performed on the Illumina MiSeq platform by the Centre for Genomic Research in the University of Liverpool (Liverpool, United Kingdom). The sequencing data from this study are available through the NCBI database under the project accession number PRJNA616223 with sample information available (Supplementary Table 1).
Bioinformatics and Statistical Analysis
Abundance tables were generated by constructing amplicon sequencing variants (ASVs) using the Qiime2 pipeline and the DADA2 algorithm () with details given at1. A total of 2 175 ASVs from n = 24 samples were identified, with summary statistics for reads per sample as follows: (1st Quantile: 30 203; Median: 35 630; Mean: 36 183; 3rd Quantile: 40 885; Maximum: 49 488). Within the workflow, qiime feature-classifier was used to classify the ASVs against SILVA SSU Ref NR database release v.132, and then qiime phylogeny align-to-tree-mafft-fasttree generated the rooted phylogenetic tree. The biom file for the ASVs was generated by combining the abundance table with taxonomy information using biom utility available in qiime2 workflow. In addition, we have removed contaminants such as chloroplasts and mitochondria as is recommended in taxonomy-based filtering of artifacts in the Qiime2 workflow2.
Next Picrust2 (), and its qiime2 plugin3 using the parameters –p-hsp-method pic –p-max-nsti 2 in qiime picrust2 full-pipeline, was used to find KEGG enzymes and MetaCyc pathway predictions. Although prediction process is highly dependent on the number of pathways available for the reference genomes, with Picrust2, it is now possible to generate accurate and putative metabolic maps at community level by virtue of its comprehensive database (∼20,000 genomes), which was not possible for its predecessor, Picrust1 (only 2 011 genomes). The algorithm consistently predicts pathways that have greater than 0.8 correlation with the actual pathways observed using shotgun metagenomic equivalents as highlighted by the authors (). Indeed, only a few (∼30 ASVs) from this study were not observed in the Picrust2 reference database. This increases our confidence on the prediction by virtue of high coverage.
For ensuing statistical analysis, the reads were rarefied to the sample with the minimum number of sample reads (23 510). This yielded a 24 (sample) × 1 829 (ASV) abundance table. Statistical analyses were performed in R (v.3.4.4) using the combined data generated from the bioinformatics as well as meta-data associated with the study (details available in Supplementary Material).
Null-Modelling
A multi-phasic null-model approach was used for comprehensive, quantitative, insights into the underlying ecological mechanisms driving community assembly. To apply null modelling techniques, we have used the full ASV table as obtained after preprocessing including removal of contaminants such as mitochondria and chloroplasts as is typical in amplicon workflows. Next, we wanted to understand the influence of the environment on microbial community assemblage. For this purpose, we used Nearest Taxon Index (NTI) to explore phylogenetic dispersion in the data. NTI is preferred because of presence of significant phylogenetic signal across short phylogenetic distances (Wang et al., 2013). Additionally, NTI is useful when phylogenetic signal cannot be measures, as was the case in this study due to lack of substantial trait data. Therefore, NTI helped determine whether the community was structured due to strong environmental pressure (local clustering in the phylogenetic tree). The NTI was calculated using mntd() and ses.mntd() functions from the picante package (). NTI represents the negative of the output from ses.mntd(). Additionally, it quantifies the number of standard deviations that separate the observed values from the mean of the null distribution (999 randomisation using null.model-“richness” in the ses.mntd() function and only considers taxa as either present or absent regardless of their relative abundance). Positive values indicate that species co-occur with more closely related species more frequently than expected by chance, with negative values suggesting otherwise. NTI measures tip-level divergences (putting more emphasis on terminal clades and is akin to “local” clustering) in phylogeny. For NTI, values > + 2 indicate strong environmental pressure, and values <−2 indicate strong competition among species as the driver of community structure.
Next, a stochasticity ratio calculation was implemented for which Jaccard (incidence-based) and Ružička (abundance-based) metrics were applied to determine the normalised stochasticity ratio (NST) based on author recommendations (). Taxa-Richness constraints of proportional-proportional (P-P) and proportional-fixed (P-F) were applied for each metric. To obtain significance for NST between treatments, we have used permutational multivariate ANOVA (henceforth referred to as PANOVA), as recommended by the authors ().
Next, the quantitative process estimates (QPE) method was used. This is based on an ecological framework that describes assembly processes in terms of selection (variable or homogenous), dispersal (dispersal limitation, or homogenising dispersal) or “undominated” mechanisms (Vellend, 2010; Stegen et al., 2015). Variable selection gives rise to high compositional differences in community structure due to different selective environmental conditions, while homogenous selection is when unchanging environmental conditions result in consistent selective pressure. Dispersal processes refer to the movement of organisms throughout space. High rates of dispersal result in similar communities, referred to as homogenising dispersal. Conversely, dispersal limitation increases differences in community composition resulting in more dissimilar communities. Conceptually, it occurs when low dispersal rates result in a high community turnover (Stegen et al., 2015). It is dispersal limitation which drives ecological drift (stochastic assembly).
The advantage of QPE is that it considers abundances as well as phylogeny of the ASVs. Following the previously described method (Stegen et al., 2015; ), deviation from the observed βMNTD (β-mean-nearest-taxon-distance) and the mean of the null distribution was evaluated using βNTI (β-nearest-taxon-index). When the observed value of βMNTD deviated significantly from the null expectation, the community was assembled by variable (βNTI >+2) or homogenous (βNTI <−2) selection processes. If the difference was not significant, the observed differences in phylogenetic composition were considered to be the result of dispersal mechanisms enabling ecological drift. These were differentiated using the abundance-based βRC and a Bray-Curtis dissimilarity metric for beta diversity. If the βRCbray > + 0.95, assembly was by dispersal limitation coupled with drift; if βRCbray < −0.95 then homogenising dispersal mechanisms contributed to community assembly; and if βRCbray was between −0.95 and + 0.95, community turnover was due to undominated mechanisms (i.e., dominated neither by dispersal nor selection processes).
Finally, we applied the competitive lottery model for clade-based community assembly (Verster and Borenstein, 2018). Briefly, the model describes the abundance of the most prevalent ASVs in the samples according to a group-based competitive lottery schema. Groups/clades are categorised phylogenetically, in our case at the family-level. It assumes that phylogenetically similar ASVs will have a similar gene content, metabolism, and preferential niche space—all resulting in high levels of competition. Lottery “winners” were defined as a member (ASV) that captures >90% of the groups abundance, and were identified. Families which included winners were plotted based on winner prevalence (the fraction of samples which include a winner ASV for that family) and winner diversity (the frequency that each ASV occurs as the winner in the samples from which winners are observed). A low diversity suggested that the same ASV was dominating that family in all samples while a high diversity suggested a more even spread of ASVs as winners in that group.
Results
In this study we sampled floating and settled granules from three replicated, low-temperature, laboratory-scale UASB bioreactors from one timepoint, when all three reactors were experiencing the flotation phenomenon. Both floating and settled biomass was compared in terms of density, settling velocity, LCFA composition and aggregate size. Additionally, the makeup of the microbial community was analysed by both qPCR and sequencing the 16S rRNA genes from DNA and cDNA of both floating and settled granules. A comprehensive range of multivariate analysis and null-models provided a novel way of understanding the differences between floating and settled microbial communities.
Characterisation of Floating and Settled Granules
Overall, floating and settled fractions had similar size distributions. Settled granules ranged in diameter from 0.2 (the smallest granules sampled) to 4.12 mm. Floating granules were only slightly smaller ranging from 0.2 to 3.89 mm. The mean granule diameter for floating and settled biomass was statistically similar (p = 0.73; Supplementary Figure 2). Settled granules displayed significantly higher settling velocities (p < 0.0001) and densities (p < 0.0001; Figures 1A,B). Total LCFA analysis showed large differences between floating and settled biomass, but a high standard deviation for the floating granules (±19.7) demonstrated high variability in LCFA concentrations from that category (Figure 1C), resulting in no significant differences (p = 0.13). Statistically significant differences were observed, however, for two of the individual LCFA (stearic acid, p = 0.05; oleic acid, p = 0.05), with higher concentrations in the settled biomass.
FIGURE 1
DNA-based gene copy numbers revealed significantly higher numbers of bacteria than archaea in both settled (p < 0.005) and floating (p < 0.05) biomass (Figure 1D). However, they did not differ statistically between floating and settled biomass. Gene copy numbers from cDNA showed similar gene copy numbers for bacteria in floating and settled granules, but detected significantly (p = 0.0001) higher numbers of archaea in the floating granules. Indeed, archaeal means approached 1.5 × 1010 and 3.8 × 109 for floating and settled biomass, respectively (Figure 1E).
Microbiome of Floating and Settled Granules
The rarefied richness (numbers of ASVs) was statistically similar between floating and settled granules (Figure 2A), while the balance of the community, measured as Shannon Entropy, was significantly lower in settled granules. This was observed in both DNA (p < 0.05) and cDNA (p < 0.01) samples (Figure 2B and Supplementary Figure 3). This suggests that while the distribution of the community was different, the actual number of observed taxa remained stable.
FIGURE 2
Significant differences in community composition were observed between floating and settled samples for both the total (DNA-based) and active (cDNA-based) communities using weighted UniFrac distances (Figure 2C; see S4 for unweighted UniFrac). Indeed, PERMANOVA confirmed that the nucleic acid and position of the biomass (i.e., whether it was floating or settled) significantly contribute to (p = 0.001) the variance between categories. Additional differences in diversity were observed between the total and active communities, which clustered separately (Figure 2C).
The makeup of the most abundant (top-25) taxa in the microbiome across all samples showed several interesting differences between floating and settled biomass (Figure 3). Notably, Lactococcus was relatively more abundant in settled biomass as opposed to floating. However, the relative abundance of fermentative bacteria, in general, fluctuated between the two categories. Anaerolineaceae, Lentimicrobiaceae, and Arcobacter were more relatively abundant in the floating granules compared to settled. Heat-trees identified entire clades within the microbiome which were enriched between categories (Figure 4). Comparison of the active communities (cDNA) of floating and settled biomass revealed that while the settled biomass had enriched groups of Methanobacterium, Clostridia, Actinobacteria, and Alphaproteobacteria, the floating biomass contained active clades of Methanomicrobiales, Bacteroidetes, Campylobacteria, and Gammaproteobacteria. Additionally, differences between DNA and cDNA samples were observed. Notably, the Euryarchaeota and Proteobacteria were enriched in the cDNA more than the DNA suggesting that they have a strong, potentially active role in the community. However, while many groups were dynamic, a core microbiome (Supplementary Figure 5) containing multiple taxa across all of the critical trophic groups required for complete AD persisted (spanning hydrolysis, acidogenesis, acetogenesis, and methanogenesis). This highlights the functional redundancy of the AD microbiome.
FIGURE 3
FIGURE 4
Subset analysis was used to identify a minimal group of ASVs, which in several combinations, could statistically explain the observed differences in community structure. Remarkably, only five ASVs (Lactococcus, Methanosaeta, Bacteroidetes vadinHA17, Arcobacter butzleri, and Lentimicrobiaceae), grouped in four different combinations (Figure 5A), explained between 26 and 30% of variation between floating and settled granules. A focus on the effects these ASVs have on the community structure and how they differed between settled and floating biomass could help explain flotation. For example, heat tree analysis of only these five ASVs revealed that when comparing cDNA of floating and settled biomass, Methanosaeta, a methane-producing archaeon, was enriched in the floating biomass. Additionally, Lactococcus, which was previously identified as a bacterium of interest (based on relative abundance; Figure 3), was further linked to the settled biomass (Figure 5).
FIGURE 5
Regressions of environmental/physico-chemical parameters against various diversity measures (Supplementary Figure 6) revealed that biomass position consistently had a significant, and complimentary influence on community structure—floating granules were consistently positively correlated with increased diversity. Notably, archaeal gene copy numbers were also positively correlated with diversity.
Putative functionality, based on Picrust2 algorithms was assessed. Differential pathway analysis was used to identify pathways and kegg orthologues which were significantly enriched between floating and settled biomass (Supplementary Figures 7–9). DNA and cDNA from each position identified 14 such pathways for each nucleic acid type (Supplementary Figures 7, 8). Putative differential pathways based on the DNA, were nearly exclusively (bar one) upregulated in the floating biomass (Supplementary Figure 8). Notably, the mevalonate pathway III, exclusive to the archaea, was putatively enriched in floating biomass (cDNA-based analysis), while pathways pertaining to lactose, or lactate were enriched in the settled biomass.
Community Assembly Processes
Through NTI analysis, strong environmental pressure was observed for all four categories although no significant changes between floating and settled biomass was detected (Figure 6A). In general, the observed values were lower in cDNA than in DNA, with significant (p < 0.01) differences in nucleic acid type for the settled biomass.
FIGURE 6
Stochasticity was quantified using the normalised stochasticity ratio (Figure 6B). In this case floating and settled categories showed differences in the percent stochasticity. According to NST (P-F), stochasticity accounted for 47–66% of assembly for floating granules, and 38–56% for settled communities. However, none of the differences between floating and settled biomass were significant (p = 0.09; PANOVA). The combination of both of these approaches showed that stochastic processes were relevant, but not the only ecological mechanism shaping the microbiome.
To further quantify and identify assembly processes, we used the QPE approach (Figure 6C) which is based on an ecological framework defined by selection (deterministic processes), dispersal (stochastic processes resulting in ecological drift), and undominated mechanisms. Here we observed that undominated processes (neither selection driven, nor dispersal driven) accounted for 13–33% of community turnover. Variable selection played no role, which is explained by the identical environmental conditions of the samples (same temperature, pH, salinity, etc.). Homogeneous selection was the dominant assembly mechanism in all categories, ranging from 40 to 80%. Clear differences in assembly processes were observed between settled and floating categories. Dispersal limitation from cDNA samples increased in relative importance from 0 to 20% between settled and floating biomass, respectively.
Finally, the lottery model used a different lens to assess clade-based community assembly. Lottery “winners” were identified (those making up >90% of community abundance within their clade) and were plotted based on winner prevalence and diversity (Supplementary Figure 11; genus-level analysis in Supplementary Figure 12). Notably, most families showing lottery-behaviour had low diversities, indicating that only one ASV was dominating as the winner. Two families showed intense lottery-like behaviour—having both high prevalence and diversity, indicating that winners were found in nearly all samples and that several different ASVs were selected as the winners. Most interesting, however, is the difference in winner behaviour between floating and settled samples (Figure 7). The lottery tree highlights winning clades to their most refined taxonomic level, also identifying clades with no lottery behaviour. Clades with listed groupings highlight those where lottery behaviour is changing between groups (i.e., where not all groups were identified for that clade). For example, Smithella was a winner in all groups, and thus lottery assembly did not differ due to flotation (no groups listed). Spirochaetaceae, however, was only a winner in settled biomass (cDNA) and lost its winner status during flotation. Conversely, Arcobacter butzerli gained winner status in the “floating microbiome.”
FIGURE 7
Discussion
Ecological Mechanisms Underpinning Granule Flotation
The combination of null models used in this study revealed that both stochastic and deterministic processes simultaneously drive community assembly—an idea that could reconcile conflicting reports from built environments (Ofiţeru et al., 2010; Zhou et al., 2013; Vanwonterghem et al., 2014; ; ). QPE analysis took this one step further to classify and quantify these processes in terms of selection and dispersal. Interestingly, homogenous selection was the most dominant ecological principle for all categories—where selection refers to ecological fitness and ongoing abiotic and biotic interactions shaping the community. The lottery model additionally suggested that these selection processes may have been governed by intense competition between species.
According to QPE, dispersal processes have a stronger role to play in floating granules. In view of homogeneous selection being a dominant process, as well as dispersal limitation, which was observed in floating samples, we propose the following hypotheses: (i) a lack of connectivity between positions (floating and settled) may result in high compositional turnover due to dispersal limitation and drift; or, (ii) that environmental selection pressure (non-optimal temperatures, pH fluctuations, mixing strategies, etc.) may make it harder for fermenters such as Lactococcus to colonise, even when there is connectivity and that these uncolonised granules may eventually float. Moreover, these findings add to the growing series of literature which identifies dispersal limitation as an important factor driving community turnover across various ecosystems (Stegen et al., 2013, 2015; ; Vass et al., 2020). It should be noted that typically a lack of phylogenetic signal leads to non-significant βNTI (required for QPE analysis). The observed significant values of βNTI supports the post hoc assumption of phylogenetic signal across short distances, which in turn, increases confidence in the results.
Finally, the lottery model for clade-based assembly identified lottery “winners” across a wide range of the microbiome. Interestingly, clades displaying winner behaviour generally showed low winner diversity indicating that within those clades only one ASV was dominating as the winner. Such organisms are deemed to be out-contending other members of their clade, fiercely competing for available nutrients within their niche space. The strongest lottery-like groups, however, have a high diversity and high winner prevalence (Verster and Borenstein, 2018), indicating a more variable selection of ASVs as winners within the clade. In the context of flotation, changing lottery behaviour gives key insights into how the “floating microbiome” is shaped. Several winners were present in both floating and settled aggregates. Such winners can be thought of as strong competitors regardless of flotation, and are perhaps more resilient to such eco-system disturbances. Winners, however, which are only present in either floating or settled granules indicated changing clade-based assembly mechanisms. Spirochaetaceae was a lottery winner (cDNA) in settled granules, but lost its winner status in the floating biomass. This suggests that the winners within Spirochaetaceae were not as fit in the new floating biomass. Conversely, some members of other families, such as Acidaminococcaceae, found their niche, out-competing the rest of the clade in the floating biomass. Such groups, in particular, give rise to the distinct “floating microbiome.”
A Distinct “Floating Microbiome”
The microbial communities of floating and settled granules were significantly different from one another. It is possible that the microbial community may play a contributing role in causing the observed reduction in density. Such ideas have been previously explored, for example, the role of methanogenic archaea (Saiki et al., 2002), specifically regarding the acetoclastic genus Methanosaeta. Here, Methanosaeta was identified during subset analysis, as one of only five ASVs, which in several combinations, statistically explain variances in community structure. Notably, cDNA analysis identified the enrichment of Methanosaeta in floating granules (Figure 5), compared to settled. Furthermore, exclusively archaeal pathways, such as the mevalonate pathway—responsible for the synthesis of isoprenoid compounds which are found in archaeal lipid membranes, fundamentally distinguishing them from bacteria—were putatively enriched in the floating granules. This pathway is presumed to be widespread amongst archaea, but has been identified amongst several methanogens (; Yoshida et al., 2020). Finally, archaeal gene copy numbers were significantly higher in floating granules according to cDNA analysis. The combination of these analyses suggests that Methanosaeta may be linked to flotation.
There are two likely mechanisms by which Methanosaeta, and other gas-producing organisms, may contribute to flotation. Both relate to the position such organisms occupy within the biofilm (i.e., on the surface, or in the interior). The primary hypothesis relates to the density; it is possible that granules having relatively higher numbers and/or proportions of active biogas producers in the biofilm interior were generating more biogas. These granules were therefore more likely to have trapped gas pockets, decreasing the biofilm density and leading to flotation. An alternative hypothesis suggests that although methanogenic archaea (primary gas-producing organisms in AD processes) are generally accepted to be located in the core of the granule (; Sekiguchi et al., 1999) they can also occupy the surface layers (Saiki et al., 2002). Previous work has explored how a hydrophilic coating of extracellular polymeric substances (EPS) can protect granules from flotation (). However, in our case, if hydrophobic Methanosaeta (; ) occupied biofilm surface, they may contribute to flotation.
Flotation Linked to Density and Not LCFA
Flotation is a serious problem for AD, resulting in the loss of active biomass, reducing system capacity (Yoda and Nishimura, 1997; ) and in severe cases, resulting in process failure (). The phenomenon is complex—and likely to be influenced by several factors. Flotation has been observed in many cases and many bioreactor types, with several reported causes including reduced biofilm density (; ; , ) sometimes attributed to trapped biogas (; Yoda and Nishimura, 1997; ; ; ), the influence of the microbiome (; ; Saiki et al., 2002; ) or the influence of LCFA (Rinzema et al., 1994; ; ; Singh et al., 2019). Specifically, several studies have attributed flotation to LCFA accumulation, especially when treating lipid-rich wastewaters such as dairy wastewater (Rinzema et al., 1994; ; Vidal et al., 2000; ; Singh et al., 2019; ). It is generally accepted that due to inherent hydrophobicity, the sorption of LCFA to the granule exterior contributes to flotation (). The results from our study indicate, however, that LCFA were unlikely to be the primary cause of flotation as more LCFA were found in the settled biomass than in the floating biomass. The findings from the physico-chemical characterisation revealed only that density and settling velocity were significantly reduced in the floating biomass, as other studies have previously reported (; Yoda and Nishimura, 1997; ; , ). This could likewise be due to trapped biogas in the interior of the biofilm.
Implications for Low-Temperature System Management
Flotation likely has several root causes, but always results in process instabilities including, but not limited to, biomass washout. The microbial community structure of floating and settled granules from this study differed significantly, as did the putative functional diversity and capacity at each position. This suggests that if floating biomass is preferentially washed out of the system, several groups, with particular functions may over time, decline in abundance. Namely, Anaerolineaceae, Arcobacter, Lentimicrobiaceae, along with high numbers of methanogens, were all identified in the floating biomass. Anaerobic microbiomes are known to be functionally redundant (; Zhu et al., 2020) however, the sustained loss of multiple key, and potentially active groups may still negatively impact system performance, especially for sensitive, low-temperature operation (). For example, Anaerolineaceae are a family of metabolically diverse bacteria, capable of semi-syntrophic primary, and secondary fermentations (Sekiguchi et al., 2003; , ) and may further be important for granular structure (Yamada et al., 2005). They have been previously identified as key players during low-temperature AD () and their reduction could disrupt low-temperature AD eco-system function. Finally, for systems optimising carbon transformations culminating in the production of methane-based biogas, decreased numbers of methanogenic archaea could have severe implications and a devastating impact on system performance (). Whilst this study was performed on laboratory-scale bioreactors, microbial ecology of full-scale systems is slightly more complex with respect to substrate composition, strength and loading rate, and feeding and heating regimes, as per our previous study (). Therefore, it is plausible to expect full-scale ecology to be slightly more dominated by selection processes than under controlled laboratory conditions.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA616223.
Author contributions
PM, VO’F, and FA designed the study. PM performed the qPCRs and prepared the sequencing libraries. AT collaborated with UI on the bioinformatics and statistical analysis and UI wrote the scripts to generate the raw figures. AT finalised figures and contributed to application of ecological theory. Results were interpreted by PM, AT, CN, TM, VO’F, and UI. AT drafted the manuscript. UI, FA, and VO’ revised the document. All authors approved the manuscript and agreed for accountability of the work therein.
Funding
This work was financially supported by the grants from the Higher Education Authority (HEA) of Ireland through: the Programme for Research at Third Level Institutions, Cycle 5 (PRTLI-5), co-funded by the European Regional Development Fund (ERDF); the Enterprise Ireland Technology Centres Programme (TC/2014/0016) and Science Foundation Ireland (14/IA/2371 and 16/RC/3889). UI was supported by the NERC Independent Research Fellowship (NE/L011956/1).
Acknowledgments
We would like to thank Carbery Milk Products for providing the anaerobic granules used in this study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.666584/full#supplementary-material
Footnotes
1.^https://github.com/umerijaz/tutorials/blob/master/qiime2_tutorial.md
References
1
AliM.WangZ.SalamK.HariA. R.PronkM.Van LoosdrechtM. C. M.et al (2019). Importance of species sorting and immigration on the bacterial assembly of different-sized aggregates in a full-scale aerobic granular sludge plant.Environ. Sci. Technol.538291–8301. 10.1021/acs.est.8b07303
2
AlphenaarP. A. (1994). Anaerobic granular sludge: characterization, and factors affecting its functioning.Wageningen: Landbouwuniversiteit te Wageningen.
3
AlvesM. M.PereiraM. A.SousaD. Z.CavaleiroA. J.PicavetM.SmidtH.et al (2009). Waste lipids to energy: how to optimize methane production from long-chain fatty acids (LCFA).Microb. Biotechnol.2538–550. 10.1111/j.1751-7915.2009.00100.x
4
APHA. (2005). Standard methods for the examination of water and wastewater. 21st ed.New York: American Public Health Association.
5
BatstoneD. J.KellerJ. (2001). Variation of bulk properties of anaerobic granules with wastewater type.Water Res.351723–1729. 10.1016/S0043-1354(00)00446-2
6
BolyenE.RideoutJ. R.DillonM. R.BokulichN. A.AbnetC. C.Al-GhalithG. A.et al (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2.Nat. Biotechnol.37852–857.
7
BottosE. M.KennedyD. W.RomeroE. B.FanslerS. J.BrownJ. M.BramerL. M.et al (2018). Dispersal limitation and thermodynamic constraints govern spatial structure of permafrost microbial communities.FEMS Microbiol. Ecol.94:fiy110. 10.1093/femsec/fiy110
8
CampanaroS.TreuL.Rodriguez-RL. M.KovalovszkiA.ZielsR. M.MausI.et al (2019). The anaerobic digestion microbiome: a collection of 1600 metagenome-assembled genomes shows high species diversity related to methane production.bioRxiv2019:680553. 10.1101/680553
9
CamposJ. L.del RíoA. V.PedrousoA.RauxP.GiustinianovichcE. A.Mosquera-CorralA. (2017). Granular biomass floatation: a simple kinetic/stoichiometric explanation.Chem. Eng. J.31163–71.
10
CaporasoJ. G.LauberC. L.WaltersW. A.Berg-LyonsD.LozuponeC. A.TurnbaughP. J.et al (2011). Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample.Proc. Natl. Acad. Sci.1084516–4522. 10.1073/pnas.1000080107
11
CardinaleB. J.DuffyJ. E.GonzalezA.HooperD. U.PerringsC.VenailP.et al (2012). Biodiversity loss and its impact on humanity.Nature48659–67.
12
ChaseJ. M.KraftN. J. B.SmithK. G.VellendM.InouyeB. D. (2011). Using null models to disentangle variation in community dissimilarity from variation in α-diversity.Ecosphere21–11.
13
ChenJ.JiQ.ZhengP.ChenT.WangC.MahmoodQ. (2010). Floatation and control of granular sludge in a high-rate anammox reactor.Water Res.443321–3328.
14
ConnellyS.ShinS. G.DillonR. J.IjazU. Z.QuinceC.SloanW. T.et al (2017). Bioreactor scalability: laboratory-scale bioreactor design influences performance, ecology, and community physiology in expanded granular sludge bed bioreactors.Front. Microbiol.8:664. 10.3389/fmicb.2017.00664
15
Cuervo-LópezF. M.MartinezF.Gutiérrez-RojasM.NoyolaR. A.GómezJ. (1999). Effect of nitrogen loading rate and carbon source on denitrification and sludge settleability in upflow anaerobic sludge blanket (UASB) reactors.Water Sci. Technol.40123–130.
16
DaffonchioD.ThaveesriJ.VerstraeteW. (1995). Contact angle measurement and cell hydrophobicity of granular sludge from upflow anaerobic sludge bed reactors.Appl. Environ. Microbiol.613676–3680.
17
de BeerD.O’FlahartyV.ThaveesriJ.LensP.VerstraeteW. (1996). Distribution of extracellular polysaccharides and flotation of anaerobic sludge.Appl. Microbiol. Biotechnol.46197–201.
18
DemirelB.SchererP. (2008). The roles of acetotrophic and hydrogenotrophic methanogens during anaerobic conversion of biomass to methane: a review.Rev. Environ. Sci. Bio/Technol.7173–190.
19
DouglasG. M.MaffeiV. J.ZaneveldJ.YurgelS. N.BrownJ. R.TaylorC. M.et al (2019). PICRUSt2: an improved and extensible approach for metagenome inference.BioRxiv2019:672295.
20
EftaxiasA.DiamantisV.MichailidisC.StamatelatouK.AivasidisA. (2020). Comparison of anaerobic digesters performance treating palmitic, stearic and oleic acid: determination of the LCFA kinetic constants using ADM1.Bioprocess Biosyst. Eng.431329–1338.
21
FritscheW.HofrichterM. (2000). “Aerobic Degradation by Microorganisms,” in Environmental Processes - soil decontamination, ed.KleinJ. (Weinheim, Germany: Wiley-VCH), 146–155.
22
GastonK. J. (2000). Global patterns in biodiversity.Nature405220–227.
23
GriffithsR. I.WhiteleyA. S.O’DonnellA. G.BaileyM. J. (2000). Rapid method for coextraction of DNA and RNA from natural environments for analysis of ribosomal DNA- and rRNA-based microbial community composition.Appl. Environ. Microbiol.665488–5491.
24
HalalshehM.KoppesJ.Den ElzenJ.ZeemanG.FayyadM.LettingaG. (2005). Effect of SRT and temperature on biological conversions and the related scum-forming potential.Water Res.392475–2482.
25
HwuC.-S.TsengS.-K.YuanC.-Y.KulikZ.LettingaG. (1998). Biosorption of long-chain fatty acids in UASB treatment process.Water Res.321571–1579.
26
KazievaE.YamamotoY.TajimaY.YokoyamaK.KatashkinaJ.NishioY. (2017). Characterization of feedback-resistant mevalonate kinases from the methanogenic archaeons Methanosaeta concilii and Methanocella paludicola.Microbiology1631283–1291. 10.1099/mic.0.000510
27
KeatingC.HughesD.MahonyT.CysneirosD.IjazU. Z.SmithC. J.et al (2018). Cold adaptation and replicable microbial community development during long-term low-temperature anaerobic digestion treatment of synthetic sewage.FEMS Microbiol. Ecol.94:fiy095. 10.1093/femsec/fiy095
28
KembelS. W.CowanP. D.HelmusM. R.CornwellW. K.MorlonH.AckerlyD. D.et al (2010). Picante: R tools for integrating phylogenies and ecology.Bioinformatics261463–1464. 10.1093/bioinformatics/btq166
29
LettingaG. (1995). Anaerobic digestion and wastewater treatment systems.Antonie Van Leeuwenhoek673–28. 10.1007/bf00872193
30
LettingaG.RoersmaR.GrinP. (1983). Anaerobic treatment of raw domestic sewage at ambient temperatures using a granular bed UASB reactor.Biotechnol. Bioeng.251701–1723. 10.1002/bit.260250703
31
LeventhalG. E.BoixC.KuechlerU.EnkeT. N.SliwerskaE.HolligerC.et al (2018). Strain-level diversity drives alternative community types in millimetre-scale granular biofilms.Nat. Microbiol.3:1295. 10.1038/s41564-018-0242-3
32
LiJ.HuB.ZhengP.QaisarM.MeiL. (2008). Filamentous granular sludge bulking in a laboratory scale UASB reactor.Bioresour. Technol.993431–3438. 10.1016/j.biortech.2007.08.005
33
LiW.ZhengP.JiJ.ZhangM.GuoJ.ZhangJ.et al (2014). Floatation of granular sludge and its mechanism: a key approach for high-rate denitrifying reactor.Bioresour. Technol.152414–419. 10.1016/j.biortech.2013.11.056
34
LiuY.TayJ. H. (2004). State of the art of biogranulation technology for wastewater treatment.Biotechnol. Adv.22533–563. 10.1016/j.biotechadv.2004.05.001
35
LoceyK. J.LennonJ. T. (2016). Scaling laws predict global microbial diversity.Proc. Natl. Acad. Sci.1135970–5975. 10.1073/pnas.1521291113
36
LuH.ZhengP.JiQ.ZhangH.JiJ.WangL.et al (2012). The structure, density and settlability of anammox granular sludge in high-rate reactors.Bioresour. Technol.123312–317. 10.1016/j.biortech.2012.07.003
37
LuX.ZhenG.EstradaA. L.ChenM.NiJ.HojoT.et al (2015). Operation performance and granule characterization of upflow anaerobic sludge blanket (UASB) reactor treating wastewater with starch as the sole carbon source.Bioresour. Technol.180264–273. 10.1016/j.biortech.2015.01.010
38
MacLeodF. A.GuiotS. R.CostertonJ. W. (1990). Layered structure of bacterial aggregates produced in an upflow anaerobic sludge bed and filter reactor.Appl. Environ. Microbiol.561598–1607. 10.1128/aem.56.6.1598-1607.1990
39
McAteerP. G.TregoA. C.ThornC.MahonyT.AbramF.O’FlahertyV. (2020). Reactor configuration influences microbial community structure during high-rate, low-temperature anaerobic treatment of dairy wastewater.Bioresour. Technol.2020:123221. 10.1016/j.biortech.2020.123221
40
NarihiroT.NobuM. K.KimN.KamagataY.LiuW. (2015). The nexus of syntrophy-associated microbiota in anaerobic digestion revealed by long-term enrichment and community survey.Environ. Microbiol.171707–1720. 10.1111/1462-2920.12616
41
NarihiroT.TeradaT.OhashiA.KamagataY.NakamuraK.SekiguchiY. (2012). Quantitative detection of previously characterized syntrophic bacteria in anaerobic wastewater treatment systems by sequence-specific rRNA cleavage method.Water Res.462167–2175. 10.1016/j.watres.2012.01.034
42
NevesL.PereiraM. A.MotaM.AlvesM. M. (2009). Detection and quantification of long chain fatty acids in liquid and solid samples and its relevance to understand anaerobic digestion of lipids.Bioresour. Technol.10091–96. 10.1016/j.biortech.2008.06.018
43
NingD.DengY.TiedjeJ. M.ZhouJ. (2019). A general framework for quantitatively assessing ecological stochasticity.Proc. Natl. Acad. Sci.11616892–16898. 10.1073/pnas.1904623116
44
OfiţeruI. D.LunnM.CurtisT. P.WellsG. F.CriddleC. S.FrancisC. A.et al (2010). Combined niche and neutral effects in a microbial wastewater treatment community.Proc. Natl. Acad. Sci.10715345–15350. 10.1073/pnas.1000604107
45
PresleyS. J.HigginsC. L.WilligM. R. (2010). A comprehensive framework for the evaluation of metacommunity structure.Oikos119908–917. 10.1111/j.1600-0706.2010.18544.x
46
RinzemaA.BooneM.van KnippenbergK.LettingaG. (1994). Bactericidal effect of long chain fatty acids in anaerobic digestion.Water Environ. Res.6640–49. 10.2175/wer.66.1.7
47
SaikiY.IwabuchiC.KatamiA.KitagawaY. (2002). Microbial analyses by fluorescence in situ hybridization of well-settled granular sludge in brewery wastewater treatment plants.J. Biosci. Bioeng.93601–606. 10.1016/s1389-1723(02)80244-4
48
SekiguchiY.KamagataY.NakamuraK.OhashiA.HaradaH. (1999). Fluorescence in situ hybridization using 16S rRNA-targeted oligonucleotides reveals localization of methanogens and selected uncultured bacteria in mesophilic and thermophilic sludge granules.Appl. Environ. Microbiol.651280–1288. 10.1128/aem.65.3.1280-1288.1999
49
SekiguchiY.YamadaT.HanadaS.OhashiA.HaradaH.KamagataY. (2003). 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.531843–1851. 10.1099/ijs.0.02699-0
50
SheltonD. R.TiedjeJ. M. (1984). Isolation and partial characterization of bacteria in an anaerobic consortium that mineralizes 3-chlorobenzoic acid.Appl. Environ. Microbiol.48840–848. 10.1128/aem.48.4.840-848.1984
51
SinghS.Rinta-KantoJ. M.KettunenR.TolvanenH.LensP.CollinsG.et al (2019). Anaerobic treatment of LCFA-containing synthetic dairy wastewater at 20°C: Process performance and microbial community dynamics.Sci. Total Environ.691960–968. 10.1016/j.scitotenv.2019.07.136
52
StegenJ. C.LinX.FredricksonJ. K.ChenX.KennedyD. W.MurrayC. J.et al (2013). Quantifying community assembly processes and identifying features that impose them.ISME J.72069–2079. 10.1038/ismej.2013.93
53
StegenJ. C.LinX.FredricksonJ. K.KonopkaA. E. (2015). Estimating and mapping ecological processes influencing microbial community assembly.Front. Microbiol.6:370. 10.3389/fmicb.2015.00370
54
StegenJ. C.LinX.KonopkaA. E.FredricksonJ. K. (2012). Stochastic and deterministic assembly processes in subsurface microbial communities.Isme J.6:1653. 10.1038/ismej.2012.22
55
SuzukiM. T.TaylorL. T.DeLongE. F. (2000). Quantitative analysis of small-subunit rRNA genes in mixed microbial populations via 5′-nuclease assays.Appl. Environ. Microbiol.664605–4614. 10.1128/aem.66.11.4605-4614.2000
56
TregoA. C.MillsS.CollinsG. (2020). Granular biofilms: function, application, and new trends as model microbial communities.Crit. Rev. Environ. Sci. Technol.20201–24. 10.1080/10643389.2020.1769433
57
TuckerC. M.ShoemakerL. G.DaviesK. F.NemergutD. R.MelbourneB. A. (2016). Differentiating between niche and neutral assembly in metacommunities using null models of β-diversity.Oikos125778–789. 10.1111/oik.02803
58
VanwonterghemI.JensenP. D.DennisP. G.HugenholtzP.RabaeyK.TysonG. W. (2014). Deterministic processes guide long-term synchronised population dynamics in replicate anaerobic digesters.ISME J.82015–2028. 10.1038/ismej.2014.50
59
VassM.SzékelyA. J.LindströmE. S.LangenhederS. (2020). Using null models to compare bacterial and microeukaryotic metacommunity assembly under shifting environmental conditions.Sci. Rep.102455. 10.1038/s41598-020-59182-1
60
VellendM. (2010). Conceptual synthesis in community ecology.Q. Rev. Biol.85183–206. 10.1086/652373
61
VersterA. J.BorensteinE. (2018). Competitive lottery-based assembly of selected clades in the human gut microbiome.Microbiome6:186. 10.1186/s40168-018-0571-8
62
VidalG.CarvalhoA.MendezR.LemaJ. M. (2000). Influence of the content in fats and proteins on the anaerobic biodegradability of dairy wastewaters.Bioresour. Technol.74231–239. 10.1016/s0960-8524(00)00015-8
63
WangB.WuD.ZhangX.MackeyH. R.ChenG.-H. (2018). Sludge flotation, its causes and control in granular sludge upflow reactors.Appl. Microbiol. Biotechnol.1026383–6392. 10.1007/s00253-018-9131-1
64
WangJ.ShenJ. I.WuY.TuC.SoininenJ.StegenJ. C.et al (2013). Phylogenetic beta diversity in bacterial assemblages across ecosystems: deterministic versus stochastic processes.ISME J.71310–1321. 10.1038/ismej.2013.30
65
YamadaT.SekiguchiY.ImachiH.KamagataY.OhashiA.HaradaH. (2005). Diversity, localization, and physiological properties of filamentous microbes belonging to chloroflexi subphylum I in mesophilic and thermophilic methanogenic sludge granules.Appl. Environ. Microbiol.717493–7503. 10.1128/aem.71.11.7493-7503.2005
66
YodaM.NishimuraS. (1997). Controlling granular sludge floatation in UASB reactors.Water Sci. Technol.36:165. 10.2166/wst.1997.0588
67
YoshidaR.YoshimuraT.HemmiH. (2020). Reconstruction of the “archaeal” mevalonate pathway from the methanogenic archaeon &It;span class="named-content genus-species" id="named-content-1">methanosarcina mazei&It;/span> in &It;span class="named-content genus-speci.Appl. Environ. Microbiol.86e2889–19. 10.1128/AEM.02889-19
68
YuY.LeeC.HwangS. (2005). Analysis of community structures in anaerobic processes using a quantitative real-time PCR method.Water Sci. Technol.5285–91. 10.2166/wst.2005.0502
69
ZhangX.JohnstonE. R.LiuW.LiL.HanX. (2016). Environmental changes affect the assembly of soil bacterial community primarily by mediating stochastic processes.Glob. Chang. Biol.22198–207. 10.1111/gcb.13080
70
ZhouJ.LiuW.DengY.JiangY.-H.XueK.HeZ.et al (2013). Stochastic assembly leads to alternative communities with distinct functions in a bioreactor microbial community.mBio4:e00584–12.
71
ZhouJ.NingD. (2017). Stochastic community assembly: does it matter in microbial ecology?Microbiol. Mol. Biol. Rev.81:e00002–e17. 10.1128/MMBR.00002-17
72
ZhuX.CampanaroS.TreuL.SeshadriR.IvanovaN.KougiasP. G.et al (2020). Metabolic dependencies govern microbial syntrophies during methanogenesis in an anaerobic digestion ecosystem.Microbiome8:22. 10.1186/s40168-019-0780-9
Summary
Keywords
anaerobic digestion, community assembly, low-temperature anaerobic digestion, sludge flotation, dairy wastewater
Citation
Trego AC, McAteer PG, Nzeteu C, Mahony T, Abram F, Ijaz UZ and O’Flaherty V (2021) Combined Stochastic and Deterministic Processes Drive Community Assembly of Anaerobic Microbiomes During Granule Flotation. Front. Microbiol. 12:666584. doi: 10.3389/fmicb.2021.666584
Received
10 February 2021
Accepted
20 April 2021
Published
14 May 2021
Volume
12 - 2021
Edited by
Sabine Kleinsteuber, Helmholtz Centre for Environmental Research (UFZ), Germany
Reviewed by
Oskar Modin, Chalmers University of Technology, Sweden; Seung Gu Shin, Gyeongnam National University of Science and Technology, South Korea; James C. Stegen, Pacific Northwest National Laboratory (DOE), United States
Updates
Copyright
© 2021 Trego, McAteer, Nzeteu, Mahony, Abram, Ijaz and O’Flaherty.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Umer Zeeshan Ijaz, umer.ijaz@glasgow.ac.ukVincent O’Flaherty, vincent.oflaherty@nuigalway.ie
†These authors have contributed equally to this work
This article was submitted to Microbiotechnology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.