Metagenomic and metatranscriptomic analysis of microbial community structure and gene expression of activated sludge

A hundred years of activated sludge: time for a rethink

Biological wastewater treatment plants (BWWTPs) based on the activated sludge (AS) process have dramatically improved worldwide water sanitation despite increased urbanization and industrialization. However, current AS-based operations are considered economically and environmentally unsustainable. In this Perspective, we discuss our current understanding of microbial populations and their metabolic transformations in AS-based BWWTPs in view of developing more sustainable processes in the future. In particular, much has been learned over the course of the past 25 years about specialized microorganisms, which could be more comprehensively leveraged to recover energy and/or nutrients from wastewater streams. To achieve this, we propose a bottom-up design approach, focused around the concept of a "wastewater biorefinery column", which would rely on the engineering of distinct ecological niches into a BWWTP in order to guarantee the targeted enrichment of specific organismal groups which in turn will allow the harvest of high-value resources from wastewater. This concept could be seen as a possible grand challenge to microbial ecologists and engineers alike at the centenary of the discovery of the AS process.

Analysis of microbial diversity by pyrosequencing the small-subunit ribosomal RNA without PCR amplification

To avoid the biases associated with PCR amplification in analysis of microbial communities, a new method has been tested for direct sequencing of the cDNA of full-length Small-subunit Ribosomal RNA withOut specific PCR amplification (SROP). In silico analysis of the SROP method demonstrated that more than 99 % of the SROP sequences could be correctly annotated. Two environmental samples (activated sludge and anaerobic sludge) with complex microbial communities were used for comparison in this study. The SROP results demonstrated that the families Rhodocyclaceae and Nitrosomonadaceae in activated sludge and the phyla Synergistetes and Spirochaetes in anaerobic sludge were underestimated by PCR-based detection. One third of the 16S ribosomal RNA (rRNA) sequences obtained by the SROP method covered the V3 amplicon region, and they are suitable for phylogenetic and diversity index analyses. The microbial diversity index calculated from the rRNA sequences by the SROP was much higher than that calculated by conventional PCR, particularly for the anaerobic sludge. The metatranscriptome-based SROP method will contribute to our better understanding of the diversity of complex microbial communities.

Mining and assessment of catabolic pathways in the metagenome of a common effluent treatment plant to induce the degradative capacity of biomass

Metagenome analysis was used to understand the microbial community in activated sludge treating industrial wastewaters at a Common Effluent Treatment Plant (CETP) in South India. The taxonomic profile mapped onto National Center for Biotechnology Information (NCBI) taxonomy using MEtaGenome ANalyzer (MEGAN), demonstrated that the most abundant domain belonged to prokaryotes, dominated by bacteria. Bacteria representing nine phyla were identified from the sequence data including representatives from two new phyla, Synergistetes and Elusimicrobia. Functional analysis of the metagenome, with specific reference to the metabolism of aromatic compounds, revealed the dominance of genes of the central meta-cleavage pathway. This information was used to improve the degradative efficiency in the wastewater treatment plant. A pilot scale plant was set up with 200L of activated sludge using salicylate induced sludge and results demonstrated 52% removal in chemical oxygen demand (COD) against non-induced biomass.

Interactions in the microbiome: communities of organisms and communities of genes

A central challenge in microbial community ecology is the delineation of appropriate units of biodiversity, which can be taxonomic, phylogenetic, or functional in nature. The term 'community' is applied ambiguously; in some cases, the term refers simply to a set of observed entities, while in other cases, it requires that these entities interact with one another. Microorganisms can rapidly gain and lose genes, potentially decoupling community roles from taxonomic and phylogenetic groupings. Trait-based approaches offer a useful alternative, but many traits can be defined based on gene functions, metabolic modules, and genomic properties, and the optimal set of traits to choose is often not obvious. An analysis that considers taxon assignment and traits in concert may be ideal, with the strengths of each approach offsetting the weaknesses of the other. Individual genes also merit consideration as entities in an ecological analysis, with characteristics such as diversity, turnover, and interactions modeled using genes rather than organisms as entities. We identify some promising avenues of research that are likely to yield a deeper understanding of microbial communities that shift from observation-based questions of 'Who is there?' and 'What are they doing?' to the mechanistically driven question of 'How will they respond?'

Metagenomic approaches for exploiting uncultivated bacteria as a resource for novel biosynthetic enzymology

Most biologically active microbial natural products are known from strains that can be isolated and cultivated in the laboratory. However, the genomics era has revealed that cultured bacteria represent a mere fraction of total estimated bacterial biodiversity. With the development of community genomics, termed metagenomics, the uncultivated majority became accessible for functional analysis. Through metagenomic studies, novel biocatalysts and biosynthetic pathways are being discovered at a pace previously not possible using traditional molecular biology techniques. Additionally, the study of uncultivated bacteria has provided valuable insights into previously overlooked biocatalysts from cultured strains. This perspective highlights recent discoveries from metagenomics of uncultivated bacteria and discusses the impact of those findings on the field of natural products.

Comparative analysis of the enzyme activities and the bacterial community structure based on the aeration source supplied to an MBR to treat urban wastewater

A comparative analysis was performed in a pilot-scale membrane bioreactor (MBR) treating urban wastewater supplied with either pure oxygen (O2) or air, to assess the influence of each aeration source on the diversity and activity of the bacterial communities in the sludge. The MBR was operated in three experimental stages with different concentrations of volatile suspended solids (VSS) and temperature, and under both aeration conditions. α-Glucosidases, proteases, esterases and phosphatases were tested as markers of organic matter removal in the sludge, and the diversity of the bacterial community was analysed by fingerprinting (temperature-gradient gel electrophoresis of partially-amplified 16S-rRNA genes). Redundancy analysis (RDA) revealed that temperature and VSS concentration were the only factors that significantly influenced the levels of enzyme activities and the values of both the Shannon-Wiener diversity index (H') and the functional organisation index (Fo), while the bacterial community structure experienced significant changes depending on the aeration source supplied in each experimental stage.

Efficiency of RNA extraction from selected bacteria in the context of biogas production and metatranscriptomics

Understanding the microbial population in anaerobic digestion is an essential task to increase efficient substrate use and process stability. The metabolic state, represented e.g. by the transcriptome, of a fermenting system can help to find markers for monitoring industrial biogas production to prevent failures or to model the whole process. Advances in next-generation sequencing make transcriptomes accessible for large-scale analyses. In order to analyze the metatranscriptome of a mixed-species sample, isolation of high-quality RNA is the first step. However, different extraction methods may yield different efficiencies in different species. Especially in mixed-species environmental samples, unbiased isolation of transcripts is important for meaningful conclusions. We applied five different RNA-extraction protocols to nine taxonomic diverse bacterial species. Chosen methods are based on various lysis and extraction principles. We found that the extraction efficiency of different methods depends strongly on the target organism. RNA isolation of gram-positive bacteria was characterized by low yield whilst from gram-negative species higher concentrations can be obtained. Transferring our results to mixed-species investigations, such as metatranscriptomics with biofilms or biogas plants, leads to the conclusion that particular microorganisms might be over- or underrepresented depending on the method applied. Special care must be taken when using such metatranscriptomics data for, e.g. process modeling.

Metagenomic profiling of antibiotic resistance genes and mobile genetic elements in a tannery wastewater treatment plant

Antibiotics are often used to prevent sickness and improve production in animal agriculture, and the residues in animal bodies may enter tannery wastewater during leather production. This study aimed to use Illumina high-throughput sequencing to investigate the occurrence, diversity and abundance of antibiotic resistance genes (ARGs) and mobile genetic elements (MGEs) in aerobic and anaerobic sludge of a full-scale tannery wastewater treatment plant (WWTP). Metagenomic analysis showed that Proteobacteria, Firmicutes, Bacteroidetes and Actinobacteria dominated in the WWTP, but the relative abundance of archaea in anaerobic sludge was higher than in aerobic sludge. Sequencing reads from aerobic and anaerobic sludge revealed differences in the abundance of functional genes between both microbial communities. Genes coding for antibiotic resistance were identified in both communities. BLAST analysis against Antibiotic Resistance Genes Database (ARDB) further revealed that aerobic and anaerobic sludge contained various ARGs with high abundance, among which sulfonamide resistance gene sul1 had the highest abundance, occupying over 20% of the total ARGs reads. Tetracycline resistance genes (tet) were highly rich in the anaerobic sludge, among which tet33 had the highest abundance, but was absent in aerobic sludge. Over 70 types of insertion sequences were detected in each sludge sample, and class 1 integrase genes were prevalent in the WWTP. The results highlighted prevalence of ARGs and MGEs in tannery WWTPs, which may deserve more public health concerns.

Mining RNA-seq data for infections and contaminations

RNA sequencing (RNA-seq) provides novel opportunities for transcriptomic studies at nucleotide resolution, including transcriptomics of viruses or microbes infecting a cell. However, standard approaches for mapping the resulting sequencing reads generally ignore alternative sources of expression other than the host cell and are little equipped to address the problems arising from redundancies and gaps among sequenced microbe and virus genomes. We show that screening of sequencing reads for contaminations and infections can be performed easily using ContextMap, our recently developed mapping software. Based on mapping-derived statistics, mapping confidence, similarities and misidentifications (e.g. due to missing genome sequences) of species/strains can be assessed. Performance of our approach is evaluated on three real-life sequencing data sets and compared to state-of-the-art metagenomics tools. In particular, ContextMap vastly outperformed GASiC and GRAMMy in terms of runtime. In contrast to MEGAN4, it was capable of providing individual read mappings to species and resolving non-unique mappings, thus allowing the identification of misalignments caused by sequence similarities between genomes and missing genome sequences. Our study illustrates the importance and potentials of routinely mining RNA-seq experiments for infections or contaminations by microbes and viruses. By using ContextMap, gene expression of infecting agents can be analyzed and novel insights in infection processes and tumorigenesis can be obtained.

Using high-throughput sequencing to assess the impacts of treated and untreated wastewater discharge on prokaryotic communities in an urban river

In many megacities wastewater is an important source of surface water, particularly during drought periods. While changes in surface water chemistry associated with effluent inflow have generally been well-studied, few data have been collected on the effects to prokaryotic communities. The objective of this study was to explore the impacts of treated and untreated wastewater discharges on prokaryotic community in an urban river. High-throughput sequencing was conducted for analyzing the prokaryotic community composition and function in river water, treated wastewater and untreated wastewater. Results revealed that the prokaryotic community compositions in the upstream river reach were dominated by treated wastewater discharge. In the middle- and downstream river reaches, untreated effluent volumes are higher, thus affecting the structure of the prokaryotic community, promoting a rise in Cyanobacteria and Thaumarchaeota. Function annotation revealed a number of genes associated with xenobiotic metabolism and human diseases were observed in river and wastewater samples, suggesting wastewater discharge to river may pose a risk to human health. Quantitative real-time PCR results revealed that the treated and untreated wastewater discharges also affected the abundance of ammonia oxidation bacteria (AOB) and ammonia oxidation archaea (AOA) in river.

Industrial activated sludge exhibit unique bacterial community composition at high taxonomic ranks

Biological degradation of domestic and industrial wastewater by activated sludge depends on a common process of separation of the diverse self-assembled and self-sustained microbial flocs from the treated wastewater. Previous surveys of bacterial communities indicated the presence of a common core of bacterial phyla in municipal activated sludge, an observation consistent with the concept of ecological coherence of high taxonomic ranks. The aim of this work was to test whether this critical feature brings about a common pattern of abundance distribution of high bacterial taxa in industrial and domestic activated sludge, and to relate the bacterial community structure of industrial activated sludge with relevant operational parameters. We have applied 454 pyrosequencing of 16S rRNA genes to evaluate bacterial communities in full-scale biological wastewater treatment plants sampled at different times, including seven systems treating wastewater from different industries and one plant that treats domestic wastewater, and compared our datasets with the data from municipal wastewater treatment plants obtained by three different laboratories. We observed that each industrial activated sludge system exhibited a unique bacterial community composition, which is clearly distinct from the common profile of bacterial phyla or classes observed in municipal plants. The influence of process parameters on the bacterial community structure was evaluated using constrained analysis of principal coordinates (CAP). Part of the differences in the bacterial community structure between industrial wastewater treatment systems were explained by dissolved oxygen and pH. Despite the ecological relevance of floc formation for the assembly of bacterial communities in activated sludge, the wastewater characteristics are likely to be the major determinant that drives bacterial composition at high taxonomic ranks.

Abundance and diversity of ammonia-oxidizing archaea in response to various habitats in Pearl River Delta of China, a subtropical maritime zone

Ammonia-oxidizing archaea (AOA) are widely considered key to ammonia oxidation in various environments. However, little work has been conducted to simultaneously investigate the abundance and diversity of AOA as well as correlations between archaeal amoA genotypes and environmental parameters of different ecosystems at one district. To understand the abundance, diversity, and distribution of AOA in Pearl River Delta of China in response to various habitats, the archaeal amoA genes in soil, marine, river, lake, hot spring and wastewater treatment plant (WWTP) samples were investigated using real-time fluorescent quantitative PCR and clone libraries. Our analyses indicated that the diversity of AOA in various habitats was different and could be clustered into five major clades, i.e., estuary sediment, marine water/sediment, soil, hot spring and Cluster 1. Phylogenetic analyses revealed that the structure of AOA communities in similar ecological habitats exhibited strong relation. The canonical correspondence method indicated that the AOA community structure was strongly correlated to temperature, pH, total organic carbon, total nitrogen and dissolved oxygen variables. Assessing AOA amoA gene copy numbers, ranging from 6.84 x 10(6) to 9.45 x 10(7) copies/g in dry soil/sediment, and 6.06 x 10(6) to 2.41 x 10(7) copies/L in water samples, were higher than ammonia-oxidizing bacteria (AOB) by 1-2 orders of magnitude. However, AOA amoA copy numbers were much lower than AOB in WWTP activated sludge samples. Overall, these studies suggested that AOA may be a major contributor to ammonia oxidation in natural habitats but play a minor role in highly aerated activated sludge. The result also showed the ratio of AOA to AOB amoA gene abundance was positively correlated with temperature and less correlated with other environmental parameters. New data from our study provide increasing evidence for the relative abundance and diversity of ammonia-oxidizing archaea in the global nitrogen cycle.

Human health implications of clinically relevant bacteria in wastewater habitats

The objective of this review is to reflect on the multiple roles of bacteria in wastewater habitats with particular emphasis on their harmful potential for human health. Indigenous bacteria promote a series of biochemical and metabolic transformations indispensable to achieve wastewater treatment. Some of these bacteria may be pathogenic or harbour antibiotic resistance or virulence genes harmful for human health. Several chemical contaminants (heavy metals, disinfectants and antibiotics) may select these bacteria or their genes. Worldwide studies show that treated wastewater contain antibiotic resistant bacteria or genes encoding virulence or antimicrobial resistance, evidencing that treatment processes may fail to remove efficiently these bio-pollutants. The contamination of the surrounding environment, such as rivers or lakes receiving such effluents, is also documented in several studies. The current state of the art suggests that only some of antibiotic resistance and virulence potential in wastewater is known. Moreover, wastewater habitats may favour the evolution and dissemination of new resistance and virulence genes and the emergence of new pathogens. For these reasons, additional research is needed in order to obtain a more detailed assessment of the long-term effects of wastewater discharges. In particular, it is important to measure the human and environmental health risks associated with wastewater reuse.

Assessing bacterial diversity in a seawater-processing wastewater treatment plant by 454-pyrosequencing of the 16S rRNA and amoA genes

The bacterial community composition of activated sludge from a wastewater treatment plant (Almería, Spain) with the particularity of using seawater was investigated by applying 454-pyrosequencing. The results showed that Deinococcus-Thermus, Proteobacteria, Chloroflexi and Bacteroidetes were the most abundant retrieved sequences, while other groups, such as Actinobacteria, Chlorobi, Deferribacteres, Firmicutes, Planctomycetes, Spirochaetes and Verrumicrobia were reported at lower proportions. Rarefaction analysis showed that very likely the diversity is higher than what could be described despite most of the unknown microorganisms probably correspond to rare diversity. Furthermore, the majority of taxa could not be classified at the genus level and likely represent novel members of these groups. Additionally, the nitrifiers in the sludge were characterized by pyrosequencing the amoA gene. In contrast, the nitrifying bacterial community, dominated by the genera Nitrosomonas, showed a low diversity and rarefaction curves exhibited saturation. These results suggest that only a few populations of low abundant but specialized bacteria are responsible for removal of ammonia in these saline wastewater systems.

Construction of customized sub-databases from NCBI-nr database for rapid annotation of huge metagenomic datasets using a combined BLAST and MEGAN approach

We developed a fast method to construct local sub-databases from the NCBI-nr database for the quick similarity search and annotation of huge metagenomic datasets based on BLAST-MEGAN approach. A three-step sub-database annotation pipeline (SAP) was further proposed to conduct the annotation in a much more time-efficient way which required far less computational capacity than the direct NCBI-nr database BLAST-MEGAN approach. The 1^st BLAST of SAP was conducted using the original metagenomic dataset against the constructed sub-database for a quick screening of candidate target sequences. Then, the candidate target sequences identified in the 1^st BLAST were subjected to the 2^nd BLAST against the whole NCBI-nr database. The BLAST results were finally annotated using MEGAN to filter out those mistakenly selected sequences in the 1^st BLAST to guarantee the accuracy of the results. Based on the tests conducted in this study, SAP achieved a speedup of ∼150–385 times at the BLAST e-value of 1e–5, compared to the direct BLAST against NCBI-nr database. The annotation results of SAP are exactly in agreement with those of the direct NCBI-nr database BLAST-MEGAN approach, which is very time-consuming and computationally intensive. Selecting rigorous thresholds (e.g. e-value of 1e–10) would further accelerate SAP process. The SAP pipeline may also be coupled with novel similarity search tools (e.g. RAPsearch) other than BLAST to achieve even faster annotation of huge metagenomic datasets. Above all, this sub-database construction method and SAP pipeline provides a new time-efficient and convenient annotation similarity search strategy for laboratories without access to high performance computing facilities. SAP also offers a solution to high performance computing facilities for the processing of more similarity search tasks.

Diversity and isolation of rare actinomycetes: an overview

A renewed interest in the development of new antimicrobial agents is urgently needed to combat the increasing number of antibiotic-resistant strains of pathogenic microorganisms. Actinomycetes continue to be the mainstream supplier of antibiotics used in industry. The likelihood of discovering a new compound with novel chemical structure can be increased with intensive efforts in isolating and screening of rare genera of microorganisms to include in natural-product-screening collections. An unexpected variety of rare actinomycetes is now being isolated worldwide from previously uninvestigated diverse natural habitats, using different selective isolation methods. These isolation efforts include methods to enhance growth (enrichment) of rare actinomycetes, and eliminate unwanted microorganisms (pretreatment). To speed up the strain isolation process, knowledge about the distribution of such unexploited groups of microorganisms must also be augmented. This is a summary of using these microorganisms as new potential biological resources, and a review of almost all of the selective isolation methods, including pretreatment and enrichment techniques that have been developed to date for the isolation of rare actinomycetes.