Previous Article | Next Article ![]()
Applied and Environmental Microbiology, June 2005, p. 2925-2933, Vol. 71, No. 6
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.6.2925-2933.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Laboratory of Applied Microbiology, Marine Biotechnology Institute, Heita, Kamaishi, Iwate 026-0001,1 Graduate School of Natural Science and Technology, Kanazawa University, Kakuma, Kanazawa, Ishikawa 920-1192, Japan2
Received 3 October 2004/ Accepted 22 December 2004
|
|
|---|
-Proteobacteria were the most abundant, and those affiliated with
-Proteobacteria were the second-most abundant. In contrast, the Kamaishi library was occupied by phylotypes affiliated with Planctomycetes,
-Proteobacteria,
-Proteobacteria, and Crenarchaeota. In the
-Proteobacteria, many Yamada phylotypes were related to free-living and symbiotic sulfur oxidizers, whereas the Kamaishi phylotype was related to the genus Pseudomonas. These results allowed us to hypothesize that sulfate-reducing and sulfur-oxidizing bacteria have become abundant in the Yamada sediment. This hypothesis was supported by quantitative competitive PCR (qcPCR) with group-specific primers. The qcPCR also suggested that organisms closely related to Desulfotalea in the Desulfobulbaceae were the major sulfate-reducing bacteria in these sediments. In addition, potential sulfate reduction and sulfur oxidation rates in the sediment samples were determined, indicating that the sulfur cycle has become active in the Yamada sediment beneath the areas of intensive shellfish aquaculture. |
|
|---|
-Proteobacteria,
-Proteobacteria, Bacteroidetes, and Planctomycetes are abundant and ecologically important. Information obtained in these studies has been combined with results of activity measurements (19, 20, 34), suggesting that sulfate reduction is the most important electron-accepting process for organic matter decomposition in coastal and continental shelf sediments. Even though much information is available that gives a general view of phylogenetic makeup and activity of prokaryotes in marine sediment, relatively little is known about the formation and transition of prokaryotic communities in response to changes in environmental conditions (27, 30, 33). Rosselló-Mora et al. have examined the response of the sediment microbial community to cyanobacterial biomass input under laboratory simulation conditions (33). A field study on sediments of eutrophic bays suggested that the sediment microbial community changes in response to eutrophication (30). Considering that prokaryotic populations rapidly respond to environmental perturbations and are reliably studied by application of molecular ecological tools, changes in sediment prokaryotes will be useful indices for assessing the state of coastal environments.
The Sanriku coast is situated in northeastern Japan (Fig. 1) and is famous for the complex Rias coastline. The Rias coastline is punctuated with numerous small, semienclosed bays, many of which have been used as fishing ports. In addition, some bays have also been utilized for intensive commercial shellfish aquaculture (mainly oysters and scallops) and seaweeds for over 30 years, which has become one of the important industries in the Sanriku region. Consequently, sustainable usage of the coastline has currently become an important issue for the Sanriku community, as is the case with other coastal communities.
![]() View larger version (41K): [in a new window] |
FIG. 1. Map of study sites.
|
|
|
|---|
Although the Yamada bay is larger in area (30 km2) than the Kamaishi bay (10 km2), they have similar enclosure indices (EI; 1.43 for the Yamada bay and 1.28 for the Kamaishi bay). An EI is estimated according to the equation
, where S is the bay area (in square kilometers), W is the width at the mouth of a bay (in kilometers), D1 is the depth at the center of a bay (in meters), and D2 is the depth at the mouth of a bay (in meters) (http://www.emecs.or.jp/closedsea-jp/closedsea-jp.htm). Since this index is correlated with the water exchange capacity, the index has been considered important for evaluating influences of biological activities in seawater to surrounding regions, including sediments. Water temperature in both bays ranged from 3°C (February and March) to 22°C (August and September). The degree of cleanliness of the seawater in these two bays is reported to be similar, as judged by chemical oxygen demand (COD) (<3 mg liter1), total nitrogen (<0.4 mg liter1) and total phosphorus (<0.04 mg liter1) (http://www.pref.iwate.jp/
hp0318/index.html). There has been no history of eutrophication in these bays.
Figure 1 presents sampling points in these bays and the aquaculture area in the Yamada bay. As shown in this figure, sampling points Y-A, Y-B, and Y-C in the Yamada bay were situated beneath shellfish beds, while Y-D was situated at the mouth of the Orikasa river and out of the aquaculture zone. The water depths of the sampling sites were between 5 and 20 m (Table 1). Sediments in these two bays are composed of fine sand and silt.
|
View this table: [in a new window] |
TABLE 1. Characteristics of marine sediment samples taken from the Yamada and Kamaishi baysa
|
The concentration of total sulfide (TS) in a sediment sample was measured using a Gastech detector (Gastech Corp., Kanagawa, Japan) according to the manufacturer's instructions. To determine the sulfate concentration in total sediment, a sediment sample was washed with 5 volumes of distilled water by being vigorously shaken. The washings were then pooled, and the sulfate concentration was determined by ion chromatography with an IA-100 ion analyzer (DKK Toa, Tokyo, Japan). For measurements of COD in sediment, approximately 30 mg of sediment was suspended in 2 ml of water and subjected to a COD analyzer (Hach, Loveland, CO). Total direct counts (TDCs) of microorganisms were determined by using fluorescence microscopy, after microorganisms were detached from sediment particles by gentle sonication (12) and stained with 4',6'-diamidino-2-phenylindole (DAPI). These analyses were conducted within 6 h after sampling.
DNA extraction and purification.
For extraction of DNA, sediment samples were frozen within 2 h after sampling and stored in a freezer before use. DNA was extracted from 5 g of sediment according to a method described previously (48). The extraction procedure included proteinase K treatment, hot detergent treatment, three cycles of freezing/thawing, and phenol-chloroform extraction (37). The extracted nucleic acids were recovered by ethanol precipitation (37) and subjected to purification by being passed through a hydroxyl apatite column (39). The nucleic acids were finally treated with RNase (37). The quality and quantity of DNA were checked by measuring its UV absorption spectrum (37).
PCR, cloning, and sequencing of 16S rRNA gene.
PCR amplification of 16S rRNA genes was performed using B341f or U515f as a forward primer and U1492r as a reverse primer (Table 2). Amplification was performed with a Progene thermal cycler (Techne, Cambridge, England) by using a 50-µl mixture containing 1.25 U of Taq DNA polymerase (Amplitaq Gold; Applied Biosystems, Foster City, CA), 10 mM Tris-HCl (pH 8.3), 50 mM KCl, 1.5 mM MgCl2, 0.001% (wt vol) gelatin, each deoxynucleoside triphosphate at a concentration of 200 µM, 50 pmol of each primer, and 10 ng of DNA. The PCR conditions used were as follows: 10 min of activation of the polymerase at 94°C; 30 cycles, each consisting of 1 min at 94°C, 1 min at 50°C, and 2 min at 72°C; and finally 10 min of extension at 72°C.
|
View this table: [in a new window] |
TABLE 2. PCR primers used in this study
|
Nucleotide sequence analysis.
Database searches for related 16S rRNA gene sequences were conducted using the GenBank database. The profile alignment technique of ClustalW, version 1.7 (41), was used to align the sequences. Each alignment was refined by visual inspection, and secondary structures were considered for the refinement analysis (18). A phylogenetic tree was constructed by the neighbor-joining method (36) and the maximum-likelihood method (13). Nucleotide positions at which any sequence had a gap or an ambiguous base were not included in the phylogenetic calculations. Checks for chimeric sequences were conducted by the chimera check program in the Ribosomal Database Project database (25), and possible chimeras were excluded from further analyses.
Quantitative competitive PCR (qcPCR).
Primers used for qcPCR are listed in Tables 2, and qcPCR assays are summarized in Table 3. Competitor fragments were produced using the competitive DNA construction kit (Takara, Ohtsu, Japan). The composition of PCR was as described above, except for the competitor fragment being added at a known copy number. The PCR conditions were as described above except that 35 cycles of amplification were performed and the annealing temperatures were as described in Table 3. Two microliters of the PCR product was analyzed by electrophoresis through 1.5% (wt vol) agarose gels with TBE buffer (37), and the gels were photographed after being stained with SYBR Gold (FMC Bioproducts, Vallensbaek Strand, Denmark). The band intensities of the target and competitor fragments were quantified by using the Multianalyst program supplied with Gel Doc 2000 (Bio-Rad, Hercules, CA). At least three PCR assays with different concentrations of a competitor (generally, decimal dilutions) were conducted for estimating a copy number of a target fragment. A relative amplification efficiency of a target to a competitor was estimated by PCR with known numbers of target and competitor fragments. The copy number of the target was estimated by considering the band intensities, length of the fragment, relative amplification efficiency of the target to the competitor, and copy number of the competitor as described by Lee et al. (22).
|
View this table: [in a new window] |
TABLE 3. Summary of qcPCR assays
|
Potential sulfur oxidation rate (SOR).
Incubation for determining SOR was initiated within 3 h after sampling. The sediment sample (20 g) was suspended in an equal weight of sulfate-depleted artificial seawater composed per liter of NaCl, 28.5 g; MgCl2 · 6H2O, 5.4 g; CaCl2 · 2H2O, 1.47 g; KCl, 0.725 g; SrCl2 · 6H2O, 0.024 g; NaBr, 0.084 g; H3BO3, 0.027 g; NaF, 2.9 mg; and KI, 0.1 mg (pH 7.2). This slurry was transferred to a bottle (100 ml) under air and supplemented with 0.4 ml of a Na2S · 9H2O solution (0.2 M; pH 7.0) and 0.4 ml of a 1,8-dihydroxyanthraquinone solution (1 mM in acetone; an inhibitor for sulfate-reducing bacteria [SRB]) (8). The bottle was gently shaken on a rotary shaker (50 rpm) at 20°C for 2 days, and the sulfate concentration was periodically measured by ion chromatography. Autoclaved sediment was incubated under the same conditions, which served as the abiotic control for estimating rates for chemical oxidation. The SOR was estimated by subtracting the chemical oxidation rate from the total oxidation rate.
Statistics.
Data were statistically analyzed by the Student t test (P = 0.05 unless otherwise stated) using the Koneko program (Kodansha Scientific, Tokyo, Japan).
Nucleotide sequence accession numbers.
The nucleotide sequences reported in this paper have been deposited in the DDBJ, EMBL, and National Center for Biotechnology Information nucleotide sequence databases under accession numbers AB116389 to AB116513.
|
|
|---|
16S rRNA gene clone library analysis.
Two sets of PCR primers (B341f/U1492r and U515f/U1492r) were used for amplification of prokaryotic 16S rRNA gene fragments, each of which successfully amplified a single fragment with the expected size from the sediment DNA. The16S rRNA gene clone libraries were constructed for the Kamaishi sediment obtained in August 2001 (the KM library) and Yamada sediments obtained in March 2001 (the Y3 library) and August 2001 (the Y8 library), each of which contained almost equal numbers of 16S rRNA gene clones derived from B341f/U1492r and U515f/U1492r-amplified fragments. Complete sequences of a total of 219 clones were determined (66 from the Y3 library, 79 from the Y8 library, and 76 from the K8 library), which produced 125 different phylotypes (a unique clone or group of clones with sequence similarity of >97%) (Table 4). Among them, 81 phylotypes (approximately 65%) comprised single clones, suggesting that the 16S rRNA gene clone libraries for the marine sediment prokaryotes were quite diverse.
|
View this table: [in a new window] |
TABLE 4. Phylotypes obtained in this study
|
-proteobacterial phylotypes and reference sequences. The phylogenetic relations appearing in Fig. 2 were confirmed by the maximum-likelihood algorithm (data not shown). Phylogenetic distributions of 16S rRNA gene clones in each of the Y3, Y8, and KM libraries at the phylum or subphylum level are presented in Fig. 3, in which clone distribution in the Antarctic library (4) is also shown for comparison. Although many phylotypes were present only in one library (Table 4), grouping analyses as shown in Fig. 3 revealed similarity in clone distribution between the Y3 and Y8 libraries, which was different from the KM library. In the Y3 and Y8 libraries, phylotypes affiliated with
-Proteobacteria were the most abundant, and those with
-Proteobacteria were the second. In contrast, the Kamaishi library was occupied evenly by phylotypes affiliated with Planctomycetes,
-Proteobacteria,
-Proteobacteria, and Crenarchaeota. Among them, those affiliated with Planctomycetes were the most abundant (18% in the KM library). In addition, the Kamaishi library contained many phylotypes affiliated with candidate divisions OP8 and OP9.
![]() View larger version (48K): [in a new window] |
FIG. 2. Neighbor-joining trees showing phylogenetic relationships among phylotypes obtained in this study and reference 16S rRNA gene sequences retrieved from the nucleotide sequence databases for -Proteobacteria (A) and -Proteobacteria (B). The bold vertical bar to the right of the tree in panel A represents the range of Thiotrichales, while the one to the right of the tree in panel B represents Desulfobulbaceae. The numbers at the branch nodes are bootstrap values per 100 trials; only values of >50 are shown. Accession numbers for the reference sequences are shown in parentheses.
|
![]() View larger version (26K): [in a new window] |
FIG. 3. Phylogenetic distribution of clones in each library. Y3, the Yamada March library; Y8, the Yamada August library; KM, the Kamaishi library (in August); Antarctic sediment, a 16S clone library for Antarctic continental shelf sediment (4).
|
-Proteobacteria were affiliated with Thiotrichales (15) and related to free-living and symbiotic sulfur-oxidizing bacteria (SOB), whereas the Kamaishi phylotype (KY142) was related to the genus Pseudomonas. In
-Proteobacteria (Table 4 and Fig. 2B), phylotype KY184 was the most abundant and included 10 clones that are distributed in all the three libraries; phylogenetic analysis indicated that this phylotype was closely related to the genus Desulofotalea (21), which includes psychrophilic SRB. Another phylotype found abundantly in the libraries was K82 (containing nine clones in the KM library); this phylotype and three other phylotypes (Y67, Y170, and Y162) were affiliated with Crenarchaeota and related to 16S rRNA gene clones in marine benthic group C (i.e., CRA9-27cm) (43). In addition, a significant number of phylotypes were related to the genus Pirellula in the Planctomycetes (38) (Table 4).
Quantification of rRNA gene copies by qcPCR.
The 16S rRNA gene clone library analysis detected differences in clone composition in libraries for the Kamaishi sediment (KM) and the Yamada sediment (Y3 and Y8) libraries; particularly, phylotypes related to potential sulfur oxidizers and those affiliated with
-Proteobacteria (potential sulfate reducers) were detected in abundance in the Y3 and Y8 libraries (Table 4; Fig. 2 and 3). To further evaluate these differences, qcPCR assays with group-specific primers were developed and are summarized in Table 2. The utility and limitation of qcPCR assays for quantifying microbial populations in natural ecosystems have been described previously (22, 44). To check the specificities of these PCR primers and optimize PCR conditions, we carried out PCR experiments using Escherichia coli cells harboring sediment rRNA gene clones (i.e., 125 phylotypes). The results confirmed that these PCR assays could amplify fragments with the expected sizes only from those clones harboring the target phylotypes as listed in Table 3. Primers B341f and U533r were used in the BAC qcPCR (Table 3) for quantifying the total bacterial population, and ratios in copy number (relative abundances) of the target phylotypes to the total bacterial population (the bacterial artificial chromosome copy number) were estimated. This method can prevent a possible incorrect estimate of the copy number resulting from variation in the DNA extraction efficiency.
Results of qcPCR for sediments obtained in July 2004 are presented in Fig. 4. The variations in relative abundance of SOB, SRB, and SVA (Table 3) among the Kamaishi sediment samples (K-A, K-B, and K-C) were small, as were the variations among the Yamada sediment samples in the aquaculture zone (Y-A, Y-B, and Y-C). However, the in-site variations in the DSV and DSS data were large. Comparison of the K-mean with the Y-mean in Fig. 4 shows that the SOB, SRB, and SVA values of the Y-mean were significantly higher than those of the K-mean. Similar trends in qcPCR were observed for sediments obtained in 2003 (data not shown).
![]() View larger version (34K): [in a new window] |
FIG. 4. qcPCR assays for quantifying rRNA gene copies of phylotypes related to SOB and SRB. Ratios of specific rRNA gene copies to the total bacterial rRNA gene copies are shown. The names of qcPCR assays refer to Table 3. Data are means (n = 3 samples for each sediment sample and n = 9 samples for K-mean and Y-mean), and error bars represent standard deviations.
|
![]() View larger version (27K): [in a new window] |
FIG. 5. SOR and SRR of marine sediment samples. Data are means (n = 3 samples for each sediment sample and n = 9 samples for K-mean and Y-mean), and error bars represent standard deviations.
|
|
|
|---|
-Proteobacteria and those related to potential SRB in
-Proteobacteria were more abundant in the Yamada sediment than those in the Kamaishi sediment. This finding was supported by several different lines of evidence, including the rRNA gene clone library analysis (Table 4; Fig. 2 and 3), the qcPCR assays (Fig. 4), and the activity measurements (Fig. 5). In addition, these differences were also observed between within (Y-A, Y-B, and Y-C) and outside (Y-D) the aquaculture zone. These results indicate that the sulfur cycle has become active in the Yamada sediment beneath the shellfish bed. The 16S rRNA gene clone library analysis (Table 4) suggests that sediment prokaryotes in the Yamada and Kamaishi bays are quite diverse (125 phylotypes from 219 clones sequenced). Similar levels of prokaryotic diversity in the marine sediment have been reported previously (4, 31). Ravenschlag et al. analyzed bacterial populations in permanently cold marine sediment in the Arctic Ocean, showing that 140 different ARDRA (amplified rRNA gene restriction analysis) patterns were obtained from 353 clones in a 16S rRNA gene library (31). Bowman and McCuaig found that 496 phylotypes were present in a subsample of clones (1,046 clones) in a 16S rRNA gene library constructed from Antarctic marine sediments (4). A rarefraction analysis conducted for the Arctic sediment library (31) has, however, suggested that the prokaryotic diversity in the marine sediment can only be partially assessed at this scale of analysis (with several hundred clones). Nevertheless, when the results of clone library analyses are compared with those of quantitative community analyses, e.g., FISH (32), rRNA slot blot hybridization (3, 32), and qcPCR (this study), it can be seen that the clone libraries contained 16S rRNA gene types representing numerically major groups of microbial populations. It is therefore possible to obtain a fundamental insight into the major prokaryotic populations by clone-library analysis on a similar scale.
In our clone libraries, KY184 was the most abundant
-proteobacterial phylotype (Table 4). This phylotype was closely related to the genus Desulfotalea that had been isolated from permanently cold Arctic marine sediments off the coast of Svalbard (21) and has been identified as psychrophilic SRB capable of incomplete oxidization of fatty acids (21). Related sequences have also been cloned from the Svalbard sediments (31), Antarctic shelf sediment (4), deep-sea sediments (23), and coastal sediments off Japan (42), indicating that this group of SRB is widely distributed in marine sediment. Our qcPCR assay using the SVA system (using a primer with the nucleotide sequence of probe Sval428) (35) estimated the relative abundance of this group (the SVA group) in the total bacterial population in the Sanriku sediments to be 1 to 19%. In addition, its relative abundance in the total SRB population was estimated to be 18 to 61%, indicating that it is the major sulfate-reducing population. A previous study has used the rRNA blot analysis for estimating this group of SRB in the Svalbard sediments (35), showing that they made up from 1.4 to 20.9% of the total eubacterial population. In the present study, it was also shown that the relative abundance of the SVA group (Fig. 4) corresponded to the SRR determined for the Yamada and Kamaishi sediments (Fig. 5). Based on these data, we concluded that SVA group SRB are the important sulfate-reducing population in the Sanriku coast sediment and suggest that their importance is not restricted to permanently cold marine sediments.
The qcPCR (Fig. 4) found that potential SOB made up from 9 to 16% of the total bacterial population in the Yamada bay sediments, while they only made up approximately 3 to 4% in the Kamaishi sediment. These values are compatible to estimates for SOB in the arctic Svalbard sediments (approximately 2% as determined by FISH using probe GAM660 and approximately 13% as determined by the rRNA blot analysis using probe GAM660) (32), while they are lower than the estimates for SOB in Antarctic shelf sediments (22 to 27% as determined by rRNA blot analysis with GAM660) (3). These data show that sulfur-oxidizing
-Proteobacteria detectable with probe GAM660 are important members of the prokaryotic community in the marine sediment, while the variation in their abundance can be ascribed to the use of different methods for the detection and estimation (47) and/or actual differences in the community structure. We should also consider a possibility that some SOB in marine sediments escaped detection by the SOB qcPCR assay. For instance, 16S rRNA gene sequences of Thiomicrospira spp. have one to several mismatches to GAM660; this genus represents aerobic lithotrophic SOB inhabiting marine sediments (5, 6). In addition, phylotypes Y121, Y189, and Y218 that have one or two mismatches when compared with GAM660 formed a coherent cluster together with the GAM660 match phylotypes (Y56, Y74, Y169, Y182, andY222) and SOB-related reference sequences (Fig. 2A). Further studies are thus needed to understand the total population structure of SOB in marine sediments.
A previous study showed that mussel aquaculture affected the macrobenthic community (40), demonstrating that dominant species were shifted from suspension feeders to deposit feeders. Such a difference in macrobenthos populations has not been observed in the Kamaishi and Yamada bays; the major benthos in these bays have been reported to be lugworms affiliated with Polychaeta (http://www.sui.pref.iwate.jp/
sui/frame2.htm). Other studies analyzed changes in compositions of biological molecules, e.g., phospholipids (28) and carbohydrates (9), in marine sediments in response to eutrophication. The results suggested that these biological signatures can be used as indices for the trophic states of the coastal environment. In the present study, we detected differences in the number and activity of sulfate reducers and sulfur oxidizers in marine sediments taken from two different bays, suggesting that these bacteria can serve as indices for assessing the organic load to the marine sediment. To further support this idea, comparative ecological surveys should be done, in which sediment samples will be taken from several different Sanriku bays with different levels of aquaculture and analyzed by the methods employed in the present study. In addition, it will also be necessary to examine contribution of other types of terminal electron-accepting processes (e.g., metal reduction) to organic matter decomposition. If this idea is fully supported by these analyses, sediment bacteria involved in the sulfur cycle (e.g., SVA PCR-detectable organisms) can be regarded as important indicators for assessing effects of aquaculture on the coastal environment in the Sanriku region.
This work was supported by the New Energy and Industrial Technology Development Organization (NEDO).
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»