In biofilms, microbial activities form gradients of substrates and electron acceptors, creating a complex landscape of microhabitats, often resulting in structured localization of the microbial populations present. To understand the dynamic interplay between and within these populations, quantitative measurements and statistical analysis of their localization patterns within the biofilms are necessary, and adequate automated tools for such analyses are needed. We have designed and applied new methods for fluorescencein situhybridization (FISH) and digital image analysis of directionally dependent (anisotropic) multispecies biofilms. A sequential-FISH approach allowed multiple populations to be detected in a biofilm sample. This was combined with an automated tool for vertical-distribution analysis by generatingin silicobiofilm slices and the recently developed Inflate algorithm for coaggregation analysis of microbial populations in anisotropic biofilms. As a proof of principle, we show distinct stratification patterns of the ammonia oxidizersNitrosomonas oligotrophasubclusters I and II and the nitrite oxidizerNitrospirasublineage I in three different types of wastewater biofilms, suggesting niche differentiation between theN. oligotrophasubclusters, which could explain their coexistence in the same biofilms. Coaggregation analysis showed thatN. oligotrophasubcluster II aggregated closer toNitrospirathan didN. oligotrophasubcluster I in a pilot plant nitrifying trickling filter (NTF) and a moving-bed biofilm reactor (MBBR), but not in a full-scale NTF, indicating important ecophysiological differences between these phylogenetically closely related subclusters. By using high-resolution quantitative methods applicable to any multispecies biofilm in general, the ecological interactions of these complex ecosystems can be understood in more detail.

Moving bed biofilm reactors (MBBRs) are increasingly used for nitrogen removal with nitritation-anaerobic ammonium oxidation (anammox) processes in wastewater treatment. Carriers provide protected surfaces where ammonia oxidizing bacteria (AOB) and anammox bacteria form complex biofilms. However, the knowledge about the organization of microbial communities in MBBR biofilms is sparse. We used new cryosectioning and imaging methods for fluorescence in situ hybridization (FISH) to study the structure of biofilms retrieved from carriers in a nitritation-anammox MBBR. The dimensions of the carrier compartments and the biofilm cryosections after FISH showed good correlation, indicating little disturbance of biofilm samples by the treatment. FISH showed that Nitrosomonas europaea/eutropha-related cells dominated the AOB and Candidatus Brocadia fulgida-related cells dominated the anammox guild. New carriers were initially colonized by AOB, followed by anammox bacteria proliferating in the deeper biofilm layers, probably in anaerobic microhabitats created by AOB activity. Mature biofilms showed a pronounced three-dimensional stratification where AOB dominated closer to the biofilm-water interface, whereas anammox were dominant deeper into the carrier space and towards the walls. Our results suggest that current mathematical models may be oversimplifying these three-dimensional systems and unless the multidimensionality of these systems is considered, models may result in suboptimal design of MBBR carriers.

Background: High-throughput amplicon sequencing of marker genes, such as the 16S rRNA gene in Bacteria and Archaea, provides a wealth of information about the composition of microbial communities. To quantify differences between samples and draw conclusions about factors affecting community assembly, dissimilarity indices are typically used. However, results are subject to several biases, and data interpretation can be challenging. The Jaccard and Bray-Curtis indices, which are often used to quantify taxonomic dissimilarity, are not necessarily the most logical choices. Instead, we argue that Hill-based indices, which make it possible to systematically investigate the impact of relative abundance on dissimilarity, should be used for robust analysis of data. In combination with a null model, mechanisms of microbial community assembly can be analyzed. Here, we also introduce a new software, qdiv, which enables rapid calculations of Hill-based dissimilarity indices in combination with null models. Results: Using amplicon sequencing data from two experimental systems, aerobic granular sludge (AGS) reactors and microbial fuel cells (MFC), we show that the choice of dissimilarity index can have considerable impact on results and conclusions. High dissimilarity between replicates because of random sampling effects make incidence-based indices less suited for identifying differences between groups of samples. Determining a consensus table based on count tables generated with different bioinformatic pipelines reduced the number of low-abundant, potentially spurious amplicon sequence variants (ASVs) in the data sets, which led to lower dissimilarity between replicates. Analysis with a combination of Hill-based indices and a null model allowed us to show that different ecological mechanisms acted on different fractions of the microbial communities in the experimental systems. Conclusions: Hill-based indices provide a rational framework for analysis of dissimilarity between microbial community samples. In combination with a null model, the effects of deterministic and stochastic community assembly factors on taxa of different relative abundances can be systematically investigated. Calculations of Hill-based dissimilarity indices in combination with a null model can be done in qdiv, which is freely available as a Python package (https://github.com/omvatten/qdiv). In qdiv, a consensus table can also be determined from several count tables generated with different bioinformatic pipelines.

Granulation of activated sludge is an increasingly important area within the field of wastewater treatment. Granulation is usually achieved by high hydraulic selection pressure, which results in the wash-out of slow settling particles. The effect of the harsh wash-out conditions on the granular sludge ecosystem is not yet fully understood, but different bacterial groups may be affected to varying degrees. In this study, we used high-throughput amplicon sequencing to follow the community composition in granular sludge reactors for 12 weeks, both in the granular phase and the suspended phase (effluent). The microbiome of the washed out biomass was similar but not identical to the microbiome of the granular biomass. Certain taxa (e.g. Flavobacterium spp. and Bdellovibrio spp.) had significantly (p < 0.05) higher relative abundance in the granules compared to the effluent. Fluorescence in situ hybridization images indicated that these taxa were mainly located in the interior of granules and therefore protected from erosion. Other taxa (e.g. Meganema sp. and Zooglea sp.) had significantly lower relative abundance in the granules compared to the effluent, and appeared to be mainly located on the surface of granules and therefore subject to erosion. Despite being washed out, these taxa were among the most abundant members of the granular sludge communities and were likely growing fast in the reactors. The ratio between relative abundance in the granular biomass and in the effluent did not predict temporal variation of the taxa in the reactors, but it did appear to predict the spatial location of the taxa in the granules.

The activated sludge process is one of the most widely used methods for treatment of wastewater and the microbial community composition in the sludge is important for the process operation. While the bacterial communities have been characterized in various activated sludge systems little is known about archaeal communities in activated sludge. The diversity and dynamics of the Archaea community in a full-scale activated sludge wastewater treatment plant were investigated by fluorescence in situ hybridization, terminal restriction fragment length polymorphism analysis and cloning and sequencing of 16S rRNA genes. The Archaea community was dominated by Methanosaeta-like species. During a 15 month period major changes in the community composition were only observed twice despite seasonal variations in environmental and operating conditions. Water temperature appeared to be the process parameter that affected the community composition the most. Several terminal restriction fragments also showed strong correlations with sludge properties and effluent water properties. The Archaea were estimated to make up 1.6% of total cell numbers in the activated sludge and were present both as single cells and colonies of varying sizes. The results presented here show that Archaea can constitute a constant and integral part of the activated sludge and that it can therefore be useful to include Archaea in future studies of microbial communities in activated sludge.

The aerobic granular sludge process is a promising technology for the removal of nutrients and organic contaminants from wastewater. However, a large amount of the sludge is often washed out during the start-up of granular reactors, which results in reduced process performance and a protracted start-up phase. In this study, the possibility of a rapid start-up of the nitrification process through a stepwise decrease of the settling time was investigated, and the bacterial population dynamics in two lab-scale sequencing batch reactors were studied. The results demonstrated that the stepwise decrease of the settling time enabled fast granulation and rapid start-up of the process. Small cores of granules were already observed after 10 days of operation, and the biomass was dominated by granules after 28 days. The removal of organic matter and ammonium was >95% after one day and 14 days, respectively. The bacterial community composition changed rapidly during the first 21 days, resulting in strongly reduced richness and evenness. The diversity increased at a later stage, and the bacterial community continued changing, albeit at a slower pace. The rate of the stepwise decrease in settling time strongly affected the abundance of nitrifying organisms, but not the general composition of the bacterial community. The results of this study support the idea that a stepwise decrease of the settling time is a successful strategy for the rapid start-up of aerobic granular sludge reactors.