Estimating ancestral distributions of lineages with uncertain sister groups： a statistical approach to Dispersal-Vicariance Analysis and a case using Aesculus L. （Sapindaceae） including fossils

We propose a simple statistical approach for using Dispersal-Vicariance Analysis （DIVA） software to infer biogeographic histories without fully bifurcating trees. In this approach, ancestral ranges are first optimized for a sample of Bayesian trees. The probability P of an ancestral range r at a node is then calculated as P（rY） = ∑t^n=1 F（rY）t Pt where Y is a node, and F（rY） is the frequency of range r among all the optimal solutions resulting from DIVA optimization at node Y, t is one of n topologies optimized, and Pt is the probability of topology t. Node Y is a hypothesized ancestor shared by a specific crown lineage and the sister of that lineage ＂x＂, where x may vary due to phylogenetic uncertainty （polytomies and nodes with posterior probability 〈 100%）. Using this method, the ancestral distribution at Y can be estimated to provide inference of the geographic origins of the specific crown group of interest. This approach takes into account phylogenetic uncertainty as well as uncertainty from DIVA optimization. It is an extension of the previously described method called Bayes-DIVA, which pairs Bayesian phylogenetic analysis with biogeographic analysis using DIVA. Further, we show that the probability P of an ancestral range at Y calculated using this method does not equate to pp＊F（rY） on the Bayesian consensus tree when both variables are 〈 100%, where pp is the posterior probability and F（rY） is the frequency of range r for the node containing the specific crown group. We tested our DIVA-Bayes approach using Aesculus L., which has major lineages unresolved as a polytomy. We inferred the most probable geographic origins of the five traditional sections of Aesculus and ofAesculus californica Nutt. and examined range subdivisions at parental nodes of these lineages. Additionally, we used the DIVA-Bayes data from Aesculus to quantify the effects on biogeographic inference of including two wildcard fossil taxa in phylogenetic analysis. Our analysis resolved the geographic     

Inferring the biogeographic origins of inter‐continental disjunct endemics using a Bayes‐DIVA approach

The arcto‐Tertiary relictual flora is comprised of many genera that occur non‐contiguously in the temperate zones of eastern Asia, Europe, eastern North America, and western North America. Within each distributional area, species are typically endemic and may thus be widely separated from closely related species within the other areas. It is widely accepted that this common pattern of distribution resulted from of the fragmentation of a once more‐continuous arcto‐Tertiary forest. The historical biogeographic events leading to the present‐day disjunction have often been investigated using a phylogenetic approach. Limitations to these previous studies have included phylogenetic uncertainty and uncertainty in ancestral range reconstructions. However, the recently described Bayes‐DIVA method handles both types of uncertainty. Thus, we used Bayes‐DIVA analysis to reconstruct the stem lineage distributions for 185 endemic lineages from 23 disjunct genera representing 17 vascular plant families. In particular, we asked whether endemic lineages within each of the four distributional areas more often evolved from (1) widespread ancestors, (2) ancestors dispersed from other areas, or (3) endemic ancestors. We also considered which of these three biogeographic mechanisms may best explain the origins of arcto‐Tertiary disjunct endemics in the neotropics. Our results show that eastern Asian endemics more often evolved from endemic ancestors compared to endemics in Europe and eastern and western North America. Present‐day endemic lineages in the latter areas more often arose from widespread ancestors. Our results also provide anecdotal evidence for the importance of dispersal in the biogeographic origins of arcto‐Tertiary species endemic in the neotropics.     

From Africa via Europe to South America: migrational route of a species‐rich genus of Neotropical lowland rain forest trees (Guatteria,Annonaceae)

Aim Several recent studies have suggested that a substantial portion of today’s plant diversity in the Neotropics has resulted from the dispersal of taxa into that region rather than by vicariance. In general, three routes have been documented for the dispersal of taxa onto the South American continent: (1) via the North Atlantic Land Bridge, (2) via the Bering Land Bridge, or (3) from Africa directly onto the continent. Here a species‐rich genus of Neotropical lowland rain forest trees (Guatteria, Annonaceae) is used as a model to investigate these three hypotheses. Location The Neotropics. Methods The phylogenetic relationships within the long‐branch clade of Annonaceae were reconstructed (using maximum parsimony, maximum likelihood and Bayesian inference) in order to gain insight in the phylogenetic position of Guatteria. Furthermore, Bayesian molecular dating and Bayesian dispersal–vicariance (Bayes‐DIVA) analyses were undertaken. Results Most of the relationships within the long‐branch clade of Annonaceae were reconstructed and had high support. However, the relationship between the Duguetia clade, the Xylopia–Artabotrys clade and Guatteria remained unclear. The stem node age estimate of Guatteria ranged between 49.2 and 51.3 Ma, whereas the crown node age estimate ranged between 11.4 and 17.8 Ma. For the ancestral area of Guatteria and its sister group, the area North America–Africa was reconstructed in 99% of 10,000 DIVA analyses, while South America–North America was found just 1% of the time. Main conclusions The estimated stem to crown node ages of Guatteria in combination with the Bayes‐DIVA analyses imply a scenario congruent with an African origin followed by dispersal across the North Atlantic Land Bridge in the early to middle Eocene and further dispersal into North and Central America (and ultimately South America) in the Miocene. The phylogenetically and morphologically isolated position of the genus is probably due to extinction of the North American and European stem lineages in the Tertiary.     

An evaluation of new parsimony‐based versus parametric inference methods in biogeography: a case study using the globally distributed plant family Sapindaceae

Aim Recently developed parametric methods in historical biogeography allow researchers to integrate temporal and palaeogeographical information into the reconstruction of biogeographical scenarios, thus overcoming a known bias of parsimony‐based approaches. Here, we compare a parametric method, dispersal–extinction–cladogenesis (DEC), against a parsimony‐based method, dispersal–vicariance analysis (DIVA), which does not incorporate branch lengths but accounts for phylogenetic uncertainty through a Bayesian empirical approach (Bayes‐DIVA). We analyse the benefits and limitations of each method using the cosmopolitan plant family Sapindaceae as a case study. Location World‐wide. Methods Phylogenetic relationships were estimated by Bayesian inference on a large dataset representing generic diversity within Sapindaceae. Lineage divergence times were estimated by penalized likelihood over a sample of trees from the posterior distribution of the phylogeny to account for dating uncertainty in biogeographical reconstructions. We compared biogeographical scenarios between Bayes‐DIVA and two different DEC models: one with no geological constraints and another that employed a stratified palaeogeographical model in which dispersal rates were scaled according to area connectivity across four time slices, reflecting the changing continental configuration over the last 110 million years. Results Despite differences in the underlying biogeographical model, Bayes‐DIVA and DEC inferred similar biogeographical scenarios. The main differences were: (1) in the timing of dispersal events – which in Bayes‐DIVA sometimes conflicts with palaeogeographical information, and (2) in the lower frequency of terminal dispersal events inferred by DEC. Uncertainty in divergence time estimations influenced both the inference of ancestral ranges and the decisiveness with which an area can be assigned to a node. Main conclusions By considering lineage divergence times, the DEC method gives more accurate reconstructions that are in agreement with palaeogeographical evidence. In contrast, Bayes‐DIVA showed the highest decisiveness in unequivocally reconstructing ancestral ranges, probably reflecting its ability to integrate phylogenetic uncertainty. Care should be taken in defining the palaeogeographical model in DEC because of the possibility of overestimating the frequency of extinction events, or of inferring ancestral ranges that are outside the extant species ranges, owing to dispersal constraints enforced by the model. The wide‐spanning spatial and temporal model proposed here could prove useful for testing large‐scale biogeographical patterns in plants.     

A comparative study in ancestral range reconstruction methods: retracing the uncertain histories of insular lineages

Island systems have long been useful models for understanding lineage diversification in a geographic context, especially pertaining to the importance of dispersal in the origin of new clades. Here we use a well-resolved phylogeny of the flowering plant genus Cyrtandra (Gesneriaceae) from the Pacific Islands to compare four methods of inferring ancestral geographic ranges in islands: two developed for character-state reconstruction that allow only single-island ranges and do not explicitly associate speciation with range evolution (Fitch parsimony [FP; parsimony-based] and stochastic mapping [SM; likelihood-based]) and two methods developed specifically for ancestral range reconstruction, in which widespread ranges (spanning islands) are integral to inferences about speciation scenarios (dispersal-vicariance analysis [DIVA; parsimony-based] and dispersal-extinction-cladogenesis [DEC; likelihood-based]). The methods yield conflicting results, which we interpret in light of their respective assumptions. FP exhibits the least power to unequivocally reconstruct ranges, likely due to a combination of having flat (uninformative) transition costs and not using branch length information. SM reconstructions generally agree with a prior hypothesis about dispersal-driven speciation across the Pacific, despite the conceptual mismatch between its character-based model and this mode of range evolution. In contrast with narrow extant ranges for species of Cyrtandra, DIVA reconstructs broad ancestral ranges at many nodes. DIVA results also conflict with geological information on island ages; we attribute these conflicts to the parsimony criterion not considering branch lengths or time, as well as vicariance being the sole means of divergence for widespread ancestors. DEC analyses incorporated geological information on island ages and allowed prior hypotheses about range size and dispersal rates to be evaluated in a likelihood framework and gave more nuanced inferences about range evolution and the geography of speciation than other methods tested. However, ancestral ranges at several nodes could not be conclusively resolved, due possibly to uncertainty in the phylogeny or the relative complexity of the underlying model. Of the methods tested, SM and DEC both converge on plausible hypotheses for area range histories in Cyrtandra, due in part to the consideration of branch lengths and/or timing of events. We suggest that DEC model-based methods for ancestral range inference could be improved by adopting a Bayesian SM approach, in which stochastic sampling of complete geographic histories could be integrated over alternative phylogenetic topologies. Likelihood-based estimates of ancestral ranges for Cyrtandra suggest a major dispersal route into the Pacific through the islands of Fiji and Samoa, motivating future biogeographic investigation of this poorly known region.     

Accounting for phylogenetic uncertainty in biogeography: a Bayesian approach to dispersal-vicariance analysis of the thrushes (Aves: Turdus)

The phylogeny of the thrushes (Aves: Turdus) has been difficult to reconstruct due to short internal branches and lack of node support for certain parts of the tree. Reconstructing the biogeographic history of this group is further complicated by the fact that current implementations of biogeographic methods, such as dispersal-vicariance analysis (DIVA; Ronquist, 1997), require a fully resolved tree. Here, we apply a Bayesian approach to dispersal-vicariance analysis that accounts for phylogenetic uncertainty and allows a more accurate analysis of the biogeographic history of lineages. Specifically, ancestral area reconstructions can be presented as marginal distributions, thus displaying the underlying topological uncertainty. Moreover, if there are multiple optimal solutions for a single node on a certain tree, integrating over the posterior distribution of trees often reveals a preference for a narrower set of solutions. We find that despite the uncertainty in tree topology, ancestral area reconstructions indicate that the Turdus clade originated in the eastern Palearctic during the Late Miocene. This was followed by an early dispersal to Africa from where a worldwide radiation took place. The uncertainty in tree topology and short branch lengths seems to indicate that this radiation took place within a limited time span during the Late Pliocene. The results support the role of Africa as a probable source area for intercontinental dispersals as suggested for other passerine groups, including basal diversification within the songbird tree.     

Deep-sea squat lobster biogeography (Munidopsidae: Leiogalathea) unveils Tethyan vicariance and evolutionary patterns shared by shallow-water relatives

The ecology, abundance and diversity of galatheoid squat lobsters make them an ideal group to study deep-sea diversification processes. Here, we reconstructed the evolutionary and biogeographic history of Leiogalathea, a genus of circum-tropical deep-sea squat lobsters, in order to compare patterns and processes that have affected shallow-water and deep-sea squat lobster species. We first built a multilocus phylogeny and a calibrated species tree with a relaxed clock using StarBEAST2 to reconstruct evolutionary relationships and divergence times among Leiogalathea species. We used BioGeoBEARS and a DEC model, implemented in RevBayes, to reconstruct ancestral distribution ranges and the biogeographic history of the genus. Our results showed that Leiogalathea is monophyletic and comprises four main lineages; morphological homogeneity is common within and between clades, except in one; the reconstructed ancestral range of the genus is in the Atlantic and Indian oceans (Tethys). They also revealed the divergence of the Atlantic species around 25 million years ago (Ma), intense cladogenesis 15–25 Ma and low levels of speciation over the last 5 million years (Myr). The four Leiogalathea lineages showed similar patterns of speciation: allopatric speciation followed by range expansion and subsequent stasis. Leiogalathea started diversifying during the Oligocene, likely in the Tethyan. The Atlantic lineage then split from its Indo-Pacific sister group due to vicariance driven by closure of the Tethys Seaway. The Atlantic lineage is less speciose compared with the Indo-Pacific lineages, with the Tropical Southwestern Pacific being the current centre of diversity. Leiogalathea diversification coincided with cladogenetic peaks in shallow-water genera, indicating that historical biogeographic events similarly shaped the diversification and distribution of both deep-sea and shallow-water squat lobsters.     

A likelihood framework for inferring the evolution of geographic range on phylogenetic trees

Abstract At a time when historical biogeography appears to be again expanding its scope after a period of focusing primarily on discerning area relationships using cladograms, new inference methods are needed to bring more kinds of data to bear on questions about the geographic history of lineages. Here we describe a likelihood framework for inferring the evolution of geographic range on phylogenies that models lineage dispersal and local extinction in a set of discrete areas as stochastic events in continuous time. Unlike existing methods for estimating ancestral areas, such as dispersal‐vicariance analysis, this approach incorporates information on the timing of both lineage divergences and the availability of connections between areas (dispersal routes). Monte Carlo methods are used to estimate branch‐specific transition probabilities for geographic ranges, enabling the likelihood of the data (observed species distributions) to be evaluated for a given phylogeny and parameterized paleogeographic model. We demonstrate how the method can be used to address two biogeographic questions: What were the ancestral geographic ranges on a phylogenetic tree? How were those ancestral ranges affected by speciation and inherited by the daughter lineages at cladogenesis events? For illustration we use hypothetical examples and an analysis of a Northern Hemisphere plant clade (Cercis), comparing and contrasting inferences to those obtained from dispersal‐vicariance analysis. Although the particular model we implement is somewhat simplistic, the framework itself is flexible and could readily be modified to incorporate additional sources of information and also be extended to address other aspects of historical biogeography.     

ORIGIN AND BIOGEOGRAPHY OF AESCULUS L. (HIPPOCASTANACEAE): A MOLECULAR PHYLOGENETIC PERSPECTIVE

Sequences of chloroplast gene matK and internal transcribed spacers of nuclear ribosomal RNA genes were used for phylogenetic analyses of Aesculus, a genus currently distributed in eastern Asia, eastern and western North America, and southeastern Europe. Phylogenetic relationships inferred from these molecular data are highly correlated with the geographic distributions of species. The identified lineages closely correspond to the five sections previously recognized on the basis of morphology. Ancestral character-state reconstruction, a molecular clock, and fossil evidence were used to infer the origin and biogeographic history of the genus within a phylogenetic framework. Based on the molecular phylogenetic reconstruction of the genus, sequence divergence, and paleontological evidence, we infer that the genus originated during the transition from the Cretaceous to the Tertiary (~65 M.Y.B.P.) at a high latitude in eastern Asia and spread into North America and Europe as an element of the “boreotropical flora”; the current disjunct distribution of the genus resulted from geological and climatic changes during the Tertiary.     

Molecular phylogeography of Troglophilus cave crickets (Orthoptera,Rhaphidophoridae): A combination of vicariance and dispersal drove diversification in the East Mediterranean region

In this study, we investigated the molecular phylogenetic divergence and historical biogeography of cave crickets belonging to the genus Troglophilus (Orthoptera, Rhaphidophoridae) from caves in eastern Mediterranean and Anatolia regions. Three mitochondrial DNA genes (COI, 12S rDNA, and 16S rDNA) and two nuclear ones (18S rDNA and 28S rDNA) were amplified and partially sequenced to reconstruct phylogenetic relationships among most of the known Troglophilus species. Results showed a well‐resolved phylogeny with three main clades representing the Balkan, the Anatolian, and the Cycladian–Cretan lineages. Based on Bayesian analyses, we applied a relaxed molecular clock model to estimate the divergence times between these three lineages. Dating estimates indicate that radiation of the ingroup might have been triggered by the opening of the Mid‐Aegean trench, while the uplift of the Anatolian Plateau in Turkey and the changes of relief, emergence, and disappearance of orographic and hydrographical barriers in the Balkan Peninsula are potential paleogeographic events responsible for the initial diversification of the genus Troglophilus. A possible biogeographic scenario, reconstructed using S‐DIVA with RASP software, suggested that the current distribution of Troglophilus species can be explained by a combination of both dispersal and vicariance events that occurred in particular in the ancestral populations. The radiation of Troglophilus species likely started from the Aegean and proceeded eastward to Anatolia and westward to the Balkan region. Results are additionally compared to those available for Dolichopoda, the only other representative genus of Rhaphidophoridae present in the Mediterranean area.     

Different filtering strategies of genotyping‐by‐sequencing data provide complementary resolutions of species boundaries and relationships in a clade of sexually deceptive orchids

Ongoing hybridization and retained ancestral polymorphism in rapidly radiating lineages could mask recent cladogenetic events. This presents a challenge for the application of molecular phylogenetic methods to resolve differences between closely related taxa. We reanalyzed published genotyping‐by‐sequencing (GBS) data to infer the phylogeny of four species within the Ophrys sphegodes complex, a recently radiated clade of orchids. We used different data filtering approaches to detect different signals contained in the dataset generated by GBS and estimated their effects on maximum likelihood trees, global F_ST and bootstrap support values. We obtained a maximum likelihood tree with high bootstrap support, separating the species by using a large dataset based on loci shared by at least 30% of accessions. Bootstrap and F_ST values progressively decreased when filtering for loci shared by a higher number of accessions. However, when filtering more stringently to retain homozygous and organellar loci, we identified two main clades. These clades group individuals independently from their a priori species assignment, but were associated with two organellar haplotype clusters. We infer that a less stringent filtering preferentially selects for rapidly evolving lineage‐specific loci, which might better delimit lineages. In contrast, when using homozygous/organellar DNA loci the signature of a putative hybridization event in the lineage prevails over the most recent phylogenetic signal. These results show that using differing filtering strategies on GBS data could dissect the organellar and nuclear DNA phylogenetic signal and yield novel insights into relationships between closely related species.     

Reconstructing ancestral ranges in historical biogeography: properties and prospects

Recent years have witnessed a proliferation of quantitative methods for biogeographic inference. In particular, novel parametric approaches represent exciting new opportunities for the study of range evolution. Here, we review a selection of current methods for biogeographic analysis and discuss their respective properties. These methods include generalized parsimony approaches, weighted ancestral area analysis, dispersal-vicariance analysis, the dispersal--extinction--cladogenesis model and other maximum likelihood approaches, and Bayesian stochastic mapping of ancestral ranges, including a novel approach to inferring range evolution in the context of island biogeography. Some of these methods were developed specifically for problems of ancestral range reconstruction, whereas others were designed for more general problems of character state reconstruction and subsequently applied to the study of ancestral ranges. Methods for reconstructing ancestral history on a phylogenetic tree differ not only in the types of ancestral range states that are allowed, but also in the various historical events that may change the ancestral ranges. We explore how the form of allowed ancestral ranges and allowed transitions can both affect the outcome of ancestral range estimation. Finally, we mention some promising avenues for future work in the development of model-based approaches to biogeographic analysis.     

Dynamic diversification history with rate upshifts in Holarctic bell‐flowers (Campanula and allies)

Campanula s.l. is one of the most speciose flowering plant lineages of the Holarctic (ca. 600 species). In the present study we sequenced three regions of the plastid genome (petD, rpl16 and trnK/matK) across a broad sample of Campanula s.l., which markedly improved phylogenetic resolution and statistical support compared to previous studies. Based on this robust phylogenetic hypothesis we estimated divergence times using BEAST, diversification rate shifts using Bayesian Analysis of Macroevolutionary Mixture (BAMM) and TreePar, and ancestral ranges using Biogeography with Bayesian (and likelihood) Evolutionary Analyses in R. Campanula s.l. is estimated to have originated during the Early Eocene but the major diversification events occurred between the Late Oligocene and Middle Miocene. Two upward diversification rate shifts were revealed by BAMM, specific to the crown nodes of two Campanula clades: CAM17, a mostly South European‐SW Asian lineage originating during the Middle Miocene and containing nearly half of all known Campanula species; and CAM15B, a SW Asian–Sino‐Himalayan lineage of nine species originating in the early Pleistocene. The dynamic diversification history of Campanula and the inferred rate shifts are discussed in a geo‐historical context.     

The biogeography of tropical reef fishes: endemism and provinciality through time

The largest marine biodiversity hotspot straddles the Indian and Pacific Oceans, driven by taxa associated with tropical coral reefs. Centred on the Indo‐Australian Archipelago (IAA), this biodiversity hotspot forms the ‘bullseye’ of a steep gradient in species richness from this centre to the periphery of the vast Indo‐Pacific region. Complex patterns of endemism, wide‐ranging species and assemblage differences have obscured our understanding of the genesis of this biodiversity pattern and its maintenance across two‐thirds of the world's oceans. But time‐calibrated molecular phylogenies coupled with ancestral biogeographic estimates have provided a valuable framework in which to examine the origins of coral reef fish biodiversity across the tropics. Herein, we examine phylogenetic and biogeographic data for coral reef fishes to highlight temporal patterns of marine endemism and tropical provinciality. The ages and distribution of endemic lineages have often been used to identify areas of species creation and demise in the marine tropics and discriminate among multiple hypotheses regarding the origins of biodiversity in the IAA. Despite a general under‐sampling of endemic fishes in phylogenetic studies, the majority of locations today contain a mixture of potential paleo‐ and neo‐endemic fishes, pointing to multiple historical processes involved in the origin and maintenance of the IAA biodiversity hotspot. Increased precision and sampling of geographic ranges for reef fishes has permitted the division of discrete realms, regions and provinces across the tropics. Yet, such metrics are only beginning to integrate phylogenetic relatedness and ancestral biogeography. Here, we integrate phylogenetic diversity with ancestral biogeographic estimation of lineages to show how assemblage structure and tropical provinciality has changed through time.     

Fossil biogeography: a new model to infer dispersal,extinction and sampling from palaeontological data

Methods in historical biogeography have revolutionized our ability to infer the evolution of ancestral geographical ranges from phylogenies of extant taxa, the rates of dispersals, and biotic connectivity among areas. However, extant taxa are likely to provide limited and potentially biased information about past biogeographic processes, due to extinction, asymmetrical dispersals and variable connectivity among areas. Fossil data hold considerable information about past distribution of lineages, but suffer from largely incomplete sampling. Here we present a new dispersal–extinction–sampling (DES) model, which estimates biogeographic parameters using fossil occurrences instead of phylogenetic trees. The model estimates dispersal and extinction rates while explicitly accounting for the incompleteness of the fossil record. Rates can vary between areas and through time, thus providing the opportunity to assess complex scenarios of biogeographic evolution. We implement the DES model in a Bayesian framework and demonstrate through simulations that it can accurately infer all the relevant parameters. We demonstrate the use of our model by analysing the Cenozoic fossil record of land plants and inferring dispersal and extinction rates across Eurasia and North America. Our results show that biogeographic range evolution is not a time-homogeneous process, as assumed in most phylogenetic analyses, but varies through time and between areas. In our empirical assessment, this is shown by the striking predominance of plant dispersals from Eurasia into North America during the Eocene climatic cooling, followed by a shift in the opposite direction, and finally, a balance in biotic interchange since the middle Miocene. We conclude by discussing the potential of fossil-based analyses to test biogeographic hypotheses and improve phylogenetic methods in historical biogeography.     

Diversification and asymmetrical gene flow across time and space: lineage sorting and hybridization in polytypic barking frogs

Young species complexes that are widespread across ecologically disparate regions offer important insights into the process of speciation because of their relevance to how local adaptation and gene flow influence diversification. We used mitochondrial DNA and up to 28 152 genomewide single nucleotide polymorphisms from polytypic barking frogs (Craugastor augusti complex) to infer phylogenetic relationships and test for the signature of introgressive hybridization among diverging lineages. Our phylogenetic reconstructions suggest (i) a rapid Pliocene–Pleistocene radiation that produced at least nine distinct lineages and (ii) that geographic features of the arid Central Mexican Plateau contributed to two independent northward expansions. Despite clear lineage differentiation (many private alleles and high between‐lineage F_ST scores), D‐statistic tests, which differentiate introgression from ancestral polymorphism, allowed us to identify two putative instances of reticulate gene flow. Partitioned D‐statistics provided evidence that these events occurred in the same direction between clades but at different points in time. After correcting for geographic distance, we found that lineages involved in hybrid gene flow interactions had higher levels of genetic variation than independently evolving lineages. These findings suggest that the nature of hybrid compatibility can be conserved overlong periods of evolutionary time and that hybridization between diverging lineages may contribute to standing levels of genetic variation.     

Systematic revision,molecular phylogeny and biogeography of the antlion tribe Acanthoplectrini (Neuroptera: Myrmeleontidae: Dendroleontinae), with emphasis on the Oriental lineage

The tribe Acanthoplectrini (Myrmeleontidae: Dendroleontinae) includes a group of antlion genera widely distributed across the Australasian and Oriental regions. The intergeneric and interspecific relationships between or within the Australian and Oriental lineages of this tribe as well as their historical biogeography remain largely unexplored. Here, we present a molecular phylogenetic and biogeographic analyses of Acanthoplectrini to infer the diversification history of this tribe, with emphasis on the Oriental lineage. Both the Oriental and Australian lineages are monophyletic and recovered as sister groups. Ancestral area reconstruction suggests that the ancestor of Acanthoplectrini might have been once widely distributed from Indochina to Australia and then split into the Oriental and Australian lineages during the early-Miocene. Our analyses recovered northeastern Indochina and south China as the ancestral range of the Oriental Acanthoplectrini. During the mid-Miocene to the mid-Pliocene, orographic events such as the rising of mountain ranges (including the Himalayas) and the formation of major islands in southeastern Asia triggered several dispersal and vicariance events in the Oriental Acanthoplectrini, driving their speciation. We revise the classification of the Oriental Acanthoplectrini, establishing the new genus Paralayahima gen. n. , which is recovered sister to Layahima Navás. Moreover, we describe four new species of Layahima, Layahima aspoeckorum sp. n. , Layahima monba sp. n. , Layahima lhoba sp. n. and Layahima xinliae sp. n. , and we reinstate two previously synonymized species, Layahima melanocoris (Yang) stat. rev. and comb. n. and Layahima nebulosa Navás stat. rev.     

Museum archives revisited: Central Asiatic hawkmoths reveal exceptionally high late Pliocene species diversification (Lepidoptera,Sphingidae)

Three high‐elevation Hyles species of Central Asia have proven difficult to sample, and thus, only a limited number of specimens are available for study. Ancient DNA techniques were applied to sequence two mitochondrial genes from ‘historic’ museum specimens of Hyles gallii, Hyles renneri and Hyles salangensis to elucidate the phylogenetic relationships of these species. This approach enabled us to include the holotypes and paratypes. The status of H. salangensis as a species endemic to a mountain range north of Kabul in Afghanistan is confirmed by this study. It is most closely related to Hyles nicaea and H. gallii, and quite distant from the clade comprising the species from Hyles vespertilio to Hyles tithymali, despite this group and H. salangensis both completely lacking an arolium on the pretarsus. Our results show that the samples assigned to H. renneri and Hyles livornica tatsienluica are conspecific and so we reinstate Hyles tatsienluica stat. nov. as the valid name for this species and synonymize H. renneri syn. nov. with it. This study shows that the distribution range postulated for H. tatsienluica extends from Nepal well into the mountains of south‐western China. The distribution ranges of H. livornica and H. tatsienluica overlap. The study confirms the previously proposed synonymies of Hyles nepalensis, Hyles gallii intermedia and Hyles gallii tibetanica with H. gallii. Extensive species sampling (over 80% of Hyles species) in this study allowed additional analyses. The dated phylogeny reveals the global Hyles hawkmoth radiation to be much more recent than previously thought: it began in the Late Miocene and culminated in a Pleistocene burst of diversification in the Northern Hemisphere. Ancestral ranges of basal nodes were reconstructed as highly equivocal, but the Neotropics has the highest probability in the two oldest nodes. Although the origins of the Madagascan and Australian species also remain ambiguous, a large crown clade of fifteen species was reconstructed to have originated in the Palaearctic. The wide distribution ranges of the two migratory species, H. livornica and H. gallii, appear to blur any traces of the biogeographic origin of the clades containing these species. Specialization in larval host plant use onto particular plant families from the ancestral condition ‘polyphagous’ may have led to an increased rate of speciation and phylogenetic diversification in three subgroups of Hyles (the Hawaiian clade, the Hyles centralasiae group and the Hyles euphorbiae complex).     

Macaronesia: a source of hidden genetic diversity for post‐glacial recolonization of western Europe in the leafy liverwort Radula lindenbergiana

Aim Bryophytes exhibit apparently low rates of endemism in Macaronesia and differ from angiosperms in their diversity patterns by the widespread occurrence of endemics within and among archipelagos. This paper investigates the phylogeography of the leafy liverwort Radula lindenbergiana to determine: (1) whether or not morphologically cryptic diversification has occurred in Macaronesia, and (2) the relationships between Macaronesian and continental populations. Location Macaronesia, Europe, Africa. Methods Eighty‐four samples were collected across the species’ distribution range and sequenced at four chloroplast DNA (cpDNA) loci (atpB–rbcL, trnG, trnL and rps4). Phylogenetic reconstructions and Bayesian ancestral area reconstructions were used in combination with population genetics statistics (H, N_ST, F_ST) to describe the pattern of present genetic diversity in R. lindenbergiana and infer its biogeographic history. Results Patterns of genetic diversity in R. lindenbergiana exhibit a striking westwards gradient, wherein haplotype (0.90) and nucleotide (0.0038 ± 0.0019) diversity peak in Macaronesia, with a substantial endemic component. We found 20.9% of the genetic variance between biogeographic regions, and most pairwise F_ST comparisons between regions are significantly different from zero. The global N_ST (0.78) is significantly higher than the global F_ST (0.20), providing evidence for the presence of phylogeographic signal in the data. Ancestral area reconstructions suggest that the haplotypes currently found in western Europe share a Macaronesian common ancestor. Main conclusions The haplotype diversification exhibited by R. lindenbergiana in Macaronesia is comparable to that reported for many angiosperm groups at the species level. The apparent lack of radiation among Macaronesian bryophytes may thus reflect the reduced morphology of bryophytes in comparison with angiosperms. The high diversity found among Macaronesian haplotypes, especially in Madeira and the Canary Islands, and the significant N_ST/F_ST ratio between Macaronesia and all the other biogeographic regions (an indication that mutation rate exceeds dispersal rates) suggest that Macaronesian archipelagos could have served as a refugium during the Quaternary glaciations. Many haplotypes currently found in Europe share a Macaronesian common ancestor, and this further suggests that Macaronesia might have played a key role in the back‐colonization of the continent.     

Diversification of Morpho butterflies (Lepidoptera,Nymphalidae): a re‐evaluation of morphological characters and new insight from DNA sequence data

This study compiles previously published morphological, colour and behavioural characters and includes new DNA sequence data for eight markers (one mitochondrial and seven nuclear) to re‐evaluate phylogenetic relationships and estimate times of divergence for Morpho butterflies using parsimony and Bayesian methods. We note an effect of missing data on phylogenetic inference and calculations of Partitioned Bremer Support. Morphology and DNA trees were moderately congruent, and the combined analyses of all data included elements of both sources. Both morphology and DNA support the monophyly of Morpho and the early separation of the sister pair M. marcus plus M. eugenia, but trees from different data sources are congruent mostly at derived nodes, and differ at several internal nodes. The analyses of combined data indicate that Morpho is composed of four clades each of which include one or more previously proposed subgenera. The subgenera Pessonia and Morpho were not monophyletic, and to address this issue we propose that Pessonia, syn.nov. be subsumed within Morpho. The ancestor of Morpho probably arose during the Oligocene, and most diversification seems to have occurred during the late Miocene. S‐DIVA analysis suggests eastern Andean region as the ancestral area for Morpho, and that the South American Atlantic Forest was colonized multiple times.