Comparative analysis reveals that polyploidy does not decelerate diversification in fish

While the proliferation of the species‐rich teleost fish has been ascribed to an ancient genome duplication event at the base of this group, the broader impact of polyploidy on fish evolution and diversification remains poorly understood. Here, we investigate the association between polyploidy and diversification in several fish lineages: the sturgeons (Acipenseridae: Acipenseriformes), the botiid loaches (Botiidae: Cypriniformes), Cyprininae fishes (Cyprinidae: Cypriniformes) and the salmonids (Salmonidae: Salmoniformes). Using likelihood‐based evolutionary methodologies, we co‐estimate speciation and extinction rates associated with polyploid vs. diploid fish lineages. Family‐level analysis of Acipenseridae and Botiidae revealed no significant difference in diversification rates between polyploid and diploid relatives, while analysis of the subfamily Cyprininae revealed higher polyploid diversification. Additionally, order‐level analysis of the polyploid Salmoniformes and its diploid sister clade, the Esociformes, did not support a significantly different net diversification rate between the two groups. Taken together, our results suggest that polyploidy is generally not associated with decreased diversification in fish – a pattern that stands in contrast to that previously observed in plants. While there are notable differences in the time frame examined in the two studies, our results suggest that polyploidy is associated with different diversification patterns in these two major branches of the eukaryote tree of life.     

Molecular evidence for the monophyly of East Asian groups of Cyprinidae (Teleostei: Cypriniformes) derived from the nuclear recombination activating gene 2 sequences

The family Cyprinidae is one of the largest families of fishes in the world and a well-known component of the East Asian freshwater fish fauna. However, the phylogenetic relationships among cyprinids are still poorly understood despite much effort paid on the cyprinid molecular phylogenetics. Original nucleotide sequence data of the nuclear recombination activating gene 2 were collected from 109 cyprinid species and four non-cyprinid cypriniform outgroup taxa and used to infer the cyprinid phylogenetic relationships and to estimate node divergence times. Phylogenetic reconstructions using maximum parsimony, maximum likelihood, and Bayesian analysis retrieved the same clades, only branching order within these clades varied slightly between trees. Although the morphological diversity is remarkable, the endemic cyprinid taxa in East Asia emerged as a monophyletic clade referred to as Xenocypridini. The monophyly for the subfamilies including Cyprininae and Leuciscinae, as well as the tribes including Labeonini, Gobionini, Acheilognathini, and Leuciscini, was also well resolved with high nodal support. Analysis of the RAG2 gene supported the following cyprinid molecular phylogeny: the Danioninae is the most basal subfamily within the family Cyprinidae and the Cyprininae is the sister group of the Leuciscinae. The divergence times were estimated for the nodes corresponding to the principal clades within the Cyprinidae. The family Cyprinidae appears to have originated in the mid-Eocene in Asia, with the cladogenic event of the key basal group Danioninae occurring in the early Oligocene (about 31-30 MYA), and the origins of the two subfamilies, Cyprininae and Leuciscinae, occurring in the mid-Oligocene (around 26 MYA).     

Ecological and geological processes impacting speciation modes drive the formation of wide-range disjunctions within tribe Putorieae (Rubiaceae)

Wide-range geographically discontinuous distributions have long intrigued scientists. We explore the role of ecology, geology, and dispersal in the formation of these large-scale disjunctions, using the angiosperm tribe Putorieae (Rubiaceae) as a case study. From DNA sequences of nuclear ITS and six plastid markers, we inferred a phylogeny with 65% of all known Putorieae species. Divergence times, ancestral ranges, and diversification rate shifts were then estimated using Bayesian inference. We further explored species climatic tolerances and performed ancestral niche reconstruction to discriminate among alternative speciation modes, including geographical and ecological vicariance, and ecogeographical, ecological, and dispersal-mediated speciation. As a result, we identified seven major clades in Putorieae, some of which exhibit striking geographical disjunctions, matching the Rand Flora pattern, with sister species in the Canary Islands andeastern and southern Africa. Initial diversification within the tribe occurred in the early Miocene, coincident with a period of climate warming; however, most clades diverged within the last 10 Myr. Aridification and high extinction rates, coupled with ecological vicariance, explain the oldest disjunctions. Adaptation to new environmental conditions, after allopatry, is observed in several clades. Dispersal, either long-distance or via corridors made available by mountain uplift, is behind the most recent disjunctions. Some of these events were followed by ecological speciation and rapid diversification, with species becoming adapted to xeric or increasingly colder continental climates. We show that an integrative approach may help discriminate among speciation modes invoked to explain disjunctions at macroevolutionary time scales, even when extinction has erased the signature of past events.     

The legacy of a vanished sea: a high level of diversification within a European freshwater amphipod species complex driven by 15 My of Paratethys regression

The formation of continental Europe in the Neogene was due to the regression of the Tethys Ocean and of the Paratethys Sea. The dynamic geology of the area and repetitious transitions between marine and freshwater conditions presented opportunities for the colonization of newly emerging hydrological networks and diversification of aquatic biota. Implementing mitochondrial and nuclear markers in conjunction with a large‐scale sampling strategy, we investigated the impact of this spatiotemporal framework on the evolutionary history of a freshwater crustacean morphospecies. The Gammarus balcanicus species complex is widely distributed in the area previously occupied by the Paratethys Sea. Our results revealed its high diversification and polyphyly in relation to a number of other morphospecies. The distribution of the studied amphipod is generally characterized by very high local endemism and divergence. The Bayesian time‐calibrated reconstruction of phylogeny and geographical distribution of ancestral nodes indicates that this species complex started to diversify in the Early Miocene in the central Balkans, partially in the shallow epicontinental sea. It is possible that there were several episodes of inland water colonization by local brackish water lineages. Subsequent diversification within clades and spread to new areas could have been induced by Alpine orogeny in the Miocene/Pliocene and, finally, by Pleistocene glaciations. The present distribution of clades, in many cases, still reflects Miocene palaeogeography of the area. Our results point out that investigations of the historical aspect of cryptic diversity in other taxa may help in a general understanding of the origins of freshwater invertebrate fauna of Europe.     

The evolutionary history of Melianthus (Melianthaceae)

The evolutionary origins of the morphological and taxonomic diversity of angiosperms is poorly known. We used the genus Melianthus to explore the diversification of the southern African flora. Melianthus comprises eight species, and a phylogeny based on one nuclear and two plastid genes, as well as a morphological data set, confirmed that the genus is monophyletic. The two earliest diverging lineages are found in relatively mesic habitats, whereas the two terminal clades (an eastern and a western clade), each with three species, favor more arid habitats. The eastern clade is largely restricted to the summer-rainfall parts of southern Africa, and the western clade is found in winter-rainfall region. Molecular dating indicates a mid-Tertiary origin of the genus, with diversification of the eastern and western clades coincident with the Late Miocene-Pliocene uplift of the Escarpment mountains and the establishment of summer aridity along the west coast. The remarkably complex flowers are indicative of sunbird pollination, but many smaller birds can also visit. Speciation may be the consequence of allopatric divergence into edaphic-climatic niches. Divergence in flower and inflorescence morphology might be in response to the divergent pressures for nectar conservation in arid regions coupled with the need for signaling to avian pollinators in generally shrubby vegetation.     

Molecular phylogeny of congeneric monogenean parasites (Dactylogyrus): a case of intrahost speciation

Dactylogyrus species (Dactylogyridae: Monogenea) are a group of monogenean gill parasites that are highly specific to freshwater fish of the family Cyprinidae. Dactylogyrus species were sampled from 19 cyprinids and one percid collected in Europe. Using partial 18S rDNA and ITS1 sequences, a phylogeny of 51 Dactylogyrus species was reconstructed to investigate the patterns of parasite speciation and diversification. Three main Dactylogyrus lineages were recognized from all phylogenetic trees, that is, analysis of 18S rDNA alone and combined 18SrDNA and ITS1. The first lineage associates the Dactylogyrus species of Cyprinus carpio and Carassius auratus of the Cyprininae; the second associates Dactylogyrus species of the Gobioninae, Pseudorasbora parva of the Rasborinae, and Ctenopharyngodon idella of the Cyprininae; and the third associates Dactylogyrus species of the Leuciscinae and Alburninae and Barbus barbus of the Cyprininae. Our results suggest that the genus Dactylogyrus is of quite recent origin and that these three lineages separated from each other in a very short period of time. Host subfamily mapping onto the parasite tree inferred from analysis of the combined dataset showed that the Cyprininae could be plesiomorphic hosts for Dactylogyrus. Dactylogyrus parasites would have secondarily colonized the Percidae and representatives of the Leuciscinae, Alburninae, Gobioninae, and Rasborinae. Comparison of host and parasite phylogenetic relationships indicated that a very high number of parasite duplications occurred within two of the three Dactylogyrus lineages. Dactylogyrus diversification can be mainly explained by sympatric intrahost speciation events that seem to be correlated to strict host specificity. Moreover, the present study shows that the congeneric parasites speciating within one host tend to occupy niches within hosts differing at least in one niche parameter.     

Ecological opportunity alters the timing and shape of adaptive radiation

The uneven distribution of diversity is a conspicuous phenomenon across the tree of life. Ecological opportunity is a prominent catalyst of adaptive radiation and therefore may alter patterns of diversification. We evaluated the distribution of shifts in diversification rates across the cichlid phylogeny and the distribution of major clades across phylogenetic space. We also tested if ecological opportunity influenced these patterns. Colonization‐associated ecological opportunity altered the tempo and mode of diversification during the adaptive radiation of cichlid fishes. Clades that arose following colonization events diversified faster than other clades. Speciation rate shifts were nonrandomly distributed across the phylogeny such that they were disproportionally concentrated around nodes that corresponded with colonization events (i.e., of continents, river basins, or lakes). Young clades tend to expand faster than older clades; however, colonization‐associated ecological opportunity accentuated this pattern. There was an interaction between clade age and ecological opportunity that explained the trajectory of clades through phylogenetic space over time. Our results indicate that ecological opportunities afforded by continental and ecosystem‐scale colonization events explain the dramatic speciation rate heterogeneity and phylogenetic imbalance that arose during the evolutionary history of cichlid fishes.     

Diversification in the Andes: age and origins of South American Heliotropium lineages (Heliotropiaceae, Boraginales)

The uplift of the Andes was a major factor for plant diversification in South America and had significant effects on the climatic patterns at the continental scale. It was crucial for the formation of the arid environments in south-eastern and western South America. However, both the timing of the major stages of the Andean uplift and the onset of aridity in western South America remain controversial. In this paper we examine the hypothesis that the Andean South American groups of Heliotropium originated and diversified in response to Andean orogeny during the late Miocene and a the subsequent development of aridity. To this end, we estimate divergence times and likely biogeographical origins of the major clades in the phylogeny of Heliotropium, using both Bayesian and likelihood methods. Divergence times of all Andean clades in Heliotropium are estimated to be of late Miocene or Pliocene ages. At least three independent Andean diversification events can be recognized within Heliotropium. Timing of the diversification in the Andean lineages Heliotropium sects.Heliothamnus, Cochranea, Heliotrophytum, Hypsogenia, Plagiomeris, Platygyne clearly correspond to a rapid, late Miocene uplift of the Andes and a Pliocene development of arid environments in South America.     

Biogeographic patterns of polyploid species for the angiosperm flora in China

Exploring the frequency and distribution pattern of polyploid species in geographic parameters is of significance in understanding the mechanisms of polyploid speciation and evolutionary drivers of biodiversity. We here explored polyploid and paleopolyploid incidence frequency in a scale of 100 × 100 km grids in China. We found 33% of angiosperm species are polyploidy in China, and 23% of polyploid speciation. Western China and eastern China showed a significantly different polyploid and paleopolyploid frequency, with an evolutionary cradle of polyploid angiosperms in the Qinghai–Tibetan plateau. Herbaceous species exhibited higher polyploid frequency but lower paleopolyploid frequency than woody species, indicating the former experienced more rapid differentiation and speciation than the latter. Our results indicate that western China is an evolutionary cradle for polyploid angiosperms where harsh environment facilitates the establishment and survival of polyploids, while polyploid lineages tend to rediploidize to be diploids with sufficient time in suitable environment.     

Population genetic structure and geographic differentiation in the hot spring snake Thermophis baileyi (Serpentes, Colubridae): indications for glacial refuges in southern-central Tibet

The endangered hot spring snake (Thermophis baileyi) is endemic to the Tibetan Plateau and restricted to a few sites at high altitudes (>3500 m a.s.l.). Its strong preference for habitats with hot springs might be an early adaptation to the cooling climate during the plateau uplift. Some of these thermal sites may have been free of ice during the last glacial maximum (LGM) serving as refuges for the snake. To test this and to examine the population genetic structure data of 12 microsatellites and mitochondrial DNA were obtained from 153 individuals from 12 sites across the plateau. Four of the microsatellite loci were newly developed for this study. The results suggest that T. baileyi has at least two genetically diverse clades in Tibet, which developed during the Pleistocene and expanded after the LGM. Accordingly, the existence of separate glacial refuges on the central plateau can be assumed. Analyses of the genetic variation indicated a high level of geographic differentiation and population structure on a regional as well as on a range-wide scale. The study shows that, apart from the phylogeographic signatures, the diversification of current Thermophis populations is caused by (a) a limited dispersal due to mountain ranges, (b) a strong preference for hot springs and (c) the insular distribution of suitable habitats on the plateau. The present data should be taken into account in necessary monitoring and management plans for T. baileyi.     

Phylogenetic relationships in genus Niviventer (Rodentia: Muridae) in China inferred from complete mitochondrial cytochrome b gene

Chinese species of the genus Niviventer, predominantly distributed in the southeastern Tibetan Plateau and in Taiwan, are a diverse group and have not yet received a thorough molecular phylogenetic analysis. Here, we reconstructed the phylogenetic relationships of 32 specimens representing nine Chinese species of Niviventer, based on sequences of the complete mitochondrial cytochrome b gene. Maximum parsimony, maximum likelihood and Bayesian analysis resulted in three consistent trees, each supported by high bootstrap values. The results showed that the Niviventer species included here are monophyletic. The nine species were classified into three distinct clades: clade A with Niviventer brahma, N. confucianus, N. coxingi, N. culturatus, N. eha and N. fulvescens; clade B with N. andersoni and N. excelsior; clade C with N. cremoriventer. Our results also suggested that N. culturatus should be a valid species rather than a subspecies of N. confucianus. Divergence times among species were calibrated according to the middle-late Pleistocene (1.2-0.13 Mya) fossil records of N. confucianus. The results demonstrated that the first radiation event of the genus Niviventer occurred in early Pleistocene (about 1.66 Mya), followed by the divergence of clades A and B at about 1.46 Mya. Most of the extant Niviventer species appeared during early to middle Pleistocene (about 1.29-0.67 Mya). These divergence times are coincidental with the last uplift events of the Tibetan Plateau, Kun-Huang movement, Pleistocene glaciations and the vicariant formation of Taiwan Strait. Consequently geographical events and Pleistocene glaciations have played a great role in the diversification of Niviventer.     

Low-copy nuclear DNA sequences reveal a predominance of allopolyploids in a New Zealand Asplenium fern complex

Recent generalisations about polyploidy in plants have been largely based on studies of angiosperms. A compelling group to compare with angiosperms is ferns, because of their high polyploidy. The bi-parental inheritance of nuclear DNA sequence markers makes them advantageous for investigating polyploid complexes, but few such markers have been available for ferns. We have used DNA sequences from the low-copy nuclear LFY locus to study an Asplenium polyploid complex. The New Zealand species of this Austral group comprise seven tetraploids and eight octoploids. LFY sequences indicate that allopolyploidy is much more predominant than previously thought, being implicated in the origins of seven of the octoploids. One of the tetraploids has had a central role, being a progenitor for five of the octoploids. All of the octoploids appear to have relatively recent origins, with the dynamic environmental conditions of the Pleistocene possibly playing a role in their formation and/or establishment.     

How important is it to consider lineage diversification heterogeneity in macroevolutionary studies? Lessons from the lizard family Liolaemidae

Macroevolutionary and biogeographical studies commonly apply multiple models to test state-dependent diversification. These models track the association between states of interest along a phylogeny, although many of them do not consider whether different clades might be evolving under different evolutionary drivers. Yet, they are still commonly applied to empirical studies without careful consideration of possible lineage diversification heterogeneity along the phylogenetic tree. A recent biogeographic study has suggested that orogenic uplift of the southern Andes has acted as a species pump, driving diversification of the lizard family Liolaemidae (307 described species), native to temperate southern South America. Here, we argue against the Andean uplift as main driver of evolution in this group. We show that there is a clear pattern of heterogeneous diversification in the Liolaemidae, which biases state- and environment-dependent analyses in, respectively, the GeoSSE and RPANDA programs. We show here that there are two shifts to accelerated speciation rates involving two clades that have both been classified as having “Andean” distributions. We incorporated the Geographic Hidden-State Speciation and Extinction model (GeoHiSSE) to accommodate unrelated diversification shifts, and also re-analyzed the data in RPANDA program after splitting biologically distinct clades for separate analyses, as well as including a more appropriate set of models. We demonstrate that the “Andean uplift” hypothesis is not supported when the heterogeneous diversification histories among these lizards is considered. We use the Liolaemidae as an ideal system to demonstrate potential risks of ignoring clade-specific differences in diversification patterns in macroevolutionary studies. We also implemented simulations to show that, in agreement with previous findings, the HiSSE approach can effectively and substantially reduce the level of distribution-dependent models receiving the highest AIC weights in such scenarios. However, we still find a relatively high rate (15%) of distribution-dependent models receiving the highest AIC weights, and provide recommendations related to the set of models included in the analyses that reduce these rates by half. Finally, we demonstrate that trees including clades following different dependent-drivers affect RPANDA analyses by producing different outcomes, ranging from partially correct models to completely misleading results. We provide recommendations for the implementation of both programs.     

Origins and introgression of polyploid species in Mentzelia section Trachyphytum (Loasaceae)

? Premise of the study: Polyploid speciation has been important in plant evolution. However, the conditions that favor the origination and persistence of polyploids are still not well understood. Here, we examine origins of 16 polyploid species in Mentzelia section Trachyphytum. ? Methods: We used phylogeny reconstructions based on DNA sequences from plastid regions and the nuclear gene isocitrate dehydrogenase (idh) to construct hypotheses of introgression and polyploidization. ? Key results: Molecular data suggest that homoploid hybridization has been surprisingly common in Trachyphytum. Diploid species had unequal involvement in polyploid origins, but most polyploid taxa had allopolyploid origins from extant progenitors. A few polyploids with extreme phenotypes did not appear to have extant progenitors. We infer that the progenitors of these species were derived from extinct diploid lineages or ancestral lineages of multiple extant diploids. In agreement with other recent studies, we recovered molecular evidence of multiple phylogenetically distinct origins for several polyploid taxa, including the widespread octoploid M. albicaulis. ? Conclusions: Evidence of high levels of introgression and allopolyploidy suggests that hybridization has played an important role in the evolution of Trachyphytum. Although idh sequences exhibited complicated evolution, including gene duplication, deletion, and recombination, they provided a higher percentage of informative characters for phylogeny reconstruction than the most variable plastid regions, allowing tests of hypotheses regarding polyploid origins. Given the necessity for rapidly evolving low-copy nuclear genes, researchers studying hybridization and polyploidy may increasingly turn to complex sequence data.     

Molecular phylogenetic relationship of snow finch complex (genera Montifringilla, Pyrgilauda, and Onychostruthus) from the Tibetan plateau

The snow finch complex (Montifringilla, Pyrgilauda, and Onychostruthus) has its center of distribution on the Tibetan plateau, with six out of seven species in the genera occurring there. Phylogenetic relationships among these six species of three genera have been studied based on DNA sequence data obtained from the mitochondrial cytochrome b gene and the nuclear myoglobin gene. The results support monophyly of the snow finch complex group and three major evolutionary lineages are recognized. The first clade consists of ruficollis, blanfordi, and davidiana. These three taxa are sometimes placed in their own genus, Pyrgilauda, and the DNA data supports this. The three taxa nivalis, henrici, and adamsi have traditionally been placed in the genus Montifringilla, and they group together strongly in the present analysis. The results further suggest that nivalis and adamsi are more closely related to each other than are nivalis and henrici, despite that the latter two are often regarded as conspecific. The third distinct lineage within the snow finch complex consists of taczanowskii, which has been placed its own genus, Onychostruthus. This taxon has a basal position in the phylogenetic tree and is sister to all other snow finches. We estimated that taczanowskii split from the other taxa between 2 and 2.5 mya, i.e., about the time for the most recent uplift of the Tibetan plateau, "the Tibet movement", 3.6-1.7 mya. Cladogenesis within the Montifringilla and Pyrgilauda clades seems to be contemporary with the second phase of "Tibet movement" at 2.5 mya and the third phase at 1.7 mya and "Kunhuang movement" in 1.5-0.6 mya. The dramatic climatic and ecological changes following from the uplift of the Tibetan plateau, together with the cyclic contraction and expansion of suitable habitats during the Pleistocene, are probably the most important factors for the cladogenesis in snow finch complex.     

Rapid radiation of Rheum (Polygonaceae) and parallel evolution of morphological traits

In this study, we examined diversification history of Rheum and tested the hypothesis that morphological traits related to plant 'body-plans' evolved in parallel in this genus. We sequenced eight chloroplast DNA fragments (representing more than 8000 bps of the sequence for each species) of 34 species from the genus and 13 species from closely related genera. Phylogenetic analyses indicate that all species of Rheum form a monophyletic lineage sister to the two-species genus Oxyria, indicating that radiative diversifications have occurred in its evolutionary history. Our dating analyses further suggest that these radiations largely coincided with the extensive uplifts of the Qinghai-Tibetan Plateau (QTP). Ancestral state reconstruction and likelihood sensitivity tests strongly indicate that both decumbent and 'glasshouse-like' body-plan traits evolved in parallel in different clades. Our findings highlight the importance of the uplift of the QTP in promoting species diversification and the parallel evolution of morphological traits that are putatively adaptive during such an evolutionary history.     

Multiple origins of polyploidy and the geographic structure of Heuchera grossulariifolia

Multiple origins of polyploidy from an ancestral diploid plant species were investigated using restriction site polymorphism and sequence variation in the chloroplast DNA (cpDNA) of Heuchera grossulariifolia (Saxifragaceae). Phylogenetic analysis indicated that autopolyploidy has arisen at least twice in the evolutionary history of this species and potentially up to as many as seven times. These results suggest a greater range of independent polyploid origins as compared to a previous study of H. grossulariifolia using cpDNA restriction sites that indicated a minimum of three independent origins. Moreover, most polyploid populations did not contain cpDNA haplotypes from a single origin, but rather combined haplotypes from at least two polyploid origins. Past migration among polyploid populations of independent origin or localized polyploid formation may explain the distribution of polyploid haplotypes within and among populations. The analysis also revealed a discrepancy between relatedness and geographical location. In nearly all sympatric populations of diploids and polyploids, polyploids had the same cpDNA haplotypes as diploids from a geographically remote population. This geographical discordance has several possible explanations, including small sample sizes, extinction of parental diploid haplotypes, chloroplast introgression, and homoplasy in the cpDNA sequence data. We conclude that the recurrent formation of polyploids is an important evolutionary mechanism in the diversification of H. grossulariifolia .     

Exploring taxonomic diversity and biogeography of the family Nemacheilinae (Cypriniformes)

Nemacheilidae, in the superfamily Cobitoidea, is comprised of many of morphologically similar fish species that occur in Eurasian water bodies. This large group shows inconsistencies between traditional morphological taxonomy and molecular phylogenetic data. We used mitochondrial genomes, recombinase‐activating gene proteins 1 (RAG1) and the mitochondrial cytochrome c oxidase I gene (COI) to study the phylogenetic relationships among Nemacheilidae species using Bayesian inference and maximum likelihood approaches. Phylogenetic analyses based on mitogenomes provided support for two clades (I and II). The mitogenomes, RAG1, and COI results indicated that several species and genera were not consistent with the traditional morphological subdivisions. The two clades inferred from mitogenomes showed clear geographical patterns. The Tibetan Plateau, Hengduan Mountains, and the Iran Plateau may act as a barrier dividing the clades. The estimated timing of clades separation (36.05 million years ago) coincides with the first uplift of the Tibetan Plateau. We conclude that the geological history of the Tibetan Plateau played a role in the diversification and distribution of the Nemacheilidae taxa. These results provided a phylogenetic framework for future studies of this complex group.     

A complex biogeographic history of diversification in Neotropical lancehead pitvipers (Serpentes,Viperidae)

Based on the literature, we had predicted that the diversification within the Neotropical snake genus Bothrops occurred along a latitudinal gradient from north to south, with diversification into unoccupied niches through ecological opportunity, not correlated with geoclimatic events. Using a dated phylogeny and estimating likelihoods of ancestral states at cladogenesis events, we reconstructed ancestral areas and assessed major events of the diversification of Bothrops clades, and we also discuss systematic implications for this group. Based on the phylogeny we produced, B. lojanus was not considered as part of the genus Bothrops since the results recovered this species nested within the Bothrocophias clade. We infer that the diversification of the Miocene Bothrops pictus and Bothrops alternatus clades may be related to the uplift of the western slopes of the Andes and the Argentinian Patagonian Andes, respectively. The Pliocene Bothrops taeniatus and Bothrops osbornei clades may be related to the uplift of the eastern and northern Andes, respectively. The Plio-Pleistocene Bothrops neuwiedi clade may be related to the habitat transitions from a warmer and forested environment to a cooler and open landscape; the Bothrops jararaca (i.e. island endemic species) and Bothrops lanceolatus clades to over-water dispersal with island speciation; and Bothrops atrox clade to the appearance of the Panamanian land bridge. We found that a multitemporal and multidirectional history of diversification may be correlated with geoclimatic and dispersalist events. We argue that the vacant niche hypothesis by itself does not explain Bothrops diversification.