Population genetic structure of the endangered limpet Cymbula nigra in a temperate Northern hemisphere region: influence of palaeoclimatic events?

Transcription

1 Marine Ecology. ISSN SHORT COMMUNICATION Population genetic structure of the endangered limpet Cymbula nigra in a temperate Northern hemisphere region: influence of palaeoclimatic events? Free Espinosa 1, Tomoyuki Nakano 2, Jose M. Guerra-García 1 & Jose C. García-Gómez 1 1 Laboratorio de Biologia Marina, Universidad de Sevilla, Seville, Spain 2 Department of Geology and Palaeontology, National Museum of Nature and Science, Tokyo, Japan Keywords Biogeography; Cymbula nigra; glaciations; limpets; Mediterranean. Correspondence Free Espinosa, Laboratorio de Biología Marina, Universidad de Sevilla, Avda. Reina Mercedes 6, 41012, Seville, Spain. free@us.es Accepted: 8 September 2010 doi: /j x Abstract Cymbula nigra (da Costa L, 1771) is the largest limpet in Europe and it is at present considered an endangered species. Nevertheless, its biology is poorly understood and there is a lack of knowledge about population genetic structure. A total of 61 C. nigra specimens were sampled and studied for a fragment of the mitochondrial gene cytochrome c oxidase subunit I. The low genetic diversity, the absence of geographic structure, the significant values reported by measures to test for a recent population expansion and a star-shaped phylogeny pattern are here interpreted as forming a hypothesis concerning the influence of palaeoclimatic events in the genetic pool of limpets, considering also the information available on different limpet species. Introduction The Family Patellidae comprises four extant genera Patella, Cymbula, Helcion and Scutellastra. Most species of Patellidae are found in the Eastern Atlantic, but Scutellastra has a wide distribution ranging from Southern Africa to the Indo-Pacific. In the Mediterranean Sea, only two genera, Cymbula and Patella, are recognized. Cymbula nigra (da Costa, 1771) (see Barea et al. 2008) is the only Cymbula species present in the Mediterranean, although it has a Western African origin (Christiaens 1973). Commonly known as black limpet, C. nigra is mainly distributed in the southernmost coasts of the Iberian Peninsula and North Africa to Angola (Gofas et al. 1988). It was catalogued as an endangered and threatened species (Annex II) by the Barcelona Convention (1993), as strictly protected by the Berna Convention (1995) and as vulnerable by the Andalusian Red List of Threatened Invertebrate Species (Moreno & Arroyo 2008). Little is known about the biology of this species but it seems to be relatively abundant in artificial breakwaters from the Strait of Gibraltar, showing a preference for smooth rocky surfaces (Rivera-Ingraham et al. 2010). Its larval stage and dispersal ability have not been studied in detail. However, Renault & Moueza (1971) reported a reproductive strategy similar to that of Patella caerulea and Patella vulgata (protandrous hermaphroditism), and the duration of the larval period could therefore be similar to these species, at around 10 days (see Dodd 1957). The genetic structure of C. nigra has only been considered in broader studies dealing with limpet phylogenies (see Nakano & Ozawa 2004, 2007; Sá-Pinto et al. 2005, 2008), and little is known about the genetic structure of its populations. The aim of the present study was to explore the genetic relationships of populations of C. nigra in a temperate area of the Northern hemisphere. Material and Methods Sampling sites, number of individuals and codes are shown in Table 1 and Fig. 1. We determined the nucleotide sequence of a 658-bp fragment of the mitochondrial gene cytochrome c oxidase subunit I (hereafter COI) using the primer pair described by Folmer et al. (1994). For details on the experimental protocol see Nakano & Espinosa (2010). Marine Ecology (2010) 1 5 ª 2010 Blackwell Verlag GmbH 1

2 Genetic patterns of Cymbula nigra Espinosa, Nakano, Guerra-García & García-Gómez DNA sequences were aligned using MACCLADE 4.03 (Maddison & Maddison 2002) referring to translated amino acid sequences. Third-codon positions were retained. All haplotypes of 61 individuals were defined using the Redundant Taxa option in MACCLADE 4.03 (Maddison & Maddison 2002), and a minimum spanning tree (MST) (Kruskal 1956; Prim 1957) was constructed Table 1. Sampling localities for Cymbula nigra. Locality Code Station number No. of individuals Benzu (Spain) BE Chorrillo (Spain) CH Cabo Negro (Morocco) CN Cap Malabata (Morocco) CM Asilah (Morocco) A Larache (Morocco) L Rabat (Morocco) R Fatna (Morocco) F Alcazarquivir (Morocco) AL Cap Spartel (Morocco) CS Puerto Saladillo (Spain) PS Estepona (Spain) E TOTAL 61 8 Europe Atlantic Ocean Africa Spain Morocco km 1, 2 3 No. of haplotypes Mediterranean Sea Fig. 1. Map of the sampling localities for Cymbula nigra. Numbers correspond to station number in Table 1. using ARLEQUIN version 2.0 (Schneider et al. 2000). The MST is computed from the matrix of pairwise distances calculated between all pairs of haplotypes using a modification of the algorithm described in Rohlf (1973). The MST is convertible to a minimum spanning network (MSN). Several population parameters were also calculated in ARLEQUIN version 2.0. Haplotype (h) and nucleotide (p) diversity values were estimated according to Nei & Tajima (1981) and Nei (1987). Tajima s D (Tajima 1989) and Fu s F s (Fu 1996) measures were estimated to test for evidence of population expansion and F st fixation index to test for geographical structure. All sequences of Cymbula nigra determined in this study have been deposited in DDBJ and GenBank under accession numbers AB to AB Results and Discussion The MSN obtained for the C. nigra samples shows a starshaped phylogeny with eight different haplotypes (Table 1, Fig. 2). Note that 46 of the 61 individuals analyzed (c. 75% of individuals) had the same haplotype. Values of genetic diversity are low (h = 0.415; p = ), whereas the analyses used to test for evidence of recent population expansion are significantly negative (Tajima s D: -1.84; Fu s F s : -5.82). The fixation index is also very low (F st = ), indicating the absence of geographic structure. [Correction added after online publication 6th January Values of... to...geographic structure inserted]. The minimum spanning network showed a star-shaped phylogeny. According to Rogers (1995) and Avise (2000) this is a signature of species that have recently expanded in size from small or modest numbers of founders, or have experienced recent colonization events (bottleneck and founder effects). Compared with data on other limpet genetic structure (Espinosa & Ozawa 2006 in Patella ferruginea and Patella caerulea; data from Sá-Pinto et al and Mauro et al reported in GenBank either in P. ferruginea and P. caerulea or in P. rustica and P. ulyssiponensis), all patellid species from the temperate Northern Hemisphere show a similar pattern of a star-shaped phylogeny. Furthermore, measures of genetic variability may increase in a spatially structured population (Amos & Harwood 1998). On the other hand, the measures may show quite a uniform structure if the population has expanded rapidly in size following a bottleneck, even if there is some spatial structuring of the population (Walker et al. 1998). The postulated expansion was further supported by the significant negative Tajima s D and Fu s F s values, as well as the absence of geographic structure among populations (F st < 0.05), taking into account that F st values of may indicate little genetic differentiation (Hartl & Clark 2 Marine Ecology (2010) 1 5 ª 2010 Blackwell Verlag GmbH

3 Espinosa, Nakano, Guerra-García & García-Gómez Genetic patterns of Cymbula nigra Hap 8 CS2 Hap 6 R3 Hap 1 (46) Hap 3 CN5 BE1,BE3,BE4,CH1,CH3,CH 4,CH5,CN1,CN2,CN3,CN4,C M2,CM3,CM5,A1,A2,A3,A4, A5,L1,L3,L5,R2,R4,F1,F2,F4,F5,F7,AL1,AL2,AL3,AL4,A L5,AL6,CS1,CS3,PS1,PS2,P S4,PS5,E1,E2,E3,E4,E5 Hap 2 (9) BE2,BE5,CH 2,CM1,L2,L4, F3,F6,PS3 Hap 5 R1 Hap 4 CM4 Hap 7 R5 Fig. 2. Minimum spanning network of the haplotypes found for Cymbula nigra. For designation of individuals see Table 1. Black spots indicate mutations. Numbers in brackets indicate the total number of individuals with the same haplotype. 1989). These results and the presence of shared haplotypes in almost all the sampled populations support the notion of a population expansion of the species related to the colonization of the present distribution area. Several environmental fluctuations could be responsible for the genetic pattern found. On the one hand, during the Messinian salinity crisis ( Ma) the Mediterranean Sea nearly dried up (Bover & Alcover 2000; Bianchi & Morri 2000; Taviani 2003) and many species recolonized the Mediterranean basin through the Gibraltar Strait when the sea level rose again, probably causing the founder effects. On the other hand, fluctuations produced by Quaternary glacial interglacial periods in Europe (see Zazo et al. 2003) may have caused important decreases in population size of patellids and generated bottleneck events. The sea surface temperature (SST) in the Mediterranean Sea ranged from 2.7 to 6.5 C (winter summer) to C during the Quaternary (González-Donoso et al. 2000). For instance, during the last glacial maximum (28 15 kya) the SST was 8 9 C mean of the whole year, around 8 C colder than present values (Melki et al. 2009). These events produced an abrupt modification of the environmental characteristics over relatively short periods and initiated great marine regressions and weak transgressions, changing the shoreline, water air temperatures and habitats, and affecting the settled species. The low levels of genetic variation are more common among intertidal species that live relatively high on the shore, where exposure times to cold stress in air are longer than for species that live at greater depths (Marko 2004). This is the case for the different limpet species in Europe and nearby North African coasts. Furthermore, several papers have reported the relationship between palaeo-environmental changes and genetic patterns in marine species. For instance, the Chinese shrimp Feneropenaeus chinensis may have undergone a genetic bottleneck during the late Pleistocene global glacial period (Li et al. 2009), whereas Earl et al. (2010) have suggested that Pleistocene Holocene marine transgressions produced a star-shaped phylogeny in the tidewater goby Eucyclogobius newberry. As Hewitt (1996) pointed out, several species must have gone through a series of bottlenecks and expansions in population size and range depending on climatic conditions. Lately, several papers have documented star phylogenies related to post-glacial colonization events from southern glacial refugia (see Marko et al and references herein) or have suggested that a recent population expansion has taken place in Atlantic and Mediterranean basins, most probably following the Messinian salinity crisis in the sea urchin Paracentrotus lividus (Calderón et al. 2008). Therefore, the present genetic pattern in Cymbula nigra could have originated from founder events from tropical areas of West Africa, instead of being influenced by glacial periods in Europe. From a conservation perspective, the low number of haplotypes in the populations analyzed could represent an additional threat to this endangered species. It would be desirable to promote interbreeding among specimens from different populations to increase genetic variability, as this may enhance the resilience of the species against environmental changes. Further studies on an extensive number of samples from tropical West Africa would be Marine Ecology (2010) 1 5 ª 2010 Blackwell Verlag GmbH 3

4 Genetic patterns of Cymbula nigra Espinosa, Nakano, Guerra-García & García-Gómez required to determine the genetic structure of the species in areas not influenced by glacial cycles and to test the hypothesis of a greater genetic diversity in these zones. Acknowledgements We express our thanks to Dr Emilio Rolán and Dr José Templado (Museo Nacional de Ciencias Naturales, Spain) for contributing specimens studied in this work. Our gratitude to Autoridad Portuaria de Ceuta and Grant-in- Aid for JSPS Fellows no to T. Nakano from the Japan Society for Promotion of Science for Financial Support. All material has been collected under appropriate collection protocols and approved ethics guidelines. Thanks to Consejeria de Medio Ambiente-Junta de Andalucía for providing facilities to obtain samples from Alboran Island as well as to Obimasa-Consejeria de Medio Ambiente-Ciudad Autónoma de Ceuta. We also are grateful to Dr Gonzalo Giribet and two anonymous referees who have greatly improved the original version of the manuscript. References Amos W., Harwood J. (1998) Factors affecting levels of genetic diversity in natural populations. Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences, 353, Avise J.C. (2000) Phylogeography. The History and Formation of Species. Harvard University Press, Cambridge, MA. Barea J.M., Ballesteros E., Moreno D. (2008) Libro Rojo de los Invertebrados de Andalucía. Consejería de Medio Ambiente, Junta de Andalucía, Seville. Bianchi C.N., Morri C. (2000) Marine biodiversity of the Mediterranean Sea: situation, problems and prospects for future research. Mar. Pollut. Bull., 40 (5), Bover P., Alcover J.A. (2000) Extreme insular evolution in Myotragus balearicus Bate, 1909 (Artiodactyla, Caprinae). Tropics, 10, Calderón I., Giribet G., Turon X. (2008) Two markers and one history: phylogeography of the edible common sea urchin Paracentrotus lividus in the Lusitanian region. Marine Biology, 154, Christiaens J. (1973) Révision du genre Patella (Mollusca, Gastropoda). Bulletin du Musée National d Histoire Naturelle, 182, Dodd J.M. (1957) Artificial fertilisation, larval development and metamorphosis in Patella vulgata L. and Patella caerulea L. Publicazione della Stazione Zoologica di Napoli, 29, Earl D.A., Louis K.D., Bardeleben C., Swift C.C., Jacobs D.K. (2010) Rangewide microsatellite phylogeography of the endangered tidewater goby, Eucyclogobius newberry (Teleostei: Gobiidae), a genetically subdivided coastal fish with limited marine dispersal. Conservation Genetics, 11, Espinosa F., Ozawa T. (2006) Population genetics of the endangered limpet Patella ferruginea (Gastropoda: Patellidae): taxonomic, conservation and evolutionary considerations. Journal of Zoological Systematics and Evolutionary Research, 44, Folmer O., Black M., Hoeh W., Lutz R., Vrijenhoek R. (1994) DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Molecular Marine Biology and Biotechnology, 3, Fu Y.S. (1996) New statistical test of neutrality for DNA samples from a population. Genetics, 143, Gofas S., Pinto Alonso J., Brandao M. (1988) Conchas e moluscos de Angola. Universidade Agostinho Neto, Elf Aquitaine Angola. González-Donoso J.M., Serrano F., Linares D. (2000) Sea surface temperature during the Quaternary at ODP sites 976 and 975 (western Mediterranean). Palaeogeography, Palaeoclimatology, Palaeoecology, 162, Hartl D.L., Clark A.G. (1989) Principles of Population Genetics, 2nd edn. Sinauer Associates, Sunderland, MA. Hewitt G.M. (1996) Some genetic consequences of ice ages, and their role in divergence and speciation. Biological Journal of the Linnean Society, 58, Kruskal J.B. (1956) On the shortest spanning subtree of a graph and the traveling salesman problem. Proceedings of the Mathematical Society, 7, Li Y.L., Kong X.Y., Yu Z.N., Kong J., Ma S., Chen L.M. (2009) Genetic diversity and historical demography of Chinese shrimp Feneropenaeus chinensis in Yellow Sea and Bohai Sea based on mitochondrial DNA analysis. African Journal of Biotechnology, 8, Maddison D.R., Maddison W.P. (2002) MacClade 4: Analysis of Phylogeny and Character Evolution, Version Sinauer Associations, Sunderland, MA. Marko P.B. (2004) What s larvae got to do with it? Disparate patterns of post-glacial population structure in two benthic marine gastropods with identical dispersal potential. Molecular Ecology, 13, Marko P.B., Hoffman J.M., Emme S.A., McGovern T.M., Keever C.C., Cox L.N. (2010) The Expansion-Contraction model of Pleistocene biogeography: rocky shores suffer a sea change? Molecular Ecology, 19, Mauro A., Arculeo M., Parrinello N. (2003) Morphological and molecular tools in identifying the Mediterranean limpets Patella caerulea, Patella aspera and Patella rustica. Journal of Experimental Marine Biology and Ecology, 295, Melki T., Kallel N., Jorissen F.J., Guichard F., Dennielou B., Berne S., Labeyrie L., Fontugne M. (2009) Abrupt climate change, sea surface salinity and paleoproductivity in the western Mediterranean Sea (Gulf of Lion) during the last 28 kyr. Palaeogeography, Palaeoclimatology, Palaeoecology, 279, Marine Ecology (2010) 1 5 ª 2010 Blackwell Verlag GmbH

5 Espinosa, Nakano, Guerra-García & García-Gómez Genetic patterns of Cymbula nigra Moreno D., Arroyo M.C. (2008) Cymbula nigra (da Costa, 1771). In: Barea-Azcón J.M., Ballesteros E., Moreno D. (Eds), Libro Rojo de los Invertebrados de Andalucía. Vol. 1. Consejería de Medio Ambiente, Junta de Andalucía, Seville: Nakano T., Espinosa F. (2010) New alien species in the Atlantic Ocean? Marine Biodiversity Records. doi: / S Nakano T., Ozawa T. (2004) Phylogeny and historical biogeography of limpets of the order Patellogastropoda based on mitochondrial DNA sequences. Journal of Molluscan Studies, 70, Nakano T., Ozawa T. (2007) Worldwide phylogeography of limpets of the Order Patellogastropoda: molecular, morphological and paleontological evidence. Journal of Molluscan Studies, 73, Nei M. (1987) Molecular Evolutionary Genetics. Columbia University Press, New York. Nei M., Tajima F. (1981) DNA polymorphism detectable by restriction endonucleases. Genetics, 97, Prim R.C. (1957) Shortest connection networks and some generalizations. Bell System Technical Journal, 36, Renault L., Moueza M. (1971) Contribution a l etude de Patella safiana Lamarck. Haliotis, 1, Rivera-Ingraham G., Espinosa F., García-Gómez J.C. (2010) Niche differentiation between juvenile and adult black limpets (Cymbula nigra). Journal of the Marine Biological Association of the United Kingdom. doi: / S Rogers A.R. (1995) Genetic evidence for a Pleistocene population explosion. Evolution, 49, Rohlf F.J. (1973) Algorithm 76. Hierarchical clustering using the minimum spanning tree. Computer Journal, 16, Sá-Pinto A., Branco M., Harris D.J., Alexandrino P. (2005) Phylogeny and phylogeography of the genus Patella based on mitochondrial DNA sequence data. Journal of Experimental Marine Biology and Ecology, 325, Sá-Pinto A., Branco M., Sayanda D., Alexandrino P. (2008) Patterns of colonization, evolution and gene flow in the species of the genus Patella in the Macaronesian Islands. Molecular Ecology, 17, Schneider S., Roessli D., Excoffier L. (2000) ARLEQUIN ver 2.000: A Software for Population Genetics Data Analysis. University of Geneva, Geneva. Tajima F. (1989) Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics, 123, Taviani M. (2003) Shaping the biogeography of the Mediterranean basin: one geologist s perspective. Biogeographia, 24, Walker D.E., Moler P.E., Buhlmann K.A., Avise J.C. (1998) Phylogeographic uniformity in mitochondrial DNA of the snapping turtle (Chelydra serpentina). Animal Conservation, 1, Zazo C., Goy J.L., Dabrio C.J., Bardají T., Marcel C.H., Ghaleb B., González-Delgado J.A., Soler V. (2003) Pleistocene raised marine terraces of the Spanish Mediterranean and Atlantic coasts: records of coastal uplift, sea-level highstands and climate changes. Marine Geology, 194, Marine Ecology (2010) 1 5 ª 2010 Blackwell Verlag GmbH 5