Extensive profiling of a complex microbial community by high-throughput sequencing

Complex microbial communities remain poorly characterized despite their ubiquity and importance to human and animal health, agriculture, and industry. Attempts to describe microbial communities by either traditional microbiological methods or molecular methods have been limited in both scale and precision. The availability of genomics technologies offers an unprecedented opportunity to conduct more comprehensive characterizations of microbial communities. Here we describe the application of an established molecular diagnostic method based on the chaperonin-60 sequence, in combination with high-throughput sequencing, to the profiling of a microbial community: the pig intestinal microbial community. Four libraries of cloned cpn60 sequences were generated by two genomic DNA extraction procedures in combination with two PCR protocols. A total of 1,125 cloned cpn60 sequences from the four libraries were sequenced. Among the 1,125 cloned cpn60 sequences, we identified 398 different nucleotide sequences encoding 280 unique peptide sequences. Pairwise comparisons of the 398 unique nucleotide sequences revealed a high degree of sequence diversity within the library. Identification of the likely taxonomic origins of cloned sequences ranged from imprecise, with clones assigned to a taxonomic subclass, to precise, for cloned sequences with 100% DNA sequence identity with a species in our reference database. The compositions of the four libraries were compared and differences related to library construction parameters were observed. Our results indicate that this method is an alternative to 16S rRNA sequence-based studies which can be scaled up for the purpose of performing a potentially comprehensive assessment of a given microbial community or for comparative studies.     

Use of Subtractive Hybridization To Design Habitat-Based Oligonucleotide Probes for Investigation of Natural Bacterial Communities

We describe a rapid oligonucleotide probe design strategy based on subtractive hybridization which yields probes for 16S rRNA or rRNA genes of individual members of microbial communities that are specific within the context of those communities. This strategy circumvents the need to sequence many similar or identical clones of dominant members of a community. Radioactively labeled subfragments of a cloned 16S rRNA gene sequence for which a probe is required (target) were hybridized with biotinylated total 16S ribosomal DNA (rDNA) amplified from the microbial community, and the hybrids formed were subsequently discarded. The remaining enriched fragments were used to screen a library consisting of cloned subfragments of the target sequence by colony hybridization in order to identify the variable regions of the 16S rRNA gene with the required specificity. The sequencing of random clones in one 16S rDNA library demonstrated that only those clones with 100% sequence identity with the probe fragment were detected by it. Moreover, sequencing of other, randomly selected, probe-positive clones revealed 100% sequence identity with the probe. Probes developed in this way tended to correspond to more variable regions of the 16S rRNA if the target sequences were similar to the sequences of other clones in the library and to less variable regions if the target sequences were phylogenetically isolated within the clone library. Although the absolute specificity of the latter probes, as assessed by comparison with available database sequences, was lower than the absolute specificity of the probes from the more variable regions, they were specific within the context of the environmental samples from which they were derived.     

Prospecting for novel biocatalysts in a soil metagenome

The metagenomes of complex microbial communities are rich sources of novel biocatalysts. We exploited the metagenome of a mixed microbial population for isolation of more than 15 different genes encoding novel biocatalysts by using a combined cultivation and direct cloning strategy. A 16S rRNA sequence analysis revealed the presence of hitherto uncultured microbes closely related to the genera Pseudomonas, Agrobacterium, Xanthomonas, Microbulbifer, and Janthinobacterium. Total genomic DNA from this bacterial community was used to construct cosmid DNA libraries, which were functionally searched for novel enzymes of biotechnological value. Our searches in combination with cosmid sequencing resulted in identification of four clones encoding 12 putative agarase genes, most of which were organized in clusters consisting of two or three genes. Interestingly, nine of these agarase genes probably originated from gene duplications. Furthermore, we identified by DNA sequencing several other biocatalyst-encoding genes, including genes encoding a putative stereoselective amidase (amiA), two cellulases (gnuB and uvs080), an alpha-amylase (amyA), a 1,4-alpha-glucan branching enzyme (amyB), and two pectate lyases (pelA and uvs119). Also, a conserved cluster of two lipase genes was identified, which was linked to genes encoding a type I secretion system. The novel gene aguB was overexpressed in Escherichia coli, and the enzyme activities were determined. Finally, we describe more than 162 kb of DNA sequence that provides a strong platform for further characterization of this microbial consortium.     

Molecular analysis of shower curtain biofilm microbes

Households provide environments that encourage the formation of microbial communities, often as biofilms. Such biofilms constitute potential reservoirs for pathogens, particularly for immune-compromised individuals. One household environment that potentially accumulates microbial biofilms is that provided by vinyl shower curtains. Over time, vinyl shower curtains accumulate films, commonly referred to as "soap scum," which microscopy reveals are constituted of lush microbial biofilms. To determine the kinds of microbes that constitute shower curtain biofilms and thereby to identify potential opportunistic pathogens, we conducted an analysis of rRNA genes obtained by PCR from four vinyl shower curtains from different households. Each of the shower curtain communities was highly complex. No sequence was identical to one in the databases, and no identical sequences were encountered in the different communities. However, the sequences generally represented similar phylogenetic kinds of organisms. Particularly abundant sequences represented members of the alpha-group of proteobacteria, mainly Sphingomonas spp. and Methylobacterium spp. Both of these genera are known to include opportunistic pathogens, and several of the sequences obtained from the environmental DNA samples were closely related to known pathogens. Such organisms have also been linked to biofilm formation associated with water reservoirs and conduits. In addition, the study detected many other kinds of organisms at lower abundances. These results show that shower curtains are a potential source of opportunistic pathogens associated with biofilms. Frequent cleaning or disposal of shower curtains is indicated, particularly in households with immune-compromised individuals.     

Molecular Analysis of Shower Curtain Biofilm Microbes

Households provide environments that encourage the formation of microbial communities, often as biofilms. Such biofilms constitute potential reservoirs for pathogens, particularly for immune-compromised individuals. One household environment that potentially accumulates microbial biofilms is that provided by vinyl shower curtains. Over time, vinyl shower curtains accumulate films, commonly referred to as “soap scum,” which microscopy reveals are constituted of lush microbial biofilms. To determine the kinds of microbes that constitute shower curtain biofilms and thereby to identify potential opportunistic pathogens, we conducted an analysis of rRNA genes obtained by PCR from four vinyl shower curtains from different households. Each of the shower curtain communities was highly complex. No sequence was identical to one in the databases, and no identical sequences were encountered in the different communities. However, the sequences generally represented similar phylogenetic kinds of organisms. Particularly abundant sequences represented members of the α-group of proteobacteria, mainly Sphingomonas spp. and Methylobacterium spp. Both of these genera are known to include opportunistic pathogens, and several of the sequences obtained from the environmental DNA samples were closely related to known pathogens. Such organisms have also been linked to biofilm formation associated with water reservoirs and conduits. In addition, the study detected many other kinds of organisms at lower abundances. These results show that shower curtains are a potential source of opportunistic pathogens associated with biofilms. Frequent cleaning or disposal of shower curtains is indicated, particularly in households with immune-compromised individuals.     

Diversity of cyanobacterial species and phylotypes in biofilms from the littoral zone of Lake Baikal

The majority of naturally occurring biofilms contain numerous microorganisms that have not yet been cultured. Additionally, there is little information available regarding the genetic structure and species diversity of these communities. Therefore, we characterised the species diversity, structure and metagenome of biofilms grown on stones and steel plates in the littoral zone of Lake Baikal (East Siberia, Russia) by applying three different approaches. First, light microscopy enabled identification of the species diversity of biofilm-forming cyanobacteria on different substrates with the dominance of Rivularia rufescens, Tolypothrix limbata, Chamaesiphon fuscus, Ch. subglobosus, and Heteroleibleinia pusilla. Additionally, scanning electron microscopy was used to show the spatial structure of biofilms. Finally, sequence analysis of 30,660 16S rRNA clones indicated a high diversity within the biofilm communities, with the majority of the microbes being closely related to Cyanobacteria (8–46% sequences), Proteobacteria (14–43%), and Bacteroidetes (10–41%). Rivularia sp., Pseudanabaena sp., and Chamaesiphon spp. were the dominant cyanobacterial phylotypes.     

Turkey fecal microbial community structure and functional gene diversity revealed by 16S rRNA gene and metagenomic sequences

The primary goal of this study was to better understand the microbial composition and functional genetic diversity associated with turkey fecal communities. To achieve this, 16S rRNA gene and metagenomic clone libraries were sequenced from turkey fecal samples. The analysis of 382 16S rRNA gene sequences showed that the most abundant bacteria were closely related to Lactobacillales (47%), Bacillales (31%), and Clostridiales (11%). Actinomycetales, Enterobacteriales, and Bacteroidales sequences were also identified, but represented a smaller part of the community. The analysis of 379 metagenomic sequences showed that most clones were similar to bacterial protein sequences (58%). Bacteriophage (10%) and avian viruses (3%) sequences were also represented. Of all metagenomic clones potentially encoding for bacterial proteins, most were similar to low G+C Gram-positive bacterial proteins, particularly from Lactobacillales (50%), Bacillales (11%), and Clostridiales (8%). Bioinformatic analyses suggested the presence of genes encoding for membrane proteins, lipoproteins, hydrolases, and functional genes associated with the metabolism of nitrogen and sulfur containing compounds. The results from this study further confirmed the predominance of Firmicutes in the avian gut and highlight the value of coupling 16S rRNA gene and metagenomic sequencing data analysis to study the microbial composition of avian fecal microbial communities.     

A glimpse under the rim – the composition of microbial biofilm communities in domestic toilets

Aim: To determine the microbial composition of biofilms in domestic toilets by molecular means. Methods and Results: Genomic DNA was extracted from six biofilm samples originating from households around Düsseldorf, Germany. While no archaeal 16S rRNA or fungal ITS genes were detected by PCR, fingerprinting of bacterial 16S rRNA genes revealed a diverse community in all samples. These communities also differed considerably between the six biofilms. Using the Ribosomal Database Project (RDP) classifier tool, 275 cloned 16S rRNA gene sequences were assigned to 11 bacterial phyla and 104 bacterial genera. Only 15 genera (representing 121 sequences affiliated with Acidobacteria, Actinobacteria, Bacteroidetes, Planctomycetes and Proteobacteria) occurred in at least half of the samples or contributed at least 10% of the sequences in a single biofilm. These sequences were defined as ‘typical’ for toilet biofilms, and they were examined in more detail. On a 97% sequence similarity level, these sequences represented 56 species. Twelve of these were closely related to well‐described bacterial species, and only two of them were categorized as belonging to risk group 2. No 16S rRNA genes of typical faecal bacteria were detected in any sample. Virtually all ‘typical’ clones were found to be closely related to bacteria or to sequences obtained from environmental sources, implicating that the flushing water is the main source of recruitment. Conclusion: In view of the great diversity of mostly yet‐uncultured bacteria and the considerable differences between individual toilets, very general strategies appear to be most suited for the removal and prevention of toilet biofilms. Significance and Impact of the Study: For the first time, a molecular fingerprinting and cloning approach was used to monitor the species composition in biofilm samples taken from domestic toilets. Knowledge about the microbial composition of biofilms in domestic toilets is a prerequisite for developing and evaluating strategies for their removal and prevention.     

Serial analysis of V6 ribosomal sequence tags (SARST-V6): a method for efficient, high-throughput analysis of microbial community composition

Serial analysis of ribosomal sequence tags (SARST) is a novel technique for characterizing microbial community composition. The SARST method captures sequence information from concatemers of short 16S rDNA polymerase chain reaction (PCR) amplicons from complex populations of DNA. Here, we describe a similar method, serial analysis of V6 ribosomal sequence tags (SARST-V6), which targets the V6 hypervariable region of bacterial 16S rRNA genes. The SARST-V6 technique exploits internal primer sequences to generate compatible restriction digest overhangs, thereby improving upon the efficiency of SARST. Serial analysis of V6 ribosomal sequence tags of bacterial community composition in hydrothermal marine sediments from Guaymas Basin resembled results of cloning and sequencing of single, full-length PCR products from ribosomal RNA genes of the same microbial community. Both methods identified the same major bacterial groups, but only SARST-V6 recovered thermodesulfobacteria and gamma-proteobacteria sequences, while only full-length PCR product cloning recovered candidate division OP11 se-quences. There were differences in the relative frequencies of some phylotypes. The disparities reflect differences in the amplicon pool obtained during initial amplification that may result from different primer affinities or DNA degradation. These results demonstrate the utility of SARST-V6 in collecting taxonomically informative data for high-throughput analysis of microbial communities.     

Phylogenetic analyses of ribosomal DNA-containing bacterioplankton genome fragments from a 4000 m vertical profile in the North Pacific Subtropical Gyre

High-throughput identification of rRNA gene-containing clones in large insert metagenomic libraries is difficult, because of the high background of host ribosomal RNA (rRNA) and rRNA genes. To address this challenge, a membrane hybridization method was developed to identify all bacterial small subunit rRNA-containing fosmid clones of microbial community DNA from seven different depths in the North Pacific Subtropical Gyre. Out of 101,376 clones screened, 751 rDNA-containing clones were identified that grouped in ∼60 different clades. Several rare sequences only remotely related to known groups were detected, including a Wolbachia -related sequence containing a putative intron or intervening sequence, as well as seven sequences from Order Myxococcales not previously detected in pelagic habitats. Stratified, depth-specific population structure was evident within both cultured and uncultured lineages. Conversely, some eurybathyal members of the genera Alcanivorax and Rhizobium shared identical small subunit ribosomal DNA sequences that were distributed from surface waters to the 4000 m depth. Comparison with similar analyses in Monterey Bay microbial communities revealed previously recognized, as well as some distinctive, depth-stratified partitioning that distinguished coastal from open ocean bacterioplankton populations. While some bias was evident in fosmid clone recovery in a few particular lineages, the overall phylogenetic group recovery and distributions were consistent with previous studies, as well as with direct shotgun sequence data from the same source DNA.     

Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies

The advent of next generation sequencing has coincided with a growth in interest in using these approaches to better understand the role of the structure and function of the microbial communities in human, animal, and environmental health. Yet, use of next generation sequencing to perform 16S rRNA gene sequence surveys has resulted in considerable controversy surrounding the effects of sequencing errors on downstream analyses. We analyzed 2.7×10(6) reads distributed among 90 identical mock community samples, which were collections of genomic DNA from 21 different species with known 16S rRNA gene sequences; we observed an average error rate of 0.0060. To improve this error rate, we evaluated numerous methods of identifying bad sequence reads, identifying regions within reads of poor quality, and correcting base calls and were able to reduce the overall error rate to 0.0002. Implementation of the PyroNoise algorithm provided the best combination of error rate, sequence length, and number of sequences. Perhaps more problematic than sequencing errors was the presence of chimeras generated during PCR. Because we knew the true sequences within the mock community and the chimeras they could form, we identified 8% of the raw sequence reads as chimeric. After quality filtering the raw sequences and using the Uchime chimera detection program, the overall chimera rate decreased to 1%. The chimeras that could not be detected were largely responsible for the identification of spurious operational taxonomic units (OTUs) and genus-level phylotypes. The number of spurious OTUs and phylotypes increased with sequencing effort indicating that comparison of communities should be made using an equal number of sequences. Finally, we applied our improved quality-filtering pipeline to several benchmarking studies and observed that even with our stringent data curation pipeline, biases in the data generation pipeline and batch effects were observed that could potentially confound the interpretation of microbial community data.     

Design and Experimental Application of a Novel Non-Degenerate Universal Primer Set that Amplifies Prokaryotic 16S rRNA Genes with a Low Possibility to Amplify Eukaryotic rRNA Genes

The deep sequencing of 16S rRNA genes amplified by universal primers has revolutionized our understanding of microbial communities by allowing the characterization of the diversity of the uncultured majority. However, some universal primers also amplify eukaryotic rRNA genes, leading to a decrease in the efficiency of sequencing of prokaryotic 16S rRNA genes with possible mischaracterization of the diversity in the microbial community. In this study, we compared 16S rRNA gene sequences from genome-sequenced strains and identified candidates for non-degenerate universal primers that could be used for the amplification of prokaryotic 16S rRNA genes. The 50 identified candidates were investigated to calculate their coverage for prokaryotic and eukaryotic rRNA genes, including those from uncultured taxa and eukaryotic organelles, and a novel universal primer set, 342F-806R, covering many prokaryotic, but not eukaryotic, rRNA genes was identified. This primer set was validated by the amplification of 16S rRNA genes from a soil metagenomic sample and subsequent pyrosequencing using the Roche 454 platform. The same sample was also used for pyrosequencing of the amplicons by employing a commonly used primer set, 338F-533R, and for shotgun metagenomic sequencing using the Illumina platform. Our comparison of the taxonomic compositions inferred by the three sequencing experiments indicated that the non-degenerate 342F-806R primer set can characterize the taxonomic composition of the microbial community without substantial bias, and is highly expected to be applicable to the analysis of a wide variety of microbial communities.     

VITCOMIC2: visualization tool for the phylogenetic composition of microbial communities based on 16S rRNA gene amplicons and metagenomic shotgun sequencing

Background

The 16S rRNA gene-based amplicon sequencing analysis is widely used to determine the taxonomic composition of microbial communities. Once the taxonomic composition of each community is obtained, evolutionary relationships among taxa are inferred by a phylogenetic tree. Thus, the combined representation of taxonomic composition and phylogenetic relationships among taxa is a powerful method for understanding microbial community structure; however, applying phylogenetic tree-based representation with information on the abundance of thousands or more taxa in each community is a difficult task. For this purpose, we previously developed the tool VITCOMIC (VIsualization tool for Taxonomic COmpositions of MIcrobial Community), which is based on the genome-sequenced microbes’ phylogenetic information. Here, we introduce VITCOMIC2, which incorporates substantive improvements over VITCOMIC that were necessary to address several issues associated with 16S rRNA gene-based analysis of microbial communities.

Results

We developed VITCOMIC2 to provide (i) sequence identity searches against broad reference taxa including uncultured taxa; (ii) normalization of 16S rRNA gene copy number differences among taxa; (iii) rapid sequence identity searches by applying the graphics processing unit-based sequence identity search tool CLAST; (iv) accurate taxonomic composition inference and nearly full-length 16S rRNA gene sequence reconstructions for metagenomic shotgun sequencing; and (v) an interactive user interface for simultaneous representation of the taxonomic composition of microbial communities and phylogenetic relationships among taxa. We validated the accuracy of processes (ii) and (iv) by using metagenomic shotgun sequencing data from a mock microbial community.

Conclusions

The improvements incorporated into VITCOMIC2 enable users to acquire an intuitive understanding of microbial community composition based on the 16S rRNA gene sequence data obtained from both metagenomic shotgun and amplicon sequencing.

     

Genetic diversity of bacterial strains isolated from soils, contaminated with polycyclic aromatic hydrocarbons, by 16S rRNA gene sequencing and amplified fragment length polymorphism fingerprinting

In order to study microbial diversity in a polycyclic aromatic hydrocarbon-impacted soil, 14 bacterial strains were analyzed by 16S rRNA gene sequencing and amplified fragment length polymorphism (AFLP) analysis. Bacterial strains isolated from two different hydrocarbon-polluted sites were identified to the species level by 16S rRNA full-gene sequencing using MicroSeq 16S rRNA gene sequencing. Their genome was subsequently analyzed by high-resolution genotyping with AFLP analysis, in order to monitor species variability and to differentiate closely related strains. Cluster analysis based on AFLP fingerprinting showed intra-specific polymorphism, even among strains with 100% 16S rRNA gene sequence identity. The results show that AFLP is a powerful, highly reproducible and discriminatory tool for revealing genetic relationships in bacterial populations. The ability to differentiate and track related closely microbes is fundamental for studying structure and dynamics of microbial communities in contaminated ecosystems.     

Phylogenetic screening of a bacterial,metagenomic library using homing endonuclease restriction and marker insertion

Metagenomics provides access to the uncultured majority of the microbial world. The approaches employed in this field have, however, had limited success in linking functional genes to the taxonomic or phylogenetic origin of the organism they belong to. Here we present an efficient strategy to recover environmental DNA fragments that contain phylogenetic marker genes from metagenomic libraries. Our method involves the cleavage of 23S ribsosmal RNA (rRNA) genes within pooled library clones by the homing endonuclease I-CeuI followed by the insertion and selection of an antibiotic resistance cassette. This approach was applied to screen a library of 6500 fosmid clones derived from the microbial community associated with the sponge Cymbastela concentrica. Several fosmid clones were recovered after the screen and detailed phylogenetic and taxonomic assignment based on the rRNA gene showed that they belong to previously unknown organisms. In addition, compositional features of these fosmid clones were used to classify and taxonomically assign a dataset of environmental shotgun sequences. Our approach represents a valuable tool for the analysis of rapidly increasing, environmental DNA sequencing information.     

Identification and characterization of metagenomic fragments from tidal flat sediment

Phylogenetic surveys based on cultivation-independent methods have revealed that tidal flat sediments are environments with extensive microbial diversity. Since most of prokaryotes in nature cannot be easily cultivated under general laboratory conditions, our knowledge on prokaryotic dwellers in tidal flat sediment is mainly based on the analysis of metagenomes. Microbial community analysis based on the 16S rRNA gene and other phylogenetic markers has been widely used to provide important information on the role of microorganisms, but it is basically an indirect means, compared with direct sequencing of metagenomic DNAs. In this study, we applied a sequence-based metagenomic approach to characterize uncultivated prokaryotes from tidal flat sediment. Two large-insert genomic libraries based on fosmid were constructed from tidal flat metagenomic DNA. A survey based on end-sequencing of selected fosmid clones resulted in the identification of clones containing 274 bacterial and 16 archaeal homologs in which majority were of proteobacterial origins. Two fosmid clones containing large metagenomic DNAs were completely sequenced using the shotgun method. Both DNA inserts contained more than 20 genes encoding putative proteins which implied their ecological roles in tidal flat sediment. Phylogenetic analyses of evolutionary conserved proteins indicate that these clones are not closely related to known prokaryotes whose genome sequence is known, and genes in tidal flat may be subjected to extensive lateral gene transfer, notably between domains Bacteria and Archaea. This is the first report demonstrating that direct sequencing of metagenomic gene library is useful in underpinning the genetic makeup and functional roles of prokaryotes in tidal flat sediments.     

Prospecting for Novel Biocatalysts in a Soil Metagenome

The metagenomes of complex microbial communities are rich sources of novel biocatalysts. We exploited the metagenome of a mixed microbial population for isolation of more than 15 different genes encoding novel biocatalysts by using a combined cultivation and direct cloning strategy. A 16S rRNA sequence analysis revealed the presence of hitherto uncultured microbes closely related to the genera Pseudomonas, Agrobacterium, Xanthomonas, Microbulbifer, and Janthinobacterium. Total genomic DNA from this bacterial community was used to construct cosmid DNA libraries, which were functionally searched for novel enzymes of biotechnological value. Our searches in combination with cosmid sequencing resulted in identification of four clones encoding 12 putative agarase genes, most of which were organized in clusters consisting of two or three genes. Interestingly, nine of these agarase genes probably originated from gene duplications. Furthermore, we identified by DNA sequencing several other biocatalyst-encoding genes, including genes encoding a putative stereoselective amidase (amiA), two cellulases (gnuB and uvs080), an α-amylase (amyA), a 1,4-α-glucan branching enzyme (amyB), and two pectate lyases (pelA and uvs119). Also, a conserved cluster of two lipase genes was identified, which was linked to genes encoding a type I secretion system. The novel gene aguB was overexpressed in Escherichia coli, and the enzyme activities were determined. Finally, we describe more than 162 kb of DNA sequence that provides a strong platform for further characterization of this microbial consortium.     

Molecular characterizations of microbial communities fouling painted and unpainted concrete structures

This study reports the use of culture-independent and culture-dependent approaches to identify naturally occurring communities of Bacteria and Fungi fouling the surfaces of concrete structures with and without an acrylic paint coating in Georgia, USA. Genomic DNA was extracted from four different sites and PCR amplification of bacterial ribosomal RNA (16S rRNA) genes and the internal transcribed spacer (ITS) region of fungal rRNA genes was conducted. Bacterial and fungal community composition was determined by restriction analysis of amplified DNA of eight clone libraries and sequencing. Five bacterial phyla were observed, and representatives of the phylum Cyanobacteria and the classes Betaproteobacteria and Gammaproteobacteria dominated the bacterial clone libraries. The ITS region of rRNA gene sequences revealed the dominant phylotypes in the fungal clone libraries to be most closely related to Alternaria, Cladosporium, Epicoccum and Udeniomyces. The majority of these fungal genera could be cultured from the sites and successfully used to foul concrete in laboratory-based experiments. While the fungal sequences were most closely related to cultured isolates, the vast majority of bacterial sequences in the libraries were related to uncultured environmental clones. Results show phylogenetically distinct microbial populations occurring at the four sites.