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.     

A young clade repeating an old pattern: diversity in Nothonotus darters (Teleostei: Percidae) endemic to the Cumberland River

Hypotheses of diversification in eastern North American freshwater fishes have focused primarily on allopatric distributions of species between disjunct highland areas and major river systems. However, these hypotheses do not fully explain the rich diversity of species within highland regions and river systems. Relatively old diversification events at small geographic scales have been observed in the Barcheek Darter subclade that occurs in the Cumberland River drainage (CRD) in Kentucky and Tennessee, United States of America, but it is unknown if this pattern is consistent in other darter subclades. We explored phylogeographic diversity in two species of Nothonotus darters, N. microlepidus and N. sanguifluus, endemic to the CRD to compare phylogenetic patterns between Barcheek Darters and species of Nothonotus. We collected sequence data for a mitochondrial gene (cytb) and three nuclear genes (MLL, S7 and RAG1) from 19 N. microlepidus and 35 N. sanguifluus specimens. Gene trees were estimated using maximum likelihood and Bayesian methods, and a 'species tree' was inferred using a Bayesian method. These trees indicate that species diversity in Nothonotus is underestimated. Five distinct lineages were evident, despite retained ancestral polymorphism and unsampled extirpated populations. Comparison of chronograms for Barcheek Darters and Nothonotus revealed that microendemism resulting from species diversification at small geographic scales in the CRD is a consistent pattern in both old and young darter subclades. Our analyses reveal that geographic isolating mechanisms that result in similar phylogeographic patterns in the CRD are persistent through long expanses of evolutionary time.     

Genetic structure,species limits and evolution of the parasitoid wasp genus Stenocorse (Braconidae: Doryctinae) based on nuclear 3RAD and mitochondrial data

Stenocorse (Braconidae: Doryctinae) is an ectoparasitoid wasp genus widely distributed from southern U.S.A. to central Argentina, which is known to attack over 60 bruchine beetle species that feed on a variety of legumes. This genus currently comprises a single described species, S. bruchivora (Crawford, 1910), though it is probably composed of additional, morphologically similar species. In this study we investigated genetic structure, species limits and evolutionary relationships of sampled populations of Stenocorse obtained along its known geographic distribution, with an emphasis on the Mexican territory. For this, we generated DNA sequences from both the mitochondrial (mt) barcoding locus and nuclear genome data obtained with the 3RAD technique. The Bayesian phylogenetic reconstructions based on the two datasets were mostly congruent, both showing clear geographic structure. Four clusters were recovered with nuclear data, which were also congruent with the geographic distribution of the samples and the species delimitation analyses. We delimited six evolutionary lineages with the two data sources, and three additional species were also delimited with the COI locus alone. Stenocorse probably originated during the mid- to early Miocene, and the diversification events that led to its current species distribution probably occurred between the early Pliocene and late Pleistocene. Our species delineation analyses, phylogenetic reconstructions and rearing records suggest that the species of Stenocorse have generalist feeding habits, contrary to what was expected due to host plant specialization of the bruchine species and the cascading effect that this could have in the species diversification of their parasitoids.     

Molecular phylogeography of Dolichopoda cave crickets (Orthoptera, Rhaphidophoridae): a scenario suggested by mitochondrial DNA

This study focuses on the phylogenetic relationships among a number of West-Mediterranean cave crickets species belonging to Dolichopoda; primarily a Mediterranean genus, distributed from eastern Pyrenees to Caucasus. In this paper, 11 Dolichopoda species from the French Pyrenees (D. linderi), the island of Corsica (D. bormansi and D. cyrnensis), and northern, central, and southern Italy (D. ligustica, D. schiavazzii, D. aegilion, D. baccettii, D. laetitiae, D. geniculata, D. capreensis, and D. palpata) were studied. Two more species, one from the Caucasus, D. euxina, and one from Greece, D. remyi, were also included in the analyses, together with more distant species within the same family to be used as outgroups. Fifteen hundred base pairs of mitochondrial DNA, corresponding to the small subunit of the ribosomal RNA (16S rRNA) and to the subunit I of the cytochrome oxidase I (COI), were sequenced in order to clarify the phylogenetic relationships and biogeography of this group of Mediterranean cave crickets. The molecular data are congruent with a phylogeographic pattern; with the geographically close species also the most related ones. Based on mtDNA divergence, the present-day distribution of genetic diversity seems to have been impacted by climatic events due to glacial and interglacial cycles that have characterized the Pleistocene era.     

Diversification of subgenus Calathus (Coleoptera: Carabidae) in the Mediterranean region – glacial refugia and taxon pulses

Aim To investigate the effects of Pleistocene climatic variations on the diversification rate of the subgenus Calathus (Coleoptera: Carabidae), and to estimate the role of vicariance and dispersal for explaining current distributional patterns. Location Western Palaearctic Region, particularly the Mediterranean Basin. Methods Fragments of the mitochondrial cox1–cox2 and the nuclear 28S and EF1α genes were analysed by Bayesian inference. Lineage divergence times were estimated using a Bayesian relaxed molecular clock. Three diversification rate analyses were conducted, namely gamma (γ)‐statistic, birth–death likelihood (BDL) test and survival analyses, in order to test departures from a constant rate model of diversification. A Bayesian approach to dispersal–vicariance analysis was developed to reconstruct the most probable ancestral area of subgenus Calathus and subsequent events of dispersal and colonization. Results A constant rate of speciation events from the late Miocene onwards was found for the subgenus Calathus, whereas recent Pleistocene climatic oscillations played an important role only in shaping intraspecific diversity. Overall diversification patterns for the subgenus are best explained by at least four westward dispersal events from the eastern Mediterranean Basin. Three distinct phylogroups were found for the widely distributed Calathus fuscipes. Incongruence between mitochondrial and nuclear loci was found for a number of species. Main conclusions Diversification analyses suggest either a constant rate of diversification (BDL analysis) or a decrease in diversification rates for the subgenus (survival or γ‐statistics analyses), but not an increase related to the effects of glaciation cycles. Diversification patterns in the subgenus Calathus agree with predictions of the taxon pulse model. From the middle Miocene onwards the Anatolian Peninsula was possibly the main centre of diversification, with successive dispersal events towards the western Mediterranean Basin. Range expansion and secondary contact zones are postulated between members of different phylogroups in C. fuscipes.     

Phylogeny, biogeography and evolution of clutch size in South American lizards of the genus Kentropyx (Squamata: Teiidae)

The lizard genus Kentropyx (Squamata: Teiidae) comprises nine species, which have been placed in three species groups (calcarata group, associated to forests ecosystems; paulensis and striata groups, associated to open ecosystems). We reconstructed phylogenetic relationships of Kentropyx based on morphology (pholidosis and coloration) and mitochondrial DNA data (12S and 16S), using maximum parsimony and Bayesian methods, and evaluated biogeographic scenarios based on ancestral areas analyses and molecular dating by Bayesian methods. Additionally, we tested the life‐history hypothesis that species of Kentropyx inhabiting open ecosystems (under seasonal environments) produce larger clutches with smaller eggs and that species inhabiting forest ecosystems (under aseasonal conditions) produce clutches with fewer and larger eggs, using Stearns’ phylogenetic‐subtraction method and canonical phylogenetic ordination to take in to account the effects of phylogeny. Our results showed that Kentropyx comprises three monophyletic groups, with K. striata occupying a basal position in opposition to previous suggestions of relationships. Additionally, Bayesian analysis of divergence time showed that Kentropyx may have originated at the Tertiary (Eocene/Oligocene) and the ‘Pleistocene Refuge Hypothesis’ may not explain the species diversification. Based on ancestral reconstruction and molecular dating, we argued that a savanna ancestor is more likely and that historical events during the Tertiary of South America promoted the differentiation of the genus, coupled with recent Quaternary events that were important as dispersion routes and for the diversification at populational levels. Clutch size and egg volume were not significantly different between major clades and ecosystems of occurrence, even accounting for the phylogenetic effects. Finally, we argue that phylogenetic constraints and phylogenetic inertia might be playing essential roles in life history evolution of Kentropyx.     

Historical biogeography of the arid‐adapted velvet ant Sphaeropthalma arota (Hymenoptera: Mutillidae) reveals cryptic species

Aim To investigate the phylogeographic patterns among populations of the wide‐ranging velvet ant Sphaeropthalma arota to (1) examine the biogeographic patterns within this species, and (2) associate major genetic divergences to historical vicariance events in order to gain insight into the processes that drove diversification in arid‐adapted organisms. Location Western North America. Methods Phylogenetic relationships were determined by analysing the two ribosomal DNA internal transcribed spacers (ITS1 and ITS2) in Mr Bayes . Divergence dates were estimated for major nodes using two different molecular dating analyses: a penalized likelihood approach to rate smoothing using the program r8s , and a Bayesian MCMC averaging approach to rate smoothing using the program beast . Both analyses were calibrated using fossils from Dominican amber. Haplotype networks were estimated using tcs . Ecological niche models (ENMs) were developed using six bioclimatic variables from the WorldClim data set in the program Maxent . Results The phylogenetic reconstructions indicate that S. arota can be split into four deeply divergent lineages that probably represent distinct species. No morphological characters were found that can be used to identify these four species, making the S. arota species complex the first documented cryptic species complex in the family Mutillidae. Ecological niche models provide estimated distributions, which indicate that each of the four species inhabits a distinct niche. Divergence date estimates suggest that major diversification events occurred in the late Neogene. Main conclusions Sphaeropthalma arota is composed of four genetically distinct species that cannot be distinguished morphologically based on current methods. We suggest that the members of this group be identified as the S. arota species complex. Major diversification events in this species complex can be linked to late Neogene mountain building and aridification events, specifically the uplift of the mountain ranges in southern California and the expansion of the Bouse Sea.     

Diversification of montane species via elevation shifts: the case of the Kaçkar cricket Phonochorion (Orthoptera)

Localized up‐down altitudinal shifts and subsequent isolation–admixture of montane species in response to glacial cycles has been proposed as a mechanism for the high diversity along Anatolian mountains. However, specific predictions of the proposed mechanism (the elevation shift model) have yet to be tested. Here, we provide a first assessment of this model for promoting inter‐ and intraspecific genetic diversity in the bush‐cricket genus Phonochorion endemic to the West Lesser Caucasus hotspot. Mitochondrial genes were analysed by Bayesian Markov Chain Monte Carlo inferences and coalescent simulations. Timing of diversification was estimated using a multispecies coalescent model. Divergence with gene flow was tested using an isolation with migration model. Population genetic parameters and genetic structuring were determined using Bayesian coalescent methods and spatial analysis. Demographic history was assessed using mismatch distributions and extended Bayesian skyline plots. Speciation events corresponded both to the Miocene and Pleistocene while intraspecific divergence was Pleistocene based. There was evidence for moderate levels of gene flow between species during diversification; however, incomplete lineage sorting could explain the data as well as gene flow. Overall diversification patterns within the genus Phonochorion agree with the predictions of the elevations shift model. Genetic patterns of diversification were driven mainly by Pleistocene glacial cycles and reflected the nature and distribution of sky islands. There was also some albeit weak evidence of demographic expansions coinciding with glacial cooling. However, evidence for divergence with gene flow was inconclusive.     

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.     

Out of South America: phylogeny of non‐biting midges in the genus Labrundinia suggests multiple dispersal events to Central and North America

Non‐biting midges of the genus Labrundinia (Chironomidae: Tanypodinae) are minute dipterans with immature stages living in a variety of unpolluted water bodies, from small streams and ponds to lakes and bays. Extensively recorded in ecological studies, the genus comprises 39 species, all except one described from areas outside the Palearctic region. Internal structure among Labrundinia species was postulated by S. S. Roback, who recognized four species groups based on morphological characters of immature stages. We examined phylogenetic relationships among known Labrundinia species using partial DNA sequences of the nuclear protein‐coding gene CAD and morphological characters. Both analyses with Bayesian inference and parsimony methods recovered the monophyly of Labrundinia, strongly supported by five morphological synapomorphies. Internal relationships within the genus partly supported Roback's species groups with the addition of later described species. Biogeographical inferences were obtained by applying Bayesian binary MCMC (BBM) analysis and favoured a scenario where Labrundinia had its initial diversification in the Neotropical region and that current presence in the Nearctic region and southern South America is due to subsequent dispersal.     

Phylogeny and biogeography of the sandalwoods (Santalum, Santalaceae): repeated dispersals throughout the Pacific

Results of the first genus-wide phylogenetic analysis for Santalum (Santalaceae), using a combination of 18S-26S nuclear ribosomal (ITS, ETS) and chloroplast (3' trnK intron) DNA sequences, provide new perspectives on relationships and biogeographic patterns among the widespread and economically important sandalwoods. Congruent trees based on maximum parsimony, maximum likelihood, and Bayesian methods support an origin of Santalum in Australia and at least five putatively bird-mediated, long-distance dispersal events out of Australia, with two colonizations of Melanesia, two of the Hawaiian Islands, and one of the Juan Fernandez Islands. The phylogenetic data also provide the best available evidence for plant dispersal out of the Hawaiian Islands to the Bonin Islands and eastern Polynesia. Inability to reject rate constancy of Santalum ITS evolution and use of fossil-based calibrations yielded estimates for timing of speciation and colonization events in the Pacific, with dates of 1.0-1.5 million yr ago (Ma) and 0.4-0.6 Ma for onset of diversification of the two Hawaiian lineages. The results indicate that the previously recognized sections Polynesica, Santalum, and Solenantha, the widespread Australian species S. lanceolatum, and the Hawaiian species S. freycinetianum are not monophyletic and need taxonomic revision, which is currently being pursued.     

Diversification of sympatric Sapromyza (Diptera: Lauxaniidae) from Madeira: six morphological species but only four mtDNA lineages

A series of recent studies on speciation of insects within the Canary Islands have indicated considerable within-island diversification, similar to that described in the Hawaiian islands. Little work has yet been carried out on the neighboring Madeiran archipelago, which is also volcanic. This study examines relationships among all known Lauxaniid flies of the genus Sapromyza from Madeira (including six newly described morphological species) based on mitochondrial gene trees constructed from cytochrome c oxidase (subunit I) and 16S rRNA partial sequences. Phylogenies based on maximum likelihood distances, a Bayesian method based on Markov chain Monte Carlo sampling from the posterior probability distribution, and maximum parsimony show that eight of the nine Madeiran species comprise a single monophyletic group. This clade is also split into two subclades representing black- and yellow/orange-bodied forms. The latter mtDNA clade corresponds to only two species (Sapromyza imitans and Sapromyza indigena) which are not reciprocally monophyletic. Monophyly is strongly supported within four of the six black-bodied species but not for the species pair (Sapromyza inconspicua, Sapromyza laurisilvae). We discuss the double occurrence (at least) of introgressive hybridization/incomplete lineage sorting within this group and suggest that recent speciation is the most likely explanation. The remaining species on the island, Sapromyza madeirensis, is very divergent from the aforementioned group, occupying a more basal position in the tree than the other Atlantic island and continental Sapromyza that were included in the analysis. At least two speciation events for Madeiran Sapromyza appear to correspond to quite ancient periods relative to the age of the island, while others are more recent. This suggests that a combination of island colonization and within-island sympatric and/or vicariance-mediated speciation may explain the observed diversity.     

The highly conserved spermatophyte cell wall DUF642 protein family: phylogeny and first evidence of interaction with cell wall polysaccharides in vitro

The evolution of spermatophyte plants involved fundamental changes in cell wall structure and function which resulted from diversification of carbohydrates and proteins. Cell wall proteomic analyses identified a novel family of proteins of yet unknown function, the DUF642 (Domain of Unknown Function 642) proteins. To investigate the evolution of the DUF642 gene family, 154 gene sequences from 24 plant species were analyzed, and phylogenetic inferences were conducted using the Maximum Likelihood and Bayesian Inference methods. Orthologous genes were detected in spermatophyte species and absent in non-seed known plant genomes. Protein sequences shared conserved motifs that defined the signature of the family. Distribution of conserved motifs indicated an ancestral intragenic duplication event. Gene phylogeny documented paleoduplication events originating three or four clades, depending on root position. When based on mid-point rooting, it retrieved four monophyletic clades: A, B, C, and D. A glycosylphosphatidylinositol (GPI)-anchor site and one or two galactose-binding domains-like (GBDLs) could be predicted for some DUF642 proteins. The B, C, and D clades grouped the predicted GPI-anchored proteins. First evidence of in vitro interaction of a DUF642 protein with a cell wall polysaccharide fraction is provided. A competition assay with cellulose prevented this interaction. The degree of diversification and the conservation of the family suggested that DUF642 proteins are key components in seed plant evolution.     

The role of geography,ecology, and hybridization in the evolutionary history of Canary Island Descurainia

Premise

Oceanic islands offer the opportunity to understand evolutionary processes underlying rapid diversification. Along with geographic isolation and ecological shifts, a growing body of genomic evidence has suggested that hybridization can play an important role in island evolution. Here we use genotyping-by-sequencing (GBS) to understand the roles of hybridization, ecology, and geographic isolation in the radiation of Canary Island Descurainia (Brassicaceae).

Methods

We carried out GBS for multiple individuals of all Canary Island species and two outgroups. Phylogenetic analyses of the GBS data were performed using both supermatrix and gene tree approaches and hybridization events were examined using D-statistics and Approximate Bayesian Computation. Climatic data were analyzed to examine the relationship between ecology and diversification.

Results

Analysis of the supermatrix data set resulted in a fully resolved phylogeny. Species networks suggest a hybridization event has occurred for D. gilva, with these results being supported by Approximate Bayesian Computation analysis. Strong phylogenetic signals for temperature and precipitation indicate one major ecological shift within Canary Island Descurainia.

Conclusions

Inter-island dispersal played a significant role in the diversification of Descurainia, with evidence of only one major shift in climate preferences. Despite weak reproductive barriers and the occurrence of hybrids, hybridization appears to have played only a limited role in the diversification of the group with a single instance detected. The results highlight the need to use phylogenetic network approaches that can simultaneously accommodate incomplete lineage sorting and gene flow when studying groups prone to hybridization; patterns that might otherwise be obscured in species trees.     

Plate tectonics and biogeographical patterns of the Pseudophoxinus (Pisces: Cypriniformes) species complex of central Anatolia, Turkey

We investigated the phylogenetic relationships of Pseudophoxinus (Cyprinidae: Leuciscinae) species from central Anatolia, Turkey to test the hypothesis of geographic speciation driven by early Pliocene orogenic events. We analyzed 1141 aligned base pairs of the complete cytochrome b mitochondrial gene. Phylogenetic relationships reconstructed by maximum likelihood, Bayesian likelihood, and maximum parsimony methods are identical, and generally well supported. Species and clades are restricted to geologically well-defined units, and are deeply divergent from each other. The basal diversification of central Anatolian Pseudophoxinus is estimated to have occurred approximately 15 million years ago. Our results are in agreement with a previous study of the Anatolian fish genus Aphanius that also shows a diversification pattern driven by the Pliocene orogenic events. The distribution of clades of Aphanius and Pseudophoxinus overlap, and areas of distribution comprise the same geological units. The geological history of Anatolia is likely to have had a major impact on the diversification history of many taxa occupying central Anatolia; many of these taxa are likely to be still unrecognized as distinct.     

Two pulses of diversification across the Isthmus of Tehuantepec in a montane Mexican bird fauna

Understanding the evolutionary history of the species in a particular region provides insights into how that fauna was formed. Of particular interest to biogeographers is examining the impact a geographical barrier had in generating temporal genetic diversity among codistributed species. We examined the impact a major New World barrier, the Isthmus of Tehuantepec (IT) in southern Mexico, had on a regional bird fauna. Specifically, genetic data from 10 montane-forest bird taxa were analysed using approximate Bayesian computation (ABC) to test the hypothesis of simultaneous intraspecific diversification at the IT. Because effective population size (N_e) has the greatest impact on coalescent times, thereby affecting tests of divergence among codistributed taxa, we chose priors for both current and ancestral N_e using empirical estimates of theta. The ABC method detected two discrete diversification events. Subsequent analysis with the number of diversification events constrained to two suggests that four taxa diverged in an older event, with the remaining six diverging more recently. Application of a range of mutation rates from 2.0 to 5.0% Myr⁻¹ places both events within the Pleistocene or Late Pliocene, suggesting that fluctuations in montane habitat induced by climate cycles and a late Pliocene seaway may have fractured this montane bird fauna. The results presented here suggest this avian fauna responded in a relatively concerted fashion over the last several million years.     

Species tree phylogeny,character evolution,and biogeography of the Patagonian genus <Emphasis Type="Italic">Anarthrophyllum</Emphasis> Benth. (Fabaceae)

Geologic events promoting the aridization of southern South America contributed to lineage divergences and species differentiation through geographic (allopatric divergence) and biotic and abiotic factors (ecological divergence). For the genus Anarthrophyllum, which is distributed in arid and semi-arid regions of Patagonia, we assessed how these factors affected species diversification and reconstructed its possible biogeographic history in South American arid environments. Sequences were obtained from two molecular markers: the ITS nuclear region and the trnS-trnG plastid region. Using Parsimony, Maximum likelihood and Bayesian inference individual gene trees were reconstructed, and a species tree was obtained using multi-species coalescent analysis. Divergence times among species were estimated using secondary calibrations. Flexible Bayesian models and stochastic character mapping were used to elucidate ancestral geographic distributions and the evolution of the floral and vegetative phenotypes in the genus. Gene trees and species tree analyses strongly support Anarthrophyllum as monophyletic; all analyses consistently retrieved three well-supported main clades: High Andean Clade, Patagonian Clade 1, and Patagonian Clade 2. Main diversification events occurred concomitant with the Andean uplift and steppe aridization; the Andean mountain range possibly acted as a species barrier for the High Andean Clade. Vegetative traits showed adaptations to harsh climates in some clades, while pollinator-related floral features were associated with independent diversification in bee- and bird-pollinated clades within both Patagonian Clades. In conclusion, evolutionary and biogeographic history of Anarthrophyllum resulted from the action of ecological, historical, and geographic factors that acted either alternatively or simultaneously.     

A phylogenetic perspective on the evolution of Mediterranean teleost fishes

The Mediterranean Sea is a highly diverse, highly studied, and highly impacted biogeographic region, yet no phylogenetic reconstruction of fish diversity in this area has been published to date. Here, we infer the timing and geographic origins of Mediterranean teleost species diversity using nucleotide sequences collected from GenBank. We assembled a DNA supermatrix composed of four mitochondrial genes (12S ribosomal DNA, 16S ribosomal DNA, cytochrome c oxidase subunit I and cytochrome b) and two nuclear genes (rhodopsin and recombination activating gene I), including 62% of Mediterranean teleost species plus 9 outgroups. Maximum likelihood and Bayesian phylogenetic and dating analyses were calibrated using 20 fossil constraints. An additional 124 species were grafted onto the chronogram according to their taxonomic affinity, checking for the effects of taxonomic coverage in subsequent diversification analyses. We then interpreted the time-line of teleost diversification in light of Mediterranean historical biogeography, distinguishing non-endemic natives, endemics and exotic species. Results show that the major Mediterranean orders are of Cretaceous origin, specifically ~100-80 Mya, and most Perciformes families originated 80-50 Mya. Two important clade origin events were detected. The first at 100-80 Mya, affected native and exotic species, and reflects a global diversification period at a time when the Mediterranean Sea did not yet exist. The second occurred during the last 50 Mya, and is noticeable among endemic and native species, but not among exotic species. This period corresponds to isolation of the Mediterranean from Indo-Pacific waters before the Messinian salinity crisis. The Mediterranean fish fauna illustrates well the assembly of regional faunas through origination and immigration, where dispersal and isolation have shaped the emergence of a biodiversity hotspot.