Impacts of Ocean Acidification on Early Life-history Stages and Settlement of the Coral-eating Sea Star Acanthaster Planci

Coral reefs are marine biodiversity hotspots, but their existence is threatened by global change and local pressures such as land-runoff and overfishing. Population explosions of coral-eating crown of thorns sea stars (COTS) are a major contributor to recent decline in coral cover on the Great Barrier Reef. Here, we investigate how projected near-future ocean acidification (OA) conditions can affect early life history stages of COTS, by investigating important milestones including sperm motility, fertilisation rates, and larval development and settlement. OA (increased pCO 2 to 900–1200 matm pCO 2) significantly reduced sperm motility and, to a lesser extent, velocity, which strongly reduced fertilization rates at environmentally relevant sperm concentrations. Normal development of 10 d old larvae was significantly lower under elevated pCO 2 but larval size was not significantly different between treatments. Settlement of COTS larvae was significantly reduced on crustose coralline algae (known settlement inducers of COTS) that had been exposed to OA conditions for 85 d prior to settlement assays. Effect size analyses illustrated that reduced settlement may be the largest bottleneck for overall juvenile production. Results indicate that reductions in fertilisation and settlement success alone would reduce COTS population replenishment by over 50%. However, it is unlikely that this effect is sufficient to provide respite for corals from other negative anthropogenic impacts and direct stress from OA and warming on corals. Copyright: ß 2013 Uthicke et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Funding: The study was funded from appropriation funding of the Australian Institute of Marine Science. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


Introduction
Carbon dioxide (CO 2 ) concentrations in the atmosphere have increased by 40% over the past 250 years and are likely to double by the end of this century [1]. Increased atmospheric CO 2 leads to increased sea surface temperatures and a reduction in ocean pH, decreased carbonate and increased dissolved inorganic carbon (DIC) concentrations [2]. While some marine primary producers such as seagrasses, phytoplankton and macroalgae may benefit from increased DIC [3,4], many calcifying organisms exhibit reduced calcification due to lower saturation states of carbonate (e.g. [5][6][7]). Coral reefs in particular are threatened both by increased temperature (more frequent coral bleaching events) and ocean acidification (OA) [8]. Although susceptibility of corals to increased pCO 2 varies with species [9], it appears that many structurally complex corals will be lost, leading to a decline in habitat available to a variety of other species and changes to ecosystem structure and function [10].
In addition to the global pressures OA and climate warming, corals are under simultaneous threat from a variety of local stressors such as nutrient runoff, overfishing and cyclones that also dramatically shape the health of the system [11,12]. Episodic outbreaks of the corallivorous crown of thorns seastars (COTS, Acanthaster planci) reduce coral cover on many reefs in the Indo-Pacific region [13] and are a major contributor to recent coral decline in the Great Barrier Reef (GBR) [12]. The underlying causes of these outbreaks have long been debated [13]; however, higher larval survivorship caused by greater food (phytoplankton) availability, driven by agricultural land runoff, is currently the most widely accepted hypothesis to explain primary COTS outbreaks [14].
Sexual reproduction and early life history stages of marine invertebrates can be highly vulnerable to OA [15,16], with fertilization affected in some species but not others [17][18][19][20]. Calcifying larvae of echinoids (sea urchins) are vulnerable to reduced pH/increased pCO 2 with reduced calcification, abnormal growth and lower survival rates [15,21,22]. Similar results have also been reported for non-calcifying asteroid (sea star) larvae [17,23]. Reductions in settlement success due to OA have not been investigated in echinoderms, but have been documented in corals [24][25][26][27]. Here, we examine the effects of OA on sperm motility, fertilisation, larval development and larval settlement in the corralivorous sea star Acanthaster planci. This is the first study to comprehensively investigate several of these life history stages and also include settlement for an important marine species, and to provide a minimum estimate of the effects on population replenishment.

Results
Fertilisation rates depended on sperm concentrations under three different pH/pCO 2 conditions (Figure 1, see Table 1 for carbon chemistry parameters for all experiments). Functions for the three pH/pCO2 conditions tested had a similar slope (Table 2), but inflection points (ie the sperm concentration at which 50% fertilisation occurs) were significantly different among treatments, and increased with increasing pCO 2 ( Table 2, 'right-shift' in Figure 1). The percentage of fertilisation was reduced by .7% at pH NBS 7.9 (pCO 2 : 877 matm) and 25% at pH NBS 7.7 (pCO 2 : 1658 matm) across sperm concentrations spanning more than 3 orders of magnitude (,10 4 to 10 7 sperm ml 21 ). This reduction in fertilisation success is coincident with a significant reduction in sperm velocity and the percentage of motile (i.e. moving) sperm at elevated pCO 2 , with significant variation among individual males (Table 3, Figure 2). The coelomic fluid surrounding gonads of A. planci had a mean pH NBS of ,7.49 (N = 5, 95% CI: 7.39-7.61).
By 10 d, the majority of larvae in all treatments had reached late bipinnaria to early brachiolaria stages. The percentage of normal developing larvae declined by ,13% at pH NBS 7.6 compared with controls ( Figure 3) and this difference was significant (ANOVA, F 2,12 = 12.34, p = 0.0012; Tukey-Kramer posthoc test, p,0.05). Larval morphometry (e.g., length, stomach width) was not significantly different between treatments ( Figure 3; Table 4).
By 17 d, larvae were advanced brachiolaria stage, and a high proportion (75.4%) had a well-developed rudiment and were competent to settle. Two settlement experiments were conducted in parallel, achieving typically [28] moderate levels of settlement success of 10-20% in ambient conditions. In the first experiment settlement substrata known to induce settlement in COTS (CCAcrustose coralline red algae, and/or biofilms) that had developed under four pH conditions for 85 d were used. In the second experiment CCA was only exposed to ambient conditions prior to the settlement tests which were then conducted under the four different pH/pCO 2 conditions. There was no difference in the settlement success between CCA and biofilms (Table 5). In contrast, differences between pH/pCO 2 treatments were significant, with settlement on pre-exposed CCA and biofilms declining by about 50% from the pre-industrial treatment (pH NBS 8.1, pCO 2 = 350 matm) to the elevated pCO 2 (pH NBS 7.8, pCO 2 = 800 matm) treatment ( Figure 4). Settlement success for larvae offered CCA cultured under present day pCO 2 conditions but in seawater of the four pH levels was not significantly different (ANOVA, Figure 4); indicating that changes   Table 3). B) Differences in sperm velocity and percentage of motile sperm between individual males (Factor Male, see ANOVA in Table 3. Error bars represent 1 standard error. Averages with different letters are significantly different (p,0.05, Tukey-Kramer posthoc tests). doi:10.1371/journal.pone.0082938.g002 in settlement response was not due to the seawater pH/pCO 2 directly.
To evaluate the importance of individual early life history parameters for overall recruitment success, we calculated effect sizes and Bayesian confidence intervals for each parameter ( Figure 5). Environmentally relevant sperm concentrations were estimated using a sperm diffusion model and parameters given in Babcock et al. [29]. Diffusion model results indicated that eggs released from females 0.5-10 m downstream of 4 spawning males experience sperm concentrations in the range of 10 3 -10 4 ml -1 .
Under these conditions, fertilisation at pH NBS 7.9 and 7.6 was reduced by 29% and 75%, respectively ( Figure 5). Comparison of the effect sizes illustrated that the effect of pH/pCO 2 on fertilisation was amongst the most distinct of all parameters measured. While both the reduction in the percentage of motile sperm and sperm swimming speed influence fertilization, effect sizes suggest that motility (16 and 30% reduction) may be the driving force behind reduced fertilization success at low pH. Effects on sperm velocity were comparatively small (5 and 6%). OA effects on larval development and growth were subtle (11-  18% changes in the pH NBS 7.8 treatment) and it is uncertain if these effects are of sufficient importance for overall development to affect recruitment rates and population maintenance. However, slower development will lead to an extended period in the plankton, and daily mortality rates for echinoderm larvae due to predation and other factors can be 5-15% [30]. With an average reduction of settlement on biofilms and CCA of 36%, effects on settlement were roughly in the same range as for fertilisation. Even if only the effects of pH/pCO 2 on fertilisation and settlement are taken into account, overall recruitment under near future OA conditions may be reduced by <50% (1002(100*(120.29)*(120.36)) = 54.5%).  Container is a random/nested factor representing measurements taken from four independent culture containers. All data were log-transformed prior to analysis. doi:10.1371/journal.pone.0082938.t004  Table 1

Discussion
Near-future OA has the potential to severely impact population maintenance and growth of COTS by affecting several early life history stages including fertilisation and recruitment success. In the present study, depressed pH/elevated pCO 2 affected the percentage of motile sperm more distinctly than velocity, similar to previous experiments on echinoids [19,20]. Early work by Mortensen [31] already described that asteroid sperm requires pH elevation for activation. It is thought that decreased pH in the gonad suppresses activation of sperm mitochondria in echinoderms. Thus, sperm require elevated pH, such as that of the surrounding seawater, for activation of mitochondrial processes [18]. For example, the gonad of the sea star Patiria pectinifera has a pH of 6.4, and elevation of pH activates sperm [32]. Similarly, in the present study, the coelomic fluid of A. planci had a pH NBS of ,7.5, indicating sperm is also stored inactive at depressed pH. Previous studies on echinoids also show reductions the percentage of motile sperm at decreased pH and ultimately reproductive success [18][19][20]. A second mechanism for reduced fertilisation success at depressed pH involves reduced efficiency of the fast sperm block for polyspermia [18], although this would be evident as a drop off in fertilisation at higher sperm concentrations, an outcome we did not observe under the sperm concentration tested. Significant variation in sperm velocity and motility amongst males was also observed in sea urchins [20], and may provide a mechanism for future selection and potential adaptation of the species to increased pCO 2 . Recent population genomic [33] and quantitative genetic [34][35][36] studies on sea urchins suggest a potential for adaptation to OA through selection in that class of the Echinoderms. Whether this is also the case for COTS should be subject to further investigation.
Ocean acidification significantly reduces the growth of calcifying echinoderm larvae (e.g. [15,21,22,37,38]). Interestingly, noncalcifying asteroid larvae, such as A. planci, are also sensitive to OA; for example, increased mortality and a stunted growth were recently reported in a temperate and an Antarctic asteroid species [17,23]. The underlying mechanism(s) responsible for OA-induced mortality and/or depressed growth of non-calcifying larvae is unclear, but it is possible that direct teratogenic effects of decreased pH stunted growth or that increased metabolism associated with acid/base regulation reduces energy availability, or 'scope' for growth [39] [40].
The present study is the first to experimentally demonstrate effects of OA on settlement of invertebrates other than corals. Both temperature and OA can alter biofilm and/or crustose coralline algal (CCA) communities that are critical for inducing settlement Table 5. Analysis of variance (type III sums of squares) for settlement of COTS larvae on CCA and biofilms exposed for 85 d in water with pre-industrial, present day control or two near future pCO 2 conditions (see Table 1).  of coral larvae [24,26,27,41]. Microbial studies demonstrated that OA alters the chemistry and bacterial community composition of settlement substrata [24,42], but it is unclear which factor is ultimately responsible for reductions in larval settlement. Larval settlement success observed here on CCA not exposed to decreased pH was in the same range as those measured previously on the same algae [28]. That earlier study also demonstrated that addition of antibiotics reduced settlement, suggesting bacterial biofilms as possible agents of settlement induction. This is consistent with our observation that settlement success was similar on CCA and biofilms. Our second settlement experiment illustrated that reduced larval settlement under OA conditions (depressed pH/elevated pCO 2 ) is not simply the direct result of sea water chemistry. These findings are consistent with those for coral larvae [24,27]. Detecting similar trends in larval settlement in representatives from different tropical phyla, and given similar cues are important for the recruitment of urchin [43] and abalone [44] larvae, reduced settlement through OA altering settlement cues could be an ecologically important phenomenon also impacting polar and temperate species. It is intriguing that projected ocean acidification scenarios have similar effects on early life history stages of prey (coral) and predator (COTS), despite representing widely divergent phylogenetic groups. Fertilisation curves derived here for COTS are similar to those observed for corals [26,45], and settlement experiments yield similar results (see above). Our model calculations suggests that population replenishment of COTS could be reduced by at least 50% if pCO 2 increases as expected by 2100. This estimate may be conservative, as it does not include additional unknown impacts of OA on juveniles and adults. However, the estimate also assumes populations are recruitmentlimited as assumed for most open marine populations [46]. In addition, this estimate does not account for potential acclimation/ adaptation in COTS gametes and developmental stages that may reduce effects of OA (see above). Although this figure is high, it must be interpreted in the context of other environmental parameters known to influence population development in this species. For example, a small increase in food availability for COTS larvae caused by terrestrial runoff can cause an 8-fold increase in larval survivorship [14]. The same study also showed that higher food availability can increase larval size whereas OA had no detectable effect on larval size in the present study. Together, these studies indicate that while OA may negatively affect population growth of this coral predator, improved understanding of the drivers of these outbreaks and potential management strategies are essential to control future outbreaks across the Indo-Pacific region. In addition, it is unlikely that reductions in COTS population size would take effect before corals themselves are impacted by OA or increased sea surface temperature [8,9].

Ethics statement
All experiments were conducted in accordance with Australian laws and specimens collections were approved by the Great Barrier Reef Martine Park Authority (Permit No. G12/35236.1).
All experiments were conducted at 28uC which represents seawater temperatures in the source area of the A. planci (Green Island, Northern Great Barrier Reef, 16u46.59S, 145u59.39E) around spawning time (,November to January).

Fertilisation assay
Ovaries and testes were dissected from the base of the animals' arms. Testes spontaneously released sperm which was collected dry. Post-vitellogenic oozytes were obtained by treating dissected ovaries in 10 25 M 1-methyladenine in FSW [47].To represent a population response, fertilisation assays were conducted with the combined sperm from four males and oocytes from five females.
Experiments were conducted in 20 ml vials containing filtered (5 mm) natural sea water (salinity: 35.5) at the respective treatment pH. Dry sperm (0.1 ml) was diluted in 10.9 ml treatment water in each of 10 replicate vials for each pCO 2 treatment to achieve an approximate sperm concentration of 10 7 -10 8 sperm ml 21 . A 1 ml sample was taken from each vial for sperm counts in a haemocytometer. We then diluted 1 ml of each replicate with 9 ml of water in six serial dilutions. To accurately fit statistical models, we calculated the sperm concentration in each vial as the measured concentration of the corresponding undiluted replicate multiplied by the dilution factor. We added 1 ml of oocyte suspension to each vial to achieve a final volume of 10 ml, and a concentration of 100 oocytes ml 21 . After 10 minutes 2 ml were sampled from each vial, placed in an EppendorfH tube and development stopped with 1 drop of 7% formalin. 50 eggs in each replicate were checked, the number of eggs with fertilization envelopes amongst those counted, and the percentage of fertilization calculated. Ambient pH NBS (8.1), pH NBS 7.9 and 7.6 were used as treatments representing present day and future (2100 to 2300) pH and pCO 2 concentrations under different scenarios [48]. pH was measured using a temperature-corrected pH meter (OAKTON, USA; pH probe: EUTECH, USA), and controlled using a Tris standard (Dixon, Batch 5). Treatment sea water was prepared by slowly bubbling pure CO 2 through a 10 l container of filtered seawater until the nominal target pH was reached. Stability of the pH was measured over the following 15 min, and for each treatment samples for alkalinity and dissolved organic carbon (DIC) analysis (250) were fixed with 125 ml of saturated mercuric chloride. Alkalinity/DIC Samples were analysed at 24uC using a VINDTA 3C titrator (Marianda, Germany). Calibration was conducted using Certified Reference seawaters (A. G. Dickson, Scripps Institute of Oceanography, Dixon, Batch 106).
We measured the pH of the coelomic fluid around the gonads of five A. planci individuals by carefully removing fluid around the gonads using a syringe. pH was measured directly after fluid removal with a pH probe as described above.

Sperm velocity and motility, and pH of the coelomic fluid
Sperm point-to-point velocity (Velocity Curvi-Linear) and percentage of motile sperm (referred to as motility) were measured from 6 male A. planci, using similar techniques as were previously used for echinoids [20]. Sperm from each individual was kept at 28uC and assays conducted within 30 min. Control and treatment water was prepared in the same way as for the fertilisation assay. For each dilution, 2 ml of dry sperm were diluted with 4 ml of treatment water in a separate 2 ml scintillation vial.
Sperm dilutions were taken up in 0.3 mm capillary slides for microscopy. All slides were thoroughly washed and rinsed 3 times and coated in bovine serum albumin (1%) to avoid sperm sticking to the slides. We used a Zeis axioscope with 200X magnification for microscopic filming. The camera (Pixilink PL-B623) was set to take 25 frames per second over a two second period. Focus was set mid-plane to minimize wall effects on swimming speeds. From each slide at least 5 videos were taken to film a minimum of 200 sperm. Three slides from independent sperm dilutions were taken for each male x treatment combination.
For velocity analysis, we used ImageJ with the CASA (computer assisted sperm analysis) plugin. Motility was scored visually from the first 5 digital frames of each video counting all mobile and immobile sperm.

Larval development
Larval development experiments were conducted in a flowthrough seawater system (,180 L h 21 ) with UV sterilised and filtered (1 mm) water. The rearing containers were 100 ml plastic jars with a window that maintained a constant volume of 30 ml. A 45 mm mesh screen was glued in front of the widow to retain larvae. Experimental pH was regulated by injection of CO 2 into the seawater reservoirs using an automatic CO 2 injection system. The pH in sections of the system was regulated according to water chemistry conditions in the rearing containers with two pH controllers (Tunze), set at pH NBS , 7.6 and pH NBS 7.8 with a third section allowed to track ambient pH.
Eggs from 4 females were pooled and placed in each rearing containers (N = 4) and allowed to acclimatise for a 15 min. Sperm of three males were pooled and added to achieve a sperm to egg ratio of 200:1; 5610 3 sperm ml 21 . The flow of seawater was turned off for 10 min to allow fertilisation. From day three larvae were fed three times a day with 1-15610 6 cells of Proteomonas sulcata for each container. At day 10 the larvae were harvested. The first 30 larvae were haphazardly sampled and scored for normal development. Normally developed larvae were photographed and their total length along the midline from posterior to the anterior, their maximum width and the length of the stomach from the end of the oesophageal tube to the beginning of the hind gut was measured.
For settlement experiments, larvae produced from the same males and females were reared in 300 l aerated containers (,5 larvae ml 21 ). Larvae were fed daily with P. sulcata at approximately 10 5 cells ml 21 .

Larval settlement
Crustose coralline red algae (CCA) of the species Hydrolithon onkodes were collected from Davies Reef (18u49.29S, 147u37.99E) in August 2012. That species was chosen because it was identified as a 'refuge substratum' for juvenile A. planci [49]; although one deep water species can induce higher settlement rates [28]. CCA chips (,1 cm diameter) were embedded in non-toxic under water glue (Mr. Sticky'sH, Fair Oaks, CA) to avoid open carbonate surfaces and glued to PVC slides. We also prepared slides only with the glue to test for potential toxic effects. Slides were mounted in custom perspex holders which were held in place on aquarium walls using magnets. The aquarium system used was similar to the system described for larval experiments. However, pH and temperature were monitored continuously (30 sec sampling rate) with ISFET type pH probes (Endress Hauser CPS-471D). Fresh filtered seawater (0.4 mm) was added to the system at a flow rate of 1 l min 21 for a 200% replacement rate per day. Flow rates in each experimental tank were 12 l min 21 . In addition to a present day (pH NBS 8.1 target), medium (future pH NBS 8.0 target) and low (future pH NBS 7.9 target) acidification treatments, this experiment included a pre-industrial treatment (past pH NBS 8.25 target). Acidified treatments were achieved by bubbling CO 2 into sump tanks with solenoid valves (SMC pneumatics) and pH setpoints, while the pre-industrial treatment was achieved by passing a stream of air through 2 sodalime canisters and mixing the low CO 2 scrubbed air with the incoming seawater in a counter current exchange tower prior to flowing into each experimental tank. Temperatures were controlled with a heater chiller unit (EvoHeat DHP40) creating a recirculating water bath for temperature control. Seasonal temperature profiles were maintained with an average of 26.1uC60.2 SD (Max = 27.9uC & Min = 24.0uC).
Settlement assays we conducted in 6 well tissue culture dishes (TCD). Each well contained 10 ml of FSW of the respective treatment. CCA were removed from the slide holders and carefully placed into the wells. We added 20-30 competent larvae to each well. Competent COTS larvae are late stage brachiolaria with a well-developed rudiment. Results from some wells which had a lower count of competent larvae at the end of the experiments were omitted from the analysis. The first experiment contained CCA pre-incubated in one of the 4 pCO 2 treatments for 85 d and we ran 18 replicates for each treatment. In addition to the CCA we used biofilms grown on the glue for the same amount of time as additional settlement substratum. Six replicates were employed for each treatment.
To distinguish effects caused by changes in cues the CCA and biofilms growing under different pCO 2 concentrations from those potentially caused by the differences in seawater chemistry, a second type of settlement experiment was conducted. In this experiment, all CCA (N = 12) used were grown under ''present day'' conditions. During the settlement period, these CCA and the larvae were inserted into the wells together with water from the 4 different treatments. Both types of experiments were run in parallel with larvae of the same age and maturity. In order to minimise re-equilibration of the CO 2 in treatment water to present day conditions, 6 well plates were closed with a lid, and floated individually in closed 400 ml plastic containers containing 200 ml of treatment water.

Statistics
Fertilisation rates at different sperm concentration were modelled using a two-parameter logistic model (asymptote fixed to 100), with log-transformed sperm concentrations. Standard errors for the predicted curves where estimated using the delta method, a method used in statistics to derive approximate probability distributions. These standard errors were used to calculate confidence intervals.
A total of 38,121 sperm were analysed for velocity and motility. We used the average velocity per replicate within each slide for statistical analysis. Mixed model Analyses of Variance (ANOVA) was used to analyse both sperm velocity and motility. pH/pCO 2 treatments and Males were treated as fixed factors, whereas replicate dilutions (slide) were treated as a random/nested factor. Sperm velocity was log transformed and motility arcsine-square root-transformed prior to analysis.
Data for development were analysed with a one-factor (% normal development, arcsine-square-root transformed) or mixed model (larval size, rearing container used as nested factor) ANOVA.
Data for larval settlement (% settlement of competent larvae) were arcsine-square root-transformed and subjected to two-factor ANOVA with Treatment (4 different pH/pCO 2 levels) and CCA vs biofilm as fixed factors. Type III sums of squares were used for this ANOVA because removal of some replicates with low numbers of competent larvae (see above) and unequal sample size between CCA and biofilm-controls lead to an unbalanced design. One way ANOVA (also Type III) was used to for the water-control experiment. Tukey-Kramer tests were choses as post-hoc tests for all ANOVA if indicated. Transformed data for all ANOVA conformed to normality and equal variance assumptions.
To compare the effect size of all individual parameters, the change in response was expressed as a precent change compared to control (present day) conditions. Based on the standard deviations of the control and respective treatment, we calculated a Bayesian 95% confidence interval for the reduction using Markov Chain Monte Carlo sampling. After a burn-in period of 2000 steps, 5000 steps were sampled and three parallel chains run.
Statistical analyses were conducted in R [50], and mixed model ANOVAs conducted in NCSS [51].