Contribution to Zoology, 75 (3/4) – 2006
Phenotypic plasticity revealed by molecular studies on reef corals of Fungia (Cycloseris) spp. (Scleractinia: Fungiidae) near river outlets
Adriaan Gittenberger , Bert W. Hoeksema
Keywords: eco-phenotype; plasticity; allozymes; river outlets; mushroom corals; Fungiidae; Cycloseris; Indonesia
On a patch reef off Makassar, Sulawesi, Indonesia, corals identified as Fungia (Cycloseris) costulata, Fungia (Cycloseris) tenuis and Fungia (Cycloseris) cf costulata were collected down to a maximum depth of 10 m. The corals lived sympatrically. Mushroom coral clones resulting from fragmentation can be recognized by their equal coloration and close proximity. Therefore, to ensure that no clones were collected, corals of dissimilar colors were selected at a mutual distance of 5 m. The corals were kept alive in two 30 liter sea-water aquariums with an air-pump. They were photographed in detail. Using allozyme electrophoresis in a laboratory close to the field area, it was tested whether the separate coral morphs should be considered three species. Eventually it was concluded that there are only two species, i.e. F. (C.) costulata and F. (C.) tenuis, of which F. (C.) costulata has two distinct morphs, one of which may be an eco-phenotype occurring on reefs off river outlets or inside estuaries.
Intraspecific variation in scleractinian corals is a classic problem in their taxonomy, both in regard to recent and fossil species (Best et al., 1999; Knowlton and Budd, 2001). Habitat-induced variability has been observed in coral species distributed along depth ranges, in which the specimens from the deeper sites used to be flatter than those from shallower places. The former ones exposing more surface area in order to compensate for less light penetration at greater depths (Wijsman-Best, 1972, 1974; Hoeksema, 1993). Other environmental factors of importance in coral shape plasticity may be sedimentation, salinity and water temperature, turbulence and flow; in addition to coral morphology also the pigmentation in the soft tissue may be affected (Bruno and Edmunds, 1997; Todd et al., 2002a, 2002b, 2004b, 2004c). Coral damage and subsequent regeneration also have specific effects on shape, which may impede taxonomy and easy identification (Hoeksema, 1989, 1991b, 1993; Oren et al., 1997; Nagelkerken and Bak, 1998).
For no obvious reason, most of the research on intraspecific coral variation has been performed on massive faviids, both in the Atlantic (Dustan, 1975; Foster, 1977, 1979; Lasker, 1981; Graus and Macintyre, 1982; Dodge, 1992[b; Beltran-Torres and Carricart-Gavinet, 1993; Amaral, 1994; Manica and Carter, 2000) and the Indo-Pacific (Wijsman-Best, 1972, 1974; Miller, 1994; Oren et al., 1997; Todd et al., 2001, 2002a, 2002b, 2004a, 2004b, 2004c). Mushroom corals (Fungiidae), which are endemic to the Indo-Pacific, have mainly been studied in relation with light penetration, sedimentation, and traumatic damage (Hoeksema, 1989, 1991b, 1993; Hoeksema and Moka, 1989). A physiological difference with regard to vulnerability to elevated temperature has also been shown among individuals of the same mushroom coral species within close range of each other (Hoeksema, 1991a).
During earlier taxonomic and morphological studies on mushroom corals three different morphs of Fungia (Cycloseris) spp. were distinguished (Hoeksema, 1989; Hoeksema and Moka, 1989). Two of these were considered separate species, viz. Fungia (Cycloseris) costulata Ortmann, 1889 (figs. 1-2), and F. (C.) tenuis Dana, 1864 (figs. 7-8); the third one was seen as a morph of the first (figs. 3-6; Hoeksema and Moka, 1989: fig. 12). Fungia (C.) tenuis has much rougher costae and usually a slightly different coloration (dark brown stomatal ends) as compared to F. (C.) costulata. The two are usually observed on the same reefs. However, F. (C.) costulata may also be present on reefs that are more nearshore, in more sediment-rich water, and they may occur deeper when found on the same reef as F. (C.) tenuis. The alleged separate morph of F. (C.) costulata is thinner and shows an evenly brown-olive green color (figs. 3-6). It occurs on nearshore reefs, either on patch reefs near river outlets or inside deep bays and estuaries. The three morphs cannot always be distinguished easily. Therefore, in order to investigate the taxonomic implications of the morphological differences, several allozymes of the morphs were compared. Allozyme electrophoresis has successfully been used to solve similar taxonomical problems (Gittenberger et al., 2001; Sanjuan et al., 1997). The three morphs will be referred to as morphs A, B and C for respectively Fungia (Cycloseris) costulata, F. (C.) tenuis and F. (C.) cf costulata.
Figs. 1-6. Upper and lower surfaces of mushroom corals. 1-2, Fungia (Cycloseris) costulata. 3-6, Fungia (Cycloseris) cf costulata. Scale = 1:1.
Material and methods
Specimens of morphs A, B and C were found sympatrically at Bone Baku reef, off Makassar, Sulawesi, Indonesia. Twelve specimens of each morph were collected within an area of about 200 m2 at depths between 2 and 12 meters. They were coded A1-12, B1-12 and C1-12. To make sure that no clones were included, only individuals that differed in polyp coloration were selected. While diving, the specimens were individually put into separate plastic bags and transported to the laboratory in a bucket with seawater. They were kept alive in two 30 liter aquariums with two air-pumps each. To reduce the pollution in the aquariums, the water was taken a few kilometers off the coast and filtered through a coffee filter. Before the allozyme electrophoresis, all corals were digitally photographed on both sides with a Fujifilm MX-2700 camera.
Each coral was taken out of the aquarium. After that half a 1.5 ml test-tube of coral tissue mixed with small pieces of skeleton was rasped of the septae with a scalpel, and 0.050 ml of homogenizing buffer (0.01M Tris, 0.001 M NaEDTA, 0.01 M Maleic acid and 0.001 M MgCl2) was added. The mixture was ground with a micro-pestle and put on ice. The damaged specimens were digitally photographed on both sides and conserved in 96% alcohol, as reference material. To create a centrifuge, the blades of a small table-ventilator were removed and the tubes were stuck to the spindle with heavy duty tape. Each sample was centrifuged for 30 seconds at maximum speed. The supernatants were extracted with a 0.100 ml pipette and added to a new tube, which was centrifuged for 30 seconds and put on ice.
Occasionally the supernatant was too slimy (highly viscose and sticky) to be extracted into a 0.100 ml pipette point. In that case, the top of this pipette point was cut off with a scalpel to enable the extraction of the “slimy” supernatant into a new tube. An additional 0.050 ml homogenization buffer was added and everything was mixed by sucking it up and down into the cut pipette point. After centrifuging this mixture at maximum speed for 30 seconds, the supernatant could be extracted with a 0.100 ml pipette point. It was added to a new tube and put on ice.
The slots of a well-plate were filled with 0.010 ml supernatant each. An applicator was used to load and apply the supernatant to a cellulose acetate gel (Hillis et al., 1987). Supernatants of all samples were run on a gel for 25 minutes and on an additional gel for 40 minutes. The electrophoresis was performed in a refrigerator at 4°C.
To test which allozyme systems in combination with which running buffers work best for Fungiidae, i.e. show polymorphic loci with a good resolution and activity, a preliminary study was done in the Netherlands. In total 10 allozyme systems were tested on two specimens of Heliofungia actiniformis (Quoy and Gaimard, 1833) and two specimens of Fungia (Verrillofungia) repanda Dana, 1846, viz. apartate aminotransferase (sAAT 18.104.22.168), alcohol dehydrogenase (ADH, 22.214.171.124), glucose dehydrogenase (GCDH, 126.96.36.199), glucose 6 phosphate dehydrogenase (G6PDH, 188.8.131.52), hexokinase (HK, 184.108.40.206), L-iditol dehydrogenase (IDDH, 220.127.116.11), isocitrate dehydrogenase (IDH, 18.104.22.168), malate dehydrogenase (MDH, 22.214.171.124), glucose-6-phosphate isomerase (GPI, 126.96.36.199), and phospho-glucomutase (PGM, 188.8.131.52). The nomenclature and IUBNC numbers are according to those of the standard of the International Union of Biochemistry (IUBNC, 1984). All allozyme systems were tested in combination with three buffers (Saccheri, 1995), i.e. TG (pH 8.5, 25mM Tris, 192 mM Glycine), TM (pH 7.8, 50 mM Tris, 20 mM Maleic acid) and P (pH 7.0, 11.6 mM Na2HPO4.2H2O, 8.4 mM NaH2PO4.H2O). The loci of the allozyme systems G6PDH, HK, MDH and GPI were polymorphic, with a good activity and a reasonable resolution using TM, TM, TG and P buffer respectively. The chemicals to test these allozyme systems were either bought in Indonesia or imported ice-packs from the Netherlands. These allozyme system-buffer combinations also gave good results during an additional test in Indonesia for specimens of Herpolitha limax (Esper, 1797), Zoopilus echinatus Dana, 1846, Fungia (Dana-fungia) fralinae Nemenzo, 1955, F. (D.) scruposa Kluzinger, 1879, F. (Pleuractis) gravis Nemenzo, 1955, F. (Verrillofungia) repanda Dana, 1846 and F. (V.) scabra Döderlein, 1901. A spider extract was used as a reference and positive control in all analyses. It showed a good activity and resolution for the allozyme system-buffer combinations described above. Unexpectedly, none of the specimens of morphs A, B or C showed any clear bands. Therefore, the allozyme system-buffer combinations were tested again for these corals. A good activity and reasonable resolution was only seen for the PGI allozyme system in combination with TG buffer. It showed one polymorphic and one homomorphic locus and was studied for the three morphs. The spider extract showed 4 bands with a high activity and good resolution and was used as a reference.
The resulting bands for each specimen were scored independently on two gels which had run for 20 and 40 minutes respectively. The bands that were scored twice were used for further analysis. The package of Swofford and Selander (1981), BIOSYS-1, was used to analyse the data. The exact probability test was used to test for Hardy-Weinberg at the polymorphic locus.
All the corals used in the experiment were investigated and identified morphologically in the field and, independently, from the photographs.
Except for the morphological differences described in the introduction, two additional characters distinguishing between A and B on the one hand, and C on the other hand were noticed.
It took about 2 minutes per specimen of the morphs A and C to scrape off sufficient tissue mixed with septal skeleton pieces to fill half a 1.5 ml tube, while it only took about 15 seconds for each specimen of B, indicating that the skeletal structure of B was weaker. Furthermore, all specimens of only the morphs A and C had slimy mucous for which the protocol had to be adjusted (see “Material and methods section”).
The frequencies of the 6 alleles found for the polymorphic locus of PGI are shown in table 1. The alleles A, D and F, accounting for 77% of all alleles scored for morph B, are not present in A and C.
Fig. 9. UPGMA dendrogram based on Rogers’ genetic distances for the three samples (cophenetic correlation = 0.99).
None of the three samples A, B and C (respectively p = 0.57, 1.00 and 0.48) was significantly (α = 0.05) deviating from Hardy-Weinberg. However, when pooling the samples, the frequencies of the alleles were significantly (p = 0.02) deviating for the sample A+B+C, almost significantly (respectively p = 0.06 and p = 0.12) for A+B and B+C, and not (p = 0.35) for A+C. This indicates that the samples A, B, C and A+C can each be considered representatives of single demes, while the samples in which morph B was pooled with the morphs A and C cannot.
Nei’s genetic distances (D), ranging between 0.02 and 0.32, and Rogers’ genetic distances, ranging between 0.087 and 0.317, are shown in table 2. The dendrogram (fig. 9) resulting from an UPGMA on Rogers’ genetic distances has a very high cophenetic correlation (0.99), indicating that it accurately reflects the pattern of genetic variation in the matrix of genetic distances (Sneath and Sokal, 1973). An UPGMA using Nei’s genetic distances gave similar results.
The fact that the allozyme system buffer combinations that work best for Fungia (Danafungia) fralinae, F. (D.) scruposa, F. (Pleuractis) gravis, F. (Verrillofungia) repanda and F. (V.) scabra do not give any clear results for the Fungia (Cycloseris) specimens, could be an indication that the allozymes of the latter taxon differ considerably from those of the former taxa. Therefore it might be more appropriate to refer to Cycloseris as a genus. DNA-analyses of Fungiidae also support this view (Gittenberger et al., 2006).
A significant deviation of Hardy Weinberg was found when the alleles of all the specimens were pooled. This indicates that they should not be considered representatives of a single panmictic population. No proof was found for a reproduction barrier between F. (C.) costulata offshore morphs and F. (C.) costulata nearshore morphs. Pooling their alleles, no deviation of Hardy Weinberg was found. They should therefore be referred to as phenotypes within F. (C.) costulata.
Nei’s and Rogers’ genetic distances and the resulting dendrogram (fig. 9) clearly show that there is very little to no gene flow between F. (C.) tenuis and F. (C.) costulata at Bone Baku reef. In total 77% of all alleles scored for F. (C.) tenuis were not present in the specimens of both forms of F. (C.) costulata, and vice versa 95% of the alleles scored for F. (C.) costulata were only accounting for 9% in F. (C.) tenuis. These results combined with the morphological differences, i.e. the roughness of the costae, the mouth coloration, skeletal strength and sliminess of the mucus, support the view that Fungia (Cycloseris) costulata Ortmann, 1889, and F. (C.) tenuis Dana, 1864, are two valid species; and that F. (C.) costulata has a nearshore ecomorph that may be related to low-salinity sea water.
We would like to thank drs P. v. Bragt from the Hoge Laboratorium School, Breda, the Netherlands, for providing the space and equipment to do the preliminary allozyme electrophoresis study in the Netherlands. Dr A. Noor is thanked for his help concerning the permits and facilities enabling research off Makassar, Indonesia. This study was supported by WOTRO (grant nr. W 82-249).
Amaral FD. 1994. Morphological variation in the reef building coral Montastrea cavernosa in Brasil. Coral Reefs13: 113-117.
Best MB, Boekschoten GJ, Oosterbaan A. 1999. Species concept and ecomorph variation in living and fossil Scleractinia. Paleonotographica Americana 54: 70-79.
Beltran-Torres AU, Carricart-Gavinet JP. 1993. Skeletal morphological variation in Montastrea cavernosa (Cnidaria: Scleractinia) at Isla Verde Coral Reef, Veracruz, Mexico. Rev. Biol. Trop. 41: 559-562.
Bruno JF, Edmunds PJ. 1997. Clonal variation for phenotypic plasticity in the coral Madracis mirabilis. Ecology 78: 2177-2190.
Dodge RE. 1982. Effects of drilling mud on the reef-building coral Montastrea annularis. Mar. Biol. 71: 141-147.
Dustan P. 1975. Growth and form in the reef-building coral Montastrea annularis. Mar. Biol. 33: 101-107.
Foster AB. 1977. Patterns of small-scale variation of skeletal morphology within the scleractinian corals Montastrea annularis and Siderastrea siderea. Proc. 3rd Int. Coral Reef Symp., Miami 2: 409-415.
Foster AB. 1979. Phenotypic plasticity in the reef corals Montastrea annularis and Siderastrea siderea. Proc. 3rd Int. Coral Reef Symp., Miami 2: 409-415. J. Exp. Mar. Biol. Ecol. 39: 25-54.
Gittenberger A, Reijnen BT, Hoeksema BW. 2006. A molecular analysis of the evolutionary history of mushroom corals (Scleractinia: Fungiidae) and its consequences for taxonomic classification. Chapter 3 in thesis, The evolutionary history of parasitic gastropods and their coral hosts in the Indo-Pacific, ISBN 90-9021100-4, Leiden: 35-56.
Gittenberger A, Vrieling K, Gittenberger E. 2001. Restricted gene flow between two alleged subspecies of Albinaria cretensis (Gastropoda, Pulmonata, Clausiliidae). Neth. J. Zool. 51 (1): 71-84.
Graus RR, Macintyre IG. 1982. Variation in the growth forms of the reef coral Montastrea annularis (Ellis and Solander): a quantitative examination of growth response to light distribution using computer simulation. Smith. Contr. Mar. Sci. 12: 441-464.
Hillis DM, Rosenfield DS, Sanches M. 1987. Allozymic variability and heterozygote deficiency within and among morfologically polymorphic populations of Liguus fasciatus (Mollusca:Pulmonata:Bulimulidae). Am. Malacol. Bull. 5: 153-157.
IUBNC. 1984. International Union of Biochemistry: Nomenclature Committee. 1984. Enzyme Nomenclature, 1984. Academic Press, Orlando, Florida.
Hoeksema BW. 1989. Taxonomy, phylogeny and biogeography of mushroom corals (Scleractinia: Fungiidae). Zool. Verh., Leiden 254: 1-295.
Hoeksema BW. 1991a. Control of bleaching in mushroom coral populations (Scleractinia: Fungiidae) in the Java Sea: stress tolerance and interference by life history strategy. Mar. Ecol. Prog. Ser. 74: 225-237.
Hoeksema BW. 1991b. Evolution of body size in mushroom corals (Scleractinia: Fungiidae) and its ecomorphological consequences. Neth. J. of Zool. 41: 122-139.
Hoeksema BW. 1993. Phenotypic corallum variability in Recent mobile reef corals. Cour. Forsch. Inst. Senckenberg 164: 263-272.
Hoeksema BW, Moka W. 1989. Species assemblages and ecomorph variation of mushroom corals (Scleractinia: Fungiidae) related to reef habitats in the Flores Sea. Neth. J. of Sea Res. 23: 149-160.
Knowlton N, Budd AF. 2001. Recognizing coral species past and present. In: J.B.C. Jackson, S. Ligard and F.K. McKinney (eds.), Evolutionary patterns: Growth, form, and tempo in the fossil record. Univ. Chicago Press, Chicago: 97-119.
Lasker HR. 1981. Phenotypic variation in the coral Montastrea cavernosa and its effect on colony energetics. Biol. Bull. 160: 292-302.
Manica A, Carter RW. 2000. Morphological and fluorescence analysis of the Montastrea annularis complex. Mar. Biol. 137: 899-906.
Miller KJ. 1994. Morphological variation in the coral genus Platygyra: environmental influences and taxonomic implications. Mar. Ecol. Prog. Ser. 110: 19-28.
Nagelkerken I, Bak RPM. 1998. Differential regeneration of artificial lesions among sympatric morphs of the Caribbean corals Porites astreoides and Stephanocoenia michelinii. Mar. Ecol. Prog. Ser. 163: 279-283.
Oren U, Benayahu Y, Loya Y. 1997. Effect of lesion size and shape of regeneration of the Red Sea coral Favia favus. Mar. Ecol. Prog. Ser. 146: 101-107.
Saccheri IJ. 1995. An experimental study of the effects of population bottlenecks on genetic variation and fitness in the butterfly Bicyclus anynana (Satyridae) Thesis. University of London, London.
Sanjuan A, Pérez-Losada M, Rolan E. 1997. Allozyme evidence for cryptic speciation in sympatric populations of Nassarius spp. (Mollusca: Gastropoda). J. Mar. Biol. Ass. U.K. 77: 773-784.
Sneath PHA, Sokal RR. 1973. Numerical taxonomy. Freeman, San Francisco.
Swofford DL, Selander RB. 1981. BIOSYS-1: A FORTRAN program for the comprehensive analysis of electrophoretic data in population genetics and systematics. J. Hered. 72: 281-283.
Todd PA, Sanderson PG, Chou LM. 2001. Morphological variation in the polyps of the scleractinian coral Favia speciosa (Dana) around Singapore. Hydrobiologia 444: 227-235.
Todd PA, Sidle RC, Chou LM. 2002a. Plastic corals from Singapore: 1. Coral Reefs 21: 391-392.
Todd PA, Sidle RC, Chou LM. 2002b. Plastic corals from Singapore: 2. Coral Reefs 21: 407-408.
Todd PA, Laddle RJ, Lewin-Koh NJI, Chou LM. 2004a. Flesh or bone? Qunatifying small-scale coral morphology using with-tissue and without-tissue techniques. Mar. Biol. 145: 323-328.
Todd, PA, Laddle RJ, Lewin-Koh NJI, Chou LM. 2004b. Genotype x environmental interactions in transplanted clones of the massive corals Favia speciosa and Diploastrea heliopora. Mar. Ecol. Prog. Ser. 271: 167-182.
Todd PA, Sidle RC, Lewin-Koh NJI. 2004c. An aquarium experiment for identifying the physical factors inducing morphological change in two massive scleractinian corals. J. Exp. Mar. Biol. Ecol. 299: 97-113.
Wijsman-Best M. 1972. Systematics and ecology of New Caledonian Faviinae (Coelenterata – Scleractinia). Contrib. Zool. 42: 3-90.
Wijsman-Best M. 1974. Habitat-induced modification of reef corals (Faviidae) and its consequences for taxonomy. Proc. 2nd Int. Coral Reef Symp., Brisbane 2: 217-228.