Statistical Approach of Functional Profiling for a Microbial Community

Background

Metagenomics is a relatively new but fast growing field within environmental biology and medical sciences. It enables researchers to understand the diversity of microbes, their functions, cooperation, and evolution in a particular ecosystem. Traditional methods in genomics and microbiology are not efficient in capturing the structure of the microbial community in an environment. Nowadays, high-throughput next-generation sequencing technologies are powerfully driving the metagenomic studies. However, there is an urgent need to develop efficient statistical methods and computational algorithms to rapidly analyze the massive metagenomic short sequencing data and to accurately detect the features/functions present in the microbial community. Although several issues about functions of metagenomes at pathways or subsystems level have been investigated, there is a lack of studies focusing on functional analysis at a low level of a hierarchical functional tree, such as SEED subsystem tree.

Results

A two-step statistical procedure (metaFunction) is proposed to detect all possible functional roles at the low level from a metagenomic sample/community. In the first step a statistical mixture model is proposed at the base of gene codons to estimate the abundances for the candidate functional roles, with sequencing error being considered. As a gene could be involved in multiple biological processes the functional assignment is therefore adjusted by utilizing an error distribution in the second step. The performance of the proposed procedure is evaluated through comprehensive simulation studies. Compared with other existing methods in metagenomic functional analysis the new approach is more accurate in assigning reads to functional roles, and therefore at more general levels. The method is also employed to analyze two real data sets.

Conclusions

metaFunction is a powerful tool in accurate profiling functions in a metagenomic sample.     

Microbial Community of a Hydrothermal Mud Vent Underneath the Deep-Sea Anoxic Brine Lake Urania (Eastern Mediterranean)

The composition of a metabolically active prokaryotic community thriving in hydrothermal mud fluids of the deep-sea hypersaline anoxic Western Urania Basin was characterized using rRNA-based phylogenetic analysis of a clone library. The physiologically active prokaryotic assemblage in this extreme environment showed a great genetic diversity. Most members of the microbial community appeared to be affiliated to yet uncultured organisms from similar ecosystems, i.e., deep-sea hypersaline basins and hydrothermal vents. The bacterial clone library was dominated by phylotypes affiliated with the epsilon-Proteobacteria subdivision recognized as an ecologically significant group of bacteria inhabiting deep-sea hydrothermal environments. Almost 18% of all bacterial clones were related to delta-Proteobacteria, suggesting that sulfate reduction is one of the dominant metabolic processes occurring in warm mud fluids. The remaining bacterial phylotypes were related to alpha- and beta-Proteobacteria, Actinobacteria, Bacteroides, Deinococcus-Thermus, KB1 and OP-11 candidate divisions. Moreover, a novel monophyletic clade, deeply branched with unaffiliated 16S rDNA clones was also retrieved from deep-sea sediments and halocline of Urania Basin. Archaeal diversity was much lower and detected phylotypes included organisms affiliated exclusively with the Euryarchaeota. More than 96% of the archaeal clones belonged to the MSBL-1 candidate order recently found in hypersaline anoxic environments, such as endoevaporitic microbial mats, Mediterranean deep-sea mud volcanoes and anoxic basins. Two phylotypes, represented by single clones were related to uncultured groups DHVE-1 and ANME-1. Thus, the hydrothermal mud of hypersaline Urania Basin seems to contain new microbial diversity. The prokaryotic community was significantly different from that occurring in the upper layers of the Urania Basin since 60% of all bacterial and 40% of all archaeal phylotypes were obtained only from mud fluids. The uniqueness of the composition of the active prokaryotic community could be explained by the complex environmental conditions at the site. The interaction of oxygenated warm mud fluids with the cold hypersaline brine of the Urania Basin seems to simultaneously select for various metabolic processes, such as aerobic and anaerobic heterotrophy, sulfide- and methane-dependent chemotrophy along with anaerobic oxidation of methane, sulfate- and metal-reduction.     

Taxonomic Profiling and Metagenome Analysis of a Microbial Community from a Habitat Contaminated with Industrial Discharges

Industrial units, manufacturing dyes, chemicals, solvents, and xenobiotic compounds, produce liquid and solid wastes, which upon conventional treatment are released in the nearby environment and thus are the major cause of pollution. Soil collected from contaminated Kharicut Canal bank (N 22°57.878′; E 072°38.478′), Ahmedabad, Gujarat, India was used for metagenomic DNA preparation to study the capabilities of intrinsic microbial community in dealing with xenobiotics. Sequencing of metagenomic DNA on the Genome Sequencer FLX System using titanium chemistry resulted in 409,782 reads accounting for 133,529,997 bases of sequence information. Taxonomic analyses and gene annotations were carried out using the bioinformatics platform Sequence Analysis and Management System for Metagenomic Datasets. Taxonomic profiling was carried out by three different complementary approaches: (a) 16S rDNA, (b) environmental gene tags, and (c) lowest common ancestor. The most abundant phylum and genus were found to be “Proteobacteria” and “Pseudomonas,” respectively. Metagenome reads were mapped on sequenced microbial genomes and the highest numbers of reads were allocated to Pseudomonas stutzeri A1501. Assignment of obtained metagenome reads to Gene Ontology terms, Clusters of Orthologous Groups of protein categories, protein family numbers, and Kyoto Encyclopedia of Genes and Genomes hits revealed genomic potential of indigenous microbial community. In total, 157,024 reads corresponded to 37,028 different KEGG hits, and amongst them, 11,574 reads corresponded to 131 different enzymes potentially involved in xenobiotic biodegradation. These enzymes were mapped on biodegradation pathways of xenobiotics to elucidate their roles in possible catalytic reactions. Consequently, information obtained from the present study will act as a baseline which, subsequently along with other “-omic” studies, will help in designing future bioremediation strategies in effluent treatment plants and environmental clean-up projects.     

Evidence for Microbial Fe(III) Reduction in Anoxic, Mining-Impacted Lake Sediments (Lake Coeur d'Alene, Idaho)

Mining-impacted sediments of Lake Coeur d'Alene, Idaho, contain more than 10% metals on a dry weight basis, approximately 80% of which is iron. Since iron (hydr)oxides adsorb toxic, ore-associated elements, such as arsenic, iron (hydr)oxide reduction may in part control the mobility and bioavailability of these elements. Geochemical and microbiological data were collected to examine the ecological role of dissimilatory Fe(III)-reducing bacteria in this habitat. The concentration of mild-acid-extractable Fe(II) increased with sediment depth up to 50 g kg⁻¹, suggesting that iron reduction has occurred recently. The maximum concentrations of dissolved Fe(II) in interstitial water (41 mg liter⁻¹) occurred 10 to 15 cm beneath the sediment-water interface, suggesting that sulfidogenesis may not be the predominant terminal electron-accepting process in this environment and that dissolved Fe(II) arises from biological reductive dissolution of iron (hydr)oxides. The concentration of sedimentary magnetite (Fe₃O₄), a common product of bacterial Fe(III) hydroxide reduction, was as much as 15.5 g kg⁻¹. Most-probable-number enrichment cultures revealed that the mean density of Fe(III)-reducing bacteria was 8.3 × 10⁵ cells g (dry weight) of sediment⁻¹. Two new strains of dissimilatory Fe(III)-reducing bacteria were isolated from surface sediments. Collectively, the results of this study support the hypothesis that dissimilatory reduction of iron has been and continues to be an important biogeochemical process in the environment examined.     

Sediment Microbial Community Structure and Mercury Methylation in Mercury-Polluted Clear Lake, California

Spatial and temporal variations in sediment microbial community structure in a eutrophic lake polluted with inorganic mercury were identified using polar lipid fatty acid (PLFA) analysis. Microbial community structure was strongly related to mercury methylation potential, sediment organic carbon content, and lake location. Pore water sulfate, total mercury concentrations, and organic matter C/N ratios showed no relationships with microbial community structure. Seasonal changes and changes potentially attributable to temperature regulation of bacterial membranes were detectable but were less important influences on sediment PLFA composition than were differences due to lake sampling location. Analysis of biomarker PLFAs characteristic of Desulfobacter and Desulfovibrio groups of sulfate-reducing bacteria suggests that Desulfobacter-like organisms are important mercury methylators in the sediments, especially in the Lower Arm of Clear Lake.     

Defining the Functional Potential and Active Community Members of a Sediment Microbial Community in a High-Arctic Hypersaline Subzero Spring

The Lost Hammer (LH) Spring is the coldest and saltiest terrestrial spring discovered to date and is characterized by perennial discharges at subzero temperatures (−5°C), hypersalinity (salinity, 24%), and reducing (≈−165 mV), microoxic, and oligotrophic conditions. It is rich in sulfates (10.0%, wt/wt), dissolved H₂S/sulfides (up to 25 ppm), ammonia (≈381 μM), and methane (11.1 g day⁻¹). To determine its total functional and genetic potential and to identify its active microbial components, we performed metagenomic analyses of the LH Spring outlet microbial community and pyrosequencing analyses of the cDNA of its 16S rRNA genes. Reads related to Cyanobacteria (19.7%), Bacteroidetes (13.3%), and Proteobacteria (6.6%) represented the dominant phyla identified among the classified sequences. Reconstruction of the enzyme pathways responsible for bacterial nitrification/denitrification/ammonification and sulfate reduction appeared nearly complete in the metagenomic data set. In the cDNA profile of the LH Spring active community, ammonia oxidizers (Thaumarchaeota), denitrifiers (Pseudomonas spp.), sulfate reducers (Desulfobulbus spp.), and other sulfur oxidizers (Thermoprotei) were present, highlighting their involvement in nitrogen and sulfur cycling. Stress response genes for adapting to cold, osmotic stress, and oxidative stress were also abundant in the metagenome. Comparison of the composition of the functional community of the LH Spring to metagenomes from other saline/subzero environments revealed a close association between the LH Spring and another Canadian high-Arctic permafrost environment, particularly in genes related to sulfur metabolism and dormancy. Overall, this study provides insights into the metabolic potential and the active microbial populations that exist in this hypersaline cryoenvironment and contributes to our understanding of microbial ecology in extreme environments.     

Metagenomic Insights into the Evolution, Function, and Complexity of the Planktonic Microbial Community of Lake Lanier, a Temperate Freshwater Ecosystem

Lake Lanier is an important freshwater lake for the southeast United States, as it represents the main source of drinking water for the Atlanta metropolitan area and is popular for recreational activities. Temperate freshwater lakes such as Lake Lanier are underrepresented among the growing number of environmental metagenomic data sets, and little is known about how functional gene content in freshwater communities relates to that of other ecosystems. To better characterize the gene content and variability of this freshwater planktonic microbial community, we sequenced several samples obtained around a strong summer storm event and during the fall water mixing using a random whole-genome shotgun (WGS) approach. Comparative metagenomics revealed that the gene content was relatively stable over time and more related to that of another freshwater lake and the surface ocean than to soil. However, the phylogenetic diversity of Lake Lanier communities was distinct from that of soil and marine communities. We identified several important genomic adaptations that account for these findings, such as the use of potassium (as opposed to sodium) osmoregulators by freshwater organisms and differences in the community average genome size. We show that the lake community is predominantly composed of sequence-discrete populations and describe a simple method to assess community complexity based on population richness and evenness and to determine the sequencing effort required to cover diversity in a sample. This study provides the first comprehensive analysis of the genetic diversity and metabolic potential of a temperate planktonic freshwater community and advances approaches for comparative metagenomics.     

Biochemical Diversity of Carboxyl Esterases and Lipases from Lake Arreo (Spain): a Metagenomic Approach

The esterases and lipases from the α/β hydrolase superfamily exhibit an enormous sequence diversity, fold plasticity, and activities. Here, we present the comprehensive sequence and biochemical analyses of seven distinct esterases and lipases from the metagenome of Lake Arreo, an evaporite karstic lake in Spain (42°46′N, 2°59′W; altitude, 655 m). Together with oligonucleotide usage patterns and BLASTP analysis, our study of esterases/lipases mined from Lake Arreo suggests that its sediment contains moderately halophilic and cold-adapted proteobacteria containing DNA fragments of distantly related plasmids or chromosomal genomic islands of plasmid and phage origins. This metagenome encodes esterases/lipases with broad substrate profiles (tested over a set of 101 structurally diverse esters) and habitat-specific characteristics, as they exhibit maximal activity at alkaline pH (8.0 to 8.5) and temperature of 16 to 40°C, and they are stimulated (1.5 to 2.2 times) by chloride ions (0.1 to 1.2 M), reflecting an adaptation to environmental conditions. Our work provides further insights into the potential significance of the Lake Arreo esterases/lipases for biotechnology processes (i.e., production of enantiomers and sugar esters), because these enzymes are salt tolerant and are active at low temperatures and against a broad range of substrates. As an example, the ability of a single protein to hydrolyze triacylglycerols, (non)halogenated alkyl and aryl esters, cinnamoyl and carbohydrate esters, lactones, and chiral epoxides to a similar extent was demonstrated.     

Microbial Community Structure and Dynamics in the Largest Natural French Lake (Lake Bourget)

We investigated the dynamics and diversity of heterotrophic bacteria, autotrophic and heterotrophic flagellates, and ciliates from March to July 2002 in the surface waters (0–50 m) of Lake Bourget. The heterotrophic bacteria consisted mainly of “small” cocci, but filaments (>2 μm), commonly considered to be grazing-resistant forms under increased nanoflagellate grazing, were also detected. These elongated cells mainly belonged to the Cytophaga-Flavobacterium (CF) cluster, and were most abundant during spring and early summer, when mixotrophic or heterotrophic flagellates were the main bacterial predators. The CF group strongly dominated fluorescent in situ hybridization–detected cells from March to June, whereas clear changes were observed in early summer when Beta-proteobacteria and Alpha-proteobacteria increased concomitantly with maximal protist grazing pressures. The analysis of protist community structure revealed that the flagellates consisted mainly of cryptomonad forms. The dynamics of Cryptomonas sp. and Dinobryon sp. suggested the potential importance of mixotrophs as consumers of bacteria. This point was verified by an experimental approach based on fluorescent microbeads to assess the potential grazing impact of all protist taxa in the epilimnion. From the results, three distinct periods in the functioning of the epilimnetic microbial loop were identified. In early spring, mixotrophic and heterotrophic flagellates constituted the main bacterivores, and were regulated by the availability of their resources mainly during April (phase 1). Once the “clear water phase” was established, the predation pressure of metazooplankton represented a strong top-down force on all microbial compartments. During this period only mixotrophic flagellates occasionally exerted a significant bacterivory pressure (phase 2). Finally, the early summer was characterized by the highest protozoan grazing impact and by a rapid shift in the carbon pathway transfer, with a fast change-over of the main predators contribution, i.e., mixotrophic, heterotrophic flagellates and ciliates in bacterial mortality. The high abundance of ciliates during this period was consistent with the high densities of resources (heterotrophic nanoflagellates, algae, bacteria) in deep layers containing the most chlorophyll. Bacteria, as ciliates, responded clearly to increasing phytoplankton abundance, and although bacterial grazing impact could vary largely, bacterial abundance seemed to be primarily bottom-up regulated (phase 3).     

Functional Metagenomic Investigations of Microbial Communities in a Shallow-Sea Hydrothermal System

Little is known about the functional capability of microbial communities in shallow-sea hydrothermal systems (water depth of <200 m). This study analyzed two high-throughput pyrosequencing metagenomic datasets from the vent and the surface water in the shallow-sea hydrothermal system offshore NE Taiwan. This system exhibited distinct geochemical parameters. Metagenomic data revealed that the vent and the surface water were predominated by Epsilonproteobacteria (Nautiliales-like organisms) and Gammaproteobacteria ( Thiomicrospira -like organisms), respectively. A significant difference in microbial carbon fixation and sulfur metabolism was found between the vent and the surface water. The chemoautotrophic microorganisms in the vent and in the surface water might possess the reverse tricarboxylic acid cycle and the Calvin−Bassham−Benson cycle for carbon fixation in response to carbon dioxide highly enriched in the environment, which is possibly fueled by geochemical energy with sulfur and hydrogen. Comparative analyses of metagenomes showed that the shallow-sea metagenomes contained some genes similar to those present in other extreme environments. This study may serve as a basis for deeply understanding the genetic network and functional capability of the microbial members of shallow-sea hydrothermal systems.     

Isolation and Characterization of Organic Phosphorus-Mineralizing Bacteria in Sediment of a Chinese Large Shallow Eutrophic Lake (Lake Taihu)

The organic phosphorus mineralizing bacteria (OPB) play an important role in phosphorus cycling in lake sediment, to which less attention has been paid. Monthly sediment samples in 2009 ending in October, together with the samples from different seasons (May, July, October, and December) in 2008, were collected at 6 sites in a Chinese large shallow eutrophic lake (Lake Taihu). The sediment OPB numbers ranged from 2.2 × 10⁶ to 1.79 × 10⁸ cells g^?1 (dry weight) at different sampling sites and in different seasons, with the average being 3.88 × 10⁷ cells g^?1 (dry weight). Its number was highest at the most polluted site and peaked in spring and summer, which can be explained by the enrichment of organic matter in sediment. Furthermore, there existed a significant positive relationship between the OPB numbers and alkaline phosphatase activities in the sediment. The 6 OPB strains isolated from the sediment were distinct in terms of their colony morphology on the yolk agar, biochemical characteristics and phosphorus release abilities. According to the 16S rDNA sequences, these OPB belong to Bacillus cereus, Stenotrophomonas maltophilia, Stenotrophomonas sp., Bacillus cereus, Xanthomonas sp., Pseudomonas sp. They were distinguished from the OPB species recorded in a Chinese small shallow eutrophic lake whose sediment organic content was remarkably higher. Taken together, organic matter shaped the OPB community not only quantitatively but also qualitatively, which in concert facilitated the enzymatic hydrolysis of organic phosphorus in lake sediment.     

Metagenomic Profiling of Microbial Composition and Antibiotic Resistance Determinants in Puget Sound

Human-health relevant impacts on marine ecosystems are increasing on both spatial and temporal scales. Traditional indicators for environmental health monitoring and microbial risk assessment have relied primarily on single species analyses and have provided only limited spatial and temporal information. More high-throughput, broad-scale approaches to evaluate these impacts are therefore needed to provide a platform for informing public health. This study uses shotgun metagenomics to survey the taxonomic composition and antibiotic resistance determinant content of surface water bacterial communities in the Puget Sound estuary. Metagenomic DNA was collected at six sites in Puget Sound in addition to one wastewater treatment plant (WWTP) that discharges into the Sound and pyrosequenced. A total of ∼550 Mbp (1.4 million reads) were obtained, 22 Mbp of which could be assembled into contigs. While the taxonomic and resistance determinant profiles across the open Sound samples were similar, unique signatures were identified when comparing these profiles across the open Sound, a nearshore marina and WWTP effluent. The open Sound was dominated by α-Proteobacteria (in particular Rhodobacterales sp.), γ-Proteobacteria and Bacteroidetes while the marina and effluent had increased abundances of Actinobacteria, β-Proteobacteria and Firmicutes. There was a significant increase in the antibiotic resistance gene signal from the open Sound to marina to WWTP effluent, suggestive of a potential link to human impacts. Mobile genetic elements associated with environmental and pathogenic bacteria were also differentially abundant across the samples. This study is the first comparative metagenomic survey of Puget Sound and provides baseline data for further assessments of community composition and antibiotic resistance determinants in the environment using next generation sequencing technologies. In addition, these genomic signals of potential human impact can be used to guide initial public health monitoring as well as more targeted and functionally-based investigations.     

The Influence of Macrophytes on Sediment Resuspension and the Effect of Associated Nutrients in a Shallow and Large Lake (Lake Taihu,China)

A yearlong campaign to examine sediment resuspension was conducted in large, shallow and eutrophic Lake Taihu, China, to investigate the influence of vegetation on sediment resuspension and its nutrient effects. The study was conducted at 6 sites located in both phytoplankton-dominated zone and macrophyte-dominated zone of the lake, lasting for a total of 13 months, with collections made at two-week intervals. Sediment resuspension in Taihu, with a two-week high average rate of 1771 g·m^-2·d^-1 and a yearly average rate of 377 g·m^-2·d^-1, is much stronger than in many other lakes worldwide, as Taihu is quite shallow and contains a long fetch. The occurrence of macrophytes, however, provided quite strong abatement of sediment resuspension, which may reduce the sediment resuspension rate up to 29-fold. The contribution of nitrogen and phosphorus to the water column from sediment resuspension was estimated as 0.34 mg·L^-1 and 0.051 mg·L^-1 in the phytoplankton-dominated zone. Sediment resuspension also largely reduced transparency and then stimulated phytoplankton growth. Therefore, sediment resuspension may be one of the most important factors delaying the recovery of eutrophic Lake Taihu, and the influence of sediment resuspension on water quality must also be taken into account by the lake managers when they determine the restoration target.     

Insights into Biodegradation Through Depth-Resolved Microbial Community Functional and Structural Profiling of a Crude-Oil Contaminant Plume

Small-scale geochemical gradients are a key feature of aquifer contaminant plumes, highlighting the need for functional and structural profiling of corresponding microbial communities on a similar scale. The purpose of this study was to characterize the microbial functional and structural diversity with depth across representative redox zones of a hydrocarbon plume and an adjacent wetland, at the Bemidji Oil Spill site. A combination of quantitative PCR, denaturing gradient gel electrophoresis, and pyrosequencing were applied to vertically sampled sediment cores. Levels of the methanogenic marker gene, methyl coenzyme-M reductase A (mcrA), increased with depth near the oil body center, but were variable with depth further downgradient. Benzoate degradation N (bzdN) hydrocarbon-degradation gene, common to facultatively anaerobic Azoarcus spp., was found at all locations, but was highest near the oil body center. Microbial community structural differences were observed across sediment cores, and bacterial classes containing known hydrocarbon degraders were found to be low in relative abundance. Depth-resolved functional and structural profiling revealed the strongest gradients in the iron-reducing zone, displaying the greatest variability with depth. This study provides important insight into biogeochemical characteristics in different regions of contaminant plumes, which will aid in improving models of contaminant fate and natural attenuation rates.     

Phylogenetic Diversity of Archaea and Bacteria in the Anoxic Zone of a Meromictic Lake (Lake Pavin, France)

The compositions of archaeal and bacterial populations at different depths (60 m [mixolimnion-chemocline interface], 70 m [chemocline-subchemocline interface], 90 m, and 92 m [the water-sediment interface]) in the anoxic zone of the water column in Lake Pavin, a freshwater permanently stratified mountain lake in France, were determined. Phylogenetic trees were constructed from sequences to assess archaeal and bacterial diversity at the four sites.     

Genetic Diversity of Microbial Eukaryotes in Anoxic Sediment of the Saline Meromictic Lake Namako-ike (Japan): On the Detection of Anaerobic or Anoxic-tolerant Lineages of Eukaryotes

Available sequence data on eukaryotic small-subunit ribosomal DNA (SSU rDNA) directly retrieved from various environments have increased recently, and the diversity of microbial eukaryotes (protists) has been shown to be much greater than previously expected. However, the molecular information accumulated to date does still not thoroughly reveal ecological distribution patterns of microbial eukaryotes. In the ongoing challenge to detect anaerobic or anoxic-tolerant lineages of eukaryotes, we directly extracted DNA from the anoxic sediment of a saline meromictic lake, constructed genetic libraries of PCR-amplified SSU rDNA, and performed phylogenetic analyses with the cloned SSU rDNA sequences. Although a few sequences could not be confidently assigned to any major eukaryotic groups in the analyses and are debatable regarding their taxonomic positions, most sequences obtained have affiliations with known major lineages of eukaryotes (Cercozoa, Alveolata, Stramenopiles, and Opisthokonta). Among these sequences, some branched with lineages predominantly composed of uncultured environmental clones retrieved from other anoxic environments, while others were closely related to those of eukaryotic parasites (e.g. Phytomyxea of Cercozoa, Gregarinea of Alveolata, and Ichthyosporea of Opisthokonta).     

Among-habitat and Temporal Variability of Selected Macroinvertebrate Based Metrics in a Mediterranean Shallow Lake (NW Spain)

According to the European Water Framework Directive, waterbodies have to be classified on the basis of their ecological status using biological quality elements, such as macroinvertebrates. This needs to take into consideration the influence of natural variation (both spatial and temporal) of reference biological communities as this may obscure the effects caused by anthropogenic disturbance. We studied the influence of among-habitat and temporal (seasonal and interannual) changes on the macroinvertebrate communities of an Iberian shallow lake and the variability of 21 measures potentially useful for bioindication purposes. Two series of data were examined: (a) macroinvertebrate samples taken on four occasions over an annual cycle were used to assess the effects of seasonality and among-habitat variability; (b) macroinvertebrate samples collected in three consecutive summers were used to assess interannual variability. Coefficients of variation, expressed as percentage, were used to quantify the effect of among-habitat and temporal variability on the selected metrics. According to our results, % Insecta, Shannon–Wiener diversity index and the qualitative taxonomic metrics (measures based on number of taxa) were robust in terms of temporal (seasonal and interannual) and among-habitat variability. Abundance ratio and some metrics based on functional feeding groups were highly variable. Therefore, qualitative taxonomic metrics may be promising tools in biomonitoring programs of Mediterranean shallow lakes due to their comparatively low variability.