Previous Article | Next Article ![]()
Applied and Environmental Microbiology, September 2008, p. 5305-5316, Vol. 74, No. 17
0099-2240/08/$08.00+0 doi:10.1128/AEM.01271-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

*
Katja Metfies,3,
Stephan Frickenhaus,1
Britta Knefelkamp,2
Karen H. Wiltshire,2 and
Linda K. Medlin4
Alfred Wegener Institute for Polar and Marine Research, Am Handelshafen 12, 27570 Bremerhaven, Germany,1 Alfred Wegener Institute for Polar and Marine Research (BAH), Kurpromenade, 27498 Helgoland, Germany,2 GKSS Research Centre, Max Planck Strasse 1, 21502 Geesthacht, Germany,3 Marine Biological Association of the United Kingdom, The Laboratory Citadel Hill, Plymouth PL1 2PB, United Kingdom4
Received 8 June 2008/ Accepted 4 July 2008
|
|
|---|
|
|
|---|
In the last decade, molecular methods have facilitated the investigation of the physiology, ecology, and distribution of this important part of the marine food web (5). In this respect, especially molecular probes have been shown to be very useful in terms of the detection and monitoring of microbial diversity, especially in the pico-sized fraction. Target genes for molecular probes are frequently the small and large subunit rRNA genes. They feature both highly conserved regions and variable regions that allow the development of probes at different taxonomic levels. The probes can be used in combination with a wide variety of hybridization-based methods (13, 14). Furthermore, the probes can be tailored as phylogenetic markers to identify phytoplankton at different taxonomic levels in a hierarchical fashion from classes down to species or strains using whole-cell and cell-free formats (21, 33, 34). The fluorescence in situ hybridization (FISH) of Prasinophyceae has shown great potential (3, 36), but the processing and quantitative analysis of FISH samples with fluorescent microscopy can be tedious, slow, and time demanding, as only one probe can be processed at a time because of the limited choice in fluorochromes. The application of molecular probes and FISH techniques in combination with flow cytometry have greatly increased the speed, detection, and accuracy of the monitoring of picoplanktonic communities (3, 34), but the limitation of fluorochrome choice still applies. However, the application of new methods, e.g., microarrays, for the detection of picoplanktonic eukaryotes with potential for high-throughput analysis of samples taken at close temporal and spatial scales can greatly facilitate the assessment of picoplanktonic species abundances, ecology, and physiology. DNA microarrays with species-specific probes have the unparalleled opportunity to detect thousands of targets in one experiment because the target is labeled, not the probes, and many probes can be used simultaneously. Although originally developed for gene expression applications, this innovative technique is being applied to species identification and has offered a promising experimental platform for microbial ecology. The most demanding challenge for the applicability of microarrays for microbial species identification is the high number of unknown environmental species and their corresponding DNA sequences that may result in unspecific signals or an oversight of species without a probe on the chip (11). However, the accuracy of the microarray can be enhanced by the application of hierarchical probes at different taxonomic levels because the detection of species is assessed by more than one positive probe signal (34). The development of a functional chip is an elaborate task because the set of probes on the microarray has to be developed such that all constituents work specifically under the same hybridization conditions (10, 34). The applicability of microarrays for the identification of phytoplankton has been demonstrated in various publications, e.g., for the detection and monitoring of harmful algae (12, 25) and marine microalgae at the class level (30, 33).
In this study, a phylochip for the detection and monitoring of the picoplanktonic prasinophytes was developed. Some of the probes were initially designed for FISH and showed specific hybridization signals, whereas other probes were designed for new groups of prasinophytes. The microarray contains 21 probes at different hierarchical levels. The specificity and discriminative potential of the probes were tested extensively under stringent conditions. Furthermore, the potential of the microarray was used to assess the occurrence of Prasinophyceae and the variation of its population composition over three annual cycles at the Helgoland Roads time series station.
|
|
|---|
![]() View larger version (51K): [in a new window] |
FIG. 1. Tree of prasinophyte diversity, showing which probes correspond to which clades. (Modified from reference 17 with permission from Elsevier.)
|
|
View this table: [in a new window] |
TABLE 1. Origin of algal cultures used in this study and their corresponding phylogenetic clades according to Guillou et al. (17)
|
|
View this table: [in a new window] |
TABLE 2. Probe sequences for the microarray
|
(iii) DNA extraction.
The template DNA was extracted from pure cultures with the DNeasy plant mini kit (Qiagen, Hilden, Germany) according to the manufacturer's protocol.
(iv) PCR amplification of the 18S rRNA gene.
For the probe specificity tests, the entire 18S gene was amplified from the target DNA (approximately 100 ng template DNA/100 µl PCR) with the universal PCR primers 1F (5'-AAC CTG GTT GAT CCT GCC AGT-3') and 1528R (5'-TGA TCC TTC TGC AGG TTC ACC TAC-3') without the polylinkers (29). The PCR protocol was 5 min at 94°C, 2 min at 94°C, 4 min at 54°C, and 2 min at 72°C for 29 cycles and an extension for 7 min at 72°C. All PCR experiments were carried out in a Mastercycler (Eppendorf, Hamburg, Germany). For the positive control in the microarray hybridization experiments, the yeast Saccharomyces cerevisiae was used. A 250-bp fragment of the TATA-box binding protein gene (TBP) of S. cerevisiae was amplified with the primers TBP-F (5'-ATG GCC GAT GAG GAA CGT TTA A-3') and TBP-R (5'-TTT TCA GAT CTA ACC TGC ACC C-3'). The TBP PCR amplification protocol was 5 min at 94°C, 1 min at 94°C, 1 min at 52°C, and 1 min at 72°C for 35 cycles and an extension for 10 min at 72°C. All PCR fragments were purified with the QIAquick PCR purification kit (Qiagen, Hilden, Germany) with modifications of the protocol to enhance the yield of PCR fragments. Step 8 (elution with elution buffer EB) was repeated with the same eluate. The DNA concentration was measured with a Nanodrop spectrophotometer (Peqlab, Erlangen, Germany).
(v) Biotin labeling of the purified PCR fragments.
In the hybridization experiments, the Biotin DecaLabel DNA labeling kit (Fermentas, St. Leon-Rot, Germany) was utilized. One labeling reaction contained 200 ng of PCR fragments and was incubated at 37°C overnight (17 to 20 h) to achieve the best biotin incorporation. The labeled PCR fragments were purified with the MinElute PCR purification kit (Qiagen, Hilden, Germany) with modifications of the manufacturer's protocol as above to enhance the yield. The DNA concentration was measured as above.
Field samples. (i) Sampling site and filtration.
The analyzed field samples were taken in the North Sea on the island of Helgoland as part of the Helgoland Roads time series (14, 39). Samples were taken every 2 weeks between February 2004 and December 2006. The Helgoland time series site is an anchorage area between the two islands of Helgoland, termed the Roads (54°11.3'N, 07°54.0'E). Picoplankton samples were obtained either by filtration of 1 to 1.5 liters of unfractionated seawater onto a 0.2-µm Isopore GTTP membrane filter (Millipore, Schwalbach, Germany) so that the entire phytoplankton community was collected on the filter or by filtration of
1.5 liters of fractionated seawater that was prefiltered through 10-µm and 3-µm Isopore TCTP membrane filters and finally collected onto a 0.2-µm Sterivex-GP filter (Millipore, Schwalbach, Germany). For DNA extraction of the latter set of samples, the 10-µm and the 0.2-µm filters were used; thus, the fraction between 10 µm and 3 µm was missing. The picoplankton are operationally defined as the <2-µm fraction, and thus, the prasinophytes were present in both sets of samples. All filters were immediately stored at –20°C.
(ii) DNA extraction.
DNA from the field samples was extracted as described previously for the extraction of genomic DNA from laboratory cultures.
(iii) PCR amplification of the 18S rRNA gene.
For sample analysis, a fragment of the 18S gene was amplified with the universal specific PCR primers 82F (5'-GAA ACT GCG AAT GAA TGG CTC-3') and 1528R (5'-TGA TCC TTC TGC AGG TTC ACC TAC-3') from the target DNA (approximately 100 ng template DNA/100 µl PCR) in a Mastercycler (Eppendorf, Hamburg, Germany). The PCR protocol of the 18S amplification and the TBP of Saccharomyces cerevisiae (positive control) and the purification and determination of the DNA concentration were identical to those for the specificity tests.
(iv) Biotin labeling of the purified PCR fragments.
Labeling of the PCR fragments was conducted as described above for the specificity test.
Chip design and hybridization. (i) Probe synthesis.
The molecular probes, including the positive and negative controls, were synthesized by Thermo Electron Corporation (Ulm, Germany) with a C6/MMT aminolink at the 5' end.
(ii) Microarray production.
Probes were spotted onto epoxy-coated Nexterion Slide E slides (Peqlab Biotechnologie GMBH, Erlangen, Germany) at a final concentration of 1 µM in 3x saline sodium citrate buffer (SSC [1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate]). We utilized the pin printer VersArray ChipWriter Pro (Bio-Rad Laboratories GmbH, München, Germany) and split pins (Point Technologies, Inc., CO). Subsequently, the slides were incubated at 60°C for 30 min in a Shake n Stack hybridization oven (Thermo Hybaid, Ulm, Germany). The microarrays were stored at –20°C. The chip contained four replicates of each probe in four independent blocks.
(iii) Standard hybridization protocol.
The hybridization solution was prepared with 1x hybridization buffer (1 M NaCl-10 mM Tris [pH 8]-0.005% Triton X-100-1 mg/ml bovine serum albumin-0.1 µg/µl herring sperm DNA) and the biotin-labeled PCR fragment at a final concentration of 11.25 ng/µl. The TBP fragment from S. cerevisiae was added as the positive control at a final concentration of 4.7 ng/µl. Blocking of the background noise was conducted by the prehybridization of the slides at 58°C for 1 h in a slide box with 50 ml 1x STT buffer (1 M NaCl-10 mM Tris [pH 8]-0.005% Triton X-100-1 mg/ml bovine serum albumin). Next, the slides were centrifuged, and the hybridization solution was incubated at 94°C for 5 min. A special coverslip, the Lifter Slip (Implen, München, Germany), was used for the hybridization. A volume of 30 µl hybridization solution was pipetted under the coverslip, and capillary action ensured even dispersal of the hybridization solution between the chip and the coverslip. The slide was placed in a humid chamber, which was constructed from a 50-ml Sarstedt tube filled with tissues moistened with hybridization solution. The hybridization was conducted at 58°C for 1 h; afterward, the slide was washed with 2x SSC-10 mM EDTA-0.05% sodium dodecyl sulfate and then 1x SSC-10 mM EDTA for 15 min and dried by centrifugation.
(ii) Staining.
The hybridized PCR fragments on the chip were stained with streptavidin-CY5 (Amersham Biosciences, Freiburg, Germany) in 1x hybridization buffer at a final concentration of 100 ng/ml. The chip was placed at room temperature for 30 min in a humid chamber and was washed afterward, twice with 2x SSC buffer for 5 min and once with 1x SSC buffer for 5 min, to remove excess staining moieties.
(iii) Scanning and quantification of microarrays.
The chip was scanned with a GenePix 4000B scanner (Molecular Devices, Sunnyvale, CA), and analysis of the obtained fluorescence signal intensities was done with GenePix 6.0 software (Molecular Devices, Sunnyvale, CA). A grid of circles was superimposed onto the scanned image to calculate the fluorescent signals and the surrounding background intensity.
(iv) Data preparation and analysis.
The hybridization data obtained from the array scanner were processed with the program PhylochipAnalyzer (32). This program implements the computation of signal-to-noise ratios (S/Ns) across the hierarchical probe set, while eliminating false positives. As we are analyzing time-dependent S/Ns for 21 probes, we must apply techniques to reduce the complexity for a structural analysis. This is done here by methods similarly applied in the analysis of differential gene expression profiles (9) as follows. Before that, all S/N values below 2, i.e., intensities below twice the background standard deviation (derived from the negative control signals), are set to zero. The samples with only zero values are removed from the data set. The resulting data are shown in Fig. 5A in color-coded form. Values above 26 are displayed in white, whereas zero values are in brown. These data are used for hierarchical clusterings (22) as displayed in Fig. 5B. Dissimilarities are computed from Euclidian distances of rows of data (for probe clustering) or columns of data (for time clustering). A complete-linkage clustering is computed from the dissimilarity matrices by standard procedures within the R package (http://www.r-project.org). In both figures, seasons are depicted in color-coded form. In Fig. 5B, this season data row is permuted, according to the order induced by time clustering (attached at the bottom).
![]() View larger version (55K): [in a new window] |
FIG. 5. (A) Color-coded representation of processed hybridization data. Top, S/Ns with the lower cutoff of 2 are in brown, the maximum value in white, and the medium values in red to yellow. Bottom, the color-coded seasons along the time axis are as follows: blue, winter; green, spring; yellow, summer; and red, autumn. The time-axis ticks are a combination of week and year (WWYY). (B) Clustered data. The data are horizontally (time) ordered by a hierarchical clustering (see text) shown at the bottom. The season color codes are ordered in the same way as described for panel A. The vertical (probe) order is given by hierarchical clustering displayed at the left.
|
|
|
|---|
![]() View larger version (29K): [in a new window] |
FIG. 2. Comparison of the signal intensities of probes Chlo01 and Chlo02 across all strains tested.
|
![]() View larger version (50K): [in a new window] |
FIG. 3. Three-dimensional plot of the signal intensities for all probes and all strains. (A) S/Ns for probes Pras03, 04, 05, 06, 01, and 07; (B) S/Ns for probes Pras08, 09, and 10; (C) S/Ns for Pras04 and genus and species probes in this clade.
|
(i) Clade VII—unnamed order containing the genus Picocystis.
Three probes were designed for this clade, each recognizing one of the three subclades in this undescribed order of prasinophytes. Given the genetic diversity among the three subclades, the groups should be described at least at the genus level. Probe Pras05 (clade VIIA) showed no cross-hybridization with nontarget strains and gave very strong specific signals for one target strain, unidentified Chlorophyceae sp. strain CCMP1205 (Fig. 3A). The other target strain identified from Fig. 1, unidentified coccoid RCC287, did not give a signal. Probe Pras06 for clade VIIB was tested only with nontarget species because no target strains were available from this clade (Fig. 3A). Here we detected only one weak cross-hybridization with Micromonas pusilla CCMP490. Probe Pras08 (clade VIIC) showed a signal for its target species Picocystis salinarum CCMP1897 (Fig. 3B). Several strains in other clades also cross-hybridized with this probe, but most signals were below the threshold of 2 (Fig. 4).
![]() View larger version (52K): [in a new window] |
FIG. 4. Summary of specificity testing for clade and genus probes.
|
(iii) Clade III, order Pseudoscourfieldiales, family Nephroselmidaceae.
The three probes designed for this clade are also family- and genus-level probes because only the genus Nephroselmis belongs to this clade. Pras10H showed one weak cross-hybridization to a nontarget strain, whereas Pras10F showed no unspecific binding with any closely related species (Fig. 3B). Pras10B showed specific results in hybridization with its target strain Nephroselmis pyriformis CCMP717 and minor cross-hybridizations with some other strains (Fig. 3B).
(iv) Clade II, Mamiellales.
Pras04 recognizes all members of this order (Fig. 3C). It showed clear signals with all seven strains tested. An unspecific signal to coccoid strain RCC287, a nontarget species, was observed. To resolve this, we sequenced this culture (data not shown) and found that it was a mixture of Bathycoccus prasinos and Micromonas pusilla (data not shown), and both contained the target sequence of probe Pras04. Thus, the culture used in this experiment is a mixture, and the signal can be regarded as correct. There are three described genera in this clade: Bathycoccus, Ostreococcus, and Micromonas. We tested genus-level probes for each of these genera. The probe for Bathycoccus, Bathy01, showed a signal for the specific cultures B. prasinos RCC378 and CCMP1898 and for unidentified coccoid RCC287 (Fig. 3C). The two Ostreococcus cultures, O. tauri RCC116 and RCC344, showed signal intensities for probe Ostreo01, but they were just above or at the threshold level, respectively (Fig. 3C). The probe Micro01 is specific for the genus Micromonas. Two cultures of M. pusilla, CCMP490 and CCMP1195, are target species for probe Micro01 and showed a signal for this probe (Fig. 3C). The signal for M. pusilla CCMP1195 was very low but above the threshold (Fig. 3C). The probes MicroA, B, and C are specific for different subclades within the Micromonas pusilla species complex (17). At least five morphologically indistinguishable different groups can be detected in this complex (43). The culture M. pusilla CCMP490, which belongs to the Micromonas clade A (17), gave a signal for MicroA and also for MicroC (Fig. 3C). One possible explanation for the mixed signals for M. pusilla CCMP490 could be that the culture is not clonal. Representatives of the different M. pusilla clades occur worldwide (43), and members of two clades could easily be present in one culture. The subclade affiliation of culture M. pusilla CCMP1195 is unknown, because it was not in the study of Guillou et al. (17). This culture showed only signals for probe MicroA (Fig. 3C). No signals were observed for MicroB.
(v) Clade I, Pyramimonadales.
The known diversity in clade I was covered by three specific probes: Pras09A1, Pras09A2, and Pras09D (Fig. 3B). These three probes recognize the described genera Pterosperma, Halosphaera, Cymbomonas, Prasinopapilla, and Pyramimonas, the last of which is paraphyletic. Specific hybridization results were observed for Pras09A1 and Pras09A2 when hybridized with their target species, and they showed no nontarget hybridization. For the third probe of this clade, probe Pras09D, no target strains were available in culture, but one weak cross-hybridization occurred with Prasinoderma coloniale CCMP1220; thus, it can be assumed to be specific.
(vi) Clade VI, order Prasinococcales.
This clade consists of two described genera, Prasinococcus and Prasinoderma, and one unidentified strain that shows sufficient genetic divergence to be another genus. Probe Pras03 for clade VI showed strong specific signals to three out of five hybridized target strains, one target strain showed a weaker signal, and another showed no signal at all (Fig. 3A). Weak cross-hybridizations with three strains from other clades (I, II, III, and V) were detected.
Field samples.
The DNA microarray was used to characterize the community composition of the Prasinophyceae at the Helgoland time series site in the North Sea site over a 3-year cycle between February 2004 and December 2006. In total, 61 samples were taken every 2 weeks within this time period and have been analyzed in this study. Overall, comparable S/Ns could be observed for all samples analyzed. However, the observation of variable signals for different targets of the same probe and inconsistencies within the hierarchies (e.g., for Chlo02 and Pras04) of the probes for some samples forced us to cut back our analysis to a rather general view. Hence, the overall signals per probe were counted as a first step in the analysis. The group with the highest number of positive signals consisted of the probes Chlo02, Pras04, Bathy01, Pras08, and Micro01. They were displayed in
40 to 70% of the samples with a positive signal (Fig. 5A). Furthermore, probes Chlo01, Pras04, and Bathy01 each showed one sample with a rather strong signal (above 26; Fig. 5A). The remaining probes showed a positive signal in only less than
25% of the samples. Some of the probes, such as Pras01 and Pras10F, did not show a positive signal at all. These data indicated that organisms belonging to the order Mamiellales, which are target species of the first group of probes, were the major contributors to the prasinophyte community at the Helgoland Roads sampling site. All other clades of the prasinophytes appeared to be minor components. The overall trend seems to be that within the Mamiellales, Bathycoccus made a larger contribution to the community composition than Micromonas. Similar results were detected with FISH (data not shown). Our data also indicated a seasonal variation in the community composition, because in the summers of the years 2004 and 2006 (approximately weeks 30 to 40), positive signals could be observed for Micromonas but not for Bathycoccus.
In order to assess the diversity within the prasinophytes, the complexity of the samples in terms of positive signals was addressed. It could be observed that the complexity (diversity) of the signals (columns) was strongly time dependent. Phases of higher complexity (diversity) for winter 2005/2006, late summer 2006, and winter 2006/2007 were observed.
Seasonal variation and succession were further assessed by grouping similar samples together. The data were reordered in time (horizontally) by the hierarchical clustering depicted as the attached tree below the matrix plot in Fig. 5B. The hierarchical approach sorted the data into blocks of similar signals. For example, at the right, a period between 1406 and 1206 showed similarity in the Bathy01 and Chlo01 signals. This phase is subdivided in two clades between 3205 and 4505, i.e., in the right clade, signals of probes Chlo02 and Pras04 appear. There were two samples of comparable high complexity (many signals), namely 3706 and 5006. Samples 0406, 0206, 0506, and 0806 are grouped at the left. They are persistently dominated by signals present in three groups: group I (Pras10B, Pras08, and Chlo02), group II (Pras09D, Pras05, and Pras03), and group III (Bathy01 and Pras04). These groups appeared as blocks in the matrix plot because of the vertical reordering, which is induced by the hierarchical clustering attached as a tree at the right in Fig. 5B. In contrast, the probes grouped in green showed much lower and fewer signals.
Despite these described similarities, overall, the clustering did not result in a clear picture that displayed a regular seasonal variation of species composition within the prasinophytes.
|
|
|---|
Specificity tests.
In the first part of this study, 21 phytoplankton probes at different hierarchical levels were tested and evaluated with 20 algal cultures in a microarray hybridization format. Part of the probe set was previously published in the context of other hybridization techniques (30, 33), and others were newly designed. The previously published probes had to be reevaluated on the microarray, because probes that work specifically in the context of other hybridization methods perform quite differently, or even not at all, in a microarray format (31). However, for the current probe set, we observed that the majority of the probes resulted in hybridization signals on the microarray. Furthermore, most of the specific probes were significantly stronger than the signals for nontarget species. The probes Pras08 and Pras10B, for clades VII and III, respectively, showed cross-hybridization to a significant number of nontarget species, e.g., all tested species belonging to the Mamiellales. Therefore, these probes may be considered unspecific. An attempt to design new, more specific probes for these clades was hampered by the diversity of the clade. Nevertheless, if the strengths of the true-positive signal and the false-positive signal are taken into consideration for the interpretation of hybridization signals, the probes could be used for the characterization of prasinophyte communities. The true-positive signal of the probes is around 5x higher than the false-positive signal. In addition, the signal of the cross-hybridization to species that belong to the Mamiellales is in the same range as the signal for the respective probes targeting this group (e.g., Pras04). Thus, if a chip displays positive signals for probes Pras04, Pras08, and Pras10B and the signal of Pras08 and Pras10B is in the same range as the signal for Pras04, it can be concluded that the signal is originating from Mamiellales. This indicates that an identification of different clades within the Prasinophyceae with the presented probe set is feasible with a standardized hybridization protocol.
In the course of the specificity testing, we observed large differences in the signal intensities of probes hybridized to different perfectly matching targets, e.g., for the Chlo01 and Chlo02 probes, positive hybridization signals could be observed for all target species. However, the signal intensity varied over a broad range. This has been observed frequently for probes that cover broad taxonomic groups (36). Five probes (Pras01, 04, 05, and 09A1 and 09A2) were specific for their target clades (V, II, VIIA, and I, respectively) and did not cross-hybridize with closely related species. The same results were observed for nearly all species and subclade probes of clade II (Bathy01, Micro01, A, B, and C). However, Ostreo01 displayed relatively low signal intensities. For a group of four probes (Pras06, 09D, and 10F and 10H for clades VIIB, II, and III, respectively), no target strains were available. They showed no significant cross-hybridization with all other prasinophyte cultures. These results indicate a high specificity for these probes. Thus, interpretation of the field sample analysis can be complicated because it is not clear if the probes result in signals at all if a respective target species is present. Specific hybridization signals for the target and cross-hybridization with closely related species were observed for probes Pras03, 07, 08, and 10B (clades VI, V, III, and VIIC, respectively).
Three cultures resulted in inconsistent signals because they did not hybridize to the corresponding perfectly matching probe. Two of them did not result in a signal at all, whereas the third one did bind to another probe (Pras04 instead of Pras05). Sequencing of the 18S rRNA gene of the targets revealed that this result was correct and that the target species did contain the probe sequence, but this does not explain why this target species was placed in another clade in the phylogenetic analysis, unless there were originally multiple isolates in this culture which have been sequentially replaced over time. This displays the potential and accuracy of the method. DNA microarrays could be suited to serve as a valuable tool to ensure the quality of a culture collection. The other two strains, which did not result in a hybridization signal, were correctly identified. In these cases, hybridization must have been hampered by other factors. The observation that most probes displayed relatively strong variation if they were hybridized to different target organisms underpins the widely accepted observation that in situ hybridization results for probes are not always consistent with their in silico predictions. Future research has to focus on the assessment of these discrepancies. Another weakness of the DNA microarrays are cross-hybridizations. They are almost impossible to avoid in a microarray hybridization format under stringent conditions (26). In this study, we evaluated probes that have been designed and optimized in respect to the hybridization temperature independently. The aim was to evaluate their performance under a single given hybridization protocol. For example, the unspecific signals for probe Pras08 could not be prevented under our hybridization conditions. It is likely that the probe may work specifically under more stringent hybridization conditions. However, it is not clear if all other probes on the microarray would keep their specificity and sensitivity under more stringent conditions. Therefore, we have to find ways to deal with cross hybridization, e.g., by the application of multiple probes for one target or hierarchical sets or by hybridization of the same sample at two different temperatures. Both approaches can prevent the misinterpretation of false-positive signals and will further improve the reliability of the microarray. Hierarchical probe sets consist of probes that target species at different taxonomic levels. In such a probe set, a signal at the species level should only be considered truly positive if all probes in the taxonomic hierarchy also show positive signals. The PhylochipAnalyzer software offers the possibility to examine all probes in a defined hierarchy and therefore represents a major advance in the data processing and interpretation of microarray experiments (32).
Field samples.
In the second part of this study, the specificity-tested DNA microarray was used to assess the community composition of the picoeukaryotic class Prasinophyceae. In total, 61 samples were analyzed that were taken every 2 weeks between February 2004 and December 2006. Overall, the results of the different samples displayed comparable hybridization signals. However, a few inconsistencies among the probes that targeted different taxonomic levels within the prasinophytes were observed. The probes from the beginning of 2004 did not show positive signals for Chlo02, but the probes at lower hierarchical levels, such as Pras04, did. It was stressed before that hierarchical probe sets should contribute to the accuracy of the results from a DNA microarray. Technically, all signals that were observed in the absence of Chlo02 should be considered negative. However, it is very unlikely that this is the case because prasinophytes have been identified to be major contributors to the picoplankton. It is also unlikely that there are no prasinophytes in these earlier samples. We assume that we have a technical problem here that is related either to the production of the microarrays or to the normalization of the data. The quality of the chip production was assessed by staining one chip from each batch with Sybr green, a fluorescent dye that is suited to stain DNA on the chips. However, because the DNA microarrays have been spotted with a split pin system, it cannot be excluded that the pin was blocked at some point during the spotting process of the respective probes. The normalization of the data poses another problem. As an internal marker that could be used for the normalization of the data in respect to variations related to the hybridization process is missing, an external control was added to the hybridization mixture. The data were normalized to this external control. Nevertheless, this approach is very susceptible to variations in the concentration of the control, e.g., because of handling.
A general analysis of the data generated with the microarrays indicates that the major constituents of the prasinophytes at the Helgoland Reed sampling site belong to the order Mamiellales. This is in accordance with a previous publication that reported the order Mamiellales to be dominant among the Chlorophyta in the English Channel (36). However, in this publication, it was described that among the Mamiellales, Micromonas pusilla was the dominant species. Our data for the Helgoland sampling site indicate that the genus Bathycoccus could be the major contributor to prasinophyte composition. The probe targeting this genus has more positive signals compared to Micromonas and Ostreococcus, which play only a minor role in the community composition of the prasinophytes at Helgoland. Furthermore, our data indicate a seasonality of the occurrence of Bathycoccus and Micromonas, because the probe for Bathycoccus did not display a signal in the summers of the years 2004 and 2006, whereas Micromonas was detected during these periods. Nevertheless, these findings should be confirmed by other methods, e.g., quantitative PCR, but they were comparable to counts observed by FISH methods (data not shown).
Clustering.
By the hierarchical clustering of the probes, similar time-dependent hybridization patterns can be identified. In contrast, the hierarchical clustering of the samples does not reveal a clear grouping when looking at the corresponding permutation of seasons, i.e., the time structure of the samples does not obviously correlate with seasons. This finding may be taken as an indication that there is no repetitive seasonal pattern of time course detectable within the hybridization patterns which, in contrast, one would expect for a community with strong interactions.
Conclusion.
A DNA microarray has been developed and assessed for the characterization of the community composition of Prasinophyceae. It is shown that the probes evaluated and tested in this study offer the potential to generate quickly a general impression of the community structure of the Prasinophyceae. However, some limitations of the technology had to be acknowledged, and the demand for improvement has to be pointed out. Our study represents one of the first applications of a low-density phylochip for the simultaneous detection and identification of different picoplanktonic species. Our analysis identified the genus Bathycoccus to be the major constituent of the prasinophyte community in the North Sea at the island of Helgoland time series site.
We thank Helga Mehl (Alfred Wegener Institute for Polar and Marine Research, Bremerhaven, Germany) for her assistance in maintaining cultures, DNA extraction, PCR preparation, and sequencing. Bank Beszteri (Alfred Wegener Institute for Polar and Marine Research, Bremerhaven, Germany) is acknowledged for valuable discussions on clustering methods. Furthermore, we thank Kristine Carstens, Gunnar Gerdts, Robert Niestoj, Silvia Peters, and Antje Wichels (Alfred Wegener Institute for Polar and Marine Research, Helgoland, Germany) for their assistance in sampling, providing samples and data, and DNA extraction.
Published ahead of print on 11 July 2008. ![]()
These authors contributed equally to this publication. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»