Population connectivity buffers genetic diversity loss in a seabird
Frontiers in Zoologyvolume 10, Article number: 28 (2013)
Ancient DNA has revolutionized conservation genetic studies as it allows monitoring of the genetic variability of species through time and predicting the impact of ecosystems’ threats on future population dynamics and viability. Meanwhile, the consequences of anthropogenic activities and climate change to island faunas, particularly seabirds, remain largely unknown. In this study, we examined temporal changes in the genetic diversity of a threatened seabird, the Cory’s shearwater (Calonectris borealis).
We analysed the mitochondrial DNA control region of ancient bone samples from the late-Holocene retrieved from the Canary archipelago (NE Atlantic) together with modern DNA sequences representative of the entire breeding range of the species. Our results show high levels of ancient genetic diversity in the Canaries comparable to that of the extant population. The temporal haplotype network further revealed rare but recurrent long-distance dispersal between ocean basins. The Bayesian demographic analyses reveal both regional and local population size expansion events, and this is in spite of the demographic decline experienced by the species over the last millennia.
Our findings suggest that population connectivity of the species has acted as a buffer of genetic losses and illustrate the use of ancient DNA to uncover such cryptic genetic events.
Over the last millennia island faunas have been progressively decimated by human exploitation, biological invasions and habitat destruction . Deterministic threats of that kind have direct mortality effects on island populations leading to short-term demographic declines. This would in turn lead to a loss of genetic diversity and an increase of genetic stochasticity in small populations, such as inbreeding and population bottlenecks, ultimately increasing the probability of extinction . This global process is particularly well exemplified among seabirds in general and albatrosses and petrels in particular, hereafter Procellariiforms . During the Holocene at least 17 Procellariiform species have vanished , and 56% of the extant species are known or suspected to be experiencing population declines . In the Canary Islands this process has led to the extinction of the dune shearwater (Puffinus holeae), the lava shearwater (P. olsoni), and an unidentified species of petrel (Pterodroma sp.) [6–8].
One of the most abundant Procellariform species inhabiting the Canary archipelago is the Atlantic Cory´s shearwater (Calonectris borealis). The breeding population of the species has been estimated as more than 270,000 pairs, however it underwent an acute decline between 1970–1990 and currently it continues to decline . This species faces a number of conservation threats globally, both at sea (i.e. fisheries by-catch) and at land (i.e. animal invasions), which have been shown to seriously impact adult survival and reproductive success (see ). In the Canary Islands the Cory’s shearwaters have been by and large exploited as a food resource by the aboriginal peoples who reached the islands sometime between 756 cal BC - 313 cal AD and has also suffered predation by introduced rats and cats for centuries . This has possibly led to a population bottleneck. Nevertheless, the impact of these threats on the genetic status of the species is unknown.
Ancient DNA (aDNA) has revolutionised biodiversity and conservation studies because it provides an excellent opportunity to monitor the genetic variability of species or populations through time . However, few studies to date have used such heterochronous' data sets for conservation purposes in seabirds [13, 14]. This is partly due to the inherent methodological difficulties of aDNA research , but also to the scarcity of fossil remains in remote oceanic archipelagos.
Using an heterochronous sampling consisting of Cory’s shearwater bone samples from late-Holocene from the Canary Islands together with current sequences representative of the entire breeding range of two sister Calonectris species, the Cory’s and the Scopoli’s (C. diomedea) shearwaters, we aim: (i) to estimate levels of genetic diversity and relationships between ancient and modern Cory’s shearwater populations, and (ii) to reconstruct the population demographic history of the Cory’s shearwater and estimate population size changes within the Canary archipelago through time.
All collected bones were morphologically identified as belonging to the Cory’s shearwater species. 14C dates show that the three samples analyzed are from the last millennium. The 2σ confidence interval of two samples fall between 1004–1530 cal yr AD being the third younger than 1499 cal AD (Additional file 1: Table S1). DNA from pooled ancient samples (N=14) resulted in 25 mtDNA distinct haplotypes (Additional file 2: Figure S1). Only three subsamples contained one single haplotype, and the remaining subsamples presented two haplotypes each. The authenticity of these 25 haplotypes is supported by: 1) negative controls at the extraction and PCR stages were uniformly clean, 2) all the fragments have been replicated, by independent PCR, at least twice and on average 3.5 times (Additional file 2: Figure S1), 3) ancient sequences obtained clearly cluster within the extant Cory’s shearwater gene pool (Figure 1 and 4) all haplotypes were unique to single subsamples, indicating that there had been no cross contamination between samples/sites.
The genetic diversity values were high and similar between the ancient and the modern population (Table 1). Most individuals presented unique haplotypes and only five of the 25 haplotypes detected in the pooled ancient samples correspond to haplotypes still present in current populations, two of them being nested within the most common modern haplotypes (Figure 1). The remaining 20 ancient haplotypes were not present in any of the 282 extant Cory’s shearwaters analysed, including the 42 newly analysed individuals from Montaña Clara [16, 17] , see Additional file 3: Table S2.
The temporal haplotype network obtained illustrates the high genetic diversity found in both ancient (Atlantic) and modern samples, and agrees with the marked genetic differentiation between the two Calonectris shearwaters sister clades, so all ancient samples clustered with Atlantic Cory’s shearwater extant haplotypes. However, there is one ancient haplotype that grouped within the Mediterranean Scopoli’s cluster (Figure 1). These findings are concordant with those found in a previous study on modern samples , that show various introgression and migration events as well as low but significant levels of interspecific gene flow [16, 17]. Our results on the ancient samples add further evidence of the rare albeit recurrent hybridization events between the two species.
The demographic analysis on the Cory’s shearwater indicated significant departures from a constant population model with the majority of credible sets containing non-zero population size changes. The Extended Bayesian Skyline Plot (EBSP) reconstruction of changes in genetic diversity through time shows a population expansion event at about 50,000 yr (Figure 2A). At a local scale and using radio-carbon rates to calibrate the demographic time-function (see Additional file 4: Supplementary methods), the EBSP analysis suggests a slight increase in genetic diversity in the Canary Archipelago over the last five millennia (Figure 2B).
Overall, these findings would not support the initial hypothesis of a loss in genetic variability in the Cory’s shearwater, in spite of its marked demographic population decline over the last millennium in the Canaries . Several factors such as the high dispersal ability, the high effective population size, and the long generation time of the species may have contributed to offset the loss caused by these anthropogenic impacts . The finding that high connectivity can buffer genetic diversity in the face of a demographic decline, has been reported in marine mammals (i.e. ), which have in common with seabirds the ability to travel long distances. This hypothesis is further supported by previous genetic and ecological data on the species that show high levels of population connectivity at both regional and local scales . While the potential for dispersal will likely help the persistence of the colony in the short-term, future viability is debatable considering the life-history traits of this species. That is, Calonectris shearwaters share the typical and extreme life history traits of Procellariiform seabirds: long-lived, low fecundity and high survival, and thereby the species is particularly susceptible to threats increasing adult mortality . In this regard, buffering effects on the genetic diversity may not be enough on a long-term basis if threats affecting the declining of populations persist. More generally, our study illustrates how aDNA can assist to better understand how populations function and evolve across time and so predict the impact of future anthropogenic and environmental changes on the seabird community, which is essential for implementing effective management strategies.
All experimental procedures on ancient samples were performed in a dedicated aDNA laboratory (IBE-PRBB, Barcelona), where no previous work with extant Calonectris shearwaters had been conducted. Ancient bone fragments were collected in Montaña Clara and Alegranza islets in the northeastern part of the Canaries (Figure 3). All samples were collected from 10 to 50 cm from the surface and separated by sampling location. Three of the best-preserved bones from Montaña Clara were dated by accelerator mass spectrometer radiocarbon analysis (AMS 14C). The rest of the bones were subsequently used for genetic analysis. Due to the fragmentary nature of the remains, the bones of different specimens could not be individualized. Subsequently, we generated fourteen subsamples (one for each of the fourteen sampling places) of pooled skeletal remains of similar weight (around 150 mg) aiming to identify a posteriori the different mitochondrial (mtDNA) haplotypes present on each subsample analysed. This pooling strategy was imposed by the limitation in size of the remains, and although it does not allow determination of the precise number of individuals analysed (only a minimum of 25, see results), it is still a useful inexpensive way to estimate the genetic diversity within a population . The fourteen pooled bones were then ground to powder and used for DNA extraction as described elsewhere . We amplified a fragment of 218 bp of the mtDNA control region (see Additional file 5: Table S3). GenBank accessions numbers KC888879– KC888904.
The modern samples included 240 mtDNA control region sequences published in previous studies from across the entire distribution range of the two sister species, the Cory’s and the Scopoli’s shearwaters [16, 17]. A subset of 42 modern samples from Montaña Clara (Lanzarote, Canary Is.) were further amplified and included in the analyses. For these samples, DNA extraction and amplification methods were performed at IBE-CMIMA laboratory (Barcelona) as described elsewhere . See Additional file 3: Table S2 for further details on the dataset. GenBank accessions numbers KC888905–KC888924.
Genetic statistics for both modern and ancient samples were performed in the program DnaSP, v. 5.10.01 . To depict the genealogical relationships between ancient and modern DNA samples and the continuity of haplotypes across time, a temporal network was constructed using TempNet . For this analysis and to facilitate data interpretation, only unique haplotypes were considered.
To investigate the demographic history of Cory’s shearwaters from the Canary archipelago, we estimated the shape of the population growth function through time by constructing Extended Bayesian Skyline Plots (EBSP) as implemented in BEAST v1.7.4  (Additional file 1: Table S1; Additional file 2: Figure S1; Additional file 3: Table S2; Additional file 4: Supplementary methods; Additional file 5: Table S3).
Illera JC, Rando JC, Richardson DS, Emerson BC: Age, origin and extinctions of the avifauna of Macaronesia: a synthesis of phylogenetic and fossil information. Quat Sci Rev. 2012, 50: 14-22.
Frankham R, Ballou JD, Briscoe DA: Introduction to conservation genetics. 2002, Cambridge: Cambridge University Press
Croxall JP, Butchart SHM, Lascelles B, Stattersfield AJ, Sullivan B, Symes A, Taylor P: Seabird conservation status, threats and priority actions: a global assessment. Bird Conserv Int. 2012, 22: 1-34. 10.1017/S0959270912000020.
Tyrberg T: Holocene avian extinctions. Holocene Extinctions. Edited by: Turvey ST. 2009, Oxford: Oxford University Press, 63-106. 352 pp
Scofield RP: Procellariform extinctions in the Holocene: threat processes and wider ecosystem-scale implications. Holocene Extinctions. Edited by: Turvey ST. 2009, Oxford: Oxford University Press, 151-166. 352 pp
Rando JC: New data of fossil birds from El Hierro (Canary Islands): Probable causes of extinction and some biogeographical considerations. Ardeola. 2002, 49: 39-49.
Rando JC, Alcover JA: Evidence for a second western Palaearctic seabird extinction during the last Millennium: the Lava Shearwater Puffinus olsoni. Ibis. 2008, 150: 188-192.
Rando JC, Alcover JA: On the extinction of Dune Shearwater (Puffinus holeae) from the Canary Islands. J Ornithol. 2010, 151: 365-369. 10.1007/s10336-009-0463-6.
BirdLife International: Birds in Europe: population estimates, trends and conservation status. 2004, Cambridge, UK: BirdLife International. (BirdLife Conservation Series No.12)
Thibault JC, Bretagnolle V, Rabouam C: Cory’s shearwater. 1997, Oxford: Oxford University Press
Rando JC, López M, Jimenez MC: Bird remains from the archaeological site of Guinea (El Hierro, Canary Islands). Inter J Osteoarch. 1997, 7: 298-302. 10.1002/(SICI)1099-1212(199707/08)7:4<298::AID-OA384>3.0.CO;2-L.
Leonard JA: Ancient DNA applications for wild life conservation. Mol Ecol. 2008, 17: 4186-4196. 10.1111/j.1365-294X.2008.03891.x.
Heupink TH, van den Hoff J, Lambert DM: King penguin population on Macquarie Island recovers ancient DNA diversity after heavy exploitation in historic times. Biol Lett. 2012, 23: 586-589.
Welch AJ, Wiley AE, James HF, Ostrom PH, Stafford TW, Fleischer RC: Ancient DNA reveals genetic stability despite demographic decline: 3,000 years of population history in the endemic Hawaiian petrel. Mol Biol Evol. 2012, 29: 3729-3740. 10.1093/molbev/mss185.
Hofreiter M, Serre D, Poinar HN, Kuch M, Pääbo S: Ancient DNA. Nat Rev Genet. 2001, 2: 353-359.
Gómez-Díaz E, González-Solís J, Peinado MA: Population structure in a highly pelagic seabird, the Cory’s shearwater Calonectris diomedea: an examination of genetics, morphology and ecology. Mar Ecol Prog Ser. 2009, 382: 197-209.
Gómez-Díaz E, González-Solís J, Peinado MA, Page RDM: Phylogeographyof the Calonectris shearwaters using molecular and morphometric data. Mol Phylogenet Evol. 2006, 41: 322-332. 10.1016/j.ympev.2006.05.006.
Friesen VL, Burg TM, McCoy KD: Mechanisms of population differentiation in seabirds. Mol Ecol. 2007, 16: 1765-1785. 10.1111/j.1365-294X.2006.03197.x.
Foote AD, Kaschner K, Schultze S, Gariloa C, Ho SYW, Higham TFG, Van der Es H, Embling C, Gregersen K, Johansson F, Willerslev E, Gilbert MTP: Ancient DNA reveals that bowhead whale lineages survived Late Pleistocene climate change and habitat shifts. Nat Commun. in press
Biology of marine birds. Edited by: Schreiber EA, Burger J. 2002, Boca Raton, FL: CRC Press
Furman BJ, Grattapaglia D, Dvorak WS, O’Malley DM: Analysis of genetic relationships of Central American and Mexican pines using RAPD markers that distinguish species. Mol Ecol. 1997, 6: 321-331. 10.1046/j.1365-294X.1997.00183.x.
Ramírez O, Illera JC, Rando JC, González-Solis J, Alcover JA, Lalueza-Fox C: Ancient DNA of the extinct lava shearwater (Puffinus olsoni) from the Canary Islands reveals incipient differentiation within the P. puffinus complex. PLoS One. 2010, 5 (12): e16072-10.1371/journal.pone.0016072.
Librado P, Rozas J: DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25: 1451-1452. 10.1093/bioinformatics/btp187.
Prost S, Anderson CNK: TempNet: a method to display statistical parsimony networks for heterochronous DNA sequence data. Methods Ecol Evol. 2011, 2: 663-667. 10.1111/j.2041-210X.2011.00129.x.
Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.
O.R. is a postdoctoral Researcher from the JAEDOC program cofounded by ESF. E.G.-D. and J.C.I were supported by a JdC and a RyC from the MICINN, Spain. This work has been founded by grants BFU2012-34157 to C.L-F and CGL2009-11278/BOS o J.G-S from the MINECO, Spain, and by a MC reintegration grant to E.G.-D. (ERG-2010-276838). We thank the Cabildo and Gobierno de Canarias and owners of Montaña Clara for permissions, Pascual Calabuig and Loly Estévez for logistic support and Teresa Darbyshire for checking the English. We thank for comments provided by Andy Foote and one anonymous referee on a previous versions of the manuscript.
The authors declare that they have no competing interests.
OR, EG-D, JG-S, JCI, JCR and CL-F contributed to the design of this research. OR, EG-D, and IO performed the experimental analyses. OR, EG-S and JCI performed the analysis. OR, EG-D, JCR, JCI, JG-S and CL-F wrote the manuscript. All authors read and approved the final manuscript.
Oscar Ramírez, Elena Gómez-Díaz contributed equally to this work.