Abundance and Distribution of Ostreococcus sp. in the San Pedro Channel, California, as Revealed by Quantitative PCR

ABSTRACT Ostreococcus is a genus of widely distributed marine phytoplankton which are picoplanktonic in size (<2 μm) and capable of rapid growth. Although Ostreococcus has been detected around the world, little quantitative information exists on its contribution to planktonic communities. We designed and implemented a genus-specific TaqMan-based quantitative PCR (qPCR) assay to investigate the dynamics and ecology of Ostreococcus at the USC Microbial Observatory (eastern North Pacific). Samples were collected from 5 m and the deep chlorophyll maximum (DCM) between September 2000 and August 2002. Ostreococcus abundance at 5 m was generally <5.0 × 103 cells ml−1, with a maximum of 8.2 × 104 cells ml−1. Ostreococcus abundance was typically higher at the DCM, with a maximum of 3.2 × 105 cells ml−1. The vertical distribution of Ostreococcus was examined in March 2005 and compared to the distribution of phototrophic picoeukaryotes (PPE) measured by flow cytometry. The largest contribution to PPE abundance by Ostreococcus was ∼70% and occurred at 30 m, near the DCM. Despite its relatively low abundance, the depth-integrated standing stock of Ostreococcus in March 2005 was ∼30 mg C m−2. Our work provides a new technique for quantifying the abundance of Ostreococcus and demonstrates the seasonal dynamics of this genus and its contribution to picoeukaryote biomass at our coastal sampling station.

Picoeukaryotes (algae and protozoa of Ͻ2 m in diameter) from marine and freshwater ecosystems have increasingly become the focus of ecological, physiological, and genomic studies. The discovery of their widespread distribution, large genetic diversity, and periodically high abundance has caused a reassessment of their importance in microbial food webs. In particular, picoeukaryotic phototrophs contribute significantly to microbial biomass and total primary productivity (39,40,70). As a consequence, phototrophic picoeukaryotes (PPE) are an important component of the base of pelagic microbial food webs. Picoeukaryotes serve as prey for nanoplanktonic (2 to 20 m) phagotrophic protists (14), whose grazing activities provide a link to higher trophic levels (57) and a mechanism for nutrient remineralization (56).
Many picoeukaryotes have only recently been described in the scientific literature. These discoveries include both phototrophic (3,16,18,31) and heterotrophic (30) taxa. The global distribution of many picoeukaryotes remains to be determined. However, Ostreococcus isolates have been documented at numerous sites around the world, including the Mediterranean Sea (22,67), Long Island Sound, New York (48), the English Channel (32,54), the Arabian Sea (9), and the North Pacific Ocean (13, 70; S. Suda et al., unpublished data).
The discovery of many picoeukaryotes has been aided by advances in instrumentation (e.g., flow cytometry [FCM]), improved culturing techniques, and the development of cultureindependent molecular assays. Culture-independent molecular tools such as environmental cloning/sequencing of rRNA genes (49) and DNA fragment-based techniques such as denaturing gradient gel electrophoresis (DGGE) (45) and terminal restriction fragment length polymorphism analysis (41) have become particularly important for revealing the presence of previously unknown taxa (21,23,30,43) and documenting the geographic ranges of both newly and previously described taxa (42,53). In addition, the accumulation of DNA sequence information in public databases has facilitated the development of taxon-directed approaches to identify and count species of interest. Fluorescent in situ hybridization (FISH) and quantitative PCR (qPCR) have proven effective for the detection and quantification of specific microbial eukaryotes from complex natural assemblages (47,71).
The use of qPCR for targeting eukaryotic microbes in environmental samples has received increased attention in recent years. A number of studies employing quantitative genetic methods for the study of bloom-forming algae that pose a threat to the health of humans or marine organisms have recently emerged (5,8,15,28,29,50,55). These approaches have detection accuracies over a large dynamic range and have substantially shorter sample processing times than more traditional techniques. Increasingly, qPCR methods have been developed for carrying out ecological studies of cryptic marine microbes.
We designed and validated a qPCR assay specific to members of the genus Ostreococcus and applied this technique to samples collected at the University of Southern California (USC) Microbial Observatory site, located in coastal waters of the eastern North Pacific. Subsurface samples from 5 m and the deep chlorophyll maximum (DCM) were analyzed for a 2-year period (September 2000 to August 2002) to reveal Ostreococcus dynamics. Our estimates constitute some of the highest abundances ever reported for Ostreococcus (up to 3.2 ϫ 10 5 cells ml Ϫ1 ) at a coastal oceanic site. Additionally, Os-treococcus comprised a substantial proportion of the total number of phototrophic picoeukaryotes in March 2005. Overall, this minute prasinophyte was almost always present and occasionally abundant at our study site off the coast of Southern California.

MATERIALS AND METHODS
Collection and processing of environmental samples. Samples were collected at approximately monthly intervals from the San Pedro Ocean Time-Series Station (latitude, 33°33ЈN; longitude, 118°24ЈW) ( Fig. 1). Water samples were collected with Niskin bottles during conductivity-temperature-depth (CTD) casts from 5 m and from the DCM between September 2000 and August 2002. The depth of the DCM was determined by examination of real-time in situ fluorometry data.
Individual water samples were prescreened through a 47-mm inline filter apparatus (Pall Corp., East Hills, NY) containing 200-m Nitex fabric (Sefar Filtration Inc., Monterey Park, CA) via gravity filtration directly from Niskin sampling bottles into carboys. All sample collection and processing were con-ducted with acid-washed (5% HCl) lab equipment. Two liters of the Ͻ200-m filtrate from each depth was filtered onto a GF/F filter (Whatman Inc., Florham Park, NJ), using a gentle vacuum (Ͻ10 mm Hg), to collect planktonic biomass. Filters were loosely rolled and placed into 15-ml Falcon tubes (BD Biosciences, San Jose, CA). Two milliliters of lysis buffer (100 mM Tris [pH 8], 40 mM EDTA [pH 8], 100 mM NaCl, 1% sodium dodecyl sulfate) was added to each tube and frozen in liquid nitrogen. All samples were stored at Ϫ80°C until processing.
A vertical profile of Ostreococcus at 10 depths through the water column was obtained during our March 2005 time-series cruise. Samples were collected at 10-m intervals from near the surface to a depth of 80 m, with one additional sample at 100 m. The purpose of this cast was to examine the vertical distribution of Ostreococcus at a finer scale than was possible with our two-depth sample archive. Samples collected from the fine-scale vertical profile were prescreened through both 200-and 80-m Nitex fabric. Size fractionation was performed to investigate whether Ostreococcus was associated with particles in the 80-to 200-m size fraction.
Additional samples were collected from unfiltered seawater for analysis of the protistan community (in particular diatoms) by light microscopy. Seawater was sampled directly from Niskin bottles and preserved with Lugol's solution in amber glass collection jars (63). Lugol's solution-preserved samples were allowed to settle overnight onto counting chambers and examined by inverted light microscopy for microplankton enumeration. qPCR sample processing. Environmental samples and Ostreococcus culturebased standards were processed as lysates to minimize sample-to-sample variability due to losses associated with traditional DNA extraction procedures. Frozen samples and standards were thawed in a 70°C water bath prior to the addition of 200 l of 0.5-mm zirconia/silica beads (BioSpec Products, Bartlesville, OK). Cells were lysed by bead beating (vortexing at highest setting for 30 to 60 s) followed by heating in a 70°C water bath for 3 to 5 min. Two additional rounds of bead beating and heating were conducted to ensure complete cell lysis. Lysates were filtered through 0.22-m, low-binding cellulose acetate syringe-tip filters (Corning Life Sciences, Acton, MA) to remove cellular debris, filter remnants, and beads. Filtered lysates were collected in sterile cryotubes (Nalge Nunc, Rochester, NY) and frozen at Ϫ20°C. qPCR probe and primer design. A TaqMan (Roche Molecular Systems, Inc.)based qPCR assay (34,35) was developed to enumerate members of the genus Ostreococcus. An 18S rRNA gene probe (Ostreo-670F; Table 1) specific to positions 670 to 697 (relative to Ostreococcus tauri [GenBank accession no. YI5814]) (17) was designed by manually searching a ClustalX (66) alignment of prasinophyte 18S rRNA gene sequences retrieved from GenBank (6). Only two Ostreococcus sequences (Y15814 [17] and AB058376 [S. Suda et al., unpublished data]) were available in GenBank at the time of probe and primer design. A full-length Ostreococcus 18S rRNA gene sequence (GenBank accession no. DQ007077) from our study site aided in the development of our qPCR probe and primers. This clone was selected for full-length sequencing (ϳ6ϫ coverage) based on its prior identification as Ostreococcus by single-pass DNA sequencing (ϳ500 bp). The full-length Ostreococcus clone was 98.7% similar to Ostreococcus tauri (Y15814) across 1,738 common nucleotide positions revealed by pairwise sequence comparison in BioEdit, version 7.0.1 (33). The three Ostreococcus sequences described above were 100% identical in the selected probe region (see Table S1 in the supplemental material). This region displayed a 6-nucleotide gap (among other mismatches) between the C and the T at a position four nucleotides from the 3Ј end of the probe relative to all other prasinophytes in an alignment of Mamiellales sequences. A recent (January 2006) BLAST search (1) indicated that 11 additional Ostreococcus 18S rRNA gene sequences had been submitted to GenBank since the time of our probe and primer design (see Table S1 in the supplemental material). The TaqMan probe showed 100% sequence identity to all additional Ostreococcus sequences, with one exception, AY425313 (32), which displayed five mismatches compared to the   Table S1 in the supplemental material). Such mismatches would likely render the Mediterranean Ostreococcus strain represented by the AY425313 sequence undetectable. In addition to Ostreococcus sequences, the probe matched a number of "uncultured environmental" 18S rRNA gene clone sequences in GenBank, suggesting the widespread distribution of this genus. Not et al. (46) designed an Ostreococcus FISH probe targeting a similar region of the 18S rRNA gene to that targeted by our TaqMan probe to successfully detect Ostreococcus in environmental samples (46). Forward and reverse qPCR primers (Table 1) were generated with Beacon Designer 2.1 (Premier Biosoft International, Palo Alto, CA), which predicted a 187-bp PCR product (see Table S1 in the supplemental material). The forward primer (Ostreo-636F) was located at nucleotide positions 636 to 656, and the reverse primer (Ostreo-822R) was located at positions 802 to 822 relative to the Ostreococcus tauri sequence (YI5814) (see Table S1 in the supplemental material). The 3Ј end of the forward primer was designed to anneal as close to the 5Ј end of the Ostreococcus probe as possible to ensure rapid hydrolysis of the probe and nearly instantaneous fluorescence detection during the combined annealing/ extension phase (11). Additionally, the forward primer was designed for specificity to Ostreococcus, thus enhancing the stringency of detection and decreasing the probability of primer hybridization to nonspecific targets. The reverse primer was not unique to Ostreococcus but was specific to phytoplankton within the order Mamiellales. Four of the 11 new Ostreococcus sequences noted above displayed either 1-or 2-bp mismatches with our Ostreococcus primers (see Table  S1 in the supplemental material). It is conceivable that such mismatches would render these strains undetectable by our method; however, it should be noted that these sequences were recovered from the Atlantic Ocean and the Mediterranean Sea (32; M. Viprey and D. Vaulot, unpublished data) and, to our knowledge, have not been detected in the Pacific Ocean. The Ostreococcus primers and probe (Ostreo-670F) were synthesized by QIAGEN (Alameda, CA). The TaqMan probe included a 6-carboxyfluorescein fluorescent reporter molecule at the 5Ј end of the oligonucleotide and Black Hole Quencher 1 (Biosearch Technologies, Inc., Novato, CA) at the 3Ј end.
Calibration of Ostreococcus qPCR method. The qPCR method was calibrated by two methods, with plasmid DNA from an Ostreococcus 18S rRNA gene clone and with cell lysates from an Ostreococcus culture. The cloned 18S rRNA gene was used to test reaction conditions and assay sensitivity, while cell lysate reactions were the primary reference for calibrating the abundance of Ostreococcus in environmental samples. A 3.2 ng l Ϫ1 stock of the cloned Ostreococcus 18S rRNA gene was serially diluted to create a standard curve spanning 8 orders of magnitude. One-microliter aliquots of each dilution were used for each qPCR (see details below).
An actively growing culture of Ostreococcus sp. (MBIC-10636) was used to establish a standard curve for calibration of the qPCR assay by direct comparison of cell numbers to the qPCR threshold cycle (C T ). Flow cytometry indicated an abundance of 3.2 ϫ 10 7 (Ϯ0.1 ϫ 10 7 ) cells ml Ϫ1 in a stock culture of Ostreococcus. Standards were prepared by adding 100, 10, and 1 ml of undiluted Ostreococcus culture to triplicate 1-liter aliquots of filtered seawater (FSW) for the three highest concentration standards. The three lowest concentration standards were constructed by adding 1 ml of serially diluted stock culture (10 Ϫ1 , 10 Ϫ2 , and 10 Ϫ3 ) to triplicate 1-liter aliquots of FSW. Ostreococcus-spiked FSW aliquots were filtered onto 47-mm GF/F filters and prepared as lysates as described above. The preparation of standards by spiking cultured cells into 1 liter of FSW simulated the processing procedure applied to natural samples. qPCR optimization and reaction conditions. Optimization reaction mixtures and plasmid standards for qPCR were prepared with 1 l of the cloned Ostreococcus 18S rRNA gene (3.2 ng l Ϫ1 ) plus 49 l of reagent master mix (see below) and run on a Bio-Rad iCycler machine (Hercules, CA). Initial testing of thermal protocols was performed with a Mg 2ϩ concentration of 2.5 mM and an annealing temperature of 55°C. Two-and three-step thermal protocols generated similar amounts of the 187-bp product, as determined by the intensity of SYBR gold (Molecular Probes, Eugene, OR)-stained DNA bands on a 1.2% agarose gel. Further optimization tests were performed using a two-step amplification protocol (described below) typical of most TaqMan assays.
The annealing/extension temperature was optimized with the gradient feature of the iCycler across a temperature range of 60 to 70°C. The 60°C reaction reached the C T ϳ2 cycles earlier than the reaction at 70°C (13.6 versus 15.4) and attained higher relative fluorescence unit (RFU) values at the plateau phase of amplification (Ͼ1,100 RFU at 60°C versus ϳ500 RFU at 70°C). This result was expected since the TaqMan probe is required to anneal to the template before the primers and to remain hybridized until hydrolysis by Taq polymerase, a condition which is favored at lower annealing/extension temperatures (10).
The optimal Mg 2ϩ concentration was determined by testing reactions prepared with 2.5 to 5.0 mM Mg 2ϩ . Average C T values decreased slightly from 12.7 to 11.4 with increasing Mg 2ϩ concentrations; however, there was no significant difference between the average C T values for reactions at 4.5 and 5.0 mM Mg 2ϩ (P Ͼ 0.10). DNA amplification reached higher RFU values with increasing Mg 2ϩ concentrations (1,100 RFU at 2.5 mM versus ϳ2,000 RFU for 4.5 and 5.0 mM). The lowest concentration of Mg 2ϩ that produced maximum RFU values and the earliest C T was 4.5 mM and was therefore selected as the optimal concentration for subsequent reactions. Higher RFU values enhance the ability to distinguish amplification over baseline fluorescence.
The Ostreococcus primer set was tested with Bio-Rad's iQ SYBR green supermix to determine the specificity of amplification, as indicated by postamplification melting-curve analysis. Melting-curve analysis indicated a single PCR product. Secondary PCR products and primer dimers were not detected in the melting-curve profile (see Fig. S1 in the supplemental material).
Sample and cell standard lysates were diluted 1:100 with sterile Milli-Q water to minimize PCR inhibition due to the presence of lysis buffer and residual cellular contents. The optimal volume of diluted lysate for use in qPCR assays was determined because increasing the lysate volume increased the reaction sensitivity only to a certain point. A lysate for this test was created by spiking 50 ml of a dense Ostreococcus culture (ϳ2.0 ϫ 10 7 cells ml Ϫ1 ) into 1 liter of FSW and processing the sample as described above. Various volumes of 1:100 diluted lysate were tested in 50-l reaction mixtures to determine the appropriate lysate volume (see Fig. S2 in the supplemental material). Lysate volumes of Ͼ10 l had no further effect on the reduction of qPCR C T values, as evaluated by unpaired t tests assuming equal variance. Specifically, the average C T values decreased significantly (P Ͻ 0.05) with increasing volumes until 15 l, which did not have a significantly different C T value from the average value obtained for 10 l (P ϭ 0.07).
The final reaction conditions in 50-l volumes included 500 nM of each primer and 250 nM of TaqMan probe. Reaction mixtures also contained the following reagents: 250 M of each deoxynucleoside triphosphate, 4.5 mM Mg 2ϩ , 1ϫ buffer B, and 2.5 units of Taq polymerase in buffer B (Promega, Madison, WI). Bovine serum albumin (A7030; Sigma) was added to each reaction to a final concentration of 300 ng l Ϫ1 to enhance the PCR (37). All reagents, samples, and standards were prepared on ice prior to thermal cycling. Ten microliters of each lysate was loaded into a well of a 96-well PCR plate, followed by the addition of 40 l of qPCR master mix. Plates were sealed with optically clear tape (Bio-Rad), centrifuged briefly to remove bubbles, and transferred to the iCycler machine.
Cell-based standards were comprised of triplicate independent lysates for each of six Ostreococcus concentrations and were run with every set of unknown samples to facilitate the conversion of C T values of unknowns to equivalent Ostreococcus abundances. Three individual reactions were set up for each sample to control for reaction variability. Samples and standards were cycled with the following thermal protocol: 95°C for 90 s (1ϫ), 70 cycles of 95°C for 15 s and 60°C for 30 s, and a hold at 4°C. Real-time data were collected during the annealing/extension step.
Cross-reactivity testing with protistan cultures. Fifty-two clonal cultures of marine and freshwater protists were grown to high abundance in 250-ml sterile culture flasks for Ostreococcus qPCR cross-reactivity testing. Marine phytoplankton cultures were grown in K medium (36) modified with 36 M PO 4 2Ϫ , while freshwater phytoplankton cultures were grown in DY-IV medium (2). Heterotrophic protists were grown in sterile seawater or freshwater amended with sterile yeast extract (final concentration, 0.005%) plus several rice grains. Abundances of cultured protists were checked by light microscopy after 1 to 2 weeks of growth to confirm a minimum abundance of ϳ1.0 ϫ 10 3 cells ml Ϫ1 . Cell cultures were filtered onto Whatman GF/F filters, lysed, and stored at Ϫ20°C. Selected culture lysates were amplified with universal eukaryote-specific PCR primers to ensure the presence of the 18S rRNA gene prior to cross-reactivity testing with Ostreococcus primers and probe. Lysates for qPCR testing were diluted 1:100 with sterile water, and 10 l of diluted lysate was used in duplicate 50-l qPCR assays. An Ostreococcus lysate served as the positive control.
Forty-eight of the 52 lysates returned negative results for duplicate qPCR amplifications (see Table S2 in the supplemental material). Two of the 52 lysates (Asterionellopsis glacialis and Paraphysomonas vestita) were positive in one of two reactions, while both duplicates of Aureoumbra lagunensis were positive. These nonspecific amplifications displayed much lower RFU values than the Ostreococcus positive control and occurred late in the thermal protocol (C T Ն 50). The three "positive" lysates were tested a second time and returned negative results for Ostreococcus cross-reactivity. The initial amplifications of nontarget taxa were attributed to aerosol contamination from the cloned Ostreococcus 18S rRNA gene, which initially was run as a secondary positive control but was omitted from later runs.
Comparison of picoeukaryotes by flow cytometry. Environmental samples were collected in triplicate from the March 2005 vertical profile for analysis by flow cytometry (described above). Samples were prescreened through 80-m Nitex fabric, preserved with formalin (final concentration, 1%), and frozen in liquid nitrogen. Analysis by FCM provided the putative abundance of Ostreococcus as a subset of total PPE.
An actively growing culture of Ostreococcus (MBIC-10636) was used to calibrate the flow cytometer for environmental detection of this genus. Ostreococcus cells were grown for approximately 1 week at 20°C on a 12-h light/dark cycle with an average irradiance of 165 microeinsteins m Ϫ2 s Ϫ1 (in water). FCM gate parameters for counting Ostreococcus were established with a four-color, duallaser FACSCalibur flow cytometer (BD Biosciences, San Jose, CA). Cultured cells formed a discrete FCM region based on forward-angle light scatter and red fluorescence (FL3 channel), detection parameters that have been used in previous studies (16,18,48). Micromonas pusilla (CCMP-487) and Micromonas sp. (a local isolate) were also analyzed by FCM to investigate potential overlap with the Ostreococcus FCM gate parameters. Micromonas and Ostreococcus are similarly sized prasinophytes that co-occur at the USC Microbial Observatory.
Nucleotide sequence accession number. An environmental Ostreococcus sp. clone from the USC Microbial Observatory (North Pacific Ocean) was sequenced in both directions to yield a full-length 18S rRNA gene construct. The sequence was deposited in GenBank and assigned the accession number DQ007077.

RESULTS
Calibration of qPCR. Initial testing of the TaqMan reaction sensitivity with the 18S rRNA gene cloned from Ostreococcus indicated linearity across 8 orders of magnitude (r 2 ϭ 0.99) ( Fig. 2A). This range corresponded to 6.5 ϫ 10 0 to 6.5 ϫ 10 8 copies of the Ostreococcus 18S rRNA gene per reaction. Standardization of the qPCR assay with the cloned 18S rRNA gene indicated a reaction sensitivity of ϳ10 copies for highly purified templates like plasmid DNA containing cloned inserts.
PCR analysis of lysates from known numbers of cultured Ostreococcus cells was used to convert C T values for unknown environmental samples directly into cell abundances. Cell lysate standards prepared from an Ostreococcus culture were analyzed on each plate with environmental samples. The cellbased standard curve displayed linearity over 5 orders of magnitude, ranging from the genomic equivalents of 16 to 160,000 cells per reaction (r 2 ϭ 0.99) (Fig. 2B). The cell-based standard curve depicted in Fig. 2B represents an average of three standard curves analyzed on three different days using the same set of lysates (stored frozen between runs). Error bars represent standard deviations of pooled determinations. This dynamic range bracketed the upper range of Ostreococcus abundance in all (diluted 1:100) environmental samples. In practice, the lowest reliably detected standard (equivalent to 16 cells per reaction) translated into the detection of 320,000 cells in 2 ml of undiluted lysate (e.g., 16 cells in 10 l of 1:100 diluted lysate, equivalent to 160 cells l Ϫ1 of undiluted lysate, ϫ 2,000 l lysis buffer ϭ 320,000 cells). Thus, the limit of detection was approximately 100 Ostreococcus cells ml Ϫ1 in natural samples for typical sample collection volumes (2 to 4 liters).
Vertical profile of Ostreococcus. Profiles of chlorophyll a, oxygen, and temperature in March 2005 revealed a mixed surface layer down to a depth of 20 m at our study site (Fig. 3). The DCM at ϳ35 m corresponded to the top of the oxycline and the middle of the main water column thermocline. This water column profile was fairly typical of profiles observed at our sampling site throughout much of the year (data not shown). The vertical profile of total PPE abundance by FCM reflected the profiles of the physical and chemical measure- below the depth of the DCM (Fig. 4). PPE abundances decreased to ϳ6.0 ϫ 10 3 cells ml Ϫ1 at 50 to 60 m and dropped precipitously to Ͻ200 cells ml Ϫ1 at 80 m and below.
Ostreococcus abundance in the Ͻ80-m size fraction was determined by qPCR and FCM and compared to the total PPE abundance (Fig. 4). Ostreococcus abundances determined by qPCR were fairly uniform from the surface to 20 m, ranging from 1.0 ϫ 10 3 to 2.0 ϫ 10 3 cells ml Ϫ1 , while a distinct subsurface maximum of 5.5 ϫ 10 3 cells ml Ϫ1 was observed at 30 m, just above the DCM, at the top of the primary water column thermocline ( Fig. 3 and 4). Ostreococcus abundances declined to Ͻ100 cells ml Ϫ1 at 60 m and below. Overall, the qPCRbased Ostreococcus abundance estimates followed the trends in the FCM-based PPE profile, with qPCR estimates generally constituting 13 to 26% of the total PPE abundance at particular depths. A notable exception occurred at 30 m, where the qPCR estimate of Ostreococcus comprised nearly 70% of the total PPE abundance, indicating that Ostreococcus can dominate PPE assemblages. Ostreococcus abundances in the 80-to 200-m size fraction were generally insignificant (Ͻ100 cells ml Ϫ1 ) at all depths tested (1 to 70 m).
In addition to measurement by qPCR, the putative abundance of Ostreococcus was measured by FCM, using gate parameters established with an Ostreococcus culture. Cytometry-based estimates were Ͻ200 cells ml Ϫ1 for depths down to 10 m and then increased to a broad peak of ϳ5.0 ϫ 10 3 cells ml Ϫ1 between 30 and 40 m before decreasing to Ͻ200 cells ml Ϫ1 at 70 m and below. Ostreococcus abundances determined by FCM differed from qPCR-based estimates in comparisons for 7 of the 10 depths (P Ͻ 0.05). Cytometry-based abundance estimates of Ostreococcus were lower than qPCR estimates for surface water samples and higher than qPCR estimates for depths below 30 m. However, cytometry and qPCR estimates were not significantly different at 10, 20, and 30 m (P ϭ 0.08, 0.53, and 0.10, respectively; ␣ ϭ 0.05).
In an effort to reconcile differences between qPCR estimates of Ostreococcus and putative abundances determined by FCM, the FCM characteristics of two Micromonas cultures were compared to those of a single Ostreococcus culture to examine the potential contribution of other small prasinophytes to the cytometric signature of Ostreococcus. Approximately 30% of the cytogram events detected with the Micromonas culture overlapped with the cytogram for Ostreococcus (Fig. 5). These results indicated a substantial overlap in Ostreococcus and Micromonas cytograms for the tested cultures. This test indicated the inability to accurately distinguish these genera by FCM and further illustrated the need for taxon-specific molecular approaches to establish species-specific abundances and distributions. Abundance and distribution of Ostreococcus through time. The depth of the DCM at the USC Microbial Observatory varied between 14 and 50 m during the 2-year study period (Fig. 6A). The shallowest depth of the DCM occurred during winter to early spring of 2001, while the deepest DCM was observed in November 2000. The overall shoaling or deepening of the DCM appeared to be more dependent on episodic events (vertical and horizontal advection) than seasonal trends. The 23-month average depth of the DCM was 32 m and varied within a relatively narrow range (Ϯ8 m [standard deviation]), which provided a stable water column feature to sample throughout the study. At times, the study site displayed relatively "openocean" characteristics (e.g., low nutrient concentrations, low planktonic biomass in net tows, and high light transmission) (data not shown), and at other times, the waters were much more coastal in nature, being visibly greener and more characteristic of nearshore (Ͻ1 km from the coast) water masses and phytoplankton communities.
Depth-integrated profiles of nitrate and phosphate (courtesy of the USC Wrigley Institute for Environmental Studies) were examined for seasonal trends to help explain the observed patterns of Ostreococcus abundance (Fig. 6B). Nutrients were expressed as depth-integrated values over the upper 40 m of the water column based on discrete measurements at 1, 10, 20, 30, and 40 m. Integrated nitrate and phosphate values were positively correlated (P Ͻ 0.01) at the USC Microbial Observatory site due to the dominance of coastal upwelling as the primary influence on nutrient delivery at this location. Two major nutrient upwelling events between March and June 2001 and February and June 2002 were revealed by the profiles (Fig. 6B).
Ostreococcus abundances determined by qPCR indicated that this genus was nearly ubiquitous at our study site, as it was detected at one or both depths on all dates except late October 2001 (Fig. 7). Month-to-month abundances of Ostreococcus during this study period were highly variable, ranging over 4 orders of magnitude at 5 m (from Ͻ100 cells ml Ϫ1 to 8.2 ϫ 10 4 cells ml Ϫ1 ; Fig. 7A) and over 5 orders of magnitude at the DCM (from Ͻ100 cells ml Ϫ1 to 3.  (Fig. 7). Two prolonged periods of low Ostreococcus abundance (Ͻ1.0 ϫ 10 3 cells ml Ϫ1 ) at 5 m were apparent, between September 2000 and January 2001 and between July and November 2001. In contrast, abundances of Ostreococcus of Ͻ1.0 ϫ 10 3 cells ml Ϫ1 at the DCM were detected on only three dates (Fig. 7B). There was no significant correlation between Ostreococcus abundance and the depth of the DCM (P Ͼ 0.05; data not shown). The average DCM depth on dates with the 11 highest Ostreococcus abundances and the 12 lowest abundances was 32 m.
A major bloom of Ostreococcus was detected at both 5 m and the DCM (27 m) in May 2001 (Fig. 7). Ostreococcus abundance increased from 3.8 ϫ 10 3 to 8.2 ϫ 10 4 cells ml Ϫ1 at 5 m and from 1.9 ϫ 10 3 to 3.2 ϫ 10 5 cells ml Ϫ1 at the DCM between April and May 2001. The Ostreococcus bloom was preceded by high abundances of diatoms in April 2001 (ϳ200 cells ml Ϫ1 at 5 m and ϳ300 cells ml Ϫ1 at the DCM). The peak in diatom abundance in April 2001 corresponded to high nitrate and phosphate concentrations in the upper 40 m of the water column (Fig. 6B). Abundances of Ostreococcus were among the lowest values detected during peak diatom blooms (Fig. 7). Nutrient upwelling events of similar magnitude occurred in April-May 2001 and May-June 2002 (Fig. 6B). Diatom abundances at 5 m and the DCM increased near the end of the study period, but Ostreococcus abundances increased only modestly at the DCM during this time frame (Fig. 7A and B).

DISCUSSION
Quantitative PCR for ecological research. Real-time qPCR analysis is increasingly the method of choice for quantifying abundances of protistan taxa in environmental samples (5,8,15,28,44,50,69,71). This approach combines extreme sensitivity with the ability to process large numbers of samples rapidly. It is conceivable that global surveys of particular taxa could be completed in a matter of days to help address the question, "Is everything everywhere?" (24)(25)(26). Although most qPCR assays have approximately the same sensitivity (detection of several to tens of target copies), we developed a TaqMan assay because the TaqMan probe provides another level of reaction specificity compared to primer-only SYBR green assays.
Several methods for relating C T values from environmental samples to gene copy numbers or cell abundance have been applied in qPCR-based studies of protistan taxa. Plasmid DNAs containing cloned target sequences have been used in a number of studies to create standard curves (28,69,71), and this is presently the only available method for calibrating qPCR detection of uncultured taxa (65). Calibration with plasmid DNA is advantageous in that exact numbers of target genes can be calculated by measuring the concentration of a DNA standard. However, a notable disadvantage when using plasmid standards is that environmental samples should be fully extracted and purified to ensure that amplification is directly comparable to the amplification of purified plasmid DNA. Complete extraction and purification may lead to high variability among replicate extracts and to spurious results due to inconsistent DNA recovery. The alternative to plasmidbased calibration involves calibration with a cultured target organism. Variations on cell-based standardizations include dilution of a single genomic DNA extract from a known number of cells (5), dilution of a cell culture of known concentration followed by multiple DNA extractions (8,15), spiking specific volumes of target organisms directly onto filters for collection and extraction (50), and spiking cultured cells into filtered seawater aliquots for subsequent collection by filtration and preparation as lysates (44). Cell-based calibration approaches (when possible) are the most direct means of relating C T values from qPCR assays to cell abundance. Also, the use of lysates prepared simultaneously from both standards and unknowns minimizes variability due to slight variances in DNA extraction efficiency and/or reagent quality.
Several caveats should be considered before the application of cell-based qPCR standardizations. Specifically, cell constituents released by lysis and the lysis buffer itself can act as PCR inhibitors leading to reductions in reaction sensitivity or accuracy. We minimized inhibition problems in our TaqMan assay by diluting all lysates 1:100 with sterile Milli-Q water and optimizing the volume of diluted lysate used in 50-l reaction mixtures (see Fig. S2 in the supplemental material). Cell-based qPCR calibrations can also be affected by the accuracy of microscopy-or FCM-based cell counts, thereby requiring replicate determinations of cell abundances in cultures. We quantified triplicate samples of live Ostreococcus cells by FCM immediately before dilutions of stock culture and subsequent collection by filtration to ensure the accuracy of standard curves. Finally, variability in the rRNA gene copy number per cell may affect the interpretation of cell-based calibrations. Recent evidence suggests that the copy number of 18S rRNA genes in Ostreococcus is relatively small, ranging from two to four copies per cell (71). However, there have been no comprehensive studies examining the variability in 18S rRNA gene copy number in Ostreococcus cells during different stages of growth. We did not attempt to control for variability in the 18S rRNA gene copy number of our cultured cells but assumed that a culture in nonsynchronous, exponential growth would provide a reasonable average copy number per cell. Furthermore, we assumed that individuals in natural populations of Ostreococcus expressed a similar range of 18S rRNA gene copy numbers as cultured cells.
Our Ostreococcus TaqMan assay was performed for 70 cycles of amplification to maximize the sensitivity of the assay. The small genome size of Ostreococcus (20,52), its low 18S rRNA gene copy number (71), and the low concentration of DNA in lysates were factors contributing to the requirement of a large number of amplification cycles. Most PCR assays are carried out using 35 or fewer amplification cycles to prevent the introduction of PCR bias to the final pool of amplicons (with overrepresentation of particular taxa) when amplifying environmental samples for cloning and sequencing (64) and to avoid the amplification of nonspecific targets, primer dimers, or contaminant DNA. Amplification of nonspecific 18S rRNA gene fragments and the formation of primer dimers were not ap- 2502 COUNTWAY AND CARON APPL. ENVIRON. MICROBIOL. parent in our qPCR assays because of the high specificity of our probe and primers for the target region (see Fig. S1 in the supplemental material). Additionally, interference of Ostreococcus-specific rRNA gene amplification by the DNAs of other protistan taxa in environmental lysates was not likely a significant factor given the results of our taxonomically diverse cross-reactivity tests (see Table S2 in the supplemental material). PCR contamination was adequately controlled, as indicated by the results of our negative controls (Milli-Q water), which remained below the amplification threshold for all sample runs (data not shown). Contamination was primarily controlled by isolating preparation steps during the qPCR setup and avoiding the use of plasmid DNAs containing 18S rRNA gene inserts in the vicinity of qPCR lysate preparation or reaction setup. Molecular detection has become an essential tool for documenting the presence and abundance of Ostreococcus because of the diminutive size of species within this genus and the lack of definitive morphological features. A variety of molecular approaches have begun to reveal its ecological dynamics and global distribution. Ostreococcus has been quantified by 18S rRNA probes (47) and detected in environmental clone libraries from a variety of marine ecosystems (7,13,32,46,54). Most recently, qPCR methods including SYBR green (71) and TaqMan (this study) have provided quantitative estimates of Ostreococcus abundances in environmental samples collected from time-series stations. Ostreococcus appears to be a significant component of PPE assemblages in some of the ecosystems where it has been detected. For example, Ostreococcus constituted ϳ10% of the clones in an 18S rRNA gene library from the USC Microbial Observatory (13).
Vertical and temporal distribution of Ostreococcus. Ostreococcus has been isolated and cultured previously from depths across the euphotic zone (68); however, our qPCR-based vertical profile from March 2005 constitutes the first discrete profile of Ostreococcus abundance through the euphotic zone at an oceanic station. This profile revealed the highest abundance of Ostreococcus near the depth of the DCM (Fig. 3 and  4). Ostreococcus comprised a substantial fraction of the PPE abundance across the euphotic zone and dominated the picoeukaryotic assemblage near the DCM. The maximum qPCRbased abundance of Ostreococcus in March 2005 reached 5.5 ϫ 10 3 cells ml Ϫ1 (Fig. 4), which was low compared to abundance estimates from time-series samples ( Fig. 7) but higher than Ostreococcus abundances from other locations (46,71). The vertical distribution of Ostreococcus in the euphotic zone (Fig.  4, open circles) is typical for phototrophic picoeukaryotes (12,39) and showed a strong correspondence to the physical structure of the water column (Fig. 3).
Ostreococcus abundances estimated by qPCR and FCM were significantly different for some depths but not significantly different for other depths (see Results) (Fig. 4). Samples for which FCM-based estimates were greater than qPCRbased estimates (e.g., samples from 40 to 70 m) were likely the result of non-Ostreococcus detections in the Ostreococcus FCM gate. The prasinophyte Micromonas was identified as a potential cause of discrepancies between qPCR-and FCM-based measurements. Micromonas and Ostreococcus regularly cooccur at other coastal sites (46,71). The overlap of Micromonas and Ostreococcus detection by FCM (ϳ30%) may have lead to an overestimation of the Ostreococcus abundance by FCM. Additionally, the lower qPCR estimates could have been due to the probe-primer mismatches noted above; however, this is extremely unlikely since none of the other strains have been detected in extensive clone libraries from the same site (P. D. Countway and D. A. Caron, unpublished data). FCMbased estimates of Ostreococcus abundance at the two shallowest depths in the vertical profile were lower than the corresponding qPCR estimates, which was possibly related to differences in the chlorophyll contents of near-surface Ostreococcus assemblages. Ostreococcus cells collected from near-surface assemblages may have been high-light-adapted ecotypes (53) with less chlorophyll per cell than the Ostreococcus culture used for FCM calibrations.
Monthly sample collections at our study site indicated that Ostreococcus abundances varied by several orders of magnitude at both 5 m and the DCM, with higher abundances typically occurring at the DCM (Fig. 7). To our knowledge, this is the first study to reveal the abundance and distribution of Ostreococcus across multiple depths for a multiyear time series. Prior investigations have reported abundances of Ostreococcus in the upper 1 to 2 m of the water column, where Ostreococcus was a lower percentage of total picoeukaryotes (46,71) than that observed in the present study. Not et al. (46) conducted a survey of prasinophytes at a time-series station in the English Channel, using FISH-tyramide signal amplification (47). Picoeukaryotes belonging to the order Mamiellales (Ͻ1.0 ϫ 10 4 cells ml Ϫ1 ) were detected in low abundance, with only 3% of the count attributed to Ostreococcus (46). Zhu et al. (71) reported the relative abundances of Ostreococcus from a timeseries study in the Mediterranean Sea. Their qPCR-based Ostreococcus detections never accounted for more than 1% of qPCR-based PPE detections (71). These low contributions contrast with our findings that Ostreococcus may constitute a significant fraction of total phototrophic picoeukaryotes, especially near the DCM (Fig. 4). However, Ostreococcus comprised a relatively low percentage of the total PPE assemblage at the shallowest depths of our study site (Fig. 4), which is consistent with previous reports. These results lead us to speculate that previous studies from other sites may have underestimated the importance of Ostreococcus, which appears to exhibit subsurface abundance maxima.
High concentrations of phototrophic picoeukaryotes (Ͼ1 ϫ 10 5 cells ml Ϫ1 ), as determined by flow cytometry, have previously been reported from nearshore locations along the Mediterranean coast of France (18,67) and Long Island, N.Y. (48). The authors of those studies reported that Ostreococcus was a significant component of PPE assemblages. Despite occasionally high abundances (Ͼ1 ϫ 10 5 cells ml Ϫ1 ), PPE from the Mediterranean site averaged 3.5 ϫ 10 4 cells ml Ϫ1 over a 27month time series at a depth of 0.5 m (67), a value similar to the interannual average DCM abundance of Ostreococcus at our study site (3.1 ϫ 10 4 cells ml Ϫ1 ) but greater than the average at 5 m (8.5 ϫ 10 3 cells ml Ϫ1 ) (Fig. 7). Ostreococcus at our study site was frequently present at abundances of Ͼ1.0 ϫ 10 4 cells ml Ϫ1 (14 of 23 samples from the DCM, but only 5 of 23 samples from 5 m) (Fig. 7).
Discrete cell-specific carbon conversion factor of 212 fg C cell Ϫ1 (70). The use of this carbon conversion factor is justified because the culture from which this value was derived (CCE9901) was isolated from a location (Scripps Pier) relatively close (ϳ130 km) to our sampling site. Furthermore, CCE9901 is a clade A Ostreococcus ecotype (32), which is the only Ostreococcus ecotype that has been detected at our site, sharing Ն99% 18S rRNA gene sequence identity with all Ostreococcus clones detected at the USC Microbial Observatory (Countway and Caron, unpublished data). The biomass of Ostreococcus reached 17.5 g C liter Ϫ1 at 5 m and 67.7 g C liter Ϫ1 at the DCM. The upper range of carbon biomass estimated for total phototrophic picoeukaryotes at Scripps Pier, a location where PPEs dominated the biomass of the picoplankton assemblage, was 33.3 g C liter Ϫ1 (70). Thus, Ostreococcus appears to be an important component of the picoeukaryote community at our coastal Pacific location, particularly at the depth of the DCM. The DCM is a common feature of the water column off the coast of Southern California and has been shown to occur near the depth of the nitracline and the 10% light level (19). It is probable that Ostreococcus plays an important role in microbial food webs within the Southern California Bight region, given its prevalence on the northern (this study) and southern (70) sides of the Bight.
Ostreococcus populations displayed tremendous resiliency at our station, showing numerous but variable peaks in abundance over 2 years (Fig. 7). Reports of growth rates for control and nutrient-amended field samples dominated by Ostreococcus range from 2 to 8 day Ϫ1 , with grazing losses in the range of 1.5 to 6.5 day Ϫ1 (27). High growth rates and grazing losses could explain the large month-to-month variability that we observed for Ostreococcus abundance. Rapid fluctuations between high and low abundances are a normal feature of picoeukaryote populations (46,67,70) and suggest that when conditions become favorable for growth (e.g., sudden nutrient availability or relief from grazing pressure), phototrophic picoeukaryote populations increase rapidly.
Two spring upwelling events were observed during our study, with one occurring between March and June 2001 and the other occurring between April and July 2002 (Fig. 6B). Nutrient concentrations in the upper 40 m of the water column increased rapidly during the early spring of 2001, while increases were more gradual and temporally shifted in 2002. Nutrient input corresponded to a threefold increase in diatom abundance at 5 m and a fivefold increase at the DCM in April 2001, while the Ostreococcus abundance did not increase until the following month (Fig. 7). Diatom abundance at the DCM returned to prebloom levels during the interval between sample collections, while Ostreococcus abundance increased Ͼ2 orders of magnitude, to 3.2 ϫ 10 5 cells ml Ϫ1 , in May 2001. A substantial increase in Ostreococcus abundance was also observed at 5 m between April and May 2001, but the diatom abundance remained unchanged. The nutrient upwelling in 2002 was of similar magnitude to the upwelling in 2001, but the highest nutrient concentrations occurred later in the year, shifting the spring diatom bloom into the summer. Ostreococcus cells did not bloom to high abundance at either depth following the upwelling in 2002. Factors controlling the abundance of Ostreococcus remain to be investigated at our study site.
The pattern of bloom succession from diatoms to other phytoplankton has been observed at many locations (4,38,58,61). Smayda and Villareal (62) suggested that an "open phytoplankton niche" develops after a typical spring diatom bloom, which would allow nondiatom phytoplankton to become relatively more abundant. This concept was extended to the "picoalgal niche" hypothesis, where particular picoplankton become dominant depending on the nutrient regime of the environment and the relative differences in grazing pressure experienced by various taxa that could occupy this niche (59,60). Picoeukaryotic phytoplankton are thought to be abundant in Thau Lagoon because of the relatively higher grazing pressures exerted upon larger-sized classes of phytoplankton and the removal of picoeukaryote grazers by cultured oysters (22,67). In general, picoeukaryotic phytoplankton in pelagic systems are thought to be more susceptible to grazing losses than larger phytoplankton because of the higher growth rates of the grazers (ciliates and flagellates) that feed upon the smallest phytoplankton (51). Many of these hypotheses remain to be adequately tested for populations of Ostreococcus; however, our newly developed method and preliminary observations may help to advance this endeavor. The sensitivity of our qPCR method will help to answer questions about the Ostreococcus dynamics at low cell abundance (prior to bloom formation), and its speed will permit the analysis of large numbers of samples.
In summary, we have developed and applied a real-time TaqMan-based qPCR assay for the detection and ecological study of Ostreococcus sp. This dynamic prasinophyte was detected at the USC Microbial Observatory on 22 of 23 sampling dates spanning a 24-month time series. Ostreococcus formed a large bloom at our study site in May 2001, reaching abundances of 8.2 ϫ 10 4 cells ml Ϫ1 at 5 m and 3.2 ϫ 10 5 cells ml Ϫ1 at the DCM. Our qPCR approach readily detected low abundances of Ostreococcus, revealing a vertical distribution profile that was relatively typical for other phytoplankton species, with an abundance maximum near the DCM. The detection of Ostreococcus at such low concentrations shows promise for investigating "prebloom" dynamics of this picoeukaryote in the future.