Effects of Predation by Protists on Prokaryotic Community Function, Structure, and Diversity in Anaerobic Granular Sludge

Unifying concepts in methanogenic, aerobic, and anammox sludge granulation

The retention of dense and well-functioning microbial biomass is crucial for effective pollutant removal in several biological wastewater treatment technologies. High solids retention is often achieved through aggregation of microbial communities into dense, spherical aggregates known as granules, which were initially discovered in the 1980s. These granules have since been widely applied in upflow anaerobic digesters for waste-to-energy conversions. Furthermore, granular biomass has been applied in aerobic wastewater treatment and anaerobic ammonium oxidation (anammox) technologies. The mechanisms underpinning the formation of methanogenic, aerobic, and anammox granules are the subject of ongoing research. Although each granule type has been extensively studied in isolation, there has been a lack of comparative studies among these granulation processes. It is likely that there are some unifying concepts that are shared by all three sludge types. Identifying these unifying concepts could allow a unified theory of granulation to be formed. Here, we review the granulation mechanisms of methanogenic, aerobic, and anammox granular sludge, highlighting several common concepts, such as the role of extracellular polymeric substances, cations, and operational parameters like upflow velocity and shear force. We have then identified some unique features of each granule type, such as different internal structures, microbial compositions, and quorum sensing systems. Finally, we propose that future research should prioritize aspects of microbial ecology, such as community assembly or interspecies interactions in individual granules during their formation and growth.

Identification and cultivation of anaerobic bacterial scavengers of dead cells

The cycle of life and death and Earth's carbon cycle(s) are intimately linked, yet how bacterial cells, one of the largest pools of biomass on Earth, are recycled back into the carbon cycle remains enigmatic. In particular, no bacteria capable of scavenging dead cells in oxygen-depleted environments have been reported thus far. In this study, we discover the first anaerobes that scavenge dead cells and the two isolated strains use distinct strategies. Based on live-cell imaging, transmission electron microscopy, and hydrolytic enzyme assays, one strain (designated CYCD) relied on cell-to-cell contact and cell invagination for degrading dead food bacteria where as the other strain (MGCD) degraded dead food bacteria via excretion of lytic extracellular enzymes. Both strains could degrade dead cells of differing taxonomy (bacteria and archaea) and differing extents of cell damage, including those without artificially inflicted physical damage. In addition, both depended on symbiotic metabolic interactions for maximizing cell degradation, representing the first cultured syntrophic Bacteroidota. We collectively revealed multiple symbiotic bacterial decomposition routes of dead prokaryotic cells, providing novel insight into the last step of the carbon cycle.

Role of Bacterial Community Composition as a Driver of the Small-Sized Phytoplankton Community Structure in a Productive Coastal System

We present here the first detailed description of the seasonal patterns in bacterial community composition (BCC) in shelf waters off the Ría de Vigo (Spain), based on monthly samplings during 2 years. Moreover, we studied the relationship between bacterial and small-sized eukaryotic community composition to identify potential biotic interactions among components of these two communities. Bacterial operational taxonomic unit (OTU) richness and diversity systematically peaked in autumn-winter, likely related to low resource availability during this period. BCC showed seasonal and vertical patterns, with Rhodobacteraceae and Flavobacteriaceae families dominating in surface waters, and SAR11 clade dominating at the base of the photic zone (30 m depth). BCC variability was significantly explained by environmental variables (e.g., temperature of water, solar radiation, or dissolved organic matter). Interestingly, a strong and significant correlation was found between BCC and small-sized eukaryotic community composition (ECC), which suggests that biotic interactions may play a major role as structuring factors of the microbial plankton in this productive area. In addition, co-occurrence network analyses revealed strong and significant, mostly positive, associations between bacteria and small-sized phytoplankton. Positive associations likely result from mutualistic relationships (e.g., between Dinophyceae and Rhodobacteraceae), while some negative correlations suggest antagonistic interactions (e.g., between Pseudo-nitzchia sp. and SAR11). These results support the key role of biotic interactions as structuring factors of the small-sized eukaryotic community, mostly driven by positive associations between small-sized phytoplankton and bacteria.

Morphology and phylogeny of two anaerobic freshwater ciliates: Brachonella comma sp. nov. and the widely-distributed but little-known caenomorphid, Ludio parvulus Penard, 1922

Hypoxic, sulfidic freshwater sediments typically support a diffuse consortium of distinctive ciliated protists, including caenomorphids, metopids, and odontostomatids among others. A recent resurgence of interest in these important members of sapropelic food webs has resulted in the description of many new species and an effort, still in its infancy, to characterize them from a morphologic, molecular, and metabolic standpoint and to determine their phylogenetic relationships. Their seemingly invariable association with prokaryotic endosymbionts and, less commonly, ectosymbionts, has become a focus for many researchers. In this report, based on morphologic and molecular data we describe a Brachonella species (Ciliophora, Metopida) new to science and analyze its phylogeny. We also provide a morphologic and molecular characterization of the smallest representative of the Caenomorphidae Poche, 1913, Ludio parvulus Penard, 1922. The phylogenetic analysis confirms the inclusion of this species in the Caenomorphidae.

Mass-immigration determines the assembly of activated sludge microbial communities

Wastewater treatment plants are engineering technologies used worldwide to protect the environment and human health. Microbial communities sustain these plants, so it is crucial to know the key factors responsible for the community assembly. We show, in contrast to existing understanding, that microbial immigration largely controls the community structure in these plants and that the fate (growth or death) of immigrating species in the plants is controlled by local factors. The community structure was quantitatively predicted by the immigrating microbial community, highlighting the need to revise the way we today understand, design, and manage microbial communities in wastewater treatment plants.

The assembly of bacterial communities in wastewater treatment plants (WWTPs) is affected by immigration via wastewater streams, but the impact and extent of bacterial immigrants are still unknown. Here, we quantify the effect of immigration at the species level in 11 Danish full-scale activated sludge (AS) plants. All plants have different source communities but have very similar process design, defining the same overall environmental growth conditions. The AS community composition in each plant was strongly reflected by the corresponding influent wastewater (IWW) microbial composition. Most species in AS across the plants were detected and quantified in the corresponding IWW, allowing us to identify their fate in the AS: growing, disappearing, or surviving. Most of the abundant species in IWW disappeared in AS, so their presence in the AS biomass was only due to continuous mass-immigration. In AS, most of the abundant growing species were present in the IWW at very low abundances. We predicted the AS species abundances from their abundance in IWW by using a partial least square regression model. Some species in AS were predicted by their own abundance in IWW, while others by multiple species abundances. Detailed analyses of functional guilds revealed different prediction patterns for different species. We show, in contrast to the present understanding, that the AS microbial communities were strongly controlled by the IWW source community and could be quantitatively predicted by taking into account immigration. This highlights a need to revise the way we understand, design, and manage the microbial communities in WWTPs.

Multi-omics analysis reveals structure and function of biofilm microbial communities in a pre-denitrification biofilter

The highly complex microbial communities in biofilm play crucial roles in the pollutant removal performance of wastewater treatment plants (WWTPs). In the present study, using multi-omics analysis, we studied microbial structure, key enzymes, functional traits, and key metabolic pathways of pre-denitrification biofilter in an urban WWTP in China. The analysis results of metagenomic and metaproteomic showed that Betaproteobacteria and Flavobacteriia were dominant in biofilms. The integrated metagenomic and metaproteomic data showed that the expression of nitrogen metabolism genes was high, and the high proportion of denitrification module indicating that denitrification was the main nitrogen removal pathway. The most abundant denitrifying bacterial genera were: Dechloromonas, Acidovorax, Bosea, Polaromonas, and Chryseobacterium. And microorganisms with denitrification potential may not be able to denitrify in the actual operation of the filter. The integrated analysis of metaproteomic and metabolomic showed that there was a correlation between biofilm microorganisms and metabolites. Metabolomic analysis indicated that metabolic profiles of biofilms varied with layer height. This study provides the first detailed microbial communities and metabolic profiles in a full-scale pre-denitrification biofilter and clarifies the mechanism of denitrification.

Food selectivity of anaerobic protists and direct evidence for methane production using carbon from prey bacteria by endosymbiotic methanogen

Anaerobic protists are major predators of prokaryotes in anaerobic ecosystems. However, little is known about the predation behavior of anaerobic protists because almost none have been cultured. In particular, these characteristics of anaerobic protists in the phyla Metamonada and Cercozoa have not been reported previously. In this study, we isolated three anaerobic protists, Cyclidium sp., Trichomitus sp., and Paracercomonas sp., from anaerobic granular sludge in an up-flow anaerobic sludge blanket reactor used to treat domestic sewage. Ingestion and digestion of food bacteria by anaerobic protists with or without endosymbiotic methanogens were demonstrated using tracer experiments with green fluorescent protein and a stable carbon isotope. These tracer experiments also demonstrated that Cyclidium sp. supplied CO₂ and hydrogen to endosymbiotic methanogens. While Cyclidium sp. and Trichomitus sp. ingested both Gram-negative and -positive bacteria, Paracercomonas sp. could only take up Gram-negative bacteria. Archaeal cells such as Methanobacterium beijingense and Methanospirillum hungatei did not support the growth of these protists. Metabolite patterns of all three protists differed and were influenced by food bacterial species. These reported growth rates, ingestion rates, food selectivity, and metabolite patterns provide important insights into the ecological roles of these protists in anaerobic ecosystems.

Factors shaping community patterns of protists and bacteria on a European scale

Factors shaping community patterns of microorganisms are controversially discussed. Physical and chemical factors certainly limit the survival of individual taxa and maintenance of diversity. In recent years, a contribution of geographic distance and dispersal barriers to distribution patterns of protists and bacteria has been demonstrated. Organismic interactions such as competition, predation and mutualism further modify community structure and maintenance of distinct taxa. Here, we address the relative importance of these different factors in shaping protists and bacterial communities on a European scale using high-throughput sequencing data obtained from lentic freshwater ecosystems. We show that community patterns of protists are similar to those of bacteria. Our results indicate that cross-domain organismic factors are important variables with a higher influence on protists as compared with bacteria. Abiotic physical and chemical factors also contributed significantly to community patterns. The contribution of these latter factors was higher for bacteria, which may reflect a stronger biogeochemical coupling. The contribution of geographical distance was similar for both microbial groups.

Thiocyanate Degradation by a Highly Enriched Culture of the Neutrophilic Halophile Thiohalobacter sp. Strain FOKN1 from Activated Sludge and Genomic Insights into Thiocyanate Metabolism

Thiocyanate (SCN^-) is harmful to a wide range of organisms, and its removal is essential for environmental protection. A neutrophilic halophile capable of thiocyanate degradation, Thiohalobacter sp. strain FOKN1, was highly enriched (relative abundance; 98.4%) from activated sludge collected from a bioreactor receiving thiocyanate-rich wastewater. The enrichment culture degraded 3.38 mM thiocyanate within 140 h, with maximum activity at pH 8.8, 37°C, and 0.18 M sodium chloride. Thiocyanate degradation was inhibited by 30 mg L^-1 phenol, but not by thiosulfate. Microbial thiocyanate degradation is catalyzed by thiocyanate dehydrogenase, while limited information is currently available on the molecular mechanisms underlying thiocyanate degradation by the thiocyanate dehydrogenase of neutrophilic halophiles. Therefore, (meta)genomic and proteomic analyses of enrichment cultures were performed to elucidate the whole genome sequence and proteome of Thiohalobacter sp. strain FOKN1. The 3.23-Mb circular Thiohalobacter sp. strain FOKN1 genome was elucidated using a PacBio RSII sequencer, and the expression of 914 proteins was identified by tandem mass spectrometry. The Thiohalobacter sp. strain FOKN1 genome had a gene encoding thiocyanate dehydrogenase, which was abundant in the proteome, suggesting that thiocyanate is degraded by thiocyanate dehydrogenase to sulfur and cyanate. The sulfur formed may be oxidized to sulfate by the sequential oxidation reactions of dissimilatory sulfite reductase, adenosine-5'-phosphosulfate reductase, and dissimilatory ATP sulfurylase. Although the Thiohalobacter sp. strain FOKN1 genome carried a gene encoding cyanate lyase, its protein expression was not detectable. The present study advances the understanding of the molecular mechanisms underlying thiocyanate degradation by the thiocyanate dehydrogenase of neutrophilic halophiles.

Temporal variation of eukaryotic community structures in UASB reactor treating domestic sewage as revealed by 18S rRNA gene sequencing

Eukaryotes are important components of ecosystems in wastewater treatment processes. However, little is known about eukaryotic community in anaerobic wastewater treatment systems. In this study, eukaryotic communities in an up flow anaerobic sludge blanket (UASB) reactor treating domestic sewage during two years of operation were investigated using V4 and V9 regions of 18S rRNA gene for amplicon sequencing. In addition, activated sludge and influent sewage samples were also analyzed and used as the references for aerobic eukaryotic community to characterize anaerobic eukaryotes. The amplicon sequence V4 and V9 libraries detected different taxonomic groups, especially from the UASB samples, suggesting that commonly used V4 and V9 primer pairs could produce a bias for eukaryotic communities analysis. Eukaryotic community structures in the UASB reactor were influenced by the immigration of eukaryotes via influent sewage but were clearly different from the influent sewage and activated sludge. Multivariate statistics indicated that protist genera Cyclidium, Platyophrya and Subulatomonas correlated with chemical oxygen demand and suspended solid concentration, and could be used as bioindicators of treatment performance. Uncultured eukaryotes groups were dominant in the UASB reactor, and their physiological roles need to be examined to understand their contributions to anaerobic processes in future studies.

Synergistic inhibition of anaerobic ammonium oxidation (anammox) activity by phenol and thiocyanate

Coke-oven wastewater discharged from the steel-manufacturing process is phenol and thiocyanate (SCN)-rich wastewater, which inhibits microbial activities in biological wastewater treatment processes. In the present study, synergistic inhibition of anaerobic ammonium oxidation (anammox) activity by phenol and SCN was examined by batch incubation and continuous operation of an anammox reactor. The comparison of anammox activities determined in the batch incubation, in which the anammox biomass was anoxically incubated with 10-250 mg L^-1 of i) phenol, ii) SCN, or iii) both phenol and SCN, showed that synergistic inhibition by phenol and SCN was greater than the inhibitions by phenol or SCN alone. The synergistic inhibition by phenol and SCN was further investigated by operating an up-flow column anammox reactor for 262 d. The removal efficiencies of NH₄⁺ and NO₂^- deteriorated when phenol and SCN concentrations in the influent increased to 16 and 32 mg L^-1, respectively, and the inhibition of anammox activity was further investigated by a¹⁵NO₂^- tracer experiment. Addition of phenol and SCN resulted in a population shift of anammox bacteria, and the dominant species changed from "Candidatus Kuenenia stuttgartiensis" to "Ca. Brocadia sinica". The relative abundance of Azoarcus and Thiobacillus 16S rRNA gene reads increased during the operation, suggesting that they were responsible for the anaerobic phenol and SCN degradation. The present study is the first to document the synergistic inhibition of anammox activity by phenol and SCN and the microbial consortia involved in the nitrogen removal as well as the phenol and SCN degradations.

Ureolytic Prokaryotes in Soil: Community Abundance and Diversity

Although the turnover of urea is a crucial process in nitrogen transformation in soil, limited information is currently available on the abundance and diversity of ureolytic prokaryotes. The abundance and diversity of the soil 16S rRNA gene and ureC (encoding a urease catalytic subunit) were examined in seven soil types using quantitative PCR and amplicon sequencing with Illumina MiSeq. The amplicon sequencing of ureC revealed that the ureolytic community was composed of phylogenetically varied prokaryotes, and we detected 363 to 1,685 species-level ureC operational taxonomic units (OTUs) per soil sample, whereas 5,984 OTUs were site-specific OTUs found in only one of the seven soil types.

Nitrogen Cycle Evaluation (NiCE) Chip for Simultaneous Analysis of Multiple N Cycle-Associated Genes

Various microorganisms play key roles in the nitrogen (N) cycle. Quantitative PCR (qPCR) and PCR amplicon sequencing of N cycle functional genes allow us to analyze the abundance and diversity of microbes responsible for N-transforming reactions in various environmental samples. However, analysis of multiple target genes can be cumbersome and expensive. PCR-independent analysis, such as metagenomics and metatranscriptomics, is useful but expensive, especially when we analyze multiple samples and try to detect N cycle functional genes present at a relatively low abundance. Here, we present the application of microfluidic qPCR chip technology to simultaneously quantify and prepare amplicon sequence libraries for multiple N cycle functional genes as well as taxon-specific 16S rRNA gene markers for many samples. This approach, named the nitrogen cycle evaluation (NiCE) chip, was evaluated by using DNA from pure and artificially mixed bacterial cultures and by comparing the results with those obtained by conventional qPCR and amplicon sequencing methods. Quantitative results obtained by the NiCE chip were comparable to those obtained by conventional qPCR. In addition, the NiCE chip was successfully applied to examine the abundance and diversity of N cycle functional genes in wastewater samples. Although nonspecific amplification was detected on the NiCE chip, this can be overcome by optimizing the primer sequences in the future. As the NiCE chip can provide a high-throughput format to quantify and prepare sequence libraries for multiple N cycle functional genes, this tool should advance our ability to explore N cycling in various samples.IMPORTANCE We report a novel approach, namely, the nitrogen cycle evaluation (NiCE) chip, by using microfluidic qPCR chip technology. By sequencing the amplicons recovered from the NiCE chip, we can assess the diversities of N cycle functional genes. The NiCE chip technology is applicable to analysis of the temporal dynamics of N cycle gene transcription in wastewater treatment bioreactors. The NiCE chip can provide a high-throughput format to quantify and prepare sequence libraries for multiple N cycle functional genes. While there is room for future improvement, this tool should significantly advance our ability to explore the N cycle in various environmental samples.

Nitrogen Cycle Evaluation (NiCE) Chip for Simultaneous Analysis of Multiple N Cycle-Associated Genes

Various microorganisms play key roles in the nitrogen (N) cycle. Quantitative PCR (qPCR) and PCR amplicon sequencing of N cycle functional genes allow us to analyze the abundance and diversity of microbes responsible for N-transforming reactions in various environmental samples. However, analysis of multiple target genes can be cumbersome and expensive. PCR-independent analysis, such as metagenomics and metatranscriptomics, is useful but expensive, especially when we analyze multiple samples and try to detect N cycle functional genes present at a relatively low abundance. Here, we present the application of microfluidic qPCR chip technology to simultaneously quantify and prepare amplicon sequence libraries for multiple N cycle functional genes as well as taxon-specific 16S rRNA gene markers for many samples. This approach, named the nitrogen cycle evaluation (NiCE) chip, was evaluated by using DNA from pure and artificially mixed bacterial cultures and by comparing the results with those obtained by conventional qPCR and amplicon sequencing methods. Quantitative results obtained by the NiCE chip were comparable to those obtained by conventional qPCR. In addition, the NiCE chip was successfully applied to examine the abundance and diversity of N cycle functional genes in wastewater samples. Although nonspecific amplification was detected on the NiCE chip, this can be overcome by optimizing the primer sequences in the future. As the NiCE chip can provide a high-throughput format to quantify and prepare sequence libraries for multiple N cycle functional genes, this tool should advance our ability to explore N cycling in various samples.IMPORTANCE We report a novel approach, namely, the nitrogen cycle evaluation (NiCE) chip, by using microfluidic qPCR chip technology. By sequencing the amplicons recovered from the NiCE chip, we can assess the diversities of N cycle functional genes. The NiCE chip technology is applicable to analysis of the temporal dynamics of N cycle gene transcription in wastewater treatment bioreactors. The NiCE chip can provide a high-throughput format to quantify and prepare sequence libraries for multiple N cycle functional genes. While there is room for future improvement, this tool should significantly advance our ability to explore the N cycle in various environmental samples.

Prokaryote species richness is positively correlated with eukaryote abundance in wastewater treatment biofilms

Biological treatment represents a key step in nutrient removal from wastewater. Until now these process has mainly been considered prokaryotic, with the interactions between prokaryotes and eukaryotes not being properly explored. We therefore investigated the co-occurrence of eukaryotes and prokaryotes in biological nitrogen removal biofilms. We found that biofilms in the nitrifying reactor contained the highest diversity and abundance of both prokaryotes and eukaryotes, with nearly three times higher prokaryote species richness than for the denitrifying reactor. The positive associations between eukaryote abundance and prokaryote diversity could potentially be explained by mutualism - and/or predator/prey interactions. Further mechanistic insight, however, is needed to determine the main diversifying mechanisms. In summary, eukaryote and prokaryote interactions seem to play a fundamental yet underexplored role in biological wastewater treatment.

SIGNIFICANCE AND IMPACT OF THE STUDY

Eukaryote and prokaryote interactions may play an important role in wastewater treatment. This study found that prokaryote species richness was nearly three times higher in the aerobe nitrification than in an anaerobe denitrification reactor, coinciding with the highest level of eukaryotes. This knowledge can be important in process control, and potentially in the development of novel approaches based on nitrate accumulating denitrifying eukaryotes.