Pyrosequencing Libraries Research Articles

Temporal shifts in dominant sulfur-oxidizing chemoautotrophic populations across the Cariaco Basin's redoxcline

A pronounced and relatively stable peak in dark carbon assimilation (DCA) has been observed within the Cariaco Basin's redoxcline over a 19-year observation period, presumably driven largely by sulfur-oxidizing chemoautotrophic bacteria (thioautotrophs). As documented in previous reports, this midwater productivity hotspot is associated with prominent biomass peaks, consisting of prokaryoplankton (Bacteria + Archaea) protists and viruses. Early in the time-series, phylogenetic studies (small subunit ribosomal RNA gene libraries) and fluorescent in situ hybridization (FISH) surveys documented that the chemoautotrophic layer was overwhelmingly populated by ε-proteobacteria related to sulfur-oxidizing symbionts from hydrothermal vent systems and β-proteobacteria. However, after May 2009, β- and ε-proteobacteria were no longer detected in significant numbers by FISH or by more exhaustive sequencing efforts (454 pyrosequencing and MiSeq iTag libraries). Sulfur oxidation gene quantification (qPCR) also confirmed that ε-proteobacteria were rare in samples collected after mid-2009 and that the chemoautotrophic layer became dominated by γ-proteobacterial sulfur-oxidizers (GSOs). Monthly hydrographic and chemical data from the CARIACO Ocean Time-Series were examined for temporal changes that might have driven this shift in the dominant thioautotrophs. Within the redoxcline (250–450 m), significant shifts occurred in water density structure and distributions of H2S, NH4+, O2, and NO2- between the β-/ε-proteobacteria-dominated and the γ-proteobacteria-dominated periods. Contrary to observations in other euxinic systems, GSO seem to have a selective advantage over other sulfur-oxidizing lieages during periods of higher and less variable H2S vertical fluxes and higher relative fluxes of H2S to O2 and NO3-. How exactly this changing seascape might have selected for one lineage over another cannot be directly determined. However, spatiotemporal variations in water density and concentration gradients altered vertical fluxes of the major reductant (H2S) and oxidants (O2, NO3-,NO2-), which we postulate provided a more geochemically stable environment favoring GSOs.

Monitoring of microbial communities in anaerobic digestion sludge for biogas optimisation

This study characterised and compared the microbial communities of anaerobic digestion (AD) sludge using three different methods – (1) Clone library; (2) Pyrosequencing; and (3) Terminal restriction fragment length polymorphism (T-RFLP). Although high-throughput sequencing techniques are becoming increasingly popular and affordable, the reliance of such techniques for frequent monitoring of microbial communities may be a financial burden for some. Furthermore, the depth of microbial analysis revealed by high-throughput sequencing may not be required for monitoring purposes. This study aims to develop a rapid, reliable and economical approach for the monitoring of microbial communities in AD sludge. A combined approach where genetic information of sequences from clone library was used to assign phylogeny to T-RFs determined experimentally was developed in this study. In order to assess the effectiveness of the combined approach, microbial communities determined by the combined approach was compared to that characterised by pyrosequencing. Results showed that both pyrosequencing and clone library methods determined the dominant bacteria phyla to be Proteobacteria, Firmicutes, Bacteroidetes, and Thermotogae. Both methods also found that sludge A and B were predominantly dominated by acetogenic methanogens followed by hydrogenotrophic methanogens. The number of OTUs detected by T-RFLP was significantly lesser than that detected by the clone library. In this study, T-RFLP analysis identified majority of the dominant species of the archaeal consortia. However, many of the more highly diverse bacteria consortia were missed. Nevertheless, the combined approach developed in this study where clone sequences from the clone library were used to assign phylogeny to T-RFs determined experimentally managed to accurately predict the same dominant microbial groups for both sludge A and sludge B, as compared to the pyrosequencing results. Results showed that the combined approach of clone library and T-RFLP accurately predicted the dominant microbial groups and thus is a reliable and more economical way to monitor the evolution of microbial systems in AD sludge.

K-shuff: A Novel Algorithm for Characterizing Structural and Compositional Diversity in Gene Libraries.

K-shuff is a new algorithm for comparing the similarity of gene sequence libraries, providing measures of the structural and compositional diversity as well as the significance of the differences between these measures. Inspired by Ripley’s K-function for spatial point pattern analysis, the Intra K-function or IKF measures the structural diversity, including both the richness and overall similarity of the sequences, within a library. The Cross K-function or CKF measures the compositional diversity between gene libraries, reflecting both the number of OTUs shared as well as the overall similarity in OTUs. A Monte Carlo testing procedure then enables statistical evaluation of both the structural and compositional diversity between gene libraries. For 16S rRNA gene libraries from complex bacterial communities such as those found in seawater, salt marsh sediments, and soils, K-shuff yields reproducible estimates of structural and compositional diversity with libraries greater than 50 sequences. Similarly, for pyrosequencing libraries generated from a glacial retreat chronosequence and Illumina® libraries generated from US homes, K-shuff required >300 and 100 sequences per sample, respectively. Power analyses demonstrated that K-shuff is sensitive to small differences in Sanger or Illumina® libraries. This extra sensitivity of K-shuff enabled examination of compositional differences at much deeper taxonomic levels, such as within abundant OTUs. This is especially useful when comparing communities that are compositionally very similar but functionally different. K-shuff will therefore prove beneficial for conventional microbiome analysis as well as specific hypothesis testing.

Conservation of the endangered Mediterranean tortoise Testudo hermanni hermanni: The contribution of population genetics and historical demography

Estimating the genetic variation and demographic trends of species in decline is of major concern in conservation genetics. This study contributes to understanding how historical and anthropogenic factors shape the distribution of current genetic diversity in one of the most endangered reptiles in Western Europe, the Hermann's tortoise Testudo hermanni hermanni. We used 17 microsatellite loci, chosen from a pyrosequencing library specifically developed for the subspecies to genotype eight populations distributed over about 30 sample localities covering almost the entire geographic distribution of the sub-species. The population genetic results reflect a very strong genetic structure and identify three major clusters among the Hermann's tortoise in the occidental Mediterranean basin: a continental cluster (Albera in Spain, Var in France and continental Italy), an insular cluster (Corsica, Sardinia and Sicily) and a cluster of mixed lineage (Minorca). Each of the eight studied populations is characterized by only one lineage except for Minorca, in which two lineages were identified. In contrast to what most empirical studies predict, the genetic diversity of the Hermann's tortoise is greater on islands than on the continent. Islands can therefore be considered as genetic sanctuaries with high conservation potential for this species, particularly in Corsica. Historical demographic patterns inferred with a generalized stepwise-mutation model (GSM) using maximum likelihood showed significant past demographic changes in only two of the eight sampled populations: a demographic bottleneck was detected in the Albera population (Spain) and a demographic expansion in Corsica (France). In contrast to what was argued in previous studies, these results indicate that the Hermann's tortoise lineage found in Corsica is autochthonous. The origin of both lineages found in Minorca remains speculative. Lastly, our study identifies the six most relevant management units (sensu Moritz, 1994) for conservation purposes on the basis that they represent a significant part of the evolutionary legacy of the species. Some conservation recommendations were proposed, in particular for the most threatened population in Albera.

Diversity of Rhodopirellula and related planctomycetes in a North Sea coastal sediment employing carB as molecular marker.

Rhodopirellula is an abundant marine member of the bacterial phylum Planctomycetes. Cultivation studies revealed the presence of several closely related Rhodopirellula species in European coastal sediments. Because the 16S rRNA gene does not provide the desired taxonomic resolution to differentiate Rhodopirellula species, we performed a comparison of the genomes of nine Rhodopirellula strains and six related planctomycetes and identified carB, coding for the large subunit of carbamoylphosphate synthetase, as a suitable molecular marker. In this study, we investigated the diversity of Rhodopirellula in coastal intertidal surface sediments of Sylt island, North Sea, using the 16S rRNA and carB genes as molecular markers. The carB clone and pyrosequencing libraries revealed the presence of 12 species of Rhodopirellula and of 66 species in closely related undescribed genera, a diversity that was not detected with a 16S rRNA gene library. This study demonstrates that the carB gene is a powerful molecular marker for detecting Rhodopirellula species in the environment and may be used for the taxonomic evaluation of new strains.

Microbial community structure accompanied with electricity production in a constructed wetland plant microbial fuel cell

This study reveals the complex structure of bacterial and archaeal communities associated with a Canna indica plant microbial fuel cell (PMFC) and its electricity production. The PMFC produced a maximum current of 105mA/m2 by utilizing rhizodeposits as the sole electron donor without any external nutrient or buffer supplements, which demonstrates the feasibility of PMFCs in practical oligotrophic conditions with low solution conductivity. The microbial diversity was significantly higher in the PMFC than non-plant controls or sediment-only controls, and pyrosequencing and clone library reveal that rhizodeposits conversion to current were carried out by syntrophic interactions between fermentative bacteria (e.g., Anaerolineaceae) and electrochemically active bacteria (e.g., Geobacter). Denitrifying bacteria and acetotrophic methanogens play a minor role in organics degradation, but abundant hydrogenotrophic methanogens and thermophilic archaea are likely main electron donor competitors.

An improved method to set significance thresholds for β diversity testing in microbial community comparisons.

Exploring the variation in microbial community diversity between locations (β diversity) is a central topic in microbial ecology. Currently, there is no consensus on how to set the significance threshold for β diversity. Here, we describe and quantify the technical components of β diversity, including those associated with the process of subsampling. These components exist for any proposed β diversity measurement procedure. Further, we introduce a strategy to set significance thresholds for β diversity of any group of microbial samples using rarefaction, invoking the notion of a meta-community. The proposed technique was applied to several in silico generated operational taxonomic unit (OTU) libraries and experimental 16S rRNA pyrosequencing libraries. The latter represented microbial communities from different biological rapid sand filters at a full-scale waterworks. We observe that β diversity, after subsampling, is inflated by intra-sample differences; this inflation is avoided in the proposed method. In addition, microbial community evenness (Gini > 0.08) strongly affects all β diversity estimations due to bias associated with rarefaction. Where published methods to test β significance often fail, the proposed meta-community-based estimator is more successful at rejecting insignificant β diversity values. Applying our approach, we reveal the heterogeneous microbial structure of biological rapid sand filters both within and across filters.

Ammonia-oxidizing bacteria and archaea in wastewater treatment plant sludge and nearby coastal sediment in an industrial area in China.

Under the increasing pressure of human activities, Hangzhou Bay has become one of the most seriously polluted waters along China's coast. Considering the excessive inorganic nitrogen detected in the bay, in this study, the impact of an effluent from a coastal industrial park on ammonia-oxidizing microorganisms (AOMs) of the receiving area was interpreted for the first time by molecular technologies. Revealed by real-time PCR, the ratio of archaeal amoA/bacterial amoA ranged from 5.68 × 10(-6) to 4.79 × 10(-5) in the activated sludge from two wastewater treatment plants (WWTPs) and 0.54-3.44 in the sediments from the effluent receiving coastal area. Analyzed by clone and pyrosequencing libraries, genus Nitrosomonas was the predominant ammonia-oxidizing bacteria (AOB), but no ammonia-oxidizing archaea (AOA) was abundant enough for sequencing in the activated sludge from the WWTPs; genus Nitrosomonas and Nitrosopumilus were the dominant AOB and AOA, respectively, in the coastal sediments. The different abundance of AOA but similar structure of AOB between the WWTPs and nearby coastal area probably indicated an anthropogenic impact on the microbial ecology in Hangzhou Bay.

Microbiota of the major South Atlantic reef building coral Mussismilia.

The Brazilian endemic scleractinian corals, genus Mussismilia, are among the main reef builders of the South Atlantic and are threatened by accelerating rates of disease. To better understand how holobiont microbial populations interact with corals during health and disease and to evaluate whether selective pressures in the holobiont or neutral assembly shape microbial composition, we have examined the microbiota structure of Mussismilia corals according to coral lineage, environment, and disease/health status. Microbiota of three Mussismilia species (Mussismilia harttii, Mussismilia hispida, and Mussismilia braziliensis) was compared using 16S rRNA pyrosequencing and clone library analysis of coral fragments. Analysis of biological triplicates per Mussismilia species and reef site allowed assessment of variability among Mussismilia species and between sites for M. braziliensis. From 173,487 V6 sequences, 6,733 coral- and 1,052 water-associated operational taxonomic units (OTUs) were observed. M. braziliensis microbiota was more similar across reefs than to other Mussismilia species microbiota from the same reef. Highly prevalent OTUs were more significantly structured by coral lineage and were enriched in Alpha- and Gammaproteobacteria. Bacterial OTUs from healthy corals were recovered from a M. braziliensis skeleton sample at twice the frequency of recovery from water or a diseased coral suggesting the skeleton is a significant habitat for microbial populations in the holobiont. Diseased corals were enriched with pathogens and opportunists (Vibrios, Bacteroidetes, Thalassomonas, and SRB). Our study examines for the first time intra- and inter-specific variability of microbiota across the genus Mussismilia. Changes in microbiota may be useful indicators of coral health and thus be a valuable tool for coral reef management and conservation.

Methane production and microbial community structure for alkaline pretreated waste activated sludge

Alkaline pretreatment was studied to analyze the influence on waste activated sludge (WAS) reduction, methane production and microbial community structure during anaerobic digestion. Methane production from alkaline pretreated sludge (A-WAS) (pH=12) increased from 251.2mL/Ld to 362.2mL/Ld with the methane content of 68.7% compared to raw sludge (R-WAS). Sludge reduction had been improved, and volatile suspended solids (VSS) removal rate and protein reduction had increased by ∼10% and ∼35%, respectively. The bacterial and methanogenic communities were analyzed using 454 pyrosequencing and clone libraries of 16S rRNA gene. Remarkable shifts were observed in microbial community structures after alkaline pretreatment, especially for Archaea. The dominant methanogenic population changed from Methanosaeta for R-WAS to Methanosarcina for A-WAS. In addition to the enhancement of solubilization and hydrolysis of anaerobic digestion of WAS, alkaline pretreatment showed significant impacts on the enrichment and syntrophic interactions between microbial communities.

Multiplexing with three-primer PCR for rapid and economical microsatellite validation.

The next generation sequencing revolution has enabled rapid discovery of genetic markers, however, development of fully functioning new markers still requires a long and costly process of marker validation. This study reports a rapid and economical approach for the validation and deployment of polymorphic microsatellite markers obtained from a 454 pyrosequencing library of Atlantic cod, Gadus morhua, Linnaeus 1758. Primers were designed from raw reads to amplify specific amplicon size ranges, allowing effective PCR multiplexing. Multiplexing was combined with a three-primer PCR approach using four universal tails to label amplicons with separate fluorochromes. A total of 192 primer pairs were tested, resulting in 73 polymorphic markers. Of these, 55 loci were combined in six multiplex panels each containing between six and eleven markers. Variability of the loci was assessed on G. morhua from the Celtic Sea (n = 46) and the Scotian Shelf (n = 46), two locations that have shown genetic differentiation in previous studies. Multilocus F(ST) between the two samples was estimated at 0.067 (P = 0.001). After three loci potentially under selection were excluded, the global F(ST) was estimated at 0.043 (P = 0.001). Our technique combines three-primer and multiplex PCR techniques, allowing simultaneous screening and validation of relatively large numbers of microsatellite loci.

Comparison of bacterial community in aqueous and oil phases of water-flooded petroleum reservoirs using pyrosequencing and clone library approaches

Bacterial communities in both aqueous and oil phases of water-flooded petroleum reservoirs were characterized by molecular analysis of bacterial 16S rRNA genes obtained from Shengli Oil Field using DNA pyrosequencing and gene clone library approaches. Metagenomic DNA was extracted from the aqueous and oil phases and subjected to polymerase chain reaction amplification with primers targeting the bacterial 16S rRNA genes. The analysis by these two methods showed that there was a large difference in bacterial diversity between the aqueous and oil phases of the reservoir fluids, especially in the reservoirs with lower water cut. At a high phylogenetic level, the predominant bacteria detected by these two approaches were identical. However, pyrosequencing allowed the detection of more rare bacterial species than the clone library method. Statistical analysis showed that the diversity of the bacterial community of the aqueous phase was lower than that of the oil phase. Phylogenetic analysis indicated that the vast majority of sequences detected in the water phase were from members of the genus Arcobacter within the Epsilonproteobacteria, which is capable of degrading the intermediates of hydrocarbon degradation such as acetate. The oil phase of reservoir fluid samples was dominated by members of the genus Pseudomonas within the Gammaproteobacteria and the genus Sphingomonas within the Alphaproteobacteria, which have the ability to degrade crude oil through adherence to hydrocarbons under aerobic conditions. In addition, many anaerobes that could degrade the component of crude oil were also found in the oil phase of reservoir fluids, mainly in the reservoir with lower water cut. These were represented by Desulfovibrio spp., Thermodesulfovibrio spp., Thermodesulforhabdus spp., Thermotoga spp., and Thermoanaerobacterium spp. This research suggested that simultaneous analysis of DNA extracted from both aqueous and oil phases can facilitate a better understanding of the bacterial communities in water-flooded petroleum reservoirs.

Comparison of the Bacterial Microbiota in a Bale of Collected Cardboard Determined by 454 Pyrosequencing and Clone Library

Biofouling, the accumulation of microorganisms, is a major problem in paper mills processing paper and cardboard. This leads to the production of lower quality recycled products. Several studies have focused on the microbial content in the paper mill and the final products. Our aim was to determine the microbial biota in a bale of collected cardboard prior to entering the paper mill. Total genomic DNA was isolated and analyzed using two different methods for comparison purposes: 454 pyrosequencing and clone library. A total of 3268 V6-V8 454 pyrosequencing reads and 322 cloned V6-V8 16S rRNA nucleotide sequences were obtained. Both methods showed the presence of three major bacterial genera: Bacillus, Solibacillus and Paenibacillus, all members of the spore-forming phylum Firmicutes. Pyrosequencing, however, revealed a richer and more diverse bacterial community than clone library. It showed the presence of additional minor Firmicute genera and of a small number of Proteobacteria. The sorting at the recycling plant, the storing, and the processing at the paper mill, the end uses, will all contribute to the bacterial microbiota present in a bale of collected cardboard as revealed here.

Highly Heterogeneous Bacterial Communities Associated with the South China Sea Reef Corals Porites lutea, Galaxea fascicularis and Acropora millepora

Coral harbor diverse and specific bacteria play significant roles in coral holobiont function. Bacteria associated with three of the common and phylogenetically divergent reef-building corals in the South China Sea, Porites lutea, Galaxea fascicularis and Acropora millepora, were investigated using 454 barcoded-pyrosequencing. Three colonies of each species were sampled, and 16S rRNA gene libraries were constructed individually. Analysis of pyrosequencing libraries showed that bacterial communities associated with the three coral species were more diverse than previous estimates based on corals from the Caribbean Sea, Indo-Pacific reefs and the Red Sea. Three candidate phyla, including BRC1, OD1 and SR1, were found for the first time in corals. Bacterial communities were separated into three groups: P. lutea and G. fascicular, A. millepora and seawater. P. lutea and G. fascicular displayed more similar bacterial communities, and bacterial communities associated with A. millepora differed from the other two coral species. The three coral species shared only 22 OTUs, which were distributed in Alphaproteobacteria, Deltaproteobacteria, Gammaproteobacteria, Chloroflexi, Actinobacteria, Acidobacteria and an unclassified bacterial group. The composition of bacterial communities within each colony of each coral species also showed variation. The relatively small common and large specific bacterial communities in these corals implies that bacterial associations may be structured by multiple factors at different scales and that corals may associate with microbes in terms of similar function, rather than identical species.

Successional development of biofilms in moving bed biofilm reactor (MBBR) systems treating municipal wastewater

Biofilm-based technologies, such as moving bed biofilm reactor (MBBR) systems, are widely used to treat wastewater. Biofilm development is important for MBBR systems as much of the microbial biomass is retained within reactors as biofilm on suspended carriers. Little is known about this process of biofilm development and the microorganisms upon which MBBRs rely. We documented successional changes in microbial communities as biofilms established in two full-scale MBBR systems treating municipal wastewater over two seasons. 16S rRNA gene-targeted pyrosequencing and clone libraries were used to describe microbial communities. These data indicate a successional process that commences with the establishment of an aerobic community dominated by Gammaproteobacteria (up to 52 % of sequences). Over time, this community shifts towards dominance by putatively anaerobic organisms including Deltaproteobacteria and Clostridiales. Significant differences were observed between the two wastewater treatment plants (WWTPs), mostly due to a large number of sequences (up to 55 %) representing Epsilonproteobacteria (mostly Arcobacter) at one site. Archaea in young biofilms included several lineages of Euryarchaeota and Crenarchaeota. In contrast, the mature biofilm consisted entirely of Methanosarcinaceae (Euryarchaeota). This study provides new insights into the community structure of developing biofilms at full-scale WWTPs and provides the basis for optimizing MBBR start-up and operational parameters.

Ectomycorrhizal-dominated boreal and tropical forests have distinct fungal communities, but analogous spatial patterns across soil horizons.

Fungi regulate key nutrient cycling processes in many forest ecosystems, but their diversity and distribution within and across ecosystems are poorly understood. Here, we examine the spatial distribution of fungi across a boreal and tropical ecosystem, focusing on ectomycorrhizal fungi. We analyzed fungal community composition across litter (organic horizons) and underlying soil horizons (0–20 cm) using 454 pyrosequencing and clone library sequencing. In both forests, we found significant clustering of fungal communities by site and soil horizons with analogous patterns detected by both sequencing technologies. Free-living saprotrophic fungi dominated the recently-shed leaf litter and ectomycorrhizal fungi dominated the underlying soil horizons. This vertical pattern of fungal segregation has also been found in temperate and European boreal forests, suggesting that these results apply broadly to ectomycorrhizal-dominated systems, including tropical rain forests. Since ectomycorrhizal and free-living saprotrophic fungi have different influences on soil carbon and nitrogen dynamics, information on the spatial distribution of these functional groups will improve our understanding of forest nutrient cycling.

Hydrocarbon-degrading bacteria enriched by the Deepwater Horizon oil spill identified by cultivation and DNA-SIP

The massive influx of crude oil into the Gulf of Mexico during the Deepwater Horizon (DWH) disaster triggered dramatic microbial community shifts in surface oil slick and deep plume waters. Previous work had shown several taxa, notably DWH Oceanospirillales, Cycloclasticus and Colwellia, were found to be enriched in these waters based on their dominance in conventional clone and pyrosequencing libraries and were thought to have had a significant role in the degradation of the oil. However, this type of community analysis data failed to provide direct evidence on the functional properties, such as hydrocarbon degradation of organisms. Using DNA-based stable-isotope probing with uniformly (13)C-labelled hydrocarbons, we identified several aliphatic (Alcanivorax, Marinobacter)- and polycyclic aromatic hydrocarbon (Alteromonas, Cycloclasticus, Colwellia)-degrading bacteria. We also isolated several strains (Alcanivorax, Alteromonas, Cycloclasticus, Halomonas, Marinobacter and Pseudoalteromonas) with demonstrable hydrocarbon-degrading qualities from surface slick and plume water samples collected during the active phase of the spill. Some of these organisms accounted for the majority of sequence reads representing their respective taxa in a pyrosequencing data set constructed from the same and additional water column samples. Hitherto, Alcanivorax was not identified in any of the previous water column studies analysing the microbial response to the spill and we discuss its failure to respond to the oil. Collectively, our data provide unequivocal evidence on the hydrocarbon-degrading qualities for some of the dominant taxa enriched in surface and plume waters during the DWH oil spill, and a more complete understanding of their role in the fate of the oil.

Characterization of thirty two microsatellite loci for three Atlanto-Mediterranean echinoderm species

Thirty two microsatellites were optimized from 454 pyrosequencing libraries for three Atlanto-Mediterranean echinoderms: Coscinasterias tenuispina, Echinaster sepositus and Arbacia lixula. We observed different frequency of microsatellite types (di-, tri-, tetra- and penta-nucleotide) throughout the genome of the species, but no significant differences were observed in allele richness among different microsatellite repeats. No loci showed linkage disequilibrium. Heterozygosity deficit and departure from Hardy–Weinberg equilibrium were observed for some loci, in two species, probably due to high levels of inbreeding. Heterozygosity excess observed in C. tenuispina could be explained by selection against homozygotes and/or outcrossing.

Characterization of bacterial diversity in drinking water by pyrosequencing

Controlling microbial contamination of drinking water is critical to public health. However, understanding of the microbial ecology of drinking water remains incomplete. Representing the first application of high-throughput sequencing in drinking water microbiology, the objective of this study is to evaluate pyrosequencing as a high-throughput technique for the characterization of bacterial diversity in drinking water in comparison with conventional clone library analysis. Pyrosequencing and clone library analysis were performed in parallel to study the bacterial community composition in drinking water samples following the concentration of microbial biomass in drinking water with ultrafiltration. Validated by clone library analysis, pyrosequencing was confirmed as a highly efficient deep-sequencing technique to characterize the bacterial diversity in drinking water. Sequences of Alphaproteobacteria and Betaproteobacteria dominated the bacterial community in drinking water with Oxalobacteraceae and Methylobacteriaceae as the most abundant bacterial families, which is consistent with the prominent abundance of these populations frequently detected in various freshwater environments where source waters originate. Bacterial populations represented by the most abundant sequences in drinking water were closely related to cultures of metabolically versatile bacterial taxa widely distributed in the environment, suggesting a potential link between environmental distribution, metabolic characteristics, and abundance in drinking water.

Exploring Pandora's Box: Potential and Pitfalls of Low Coverage Genome Surveys for Evolutionary Biology

High throughput sequencing technologies are revolutionizing genetic research. With this “rise of the machines”, genomic sequences can be obtained even for unknown genomes within a short time and for reasonable costs. This has enabled evolutionary biologists studying genetically unexplored species to identify molecular markers or genomic regions of interest (e.g. micro- and minisatellites, mitochondrial and nuclear genes) by sequencing only a fraction of the genome. However, when using such datasets from non-model species, it is possible that DNA from non-target contaminant species such as bacteria, viruses, fungi, or other eukaryotic organisms may complicate the interpretation of the results. In this study we analysed 14 genomic pyrosequencing libraries of aquatic non-model taxa from four major evolutionary lineages. We quantified the amount of suitable micro- and minisatellites, mitochondrial genomes, known nuclear genes and transposable elements and searched for contamination from various sources using bioinformatic approaches. Our results show that in all sequence libraries with estimated coverage of about 0.02–25%, many appropriate micro- and minisatellites, mitochondrial gene sequences and nuclear genes from different KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways could be identified and characterized. These can serve as markers for phylogenetic and population genetic analyses. A central finding of our study is that several genomic libraries suffered from different biases owing to non-target DNA or mobile elements. In particular, viruses, bacteria or eukaryote endosymbionts contributed significantly (up to 10%) to some of the libraries analysed. If not identified as such, genetic markers developed from high-throughput sequencing data for non-model organisms may bias evolutionary studies or fail completely in experimental tests. In conclusion, our study demonstrates the enormous potential of low-coverage genome survey sequences and suggests bioinformatic analysis workflows. The results also advise a more sophisticated filtering for problematic sequences and non-target genome sequences prior to developing markers.