Phylogenetic relationships of the extinct St Helena petrel, Pterodroma rupinarum Olson, 1975 (Procellariiformes: Procellariidae), based on ancient DNA ANDREANNA J. WELCH1,*, STORRS L. OLSON2 and ROBERT C. FLEISCHER3 1Department of Biological Sciences, University at Buffalo, Buffalo, NY 14260, USA 2Department of Vertebrate Zoology, National Museum of Natural History, Smithsonian Institution, Washington, DC 20013, USA 3Smithsonian Conservation Biology Institute, Center for Conservation and Evolutionary Genetics, National Zoological Park, Washington, DC 20013, USA Received 9 May 2013; revised 29 July 2013; accepted for publication 2 August 2013 Palaeontological studies show that three endemic procellariid seabird species became extinct on the remote island of St Helena in the South Atlantic Ocean. At least one of these, Pterodroma rupinarum Olson, 1975, is likely to have survived until human colonization of the island, although it is known only from subfossil bones. Several species of Pterodroma are distributed across the Atlantic, but the skull and bill of Pt. rupinarum were judged to be more similar to the Indo-Pacific Pterodroma rostrata group, which was recently split into the separate genus Pseudobulweria. We used ancient DNA techniques to sequence the mitochondrial cytochrome b gene of the extinct Pt. rupinarum, and conducted phylogenetic analyses to investigate the placement of this enigmatic taxon. In trees constructed using maximum likelihood and Bayesian inference, Pt. rupinarum did not group with Pseudobulweria, but instead fell within a strongly supported clade of Atlantic Pterodroma, including the endangered Black-capped [Pterodroma hasitata (Kuhl, 1820)] and Bermuda [Pterodroma cahow (Nichols & Mowbray, 1916)] petrels, as well as the Macaronesian petrels [Pterodroma madeira Mathews, 1934, Pterodroma feae (Salvadori, 1899) and Pterodroma deserta Mathews, 1934]. Pterodroma rupinarum shared a particularly close relationship with Pt. feae of the Cape Verde Islands, which is also the geographically closest species within the clade. Considering the osteological distinctiveness of Pt. rupinarum it was probably a separate species, or at least a highly diverged population that was isolated for a substantial period of time prior to its extinction. © 2013 The Linnean Society of London, Zoological Journal of the Linnean Society, 2013 doi: 10.1111/zoj.12078 ADDITIONAL KEYWORDS: cytochrome b – extinction – petrel – phylogeny – Pseudobulweria – Pterodroma – St Helena. INTRODUCTION Seabirds are highly vagile, and may forage widely over the open oceans, travelling hundreds to thousands of kilometers on a single foraging trip during the breeding season (Brooke, 2004). Their ability to disperse widely over expansive pelagic eco- systems makes them a useful indicator of marine conditions (Piatt, Sydeman & Wiese, 2007; Wiley et al., 2013), and as top marine predators they play ecologically important roles (Sánchez-Piñero & Polis, 2000; Croll et al., 2005). Unfortunately, seabirds are one of the most endangered groups of birds in the *Corresponding author. Current address: School of Biological and Biomedical Sciences, Durham University, UK; Institute of Marine Sciences, University of California, Santa Cruz, CA, USA. E-mail: andreanna05@gmail.com bs_bs_banner Zoological Journal of the Linnean Society, 2013. With 2 figures © 2013 The Linnean Society of London, Zoological Journal of the Linnean Society, 2013 1 world (Croxall et al., 2012). The family Procellariidae contains approximately 80 species of petrels and shearwaters, although the taxonomy of this group is challenging and often disputed (Brooke, 2004; Penhallurick & Wink, 2004; Rheindt & Austin, 2005). The Procellariidae ranks as the third most endan- gered seabird family (Croxall et al., 2012), and many of these species face substantial threats both at the breeding grounds as well as at sea (Weimerskirch, Brothers & Jouventin, 1997; Lewison & Crowder, 2003; Wanless et al., 2007; Fontaine, Gimenez & Bried, 2011; Barbraud et al., 2012). Despite this, investigation of the fossil record indicates that in the recent past seabird extinction has been relatively rare, especially when compared with the extinction rate of terrestrial species on oceanic islands (Steadman, 1995; Steadman, 2006). Fossils from Saint Helena, a remote island in the South Atlantic Ocean (Fig. 1), showed that three apparently endemic species of Procellariidae became extinct there, and populations of three other pro- cellariids, along with several additional species of seabirds, were extirpated (Olson, 1975). Of the three extinct species, only one, Puffinus pacificoides Olson, 1975, is likely to have disappeared prior to human arrival, probably during the Pleistocene. Another species, Bulweria bifax Olson, 1975, was much more abundant in the Pleistocene than in recent deposits, and may have been decreasing in number towards modern times. The third species, Pterodroma rupinarum Olson, 1975, was fairly common in the most recent deposits, and may have become extinct after the discovery of the island in 1502 by Portu- guese sailors (Olson, 1975). Little is known about the P. rupinarum P. arminjoniana P. cahow P. hasitata P. caribbaea P. feae P. deserta P. madeira P. mollis P. incerta † † Figure 1. Approximate recent breeding locations of Pterodroma petrels in the Atlantic Ocean. 2 A. J. WELCH ET AL. © 2013 The Linnean Society of London, Zoological Journal of the Linnean Society, 2013 avifauna of St Helena at the time of human coloniza- tion, as there appear to be no written accounts from that period; however, after permanent settlement in 1513 there are numerous accounts of large numbers of feral invasive species, including goats, pigs, rats, mice, and cats (Gosse, 1938), which would presum- ably have caused severe environmental perturbation and led to greatly increased predation on burrow- nesting seabirds. Initial work by Bourne (1956) on the taxonomic identification of a limited number of Pterodroma bones from St Helena ruled out affinity with Pterodroma arminjoniana (Giglioli & Salvadori, 1869), the only other tropical South Atlantic petrel (Fig. 1). Instead, he suggested that the fossils might belong to Macaronesian species of Pterodroma breeding on Madeira and the Cape Verde Islands in the North Atlantic. Recent genetic work on the Macaronesian petrels suggests the existence of three recently diverged species: Pterodroma feae (Salvadori, 1899) in the Cape Verde archipelago; Pterodroma deserta Mathews, 1934 on Bugio Island (first sug- gested as a distinct species by Bourne in 1955); and Pterodroma madeira Mathews, 1934 on the island of Madeira in the Madeira Archipelago (Fig. 1), although the species status of Pt. deserta is still somewhat controversial (Zino, Brown & Biscoito, 2008; Jesus et al., 2009; Gangloff et al., 2013). Bourne (1965) later suggested that the fossil Pterodroma from St Helena might belong with the Pterodroma rostrata group instead. Examination of greatly expanded material collected by Olson (1975) resulted in the naming of the fossils as a new species, Pt. rupinarum, some characters of which, particularly the skull and the angle of the bill, were considered to be more similar to the Pterodroma rostrata group than to other Pterodroma petrels. Genetic work (Bretagnolle, Attié & Pasquet, 1998) has since shown that the Pterodroma rostrata group actually repre- sents a distinct genus, Pseudobulweria, first proposed by Mathews in 1936 [type species Thalassidroma (Bulweria) macgillivrayi G.R. Gray, 1859]. The Indo-Pacific genus Pseudobulweria contains four extant species: Pseudobulweria rostrata (Peale, 1848), Pseudobulweria becki (Murphy 1928), Pseudobulweria aterrima (Bonaparte, 1857), and Pseudobulweria macgillivrayi (Gray, 1860). These are among the most endangered of procellariiforms, with three out of four of the species considered critically endangered, being known only from a handful of museum specimens, and often having been feared extinct for many years (Attié, Stahl & Bretagnolle, 1997; Shirihai, 2008; Shirihai et al., 2009; Gangloff et al., 2012). Here, we report on ancient DNA sequences from the extinct Pt. rupinarum and assess the placement of this species in a phylogeny including nearly all Pterodroma and Pseudobulweria petrels. MATERIAL AND METHODS SAMPLES AND MOLECULAR TECHNIQUES Of the thousands of available bones of Pt. rupinarum collected by Olson in 1971, ten were selected from three sites on St Helena: eight from Prosperous Bay (USNM 176607, 176608, 176610, 176613, 176615, 176623, 176627, 176641), one from Sandy Bay (USNM 176679), and one from Sugarloaf Hill, site 2 (USNM 176702). Bones that were clearly identifiable as Pt. rupinarum were selected based on size, thick- ness of bone walls, and un-weathered appearance. A mixture of complete and partial elements was selected to minimize the effects of destructive sam- pling. The age of the samples is not certainly known, but all of the deposits are thought to be late Holocene. A small piece of bone (approximately 4 mm × 6 mm) was cut from each specimen using a rotary saw with disposable blades. Each sample was ground into powder using a mortar and pestle and DNA was extracted using a phenol/chloroform procedure with centrifugal dialysis (Fleischer et al., 2000). All labora- tory work was conducted in a physically isolated facility dedicated for ancient DNA work, and strict protocols were followed to prevent contamination from modern sources. Ancient DNA was extracted in batches of six samples, each of which included an extraction control. Most phylogenetic work on procellariiform seabirds has focused on the mitochondrial cytochrome b gene (Bretagnolle et al., 1998; Nunn & Stanley, 1998; Austin, Bretagnolle & Pasquet, 2004; Gangloff et al., 2012), so we targeted that gene for our study. Polymerase chain reaction (PCR) was carried out as described in Welch et al. (2012). Briefly, PCR was conducted in 25-μL reaction volumes with 1 unit Amplitaq Gold DNA polymerase (Applied Biosystems), 2–3 μL of DNA extract, and with a total of 45 cycles. We used a series of short (120–200 bp), overlapping PCR primers designed specifically to amplify in Pterodroma (Welch et al., 2012), in Pseudobulweria and Bulweria (Table 1), and in all Procellariformes (Pyle, Welch & Fleischer, 2011). Negative controls were used to detect potential contamination. All amplifications were conducted a minimum of two times for each primer–sample combination. PCR products were cleaned using ExoSAP-IT (USB), cycle-sequenced in both directions with Big Dye Terminator 3.1 (Applied Biosystems), then purified using Sephadex G-50 fine columns (GE Healthcare Bio-Sciences), and electrophoresed in an ABI 3130 xL Genetic Analyzer (Applied Biosystems). ST HELENA PETREL, PTERODROMA RUPINARUM 3 © 2013 The Linnean Society of London, Zoological Journal of the Linnean Society, 2013 Sequences were visually inspected and edited in SEQUENCHER 4.9. DATA ANALYSIS In addition to the data we obtained for Pt. rupinarum, we downloaded representative DNA sequences from the GenBank database for procellariid taxa (Table 2), selecting sequences from all species of Pseudobulweria and nearly all species of Pterodroma (see Results). Here, we consider the three Macaronesian taxa of Pterodroma from the North Atlantic as distinct species. Only partial sequences were available for the enigmatic Ps. aterrima, Ps. becki, and Ps. macgillivrayi, so for those species multiple sequences were included. We also included multi- ple sequences for taxa of interest that breed in several different locations (e.g. Ps. rostrata and Pt. arminjoniana). Diving petrels of the genus Pelecanoides were used to root the tree (Nunn & Stanley, 1998). The cytochrome b gene is relatively well conserved in Procellariiformes, and sequences were aligned in SEQUENCHER and then manually adjusted by eye. To detect the presence of nuclear copies of mitochondrial DNA in our data set (NUMTs; Sorenson & Fleischer, 1996), the sequences were char- acterized and translated in MacClade 4.08 (Maddison & Maddison, 2008). The program jModelTest 0.1.1 (Posada, 2008) and the Akaike information criterion (AIC) were used to select the best-fitting substitution model. The GTR + I + G model was selected for the full data set as well as for data sets composed of the nucleotides in first and third codon positions, and the GTR + I model was selected for the data set composed of nucleotides in the second codon position. The transition/transversion bias was estimated using maximum likelihood in MEGA 5 (Nei & Kumar, 2000; Tamura et al., 2011). Calculations of pairwise sequence divergence using the best-fit substitution model were conducted in PAUP* (Swofford, 2002) and MEGA 5. We constructed phylogenetic trees using maximum likelihood and Bayesian inference. For both methods we first constructed trees using the full data set and then performed the analysis again, partitioning the data by codon position. Maximum-likelihood analyses were conducted with the RAxML BlackBox webser- ver (Stamatakis, Hoover & Rougemont, 2008), and branch support was estimated through 1000 boot- strap replicates. Bayesian inference was conducted using MrBayes 3.1.2 (Ronquist & Huelsenbeck, 2003; Altekar et al., 2004). For the codon-partitioned analy- ses, some parameters of the GTR substitution matrix in the partition for the first codon position had par- ticularly low effective sample size (ESS) compared with the others (e.g. some types of transversions were rare in the data set), leading to poor mixing overall; therefore, the HKY + I + G model was implemented instead. Default priors were used and all parameters except for branch lengths and topology were unlinked. We used the default heating scheme with one cold and three heated chains. Runs were conducted for 10 million generations, with the first 10% discarded as burn-in. We conducted multiple independent analyses with different, randomly generated starting seeds. We examined plots of −lnL across generations and assessed the ESS for each parameter (> 500 in all cases) in TRACER 1.5 (Drummond et al., 2012). To test whether the codon-partitioned model fit the data better than the unpartitioned model, we conducted a likelihood ratio test (LRT) for maximum-likelihood analyses, and estimated Bayes factors (BF) in TRACER for Bayesian analyses (Suchard, Weiss & Sinsheimer, 2001). We investigated the divergence time of Pt. rupinarum using a strict molecular clock gene tree analysis in BEAST 1.7.4 (Drummond et al., 2012). The analysis was performed using the unpartitioned data set with the GTR + I + G substitution model, empirical base frequencies, and a random starting tree. We used the Yule tree prior with other priors set at their default values. The substitution rate for the cytochrome b gene was set to 1.89% per million years, as estimated for procellariiform seabirds using fossil and biogeographic calibrations (Weir & Schluter, 2008). Multiple independent analyses with randomly generated starting seeds were conducted for 50 million generations, with the first 10% discarded as burn-in. As above, plots of −lnL across generations and ESS were examined in TRACER 1.5. Table 1. Polymerase chain reaction (PCR) primers devel- oped from Pseudobulweria and Bulweria sequences that were used in this study, in addition to primers developed from sequences of other procellariiform taxa (Pyle et al., 2011; Welch et al., 2012) Primer Sequence (5′→3′) aJugAF GGATCTCTCCTAGGYATCTGCCTAAT aJugAR TCGGATTAGTCAGCCGTAYTG aJugBF CAACCCTAGCCTTTTCATCC aJugBR CCGTAATAGAATCCTCGTCCA aJugCF CAAACGGAGCCTCATTYTTC aJugCR ATCCTACGAAGGCRGTTGCT aJugDF CCTGAGGGGGATTCTCAGTA aJugDR ACCTAGGGGGTTGTTTGAGC aJugEF CGCAGGACTCASMCTAATCCA aJugER TGGGAGGAATATAATTGTGAARCC aJugFF TCCACCCCTACTTCACCCTA aJugFR GGTTTGATATGRGGRGGTGT aJugGF AACTTCACCCCAGCAAACC aJugGR GGAGGAATGGGGTTAGGAAT 4 A. J. WELCH ET AL. © 2013 The Linnean Society of London, Zoological Journal of the Linnean Society, 2013 Table 2. Taxa included in phylogenetic analyses and the corresponding accession number for the mitochondrial cytochrome b sequence obtained from GenBank Taxon GenBank accession no. Bulweria bulwerii U74351 Calonectris diomedea AY139626 Calonectris leucomelas AF076045 Daption capense AF076046 Fulmarus glacialis U74348 Fulmarus glacialoides AF076055 Halobaena caerulea AF076057 Aphrodroma brevirostris NC_007174 Macronectes giganteus U48941 Macronectes halli AF076061 Pachyptila desolata AF076068 Pachyptila salvini AF076069 Pachyptila turtur AF076070 Pachyptila vittata U74349 Pagodroma nivea AF076071 Pelecanoides garnoti AF076073 Pelecanoides georgicus AF076074 Pelecanoides magellanicus AF076075 Pelecanoides urinatrix AF076076 Procellaria aequinoctialis U74350 Procellaria cinerea U48940 Procellaria parkinsoni AF076077 Procellaria westlandica AF076078 Pseudobulweria aterrima JF522101 Pseudobulweria aterrima U70483 Pseudobulweria becki JF522102 Pseudobulweria macgillivrayi JF522104 Pseudobulweria macgillivrayi JF522105 Pseudobulweria rostrata rostrata – Fiji (F) JF522103 Pseudobulweria rostrata rostrata – Marquesas (M) JF522106 Pseudobulweria rostrata rostrata – Gambier (G) U70482 Pseudobulweria rostrata trouessarti – New Caledonia (N) U70493 Pterodroma arminjoniana – Round (R) GQ328970 Pterodroma arminjoniana – Trindade (T) GQ328969 Pterodroma axillaris U74342 Pterodroma baraui U70489 Pterodroma cahow U74331 Pterodroma cookii U74345 Pterodroma deserta U74333 Pterodroma externa U74339 Pterodroma feae FJ196356 Pterodroma hasitata EU167017 Pterodroma heraldica GQ328988 Pterodroma hypoleuca AF076079 Pterodroma incerta U74335 Pterodroma inexpectata U74346 Pterodroma lessonii U74337 Pterodroma longirostris U74344 Pterodroma macroptera EU979357 Pterodroma madeira FJ196362 Pterodroma magentae U74338 Pterodroma mollis dubia U74655 Pterodroma mollis mollis U74334 ST HELENA PETREL, PTERODROMA RUPINARUM 5 © 2013 The Linnean Society of London, Zoological Journal of the Linnean Society, 2013 RESULTS DNA sequences were obtained from two samples: 194 bp from one sample from Prosperous Bay (USNM 176607) and 786 contiguous base pairs from the sample from Sandy Bay (USNM 176679). The sequence from Prosperous Bay fully overlapped with the sequence from Sandy Bay and was identi- cal. After alignment, the total sequence length for Pt. rupinarum was 949 bp. This sequence has been deposited in the GenBank database under accession number KF010837. The full alignment of comparative sequences was 1143 bp and contained 76 taxa, with an average sequence length of 1045 bp. No gaps or premature stop codons were detected. There were 478 variable sites with 21, 6, and 73% occurring in the first, second, and third codon positions, respectively. After translation, 14% of the substitutions resulted in non-synonymous changes. The nucleotide frequen- cies were A = 28.15%, T = 26.34%, C = 32.91%, and G = 12.60%, and the transition/transversion bias was estimated to be 9.42. This pattern of variation suggests that these sequences are likely to be of mitochondrial, and not nuclear, origin. Additionally, the independent amplification of multiple overlapp- ing fragments for Pt. rupinarum did not uncover any ambiguous sequences, which provides additional evi- dence for the mitochondrial origin of those sequences. Maximum likelihood and Bayesian inference pro- duced congruent topologies, with Pt. rupinarum placed within a clade of Pterodroma petrels from the Atlantic with high support. For both methods of analysis there was strong support for the codon- partitioned model over the unpartitioned model (LRT = 1254, P < 0.001; log10BF = 238), although both models also produced congruent topologies. Overall, the topology recovered here is generally in agreement with other molecular analyses of Procellariidae (Fig. 2; Nunn & Stanley, 1998; Austin et al., 2004; Penhallurick & Wink, 2004; Jesus et al., 2009; Gangloff et al., 2012). In general the cytochrome b gene provides poor resolution for deeper nodes beyond the level of genus. The branch leading to Pseudobulweria received strong support (posterior probability, PP 1.0; bootstrap support, BS 99). Within Pseudobulweria, relationships were largely unre- solved, and further analyses are particularly needed to determine subspecific relationships and taxonomic status within the Tahiti petrel (Ps. rostrata; Gangloff et al., 2012). Support was also strong for the branch leading to Pterodroma (PP 1.0; BS 100). Within Pterodroma, deeper relationships were unresolved, but some strongly supported clades were also detected. In the phylogenies, Pt. rupinarum did not group with Pseudobulweria, but instead fell within a Table 2. Continued Taxon GenBank accession no. Pterodroma neglecta U74341 Pterodroma nigripennis U74343 Pterodroma phaeopygia HQ420347 Pterodroma rupinarum KF010837 Pterodroma sandwichensis HQ420351 Pterodroma solandri U74347 Pterodroma ultima HQ420379 Puffinus assimilis AF076080 Puffinus bryani HQ589355 Puffinus bulleri AF076081 Puffinus carneipes AF076082 Puffinus creatopus AF076083 Puffinus gravis U74354 Puffinus griseus U74353 Puffinus huttoni AF076084 Puffinus lherminieri AF076085 Puffinus nativitatis AF076086 Puffinus newelli HQ589356 Puffinus opisthomelas AF076087 Puffinus pacificus AF076088 Puffinus puffinus U74355 Puffinus tenuirostris U74352 Thalassoica antarctica AF076095 6 A. J. WELCH ET AL. © 2013 The Linnean Society of London, Zoological Journal of the Linnean Society, 2013 Figure 2. Bayesian consensus tree of cytochrome b sequences from 72 procellariid taxa and four Pelecanoides out-groups. Maximum-likelihood analyses yielded a congruent topology. Bayesian posterior probability and maximum-likelihood bootstrap support are shown along the branches. Sampling locations for Pterodroma arminjoniana and Pseudobulweria rostrata are as listed in Table 2. ST HELENA PETREL, PTERODROMA RUPINARUM 7 © 2013 The Linnean Society of London, Zoological Journal of the Linnean Society, 2013 strongly supported Atlantic clade of Pterodroma (Fig. 2) that includes Pterodroma hasitata (Kuhl, 1820), Pterodroma cahow (Nichols & Mowbray, 1916), and the Macaronesian petrels, Pt. madeira, Pt. deserta, and Pt. feae (PP 1.0; BS 98). Within this clade, Pt. rupinarum was most closely related to Pt. feae (PP 0.99; BS 91), which breeds in the Cape Verde Islands, the most southerly of the Macaronesian islands and geographically closest to St Helena. Divergence time estimates from the strict molecular clock gene tree analysis suggest that Pt. rupinarum diverged from Pt. feae approximately 260 000 years ago (95% highest posterior density 126 000–418 000 years ago). The pairwise genetic dis- tance (Table 3) between Pt. rupinarum and Pt. feae was 1.10%, and between Pt. rupinarum and the next closest taxon, Pt. deserta, was 1.83%, whereas dis- tances to species of Pseudobulweria were much larger, ranging from 12.2 to 14.6%, and averaging approxi- mately 13%. The calculation of genetic distances can be influenced by missing data (Fregin et al., 2012), and therefore these distances may be slightly revised if longer sequences from Pt. rupinarum or the criti- cally endangered Pseudobulweria petrels become available. DISCUSSION PHYLOGENETIC PLACEMENT Phylogenetic analyses of the cytochrome b sequence of the extinct Pt. rupinarum strongly support its place- ment in the genus Pterodroma. It groups within a clade of petrels distributed in the North Atlantic, and is most closely related to the Macaronesian petrels, in particular Pt. feae of the Cape Verde islands. It is unlikely that the placement of Pt. rupinarum could result from contamination, postmortem DNA damage, or poor taxonomic sampling. Negative DNA extraction and PCR controls did not display evidence of contami- nation. Additionally, sequences of multiple, overlap- ping fragments from the same sample were identical, as were sequences obtained from two different indi- viduals collected at different sites under different conditions of local preservation, suggesting minimal contamination or errors from miscoding legions in the DNA. PCR amplification success is higher for shorter fragments, typical of ancient DNA. Taxonomic sam- pling is also not likely to be an issue. Using sequences available in the GenBank database all extant Pseudobulweria petrels were included in the phylog- eny. One additional extinct taxon potentially belong- ing to that genus remains to be described (Gangloff et al., 2012). Cytochrome b sequences > 600 bp were not available for seven Pacific taxa: Pterodroma alba (Gmelin, 1789), Pterodroma atrata (Mathews, 1912), Pterodroma brevipes (Peale, 1848), Pterodroma cervicalis (Salvin, 1891), Pterodroma defilippiana (Giglioli & Salvadori, 1869), Pterodroma leucoptera (Gould, 1844), and Pterodroma pycrofti Falla, 1933, none of which are thought to belong in the same clade as the Atlantic taxa that grouped with Pt. rupinarum. No sequences were available for the extinct Jamaica Table 3. Pairwise genetic distances (%) between Ptero- droma rupinarum and Pterodroma and Pseudobulweria petrels, and other selected taxa of Procellariidae Taxon Pterodroma rupinarum Pterodroma feae 1.10 Pterodroma deserta 1.83 Pterodroma madeira 2.39 Pterodroma mollis dubia 2.53 Pterodroma hasitata 3.25 Pterodroma cahow 3.32 Pterodroma mollis mollis 3.38 Pterodroma magentae 3.42 Pterodroma lessonii 4.11 Pterodroma incerta 4.25 Pterodroma macroptera 4.55 Pterodroma axillaris 7.65 Pterodroma nigripennis 7.67 Pterodroma solandri 8.01 Pterodroma externa 8.16 Pterodroma neglecta 8.94 Pterodroma baraui 9.09 Pterodroma longirostris 9.16 Pterodroma ultima 9.42 Pterodroma arminjoniana – Trindade 9.46 Pterodroma heraldica 9.54 Pterodroma sandwichensis 9.56 Pterodroma hypoleuca 9.58 Pterodroma phaeopygia 9.73 Pterodroma cookii 9.98 Pterodroma inexpectata 10.32 Pseudobulweria rostrata rostrata – Marquesas 12.15 Pseudobulweria rostrata rostrata – Fiji 12.28 Pseudobulweria aterrima 12.36 Pseudobulweria rostrata trouessarti – New Caledonia 12.71 Pseudobulweria macgillivrayi 13.11 Pseudobulweria rostrata rostrata – Gambier 13.56 Pseudobulweria becki 14.64 Bulweria bulwerii 12.72 Fulmarus glacialoides 14.86 Pachyptila desolata 13.93 Pelecanoides magellanicus 15.45 Procellaria parkinsoni 13.26 Puffinus puffinus 14.18 8 A. J. WELCH ET AL. © 2013 The Linnean Society of London, Zoological Journal of the Linnean Society, 2013 petrel, Pterodroma caribbaea Carte, 1866, which has sometimes been considered a subspecies of Pt. hasitata (see below). Analyses relying strictly on mitochondrial DNA can be biased in some circumstances. In cases of hybridi- zation, phylogenetic analysis of mitochondrial DNA will group hybrids with the maternal species. It is possible that the individuals of Pt. rupinarum we sequenced were the product of a Pterodroma × Pseudobulweria hybridization event, although this seems unlikely. There is no evidence for the existence of another species of Pterodroma on St Helena, and no species of Pseudobulweria are known from the Atlantic. Furthermore, continental landmasses as well as non-physical barriers may reduce the fre- quency of dispersal of Pseudobulweria into the Atlan- tic (Steeves, Anderson & Friesen, 2005; Friesen, Burg & McCoy, 2007). Also, hybridization is thought to be rare in procellariid seabirds in general, with only a few reported cases in the family (Hunter, 1982; Brown et al., 2010). Finally, it is unclear whether a hybridi- zation event between such divergent species would produce viable offspring. Positive selection may also act on the mitochondrial genome (Bazin, Glémin & Galtier, 2006), which could potentially influence the results of phylogenetic analyses. Additional data from the nuclear genome would help to verify the findings presented here, but the results are consistent with the geographical distributions of other species of Pterodroma in the Atlantic. Therefore, biases related to hybridization and positive selection do not seem likely. The present analysis indicates, with high support, that Pt. rupinarum belongs within a clade of Pterodroma petrels that are distributed around the North Atlantic. This group includes the endangered black-capped petrel (Pt. hasitata) of Hispaniola, Cuba, and Lesser Antilles, the vanished Jamaica petrel (Pt. caribbaea), which was once considered a subspecies of Pt. hasitata (Brooke, 2004), and was last observed in the mid 19th century, the endangered Bermuda petrel (Pt. cahow), as well as the three species of Macaronesian petrels breeding in the eastern North Atlantic. The closest taxon to Pt. rupinarum appears to be Pt. feae of the Cape Verde Islands, with Pt. deserta of Bugio Island in the Madeira Archipelago being slightly more distant. The species status of Pt. deserta has been debated. There appears to be little or no gene flow with Pt. feae, but the two taxa have diverged only very recently, within the last 40 000 years (Gangloff et al., 2013), and so the divergence process may not be complete. Pseudobulweria petrels breed in the Indian and Pacific oceans, and given this wide distribution it may not have been surprising to find a member of the genus in the Atlantic; however, the phylogenetic placement of Pt. rupinarum within Pterodroma is consistent with a biogeographic scenario in which a closely-knit clade of petrel species radiated through- out the Atlantic and Caribbean, reaching as far south as St Helena (Fig. 1). The phylogenetic placement of Pt. rupinarum may seem contradictory to earlier osteological analyses. The skull of this taxon was described as being similar to petrels of the genus Pseudobulweria, which have a marked downward slope of the bill as compared with Pterodroma petrels, where the bill attaches to the skull in a horizontal manner (Olson, 1975). When these analyses were conducted nearly 40 years ago no skeletons of Pseudobulweria were available, and therefore comparisons were made between the three- dimensional fossils and X-radiographs taken of museum study skins. Similarly, there were very few skeletons of Pterodroma species in museum collec- tions, and none of the Macaronesian species. This is slowly beginning to change. Now there are at least a good series of skeletons of several different populations of Ps. rostrata, and although very scarce and scattered in small and remote museums, there are also a few skeletons now available for the Macaronesian species of Pterodroma. Until proper comparisons of this new material can be made with the fossils of Pt. rupinarum, it remains unclear whether the apparently distinctive features of the skull and mandible are truly similar to, and presum- ably convergent with, Pseudobulweria, or whether such features may also occur in some or all of the Macaronesian species of Pterodroma. Likewise, we are unable at this point to say how much the osteology of Pt. rupinarum is similar to or differs from its apparent closest relative, Pt. feae, of the Cape Verde Islands. SPECIES STATUS Pterodroma rupinarum shows a close relationship to Pt. feae and Pt. deserta, with an uncorrected pairwise sequence divergence of 1.1–1.8%. Further work is needed to investigate genetic variation within these taxa, which may lead to a revision of the estimated genetic distance in the future; however, a distance of 1.1% is within the normal range of pairwise cytochrome b sequence divergences for acknowledged full species of petrels (Austin et al., 2004; Techow, Ryan & O’Ryan, 2009). The lowest pairwise diver- gence between species in the data set analysed here was 0.7%, which occurred between Pachyptila desolata (Gmelin, 1789) and Pachyptila salvini (Mathews, 1912), and between Puffinus creatopus Coues, 1984 and Puffinus carneipes Gould, 1844. Overall, a pairwise sequence divergence of 1.1%, in addition to osteological differences, suggests that ST HELENA PETREL, PTERODROMA RUPINARUM 9 © 2013 The Linnean Society of London, Zoological Journal of the Linnean Society, 2013 Pt. rupinarum is a distinct species from Pt. feae. At the very least it may represent a substantially differ- entiated population of Pt. feae that was isolated for a long period of time prior to its extinction, perhaps about a quarter of a million years, assuming an evolutionary rate of 1.89% per million years for the cytochrome b gene (Weir & Schluter, 2008). EXTINCTION AND CONSERVATION OF PETRELS Several species of procellariid seabirds have become extinct or extirpated since humans colonized the islands where they breed, including Puffinus spelaeus Holdaway & Worthy, 1994 and Puffinus olsoni McMinn, Jaume & Alcover, 1990 (Rando & Alcover, 2008) in the Canary Islands, and Puffinus parvus Shufeldt, 1916 (= P. boydi) (Olson, 2004) on Bermuda, as well as other species in the Pacific. In addition, many bones have been found for which the taxonomic identity remains unclear, and these may or may not represent undescribed species (Steadman, 2006). Furthermore, many populations of extant species have been extirpated as well, and these could have represented distinct but unrecognized evolutionary units. It is interesting that two extinct species occur in the same clade of Pterodroma: Pt. rupinarum on St Helena and Pt. caribbaea on Jamaica (although the extinction of the Jamaica petrel has not been confirmed). Which factors are consistently related to seabird extinctions, if any, remains unclear. Anthropogenic influences on these islands are similar to those of other oceanic islands where procellariid species have persisted. Gangloff et al. (2012) suggest that petrels that breed on gently sloping islands, as opposed to islands with steep, inaccessible cliffs, may be at increased risk for extinction. This may have been the case for Pt. rupinarum (Olson, 1975), but Pt. caribbaea was known to breed on steep cliffs (Brooke, 2004). Another factor may be that Pt. rupinarum and Pt. caribbaea are both single- island endemics, which have an increased risk of extinction because of restricted ranges, small popula- tion sizes, and lower genetic diversity (Frankham, 1997). Additionally, the characteristic philopatry of petrel species may mean that birds continually return to the same high-mortality breeding sites year after year (Igual et al., 2007). Although at least one, and potentially two, instances of extinction have occurred in this clade, there is also at least one case of suc- cessful conservation restoration. Pterodroma cahow was thought to be extinct for over 200 years until a small population was rediscovered in the early 1900s. Through intensive management, population sizes have been steadily increasing and translocated birds have begun breeding in their new colonies (Carlile et al., 2003; Gehrman, 2012; Birdlife International, 2013); therefore, although some petrel extinctions have already occurred, it may still be possible to prevent additional extinctions in the future. ACKNOWLEDGEMENTS We would like to thank Helen James, Nancy Rotzel, and members of the Center for Conservation and Evolutionary Genetics at the Smithsonian Instit- ution’s National Zoological Park for helpful discussion and logistical support. Two reviewers provided useful comments that improved this article. We would like to thank Chris Milensky, Division of Birds, National Museum of Natural History, Smithsonian Institution, for access to samples, and the Center for Computa- tional Resources at the University at Buffalo for access to the high-performance computer cluster. The Center for Conservation and Evolutionary Genetics and NSF-DEB0745604 provided funding for the laboratory work. REFERENCES Altekar G, Dwarkadas S, Huelsenbeck JP, Ronquist F. 2004. Parallel metropolis coupled Markov chain Monte Carlo for Bayesian phylogenetic inference. Bioinformatics 20: 407–415. Attié C, Stahl J-C, Bretagnolle V. 1997. New data on the endangered mascarene petrel Pseudobulweria aterrima: a third twentieth century specimen and distribution. Colonial Waterbirds 20: 406–412. Austin JJ, Bretagnolle V, Pasquet E. 2004. A global molecular phylogeny of the small Puffinus shearwaters and implications for systematics of the Little-Audobon’s shear- water complex. Auk 121: 847–864. Barbraud C, Rolland V, Jenouvrier S, Nevoux M, Delord K, Weimerskirch H. 2012. Effects of climate change and fisheries bycatch on Southern Ocean seabirds: a review. Marine Ecology Progress Series 454: 285–307. Bazin E, Glémin S, Galtier N. 2006. Population size does not influence mitochondrial genetic diversity in animals. Science 312: 570–572. Birdlife International. 2013. Species factsheet: Pterodroma cahow. Bourne WRP. 1955. The birds of the Cape Verde Islands. Ibis 97: 508–556. Bourne WRP. 1956. Notes on a skull of the genus Bulweria from St. Helena. Bulletin of the British Ornithologists’ Club 76: 79–82. Bourne WRP. 1965. The missing petrels. Bulletin of the British Ornithologists’ Club 85: 97–105. Bretagnolle V, Attié C, Pasquet E. 1998. Cytochrome-B Evidence for validity and phylogenetic relationships of Pseudobulweria and Bulweria (Procellariidae). Auk 115: 188–195. Brooke M. 2004. Albatrosses and Petrels across the world. Oxford: Oxford University Press. 10 A. J. WELCH ET AL. © 2013 The Linnean Society of London, Zoological Journal of the Linnean Society, 2013 Brown RM, Nichols RA, Faulkes CG, Jones CG, Bugoni L, Tatayah V, Gottelli D, Jordan WC. 2010. Range expansion and hybridization in Round Island petrels (Pterodroma spp.): evidence from microsatellite genotypes. Molecular Ecology 19: 3157–3170. Carlile N, Priddle D, Zino F, Natividad C, Wingate DB. 2003. A review of four successful recovery programmes for threatened sub-tropical petrels. Marine Ornithology 31: 185–192. Croll DA, Maron JL, Estes JA, Danner EM, Byrd GV. 2005. Introduced predators transform subarctic islands from grassland to tundra. Science 307: 1959–1961. Croxall JP, Butchart SHM, Lascelles B, Stattersfield AJ, Sullivan B, Symes A, Taylor P. 2012. Seabird conserva- tion status, threats and priority actions: a global assess- ment. Bird Conservation International 22: 1–34. Drummond AJ, Suchard MA, Xie D, Rambaut A. 2012. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Molecular Biology and Evolution 29: 1969–1973. Fleischer RC, Olson SL, James HF, Cooper AC. 2000. Identification of the extinct Hawaiian Eagle (Haliaeetus) by mtDNA sequence analysis. Auk 117: 1051–1056. Fontaine R, Gimenez O, Bried J. 2011. The impact of introduced predators, light-induced mortality of fledglings and poaching on the dynamics of the Cory’s shearwater (Calonectris diomedea) population from the Azores, north- eastern subtropical Atlantic. Biological Conservation 144: 1998–2011. Frankham R. 1997. Do island populations have less genetic variation than mainland populations? Heredity 78: 311–327. Fregin S, Haase M, Olsson U, Alström P. 2012. Pitfalls in comparisons of genetic distances: a case study of the avian family Acrocephalidae. Molecular Phylogenetics and Evolu- tion 62: 319–328. Friesen VL, Burg TM, McCoy KD. 2007. Mechanisms of population differentiation in seabirds. Molecular Ecology 16: 1765–1785. Gangloff B, Shirihai H, Watling D, Cruaud C, Couloux A, Tillier A, Pasquet E, Bretagnolle V. 2012. The complete phylogeny of Pseudobulweria, the most endangered seabird genus: systematics, species status and conservation impli- cations. Conservation Genetics 13: 39–52. Gangloff B, Zino F, Shirihai H, González-Solís J, Couloux A, Pasquet E, Bretagnolle V. 2013. The evolu- tion of north-east Atlantic gadfly petrels using statistical phylogeography. Molecular Ecology 22: 495–507. Gehrman E. 2012. Rare birds: the extraordinary tale of the Bermuda petrel and the man who brought it back from extinction. Boston: Beacon Press. Gosse P. 1938. St. Helena, 1502–1938. London: Cassell and Company. Holdaway RN, Worthy TH. 1994. A new fossil species of shearwater Puffinus from the late quaternary of the South Island, New Zealand, and notes on the biogeography and evolution of the Puffinus gavia superspecies. Emu 94: 201– 215. Hunter S. 1982. Inter-specific breeding in Giant petrels at South Georgia. Emu 82: 608–614. Igual JM, Forero MG, Gomez T, Oro D. 2007. Can an introduced predator trigger an evolutionary trap in a colo- nial seabird? Biological Conservation 137: 189–196. Jesus J, Menezes D, Gomes S, Oliveira P, Nogales M, Brehm A. 2009. Phylogenetic relationships of gadfly petrels Pterodroma spp. from the Northeastern Atlantic Ocean: molecular evidence for specific status of Bugio and Cape Verde petrels and implications for conservation. Bird Con- servation International 19: 199–214. Lewison RL, Crowder LB. 2003. Estimating fishery bycatch and effects on a vulnerable seabird population. Ecological Applications 13: 743–753. Maddison DR, Maddison WP. 2008. MacClade v 4.08. Sunderland, MA: Sinauer Associations. Nei M, Kumar S. 2000. Molecular evolution and phylogenetics. New York: Oxford University Press. Nunn GB, Stanley SE. 1998. Body size effects and rates of cytochrome b evolution in tube-nosed seabirds. Molecular Biology and Evolution 15: 1360–1371. Olson SL. 1975. Paleornithology of St. Helena Island, South Atlantic Ocean. Smithsonian Contributions to Paleobiology 23: 1–49. Available at: http://www.sil.si.edu/ smithsoniancontributions/paleobiology/pdf_lo/SCtP-0023.pdf Olson SL. 2004. Taxonomic review of the fossil Procellariidae (Aves: Procellariiformes) described from Bermuda by R. W. Shufeldt. Proceedings of the Biological Society of Washing- ton 117: 575–581. Penhallurick J, Wink M. 2004. Analysis of the taxonomy and nomenclature of the Procellariiformes based on com- plete nucleotide sequences of the mitochondrial cytochrome b gene. Emu 204: 125–147. Piatt JF, Sydeman WJ, Wiese F. 2007. Introduction: a modern role for seabirds as indicators. Marine Ecology Progress Series 352: 199–204. Posada D. 2008. jModelTest: phylogenetic model averaging. Molecular Biology and Evolution 25: 1253–1256. Pyle P, Welch AJ, Fleischer RC. 2011. A new species of shearwater (Puffinus) recorded from Midway Atoll, North- western Hawaiian Islands. Condor 113: 518–527. Rando JC, Alcover JA. 2008. Evidence for a second western Palaearctic seabird extinction during the last millenium: the Lava Shearwater. Ibis 150: 188–192. Rheindt FE, Austin JJ. 2005. Major analytical and concep- tual shortcomings in a recent taxonomic revision of the Procellariiformes – a reply to Penhallurick and Wink (2004). Emu 105: 181–186. Ronquist F, Huelsenbeck JP. 2003. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574. Sánchez-Piñero F, Polis GA. 2000. Bottom-up dynamics of allochthonous input: direct and indirect effects of seabirds on islands. Ecology 81: 3117–3132. Shirihai H. 2008. Rediscovery of Beck’s petrel Pseudo- bulweria becki, and other observations of tubenoses from the Bismark archipelago, Papua New Guinea. Bulletin of the British Ornithologists’ Club 128: 3–16. Shirihai H, Pym T, Kretzschmar J, Moce K, Taukei A, Watling D. 2009. First observations of Fiji Petrel ST HELENA PETREL, PTERODROMA RUPINARUM 11 © 2013 The Linnean Society of London, Zoological Journal of the Linnean Society, 2013 Pseudobulweria macgillivrayi at sea: off Gau Island, Fiji, in May 2009. Bulletin of the British Ornithologists’ Club 129: 129–148. Sorenson MD, Fleischer RC. 1996. Multiple independent transpositions of mitochondrial DNA control region sequences to the nucleus. Proceedings of the National Academy of Sciences of the United States of America 93: 15239–15243. Stamatakis A, Hoover P, Rougemont J. 2008. A rapid bootstrap algorithm for the RAxML web-servers. Systematic Biology 75: 758–771. Steadman DW. 1995. Island birds: biodiversity meets zooarchaeology. Science 267: 1123–1131. Steadman DW. 2006. Extinction and biogeography of tropical Pacific birds. Chicago: University of Chicago Press. Steeves TE, Anderson DJ, Friesen VL. 2005. The Isthmus of Panama: a major physical barrier to gene flow in a highly mobile pantropical seabird. Journal of Evolutionary Biology 18: 1000–1008. Suchard MA, Weiss RE, Sinsheimer JS. 2001. Bayesian selection of continuous-time markov chain evolutionary models. Molecular Biology and Evolution 18: 1001–1013. Swofford DL. 2002. PAUP*. Phylogenetic analysis using par- simony (*and other methods). Sunderland, MA: Sinauer Associates. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. 2011. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Molecular Biology and Evolution 28: 2731–2739. Techow NMSM, Ryan PG, O’Ryan C. 2009. Phylo- geography and taxonomy of White-chinned and Spectacled Petrels. Molecular Phylogenetics and Evolution 52: 25– 33. Wanless RM, Angel A, Cuthbert RJ, Hilton GM, Ryan PG. 2007. Can predation by invasive mice drive seabird extinctions? Biology Letters 3: 241–244. Weimerskirch H, Brothers N, Jouventin P. 1997. Popu- lation dynamics of wandering albatross Diomedea exulans and Amsterdam albatross D. amsterdamensis in the Indian Ocean and their relationships with long-line fisheries: conservation implications. Biological Conservation 79: 257– 270. Weir JT, Schluter D. 2008. Calibrating the avian molecular clock. Molecular Ecology 17: 2321–2328. Welch AJ, Fleischer RC, James HF, Wiley AE, Ostrom PH, Adams J, Duvall F, Holmes N, Hu D, Penniman J, Swindle K. 2012. Population divergence and gene flow in an endangered and highly mobile seabird. Heredity 109: 19–28. Wiley AE, Ostrom PH, Welch AJ, Fleischer RC, Gandhi H, Southon J, Stafford TW Jr, Penniman J, Hu D, Duvall F, James HF. 2013. Millennial-scale records from a wide-ranging predator show evidence of recent human impact to oceanic food webs. Proceedings of the National Academy of Sciences of the United States of America 110: 8972–8977. Zino F, Brown R, Biscoito M. 2008. The separation of Pterodroma madeira (Zino’s petrel) from Pterodroma feae (Fea’s petrel) (Aves: Procellariidae). Ibis 150: 326–334. 12 A. J. WELCH ET AL. © 2013 The Linnean Society of London, Zoological Journal of the Linnean Society, 2013