Geochemical Rate-RNA Integration Study: Ribulose-1,5-Bisphosphate Carboxylase/Oxygenase Gene Transcription and Photosynthetic Capacity of Planktonic Photoautotrophs

doi:10.1128/AEM.70.9.5459-5468.2004

Appl Environ Microbiol. 2004 September; 70(9): 5459–5468.

doi: 10.1128/AEM.70.9.5459-5468.2004

PMCID: PMC520902

Geochemical Rate-RNA Integration Study: Ribulose-1,5-Bisphosphate Carboxylase/Oxygenase Gene Transcription and Photosynthetic Capacity of Planktonic Photoautotrophs

Jorge E. Corredor,¹ Boris Wawrik,² John H. Paul,² Hiep Tran,³ Lee Kerkhof,³^* José M. López,¹ Angel Dieppa,¹ and Oswaldo Cárdenas¹

Department of Marine Sciences, University of Puerto Rico, Lajas, Puerto Rico,¹ College of Marine Science, University of South Florida, St. Petersburg, Florida,² Institute of Marine and Coastal Sciences, Rutgers University, New Brunswick, New Jersey³

^*Corresponding author. Mailing address: Institute of Marine and Coastal Sciences, Rutgers University, New Brunswick, NJ 08901. Phone: (732) 932-6555, ext. 335. Fax: (732) 932-6520. E-mail: kerkhof/at/imcs.rutgers.edu.

Received October 15, 2003; Accepted May 5, 2004.

This article has been cited by other articles in PMC.

Abstract

A pilot field experiment to assess the relationship between traditional biogeochemical rate measurements and transcriptional activity of microbial populations was carried out at the LEO 15 site off Tuckerton, N.J. Here, we report the relationship between photosynthetic capacity of autotrophic plankton and transcriptional activity of the large subunit gene (rbcL) for ribulose-1,5-bisphosphate carboxylase/oxygenase (RubisCO), the enzyme responsible for primary carbon fixation during photosynthesis. Similar diel patterns of carbon fixation and rbcL gene expression were observed in three of four time series, with maxima for photosynthetic capacity (P_max) and rbcL mRNA occurring between 10 a.m. and 1 p.m.. The lowest P_max and rbcL levels were detected between 6 p.m. and 10:30 p.m.. A significant correlation was found between P_max and form ID rbcL mRNA (R² = 0.56) and forms IA and IB (R² = 0.41 and 0.47, respectively). The correlation between the abundance of “diatom” rbcL and P_max mRNA was modest (R² = 0.49; n = 12) but improved dramatically (R² = 0.97; n = 10) upon removal of two outliers which represented afternoon samples with high P_max but lower mRNA levels. Clone libraries from reverse transcription-PCR-amplified rbcL mRNA indicated the presence of several chromophytic algae (diatoms, prymnesiophytes, and chrysophytes) and some eukaryotic green flagellates. Analogous results were obtained from amplified small rRNA sequences and secondary pigment analysis. These results suggest that diatoms were a major contributor to carbon fixation at LEO 15 at the time of sampling and that photosynthetic carbon fixation was partially controlled by transcriptional regulation of the RubisCO gene.

A primary research goal for marine environmental biogeochemistry is a thorough understanding of the molecular to global scale links and feedback mechanisms between solar irradiance, marine microbial activity, primary productivity, carbon and nitrogen cycles, and remotely sensed ocean color data. Increased emphasis is now being placed on the application of modern molecular tools to address these problems. Understanding of the genetic basis of biogeochemical processes and, more importantly, the expression of relevant genes in the marine environment is perceived as a promising route to these goals. A detailed understanding of the connection between gene expression, protein translation, and enzyme activity in environmental systems can provide significant insight into the mechanisms controlling biogeochemical transformation rates.

To this end, a pilot field experiment was carried out to test whether mRNA levels from complex, indigenous phytoplankton and bacterial populations can be related to biogeochemical rate processes mediated by these communities. The Geochemical Rate-RNA Integration Study (GRIST) incorporated several different biogeochemical measurements and molecular approaches (6). The experiment was designed to provide guiding principles on how the combination of the flux determinations and the quantification of microbial transcriptional activity can be used to better understand the regulation and response of microbial activity to different environmental parameters. Measurements have focused on C and N cycling and included primary production, bacterial production, dissolved organic matter production and uptake, inorganic and organic N assimilation, and N redox cycling in water column and sediment samples.

Concurrently, a suite of molecular approaches were used to quantify gene product abundance over two diel cycles. We employed partial-length gene probes to identify the major clades of rbcL mRNA in oceanic phytoplankton (13, 18, 23). Form IA rbcL is found in phycoerythrin-containing picocyanobacteria such as Synechococcus as well as Prochlorococcus. Form IB rbcL is found in chlorophytes and certain picocyanobacteria, whereas form ID of the enzyme is found in chromophytic phytoplankton (diatoms, prymnesiophytes, pelagophytes, bolidophytes, rhodophytes, chrysophytes, and some dinoflagellates) (22, 23). We here report on the relationship between the photosynthetic capacity of autotrophic plankton and the transcriptional activity of variants of the large subunit gene (rbcL) for ribulose-1,5-bisphosphate carboxylase/oxygenase (RubisCO), the enzyme in the first step of the Calvin-Benson-Bassham pathway of autotrophic carbon fixation (21).

MATERIALS AND METHODS

The GRIST experiment was carried out at the Rutgers University Marine Field Station at Tuckerton, N.J. (19 to 25 July 2002). Sampling was carried out at the well-characterized Long-Term Ecosystem Observatory (4). Diel periodicity was explored at 4-h intervals over two daily cycles (20 July 2002 and 22 July 2002), referred to hereafter as day 1 and day 2, respectively. Filtration for biogeochemical and molecular analyses occurred at sea aboard the R/V Arabella, and samples were flash-frozen on deck, either on dry ice or in liquid N₂.

Photosynthetic pigments and elemental analysis.

Chlorophyll a concentration of GF/F-filtered samples was determined fluorometrically after 90% acetone extraction with the method of Welschmeyer (25). High-pressure liquid chromatographic analysis of photosynthetic pigments was carried out following the procedure described by Wright et al. (26) with a C₁₈ Waters Symmetry column (3.9 by 150 mm) on a Shimadzu LC-10AT instrument equipped with an SPD-M10AV diode array detector. The high-pressure liquid chromatograph was calibrated with a pure (Sigma Chemical) chlorophyll a solution in 90% acetone. To increase resolution and allow separation of zeaxanthin from lutein, the mobile-phase program was modified as follows: 0 to 4 min, 1 ml min⁻¹, B (90:10 methanol-0.5 M aqueous ammonium acetate, pH 7.2, vol/vol), ramp to C (90:10 acetonitrile-water, vol/vol); 4 to 5.5 min, increase flow to 1.8 ml min⁻¹, 4 to 18 min, ramp to 20% C, 80% D (ethyl acetate); 18 to 21 min, 1 ml min⁻¹ ramp to C; 21 to 24 min, ramp C to B for restart. Particulate organic C and N content was determined with a CE Elantech CHN elemental analyzer.

Photosynthetic parameters.

Water samples for experimental determination of phytoplankton photosynthetic capacity representative of surface and bottom (14 m) assemblages were transported to the dockside laboratory at the Rutgers University Marine Field Station (<45-min transit time). Samples were transferred to 1-liter light-shielded, acid-washed polyethylene bottles, of which a 650-ml subsample was spiked with 0.108 mCi of [¹⁴C]bicarbonate (Amersham Pharmacia Biotech) for a 0.167-μCi initial activity. Aliquots (40 ml) of the spiked water were transferred to 40-ml borosilicate Environmental Protection Agency vials and incubated for 2 h in a photosynthetron apparatus (CHPT Mfg Co.) at constant temperature (17°C) and irradiances ranging from 0 (dark sample) to 614 μE m⁻² s⁻¹. A time-zero sample blank was immediately filtered prior to commencement of incubation. Following incubation, samples were filtered onto 25-cm GF/F glass fiber filters and treated with 250 μl of 10% HCl to drive off unfixed [¹⁴C]bicarbonate. After 24 h, 10 ml of BCS scintillation fluid (Amersham Pharmacia Biotech) were added, and sample radioactivity was determined on a Beckman L6000 liquid scintillation counter with the channels ratio mode. The resulting data were plotted in P versus E (irradiance) curves, and the photosynthetic parameters α (response to low light), P^B_max (photosynthetic capacity normalized to chlorophyll a content), and β (susceptibility to photoinhibition) were computed with the exponential formulation of Platt et al. (20).

rbcL mRNA analysis by probe hybridization.

mRNA was extracted from seawater with RNeasy spin columns (QIAGEN) as previously described (12). Briefly, either 100- or 200-ml seawater samples were treated with 0.1% diethylpyrocarbonate (Sigma Chemical Corp.) and filtered onto 25-mm, 0.45-μm HV polyvinylidene difluoride filters (Millipore Durapore). Filters were stored in liquid nitrogen in 750 μl of RLT lysis buffer (Qiagen) together with 0.2 g of baked glass beads (Biospec Products). Cell lysis was achieved by bead beating; 500 μl of lysate was then extracted following the RNeasy kit (Qiagen) protocol. Samples were split three ways. One third was digested with DNase-free RNase (Promega, Madison, Wis.), and one third was digested with RQ1 DNase (Promega). RNA was then immobilized onto Zeta-Probe charged nylon filters (Bio-Rad) by dot blotting and UV cross-linking. Duplicate samples were probed with form IA, form IB, and form ID rbcL probes previously described (7, 16). [³⁵S]UTP-labeled probes were prepared by in vitro transcription. Dot blots were analyzed with a Bio-Rad model GS363 molecular imager. Standard curves were made from opposite orientation in vitro transcripts generated from the same rbcL clones used to make the riboprobes.

To detect the proportion of form 1D hybridization signal by diatoms, a quantitative PCR assay was performed. Samples were filtered as above but the diethylpyrocarbonate treatment was omitted. Aliquots of the RNeasy eluates were amplified as described with chromophytic rbcL primers and a diatom-specific TaqMan probe (24).

Extraction of mRNA and rRNA for RT-PCR.

Two methods were used to obtain rbcL mRNA and intact ribosomes. For mRNA analyses, 200 ml of seawater was filtered from the 2:00 p.m. sample of the surface water from the first diel and filters were stored in 750 μl of RLT (Qiagen) and stored in liquid nitrogen. Samples were not treated with diethylpyrocarbonate, if they were intended for reverse transcription (RT)-PCRs. Upon extraction samples were DNase digested on the RNEasy columns for 15 min with the RNase-free DNase reagent set (Qiagen) according to the protocol provided by Qiagen.

Reverse transcriptions were performed with random hexamers, 4.7 mM MgCl, and Moloney murine leukemia virus reverse transcriptase (Promega) for 30 min at 37°C and freshly extracted RNA; 5 μl was then added to a PCR. Two primer sets were used to amplify the rbcL fragments. The form IA/B primer set produced a 615-bp fragment (forward primer: 5′-TCIGCITGRAACTAYGGTCG-3′; reverse primer: 5′-CTGAGIGGIAARAACTACGG-3′) and the form ID set produced a 554-bp fragment (forward primer: 5′-GATGATGARAAYATTAACTC-3′; reverse primer 5′-ATTTGDCCACAGTGDATACCA-3′). All sequences are stated with International Union of Pure and Applied Chemistry (IUPAC) degeneracy. PCR conditions were as follows: 1 μM final concentration for both primers, 1.5 mM MgCl₂, 2.5 mM each nucleotide, and 5 U of Taq polymerase (Promega). Cycle parameters were 3 min at 95°C, followed by 40 cycles of 1 min at 95°C, 1 min at 52°C, and 1 min 30 sec at 72°C. Cycling was followed by a 15-min 72°C elongation step.

For intact ribosome extraction from phytoplankton, a modified phenol-chloroform procedure incorporating multiple freeze-thaw cycles was used as described previously (7). Total nucleic acid pellets were resuspended in diethylpyrocarbonate-treated water, and the ribosomal small subunit was gel purified with an RNaid kit (Bio 101, Carlsbad, Calif.). The purified small subunit was diluted 10⁻⁵ to 10⁻⁶ and used as a template in an RT-PCR with an 18S rRNA eukaryotic forward primer (5′-ACC TGG TTG ATC CTG CCA G-3′ ) and a universal reverse primer (519R; 5′-ATTACC GCG GCT GCT GG-3′ ) with the Titan RT-PCR kit (Roche) containing the RNase inhibitor anti-RNase (Ambion). The cycling parameters were 50°C for 30 min, 94°C for 2.5 min, 25 cycles of 94°C for 1 min, 55°C for 2 min, and 72°C for 3.5 min, followed by a final extension at 72°C for 10 min.

Clone libraries.

rbcL genes from mRNA were RT-PCR amplified from one sampling time with the primer sets described above. Immediately after amplification, products were purified with the Quiaquick PCR purification kit (Qiagen) or with a Geneclean kit (Bio 101). Amplicons were ligated into the PCR II vector with a TA dual promoter cloning kit (Invitrogen). RubisCO clones were screened for the correct insert size by PCR and for unique restriction patterns with Sau3AI and AluI before sequencing in an Applied Biosystems model 373 sequencer by the University of Florida core sequencing laboratory (University of Florida, Gainesville, Fla.). rRNA clones were obtained from both surface and bottom samples during day 2. Clones were screened by terminal restrictions fragment length polymorphism analysis by digestion with MspI, and selected clones were sequenced on an ABI 310 as described previously (8).

Phylogenetic analysis.

RubisCO cDNA sequences were translated into amino acids and aligned with their closest match found in GenBank as determined by BLAST searching as well as sequences we previously recovered from similar environments (14, 23). Alignment was performed with Omiga 1.1 (Oxford Molecular Group, Oxford, United Kingdom) with a ClustalW pairwise weighted alignment method. Amino acid sequences were used to avoid potential biases introduced by codon usage and GC content. Alignments were inspected for obvious misalignments and exported to Mega 2.0 beta (10). Trees were built with both parsimony and the neighbor-joining method with a gamma distribution (gamma parameter = 2.0) to correct for rate heterogeneity across sites. rRNA genes were initially aligned with ClustalW and later aligned by hand. Phylogenetic trees were reconstructed by maximum-likelihood methods with fastDNAml.

Nucleotide sequence accession number.

Accession numbers for the RubisCO clones are AY356325 to AY356346. Accession numbers for the 18S rRNA clones are AY378156 to AY378164.

RESULTS

Water mass properties varied moderately during the experiment, with initial surface temperatures of 19°C increasing to 22°C during day 1 (Fig. (Fig.1).1). Day 2 surface waters increased in temperature from 22 to 24°C by the end of the experiment. Deeper water also experienced slight warming from 16 to 17°C, with upper layer stratification being apparent. Practical salinity during the sample period ranged from 31.6 (surface) to 32 (deep), with an increase of less than 0.2 in surface water at the 19:00 and 22:00 time points of day 2 (data not shown). An additional intrusion of warmer water at the sample site was also evident between 6 and 10 a.m. during day 1 at approximately 12:00. However, no samples for primary production, rbcL mRNA, or 18S rRNA were collected at these intermediate depths during the experiments.

FIG. 1.

Temperature versus depth data for day 1 (A) and day 2 (B).

Surface chlorophyll a, particulate organic C, and particulate organic N concentrations all increased from approximately 18:00 on day 1 to 19:00 on day 2, suggesting that a small phytoplankton bloom was under way (Fig. (Fig.2).2). The chlorophyll a concentration ranged between 3 and 16 mg m⁻³ during day 1 and between 2.6 and 34.7 mg m⁻³ during day 2 with deep water (14 m) chlorophyll a variation substantially lower than surface variation. Particulate organic C concentrations varied between 0.48 and 3.2 mg liter⁻¹ and particulate organic N concentrations ranged between 0.08 and 0.49 mg liter⁻¹ (Fig. (Fig.2).2). Both particulate organic C (r² = 0.86) and particulate organic N (r² = 0.90) were significantly correlated to chlorophyll a and for the most part followed a temporal pattern similar to that observed for chlorophyll a. Chlorophyll a, particulate organic C, and particulate organic N content were consistently greater near the surface than at depth during the course of the experiment.

FIG. 2.

Time course for chlorophyll a (Chl a), particulate organic C (POC), and particulate organic N (PON) at the LEO 15 node A site.

Relative to the Redfield C:N ratio of 7.1, surface particulate organic matter (8.8) was significantly depleted in N, while particulate organic matter at depth (5.9) was N rich. These differences may reflect enhanced pigment and protein synthesis in light-limited deep waters relative to greater carbohydrate and lipid synthesis in surface waters. Chlorophyll a values, as determined by HPLC, were about 65% of those reported by bulk fluorometry, indicating positive interference by other chlorophylls in the fluorometric method. Chlorophyll a was the dominant photosynthetic pigment. Fucoxanthin was the dominant accessory pigment representing 28% of the total chlorophyll a, followed by chlorophyll b (15% chlorophyll a). Fucoxanthin was significantly correlated with chlorophyll a (r² = 0.52, P < 0.01). Other pigments present in minor concentrations (<10% chlorophyll a) were, in order of abundance, chlorophyll c₂, zeaxanthin, peridinin, and chlorophyll c₁ (Table (Table11).

TABLE 1.

Concentration of photosynthetic pigments at LEO 15 during the GRIST experiment as determined by HPLC

Surface photosynthetic rates during the first day were only measured during the initial hours of the experiment when chlorophyll a variability was low. Light-saturated photosynthetic rates (P_max), the maximum achievable C fixation rates, varied little in these samples, as was the case for deep-water samples despite a greater range in phytoplankton abundance, as reflected by chlorophyll a, particulate organic C, and particulate organic N content. During the second diel experiment, P_max varied substantially (Fig. (Fig.3,3, Table Table22).

FIG. 3.

Photosynthesis-irradiance curves for surface and bottom phytoplankton populations.

TABLE 2.

Photosynthetic parameters of phytoplankton assemblages during the GRIST experiment

Photosynthetic capacity (P^B_max), the maximum achievable photosynthetic rate normalized to chlorophyll a concentration, ranged between 1.1 and 11.6 mg of C (mg of chlorophyll a)⁻¹ h⁻¹, well within the range of values reported previously for phytoplankton assemblages of temperate coastal marine waters (2). Average P^B_max was moderately higher at depth (6.9 m) than near the surface (4.8 m). The response to irradiance at low light levels (α) ranged between 0.03 and 0.25 mg of C (mg of chlorophyll a)⁻¹ μE⁻¹ m⁻² and was higher on average at depth (0.13) than near the surface (0.09). Photoinhibition was apparent only in deep-water samples.

In the surface waters, rbcL mRNA levels were highest either at the first-light sampling (dawn; day 1 surface waters) or at 10:30 a.m. (Fig. (Fig.4).4). In the subsurface, maximum transcription occurred at the 10:30 or 1 p.m. sampling. The minimum transcriptional activity occurred at about 7:00 p.m. in both surface and bottom samples. During the first diel in the surface samples, form ID rbcL mRNA was the most abundant form of rbcL mRNA during the morning peak in transcriptional activity (Fig. (Fig.4).4). During the other three diel samples (day 1 bottom water and both surface and bottom for day 2), diatom rbcL mRNA displayed the greatest dynamic range among the rbcL forms that we quantified.

FIG. 4.

Time course for rbcL mRNA abundance and light-saturated photosynthetic rate (P_max). IA, form IA rbcL mRNA; IB, form IB rbcL mRNA; ID, form ID rbcL mRNA; DRT, real-time diatom rbcL mRNA. Open squares, P^B_max. Right panels, day 1; left panels, day 2.

Taking the average rbcL mRNA levels for the daylight hours for this group of samples indicated that diatom rbcL was significantly greater than form IA or form IB, but not different from form ID (P = 0.046). Figure Figure55 shows the daylight average rbcL levels for all four time series. Bottom samples both showed a distinct increase in transcript abundance at the last sampling, 10:30 p.m. to 11:00 p.m., which was also reflected in particulate organic C and P^B_max. All three forms of rbcL mRNA quantified by hybridization were significantly correlated to each other at the 1% confidence interval (A to B: P = 0.001, A to D: P = 0.01, B to D: P = 0.001). Form IB and form ID rbcL mRNAs were also significantly correlated with diatom rbcL mRNA quantified by real-time PCR. Form IA was not (P = 0.12). P_max was significantly correlated with all forms of rbcL that we quantified (A: P = 0.01, B: P = 0.005, D: P = 0.001, diatom: P = 0.007; see Fig. Fig.9).9). P_max was least well correlated to form IA rbcL (R² = 0.41). P_max exhibited progressively tighter correlation to form IB (R² = 0.47) and to form ID rbcL mRNA (R² = 0.56). The best predictor of P_max was the cumulative rbcL mRNA signal obtained by adding form IA, IB, and ID hybridization numbers (R² = 0.58).

FIG. 5.

Relative abundance of the different forms of rbcL mRNA observed during the experiment (averages of the daylight hours).

FIG. 9.

Relationship between PCR rbcL mRNA abundance and light-saturated photosynthetic rate. The linear relationship is given by y = −6 + 95.5x (r² = 0.58; n = 16).

Cloning and sequencing of rbcL mRNA from transcriptionally active phytoplankton yielded 22 unique rbcL clones (Fig. (Fig.6).6). Table Table33 gives a list of all unique sequences detected and their closest relatives in GenBank as determined by BlastP analysis (1). Phylogenetic analysis suggested that 14 of these clones were of the form ID type and likely included four diatoms, four prymnesiophytes, four chrysophytes, a rhodophyte, and a deeply rooted sequence most closely related to Olithodiscus luteus. The remaining eight clones were of the form IB type and included two flagellated chlorophytes related to Pyramimonas, one sequence closely related to Bathycoccus, and a small clade of deeply rooted forms. No picocyanobacterium-like sequences (i.e., Synechococcus or Prochlorococcus) were detected.

FIG. 6.

Neighbor-joining tree of transcriptionally active rbcL clones obtained from the site (labeled GRIST and clone-specific suffix) and related sequences from GenBank. The clone labeled P994DY5 is from the Gulf of Mexico (25).

TABLE 3.

Closest relatives of rbcL clones in GenBank

To test whether diatoms, prymnesiophytes, or chrysophytes were actively growing in the water samples during the experiment, we created clonal libraries from intact rRNA purified from the samples with reverse transcriptase. The relationship between RNA content and growth rate has been well characterized for phytoplankton and bacteria (3, 9). A reconstructed phylogenetic tree for the 18S rRNA subunit is shown in Fig. Fig.7.7. Comparison of gene phylogenies obtained for cloned rbcL and 18S sequences revealed several important similarities and discrepancies. The GRIST CH25 rbcL clone was found to be most similar (>92% at the amino acid level, Table Table3)3) to Detonula confervacea, as demonstrated by phylogenetic analysis (Fig. (Fig.5).5). Similarly, the 18S clones G2D2T4 and C12D2T4 (Fig. (Fig.6)6) were found to be most closely related to Detonula. Two rbcL sequences were observed to be most similar to the diatom Skeletonema costatum (GRIST CH19 and CH27), but no corresponding 18S clones were observed. The GRIST CH29 clone was related to Pedinella, matching the B12D2T4 and G1D2T3 18S clones from both surface and deep samples. Sequences similar to a red alga, green algae, and one almost identical to Bathycoccus prasinos were not observed among the sequence of the corresponding 18S library. Finally, the GRIST CH36 rbcL related to E. huxleyi was similar to the D1D2T3 18S clone. These results support the observation that diatoms and prymnesiophytes were actively photosynthesizing and growing during the experiment.

FIG. 7.

Phylogenetic tree based on 18S rRNA genes obtained via reverse transcription of intact ribosomes during GRIST with maximum-likelihood methods.

DISCUSSION

Finding molecular proxies that can distinguish and quantify biogeochemically active members in a microbial population should provide a better understanding of how such populations respond to environmental parameters. Insight into the regulation and expression of photosynthetic processes is at the core of the global C cycle and is particularly pertinent to current environmental concerns. In this study, we sought to determine what photosynthetically active phytoplankton groups were present, the photosynthetic capacity of the assemblage, and the relationship of this activity to different sources of rbcL mRNA.

Diel variation of marine phytoplankton photosynthetic parameters has been well documented (11), and recent results indicate that such variability is in large part due to modulation of RubisCO transcription (19). The overriding factor in the regulation of RubisCO in vivo is light (3, 5, 14). Transcripts are known to accumulate in the light and are degraded in the dark in a number of prokaryotic and eukaryotic algae. In the environment, rbcL transcription has been investigated in several studies (13, 16, 17, 19, 23, 24, 27). Mesocosm experiments and a Langrangian study in the oligotrophic Gulf of Mexico revealed strong diel variation in form IB rbcL mRNA levels in surface waters, suggesting diel patterns in the regulation of rbcL transcription (19).

Form IB rbcL mRNA was significantly correlated with carbon fixation and was elevated during morning hours. In the water column, the highest cyano-rbcL transcripts (Synechococcus and Prochlorococcus) were observed at shallower depths, while chromo-rbcL transcripts were highest at the deep chlorophyll maximum, where chromophytic picoeukaryotes dominated, as determined by flow cytometry or epifluorescence microscopy (18). In addition to spatial segregation, a temporal separation in the peaks of chromo (form ID) and cyano (form IB) transcripts was also observed. Cyanobacterial rbcL transcript was predominantly found during morning hours, while chromophitic rbcL mRNA was most abundant in the late afternoon (16). However, both cyano rbcL mRNA and carbon fixation peaked in the early morning hours and were significantly correlated, suggesting that form IB rbcL mRNA was a good proxy for carbon fixation in these waters.

In contrast to the Gulf of Mexico and Tampa Bay studies, the Mid-Atlantic Bight (LEO 15) analysis of the transcriptionally active rbcL clones from a single time point indicated an abundance of chromophytic algae (a group that includes diatoms, prymnesiophytes, and chrysophytes), with a lesser amount of chlorophytic flagellates. These observations are in agreement with the accessory pigment analysis by HPLC and the major forms of rbcL mRNA encountered by probe hybridization (i.e., form ID and diatom-specific rbcL mRNA; Fig. Fig.44 and and5).5). This was also supported by the analysis of 18S rRNA clone libraries, which yielded sequences of similar phylogenetic affinity (in certain instances) as our rbcL clone libraries. It was also interesting to note the absence of oligotrophic oceanic picoplankton such as Synechococcus and Prochlorococcus in our rbcL clone libraries.

HPLC analysis of GRIST samples revealed predominance of pigments associated with chromophytes (fucoxanthin) and chlorophytes (chlorophyll b). Minor occurrences of peridinin and the absence of 19′ butanoyloxyfucoxanthin and 19′ hexanoyloxyfucoxanthin indicated predominance of diatoms relative to dinoflagellates and absence of prymnesiophytes and chrysophytes. The high proportions of fucoxanthin and chlorophyll b are indicative of a phytoplankton assemblage dominated by chromophytes and chlorophytes, while the low concentrations of zeaxanthin indicated lesser numbers of cyanophytes (Fig. (Fig.8).8). Thus, diatoms and chlorophytes appear to have been the predominant carbon-fixing components of the phytoplankton at the time of our experiment based on all assays: probe hybridization (rbcL mRNA), clonal analysis by RT-PCR (rbcL and 18S rRNA genes), and biochemical (HPLC) analyses.

FIG. 8.

Relative abundance of photosynthetic pigments as a percentage of chlorophyll a content.

The relationship between rbcL mRNA and P_max during the GRIST experiment may have been affected by environmental variability as advection of water masses with different physical properties (Fig. (Fig.1)1) may have resulted in the replacement of phytoplankton populations. Large changes in pigment composition associated with an increase in thermocline intensity suggest a population shift in surface waters during the morning of day 1 (Fig. (Fig.8).8). However, afternoon increases in diatom dominance as reflected in fucoxanthin concentrations did not correspond to peaks in diatom or form ID mRNA in three of four diel samples, suggesting that the transcriptional signal overshadowed the changes in community composition or in situ pigment synthesis. Total mRNA yields (total mRNA was quantified in all samples prior to dotting; data not shown) exhibited good correlations with pigment levels, particularly to chlorophyll c₁, fucoxanthin and zeaxanthin (R² = 0.68, 0.57, and 0.77, respectively; P < 0.001 for all three), but were not significantly correlated with the individual mRNA signals (form IA: P = 0.38; IB: P = 0.17; ID: P = 0.29; diatom: P = 0.627). This indicates that variations in rbcL mRNA levels were not caused by population level-dependent changes of total mRNA, but were caused by transcriptional regulation of phytoplankton rbcL mRNA.

Correlation between rbcL mRNA and P_max is, as expected, dependent on the particular type of rbcL mRNA examined (Fig. (Fig.9),9), the most robust correlation (R² = 0.56) corresponding to type ID rbcL mRNA (R² = 0.56) and the least robust (R² = 0 0.41) corresponding to type IA. Overall, correlation between rbcL hybridization data and photosynthetic carbon fixation was statistically significant at R² = 0.58 (P = 0.00064). Real-time diatom rbcL mRNA correlated moderately with P_max (R² = 0.49; n = 12). However, two afternoon (3:40 p.m. on day 1 and 4:20 p.m. on day 2) surface samples showed barely detectable levels of diatom rbcL mRNA despite the high photosynthetic capacity observed at these times. We attribute this to mRNA degradation in the presence of high levels of functional RubisCO protein. Removal of the two data points yields a remarkably tight correlation (R² = 0.97 n = 10) indicating close coupling between gene transcription and photosynthetic rate (Fig. (Fig.10)10) for most of the day.

FIG. 10.

Relationship between real-time PCR diatom rbcL mRNA abundance and light-saturated photosynthetic rate. The linear relationship, excluding the outliers marked in gray, is given by y = 13.2 + 57.1x (r² = 0.97; n = 10).

Our results suggest that transcriptional regulation may have been important for certain time series (day 2) of the GRIST experiment after initial bloom development. Correlations between rbcL mRNA and P_max were low during day 1, but biomass and photosynthetic rate variability were modest at this time. Closer correlations were observed on day 2 when light-saturated photosynthetic rates varied over one order of magnitude (13 to 160 mg of C liter⁻¹ h⁻¹). Late-night increases in both mRNA abundance and photosynthetic capacity of the deep-water assemblage constitute a novel observation.

Although transcription has been used to some degree as a surrogate for gene expression, uncertainty remains regarding the predictability of biogeochemical activity based upon transcript abundance alone. Enzymes are subject to a host of posttranscriptional controls and RubisCO control is complex, particularly in eukaryotes. Several mechanisms for metabolic regulation of RubisCO activity are known (5). These include reversible activation of the active site by CO₂, binding of noncompetitive inhibitors to the active site and binding of effectors at sites other than the active site. In general, these controls respond to the relative abundance of RubisCO substrates and to NADP/NADPH and ADP/ATP ratios. Posttranscriptional mechanisms of regulation are active in eukaryotes (13), but transcriptional regulation appears to be the overriding factor in the regulation of carbon fixation in cyanobacteria (15, 22).

In summary, a strong relationship between photosynthetic capacity and transcription of the rbcL gene by diatoms was observed in three of four diel time series at the LEO 15 site off Tuckerton, N.J. This implies that transcriptional regulation of rbcL genes can be a major factor controlling carbon fixation in this coastal environment.

Acknowledgments

GRIST was carried out under the auspices of the U.S. Department of Energy Biotechnological Investigations Ocean Margin Program (BIOMP; grants DE-FG02-00ER62978 to L.K. and J.E.C. and DE-FG02-97ER62452 to J.H.P.).

We thank Lora McGuinness for sample collection and logistics, Alvaro Cabrera for assistance in the photosynthetron incubations, Ana Lozada for analysis of chlorophyll a, particulate organic C, and particulate organic N, and the crew of the RVs Arabella and Caletta for sample collection and transport.

REFERENCES

1. Altschul, S. F., T. L. Madden, A. A. Schaffer, J. Zhang, Z. Zhang, W. Miller, and D. J. Lipman. 1997. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25:3389-3402. [PMC free article] [PubMed]

2. Côté, B., and T. Platt. 1983. Day-to-day variations in the spring-summer photosynthetic parameters of coastal marine phytoplankton. Limnol. Oceanogr. 28:320-344.

3. Dortch, Q., T. L. Roberts, J. J. R. Clayton, and S. I. Ahmed. 1983. RNA/DNA ratios and DNA concentrations as indicators of growth rate and biomass in planktonic marine organisms. Mar. Ecol. Prog. Ser. 13:61-71.

4. Glenn, S. M., M. F. Crowley, D. B. Haidvogel, and Y. T. Song. 1996. Underwater observatory captures coastal upwelling events of New Jersey. Eos Trans. Am. Geophys. Union 77:233-236.

5. Hartman, F. C., and M. R. Harpel. 1994. Structure, function, regulation, and assembly of d-ribulose-1,5-bisphosphate carboxylase/oxygenase. Annu. Rev. Biochem. 63:197-234. [PubMed]

6. Kerkhof, L., J. Corredor, J. Paul, D. Bronk, J. Lopez, and J. Cherrier. 2003. Experiment explores intercalibration of biogeochemical flux and nucleic acid measurements. Eos Trans. Am. Geophys. Union 84:167.

7. Kerkhof, L., and P. Kemp. 1999. Small ribosomal RNA content in marine proteobacteria during non-steady-state growth. FEMS Microbiol Ecol. 30:253-260. [PubMed]

8. Kerkhof, L., M. Santoro, and J. Garland. 2000. Response of soybean rhizosphere communities to human hygiene water addition as determined by community-level physiological profiling (CLPP) and terminal restriction fragment length polymorphism (TRFLP) analysis. FEMS Microbiol. Lett. 1994:95-101. [PubMed]

9. Kerkhof, L., and B. B. Ward. 1993. Comparison of nucleic acid hybridization and fluorometry for measurement of the relationship between RNA/DNA ratio and growth rate in a marine bacterium. Appl. Environ. Microbiol. 59:1303-1309. [PMC free article] [PubMed]

10. Kumar, S., K. Tamura, and M. Nei. 1993. MEGA: molecular evolutionary genetic analysis, version 2.0. The Pennsylvania State University, University Park.

11. MacCaull, W., and T. Platt. 1977. Diel variations in the photosynthetic parameters of coastal marine phytoplankton. Limnol. Oceanogr. 22:723-731.

12. Paul, J. H. 2001. Marine microbiology, vol. 30. Dowden, Hutchinson & Ross, New York, N.Y.

13. Paul, J. H., A. Alfreider, J. B. Kang, R. A. Stokes, D. Griffin, L. Campbell, and E. Ornolfsdottir. 2000. Form IA rbcLrbcL transcripts associated with a low salinity/high chlorophyll plume (′Green River') in the eastern Gulf of Mexico. Mar. Ecol. Prog. Ser. 198:1-8.

14. Paul, J. H., A. Alfreider, and B. Wawrik. 2000. Micro- and macrodiversity in rbcLrbcL sequences in ambient phytoplankton populations from southeastern Gulf of Mexico. Mar Ecol. Prog. Ser. 198:9-18.

15. Paul, J. H., J. B. Kang, and F. R. Tabita. 2000. Diel patterns of regulation of rbcLrbcL transcription in a cyanobacterium and a prymnesiophyte. Mar. Biotechnol. 2:429-436. [PubMed]

16. Paul, J. H., S. L. Pichard, J. B. Kang, G. M. F. Watson, and F. R. Tabita. 1999. Evidence for a clade-specific temporal and spatial separation in ribulose bisphosphate carboxylase gene expression in phytoplankton populations off Cape Hatteras and Bermuda. Limnol. Oceanogr. 44:12-23.

17. Pichard, S., M. E. Frisher, and J. H. Paul. 1993. Ribulose bisphosphate carboxylase gene expression in subtropical marine phytoplankton populations. Mar. Ecol. Prog. Ser. 101:55-65.

18. Pichard, S. L., L. Campbell, K. Carder, J. B. Kang, J. Patch, F. R. Tabita, and J. H. Paul. 1997. Analysis of ribulose bisphosphate carboxylase gene expression in natural phytoplankton populations by group-specific gene probing. Mar. Ecol. Prog. Ser. 149:239-253.

19. Pichard, S. L., L. Campbell, J. B. Kang, F. R. Tabita, and J. H. Paul. 1996. Regulation of ribulose bisphoshate carboxylase gene expression in natural phytoplankton communities. I. Diel rhythms. Mar. Ecol. Prog. Ser. 139:257-265.

20. Platt, T., S. Sathyendranath, and P. Ravindran. 1990. Primary production by phytoplankton: analytic solutions for daily rates per unit area of water surface. Proc. R. Soc. Lond. Ser. B 241:201-211.

21. Tabita, F. R. 1988. Molecular and cellular regulation of autotrophic carbon dioxide fixation in microorganisms. Microbiol. Rev. 52:155-189. [PMC free article] [PubMed]

22. Watson, G. M., and F. R. Tabita. 1996. Regulation, unique gene organization, and unusual primary structure of carbon fixation genes from a marine phycoerythrin-containing cyanobacterium. Plant Mol. Biol. 32:1103-1115. [PubMed]

23. Wawrik, B., J. H. Paul, L. Campbell, D. Griffin, L. Houchin, A. Fuentes-Ortega, and F. Mueller-Karger. 2003. Vertical structure of rbcL-rbcL-containing phytoplankton phylotypes associated with a coastal plume in the Gulf of Mexico. Mar. Ecol. Prog. Ser. 251:87-101.

24. Wawrik, B., J. H. Paul, and F. R. Tabita. 2002. Real-time PCR quantification of rbcLrbcL (ribulose-1,5-bisphosphate carboxylase/oxygenase) mRNA in diatoms and pelagophytes. Appl. Environ Microbiol. 68:3771-3779. [PMC free article] [PubMed]

25. Welschmeyer, N. A. 1994. Fluorometric analysis of chlorophyll aa in the presence of chlorophyll bb and pheogiments. Limnol. Oceanogr. 39:1985-1992.

26. Wright, S. W., and S. W. Jeffrey. 1997. High-resolution HPLC system for chlorophylls and carotenoids of marine phytoplankton, p. 343-360. InIn S. W. Jeffrey, R. F. C. Mantoura, and S. W. Wright (ed.), Phytoplankton pigments in oceanography. UNESCO, Paris, France.

27. Wyman, J., J. P. Zehr, and D. G. Capone. 1996. Temporal variability in nitrogenase gene expression in natural populations of the marine cyanobacterium Trichodesmium thiebautiiTrichodesmium thiebautii. Appl. Environ. Microbiol. 62:1073-1075. [PMC free article] [PubMed]

Articles from Applied and Environmental Microbiology are provided here courtesy of
American Society for Microbiology (ASM)