Evolutionary history of a vanishing radiation: isolation-dependent persistence and diversification in Pacific Island partulid tree snails
© Lee et al.; licensee BioMed Central Ltd. 2014
Received: 16 December 2013
Accepted: 16 September 2014
Published: 24 September 2014
Partulid tree snails are endemic to Pacific high islands and have experienced extraordinary rates of extinction in recent decades. Although they collectively range across a 10,000 km swath of Oceania, half of the family’s total species diversity is endemic to a single Eastern Pacific hot spot archipelago (the Society Islands) and all three partulid genera display highly distinctive distributions. Our goal was to investigate broad scale (range wide) and fine scale (within-Society Islands) molecular phylogenetic relationships of the two widespread genera, Partula and Samoana. What can such data tell us regarding the genesis of such divergent generic distribution patterns, and nominal species diversity levels across Oceania?
Museum, captive (zoo) and contemporary field specimens enabled us to genotype 54 of the ~120 recognized species, including many extinct or extirpated taxa, from 14 archipelagoes. The genera Partula and Samoana are products of very distinct diversification processes. Originating at the western edge of the familial range, the derived genus Samoana is a relatively recent arrival in the far eastern archipelagoes (Society, Austral, Marquesas) where it exhibits a stepping-stone phylogenetic pattern and has proven adept at both intra-and inter- archipelago colonization. The pronounced east—west geographic disjunction exhibited by the genus Partula stems from a much older long-distance dispersal event and its high taxonomic diversity in the Society Islands is a product of a long history of within-archipelago diversification.
The central importance of isolation for partulid lineage persistence and diversification is evident in time-calibrated phylogenetic trees that show that remote archipelagoes least impacted by continental biotas bear the oldest clades and/or the most speciose radiations. In contemporary Oceania, that isolation is being progressively undermined and these tree snails are now directly exposed to introduced continental predators throughout the family’s range. Persistence of partulids in the wild will require proactive exclusion of alien predators in at least some designated refuge islands.
The Pacific Ocean comprises a third of the Earth’s surface and contains approximately 25,000 oceanic islands . Most are clustered within hot spot and island arc volcanic archipelagoes scattered across Oceania, a tropical region extending from New Guinea to Hawaii and Easter Island . These islands represent highly distinctive evolutionary settings for terrestrial biotas: they have never been connected to continental landmasses and trans-oceanic founder speciation, followed by varying degrees of in situ cladogenesis, is the primary diversification process -.
Oceanic island biotas start out with a depauperate and disharmonic composition, i.e., suites of taxa that are prevalent in continental ecosystems, but incapable of dispersing across oceanic barriers, are missing . This facilitates the survival of relict lineages and the evolution of endemic adaptive radiations ,,, resulting in a “diversity and stability paradox” : oceanic island biotas are often species poor, but have high endemism; they appear stable (pending maintenance of biotic isolation), but are susceptible to collapse following anthropogenic disturbance. Indeed, anthropogenic habitat modification/disturbance and human-introduced continental species are the primary drivers of oceanic island extinctions . The small ranges/populations of oceanic island endemics, together with their relative lack of defensive or competitive abilities, render them exceptionally vulnerable to introduced aliens ,,-.
Although individual islands are often species poor, the cumulative effect of the “diversity and stability paradox”  is of global significance: most historical extinctions of birds , plants  and mollusks  worldwide have involved oceanic island endemics. Of particular concern in recent decades has been a major wave of Pacific Island land snail extinctions originating in a misguided biological control program that deliberately introduced a number of alien carnivorous land snails, most notably Euglandina rosea, across much of Oceania -. It is estimated that E. rosea alone has probably caused the extinction of 134 Pacific Island land snail endemics  and among the hardest impacted are members of the family Partulidae -, the subject of this study.
Evolutionary studies of oceanic island lineages share common goals such as the identification of continental sister taxa and the inference of the number, timing and directionality of colonization events among archipelagoes and islands . Unlike Galapagos finches  and Hawaiian silverswords , for instance, a convincing sister group to Partulidae has not yet been identified ; the poorly-studied New Caledonia endemic family Draparnaudiidae (6 known species) is considered the likeliest candidate ,. All available data suggest that Partulidae have had a long evolutionary history in Oceania, dispersing to new islands as older ones sink below sea level. Pilsbry and Cooke  partitioned Partulidae into 3 genera (Partula, Samoana and Eua), based primarily on male genital characters, and these have been substantially corroborated by independent molecular character sets ,,-, with the exception of one Tahitian species complex with highly variable genitalia .
Eua has the simplest genital anatomy, contains 4 species, and is restricted to the central part of the family’s range: 1 species on `Eua (Tonga), 3 in the Samoan Archipelago (Figure 1). Samoana contains approximately 23 species and has a disjunct distribution: extensively distributed throughout central and eastern archipelagoes but with a solitary western species, S. fragilis, in the Mariana Islands 5000 km from its nearest congener in Fiji (Figure 1). Partula, the most speciose genus (~92 species), exhibits a strikingly different disjunct pattern: widespread in western archipelagoes, absent from the central part of the family’s range and, excluding anthropogenic populations , restricted to a small subset of eastern islands: 1 species in Rarotonga (Cook Islands), and a remarkable 59 species ,, in a single hot spot archipelago, the Society Islands (Figure 1).
Kondo and Burch  implicitly proposed the following phylogenetic pattern of generic relationships: (Eua (Samoana, Partula)), a topology that has been corroborated by available molecular phylogenies. Land snail gene trees that have included all three genera  have placed Eua basally, a topology that positions the last common ancestor of extant Partulidae in the central part of the present-day partulid range. Partulid phylogenies rooted with Eua species ,, have consistently recovered reciprocally monophyletic Samoana and Partula genera. Goodacre and Wade’s  nuclear ribosomal phylogeny supported an eastward expansion of Samoana species from Samoa to the Society and Marquesan archipelagoes, but the Partula clade had a pronounced basal polytomy and yielded little historical insight into this genus’s disjunct distribution.
A comprehensive understanding of partulid diversification processes requires not only a range-wide perspective, but also the development of finer scale within-archipelago and within-island phylogenies. This has been best developed for the Society Island radiation where a number of workers have proposed a “progression rule” model of evolution – lineages sequentially colonize newer islands within an archipelago as they are formed  - for the well-studied Partula species of Moorea and Tahiti ,,-. In this model, neighboring older islands (Leeward group) seeded the younger islands (Windward group), first Moorea, then Tahiti, via single colonization events and subsequent speciation occurred in situ within each island ,,-. This model does not apply to the four endemic Society Island Samoana species, two of which have multi-island distributions ,,.
Summary data for partulid taxa genotyped for this study showing the fraction of species in each genus sampled, the number of islands and archipelagoes they were sourced from and the intensive sampling of Tahitian and Moorean Partula populations
Species sampled (species described)
Taxonomic designations, sampling locations, IUCN Red List Status and sources for the 54 partulid species genotyped in this study
IUCN Red list status
Eua globosa Pilsbry & Cooke, 1934
E. expansa (Pease, 1872)
E. montana (Cooke & Crampton, 1930)
E. zebrina (Gould, 1847)
Tutuila, Samoa (Am.)
Samoana fragilis (Ferussac, 1821)
S. abbreviata (Mousson, 1869)
Tutuila, Am. Samoa
S. conica (Gould, 1847)
Tutuila, Am. Samoa
S. thurstoni (Cooke & Crampton, 1930)
Ofu, Am. Samoa
S. canalis (Mousson, 1865)
S. stevensoniana (Pilsbry, 1909)
S. margaritae (Crampton & Cooke, 1953)
S. oreas (Crampton & Cooke, 1953)
S. attenuata (Pease, 1864)
S. burchi Kondo, 1973
S. diaphana Crampton & Cooke, 1953
S. bellula (Hartman, 1885)
Ua Pou, Marquesas
S. decussatula (Pfeiffer, 1850)
Tahuata/Hiva Oa, Marquesas
S. strigata (Pease, 1868)
Nuku Hiva, Marquesas
Partula calypso Semper, 1865
P. thetis Semper, 1865
Ulong /Babeldaob/Ngeruktabel, Palau
FMNH/1995, 1998, 2006
P. gibba Ferussac, 1821
P. radiolata Pfeiffer, 1846
P. emersoni Pilsbry, 1913
P. carteriensis (Quoy & Gaimard, 1832)
New Britain, Bismarcks
P. similaris Hartman, 1886
P. auraniana Hartman, 1888
P. turneri Pfeiffer, 1860
P. lirata Mousson, 1865
P. subgonochila Mousson, 1871
Alofi, Wallis & Futuna
P. assimilis Pease, 1868
P. lutea Lesson, 1831
Bora Bora, Society
P. dentifera Pfeiffer, 1853
Extinct in the Wild
P. faba (Gmelin, 1791)
Extinct in the Wild
P. hebe (Pfeiffer, 1846)
Extinct in the Wild
P. labrusca Crampton & Cooke, 1953
P. meyeri Burch, 2007
P. tristis Crampton & Cooke, 1953
Extinct in the Wild
P. turgida (Pease, 1865)
P. arguta (Pease, 1865)
P. rosea Broderip, 1832
Extinct in the Wild
P. varia Broderip, 1832
Extinct in the Wild
P. aurantia Crampton, 1932
P. exigua Crampton, 1917
P. mirabilis Crampton, 1924
Extinct in the Wild
P. mooreana Hartman, 1880
Extinct in the Wild
P. suturalis Pfeiffer, 1855
Extinct in the Wild
P. taeniata (Mörch, 1850)
P. tohiveana Crampton, 1924
Extinct in the Wild
P. affinis Pease, 1867
P. clara Pease, 1864
P. filosa Pfeiffer, 1853
P. hyalina Broderip, 1832
Mangaia/ Mauke, Cooks
Fontaine & Gargominy/2002-4
P. nodosa Pfeiffer, 1851
Extinct in the Wild
P. otaheitana (Brugière, 1789)
P. producta Pease, 1865
Results & discussion
The unexpected phylogenetic placement of the two Palauan Partula taxa with Samoana species raises a question concerning their generic status. Pilsbry and Cooke  distinguished the genera Partula and Samoana primarily on male genital characters: the latter being differentiated into distinct epiphallus and penis segments connected by an extension of the penial retractor muscle that is absent in other partulids. Kondo’s Ph.D. thesis  contains detailed diagrams of the reproductive tracts of the three endemic Palauan species (P. thetis, P. calypso and P. leucothoe, the latter was unavailable for this study) and all three unambiguously displayed the male genital anatomy characteristic of the genus Partula. These anatomical data, together with the robust, well-resolved mt gene tree topology (Figure 3) demonstrate that the genus Partula is paraphyletic and that the distinct genital anatomy of Samoana species is derived. In contrast to Kondo and Burch’s  hypothesis of intergeneric relationships - (Eua (Samoana, Partula)) - our data are consistent with (Eua (Partula1 (Partula2, Samoana))) where Partula1 represents all non-Paluan members of the genus and Partula2 represents the two genotyped Palauan species.
West–east dispersal of Samoana
The most remarkable feature of the distribution of Samoana is the presence of a single isolated species, S. fragilis, bearing the characteristic Samoana genital anatomy , in the Western Pacific (Figure 1). Genotyping S. fragilis museum specimens sampled in 1945 allowed us to put this geographic disjunction into a phylogenetic context for the first time. Our two molecular markers yielded very different levels of phylogenetic resolution. The more conserved nuclear large ribosomal gene fragment produced a shallow clade with a basal polytomy such that interrelationships among Western (including the two Palauan Partula species) and Eastern Samoana taxa were unresolved (Figure 2). The mt marker provided much better resolution (Figure 3) and we base our discussion of the evolutionary history of Samoana on this topology. It shows that the Samoana clade is anchored in the West by its exclusive sister relationship with Palauan Partula species and by the placement of the Marianas species S. fragilis at its base - the sole western species is not a derived founder lineage (Figure 3). This is a somewhat surprising result because it places the inferred evolutionary origins of the genus not in the Eastern Pacific, home to 21/22 Samoana species (Figure 1), but instead in the far western portion of the partulid range (Figure 3).
The mt gene tree topology (Figure 3) is consistent with a western origin of the genus Samoana followed by progressive eastward colonization of derived lineages. S. fragilis is sister to the Samoan congener S. thurstoni, although the posterior probability value for this node is low. A derived clade, that approximates a stepping-stone dispersal pattern from West to East, is evident for the easternmost archipelagoes with Samoa as the regional source. Three closely-related Samoan species – S. conica, S. stevensonia and S. abbreviata – are robustly sister to derived congers in the Society, Austral and Marquesan archipelagoes, the latter two forming monophyletic clades. Reconstructing the multi-archipelagic dispersal pathways of these Eastern lineages is complicated by a basal polytomy involving Austral taxa and two Society Island lineages (Figure 3). However, the sister status of the Marquesan taxa with two Society Island congeners, S. diaphana and S. burchi, is consistent with an earlier inference  of dispersal and colonization among these two archipelagoes. Our results corroborate a previous allozyme study  regarding Marquesan monophyly and adaptive radiation – Marquesan taxa exhibit a diversity of distinct shell phenotypes. Interestingly, so do the two Society island montane forest sister taxa: the thin-shelled S. diaphana (Moorea and Thaiti) and the thick-shelled Tahitian endemic S. burchi.
Although our data provide new insights into the evolution and biogeography of the genus Samoana, they raise obvious follow-on questions: if the genus originated in the West, why does almost all the extant diversity occur in Eastern archipelagoes and why are Samoana species missing from a huge swathe of the family’s range from the Marianas to Fiji (Figure 1)? We consider it unlikely that the genus extended its range across this ~5000 km gap in a single dispersal event and hypothesize that at least some intermediate archipelagoes such as the Caroline Islands, Bismarck Archipelago, Solomon Islands, Santa Cruz Islands and Vanuatu may have supported populations at some time in the past. Present day absence of the genus Samoana from at least some of these archipelagoes is probably not due to a lack of dispersal ability. Our data indicate that a Samoan lineage was capable of stepping-stone colonization of Eastern archipelagoes up to 3000 km distant, a range that ostensibly puts the Western archipelagoes of Vanuatu, Santa Cruz and the Solomon Islands within equivalent reach. The absence of the genus Samoana from much of the Western range of the family may stem from (as yet unidentified) ecological factors that have driven original populations to extinction and prevented rare long-distant migrants from reestablishing new populations.
West–east Partula disjunction – product of a single dispersal event
The ~2000 km geographic disjunction – from the Lau Islands of Fiji to Rarotonga in the Cook Islands – separating Western and Eastern members of the genus Partula (Figure 1) has yet to be put into a comprehensive phylogenetic framework. Although our sampling of the genus Partula is heavily weighted toward Eastern populations, we have genotyped 11 species from across the Western range (Marianas, Palau, Carolines, Bismarck Archipelago and the Massim Region of Papua New Guinea, Vanuatu, Wallis and Futuna, Fiji). This broad sampling allows some confidence that we have captured much of the regional phylogenetic framework, if not all of its detail. Similarly, for the Eastern Partula species, our taxon sampling is geographically comprehensive: the 26 species genotyped include the sole non-anthropogenic Cook Island species, P. assimilis, as well as representatives from the all of the salient Society Islands except for Tahaa, which shares the same lagoon as Raiatea.
Phylogenetic relationships among the non-Palauan Western Partula taxa were obscured by a basal polytomy (Figure 4). Within-archipelago monophyly was evident only for Marianas taxa and these were weakly sister to the Micronesian P. emersoni. The phylogenetic affinities of Papua New Guinean populations of P. carteriensis and P. similaris, flagged as likely anthropogenic introductions , remained uncertain. The two easternmost taxa had distinct topological placements: P. subgonochila (Wallis & Futuna) was a member of the basal polytomy whereas the Fijian P. lirata was sister to one of the Vanuatu species, P. auraniana (Torres Islands). These disparate results imply that a more comprehensive sampling of Western taxa is required to capture a high-resolution understanding of regional evolutionary history. Unfortunately, this will have to be accomplished quickly given the ongoing decline and extirpation of many Western Pacific partulid populations and species ,-.
Eastern Partula diversification
The only Eastern Partula species not endemic to the Society Islands – P. assimilis of Rarotonga (Cook Islands) – is positioned within the Society Island clade (Figure 5), weakly sister to a minor, divergent Moorean P. suturalis clade and to the Bora Bora endemic P. lutea, with which it shares a large nuclear ribosomal genotype (Figure 2) as well as allozymic affinities . We agree with Johnson et al.  that P. assimilis probably represents a founder lineage from a Society Islands source; the only know Partula transplant from that hotspot archipelago apart from anthropogenic populations of the Tahitian endemic P. hyalina.
Concerning the Society Islands topology, a number of important stem nodes within the clade were poorly supported and, consistent with parallel studies on this archipelago’s biota , there was little agreement with progression rule  expectations for the archipelago, i.e., topological congruence with the chronological sequence of island formation. Only one (Huahine) of three older Society Islands represented in the topology (Bora Bora, Raiatea, Huahine; all members of the Leeward Island subgroup) positioned basally. The oldest Society Islands with partulids, Bora Bora and Raiatea, placed in derived positions, weakly sister to a heterogeneous assemblage that included the Rarotonga species, two Moorean clades and one Tahitian clade (Figure 5).
One consistent topological difference between the Leeward and Windward Island lineages concerned their degree of within-island monophyly. Genotyped snails from each Leeward Island formed well-supported (PPS ≥97) single island-specific clades. Snails from each Windward Island formed four discrete lineages, three of which had exclusive phylogenetic relationships (PPS ≥96) with snails in the other Windward island, as described in detail in an earlier study . However, the phylogenetic affinities of the largest and most deeply branched of the four Tahitian mt clades remained enigmatic. This robust (PPS = 99) Tahitian clade positioned on a basal polytomy that encompassed both Leeward and Windward Island lineages and it lacked an identifiable sister lineage, either on Moorea, or on any of the other Society Islands (Figure 5).
The Figure 5 topology is incompatible with the prevailing speciation model for Moorean and Tahitian Partula that views all congeners on each Windward Island as the product of a single colonization event: a Leeward Island source for Moorea and a single Moorean source for the subsequent colonization of Tahiti . The model has a narrow base of empirical support: allozyme phylogenies from 24 species of Partula that have been weighted on the basis of allele frequency ,. These data overturned previous inferences – based on reproductive relationships and morphological similarities – that Tahiti was colonized by multiple Moorean lineages . However, the allozyme data were less than robust: an unweighted (UPGMA) analysis of the same allozyme dataset yielded a substantially distinct topology in which the Tahitian species were no longer monophyletic  and a later reanalysis yielded another topology consistent with back-migration from Tahiti to the southern part of Moorea . Salient DNA phylogenies, whether using nuclear ribosomal [44, this study] or mt [35, 50–52, this study] markers have failed to recover the model’s predicted topology. The mt studies, in particular, have consistently revealed the presence of multiple exclusive sister relationships among subsections of Moorean and Tahitian Partula mt treespace [35, 50–52, this study].
Despite these incongruences, the Windward Island Partula speciation model  has remained the consensus conceptual framework for speciation studies of Moorean and Tahitian Partula species: speciation is viewed as having occurred in situ within each island ,,-,. Consistent with this view, examples of Moorean/Tahitian mt polyphyly have been interpreted either as being due to convergent molecular evolution (for restriction fragment length polymorphism data ,, or to retained ancestral mt polymorphisms -, rather than being products of inter-island gene flow involving multiple discrete lineages. However, both of these non-gene flow interpretations are problematic. The large majority of mt COI nucleotide substitutions observed among Moorean and Tahitian taxa are synonymous (data not shown) and are therefore unlikely to stem from selection-driven convergent evolution. The ancestral polymorphism interpretation requires a vicariant model of inter-island genetic differentiation that is inapplicable to Moorea and Tahiti, two spatially discrete islands that have never been joined ,.
We propose a new speciation model for Tahitian Partula taxa that involves four discrete founding lineages. Three of the four have explicit phylogenetic ties to Moorean congeners, as detailed earlier , but the largest, comprising half of Tahitian mt treespace, is of undetermined provenance (Figure 5). This new model is consistent with parallel Society Island studies that have uncovered a role for multiple independent colonizations in the evolution of Tahiti’s endemic biota .
Age of the Society Islands Partula radiation
A standout feature of the biogeography of Partulidae is the markedly asymmetric distribution of alpha diversity across the family’s 10,000 km range (Figure 1). Half of the nominal species diversity is endemic to six Society Islands spanning a mere 320 km of Oceania with a quarter being endemic to a single island: the 167 km2 Raiatea . The extraordinary concentration of 59 species of Partula within this single hotspot archipelago raises questions of taxonomic equivalency (species descriptions are typically based on shell phenotype distinctions) and lineage persistence across the range. While we cannot engage meaningfully here with the fraught topic of Partula species designations ,,, we can address the issue of lineage persistence by constructing a time-calibrated phylogeny for our molecular dataset.
These time-calibrated data imply that the heightened alpha diversity characteristic of the Society Island Partula radiation, and of the island of Raiatea, resulted at least in part from a longer uninterrupted diversification timeframe than was available elsewhere in the family’s range. For this within-archipelago radiation, cladogenesis has apparently unfolded in a qualitatively different manner to that experienced by partulids elsewhere, being characterized (until very recently) by greater lineage longevity coupled with much less lineage dispersal. Within the archipelago, the hyper-speciose Raiatean radiation is correlated with that island’s age, area and relief. Despite its relatively advanced age (Figure 5), Raiatea is second only to Tahiti in area and has retained the high island profile (1017 m maximum elevation) required for autonomous rainfall generation and rain forest formation ,. Outside of the Society Island Partula radiation, the small Palauan Partula and Eua clades also show evidence for long-term persistence within modest ranges, but lineage diversification in the Western Partula and Samoana crown clades occurred across much more extensive ranges and necessarily involved multiple episodes of among-archipelago dispersal (Figures 3 and 6).
Divergent histories of Eastern Partula and Samoana radiations
Although the genera Partula and Samoana overlap in the eastern edge of the family’s range, co-occurring in the Society Islands, they appear to have experienced distinct patterns of regional cladogenesis . Partula has many more regional species (59 vs. 15 respectively), but Samoana has a greater eastern range, encompassing the Marquesas and Austral Archipelagoes as well as the Society Islands (Figure 1). They also differ in their regional phylogenetic/population genetic profiles: eastern Samoana species have much lower collective genetic diversity levels and more pronounced phylogenetic cohesiveness ,,. The BEAST topology (Figure 6) yielded new insights into the evolutionary origins of these regional distinctions, showing that the genus Partula established eastern populations much earlier (3.27 Mya) than did Samoana (1.12 Mya). This divergent chronology probably contributed to the regional generic species richness disparity, but not to their distributional disparity. Despite its much longer regional timeframe, the genus Partula has effectively been marooned on the Society Islands, establishing just one inter-archipelago founder species in Rarotonga (Cook Islands), excepting anthropogenic introductions of P. hyalina. In contrast, the later arriving Eastern Samoana lineage has proven to be a much more effective regional disperser, on both intra-archipelago (many species have multi-island distributions: 33, 35, 54, 55, 57] and inter-archipelago scales. Why this should be is unclear, although it may be relevant that many eastern Samoana species produce exceptionally sticky mucus , and this trait could plausibly increase the likelihood of rare inter-island/archipelago phoretic dispersal events, e.g., on avian vectors.
Isolation is the primary biological characteristic of oceanic islands  and our study, placing partulid species diversity into a range-wide phylogenetic context, underlines the critical importance of biotic isolation in the evolution and diversification of this endemic Pacific island tree snail family. Most partulid species diversity – 59 species of Partula and 15 species of Samoana – is endemic to the Society, Marquesas and Austral archipelagoes. These are the three most isolated archipelagoes in the family’s range, being the most remote from Sahul/Asian continental source populations (Figure 1). Prior to human settlement, they experienced the least impact from continental biotas, e.g., they lacked terrestrial mammals (including bats) and reptiles ,. This is the rarified biotic setting that supported two very distinct generic radiations unsurpassed elsewhere in the range of Partulidae; that of Partula being the result of prolonged, uninterrupted cladogenesis within a single archipelago, that of Samoana stemming from its ability to disperse and diversify among multiple archipelagoes. Human settlement and modernity/globalization have eroded marine barrier-dependent isolation all across the range of Partulidae, directly exposing these endemic tree snails to an ever growing cast of introduced continental predators. Their continued survival outside of zoos is dependent on our ability to maintain a minimum level of biotic isolation in at least some of the islands of Oceania that could serve as refuges  for surviving members of this emblematic Pacific tree snail family.
See Table 2 for sampling locations, taxonomic identity and IUCN status of partulid snails genotyped in this study. The vast majority of our samples were obtained from museum collections and/or from captive populations that stemmed from numerous individual collecting events operating over many decades. These specimens were in large part identified by their original collectors, e.g., J.B. Burch identified his extensive 1970 samples of Society Islands partulid samples in the field using Crampton’s taxonomic descriptions . Voucher specimen information is available in the supplemental information together with sampling details for the 662 genotyped partulid snails analyzed in this study (Additional file 3: Table S1).
Total genomic DNA was isolated using the E.Z.N.A.® Mollusc DNA Kit (Omega Biotech, Norcross, GA) following the manufacturer’s instructions. A 655 nucleotide fragment of the mt cytochrome oxidase I (COI) was amplified with GoTaq DNA Polymerase (Promega, Madison, WI) using the “universal” primer pairs LCO1490/HCO2198 . A 805 (aligned) nucleotide fragment of the large nuclear ribosomal subunit gene (28S) was also amplified from a subset of specimens using D23F/D6R primers . The amplification condition consisted of initial denaturation of 95°C for 2 min, 40 cycles of denaturation (95°C, 30 sec), annealing (45°C for COI and 50°C for 28S, 30 sec), extension (72°C, 1 min), and a final extension at 72°C for 5 min. Amplification products were prepared for cycle sequencing by diluting 1:5 in sterile water. Sequencing was performed in both directions, using the PCR primers, at the University of Michigan DNA Sequencing Core Facility. All DNA sequences obtained have been deposited in GenBank (Additional file 3: Table S1). The resulting chromatograms were edited by comparing both strands using Sequencher 4.8 (Gene Codes Corporation, Ann Arbor, Michigan, USA) and the edited sequences were aligned by eye in Se-Al v2.0a11 . In total, 667 partulids were typed for mt COI and 130 snails for the large nuclear ribosomal subunit gene (28S), yielding 312 and 89 distinct genotypes, respectively.
Maximum parsimony (MP), maximum likelihood (ML) and Bayesian analyses were performed using the partulid genus Eua as the outgroup . MP searches were heuristic using PAUP* 4.0b10  with 100 random stepwise additions (MaxTrees set to 10,000) and tree-bisection-reconnection (TBR) branch-swapping. Bootstrap  branch support (BS) levels were estimated with 100 replicates (MaxTrees set to 100), 10 random additions each. ML analyses were performed using PAUP* under the best-fit model of nucleotide substitution (COI: K81uf + I + G, 28S: GTR + I + G) selected by the Akaike Information Criterion in 3.7 . The neighbor-joining tree was used as the starting tree with the likelihood parameters found in Modeltest. The branch-swapping algorithm was set to nearest neighbor interchange (NNI) for COI and to TBR for 28S dataset. Bootstrap support values were estimated using a fast-heuristic search with 100 replicates. Bayesian searches were run for 100x106 generations in MrBayes 3.1.2  set for the GTR + Γ + I model. Model parameters were treated as unknown and were estimated for each analysis. Four chains were run simultaneously and trees were sampled every 10,000 cycles. Posterior probability values were estimated by generating a 50% majority rule consensus tree after the burn-in period of 2,500 using PAUP*.
Time-calibrated BEAST analyses
We chose three time calibration points by identifying within our partulid mt COI Bayesian analysis the most comprehensively sampled, robustly supported clades that were restricted to specific archipelagos or islands of known age. Unsurprisingly, all three concerned Society Islands Partula species, by far our best-sampled archipelago and genus. Each calibration was chosen to represent a particular temporal reference point from the geological history of the archipelago. The oldest, Calibration point A (Figure 6), represents the inferred age of the oldest island in the archipelago bearing partulid tree snails: Bora Bora, 4.01-3.21 Mya . The corresponding Eastern Partula clade (PPS = 100; Figure 5) includes genotypes from all individuals sampled in the Society Islands and from Rarotonga, which we consider to be a founder population. The next oldest, Calibration point B (Figure 6), represents the inferred age of Raiatea: 2.75-2.44 Mya  and the corresponding, well-supported clade (PPS = 97; Figure 5), was exclusively composed of Raiatean taxa. The youngest, Calibration point B (Figure 6), represents the inferred age of Tahiti: 1.12-0.50 Mya , the youngest Society Island bearing partulid tree snails. This was applied to the largest Tahiti-exclusive clade (PPS = 99; Figure 5), which probably represents the oldest sampled Tahitian partulid lineage.
We used BEAST 1.7.1  to estimate a time-calibrated partulid phylogeny. This involved employing an uncorrelated lognormal distribution for the molecular clock and the SRD06 Model  for codon positions. Normal distributions were used to characterize prior distributions of the clade ages. Two independent MCMC analyses were run for 50 million iterations respectively and sampled every 5000 iterations. Convergence diagnostics were conducted in Tracer v1.5  and reliable ESS values (>200) were ensured. The first 10% trees were discarded as burn-in respectively and the trees were combined in LogCombiner . A maximum credibility tree was generated from the combined trees in TreeAnnotater .
Availability of supporting data
Museum vouchers have been deposited for genotyped snails, DNA sequences have been deposited in GenBank (see Additional file 3: Table S1 for details) and phylogenetic data have been deposited in TreeBase (reference # 16092).
DÓF and TL conceived and designed the study; TL generated the sequences incrementally over a 10-year period as museum, zoo and field samples became available; TL performed the Bayesian and Maximum Likelihood phylogenetic analyses; JL and CKCC performed the BEAST analyses; DÓF drafted the manuscript and oversaw figure design; All authors read and approved the final manuscript.
This work was supported by NSF awards DEB-0425984 & OCE 0850625 to D. Ó Foighil. J. Megahan provided expert graphical assistance. Our senior colleague, J.B. Burch, generously granted us full access to his invaluable collection of 1970’s era lyophilized partulid tissues. We gladly acknowledge other colleagues who kindly provided additional specimens: A. Baldinger, P. Bouchet, T. Coote, B. Fontaine, O. Gargominy, J. Gerber, C. Hickman, B. Holland, R. Kawamoto, F. Kraus, J-Y Meyer, P. Pearce-Kelly and R. Rundell. Our thanks to the mollusk collections of the Field Museum, Bishop Museum, Carnegie Museum and Museum of Comparative Zoology for providing reference material on loan. The manuscript was improved by the critical comments of three anonymous reviewers.
- Steadman DW: Extinction & Biogeography of Tropical Pacific Birds. 2006, University of Chicago Press, ChicagoGoogle Scholar
- Paulay G, Meyer C: Diversification in the tropical Pacific: comparisons between marine and terrestrial systems and the importance of founder speciation. Integr Comp Biol. 2002, 42: 922-934. 10.1093/icb/42.5.922.PubMedView ArticleGoogle Scholar
- Paulay G: Biodiversity on oceanic islands - its origin and extinction. Am Zool. 1994, 34: 134-144.View ArticleGoogle Scholar
- Wagner WL, Funk VJ: Hawaiian Biogeography. Evolution on a Hot Spot Archipelago. 1995, Smithsonian Institution Press, Washington DCGoogle Scholar
- Cowie RH, Holland BS: Dispersal is fundamental to biogeography and the evolution of biodiversity on oceanic islands. J Biogeogr. 2006, 33: 193-198. 10.1111/j.1365-2699.2005.01383.x.View ArticleGoogle Scholar
- Haney L: Is a new paradigm emerging for oceanic island biogeography?. J Biogeogr. 2007, 34: 753-757. 10.1111/j.1365-2699.2007.01692.x.View ArticleGoogle Scholar
- Carlquist S: Island Biology. 1974, Columbia University Press, New YorkView ArticleGoogle Scholar
- Vermeij GJ: Inequality and the directionality of history. Am Nat. 1999, 153: 243-253. 10.1086/303175.View ArticleGoogle Scholar
- Gillespie RG, Roderick GK: Arthropods on islands: evolution and conservation. Ann Rev Entomol. 2002, 47: 595-632. 10.1146/annurev.ento.47.091201.145244.View ArticleGoogle Scholar
- Cronk QCB: Islands: stability, diversity, conservation. Biodiversity Conserv. 1997, 6: 477-493. 10.1023/A:1018372910025.View ArticleGoogle Scholar
- Woinarski JCZ: Biodiversity conservation in tropical forest landscapes of Oceania. Biol Cons. 2010, 143: 2385-2394. 10.1016/j.biocon.2009.12.009.View ArticleGoogle Scholar
- Pimm SL: The Balance of Nature?. 1991, University of Chicago Press, ChicagoGoogle Scholar
- Cowie RH: Pacific Island Land Snails: Relationships, Origins and Determinants of Diversity. The Origin and Evolution of Pacific Island Biotas, New Guinea to Eastern Polynesia: Patterns and Processes. Edited by: Keast A, Miller SE. 1996, SPB Publishing, Amsterdam, 347-372.Google Scholar
- Wiles GJ, Bart J, Beck RE, Aguon CF: Impacts of the brown tree snake: patterns of decline and species persistence in Guam’s avifauna. Conserv Biol. 2003, 17: 1350-1360. 10.1046/j.1523-1739.2003.01526.x.View ArticleGoogle Scholar
- Vitousek PM: Diversity and Biological Invasions of Oceanic Islands. Biodiversity. Edited by: Wilson EO. 1988, National Academy Press, Washington, 181-189.Google Scholar
- Sax DF, Gaines SD: Species invasions and extinction: the future of native biodiversity on islands. Proc Natl Acad Sci U S A. 2008, 105: 11490-11497. 10.1073/pnas.0802290105.PubMedPubMed CentralView ArticleGoogle Scholar
- Régnier C, Fontaine B, Bouchet P: Not knowing, not recording, not listing: numerous unnoticed mollusk extinctions. Conserv Biol. 2009, 23: 1214-1221. 10.1111/j.1523-1739.2009.01245.x.PubMedView ArticleGoogle Scholar
- Hadfield MG, Miller SE, Carwile AH: The decimation of endemic Hawai’ian tree snails by alien predators. Am Zool. 1993, 33: 610-622.View ArticleGoogle Scholar
- Civeyrel L, Simberloff D: A tale of two snails: is the cure worse than the disease?. Biodiversity Conserv. 1996, 5: 1231-1252. 10.1007/BF00051574.View ArticleGoogle Scholar
- Cowie RH: Non-Indigenous Land and Freshwater Molluscs in the Islands of the Pacific: Conservation Impacts and Threats. Invasive Species in the Pacific: A Technical Review and Regional Strategy. Edited by: Sherley G. 2000, South Pacific Regional Environment Programme, Samoa, 143-172.Google Scholar
- Lydeard C, Cowie RH, Ponder WF, Bouchet P, Clark SA, Cummings KS, Frest TJ, Gargominy O, Herbert DG, Hershler R, Perez KE, Roth B, Seddon M, Strong EE, Thompson FG: The global decline of nonmarine mollusks. Bioscience. 2004, 54: 321-330. 10.1641/0006-3568(2004)054[0321:TGDONM]2.0.CO;2.View ArticleGoogle Scholar
- Clarke B, Murray J, Johnson MS: The extinction of endemic species by a program of biological control. Pac Sci. 1984, 38: 97-104.Google Scholar
- Cowie RH: Evolution and extinction of Partulidae, endemic Pacific island snails. Philos Trans R Soc London B. 1992, 335: 167-191. 10.1098/rstb.1992.0017.View ArticleGoogle Scholar
- Hopper DR, Smith BD: Status of tree snails (Gastropoda: Partulidae) on Guam, with a resurvey of sites studied by H.E. Crampton in 1920. Pac Sci. 1992, 46: 77-85.Google Scholar
- Cowie RH, Cook RP: Extinction or survival: partulid tree snails in American Samoa. Biodiversity Conserv. 2001, 10: 143-159. 10.1023/A:1008950123126.View ArticleGoogle Scholar
- Coote T, Loève E: From 61 species to five: endemic tree snails of the Society Islands fall prey to an ill-judged biological control programme. Oryx. 2003, 37: 91-96. 10.1017/S0030605303000176.View ArticleGoogle Scholar
- Clench WJ: A new subspecies of Paupina with records of land mollusks from Eastern Papua and associated islands. The Nautilus. 1936, 50: 53-54.Google Scholar
- Ó Foighil D, Lee T, Slapcinsky J: Prehistoric anthropogenic introduction of partulid tree snails in Papua New Guinean archipelagoes. J Biogeography. 2011, 38: 1625-1632. 10.1111/j.1365-2699.2011.02489.x.View ArticleGoogle Scholar
- Vagvolgyi J: Body size, aerial dispersal, and origin of the Pacific land snail fauna. Syst Zool. 1975, 24: 465-488. 10.2307/2412906.View ArticleGoogle Scholar
- Kirchner C, Kratzner R, Welter-Schultes FW: Flying snails – how far can Truncatellina (Pulmonata: Vertiginidae) be blown over the sea?. J Moll Stud. 1997, 63: 479-487. 10.1093/mollus/63.4.479.View ArticleGoogle Scholar
- Wada S, Kawakami K, Chiba S: Snails can survive passage through a bird’s digestive system. J Biogeogr. 2012, 39: 69-73. 10.1111/j.1365-2699.2011.02559.x.View ArticleGoogle Scholar
- Pilsbry HA: Manual of Conchology, Structural and Systematic. Vol. 20, Caecilioides, Glessula and Partulidae. 1909–1910, Academy of Natural Sciences of Philadelphia: Conchological Department, Philadelphia, PAGoogle Scholar
- Kondo Y: Partulidae: preview of anatomical revision. Nautilus. 1968, 81: 73-77.Google Scholar
- Lee T, Burch JB, Coote T, Fontaine B, Gargominy O, Pearce-Kelly P, Ó Foighil D: Prehistoric inter-archipelago trading of Polynesian tree snails leaves a conservation legacy. Proc R Soc Lond B. 2007, 272: 2907-2914. 10.1098/rspb.2007.1009.View ArticleGoogle Scholar
- Lee T, Burch JB, Coote T, Pearce-Kelly P, Hickman C, Meyer J-Y, Ó Foighil D: Moorean tree snail survival revisited: a multi-island genealogical perspective. BMC Evol Biol. 2009, 9: 204-10.1186/1471-2148-9-204.PubMedPubMed CentralView ArticleGoogle Scholar
- Emerson BC: Evolution on oceanic islands: molecular phylogenetic approaches to understanding pattern and process. Mol Ecol. 2002, 11: 951-966. 10.1046/j.1365-294X.2002.01507.x.PubMedView ArticleGoogle Scholar
- Sato A, Tichy H, O’hUigin C, Grant PR, Klein J: On the origin of Darwin’s finches. Mol Biol Evol. 2001, 18: 299-311. 10.1093/oxfordjournals.molbev.a003806.PubMedView ArticleGoogle Scholar
- Baldwin BG, Sanderson MJ: Age and rate of diversification of the Hawaiian silversword alliance. Proc Natl Acad Sci U S A. 1998, 95: 9402-9406. 10.1073/pnas.95.16.9402.PubMedPubMed CentralView ArticleGoogle Scholar
- Wade CM, Mordan PB, Naggs F: Evolutionary relationships among the Pulmonate land snails and slugs (Pulmonata, Stylommatophora). Biol J Linn Soc. 2006, 87: 593-610. 10.1111/j.1095-8312.2006.00596.x.View ArticleGoogle Scholar
- Tillier S, Mordan PB: The anatomy and systematics of the New Caledonian land snail genus Draparnaudia Montrouzier, 1859 (Pulmonata: Orthurethra). Zool J Linn Soc. 1995, 113: 47-91.Google Scholar
- Bouchet P, Rocroi JP, Frýda J, Hausdorf B, Ponder W, Valdés Á, Warén A: Classification and nomenclator of gastropod families. Malacologia. 2005, 47: 1-397.Google Scholar
- Pilsbry HA, Cooke CM: Partulidae of Tonga and related forms. Bernice P Bishop Mus Occas Papers. 1934, 10: 3-22.Google Scholar
- Johnson MS, Murray J, Clarke B: An electrophoretic analysis of phylogeny and evolutionary rates in the genus Partula from the Society Islands. Proc R Soc Lond B. 1986, 227: 161-177. 10.1098/rspb.1986.0017.View ArticleGoogle Scholar
- Goodacre SL, Wade CM: Molecular evolutionary relationships between partulid land snails of the Pacific. Proc R Soc Lond B. 2001, 268: 1-7. 10.1098/rspb.2000.1322.View ArticleGoogle Scholar
- Lee T, Meyer JY, Burch JB, Pearce-Kelly P, Ó Foighil D: Not completely gone: two partulid tree snail species persist on the highest peak of Raiatea, French Polynesia. Oryx. 2008, 42: 615-619. 10.1017/S0030605308001427.View ArticleGoogle Scholar
- Johnson MS, Murray J, Clarke B: Evolutionary relationships and extreme genital variation in a closely related group of Partula. Malacologia. 1993, 35: 43-61.Google Scholar
- Burch JB: A new species of land snail (Stylommatophora: Partulidae) from Raiatea, French Polynesia, Oceania. Occas Pap Mus Zool Univ Mich. 2007, 740: 8-Google Scholar
- Kondo Y, Burch JB: Archaic land snails of the Pacific islands. Malac Rev. 1972, 5: 17-18.Google Scholar
- Johnson MS, Murray J, Clarke B: The ecological genetics and adaptive radiation of Partula on Moorea. Oxford Surv Evolut Biol. 1993, 9: 167-236.Google Scholar
- Murray J, Stine OC, Johnson MS: The evolution of mitochondrial DNA in Partula. Heredity. 1991, 66: 93-104. 10.1038/hdy.1991.12.View ArticleGoogle Scholar
- Goodacre SL, Wade CM: Patterns of genetic variation in Pacific island land snails: the distribution of cytochrome b lineages among Society Island Partula. Bio J Linn Soc. 2001, 73: 131-138. 10.1111/j.1095-8312.2001.tb01351.x.View ArticleGoogle Scholar
- Goodacre SL: Population structure, history and gene flow in a group of closely related land snails: genetic variation in Partula from the Society Islands of the Pacific. Mol Ecol. 2002, 11: 55-68. 10.1046/j.0962-1083.2001.01422.x.PubMedView ArticleGoogle Scholar
- Gillespie RG, Claridge EM, Goodacre SL: Biogeography of the fauna of French Polynesia: diversification within and between a series of hot spot archipelagos. Philos Trans R Soc Lond B. 2008, 363: 3335-3346. 10.1098/rstb.2008.0124.View ArticleGoogle Scholar
- Kondo Y: Samoana of the Society Islands (Pulmonata: Partulidae). Malac Rev. 1973, 6: 19-33.Google Scholar
- Johnson MS, Murray J, Clarke B: High genetic similarities and low heterozygosities in land snails of the genus Samoana from the Society Islands. Malacologia. 1986, 27: 97-106.Google Scholar
- Kondo Y: A Revision of the Family Partulidae. 1955, PhD thesis Harvard University, Department of Biology, Cambridge, MAGoogle Scholar
- Johnson MS, Murray J, Clarke B: Parallel evolution in Marquesan partulid land snails. Bio J Linn Soc. 2000, 69: 577-598. 10.1111/j.1095-8312.2000.tb01224.x.View ArticleGoogle Scholar
- Meyer JY, Jourdan H, Malau A: Recent incursions of invasive animal and plant species in the small and remote archipelago of Wallis and Futuna (South Pacific): “La Maison Brûle mais nous regardons ailleurs”. Aliens. 2008, 27: 8-9.Google Scholar
- Rundell RJ: Diversity and conservation of the land snail fauna of the western Pacific islands of Belau (Republic of Palau, Oceania). Am Malacol Bull. 2010, 28: 81-90. 10.4003/006.028.0223.View ArticleGoogle Scholar
- Fontaine B, Gargominy O, Prié V: Indigenous Land Snails. The Natural History of Santo. Edited by: Bouchet P, Guyader H, Pascal O. 2011, Museum National d’Historie Naturelle, Patrimoines Naturels, Paris, 169-179.Google Scholar
- Pelep PO, Hadfield MG: The status of the endemic snails of the genus Partula (Gastropoda: Partulidae) on Pohnpei, Federated States of Micronesia. Micronesica. 2011, 41: 253-262.Google Scholar
- Cooke CM, Crampton HE: New species of Partula. B P Bishop Mus Occ Papers. 1930, 9: 3-5.Google Scholar
- Murray J, Clarke B: The genus Partula on Moorea: speciation in progress. Proc R Soc Lond B. 1980, 211: 83-117. 10.1098/rspb.1980.0159.View ArticleGoogle Scholar
- Clarke B, Johnson MS, Murray J: Clines in the genetic distance between two species of island land snails: how “molecular leakage” can mislead us about speciation. Philos Trans R Soc Lond B. 1996, 351: 773-784. 10.1098/rstb.1996.0072.View ArticleGoogle Scholar
- Hickman CS: Origination, Evolutionary Radiation, and Extinction in the Modern World: The Story of Partulid Tree Snails. Evolution: Investigating the Evidence. Edited by: Scotchmoor J, Springer DA. 1999, Paleontological Society Special Publication, Pittsburgh, PA, 261-280.Google Scholar
- Duncan RA, McDougall I: Linear volcanism in French Polynesia. J Volcanol Geotherm Res. 1976, 1: 197-227. 10.1016/0377-0273(76)90008-1.View ArticleGoogle Scholar
- Hildenbrand A, Gillot PY, Le Roy I: Volcano-tectonic and geochemical evolution of an oceanic intra-plate volcano: Tahiti-Nui (French Polynesia). Earth Planet. Sci Lett. 2004, 217: 349-365.Google Scholar
- Lee T, Burch JB, Jung Y, Coote T, Pearce-Kelly P, Ó Foighil D: Tahitian tree snail mitochondrial clades survived recent mass-extirpation. Curr Biol. 2007, 17: R502-R503. 10.1016/j.cub.2007.05.006.PubMedView ArticleGoogle Scholar
- Meyer JY: Threat of invasive alien plants to native flora and forest of Eastern Polynesia. Pac Sci. 2004, 58: 357-375. 10.1353/psc.2004.0032.View ArticleGoogle Scholar
- Gillespie RG, Claridge EM, Roderick GK: Biodiversity dynamics in isolated island communities: interaction between natural and human-mediated processes. Mol Ecol. 2008, 17: 45-57. 10.1111/j.1365-294X.2007.03466.x.PubMedView ArticleGoogle Scholar
- Allison A: Zoogeography of Amphibians and Reptiles of New Guinea and the Pacific Region. The Origin and Evolution of Pacific Island Biotas, New Guinea to Eastern Polynesia: Patterns and Processes. Edited by: Keast A, Miller SE. 1996, SPB Academic Publishing, Amsterdam, 407-436.Google Scholar
- Flannery TF: Mammalian Zoogeography of New Guinea and the Southwestern Pacific. The Origin and Evolution of Pacific Island Biotas, New Guinea to Eastern Polynesia: Patterns and Processes. Edited by: Keast A, Miller SE. 1996, SPB Academic Publishing, Amsterdam, 399-406.Google Scholar
- Ó Foighil D: Conservation status update on Society Island Partulidae. Tentacle. 2009, 17: 30-35.Google Scholar
- Folmer O, Black M, Hoeh W, Lutz R, Vrijenhork R: DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotech. 1994, 3: 294-299.Google Scholar
- Park JK, Ó Foighil D: Sphaeriid and corbiculid clams represent separate heterodont bivalve radiations into freshwater environments. Mol Phylogenet Evol. 2000, 14: 75-88. 10.1006/mpev.1999.0691.PubMedView ArticleGoogle Scholar
- Rambaut A: SE-AL Sequence Alignment Editor, v2.0a11. 2002, University of Oxford, OxfordGoogle Scholar
- Swofford DL: PAUP*: Phylogenetic Analysis Using Parsimony (*and Other Methods). In Version 4. 2003, Sinauer Associates, Sunderland, MAGoogle Scholar
- Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985, 39: 783-791. 10.2307/2408678.View ArticleGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.PubMedView ArticleGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.PubMedView ArticleGoogle Scholar
- Uto K, Yamamoto Y, Sudo M, Uchiumi S, Ishizuka O, Kogiso T, Tsunakawa H: New K-Ar ages of the Society Islands, French Polynesia, and implications for the Society hotspot feature. Earth Planets Space. 2007, 59: 879-885. 10.1186/BF03352750.View ArticleGoogle Scholar
- Guillou H, Maury RC, Blais S, Cotton J, Legendre C, Guille G, Caroff M: Age progression along the Society hotspot chain (French Polynesia) based on new unspiked K-Ar ages. Bull Soc Geol Fr. 2005, 176: 135-150. 10.2113/176.2.135.View ArticleGoogle Scholar
- Drummond AJ, Suchard MA, Xie D, Rambaut A: Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012, 29: 1969-1973. 10.1093/molbev/mss075.PubMedPubMed CentralView ArticleGoogle Scholar
- Shapiro B, Rambaut A, Drummond AJ: Choosing appropriate substitution models for the phylogenetic analysis of protein-coding sequences. Mol Biol Evol. 2006, 23: 7-9. 10.1093/molbev/msj021.PubMedView ArticleGoogle Scholar
- Rambaut A, Drummond AJ: Tracer v. 1.5. 2009, University of Oxford, OxfordGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.