Hostname: page-component-7c8c6479df-ws8qp Total loading time: 0 Render date: 2024-03-27T23:30:34.869Z Has data issue: false hasContentIssue false

Phylogenetic placement and life history trait imputation for Grenada Dove Leptotila wellsi

Published online by Cambridge University Press:  05 August 2022

Catherine Peters*
Affiliation:
University of Chester, Biological Sciences, Parkgate Road, Chester CH1 4BJ, UK
Matthew Geary
Affiliation:
University of Chester, Biological Sciences, Parkgate Road, Chester CH1 4BJ, UK
Howard P. Nelson
Affiliation:
University of Cambridge, Pembroke Street, Cambridge, CB2 3QZ, UK
Bonnie L. Rusk
Affiliation:
Grenada Dove Conservation Programme, c/o Forestry and National Parks Department, Saint George, Grenada
Achaz Von Hardenberg
Affiliation:
University of Chester, Biological Sciences, Parkgate Road, Chester CH1 4BJ, UK
Anna Muir
Affiliation:
University of Chester, Biological Sciences, Parkgate Road, Chester CH1 4BJ, UK
*
*Author for correspondence: Catherine Peters, Email: c.peters@chester.ac.uk
Rights & Permissions [Opens in a new window]

Abstract

Phylogenetic analyses can be used to resolve taxonomic uncertainties and reconstruct a species’ evolutionary history. This can be combined with ecological data to predict missing life history traits which are important for creation of conservation management strategies. We investigated the evolutionary and life history of the ‘Critically Endangered’ Grenada Dove Leptotila wellsi by estimating its phylogenetic placement and using this new phylogeny to test the accuracy of phylogenetic comparative methods for estimating both documented and unknown life history traits. We extracted DNA from two Grenada Dove samples and obtained sequences from three mitochondrial markers: Cytochrome oxidase I (COI), NADH dehydrogenase 2 (ND2) and Cytochrome b (Cyt b); and one nuclear marker: β-Fibrinogen intron 7 (β-FIB). We present the first genetic data obtained for the Grenada Dove. Our data identify the Grey-Chested Dove Leptotila cassinii as the species which shares both a most recent common ancestor, with an estimated divergence of approximately 2.53 million years ago, and the smallest genetic distance (P = 0.0303) with the Grenada Dove. Life history trait values for the Grenada Dove predicted from our analyses using phylogenetic imputation are: clutch size = 2 (± 0.09) eggs; clutches per year = 1.4 (± 0.81); incubation time = 14.2 (± 0.75) days; hatching weight = 3.8 g (± 1.05) and single imputation: fledging age (genus median) = 15.5 days, longevity (genus median) = 8.6 years. This study contributes novel information regarding evolutionary history and life history characteristics to inform long-term conservation actions for a ‘Critically Endangered’ species.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2022. Published by Cambridge University Press on behalf of BirdLife International

Introduction

Over 12% of extant bird species are threatened with extinction (Brooks et al. Reference Brooks, Collar, Green, Marsden and Pain2008), 58.8% of which are island species (Tershy et al. Reference Tershy, Shen, Newton, Holmes and Croll2015). Natural extinction rates are higher in island populations than in mainland populations. Sensitivity to environmental change, predation by introduced predators, and small population sizes are all common characteristics among island avifaunas; therefore, island endemic species are at an increased risk of extinction (Frankham Reference Frankham1998). Taxonomic patterns in extinction risk suggest that families containing high numbers of endemic species are more susceptible to extinction (Lockwood et al. Reference Lockwood, Brooks and Mckinney2000). Resolution of unknown taxonomic and evolutionary relationships can therefore contribute to a greater understanding of extinction vulnerabilities among bird groups (Lockwood et al. Reference Lockwood, Brooks and Mckinney2000, Johnson et al. Reference Johnson, Adler and Cherry2007).

Advances in systematic biology have allowed phylogenetic information to be combined with phenotypic data to further infer the evolutionary processes, based on the assumption that phenotypic traits within groups of species will be influenced by shared ancestry (Cressler et al. Reference Cressler, Butler and King2015, Bastide et al. Reference Bastide, Solís-Lemus, Kriebel, Sparks and Ané2018). Recently, comparative phylogenetic methods for missing data imputation have been employed to predict physiological adaptations, conservation status, and functional life history traits for data deficient species (Riek and Bruggeman Reference Riek and Bruggeman2013, Swenson Reference Swenson2014, González-del-Pliego et al. Reference González-del-Pliegoy, Freckleton, Edwards, Koo, Scheffers, Pyron and Jetz2019, James et al. Reference James, Salguero-Gómez, Jones, Childs and Beckerman2021). Predictive models are particularly useful for estimation of trait values when species are rare and elusive, locations are remote, and fiscal constraints limit field surveys (Wood et al. Reference Wood, Stillman and Hilton2018, Horswill et al. Reference Horswill, Manica, Daunt, Newell, Wanless, Wood and Matthiopoulos2021). Life history traits such as generation time, longevity, male and female maturity, clutch size, and incubation and fledging times can be used to evaluate the capability of a species to adapt to climate change and habitat fragmentation as well as being essential to estimate population trends (Mace et al. Reference Mace, Collar, Gaston, Hilton-Taylor, Akçakaya, Leader-Williams, Milner-Gulland and Stuart2008, Pearson et al. Reference Pearson, Stanton, Shoemaker, Aiello-Lammens, Ersts, Horning, Fordham, Raxworthy, Ryu, Mcnees and Akçakaya2014, Storchová and Hořák Reference Storchová and Hořák2018, Cuervo and Møller Reference Cuervo and Møller2019, Horswill et al. Reference Horswill, Manica, Daunt, Newell, Wanless, Wood and Matthiopoulos2021).

Pigeons and doves (Columbiformes; Columbidae) are one of the oldest lineages of birds, inhabiting a wide range of habitats on six of the seven continents, despite which their evolutionary history is still poorly resolved (Pereira et al. Reference Pereira, Johnson, Clayton and Baker2007, Soares et al. Reference Soares, Novak, Haile, Heupink, Fjeldså, Gilbert, Poinar, Church and Shapiro2016). Neotropical columbids may form the most ancient lineage (Johnson and Clayton Reference Johnson and Clayton2000b, Shapiro et al. Reference Shapiro, Sibthorpe, Rambaut, Austin, Wragg, Bininda-Emonds, Lee and Cooper2002) but have been little studied, with fundamental taxonomy, conservation status, ecological and life-history data still largely unknown (Brooks et al. Reference Brooks, Collar, Green, Marsden and Pain2008, Devenish-Nelson et al. Reference Devenish-Nelson, Weidemann, Townsend and Nelson2019). Therefore, gathering information about Neotropical species is essential and urgent for conservation of this diverse group of birds (Latta Reference Latta2012).

The Neotropical columbid genus Leptotila comprises 11 morphologically similar species, of which one, the Grenada Dove L. wellsi, is ‘Critically Endangered’ (BirdLife International 2018). Endemic to Grenada in the eastern Caribbean Sea, the Grenada Dove is found in two seemingly isolated populations in the west and south-west of the island (Rusk et al. Reference Rusk, Temple, Wege, Wallace, Blockstein, Ledec and Sorenson2008, Rusk Reference Rusk2017). The most recent population assessment estimated that 160 ± 30 individuals remain of this species and showed a decreasing population trajectory (Rivera-Milán et al. Reference Rivera-Milán, Bertuol, Simal and Rusk2015, BirdLife International 2018). Important details about the species’ biology, such as longevity and incubation period, are currently unreported. The Grenada Dove was once considered a subspecies of the Grey-fronted Dove L. rufaxilla until reclassified as a distinct species based on morphological and acoustic evidence (Blockstein and Hardy Reference Blockstein and Hardy1989). However, its genetic distinctiveness and relationships have never been investigated. Anecdotal evidence indicates that the Grenada Dove is closely related to the Caribbean Dove L. jamaicensis due to morphological similarities (Anthony Jeremiah, Head of Forest and National Parks Department, Government of Grenada, pers. comm.).

By investigating both the evolutionary and life history of the Grenada Dove, and by estimating its phylogenetic placement, we were able to predict currently unreported life history traits. We extracted DNA from Grenada Dove samples and obtained sequences from three mitochondrial DNA markers: Cytochrome oxidase I (COI), NADH dehydrogenase 2 (ND2) and Cytochrome b (Cyt b); and one nuclear marker: β-Fibrinogen intron 7 (β-FIB). We used these markers first to infer the evolutionary history of this species amongst other Neotropical columbid species. We then used phylogenetic comparative methods for missing data imputation, based on this evolutionary history, to predict unknown life history trait values to inform conservation management decisions.

Methods

DNA Extraction and PCR

Non-invasively collected eggshell and feather samples (n = 2) were collected in 2018 and 2017, respectively. Samples were transported by airmail to the UK in June 2018 where they were cleaned with 70% ethanol and stored at −20°C. DNA barcoding was used to confirm species identification for the samples using methods outlined in Patel et al. (Reference Patel, Waugh, Millar and Lambert2010). As this is the first time Grenada dove DNA has been obtained, there is no reference available for this species. However, as it is the only member of its genus resident on Grenada, where all of our samples were collected, we conducted a BLAST® search with the resulting data and found that sequences obtained for this study were consistent with previously obtained sequences for Leptotila. We used the optimised method for DNA extraction, using the QIAGEN DNeasy® Blood and Tissue kit and the primerless PCR protocol, outlined in Peters et al. (Reference Peters, Nelson, Rusk and Muir2019).

The substrate from the primerless PCR process was used to amplify the mitochondrial markers Cyt b (882 bp), COI (613 bp) and ND2 (949 bp) and the nuclear marker β-FIB (902 bp). Primer sequences and the corresponding annealing temperatures are given in Supporting Information 1. PCR cycling parameters for all except the COI region were: initial denaturation at 95°C for 5 min, 45 cycles of 95°C for 30 s, 50–60°C (primer specific, see Appendix S1 in the online supplementary material) for 30 s, 72°C for 60 s and a final extension at 72°C for 5 min; and for the COI region: 95°C for 5 min, 10 cycles of 95°C for 30 s, 55°C for 30 s, 72°C for 60 s, 35 cycles of 95°C for 30 s, 50°C for 30 s, 72°C for 60 s and a final extension at 72°C for 5 min. Where regions could not be amplified using a single primer set, additional primers were designed using Primer3 (Untergasser et al. Reference Untergasser, Cutcutache, Koressaar, Ye, Faircloth, Remm and Rozen2012) allowing a small amplicon strategy to be used where regions were amplified in smaller overlapping sections (as per Peters et al. Reference Peters, Nelson, Rusk and Muir2019).

Prior to sequencing, samples were purified as per the manufacturer’s protocol using the QIAquick PCR Purification kit (QIAGEN Inc., Crawley). Preparation and submission of samples were carried out according to the Eurofins Genomics guidelines for the Mix2Seq kit (Eurofins Genomics, Luxembourg) where sequencing was performed using Sanger sequencing methods. As the primerless PCR technique has the potential to increase the chance of sequencing error, amplification and sequencing were performed in triplicate in order to obtain consensus data for each region (Peters et al. Reference Peters, Nelson, Rusk and Muir2019, Weber et al. Reference Weber, Stewart, Garza and Lehman2000). Sequence reconstruction and alignments were conducted using Sequencher 5.4.6 (Gene Codes Corporation 2017).

Phylogenetics

Sequences from 24 Neotropical columbids with Common Crane Grus grus as an outgroup were obtained for all markers from Genbank (Benson et al. Reference Benson, Cavanaugh, Clark, Karsch-Mizrachi, Lipman, Ostell and Sayers2013; Table 1). Phylogenetic tree construction was performed using MEGAX 10.2.4 (Kumar et al. Reference Kumar, Stecher, Li, Knyaz and Tamura2018). We used Maximum Likelihood to fit 24 different nucleotide substitution models and used Akaike information criterion (AIC) and Bayesian information criterion (BIC) to select the best models (Schwarz Reference Schwarz1978, Akaike Reference Akaike1987, Nei and Kumar Reference Nei and Kumar2000).

Table 1. List of species and their Genbank accession number for all sequences used in phylogenetic analysis.

Sources:

3 (Johnson & Clayton, Reference Johnson and Clayton2000a)

4 (Johnson & Clayton, Reference Johnson and Clayton2000b)

6 (Johnson & Weckstein, Reference Johnson and Weckstein2011)

17 (Prychitko and Moore Reference Prychitko and Moore2003)

24 (Stoeckle Reference Stoeckle2011)

25 (Sweet and Johnson Reference Sweet and Johnson2015)

The Hasegawa-Kishino-Yano (HKY) model using a discrete Gamma distribution (+G) was the best-supported model for the combined mitochondrial and nuclear DNA dataset of 3,346 bp and thus evolutionary analyses were performed using this model with 500 bootstrap replicates (Hasegawa et al. Reference Hasegawa, Kishino and Yano1985). In order to check the consistency of this tree we also built the tree using the *BEAST approach implemented in BEAST v2.6.6 (Heled and Drummond Reference Heled and Drummond2010, Bouckaert et al. Reference Bouckaert, Vaughan, Barido-Sottani, Duchêne, Fourment, Gavryushkina, Heled, Jones, Kühnert, De Maio, Matschiner, Mendes, Müller, Ogilvie, Du Plessis, Popinga, Rambaut, Siveroni, Suchard, Wu, Xie, Zhang, Stadler and Drummond2019) and the resulting species tree provided the same topology as the concatenated tree (for details see Appendix S2). As mtDNA tends to have a larger proportion of variable bases than nDNA (Allio et al. Reference Allio, Donega, Galtier and Nabholz2017), we also inferred the position of the Grenada Dove within the Leptotila genus based on the mtDNA and nDNA separately, to provide a comparison between the two marker types, using the following best-supported models as per AIC and BIC: General Time Reversible model ([+I], 38.01% sites), for a total of 2,444 bp of mtDNA (Tavaré Reference Tavaré1986), and Tamura 3-parameter model, for a total of 902 bp of nDNA (Nei and Kumar Reference Nei and Kumar2000). We also inferred the position of the Grenada Dove in relation to L. jamaicensis, using a reduced genetic dataset using the Hasegawa-Kishino-Yano model (+G, parameter = 0.2496) for a total of 2,720 bp of ND2, Cyt b and β-FIB.

Time-calibrated phylogenies were reconstructed to infer molecular-based estimates of divergence by applying the RelTime-ML method following the protocol outlined in Mello (Reference Mello2018). Fossil calibrations are widely used for molecular dating, but no fossil records are available for the taxa used in this study (Ksepka et al. Reference Ksepka, Parham, Allman, Benton, Carrano, Cranston, Donoghue, Head, Hermsen, Irmis, Joyce, Kohli, Lamm, Leehr, Patané, Polly, Phillips, Smith, Smith, Van Tuinen, Ware and Warnock2015, Peters and McClennen Reference Peters and McClennen2015). We therefore used confidence intervals of minimum and maximum boundary estimations from the TimeTree database, which provides divergence time estimates from published studies integrated with a geological timescale, to calibrate divergence time (Kumar et al. Reference Kumar, Stecher, Suleski and Hedges2017, Mello Reference Mello2018, Tamura et al. Reference Tamura, Battistuzzi, Billing-Ross, Murillo, Filipski and Kumar2012, Reference Tamura, Tao and Kumar2018).

Life history trait imputation

To predict life history traits for the Grenada Dove we used multivariate phylogenetic comparative methods for missing data imputation implemented in the Rphylopars package version 0.3.0 (Goolsby et al. Reference Goolsby, Bruggeman and Ané2017) in R version 3.6.2 (R Core Team 2020). This approach uses a multivariate procedure of phylogenetic covariance across species to predict differences in life history traits (Goolsby et al. Reference Goolsby, Bruggeman and Ané2017, James et al. 2020). We imputed phylogenetic covariance across 24 Neotropical columbid species. Twelve life history traits - clutch size, clutches per year, incubation time, hatching weight, fledging age, longevity, generation length, male and female maturity, no-sex body mass (unreported sex; Myhrvold et al. Reference Myhrvold, Baldridge, Chan, Sivam, Freeman and Ernest2015), and male and female body mass - were obtained for each of these species, where available, from primary and grey literature (Appendix S3; Myhrvold et al. Reference Myhrvold, Baldridge, Chan, Sivam, Freeman and Ernest2015, HBW Alive 2020).

Available life history trait values for the Grenada Dove were a clutch size of two eggs (Rusk, unpubl. data), an imputed generation length of 4.2 years calculated using the IUCN generation length calculator (BirdLife International 2018), a fledging age of 14 days from focal nest observations and an estimate of two or three nesting attempts per year (Rusk, unpubl. data). We imputed these traits to determine whether multivariate phylogenetic comparative methods generated similar values. We fitted evolutionary models and used AIC and BIC as well as assessing model parameters to select the best fitting model. Phylogenetic signals were calculated using a maximum-likelihood approach to quantify the relationship between life history trait values and phylogenetic signal in order to ascertain model precision per life history trait (Penone et al. Reference Penone, Davidson, Shoemaker, Marco, Brooks, Young, Graham, Costa, Di Marco, Rondinini, Brooks, Young, Graham and Costa2014, James et al. Reference James, Salguero-Gómez, Jones, Childs and Beckerman2021). We tested the accuracy of the RPhylopars method by predicting known values, which we removed individually from the dataset. We focused on the genus Columbina as this genus contained the most known data, which allowed us to sequentially remove a total of 20% of these known values for each trait with minimal decrease in the accuracy of predictions (Penone et al. Reference Penone, Davidson, Shoemaker, Marco, Brooks, Young, Graham, Costa, Di Marco, Rondinini, Brooks, Young, Graham and Costa2014). We compared the predicted values from RPhylopars to the true value to assess the accuracy of predictions. We also calculated the median values, used in some life-history imputation without explicit phylogenetic information (e.g. Bird et al. Reference Bird, Martin, Akçakaya, Gilroy, Burfield, Garnett, Symes, Taylor, Şekercioğlu and Butchart2020), for fledging age and longevity.

Results

Phylogenetics

There were no polymorphisms observed between the two Grenada Dove sequences. There were 27 nucleotide polymorphisms across 3,346 bp that were unique to the Grenada Dove and not present in any other Leptotila species used in this investigation, consisting of 21 transition mutations, two transversion mutations and one transition/transversion mutation (depending on the Leptotila species) in the mtDNA and three transitional mutations in the nDNA β-FIB gene. The phylogenetic tree constructed using concatenated sequence data - Cyt b (882 bp), COI (613 bp), ND2 (949 bp), and β-FIB (902 bp) - for Neotropical columbids indicates that the Grenada Dove falls within the monophyletic group formed by Leptotila species (Figure 1a.). The Grenada Dove shared a most recent common ancestor with its sister group containing L. cassinii and L. plumbeiceps approximately 2.53 million years ago (mya) (Figure 2a). The smallest genetic distance (P = 0.0303) where 100 nucleotides were polymorphic (Table 2) was identified between the Grenada Dove and L. cassinii. The phylogenetic tree built with the reduced data and containing L. jamaicensis (Figure 3) still shows a sister relationship between the Grenada Dove and L. cassinii and L. plumbeiceps.

Figure 1. Phylogenetic analysis for 24 Neotropical columbid species and outgroup using the Maximum Likelihood method. Grenada Dove sequences indicated by the circular symbol; a. all mitochondrial and nuclear sequence data estimated using the Hasegawa-Kishino-Yano model (+G, parameter = 0.2281); b. mtDNA using General Time Reversible model ([+I], 38.01% sites); c. β-fibrinogen using the Tamura 3-parameter model. The trees are drawn to scale, with branch lengths measured in the number of substitutions per site and nodes supported by bootstrap values of >80%.

Figure 2. Time-calibrated phylogenies for Leptotila species and outgroup using the RelTime-ML method a. all mitochondrial and nuclear sequence data estimated using the Hasegawa-Kishino-Yano model (+G, parameter = 0.1547), b. mtDNA General Time Reversible model ([+I], 35.97% sites) c. β-fibrinogen using the Tamura 3-parameter model. Number scale indicates million years ago (mya). Nodes indicate 7 estimated divergence time calibration points of estimated divergence time ranges in million years ago (mya) obtained using the TimeTree resource (Kumar et al. Reference Kumar, Stecher, Suleski and Hedges2017).

Table 2. The number of substitutions per site (p-distance) and the number of base differences (polymorphic sites) between sequences for a total of 3,346 bp positions in the final dataset. *2,720 bp positions in the final dataset.

Figure 3. Phylogenetic analysis for 25 Neotropical columbid species and outgroup using the Maximum Likelihood method. Grenada Dove sequences indicated by the circular symbol and Leptotila jamaicensis are indicated by a square using the Hasegawa-Kishino-Yano model (+G, parameter = 0.2496). Tree drawn to scale, with branch lengths measured in the number of substitutions per site and nodes supported by bootstrap values of >80%.

The mitochondrial phylogenetic trees shown in Figures 1b and 2b suggests that the Grenada Dove shared its most recent common ancestor with L. cassinii and L. plumbeiceps approximately 2.18 mya, which is consistent with the phylogenetic tree constructed with the concatenated data set in Figure 1a. The nuclear phylogenetic trees shown in Figures 1c and 2c suggest that the Grenada Dove shared its most recent common ancestor with L. megalura and L. ruffaxilla approximately 2.75 mya, dissimilar to the most recent common ancestor shown by the phylogenetic trees constructed with the combined (Figure 1a and 2a) and mtDNA (Figure 1b and 2b) datasets. Analysis of the β-FIB gene relevels an indel (Appendix S4) in L. cassinii and L. plumbeiceps sequences which is not shared by the Grenada Dove or any of the other Leptotila species included in this study.

Life history trait imputation

Strong phylogenetic signals were revealed for clutch size, no-sex body mass, and male and female body mass, whereas weak phylogenetic signal was detected for clutches per year, incubation time, hatching weight, fledging age, longevity, and generation length (Appendix S5). Model validation (Figure 4) showed lower variation relative to each trait around seven of the twelve life history trait predictions including clutch size, clutches per year, incubation time, hatching weight, fledging age, longevity, and generation length and were therefore considered accurate. However, the remaining five life history trait predictions including male and female maturity, no-sex body mass, and male and female body mass showed a high level of variation from the true value and therefore were discounted. Despite suggested trait accuracy from the model validation for fledging age and longevity, we report large standard error values in addition to weak phylogenetic signal, and as such, genus medians were considered a more accurate representation for these traits. Known values for generation length for Leptotila are themselves imputed using the IUCN generation length calculator, and thus we discounted this trait from our results. Overall, six life history trait values were predicted for the Grenada Dove with confidence (Figure 5) using multivariate and single imputation procedures: a. clutch size = 2 (±0.09 SE) eggs, b. clutches per year = 1.4 (±0.81 S.E.), c. incubation time=14.2 (±0.75 SE) days. d. hatching weight = 3.8 g (± 1.05 SE) grams, fledging age = 15.5 days (genus median) days, f. longevity = 8.6 years (genus median).

Figure 4. Plot displaying the squared difference between the eb model predicted value and the true value divided by the trait median using the genus Columbina along with Normalized Root Mean Squared Error [NRMSE (Penone et al. Reference Penone, Davidson, Shoemaker, Marco, Brooks, Young, Graham, Costa, Di Marco, Rondinini, Brooks, Young, Graham and Costa2014, James et al. Reference James, Salguero-Gómez, Jones, Childs and Beckerman2021)]. Raw data are represented by the ° and the first, second and third quartile are represented by the box plot. Mean Squared Error per life history trait.

Figure 5. Predicted trait values with upper and lower 95% confidence intervals for the Grenada Dove Leptotila wellsi in comparison with other Leptotila species for the following lift history traits using the EB model: a. clutch size, b. clutches per year, c. incubation time, d. hatching weight; and median estimates: e. fledging age, f. longevity. Species of interest highlighted red. Known trait values represented by individual grey point.

Discussion

This study reports the first genetic data obtained for the Grenada Dove and supports a sister relationship with Grey-chested Dove L. cassinii along with Grey-headed Dove L. plumbeiceps, whose ranges span Central America into Colombia (BirdLife International 2021). We provide estimations for six previously unpublished life history traits for the Grenada Dove some of which, such as clutch size and clutches per year, have been found to be correlated with extinction risk (Parlato et al. Reference Parlato, Armstrong and Innes2015).

Leptotila cassinii was the species with the smallest genetic distance from the Grenada Dove (P = 0.0303) with 100 bp nucleotide polymorphisms based on concatenated sequence data and phylogenetic reconstruction. The most recent common ancestor was shared approximately 2.53 mya between L. cassinii, L. plumbeiceps, and the Grenada Dove. However, when the mtDNA and nDNA are considered separately our results show disparity. The mtDNA shows the same relationship displayed by the concatenated data described above, but nDNA analysis revealed that the most recent common ancestor was shared between L. rufaxilla, L. megalura, and the Grenada Dove. The mito-nuclear incongruence may have arisen for a variety of reasons; a larger amount of mtDNA sequence data was used providing a stronger phylogenetic signal, the differing characteristics between these markers such as differential lineage sorting and linkage disequilibrium, all of which can lead to ambiguous patterns of variation and evolutionary inferences (Hurst and Jiggins Reference Hurst and Jiggins2005, Rubinoff and Holland Reference Rubinoff and Holland2005). Overall, we identify L. cassinii and L. plumbeiceps as the species sharing the most recent common ancestor with the Grenada Dove, however, we recommend caution due to the disparity between the mitochondrial or nuclear phylogeny. Nevertheless, our result was unexpected as L. cassinii is predominantly a Central American species and has a larger geographical distance from Grenada than L. rufaxilla, of which the Grenada Dove was previously thought to be a subspecies. We reveal a genetic difference (P = 0.04) and 126 nucleotide polymorphisms between the Grenada Dove and L. rufaxilla. Our findings support the work of Blockstein and Hardy (Reference Blockstein and Hardy1989) whose morphological analysis (showing differing morphological traits such as cinnamon underwing on primary feathers, a greater extent of white abdomen and a lesser extent of white on the tail tips) and sonographic data (no response was made by the Grenada Dove to L. rufaxilla playback presentations) led to the reclassification of the Grenada Dove as a distinct taxon. Anecdotal evidence suggested the Grenada Dove is closely related to L. jamaicensis (Anthony Jeremiah, Head of Forest and National Parks Department, Government of Grenada, pers. comm.), however our results show that the Grenada Dove shares common ancestors with four other Leptotila more recently than its shared ancestor with L. jamaicensis and with a greater genetic distance and number of nucleotide polymorphisms in comparison.

The evolutionary history of the Leptotila genus was inferred by Johnson and Weckstein (Reference Johnson and Weckstein2011) using a molecular phylogeny to investigate the role of the Central American land-bridge in avian dispersal-driven diversification events. Their results show that this genus diverged into two clades at a time that coincides with the closure of the Isthmus of Panama forming the Central American land-bridge. Their study suggests that the northern clade dispersed southwards from North America stepwise through Central America, crossing the land-bridge into South America, while the southern clade, containing L. rufaxilla and L. cassinii, dispersed in the opposite direction from South America across the Panama Isthmus (Johnson and Weckstein Reference Johnson and Weckstein2011). Our results indicate that the Grenada Dove is placed within the southern clade and not with the northern clade containing L. jamaicensis. We hypothesise that there was a shared evolutionary history with an upwards dispersal event from South America by the ancestor shared by L. rufaxilla, L. cassinii and the Grenada Dove with a range expansion to Central America and the southern islands of the Lesser Antilles. Dispersal over the barriers presented by the South American Andes and the Caribbean Sea implies reduced gene flow between the ancestral populations and eventual isolation. This is supported by the absence of the indel in the Grenada Dove nDNA, suggesting this mutation occurred after divergence from the shared common ancestor with L. cassinii and L. plumbeiceps. Our data support the hypothesis that allopatric speciation due to isolation of the Grenada Dove and exposure to differing selective pressures resulted in the evolution of Grenada Dove as a distinct species.

Our analyses were able to predict six life history values for the Grenada Dove. Values for clutch size, clutches per year and fledging age agreed with documented values from focal observations (Rusk, unpubl. data). With our model estimations of only two eggs per clutch and 1.4 clutches per year the Grenada Dove has a relatively slow breeding rate, which has been shown to be correlated with extinction risk (Robinson et al. Reference Robinson, Hau, Klasing, Wikelski, Brawn, Austin, Tarwater and Ricklefs2010, Hutchings et al. Reference Hutchings, Myers, García, Lucifora and Kuparinen2013). Given that IUCN generated a 4.2 year generation length and with a low predicted longevity (genus median = 8.6 years), these values suggest this species would be slow to recover after a major population decline (Bird et al. Reference Bird, Martin, Akçakaya, Gilroy, Burfield, Garnett, Symes, Taylor, Şekercioğlu and Butchart2020). As not all traits are phylogenetically conserved (Kamilar and Cooper Reference Kamilar and Cooper2013), we tested for phylogenetic signal and found that signal strength varied amongst traits. Clutch size and generation length both exhibited strong phylogenetic signals, as is expected owing to their correlation with body mass, which indicates phylogenetic conservatism of these life history traits (Kamilar and Cooper Reference Kamilar and Cooper2013, Calhoon et al. Reference Calhoon, Jimenez, Harper, Jurkowitz and Williams2014). All other traits exhibited low phylogenetic signal strength, which is not uncommon in traits such as these that often evolve as a response to differing environmental conditions (Kamilar and Cooper Reference Kamilar and Cooper2013, Martin et al. Reference Martin, Ton and Oteyza2018). Although our model validation has suggested that many of our predicted values are accurate, these results have to be used cautiously, and would benefit from validation from field data (Ando Reference Ando2019, Bolton et al. Reference Bolton, Van Oosterhout, Collar and Bell2016).

When comparing predicted values for life history traits for the Grenada Dove and L. cassinii (the species we reveal shares the most recent common ancestor and the smallest genetic distance) and L. rufaxilla (the species it was once considered a subspecies of) we find that a clutch size of two and more than one clutch per year were shared between all three species. The Grenada Dove and L. cassinii both having a hatching weight of over 3 g and an incubation time of around 14 days whereas L. rufaxilla has a slightly lower hatching weight of around 2 g and a slightly higher incubation time of 15 days. While there is no great difference in life history traits between studied members of the Leptotila genus, as expected, we find that the species with life history traits most similar to the Grenada Dove is L. cassinii which is the species we also reveal has the smallest genetic distance.

The more favourable conservation status of the other Leptotila species examined in this study may be attributed to their much larger ranges and availability of suitable habitat (BirdLife International 2021). We reinforce the need for increasing suitable habitat available to the Grenada Dove to reduce extinction risk, both through habitat restoration and protection, as has been outlined in the recovery and action plan for this species, as well as identification of additional habitat for establishing new populations and enabling population increase (Rusk et al. Reference Rusk, Temple, Wege, Wallace, Blockstein, Ledec and Sorenson2008).

Supplementary Materials

To view supplementary material for this article, please visit http://doi.org/10.1017/S0959270922000065.

Acknowledgements

We would like to thank the Forest and National Parks Department of the Government of Grenada for offering their support and expertise, in particular Chief Forestry Officer Anthony Jeremiah for providing us with the eggshell sample and Forest Ranger Doland Francis for providing us with the feather samples. We thank Dr Nigel Collar for his advice and support with this manuscript, Dr Ian McDowall for his support and guidance with the molecular methods, Dr Alejandro González-Voyer for his advice and support with multivariate phylogenetic comparative methods and the Chester Biological Sciences technical team for their laboratory assistance. This work was funded by the University of Chester. Ethical approval was awarded by the Faculty Research and Ethics Committee at the University of Chester. CITES import (546012/01, 567389/01) and export (Grenada IACUC 14002, 14004) permits were obtained for movement of samples.

References

Akaike, H. (1987) Factor analysis and AIC. Psychometrika 52: 371386.CrossRefGoogle Scholar
Allio, R., Donega, S., Galtier, N. and Nabholz, B. (2017) Large variation in the ratio of mitochondrial to nuclear mutation rate across animals: Implications for genetic diversity and the use of mitochondrial DNA as a molecular marker. Mol. Biol. Evol. 34: 27622772.CrossRefGoogle ScholarPubMed
Ando, H. (2019) Genetic and ecological conservation issues for oceanic island birds, revealed by a combination of the latest molecular techniques and conventional field work. Ecol. Res. 34: 255264.CrossRefGoogle Scholar
Bastide, P., Solís-Lemus, C., Kriebel, R., Sparks, K. W. and Ané, C. (2018) Phylogenetic comparative methods on phylogenetic networks with reticulations. Syst. Biol. 67: 800820.CrossRefGoogle ScholarPubMed
Benson, D.A., Cavanaugh, M., Clark, K., Karsch-Mizrachi, I., Lipman, D. J., Ostell, J. and Sayers, E. W. (2013) GenBank. Nucleic Acids. Res. 41: 3642.CrossRefGoogle ScholarPubMed
Bird, J. P., Martin, R., Akçakaya, H. R., Gilroy, J., Burfield, I. J., Garnett, S., Symes, A., Taylor, J., Şekercioğlu, Ç. H. and Butchart, S. H. M. (2020) Generation lengths of the world’s birds and their implications for extinction risk. Conserv. Biol. 34: 12521261.CrossRefGoogle ScholarPubMed
BirdLife International (2018) Leptotila wellsi. The IUCN Red List of Threatened Species. https://www.iucnredlist.org/species/22690874/131031811 (Accessed 19 February 2020).Google Scholar
BirdLife International (2021) Data Zone Species search http://datazone.birdlife.org/species/search (Accessed 6 September 2021)Google Scholar
Blockstein, D. E. and Hardy, J. W. (1989) The Grenada dove (Leptotila wellsi) is a distinct species. Auk 106: 334340.Google Scholar
Bolton, N. M., Van Oosterhout, C., Collar, N. J. and Bell, D. J. (2016) Population constraints on the Grenada Dove Leptotila wellsi: preliminary findings and proposals from south-west Grenada. Bird Conserv. Internatn. 26: 205213.CrossRefGoogle Scholar
Bouckaert, R., Vaughan, T. G., Barido-Sottani, J., Duchêne, S., Fourment, M., Gavryushkina, A., Heled, J., Jones, G., Kühnert, D., De Maio, N., Matschiner, M., Mendes, F., Müller, N. F., Ogilvie, H. A., Du Plessis, L., Popinga, A., Rambaut, A., Siveroni, I., Suchard, M. A., Wu, C. H., Xie, D., Zhang, C., Stadler, T. and Drummond, A. J. (2019) BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis. PLOS Comput. Biol. 15: 128.CrossRefGoogle ScholarPubMed
Brooks, T. M., Collar, N. J., Green, R. E., Marsden, S. J. and Pain, D. J. (2008) The science of bird conservation. Bird Conserv. Internatn. 18: 212.CrossRefGoogle Scholar
Calhoon, E. A., Jimenez, A. G., Harper, J. M., Jurkowitz, M. S. and Williams, J. B. (2014) Linkages between mitochondrial lipids and life history in temperate and tropical bird. Physiol. Biochem. Zool. 87: 265275.CrossRefGoogle Scholar
Clayton, D. H., Bush, S. E., Goates, B. M. and Johnson, K. P. (2003) Host defense reinforces host–parasite cospeciation. PNAS 100: 1569415699.CrossRefGoogle ScholarPubMed
Cressler, C., Butler, M. and King, A. (2015) Detecting adaptive evolution in phylogenetic comparative analysis using the Ornstein-Uhlenbeck model. Syst. Biol. 64: 953968CrossRefGoogle ScholarPubMed
Cuervo, J. J. and Møller, A. P. (2019) Demographic, ecological, and life-history traits associated with bird population response to landscape fragmentation in Europe. Landsc. Ecol. 9: 113.Google Scholar
Devenish-Nelson, E. S., Weidemann, D., Townsend, J. and Nelson, H. P. (2019) Patterns in island endemic forest-dependent bird research: the Caribbean as a case-study. Biodivers. Conserv. 28: 18851904.CrossRefGoogle Scholar
Ericson, P. G. P., Anderson, C. L., Britton, T., Elzanowski, A., Johansson, U. S., Källersjö, M., Ohlson, J. I., Parsons, T. J., Zuccon, D. and Mayr, G. (2006) Diversification of Neoaves: Integration of molecular sequence data and fossils. Biol. Lett. 2: 543547.CrossRefGoogle ScholarPubMed
Frankham, R. (1998) Inbreeding and extinction: Island populations. Conserv. Biol. 12: 665675.CrossRefGoogle Scholar
Gene Codes Corporation (2017) Sequencher DNA sequence analysis software from Gene Codes Corporation. Gene Codes Corporation. https://www.genecodes.com/ (Accessed 19 February 2019)Google Scholar
González-del-Pliegoy, P., Freckleton, R.P., Edwards, D.P., Koo, M. S., Scheffers, B. R., Pyron, R. A. and Jetz, W. (2019) Phylogenetic and trait-based prediction of extinction risk. Curr. Biol. 29: 15571563CrossRefGoogle Scholar
Goolsby, E. W., Bruggeman, J. and Ané, C. (2017) Rphylopars: fast multivariate phylogenetic comparative methods for missing data and within-species variation. Meth. Ecol. Evol. 8: 2227.CrossRefGoogle Scholar
Hasegawa, M., Kishino, H. and Yano, T. (1985) Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J. Mol. Evol. 22: 160174.CrossRefGoogle ScholarPubMed
HBW Alive (2020) Handbook of the birds of the world alive | HBW Alive. Barcelona: Lynx Edicions.Google Scholar
Heled, J. and Drummond, A. J. (2010) Bayesian inference of species trees from multilocus data. Mol. Biol. Evol. 27: 570580.CrossRefGoogle ScholarPubMed
Horswill, C., Manica, A., Daunt, F., Newell, M., Wanless, S., Wood, M., and Matthiopoulos, J. (2021) Improving assessments of data-limited populations using life-history theory. J. Appl. Ecol. 58: 12251236.CrossRefGoogle Scholar
Hurst, G. D. D. and Jiggins, F. M. (2005) Problems with mitochondrial DNA as a marker in population, phylogeographic and phylogenetic studies: The effects of inherited symbionts. Proc. R. Soc. B. 272: 15251534.CrossRefGoogle ScholarPubMed
Hutchings, J. A., Myers, R. A., García, VB, Lucifora, L. O. and Kuparinen, A. (2013) Life-history correlates of extinction risk and recovery potential. Ecol. Appl. 23: 515522.Google Scholar
James, T. D., Salguero-Gómez, R., Jones, O. R., Childs, D. Z. and Beckerman, A. P. (2021) Bridging gaps in demographic analysis with phylogenetic imputation. Conserv. Biol. 35: 12101221.CrossRefGoogle ScholarPubMed
Johnson, K. P., Adler, F. R. and Cherry, J. L. (2007) Genetic and phylogenetic consequences of island biogeography. Evolution 54: 387396.Google Scholar
Johnson, K. P. and Clayton, D. H. (2000a) A molecular phylogeny of the dove genus Zenaida: mitochondrial and nuclear DNA sequences. Condor 102: 864870.CrossRefGoogle Scholar
Johnson, K. P. and Clayton, D. H. (2000b) Nuclear and mitochondrial genes contain similar phylogenetic signal for pigeons and doves (Aves: Columbiformes). Mol. Phylogenet. Evol. 14: 141151.CrossRefGoogle ScholarPubMed
Johnson, K. P. and Weckstein, J. D. (2011) The Central American land bridge as an engine of diversification in New World doves. J. Biogeogr. 38: 10691076.CrossRefGoogle Scholar
Johnson, K. P., De Kort, S., Dinwoodey, K., Dinwoodey, K., Mateman, A. C., Cate, C. T., Lessells, C. M. and Clayton, D. H. (2001) A molecular phylogeny of the dove genera Streptopelia and Columba. The Auk 118: 874887.CrossRefGoogle Scholar
Kamilar, J. M. and Cooper, N. (2013) Phylogenetic signal in primate behaviour, ecology and life history. Philos. T. Roy. Soc. B. 368: 110.CrossRefGoogle ScholarPubMed
Kerr, K. C. R., Lijtmaer, D. A., Barreira, A. S., Hebert, P. D. N. and Tubaro, P. L. (2009) Probing evolutionary patterns in neotropical birds through DNA barcodes. PLoS ONE 4: 4379.CrossRefGoogle ScholarPubMed
Kerr, K. C. R., Stoeckle, M. Y., Dove, C. J., Weigt, L. A., Francis, C. M. and Hebert, P. D. N. (2006) Zenaida asiatica voucher FMNH 430791 cytochrome oxidase subunit 1 (COI) gene, partial cds; mitochondrial GenBank: DQ433271.1. NCBI. https://www.ncbi.nlm.nih.gov/nuccore/DQ433271.1 (Accessed 23 December 2019).Google Scholar
Kerr, K. C. R., Stoeckle, M. Y., Dove, C. J., Weigt, L. A., Francis, C. M. and Hebert, P. D. N. (2007) Comprehensive DNA barcode coverage of North American birds. Mol. Ecol. Notes 7: 535543.CrossRefGoogle ScholarPubMed
Krajewski, C., Sipiorski, J. T. and Anderson, F. E. (2010) Complete mitochondrial genome sequences and the phylogeny of cranes (Gruiformes: Gruidae). The Auk 127: 440452.CrossRefGoogle Scholar
Ksepka, D. T., Parham, J. F., Allman, J.F., Benton, M. J., Carrano, M. T., Cranston, K. A., Donoghue, P. C. J., Head, J. J., Hermsen, E. J., Irmis, R. B., Joyce, W. G., Kohli, M., Lamm, K. D., Leehr, D., Patané, J. L., Polly, P. D., Phillips, M. J., Smith, N. A., Smith, N. D., Van Tuinen, M., Ware, J. L. and Warnock, R. C. M. (2015) The fossil calibration database - A new resource for divergence dating. Syst. Biol. 64: 853859.CrossRefGoogle ScholarPubMed
Kumar, S., Stecher, G., Li, M., Knyaz, C. and Tamura, K. (2018) MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35: 15471549.CrossRefGoogle ScholarPubMed
Kumar, S., Stecher, G., Suleski, M. and Hedges, S. B. (2017) TimeTree: A resource for timelines, timetrees, and divergence times. Mol. Biol. Evol. 34: 3978.CrossRefGoogle ScholarPubMed
Latta, S.C. (2012) Avian research in the Caribbean: Past contributions and current priorities. J. Field Ornithol. 83: 107121.CrossRefGoogle Scholar
Li, C. H., Shi, W. and Shi, W. Y. (2015) Mitochondrial genome sequence of Egyptian swift Rock Pigeon (Columba livia breed Egyptian swift). Mitochondr. DNA 26: 479480.CrossRefGoogle ScholarPubMed
Lockwood, J. L., Brooks, T. M. and Mckinney, M. L. (2000) Taxonomic homogenization of the global avifauna. Anim. Conserv. 3: 2735.CrossRefGoogle Scholar
Mace, G. M., Collar, N. J., Gaston, K. J., Hilton-Taylor, C., Akçakaya, H. R., Leader-Williams, N., Milner-Gulland, E. J. and Stuart, S. N. (2008) Quantification of extinction risk: IUCN’s system for classifying threatened species. Conserv. Biol. 22: 14241442.CrossRefGoogle ScholarPubMed
Martin, T. E., Ton, R. and Oteyza, J. C. (2018) Adaptive influence of extrinsic and intrinsic factors on variation of incubation periods among tropical and temperate passerines. The Auk 135: 101113.CrossRefGoogle Scholar
Mello, B. (2018) Estimating timetrees with MEGA and the timetree resource. Mol. Biol. Evol. 35: 23342342.CrossRefGoogle ScholarPubMed
Miller, M. J., Bermingham, E. and Winker, K. (2008) Leptotila cassini voucher UAM 20256 NADH dehydrogenase subunit 2 (ND2) gene, complete cds; mitochondrial GenBank: FJ175707.1. NCBI. https://www.ncbi.nlm.nih.gov/nuccore/FJ175707.1 (Accessed 22 December 2019).Google Scholar
Monceau, K., Cézilly, F., Moreau, J., Motreuil, S. and Wattier, R. (2013) Colonisation and diversification of the Zenaida Dove (Zenaida aurita) in the Antilles: Phylogeography, contemporary gene flow and morphological divergence. PLoS ONE 8: 82189.CrossRefGoogle ScholarPubMed
Monceau, K., Wattier, R. A., Moreau, J., Motreuil, S. and Cezilly, F. (2012) Zenaida galapagoensis isolate galW10 cytochrome oxidase sununit I (COI) gene, partial cds; mitochondrial GenBank: JQ420133.1. NCBI. https://www.ncbi.nlm.nih.gov/nuccore/JQ420133.1 (Accessed 23 December 2019).Google Scholar
Myhrvold, N. P., Baldridge, E., Chan, B., Sivam, D., Freeman, D. L. and Ernest, S. K. M. (2015) An amniote life-history database to perform comparative analyses with birds, mammals, and reptiles. Ecology 96: 3109.CrossRefGoogle Scholar
Nei, M. and Kumar, S. (2000) Molecular evolution and phylogenetics. Oxford: Oxford University Press.Google Scholar
Pacheco, M. A., Battistuzzi, F. U., Lentino, M., Aguilar, R. F., Kumar, S. and Escalante, A. A. (2011) Evolution of modern birds revealed by mitogenomics: Timing the radiation and origin of major orders. Mol. Biol. Evol. 28:1 9271942.CrossRefGoogle ScholarPubMed
Parlato, E. H., Armstrong, D. P. and Innes, J. G. (2015) Traits influencing range contraction in New Zealand’s endemic forest birds. Oecologia 179: 319328.CrossRefGoogle ScholarPubMed
Patel, S., Waugh, J., Millar, C. D. and Lambert, D. M. (2010) Conserved primers for DNA barcoding historical and modern samples from New Zealand and Antarctic birds. Mol. Ecol. Resour. 10: 431438.CrossRefGoogle ScholarPubMed
Pearson, R. G., Stanton, J. C., Shoemaker, K. T., Aiello-Lammens, M. E., Ersts, P. J., Horning, N., Fordham, D. A., Raxworthy, C. J., Ryu, H. Y., Mcnees, J. and Akçakaya, H. R. (2014) Life history and spatial traits predict extinction risk due to climate change. Nat. Clim. Chang. 4: 217221.CrossRefGoogle Scholar
Penone, C., Davidson, A. D., Shoemaker, K. T., Marco, M. D., Brooks, T. M., Young, B. E., Graham, C. H., Costa, G. C., Di Marco, M., Rondinini, C., Brooks, T. M., Young, B. E., Graham, C. H. and Costa, G. C. (2014) Imputation of missing data in life-history trait datasets: Which approach performs the best? Meth. Ecol. Evol. 5: 961970.CrossRefGoogle Scholar
Pereira, S. L., Johnson, K. P., Clayton, D. H. and Baker, A. J (2007) Mitochondrial and nuclear DNA sequences support a cretaceous origin of Columbiformes and a dispersal-driven radiation in the paleogene. Syst. Biol. 56: 656672.CrossRefGoogle Scholar
Peters, C., Nelson, H., Rusk, B. L. and Muir, A. (2019) A novel method to optimise the utility of underused moulted plumulaceous feather samples for genetic analysis in bird conservation. Conserv. Genet. Resour. 12: 457467CrossRefGoogle Scholar
Peters, S. E. and McClennen, M. (2015) The Paleobiology Database application programming interface. Paleobiology 42: 17.CrossRefGoogle Scholar
Prychitko, T. M. and Moore, W. S. (2003) Alignment and phylogenetic analysis of β-fibrinogen intron 7 sequences among avian orders reveal conserved regions within the intron. Mol. Biol. Evol. 20: 762771.CrossRefGoogle ScholarPubMed
Qu, J. (2018) Streptopelia decaocto mitochondrion, complete genome NCBI Reference Sequence: NC_037513.1. NCBI. https://www.ncbi.nlm.nih.gov/nuccore/NC_037513.1 (Accessed 22 December 2019).Google Scholar
R Core Team (2020) R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing.Google Scholar
Riek, A. and Bruggeman, J. (2013) Estimating field metabolic rates for Australian marsupials using phylogeny. Comp. Biochem. Physiol. A. Mol. Integr. Physiol. 164: 598604.CrossRefGoogle ScholarPubMed
Rivera-Milán, F. F., Bertuol, P., Simal, F. and Rusk, B. L. (2015) Distance sampling survey and abundance estimation of the critically endangered. Condor 117: 8793.CrossRefGoogle Scholar
Robinson, W. D., Hau, M., Klasing, K. C., Wikelski, M., Brawn, J. D., Austin, S. H., Tarwater, C. E. and Ricklefs, R. E. (2010) Diversification of life histories in New World birds. The Auk 127: 253262.CrossRefGoogle Scholar
Rubinoff, D. and Holland, B. S. (2005) Between two extremes: mitochondrial DNA is neither the panacea nor the nemesis of phylogenetic and taxonomic inference. Syst. Biol. 54: 952961.CrossRefGoogle ScholarPubMed
Rusk, B. L. (2017) Long-term population monitoring of the Critically Endangered Grenada Dove (Leptotila wellsi) on Grenada, West Indies. J. Caribb. Ornithol. 30: 4956.Google Scholar
Rusk, B. L., Temple, S. A., Wege, D, Wallace, G. Blockstein, D. E., Ledec, G. and Sorenson, L. (2008) Recovery and action plan for the Grenada Dove. Unpublished draft. St. George’s, Grenada: Forestry and National Parks Department, Government of Grenada.Google Scholar
Santiago-Alarcon, D., Rodríguez-Ferraro, A., Parker, P. G. and Ricklefs, R. E. (2014) Different meal, same flavor: Cospeciation and host switching of haemosporidian parasites in some non-passerine birds. Parasites Vector 7: 286.CrossRefGoogle ScholarPubMed
Schindel, D. E., Stoeckle, M. Y., Milensky, C., Trizna, M., Schmidt, B., Gebhard, C. and Graves, G. (2011) Project description: DNA barcodes of bird species in the National Museum of Natural History, Smithsonian Institution, USA. Zookeys 152: 8792.CrossRefGoogle Scholar
Schwarz, G. (1978) Estimating the dimension of a model. Ann. Stat. 6: 461464.CrossRefGoogle Scholar
Shapiro, B., Sibthorpe, D., Rambaut, A., Austin, J., Wragg, G. M., Bininda-Emonds, O. R. P., Lee, P. L. M. and Cooper, A. (2002) Flight of the Dodo. Science 295: 1683.CrossRefGoogle ScholarPubMed
Soares, A. E. R, Novak, B. J., Haile, J., Heupink, T. H., Fjeldså, J., Gilbert, M. T. P., Poinar, H., Church, G. M. and Shapiro, B. (2016) Complete mitochondrial genomes of living and extinct pigeons revise the timing of the columbiform radiation. BMC Evol. Biol. 16: 19.CrossRefGoogle ScholarPubMed
Sorenson, M. D., Oneal, E., García-Moreno, J. and Mindell, D. P. (2003) More taxa, more characters: The hoatzin problem is still unresolved. Mol. Biol. Evol. 20: 14841498.CrossRefGoogle ScholarPubMed
Stoeckle, M. (2011) Streptopelia roseogrisea voucher CWSn/a003 cytochrome oxidase subunit 1 (COI) gene, partial cds; mitochondrial GenBank: JN801382.1. NCBI. https://www.ncbi.nlm.nih.gov/nuccore/JN801382.1 (Accessed 26 December 2019).Google Scholar
Storchová, L. and Hořák, D. (2018) Life-history characteristics of European birds. Glob. Ecol. Biogeogr. 27: 400406.CrossRefGoogle Scholar
Sweet, A. D. and Johnson, K. P. (2015) Patterns of diversification in small New World ground doves are consistent with major geologic events. The Auk 132: 300312.CrossRefGoogle Scholar
Swenson, N. G. (2014) Phylogenetic imputation of plant functional trait databases. Ecography 37: 105110.CrossRefGoogle Scholar
Tamura, K., Battistuzzi, F. U., Billing-Ross, P., Murillo, O., Filipski, A. and Kumar, S. (2012) Estimating divergence times in large molecular phylogenies. PNAS 109: 1933319338.CrossRefGoogle ScholarPubMed
Tamura, K., Tao, Q. and Kumar, S. (2018) Theoretical foundation of the reltime method for estimating divergence times from variable evolutionary rates. Mol. Biol. Evol. 35: 17701782.CrossRefGoogle ScholarPubMed
Tavaré, S. (1986) Some probabilistic and statistical problems in the analysis of DNA sequences. Lect. Math. Life Sci. 17: 5786.Google Scholar
Tavares, E. S., Gonçalves, P., Miyaki, C. Y. and Baker, A. J. (2011) DNA barcode detects high genetic structure within neotropical bird species. PLoS ONE 6: 110.CrossRefGoogle ScholarPubMed
Tershy, B. R., Shen, K. W., Newton, K. M., Holmes, N. D. and Croll, D. A. (2015) The importance of islands for the protection of biological and linguistic diversity. Bioscience 65: 592597.CrossRefGoogle Scholar
Untergasser, A., Cutcutache, I., Koressaar, T., Ye, J., Faircloth, B. C., Remm, M. and Rozen, S. G. (2012) Primer3-new capabilities and interfaces. Nucleic Acids Res. 40: 112.CrossRefGoogle ScholarPubMed
Weber, D., Stewart, B. S., Garza, J. C. and Lehman, N. (2000) An empirical genetic assessment of the severity of the northern elephant seal population bottleneck. Curr. Biol. 10: 12871290.CrossRefGoogle ScholarPubMed
Wood, K. A., Stillman, R. A. and Hilton, G. M. (2018) Conservation in a changing world needs predictive models. Anim. Conserv. 21: 8788.CrossRefGoogle Scholar
Figure 0

Table 1. List of species and their Genbank accession number for all sequences used in phylogenetic analysis.

Figure 1

Figure 1. Phylogenetic analysis for 24 Neotropical columbid species and outgroup using the Maximum Likelihood method. Grenada Dove sequences indicated by the circular symbol; a. all mitochondrial and nuclear sequence data estimated using the Hasegawa-Kishino-Yano model (+G, parameter = 0.2281); b. mtDNA using General Time Reversible model ([+I], 38.01% sites); c. β-fibrinogen using the Tamura 3-parameter model. The trees are drawn to scale, with branch lengths measured in the number of substitutions per site and nodes supported by bootstrap values of >80%.

Figure 2

Figure 2. Time-calibrated phylogenies for Leptotila species and outgroup using the RelTime-ML method a. all mitochondrial and nuclear sequence data estimated using the Hasegawa-Kishino-Yano model (+G, parameter = 0.1547), b. mtDNA General Time Reversible model ([+I], 35.97% sites) c. β-fibrinogen using the Tamura 3-parameter model. Number scale indicates million years ago (mya). Nodes indicate 7 estimated divergence time calibration points of estimated divergence time ranges in million years ago (mya) obtained using the TimeTree resource (Kumar et al.2017).

Figure 3

Table 2. The number of substitutions per site (p-distance) and the number of base differences (polymorphic sites) between sequences for a total of 3,346 bp positions in the final dataset. *2,720 bp positions in the final dataset.

Figure 4

Figure 3. Phylogenetic analysis for 25 Neotropical columbid species and outgroup using the Maximum Likelihood method. Grenada Dove sequences indicated by the circular symbol and Leptotila jamaicensis are indicated by a square using the Hasegawa-Kishino-Yano model (+G, parameter = 0.2496). Tree drawn to scale, with branch lengths measured in the number of substitutions per site and nodes supported by bootstrap values of >80%.

Figure 5

Figure 4. Plot displaying the squared difference between the eb model predicted value and the true value divided by the trait median using the genus Columbina along with Normalized Root Mean Squared Error [NRMSE (Penone et al.2014, James et al.2021)]. Raw data are represented by the ° and the first, second and third quartile are represented by the box plot. Mean Squared Error per life history trait.

Figure 6

Figure 5. Predicted trait values with upper and lower 95% confidence intervals for the Grenada Dove Leptotila wellsi in comparison with other Leptotila species for the following lift history traits using the EB model: a. clutch size, b. clutches per year, c. incubation time, d. hatching weight; and median estimates: e. fledging age, f. longevity. Species of interest highlighted red. Known trait values represented by individual grey point.

Supplementary material: File

Peters et al. supplementary material

Peters et al. supplementary material 1

Download Peters et al. supplementary material(File)
File 4.6 KB
Supplementary material: File

Peters et al. supplementary material

Peters et al. supplementary material 2

Download Peters et al. supplementary material(File)
File 213.6 KB
Supplementary material: File

Peters et al. supplementary material

Peters et al. supplementary material 3

Download Peters et al. supplementary material(File)
File 2.3 KB