Different processes lead to similar patterns: a test of codivergence and the role of sea level and climate changes in shaping a southern temperate freshwater assemblage
© Barber et al; licensee BioMed Central Ltd. 2011
Received: 29 March 2011
Accepted: 25 November 2011
Published: 25 November 2011
Understanding how freshwater assemblages have been formed and maintained is a fundamental goal in evolutionary and ecological disciplines. Here we use a historical approach to test the hypothesis of codivergence in three clades of the Chilean freshwater species assemblage. Molecular studies of freshwater crabs (Aegla: Aeglidae: Anomura) and catfish (Trichomycterus arealatus: Trichomycteridae: Teleostei) exhibited similar levels of genetic divergences of mitochondrial lineages between species of crabs and phylogroups of the catfish, suggesting a shared evolutionary history among the three clades in this species assemblage.
A phylogeny was constructed for Trichomycterus areolatus under the following best-fit molecular models of evolution GTR + I + R, HKY + I, and HKY for cytochrome b, growth hormone, and rag 1 respectively. A GTR + I + R model provided the best fit for both 28S and mitochondrial loci and was used to construct both Aegla phylogenies. Three different diversification models were observed and the three groups arose during different time periods, from 2.25 to 5.05 million years ago (Ma). Cladogenesis within Trichomycterus areolatus was initiated roughly 2.25 Ma (Late Pliocene - Early Pleistocene) some 1.7 - 2.8 million years after the basal divergences observed in both Aegla clades. These results reject the hypothesis of codivergence.
The similar genetic distances between terminal sister-lineages observed in these select taxa from the freshwater Chilean species assemblage were formed by different processes occurring over the last ~5.0 Ma. Dramatic changes in historic sea levels documented in the region appear to have independently shaped the evolutionary history of each group. Our study illustrates the important role that history plays in shaping a species assemblage and argues against assuming similar patterns equal a shared evolutionary history.
KeywordsAegla biogeography Chile cladogenesis freshwater Trichomycterus
Understanding how species assemblages have been formed and maintained is a central question in biology [1–4]. Historical processes, such as river capture and climate change, likely played significant roles in driving the formation of aquatic species assemblages [5–7], whereas ecological processes may help maintain species interactions in co-occurring species if in fact those species are interacting (i.e form a community). Understanding the relative roles these two processes have had in shaping species assemblages or communities is challenging . A logical first-step, especially in poorly studied regions, is to apply historical methods to understand how past events have shaped current patterns . In this study we apply comparative phylogeography and biogeography methods [10, 11] to test the hypotheses that historical abiotic events are responsible for generating the pattern of similar levels of genetic divergence recently observed in two diverse taxa (three clades) that are codistributed in freshwater systems of Chile.
Hypothesized outcomes on cladogenic rates caused by major geological and climatic events documented in the region.
Predictions and expectations
Establishment of freshwater systems in Chile
~ 9.0 Ma
Origins and initial diversification
Basal nodes date to this period
Dramatic increase in sea-levels
4.3 - 5.3 Ma
Increased cladogenic events
Significant shifts in cladogenic rates (Yule two rate or rate-variable birth-death)
Dramatic decrease in sea-levels
1.4 - 2.2 Ma
Decreased cladogenic events
Significant shifts in cladogenic rates (Yule two rate or rate-variable birth-death)
Late Pleistocene climate oscillations
(Driver of speciation)
1.0 - 0.1 Ma
Increase in cladogenic events
Majority of cladogenic events date to this period
Significant shifts in cladogenic rates (Yule two rate or rate-variable birth-death)
Late Pleistocene climate oscillations
1.0 - 0.1 Ma
Decline in cladogenic rates
Few cladogenic events date to this period
Decline in cladogenic rates (density dependent model [logistic or exponential] of diversification)
No change in cladogenic rates
Constant cladogenesis models: yule pure-birth or pure-birth
Dates of cladogenic events
Estimates of the ages of major nodes in the Trichomycterus areolatus and Aegla phylogenies (Figure 2).
Effective sample sizes for node estimates were over 500 (as were most parameters); thus, we assume the results are from thoroughly sampled analyses. Models selected using AIC for Trichomycterus areolatus were GTR + I + R, HKY + I, and HKY for cyt-b, gh, and rag1 respectively. The GTR + I + R model provided the best fit for both 28S and mitochondrial loci for both of the Aegla datasets. Nuclear mutation rates estimated from the analyses are as follows: Trichomycterus areolatus (rag1 = 4.48 × 10-4 [2.34 × 10-4 - 6.79 × 10-4] and gh = 1.24 × 10-3 [7.41 × 10-4 - 1.81 × 10-3]), AeglaA (28S = 1.63 × 10-3 [7.25 × 10-4 - 2.58 × 10-3], AeglaB (28S = 1.04 × 10-3 [5.38 × 10-4 - 1.60 × 10-3].
Results of maximum likelihood tests of diversification models.
Best constant model
Best variable model
pure birth (AIC = 5.82)
DDL (AIC = 4.47)
1.36 (p = 0.21)
pure birth (AIC = 12.56)
DDL (AIC = 4.77)
7.79 (p = 0.02)
pure birth (AIC = 11.79)
rvbd (AIC = 6.32)
5.47 (p = 0.03)
This study was motivated by the observation of similar levels of genetic distances between lineages in three clades that are codistributed in freshwater systems of Chile. Our goal was to determine if this pattern was the result of codivergence driven by historical abiotic events. We discuss the temporal origins of the three groups. Then we discuss results of diversification tests and mechanisms that are postulated to have shaped this assemblage of freshwater taxa.
Lineage splitting time (mean) estimates of basal nodes vary considerably (2.18 - 4.50 Ma), suggesting different processes generated initial cladogenesis (Table 2 and Figure 2). The confidence intervals of the two crab groups overlap considerably, suggesting similar temporal origins. However, different models of cladogensis were recovered for these two groups suggesting different processes have shaped these two groups (see next section). The hypothesis that the original formation of river systems in the region some 9 Ma [12, 21] played a role in shaping this fauna is rejected. However the exceptionally high sea levels that occurred from 5.0 to 4.2 Ma are consistent with the origins of both Aegla clades (Table 2 and Figure 2; ). Cladogenesis within Trichomycterus areolatus occurred relatively recently (~2.25 Ma) and is consistent with a period of low sea levels. This time period is also consistent with the second set of cladogenic events observed within AeglaB. The majority of cladogenic events within T. areolatus appear in the Pleistocene with the last two events dating to a period (< 1 Ma) of high amplitude climate oscillations (Table 2 and Figure 2). Only three cladogenic events occurred during the last 1 Ma in the two crab clades. The relatively few cladogenic events observed during the last 1 Ma reject the hypothesis that the especially high amplitude climate cycles of this period played a significant role in driving cladogensis within these taxa .
Models and mechanisms of cladogenesis
Three different models of diversification were supported rejecting the hypothesis of codivergence among taxa. The pure birth model of diversification best explains cladogenesis within Trichomycterus areolatus and suggests that the rate initiated in this taxon was maintained throughout its evolutionary history and argues against any significant impact on rates by later sea-level changes or climate oscillations. A density-dependent model of declining diversification is the best explanation of cladogenesis in AeglaA.
Declines in diversification rates can occur for several reasons, but are most notable in adaptive radiations [27, 28]. Speciation rates can decline as available niches become filled with the accumulation of species  or if an underparameterized model of molecular evolution is employed . The molecular models employed in this study are relatively complex suggesting the rate decline is not an artifact of the analysis. Alternatively, the rate decline could be due to increased extinction rates.
To our knowledge, no empirical data exist on population sizes in Aegla, thus we can only make broad assumptions on the relative sizes of populations and their susceptibility to extinction. The distributions of several species of Chilean Aegla are restricted to single or several adjoining river system suggesting small populations ; Figure 1). Restricted and fragmented distributions coupled with small isolated populations in some species  suggest that population sizes are small, possibly increasing the risk of extinction . In addition, population sizes, inferred from estimates of genetic diversity, are low and declining in one crab species, Aegla alacalufi .
Alternatively, especially high amplitude climate cycles over the last one million years may have increased extinction. Only one lineage splitting event dates to the Late Pleistocene (< 1 Ma) rejecting the hypothesis that this period was a driver of speciation. Rather the decline in diversification suggests that the Late Pleistocene climate cycles, or the lowered sea levels of ~ 2 Ma, may have inhibited speciation or increased extinction within AeglaA. A phylogeographic study of Aegla alacalufi recovered deeper divergence between non-glaciated populations than between glaciated populations . Many of the drainage systems inhabited by Chilean Aegla species were covered at least partially by ice during the last and presumably previous glacial cycles [13, 16]. It is likely then that repeated glacial cycles coupled with small population sizes in many Aegla species inhibited diversification by increasing extinction rates.
The rate-variable model best explains diversification in AeglaB and suggests that this clade experienced a more dynamic evolutionary history than the other clades. That is the onset of cladogenesis initiated approximately 5 Ma were later altered by abiotic events. A second cluster of cladogenic events begins 1.44 Ma (node H) which is a time consistent with both lowered sea levels and the onset of high amplitude glacial cycles, suggesting that these these events increased diversification. Thus in contrast to AeglaA and Trichomycterus areolatus this period is postulated to be a driver of cladogenesis in AeglaB.
Tests of codivergence
Historical approaches have been used extensively in the study of codivergence in host-parasite systems . For example, the assumption that pattern equaled process in hantavirus was assumed because the virus phylogeny formed three primary clades that coincided with the three rodent subfamilies that formed their hosts and the fact that there is a high degree of host specificity in the system . However,  were able to reject the hypothesis of codivergence between the hantavirus virus and their mammalian hosts.
In this study historical methods were employed to determine if abiotic events were possibly responsible for driving diversification in these clades, which in turn allows us to test the hypothesis of codivergence. The observation of similar levels of genetic distance between major lineags among these taxa suggests that they are the result of either shared evolutionary or ecological processes. Our study illustrates the important role that history can play in shaping current patterns and argues against assuming that the observation of similar patterns among codistributed taxa equals a shared history . It remains to be determined how, if any, ecological interactions have shaped this species assemblage.
Employing phylogeographic and phylogenetic methods to test hypotheses of codivergence among contemporary faunas is an exciting and developing research area . A fundamental challenge in this endeavor will be choosing methods that can tease apart complex evolutionary histories among taxa. In this paper we employed two different methods to test the hypothesis of codivergence among our three focal groups. Obviously there are some caveats with these methods. For example, confidence intervals on divergence dates are wide and in some cases overlap multiple abiotic events and divergence times from the other taxa. Relying on this method alone would not have allowed us to test our hypothesis. However, differences in timing of lineage splitting combined with the fact that we observed three taxon specific models of cladogenesis allows us to reject the hypothesis of codivergence. Another potential problem arises by using mutation rates that are not taxon specific. This is an issue for most comparative molecular studies due to the paucity of taxon specific rates. Because mutation rates vary widely across taxa - vertebrates versus invertebrates in our study - assuming a common rate is an unacceptable solution. Finally we recognize our relatively taxon poor phylogenies may make it difficult to rigorously compare models of diversification [27, 35]. For example missing or undiscovered taxa could cause an increase, decrease or shift in diversification rates. The previous phylogeographic and phylogenetic studies of these two taxa did recover undescribed independently-evolving lineages and these were included in this study. Despite these caveats our study illustrates the potential to test codivergence among codistributed fauna when multiple approaches are employed [34, 36–39]
Results from two different methodological approaches reject the hypothesis of codivergence among our three focal groups. Rather these representatives of the Chilean freshwater assemblage evolved through temporally-independent processes. Although this freshwater fauna is relatively depauperate additional taxa need to be studied to determine if this is a pervasive observation in freshwater systems of Chile and how ecological interactions are currently operating on the entire species assemblage. This study illustrates the importance of using historical methods to test hypotheses on how contemporary patterns have been formed and cautions against assuming that congruent patterns are the result of shared processes or ecological interactions alone.
Materials and methods
Chilean taxa and molecular data
The catfish Trichomycterus areolatus and 21 species of Aegla crabs were studied from the southern temperate freshwater river systems of Chile (Figure 1). Trichomycterus areolatus is an abundant and widespread freshwater fish with broad ecological tolerances . Data for Trichomycterus areolatus are the mitochondrial cytochrome b gene (1137 base-pairs, [bp] GenBank accession numbers FJ772091 - FJ772237) and two newly obtained nuclear loci (growth hormone [1491 bp: JN186609 - JN186714] and rag1 [992 bp: JN186409 - JN186531]) from 63 individuals comprising seven clades . Like the catfish, Chilean Aegla crabs are widely distributed and occur in a variety of freshwater habitats . Two non-sister clades of mostly Chilean Aegla species were examined [19, 23]. The first clade, hereafter referred to as AeglaA, includes nine taxa, the second, AeglaB, consists of twelve taxa . Molecular data for Aegla consisted of the cytochrome oxidase I and II genes (1598 bp: AY049985-050166) and the nuclear 28S gene (3105 bp: AY595931 - AY595973, AY596092 - AY596093), and 20 and 24 individuals from AeglaA and AeglaB, respectively [19, 20, 23].
Phylogeny and dating
The best-fitting model of molecular evolution for each locus was determined using the Akaike Information Criterion (AIC) in jModeltest . A species tree and dates of major nodes were estimated simultaneously using *BEAST  as implemented in Beast 1.6.1 . A Yule tree prior, piecewise linear and constant-root population size model and lognormal molecular clock were assumed. Analyses were run 5.0 × 107 generations with a burn-in of 10% . We assumed a mitochondrial mutation rate of 1.0 × 10-8 , and a range of 8.3 × 10-9- 1.3 × 10-8 [45, 46] for Trichomycterus areolatus and Aegla, respectively. Mutation rates for all nuclear markers were estimated, relative to the mitochondrial rate, using a normal prior. Due to uncertainty in the appropriate outgroup a midpoint root was used to root the Trichomycterus areolatus phylogeny . Two individuals of Aegla papudo (cytochrome I: AY050077 - AY050078; cytochrome II: AY050121-AY050121; 28S: AY595929 - AY595930) were used to root the Aegla phylogeny [19, 20, 23]. Mean divergence estimates in millions of years and their 95% confidence intervals were obtained for all nodes using TRACER v1.4 . If dates of major nodes or their confidence intervals overlap with abiotic events (i.e., sea level changes; Table 1), we assume those events may have been responsible for driving divergence at those nodes .
Tests of diversification models
The hypothesis of codivergence predicts the same process of cladogenesis should be observed across taxa. To this end we evaluate the fit of six explicit models of diversification to each of the three phylogenies using a maximum likelihood test implemented in LASER . In general, these models provide the means to determine if cladogenic rates were altered (i.e. by abiotic events), or remained constant (no effect), during the phylogenetic history of each group. If these taxa codiverged, due to shared historical events, we expect each to have evolved under the same model. The six models compared were: Yule pure-birth, birth-death, logistic density-dependent, exponential density-dependent, Yule two-rate, and a rate-variable birth-death model. These models capture expected patterns of cladogenesis under various historical scenarios (Table 1). For example, if changes in Pliocene sea levels altered rates, we expect to observe an instantaneous shift in rates. Similarly, we would observe a shift in rates if Late Pleistocene climate change increased lineage splitting. Under these two scenarios, we expect to observe either a Yule two-rate or rate-variable birth-death model. In contrast, if Late Pleistocene climate change decreased diversification, through increased extinction or decreased speciation, we would expect to detect either a logistic or exponential density dependent model. Finally, constant models (Yule pure-birth or birth-death) are the expected outcome if climate and sea level changes had no impact on speciation rates (pure-birth and birth-death models assume no extinction has occurred). The best constant (RC) and best rate variable (RV) model were determined using (AIC). Delta AICrc (= AICRC - AICRV) was used as a test statistic. Significance of each delta AICrc was determined using null distributions derived from 5,000 trees simulated under a constant birth-death (extinction = lineage splitting) model (two-tailed alpha = 0.05, ).
Our objective in this study was to investigate processes of cladogenesis in codistributed taxa with similar genetic distances between evolutionarily distinct lineages. Unmack et al  phylogeographic study of Trichomycterus areolatus recovered seven well-supported terminal clades with little to no geographical overlap between them. These clades had model-corrected sequence divergences between them greater than many recognized species (> 2%: e.g. [11, 50]). These observations qualify the seven clades as evolutionarily distinct lineages under several concepts including evolutionarily significant units (ESU: ) and the phylogenetic species concept . Therefore, we recognized each of the seven clades recovered in Trichomycterus areolatus as taxonomically equivalent to the currently recognized species as well as several recently discovered lineages in Aegla.
This work was funded by a grant from the U.S. National Science Foundation (NSF) PIRE program (OISE 0530267) to support collaborative research on Patagonian Biodiversity among the following institutions (listed alphabetically): Brigham Young University, Centro Nacional Patagónico (AR), Dalhousie University, Instituto Botánico Darwinion (AR), Universidad Austral de Chile, Universidad Nacional del Comahue, Universidad de Concepción, and George Washington University.
- MacArthur RH, Wilson EO: The theory of island biogeography. 1967, Princeton, N.J.,: Princeton University PressGoogle Scholar
- Hubbell SP: The unified neutral theory of biodiversity and biogeography. 2001, Princeton [N.J.]: Princeton University PressGoogle Scholar
- Wiens JJ, Donoghue MJ: Historical biogeography, ecology and species richness. Trends Ecol Evol. 2004, 19 (12): 639-644. 10.1016/j.tree.2004.09.011.View ArticlePubMedGoogle Scholar
- Lee JB, Johnson JB: Biogeography of the livebearing fish Poecilia gillii in Costa Rica: are phylogeographical breaks congruent with fish community boundaries?. Mol Ecol. 2009, 18 (19): 4088-4101. 10.1111/j.1365-294X.2009.04339.x.View ArticlePubMedGoogle Scholar
- Burridge CP, Craw D, Waters JM: River capture, range expansion, and cladogenesis: the signature of freshwater vicariance. Evolution. 2006, 60 (5): 1038-1049.View ArticlePubMedGoogle Scholar
- Hewitt G: The genetic legacy of the Quaternary ice ages. Nature. 2000, 405 (6789): 907-913. 10.1038/35016000.View ArticlePubMedGoogle Scholar
- Rosen DE: Vicariant Patterns and Historical Explanation in Biogeography. Syst Biol. 1978, 27 (2): 159-188.Google Scholar
- Cavender-Bares J, Kozak KH, Fine PVA, Kembel SW: The merging of community ecology and phylogenetic biology. Ecol Lett. 2009, 12 (7): 693-715. 10.1111/j.1461-0248.2009.01314.x.View ArticlePubMedGoogle Scholar
- Johnson MTJ, Stinchcombe JR: An emerging synthesis between community ecology and evolutionary biology. Trends Ecol Evol. 2007, 22 (5): 250-257. 10.1016/j.tree.2007.01.014.View ArticlePubMedGoogle Scholar
- Nelson GJ, Platnick NI: Systematics and biogeography: cladistics and vicariance. 1981, New York: Columbia University PressGoogle Scholar
- Avise JC: Phylogeography: the history and formation of species. 2000, Cambridge, Mass.: Harvard University PressGoogle Scholar
- Potter PE: The Mesozoic and Cenozoic paleodrainage of South America: a natural history. Journal of South American Earth Sciences. 1997, 10 (5-6): 331-344. 10.1016/S0895-9811(97)00031-X.View ArticleGoogle Scholar
- Unmack PJ, Bennin AP, Habit EM, Victoriano PF, Johnson JB: Impact of ocean barriers, topography, and glaciation on the phylogeography of the catfish Trichomycterus areolatus (Teleostei: Trichomycteridae) in Chile. Biol J Linn Soc. 2009, 97 (4): 876-892. 10.1111/j.1095-8312.2009.01224.x.View ArticleGoogle Scholar
- Beheregaray LB: Twenty years of phylogeography: the state of the field and the challenges for the Southern Hemisphere. Mol Ecol. 2008, 17 (17): 3754-3774.PubMedGoogle Scholar
- Ruzzante DE, Walde SJ, Cussac VE, Dalebout ML, Seibert J, Ortubay S, Habit E: Phylogeography of the Percichthyidae (Pisces) in Patagonia: roles of orogeny, glaciation, and volcanism. Mol Ecol. 2006, 15 (10): 2949-2968. 10.1111/j.1365-294X.2006.03010.x.View ArticlePubMedGoogle Scholar
- Xu J, Perez-Losada M, Jara CG, Crandall KA: Pleistocene glaciation leaves deep signature on the freshwater crab Aegla alacalufi in Chilean Patagonia. Mol Ecol. 2009Google Scholar
- Balian EV, Lévêque C, Segers H, Martens K: Freshwater Animal Diversity Assessment. 2008, 198:View ArticleGoogle Scholar
- Dyer B: Systematic review and biogeography of the freshwater fishes of Chile. Estudios Oceanológicos. 2000, 19: 77-98.Google Scholar
- Perez-Losada M, Bond-Buckup G, Jara CG, Crandall KA: Conservation Assessment of Southern South American Freshwater Ecoregions on the Basis of the Distribution and Genetic Diversity of Crabs from the Genus Aegla. Conserv Biol. 2009Google Scholar
- Pérez-Losada M, Jara CG, Bond-Buckup G, Crandall KA: Phylogenetic relationships among the species of Aegla (Anomura: Aeglidae) freshwater crabs from Chile. Journal of Crustacean Biology. 2002, 22 (2): 304-313. 10.1651/0278-0372(2002)022[0304:PRATSO]2.0.CO;2.View ArticleGoogle Scholar
- Lundberg JG, Marshall LG, J G, Horton B, L MMCS, F W: The stage for Neotropical fish diversification: A history of tropical South American river. 1998, Porto Alegre, Brazil: EDIPUCRSGoogle Scholar
- Perez-Losada M, Bond-Buckup G, Jara CG, Crandall KA: Molecular systematics and biogeography of the southern South american freshwater "crabs" Aegla (decapoda: Anomura: Aeglidae) using multiple heuristic tree search approaches. Syst Biol. 2004, 53 (5): 767-780. 10.1080/10635150490522331.View ArticlePubMedGoogle Scholar
- Pérez-Losada M, Jara CG, Bond-Buckup G, Crandall KA: Conservation phylogenetics of Chilean freshwater crabs Aegla (Anomura, Aeglidae): assigning priorities for aquatic habitat protection. Biol Conserv. 2002, 105 (3): 345-353. 10.1016/S0006-3207(01)00218-X.View ArticleGoogle Scholar
- Le Roux JP, Olivares DM, Nielsen SN, Smith ND, Middleton H, Fenner J, Ishman SE: Bay sedimentation as controlled by regional crustal behaviour, local tectonics and eustatic sea-level changes: Coquimbo Formation (Miocene-Pliocene), Bay of Tongoy, central Chile. Sedimentary Geology. 2006, 184 (1-2): 133-153. 10.1016/j.sedgeo.2005.09.023.View ArticleGoogle Scholar
- Lisiecki LE, Raymo ME: A Pliocene-Pleistocene stack of 57 globally distributed benthic δ18Ο records. Paleoceanography. 2005, 20 (1): 1-17.View ArticleGoogle Scholar
- Zink RM, Klicka J, Barber BR: The tempo of avian diversification during the Quaternary. Philos T Roy Soc B. 2004, 359 (1442): 215-219. 10.1098/rstb.2003.1392.View ArticleGoogle Scholar
- Phillimore AB, Price TD: Density-dependent cladogenesis in birds. Plos Biol. 2008, 6 (3): e71-10.1371/journal.pbio.0060071.View ArticlePubMedPubMed CentralGoogle Scholar
- Schluter D: The ecology of adaptive radiation/Dolph Schluter. 2000, Oxford: Oxford University PressGoogle Scholar
- McPeek MA: The ecological dynamics of clade diversification and community assembly. Am Nat. 2008, 172 (6): E270-284. 10.1086/593137.View ArticlePubMedGoogle Scholar
- Revell LJ, Harmon LJ, Glor RE: Underparameterized model of sequence evolution leads to bias in the estimation of diversification rates from molecular phylogenies. Syst Biol. 2005, 54 (6): 973-983. 10.1080/10635150500354647.View ArticlePubMedGoogle Scholar
- Newman D, Pilson D: Increased probability of extinction due to decreased genetic effective population size: experimental populations of Clarkia Pulchella. 1997, Hoboken, NJ, ETATS-UNIS: Wiley, 51:Google Scholar
- Page RDM, Charleston MA: Trees within trees: phylogeny and historical associations. Trends Ecol Evol. 1998, 13 (9): 356-359. 10.1016/S0169-5347(98)01438-4.View ArticlePubMedGoogle Scholar
- Ramsden C, Holmes EC, Charleston MA: Hantavirus Evolution in Relation to Its Rodent and Insectivore Hosts: No Evidence for Codivergence. Mol Biol Evol. 2009, 26 (1): 143-153.View ArticlePubMedGoogle Scholar
- Hickerson MJ, Carstens BC, Cavender-Bares J, Crandall KA, Graham CH, Johnson JB, Rissler L, Victoriano PF, Yoder AD: Phylogeography's past, present, and future: 10 years after Avise, 2000. Mol Phylogenet Evol. 2010, 54 (1): 291-301. 10.1016/j.ympev.2009.09.016.View ArticlePubMedGoogle Scholar
- Ricklefs RE: Estimating diversification rates from phylogenetic information. Trends Ecol Evol. 2007, 22 (11): 601-610. 10.1016/j.tree.2007.06.013.View ArticlePubMedGoogle Scholar
- Barber BR, Klicka J: Two pulses of diversification across the Isthmus of Tehuantepec in a montane Mexican bird fauna. P Roy Soc B-Biol Sci. 2010, 277 (1694): 2675-2681. 10.1098/rspb.2010.0343.View ArticleGoogle Scholar
- Carnaval AC, Hickerson MJ, Haddad CF, Rodrigues MT, Moritz C: Stability predicts genetic diversity in the Brazilian Atlantic forest hotspot. Science. 2009, 323 (5915): 785-789. 10.1126/science.1166955.View ArticlePubMedGoogle Scholar
- Templeton AR: Out of Africa again and again. Nature. 2002, 416 (6876): 45-51. 10.1038/416045a.View ArticlePubMedGoogle Scholar
- Templeton AR: Statistical hypothesis testing in intraspecific phylogeography: nested clade phylogeographical analysis vs. approximate Bayesian computation. Mol Ecol. 2009, 18 (2): 319-331. 10.1111/j.1365-294X.2008.04026.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Posada D: jModelTest: Phylogenetic Model Averaging. Mol Biol Evol. 2008, 25 (7): 1253-1256. 10.1093/molbev/msn083.View ArticlePubMedGoogle Scholar
- Heled J, Drummond AJ: Bayesian inference of species trees from multilocus data. Mol Biol Evol. 27 (3): 570-580.
- Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. Bmc Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.View ArticlePubMedPubMed CentralGoogle Scholar
- Drummond AJ, Ho SY, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. Plos Biol. 2006, 4 (5): e88-10.1371/journal.pbio.0040088.View ArticlePubMedPubMed CentralGoogle Scholar
- Near TJ, Benard MF: Rapid allopatric speciation in logperch darters (Percidae: Percina). Evolution. 2004, 58 (12): 2798-2808.View ArticlePubMedGoogle Scholar
- Knowlton N, Weigt LA, Solorzano LA, Mills DK, Bermingham E: Divergence in proteins, mitochondrial DNA, and reproductive compatibility across the isthmus of Panama. Science. 1993, 260 (5114): 1629-1632. 10.1126/science.8503007.View ArticlePubMedGoogle Scholar
- Schubart CD, Diesel R, Hedges SB: Rapid evolution to terrestrial life in Jamaican crabs. 1998, London, ROYAUME-UNI: Nature Publishing Group, 393:Google Scholar
- Rambaut A, Drummond AJ: TRACER v1.4. 2007Google Scholar
- Lemmon EM, Lemmon AR, Cannatella DC, Brown J: Geological and climatic forces driving speciation in the continentally distributed trilling chorus frogs (Pseudacris). Evolution. 2007, 61 (9): 2086-2103. 10.1111/j.1558-5646.2007.00181.x.View ArticlePubMedGoogle Scholar
- Rabosky DL: Likelihood methods for detecting temporal shifts in diversification rates. Evolution. 2006, 60 (6): 1152-1164.View ArticlePubMedGoogle Scholar
- April J, Mayden RL, Hanner RH, Bernatchez L: Genetic calibration of species diversity among North America's freshwater fishes. Proceedings of the National Academy of Sciences. 2011, 108 (26): 10602-10607. 10.1073/pnas.1016437108.View ArticleGoogle Scholar
- Moritz C: Defining 'Evolutionarily Significant Units' for conservation. Trends Ecol Evol. 1994, 9 (10): 373-375. 10.1016/0169-5347(94)90057-4.View ArticlePubMedGoogle Scholar
- Cracraft J: Speciation and its ontology: the empirical consequences of alternative species concepts for understanding patterns and processes of differentiation. Speciation and its Consequences. Edited by: Endler DOJA. 1989, Sinauer Associates, 28-59.Google Scholar
- Levinton JS: A theory of diversity equilibrium and morphological evolution. Science. 1979, 204 (4390): 335-336. 10.1126/science.204.4390.335.View ArticlePubMedGoogle 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 cited.