Document Detail

Late Cretaceous vicariance in Gondwanan amphibians.
Jump to Full Text
MedLine Citation:
PMID:  17183706     Owner:  NLM     Status:  MEDLINE    
Overseas dispersals are often invoked when Southern Hemisphere terrestrial and freshwater organism phylogenies do not fit the sequence or timing of Gondwana fragmentation. We used dispersal-vicariance analyses and molecular timetrees to show that two species-rich frog groups, Microhylidae and Natatanura, display congruent patterns of spatial and temporal diversification among Gondwanan plates in the Late Cretaceous, long after the presumed major tectonic break-up events. Because amphibians are notoriously salt-intolerant, these analogies are best explained by simultaneous vicariance, rather than by oceanic dispersal. Hence our results imply Late Cretaceous connections between most adjacent Gondwanan landmasses, an essential concept for biogeographic and palaeomap reconstructions.
Ines Van Bocxlaer; Kim Roelants; S D Biju; J Nagaraju; Franky Bossuyt
Related Documents :
17571666 - Age-60 rule: the end is in sight.
3130766 - Automatic piercing spinhaler (halermatic): a comparative study.
3565426 - Underlying cardiac lesions in adults with infective endocarditis. the changing spectrum.
17183706 - Late cretaceous vicariance in gondwanan amphibians.
20634306 - Solitaire flow-restoration device for treatment of acute ischemic stroke: safety and re...
23086786 - Right subclavian approach as a feasible alternative for transcatheter aortic valve impl...
Publication Detail:
Type:  Journal Article; Research Support, Non-U.S. Gov't     Date:  2006-12-20
Journal Detail:
Title:  PloS one     Volume:  1     ISSN:  1932-6203     ISO Abbreviation:  PLoS ONE     Publication Date:  2006  
Date Detail:
Created Date:  2007-01-08     Completed Date:  2010-03-16     Revised Date:  2013-06-06    
Medline Journal Info:
Nlm Unique ID:  101285081     Medline TA:  PLoS One     Country:  United States    
Other Details:
Languages:  eng     Pagination:  e74     Citation Subset:  IM    
Biology Department, Unit of Ecology and Systematics, Vrije Universiteit Brussel, Brussels, Belgium.
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Amphibians / classification*,  genetics*
Anura / classification,  genetics
Biological Evolution*
Models, Biological
Time Factors

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): PLoS ONE
Journal ID (publisher-id): plos
ISSN: 1932-6203
Publisher: Public Library of Science, San Francisco, USA
Article Information
Download PDF
Van Bocxlaer et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Received Day: 7 Month: 9 Year: 2006
Accepted Day: 15 Month: 9 Year: 2006
collection publication date: Year: 2006
Electronic publication date: Day: 20 Month: 12 Year: 2006
Volume: 1 Issue: 1
E-location ID: e74
ID: 1762348
Publisher Id: 06-PONE-RA-00097
DOI: 10.1371/journal.pone.0000074
PubMed Id: 17183706

Late Cretaceous Vicariance in Gondwanan Amphibians Alternate Title:Gondwanan Amphibian Vicariance
Ines Van Bocxlaer1
Kim Roelants1
S.D. Biju12
J. Nagaraju3
Franky Bossuyt1*
Michael Hofreiteredit1 Role: Academic Editor
1Biology Department, Unit of Ecology and Systematics, Vrije Universiteit BrusselBrussels, Belgium
2Centre for Environmental Management of Degraded Ecosystems (CEMDE), School of Environmental Studies, University of DelhiDelhi, India
3Laboratory of Molecular Genetics, Centre for DNA Fingerprinting and DiagnosticsHyderabad, India
Max Planck Institute for Evolutionary Anthropology, Germany
Correspondence: * To whom correspondence should be addressed. E-mail:
Contributed by footnote: Conceived and designed the experiments: KR FB JN SB IV. Performed the experiments: KR IV. Analyzed the data: KR FB IV. Contributed reagents/materials/analysis tools: KR FB JN SB IV. Wrote the paper: KR FB IV.


Ever since Alfred Wegener's theory in 1915, biologists have had vicariance at their disposal as a mechanism to logically explain transoceanic distributions in plants and animals. Meta-analyses of area cladograms [1] show that many animal phylogenies match the sequential break-up of Gondwana and thus support this hypothesis. However, there is an increasing number of Southern Hemisphere taxa for which divergence time estimates do not readily fit the temporal framework of fragmentation [2], [3]. This lack of concordance has sometimes led to the conclusion that the distribution of these organisms resulted from overseas dispersal events [e.g. 4], rather than from a vicariant history. Although transoceanic dispersal may have played an important role in several taxa [5], the possibility that terrestrial connections have existed longer than assumed under traditional geological models, has received limited attention in the dispersal-vicariance debate. Yet, if theories of prolonged connections [6], intervening landmasses [2], [3], or a completely enclosed Pacific Basin [7] in the Late Cretaceous would prove to be correct, several of the younger phylogenies might still be explained by Gondwanan vicariance.

While geology is typically used to calibrate the Tree of Life, molecular timetrees of organisms with no (or very limited) oceanic dispersal capabilities in turn have the property that they can test for the presence of associations among tectonic plates. Indeed, under this premise, divergence times between continent-scale endemic groups provide a maximum age for the actual separation of landmasses (i.e., the disappearance of a terrestrial connection). Amphibians are an extremely good model for such studies, because they are notoriously salt-intolerant and, as a rule, oceans provide an effective barrier against their transoceanic dispersal [5]. An exceptional case of short-distance transoceanic dispersal to the Comores has been shown [8], but amphibians do not occur on other oceanic islands and there is no evidence for long-distance overseas dispersal in these animals. Thus it is surprising that molecular clock analyses date the origin of two species-rich Gondwanan groups, Natatanura and Microhylidae, well in the Cretaceous [8]?[11]. If correct, the subsequent diversification of the major lineages within these clades cannot be congruent with initial Gondwanan break-up, which is traditionally depicted in the Early Cretaceous (145?100 mya). To address this controversy, and to test the hypothesis of prolonged biotic interchange between Gondwanan landmasses, we checked for phylogenetic and temporal concordance between and within Microhylidae and Natatanura, using 3.1 and 5.6 kb of sequence data for the major lineages, respectively [12].


Our Maximum likelihood (ML) and Bayesian phylogenetic trees differ from Frost et al.'s ?Amphibian Tree of Life? [13] in some important biogeographic aspects (Figure 1, Table 1). First, in Microhylidae, we recover Dyscophinae (not including the Asian Calluella) as the closest relatives of Microhylinae (i.e., a Madagascar-Eurasia/India sister-relationship) (Figure 1) instead of its alternative position as sister group of Asterophryinae. The alternative of Scaphiophryninae being the sistergroup of Microhylinae [13], although implying the same biogeographic pattern, is rejected by our data (Table 1, hypotheses 1?3). Second, our analyses do not support a sister-clade relationship between the African genus Hoplophryne and Ramanella from the Indian subcontinent (Table 1, hypothesis 4), but recover the latter in an unexpected clade endemic to the Indian Subcontinent (Figure 1A, Indian clade). Third, our analyses show that Melanobatrachinae and Microhylinae are restricted to India and India/Eurasia, respectively, and that their assumed relatives in Africa and South America represent separate lineages (Table 1, hypothesis 5 and 6). Fourth, Kalophrynus as sistergroup of all remaining microhylids, which would imply a basal position for a strictly Eurasian lineage, is also rejected by our data (Table 1, hypothesis 7). Finally, we find strong support for an Indian endemic clade composed of the genera Indirana and Micrixalus (Figure 1B). The nested position of Indirana within six hierarchic African clades, as found by Frost et al. [13] is strongly rejected (Table 1, hypotheses 8?13), identifying Africanura as an exclusively African radiation (Figure 1B). All of our alternatively recovered relationships favor simpler biogeographic scenarios than previous hypotheses and complete the evidence that amphibians experienced continent-scale endemic radiations [11] as the result of a general vicariance history with little overseas dispersal.

All hypotheses either represent long-established taxa [29]?[30] or have been proposed recently [13].

To estimate the age of early divergences in Microhylidae and Natatanura, we combined both clades in a single data matrix and performed relaxed molecular clock analyses using two independent methods and multiple combinations of calibration points [12]. Our age estimates are largely congruent with recent studies focusing on deeper anuran divergences or on natatanuran relationships [8]?[11]. Importantly, all our timetrees placed the 95% confidence intervals of early divergences in both groups in the Late Cretaceous (Figure 2A, gray block) and Early Tertiary (Table 2). These results indicate that each of the endemic groups became isolated long after the traditionally depicted break-up events [14].

For both the TK- and PL-method, the ages represent the 95% confidence intervals of the analysis that includes all calibration points except G [12]. Node numbers and colors associated with break-up correspond with figure 2. KP, Kerguelen Plateau.

We used Dispersal-Vicariance analyses [12] to reconstruct the biogeography of Late Cretaceous and Early Tertiary nodes in Natatanura and Microhylidae. This method has the advantage of not being constrained by a particular break-up sequence, and thus does not favor the traditional model of Gondwanan fragmentation a priori. Our analyses suggest six vicariance events among the following adjacent Gondwanan landmasses (Figure 2, Table 2, nodes 1, 2, 4, 5, 8, 9): Indo-Madagascar versus Africa (blue), Indo-Madagascar versus Australia-New Guinea (orange), Africa versus South America (yellow), and South America-Antarctica versus Indo-Madagascar (over the Kerguelen plateau) (green). Most of our inferred 95% credibility intervals for the divergence time estimates of these nodes show no overlap with the temporal windows of traditional plate-tectonic models, which indicate that Indo-Madagascar separated from Africa and Antarctica-Australia-New Guinea in the Early Cretaceous (165?121 mya and 130?120 mya, respectively), and Madagascar from India about 88 mya [14]. However, the strong temporal overlap between 95% age intervals of equally reconstructed splits (compare the range of identically colored intervals in Figure 2) suggests that they were shaped by the same paleogeographic events and hence, that they reflect the actual separation of the corresponding landmasses.


Under a strict interpretation of the DIVA results, the Madagascar-Eurasia/India relationships require three contemporary dispersals (Figure 2, purple) directly from Madagascar to Eurasia (i.e., transoceanic, not using India). However, accretion of landmasses, like vicariance, also allows multiple groups to simultaneously depict the same change in their distribution [15], [16]. We suggest that this contemporary, repeated pattern of area relationships is a case of ?mass coherent dispersal? [15], [16], resulting from taxa on isolated India (or Indo-Madagascar) dispersing to Eurasia after contact [17]. Under this interpretation, the Madagascar-Eurasia splits in fact reflect the vicariance between Madagascar and India (Figure 2, purple). The three extinctions needed on India to accept this scenario can be easily explained by mass extinction [12], given that the Deccan Traps flooded large parts of the subcontinent around the Cretaceous-Tertiary (KT) boundary.

Our combined phylogenetic, temporal, and biogeographic evidence in two frog clades is indicative of Late Cretaceous terrestrial connections between most adjacent Gondwanan landmasses. Although this contrasts with traditional plate-tectonic models [14], various alternative hypotheses have been proposed that could accommodate our results. First, our phylogenies suggest that African frog lineages in both groups became isolated first (i.e., an ?Africa first? scenario [2], [18]). Vicariance events in the early Late Cretaceous (Table 2, nodes 2 and 9) are consistent with the prolonged existence of a South America-Africa connection [18], as well as the postulation of a ?Central Corridor? between Africa and Indo-Madagascar until ?90 mya [6]. Second, the Kerguelen Plateau (Figure 2B) and Gunnerus Ridge, that may have connected Indo-Madagascar and Antarctica until ?80 mya [19]?[22] provide an explanation for microhylid vicariances between South America and Indo-Madagascar. Third, paleontological and paleogeographic studies have recently suggested that both the Tethys Sea and Pacific Basin may have been much smaller and more enclosed by surrounding continents than previously assumed [7], [23]. As a consequence, India may have accreted to Eurasia already in the Late Cretaceous, while it was still connected to Madagascar over the Seychelles Plateau [22], [24]. This would account for the triple, deeply nested Madagascar-Eurasia/India relationships recovered in our phylogenies.

A relaxed scenario for Gondwanan break-up helps to explain several observations in the fossil record. For example, several mammalian, dinosaur, and crocodilian taxa seem to have attained widespread distributions in Gondwana only by the Late Cretaceous [2], [18], [20]. Furthermore, several palaeontological studies have suggested close affinities of South East Asian taxa with Indian, Madagascan, and Australia-New Guinean taxa in the Late Cretaceous [25]?[27]. Although fossils inevitably underestimate the true age of divergences, these observations provide a better fit to a model that allows both prolonged biotic interchange between landmasses of Gondwanan origin, and a Late Cretaceous connection of Indo-Madagascar and Australia-New Guinea with South East Asia. Landmasses that experienced recognized mass extinctions, such as Antarctica, the Kerguelen Plateau, and India, must have played a central role in Late Cretaceous range extensions. Future palaeobiogeographical research on biota of Southern Hemispere origin will clearly benefit from key fossil evidence from these regions.

Materials and Methods

Phylogenetic relationships within Microhylidae and Natatanura were estimated using independently assembled data sets, each encompassing the major lineages within one family. The microhylid data set contains 2865 unambiguously aligned nucleotides, sampled from three nuclear protein-coding genes (Cxcr-4, Ncx-1, Rag-1) and one mitochondrial gene (16S rDNA) (Table S1). The ingroup is composed of 28 microhylid species, covering all previously and currently recognized subfamilies within this family [13], [28]?[31] (Table S2). Four natatanuran and five afrobatrachian species served as outgroup. The natatanuran data set contains 4446 unambiguously aligned nucleotides, combining the same three nuclear gene fragments mentioned above with coding regions of two additional nuclear genes (Tyr and Rhod), and an extended mitochondrial fragment (12S rDNA+tRNA(Val)+16S rDNA) (Table S3). The ingroup comprises 31 natatanuran species, while one dendrobatid, four afrobatrachian and five microhylids were selected as outgroup (Table S4). Genbank accession numbers are summarized in table S5 (Microhylidae), table S6 (Natatanura), and table S7 (outgroup species).

Phylogenetic analyses

Phylogenetic relationships were reconstructed for both data sets using Bayesian analyses and maximum likelihood (ML) searches (Figure 1). All analyses implemented the GTR+G+I model of DNA substitution, which received the best Akaike information score by Modeltest 3.06 [32]. Bayesian analyses were performed with MrBayes 3.1.2 [33], using a mixed model according to a gene-based data partition. Two runs of four MCMC chains each were executed in parallel for five million generations, with a sampling interval of 1000 generations and a burn-in corresponding to the first one million generations. Convergence of the parallel runs was confirmed by split frequency standard deviations (<0.001), and by potential scale reduction factors (PSRF; ?1.0) for all model parameters. Posterior probabilities for clades were obtained by combining the post-burn-in trees from parallel runs in a single consensus tree. Heuristic ML searches were performed with PAUP* 4.0b10 [34] and involved multiple rounds of TBR branch swapping, alternatively starting from Bayesian consensus trees or from ML trees estimated by the program Phyml 2.4.1 [35]. All searches implemented empirical nucleotide frequencies and fixed rate parameters, estimated in advance from the starting trees. Clade support under the likelihood criterion was assessed by analyzing 1000 nonparametric bootstrap replicates using Phyml. Bayesian posterior probabilities for alternative microhylid and natatanuran phylogenetic hypotheses [13], [29]?[30] were estimated by screening the post-burn-in trees sampled by MrBayes using topological constraint filters in PAUP*.

Tree calibration and dating

We performed dating analyses using two different relaxed molecular clock methods, which have recently been demonstrated to be the least sensitive to taxon sampling [36], and which have complementary advantages and limitations. Thorne & Kishino's (TK) method [37] accommodates unlinked rate variation across different loci (a ?multi-gene? approach), allows the use of time constraints on multiple divergences, and uses a Bayesian MCMC approach to approximate the posterior distribution of divergence times and rates, but uses an F84+G model (or a nested variant) for branch length estimation, and fails to incorporate phylogenetic uncertainty in the posterior distribution. Sanderson's penalized likelihood (PL) method [38] allows branch length estimation using more complex DNA substitution models (e.g., GTR+G+I) and the use of a posterior tree set to estimate credibility intervals (CI), but necessarily averages rate variation over all loci (a ?supergene? approach), and requires a time-consuming cross-validation method to determine optimal rate smoothing penalty parameters. To maximize the overall accuracy of our dating estimates, we sought to obtain an optimal phylogenetic coverage of calibration points across our tree (Figure S1, Table S8). However, to minimize the risk of over-constraining the resulting timetree, we used them only as minimum time constraints. The results of TK- and PL-analyses with different combinations of calibration points are listed in Table S9. Removal of individual time constraints in most cases resulted in highly congruent dating estimates with respect to the total set of calibration points (Figure S2), with the exception of two paleogeographic calibration points. However, these estimates resulted in a younger age for reconstructed natatanuran and microhylid vicariance events, and thus would imply an even larger discrepancy with currently accepted models for Gondwana break-up.

Biogeographic analyses

To quantify vicariance events in the early diversification of both groups, we estimated ancestral distributions for ingroup nodes using dispersal-vicariance analyses performed by the program DIVA 1.1 [39]. All ingroup taxa were coded for their present distribution across six biogeographic units (five of which represent separate fragments of Gondwana): Africa, Madagascar, India, South America, Australia-New Guinea and Eurasia (Table S10, Figures S3, S4). To simplify the analyses, we saved an extra biogeographic unit by coding the North American Gastrophryne olivacea as a South American taxon, because of its nested position within the neotropical clade Gastrophryninae.

Supporting Information Materials and Methods S1

(0.17 MB DOC)

Click here for additional data file (pone.0000074.s001.doc)

Figure S1

Phylogenetic tree used for the divergence time analyses. The natatanuran and microhylid clades are resolved according to the ML topologies obtained from their respective data sets (Figure 1). Outgroup divergences are resolved according to previous phylogenetic evidence (see Methods). Numbers at ingroup nodes are cross-referenced in Table S9, Letters indicate calibration points.

(11.84 MB TIF)

Click here for additional data file (pone.0000074.s002.png)

Figure S2

Effect of excluding individual calibration points on divergence time estimates. Letters refer to excluded calibration points and correspond to Tables S8 and S9. The bars represent mean age differences, obtained with the TK-method (orange) and the PL method (blue), compared to the analysis including all calibration points except G.

(5.34 MB TIF)

Click here for additional data file (pone.0000074.s003.png)

Figure S3

DIVA-reconstruction of ancestral distribution areas. Letter codes at internal nodes correspond to landmasses depicted on the inset globe.

(8.35 MB TIF)

Click here for additional data file (pone.0000074.s004.png)

Figure S4

DIVA-reconstruction of ancestral distribution areas, assuming three divergences on Indo-Madagascar. Reconstruction of ancestral distributions under the assumption that nodes 3, 6 and 7 represent vicariance events related to India-Madagascar break-up. Letter codes at internal nodes correspond to landmasses depicted on the inset globe.

(8.35 MB TIF)

Click here for additional data file (pone.0000074.s005.png)

Table S1

Summary of the sequence data for all sampled nuclear and mitochondrial gene fragments and the total dataset in Microhylidae.

(0.31 MB DOC)

Click here for additional data file (pone.0000074.s006.doc)

Table S2

Taxa included in the Microhylidae dataset, with voucher numbers and origin.

(0.34 MB DOC)

Click here for additional data file (pone.0000074.s007.doc)

Table S3

Summary of the sequence data for all sampled nuclear and mitochondrial gene fragments and the total dataset in Natatanura.

(0.31 MB DOC)

Click here for additional data file (pone.0000074.s008.doc)

Table S4

Taxa included in the Natatanura dataset, with voucher numbers and origin.

(0.34 MB DOC)

Click here for additional data file (pone.0000074.s009.doc)

Table S5

Taxa with GenBank accession numbers of homologous gene fragments for microhylid species.

(0.34 MB DOC)

Click here for additional data file (pone.0000074.s010.doc)

Table S6

Taxa with GenBank accession numbers of homologous gene fragments for natatanuran species.

(0.36 MB DOC)

Click here for additional data file (pone.0000074.s011.doc)

Table S7

Taxa with GenBank accession numbers of homologous gene fragments for outgroup species.

(0.37 MB DOC)

Click here for additional data file (pone.0000074.s012.doc)

Table S8

Calibration points used in this study.

(0.34 MB DOC)

Click here for additional data file (pone.0000074.s013.doc)

Table S9

Divergence time estimates (standard deviation; 95% credibility intervals) for all ingroup nodes in our dating-tree, obtained with different calibration points and dating methods.

(0.16 MB XLS)

Click here for additional data file (pone.0000074.s014.xls)

Table S10

Dating estimates (mya) for Microhylidae and Natatanura, with DIVA reconstruction of vicariance events in their early evolution. For both the TK- and PL-method, the ages represent the 95% confidence intervals of the analysis that includes all calibration points except G. Vicariance (1) represents the strict analysis; Vicariance (2) represents the analysis with the three India-Eurasia splits constrained as India-Madagascar. The two DIVA reconstructions represent the actual distribution of the frogs, while the associated continental break-up is the presumed geological unit at the moment of break-up, based on our divergence times. Node numbers and colors associated with break-up correspond with figure 2. Afr?=?Africa, Ind?=?India, SAm?=?South America, Mad?=?Madagascar, Aus?=?Australia-New Guinea, Eur?=?Eurasia, KP?=?Kerguelen Plateau.

(0.33 MB DOC)

Click here for additional data file (pone.0000074.s015.doc)


Competing Interests: The authors have declared that no competing interests exist.

Funding: Ines Van Bocxlaer is supported by the Institute for the Promotion of Innovation through Science and Technology in Flanders (IWT-Vlaanderen). Franky Bossuyt receives a postdoctoral fellowship from the Fonds voor Wetenschappelijk Onderzoek-Vlaanderen. This work was supported by research grants FWO, FWO G.0056.03, FWO G.0307.04 and VUB OZR834.

We are grateful to David J. Gower, Simon Loader, David C. Cannatella, Robert F. Inger and Steve Donnellan, David B. Wake, Miguel Vences, and Robert C. Drewes for providing tissue samples, and to Mark Wilkinson for commenting on an earlier draft of this manuscript. We thank An Mannaert and Katleen Willibal for lab assistance.

1. Sanmartin I,Ronquist F. 2004;Southern Hemisphere biogeography inferred by event-based models: plant versus animal patterns.Syst Biol 53:216–243. [pmid: 15205050]
2. Sampson SD,Witmer LM,Forster CA,Krause DW,O'Connor PM,Dodson P,Ravoavy F. 1998;Predatory Dinosaur remains from Madagascar: implications for the Cretaceous biogeography of Gondwana.Science 280:1048–1051. [pmid: 9582112]
3. Cooper A,Lalueza-Fox C,Anderson S,Rambaut A,Austin J,et al. 2001;Complete mitochondrial genome sequences of two extinct moas clarify ratite evolution.Nature 409:704–707. [pmid: 11217857]
4. Raxworthy CJ,Forstner MRJ,Nussbaum RA. 2002;Chameleon radiation by oceanic dispersal.Nature 415:784–787. [pmid: 11845207]
5. de Queiroz A. 2005;The resurrection of oceanic dispersal in historical biogeography.TREE 20:68–73. [pmid: 16701345]
6. Chatterjee S,Scotese CR. 1999;The breakup of Gondwana and the evolution and biogeography of the Indian plate.Proc Indian Natl Sci Acad 65A:397–425.
7. McCarthy D. 2005;Biogeographical and geological evidence for a smaller, completely-enclosed Pacific Basin in the Late Cretaceous.J Biogeogr 32:2161–2177.
8. Vences M,Vieites DR,Glaw F,Brinkmann H,Kosuch J,et al. 2003;Multiple overseas dispersal in amphibians.Proc R Soc Lond B 270:2435–2442.
9. San Mauro D,Vences M,Alcobendas M,Zardoya R,Meyer A. 2005;Initial diversification of living amphibians predated the breakup of Pangaea.Am Nat 165:590–599. [pmid: 15795855]
10. van der Meijden A,Vences M,Hoegg S,Meyer A. 2005;A previously unrecognized radiation of ranid frogs in Southern Africa revealed by nuclear and mitochondrial DNA sequences.Mol Phylogenet Evol 37:674–685. [pmid: 15975829]
11. Bossuyt F,Brown RM,Hillis DM,Cannatella DC,Milinkovitch MC. 2006;Phylogeny and biogeography of a cosmopolitan frog radiation: Late Cretaceous diversification resulted in continent-scale endemism in the family Ranidae.Syst Biol 55:579–594. [pmid: 16857652]
12. Materials and Methods are available as supporting material on PLoS One Online
13. Frost DR,Grant T,Faivovich J,Bain RH,Haas A,et al. 2006;The amphibian Tree of Life.Bull Am Mus Nat Hist 297:1–370.
14. Smith AG,Smith DG,Funnell BM. 1994Atlas of Mesozoic and Cenozoic coastlines. Cambridge, UK: Cambridge University Press; :99 p.
15. Upchurch P,Hunn CA. 2002;?Time?: the neglected dimension in cladistic biogeography?Geobios 35:277–286.
16. Turner AH. 2004;Crocodyliform biogeography during the Cretaceous: evidence of Gondwanan vicariance from biogeographical analysis.Proc Biol Sci 271:2003–2009. [pmid: 15451689]
17. Bossuyt F,Milinkovitch MC. 2001;Amphibians as indicators of early tertiary ?out-of-India? dispersal of vertebrates.Science 292:93–95. [pmid: 11292870]
18. Sereno PC,Wilson JA,Conrad JL. 2004;New dinosaurs link southern landmasses in the Mid-Cretaceous.Proc R Soc Lond B. 271:1325–1330.
19. Karanth KP. 2006;Out-of-India Gondwanan origin of some tropical Asian biota.Current Science 90:789–792.
20. Krause DW,Prasad GVR,von Koenigswald W,Sahni A,Grine FE. 1997;Cosmopolitanism among Late Cretaceous Gondwanan mammals.Nature 390:504–507.
21. Hay WW,DeConto RM,Wold CN,Willson KM,Voigt S,et al. Barrera E,Johnson C1999;An alternative global Cretaceous paleogeography,Evolution of the Cretaceous ocean-climate system. Geological Society of America Special Paper 332:1–48.
22. Rage JC. 2003;Relationships of the Malagasy fauna during the Late Cretaceous: Northern or Southern routes?Acta Palaeontol Pol 48:661–662.
23. Briggs JC. 2003;The biogeographic and tectonic history of India.J Biogeogr 30:381–388.
24. Patriat P,Segoufin J. 1988;Reconstruction of the Central Indian Ocean.Tectonophysics 155:211–234.
25. Prasad V,Str?mberg CAE,Alimohammadian H,Sahni A. 2005;Dinosaur coprolites and the early evolution of grasses and grazers.Science. 310:1177–1180. [pmid: 16293759]
26. Azuma Y,Curie PJ. 2000;A new carnosaur (Dinosauria: Theropoda) from the Lower Cretaceous of Japan.Can J Earth Sci 37:1735–1753.
27. Hill RV,Witmer LM,Norell MA. 2003;A new specimen of Pinacosaurus grangeri (Dinosauria: Ornithischia) from the Late Cretaceous of Mongolia: ontogeny and phylogeny of Ankylosaurs.American Museum Novitates 3395:1–30.
28. Dubois A. 2005;Amphibia Mundi. 1.1. An ergotaxonomy of recent amphibians.Alytes 23:1–24.
29. Duellman WE,Trueb L. 1994Biology of amphibians. Baltimore, USA: Johns Hopkins University Press; :670.
30. Parker HW. 1934A monograph of the frogs of the family Microhylidae, London: British Museum (Natural History); :1–203.
31. Savage JM. 1973The geographic distribution of frogs: patterns and predictions.Evolutionary Biology of the Anurans: contemporary research on major problems Vial JLColumbia: Univ. Missouri Press; :351–445.
32. Posada D,Crandall KA. 1998;MODELTEST: testing the model of DNA substitution.Bioinformatics 14:817–818. [pmid: 9918953]
33. Ronquist F,Huelsenbeck JP. 2003;MrBayes 3: Bayesian phylogenetic inference under mixed models.Bioinformatics 19:1572–1574. [pmid: 12912839]
34. Swofford DL. 2003PAUP*. Phylogenetic analyses using parsimony (* and other methods).. version 4.0 edSunderland, Massachusetts: Sinauer Associates;
35. Guindon S,Gascuel O. 2003;A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood.Syst Biol 52:696–704. [pmid: 14530136]
36. Linder HP,Hardy CR,Rutschmann F. 2005;Taxon sampling effects in molecular clock dating: an example from the African Restionaceae.Mol Phylogenet Evol 35:569–582. [pmid: 15878126]
37. Thorne JL,Kishino H. 2002;Divergence time and evolutionary rate estimation with multilocus data.Syst Biol 51:689–702. [pmid: 12396584]
38. Sanderson MJ. 2002;Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood approach.Mol Biol Evol 19:101–109. [pmid: 11752195]
39. Ronquist F. 1997;Dispersal-vicariance analysis: a new approach to the quantification of historical biogeography.Syst Biol 46:195–203.


[Figure ID: pone-0000074-g001]
doi: 10.1371/journal.pone.0000074.g001.
Figure 1  Maximum likelihood phylograms for (A) Microhylidae (-lnL?=?24126.41696) and (B) Natatanura (-lnL?=?53106.71363).

Branch support is indicated as follows: black circles: ML bootstrap support (BS)?75% and Bayesian posterior probability (PP)?0.95; triangle pointing to the right: BS<75%, PP?0.95; triangle pointing to the left: BS?75%, PP<0.95; white circles: BS<75% and PP<0.95.

Higher taxon names follow Frost et al. [13], with the exception of removal of Hoplophryne from Melanobatrachinae, and recognition of Ranixalidae.

[Figure ID: pone-0000074-g002]
doi: 10.1371/journal.pone.0000074.g002.
Figure 2  Late Cretaceous vicariance in Microhylidae and Natatanura.

(A) Molecular timetree (TK method, all calibration points except G [12]).

Horizontal colored bars and lines at internal nodes (Standard deviation and 95% credibility interval, respectively) indicate vicariance events reconstructed by DIVA-analyses, and interpreted as follows: orange: Australia <?> Indo-Madagascar; yellow: Africa <?> South America; blue: Africa <?> Indo-Madagascar; purple: Madagascar <?> India (Seychelles); green: S. America-Antarctica <?> Indo-Madagascar (the intervening Kerguelen Plateau being involved).

The latter splits in our timetree are interpreted as vicariance events between the Kerguelen plateau and Antarctica or Indo-Madagascar [12].

The branches denoting the latest colonization of Eurasia, as reconstructed by DIVA, are indicated by an asterisk.

Numbers at terminals correspond to taxon numbers in figure 1. (B) Late Cretaceous Gondwana, with indication of corresponding geological break-ups.

Abbreviations: AF?=?Africa, MA?=?Madagascar, IN?=?India, EU?=?Eurasia, SA?=?South America, AN?=?Antarctica, AU?=?Australia-New Guinea, KP?=?Kerguelen Plateau.

[TableWrap ID: pone-0000074-t001] doi: 10.1371/journal.pone.0000074.t001.
Table 1  Bayesian test for alternative hypotheses (and their associated biogeographic implication) of relationships in Microhylidae and Natatanura.
Alternative Hypothesis Reference Biogeographic implication p
(1) Monophyly of Dyscophinae (including Calluella) [29], [30] Dyscophinae not endemic to Madagascar < 0.001
(2) Dyscophinae sistergroup of Asterophryinae [13] Madagascar - Australia/New Guinea link < 0.001
(3) Scaphiophryninae sistergroup of Microhylinae [13] alternative Madagascar - India/Eurasia link < 0.001
(4) (Ramanella, Hoplophryne) clade [13] India - Africa link < 0.001
(5) Monophyly of Melanobatrachinae (Hoplophryne + Melanobatrachus) [13] Melanobatrachinae not endemic to India < 0.001
(6) Monophyly of Microhylinae s.l. (including South American taxa) [29], [30] Microhylinae not endemic to Eurasia/India < 0.001
(7) Kalophrynus sistergroup of all remaining Microhylidae [13] Early divergence of Eurasian lineage < 0.001
(8) Monophyly of Africanura [13] Africanura not endemic to Africa 0.004
(9) Monophyly of Pyxicephaloidea [13] Pyxicephaloidea not endemic to Africa < 0.001
(10) Monophyly of Petropedetidae [13] Petropedetidae not endemic to Africa < 0.001
(11) Monophyly of (Petropedetes + Arthroleptides + Indirana) [13] Petropedetidae not endemic to Africa < 0.001
(12) Monophyly of Telmatobatrachia [13] no initial isolation of African lineages 0.001
(13) Monophyly of Ametrobatrachia [13] no initial isolation of African lineages < 0.001

[TableWrap ID: pone-0000074-t002] doi: 10.1371/journal.pone.0000074.t002.
Table 2  Dating estimates (mya) for vicariance events in the early evolution of Microhylidae and Natatanura.
Node Age (TK) Age (PL) continental break-up associated with vicariance
9 [74.2, 107.7] [61.7, 75.0] Africa <?> Indo-Madagascar (blue)
2 [74.2, 97.1] [68.2, 74.9] Africa <?> South America (yellow)
Africa <?> Indo-Madagascar (blue)
1 [65.4, 88.1] [51.7, 69.1] S. America-Antarctica <?> Indo-Madagascar (disappearance KP) (green)
4 [70.4, 91.4] [66.8, 71.4] S. America-Antarctica <?> Indo-Madagascar (disappearance KP) (green)
8 [63.1, 95.0] [52.3, 68.5] Australia <?> Indo-Madagascar (orange)
5 [66.7, 84.5] [65.2, 68.5] Australia <?> Indo-Madagascar (orange)
3 [48.7, 75.8] [49.6, 66.1] Madagascar <?> India (Seychelles) (purple)
6 [65.1, 77.6] [65.0, 65.0] Madagascar <?> India (Seychelles) (purple)
7 [56.0, 86.2] [50.5, 62.8] Madagascar <?> India (Seychelles) (purple)

Article Categories:
  • Research Article
Article Categories:
  • Evolutionary Biology

Previous Document:  The syntax and meaning of wild gibbon songs.
Next Document:  Control of canalization and evolvability by Hsp90.