PyRate: a new program to estimate speciation and extinction rates from incomplete fossil data

Despite the advancement of phylogenetic methods to estimate speciation and extinction rates, their power can be limited under variable rates, in particular for clades with high extinction rates and small number of extant species. Fossil data can provide a powerful alternative source of information to investigate diversification processes. Here, we present PyRate, a computer program to estimate speciation and extinction rates and their temporal dynamics from fossil occurrence data. The rates are inferred in a Bayesian framework and are comparable to those estimated from phylogenetic trees. We describe how PyRate can be used to explore different models of diversification. In addition to the diversification rates, it provides estimates of the parameters of the preservation process (fossilization and sampling) and the times of speciation and extinction of each species in the data set. Moreover, we develop a new birth-death model to correlate the variation of speciation/extinction rates with changes of a continuous trait. Finally, we demonstrate the use of Bayes factors for model selection and show how the posterior estimates of a PyRate analysis can be used to generate calibration densities for Bayesian molecular clock analysis. PyRate is an open-source command-line Python program available at http://sourceforge.net/projects/pyrate/.

Biogeographic analysis of Crocidiinae (Diptera, Bombyliidae): finding congruence among morphological, molecular, fossil and paleogeographical data

Biogeographic studies dealing with Bombyliidae are rare in the literature and no information is available on its origin and early diversification. In this study, we found evidence from molecular phylogeny and from fossil record supporting a Middle Jurassic origin of the Bombylioidea, taken as a starting point to discuss the biogeography and diversification of Crocidiinae. Based on a previously published phylogenetic hypothesis, we performed a Brooks Parsimony Analysis (BPA) to discuss the biogeographical history of Crocidiinae lineages. This subfamily is mostly distributed over arid areas of the early components of the Gondwanaland: Chile and southern Africa, but also in southwestern Palaearctic and southwestern Nearctic. The vicariant events affecting the Crocidiinae biogeography at the generic level seems to be related to the sequential separation of a Laurasian clade from a Gondwanan clade followed by the splitting of the latter into smaller components. This also leads to a hypothesis of origin of the Crocidiinae in the Middle Jurassic, the same period in which other bombyliid lineages are supposed to have arisen and irradiated.

Fossil resource depletion and climate change emissions: the role of physical constrictions

[spa] En lo que concierne al cambio climático, los pronósticos de cercanos picos de combustible fósiles parecen buenas noticias pues la mayoría de las emisiones proceden de la quema de combustibles fósiles. Sin embargo, esto podría resultar engañoso de confirmarse las enormes estimaciones de reservas de carbón pues puede divisarse un intercambio de combustible fósiles con baja concentración de carbono (petróleo y gas) por otros de mayor (carbón). Ciñéndonos a esta hipótesis desarrollamos escenarios donde tan pronto el petróleo y el gas natural alcanzan su cénit la extracción de carbón crece lo necesario para compensar el descenso de los primeros. Estimamos las emisiones que se deriva de tales supuestos y las comparamos con el peor escenario del IPCC. Si bien dicho escenario parece improbable concluimos que los picos de petróleo y gas no son suficientes para evitar peligrosas sendas de gases de efecto invernadero. Las concentraciones de CO2 halladas superan con creces las 450 ppm sin signos de remisión.

The fossil record and evolution of freshwater plants: A review

Palaeobotany applied to freshwater plants is an emerging field of palaeontology. Hydrophytic plants reveal evolutionary trends of their own, clearly distinct from those of the terrestrial and marine flora. During the Precambrian, two groups stand out in the fossil record of freshwater plants: the Cyanobacteria (stromatolites) in benthic environments and the prasinophytes (leiosphaeridian acritarchs) in transitional planktonic environments. During the Palaeozoic, green algae (Chlorococcales, Zygnematales, charophytes and some extinct groups) radiated and developed the widest range of morphostructural patterns known for these groups. Between the Permian and Early Cretaceous, charophytes dominated macrophytic associations, with the consequence that over tens of millions of years, freshwater flora bypassed the dominance of vascular plants on land. During the Early Cretaceous, global extension of the freshwater environments is associated with diversification of the flora, including new charophyte families and the appearance of aquatic angiosperms and ferns for the first time. Mesozoic planktonic assemblages retained their ancestral composition that was dominated by coenobial Chlorococcales, until the appearance of freshwater dinoflagellates in the Early Cretaceous. In the Late Cretaceous, freshwater angiosperms dominated almost all macrophytic communities worldwide. The Tertiary was characterised by the diversification of additional angiosperm and aquatic fern lineages, which resulted in the first differentiation of aquatic plant biogeoprovinces. Phytoplankton also diversified during the Eocene with the development of freshwater diatoms and chrysophytes. Diatoms, which were exclusively marine during tens of millions of years, were dominant over the Chlorococcales during Neogene and in later assemblages. During the Quaternary, aquatic plant communities suffered from the effects of eutrophication, paludification and acidification, which were the result of the combined impact of glaciation and anthropogenic disturbance.

Unusual concentration of Early Albian arthropod-bearing amber in the Basque-Cantabrian Basin (El Soplao, Cantabria, Spain): Palaeoenvironmental and paleobiological implications.

The El Soplao site is a recently-discovered Early Albian locality of the Basque-Cantabrian Basin (northern Spain) that has yielded a number of amber pieces with abundant bioinclusions. The amber-bearing deposit occurs in a non-marine to transitional marine siliciclastic unit (Las Peñosas Formation) that is interleaved within a regressive-transgressive, carbonate-dominated Lower Aptian-Upper Albian marine sequence. The Las Peñosas Formation corresponds to the regressive stage of this sequence and in its turn it splits into two smaller regressive-transgressive cycles. The coal and amber-bearing deposits occur in deltaic-estuarine environments developed during the maximum regressive episodes of these smaller regressive-transgressive cycles. The El Soplao amber shows Fourier Transform Infrared Spectroscopy spectra similar to other Spanish Cretaceous ambers and it is characterized by the profusion of sub-aerial, stalactite-like flows. Well-preserved plant cuticles assigned to the conifer genera Frenelopsis and Mirovia are abundant in the beds associated with amber. Leaves of the ginkgoalean genera Nehvizdya and Pseudotorellia also occur occasionally. Bioinclusions mainly consist of fossil insects of the orders Blattaria, Hemiptera, Thysanoptera, Raphidioptera, Neuroptera, Coleoptera, Hymenoptera and Diptera, although some spiders and spider webs have been observed as well. Some insects belong to groups scarce in the fossil record, such as a new morphotype of the wasp Archaeromma (of the family Mymarommatidae) and the biting midge Lebanoculicoides (of the monogeneric subfamily Lebanoculicoidinae). This new amber locality constitutes a very significant finding that will contribute to improving the knowledge and comprehension of the Albian non-marine paleoarthropod fauna.

Fossil resource depletion and climate change emissions: the role of physical constrictions

[spa] En lo que concierne al cambio climático, los pronósticos de cercanos picos de combustible fósiles parecen buenas noticias pues la mayoría de las emisiones proceden de la quema de combustibles fósiles. Sin embargo, esto podría resultar engañoso de confirmarse las enormes estimaciones de reservas de carbón pues puede divisarse un intercambio de combustible fósiles con baja concentración de carbono (petróleo y gas) por otros de mayor (carbón). Ciñéndonos a esta hipótesis desarrollamos escenarios donde tan pronto el petróleo y el gas natural alcanzan su cénit la extracción de carbón crece lo necesario para compensar el descenso de los primeros. Estimamos las emisiones que se deriva de tales supuestos y las comparamos con el peor escenario del IPCC. Si bien dicho escenario parece improbable concluimos que los picos de petróleo y gas no son suficientes para evitar peligrosas sendas de gases de efecto invernadero. Las concentraciones de CO2 halladas superan con creces las 450 ppm sin signos de remisión.

Bayesian estimation of speciation and extinction from incomplete fossil occurrence data.

The temporal dynamics of species diversity are shaped by variations in the rates of speciation and extinction, and there is a long history of inferring these rates using first and last appearances of taxa in the fossil record. Understanding diversity dynamics critically depends on unbiased estimates of the unobserved times of speciation and extinction for all lineages, but the inference of these parameters is challenging due to the complex nature of the available data. Here, we present a new probabilistic framework to jointly estimate species-specific times of speciation and extinction and the rates of the underlying birth-death process based on the fossil record. The rates are allowed to vary through time independently of each other, and the probability of preservation and sampling is explicitly incorporated in the model to estimate the true lifespan of each lineage. We implement a Bayesian algorithm to assess the presence of rate shifts by exploring alternative diversification models. Tests on a range of simulated data sets reveal the accuracy and robustness of our approach against violations of the underlying assumptions and various degrees of data incompleteness. Finally, we demonstrate the application of our method with the diversification of the mammal family Rhinocerotidae and reveal a complex history of repeated and independent temporal shifts of both speciation and extinction rates, leading to the expansion and subsequent decline of the group. The estimated parameters of the birth-death process implemented here are directly comparable with those obtained from dated molecular phylogenies. Thus, our model represents a step towards integrating phylogenetic and fossil information to infer macroevolutionary processes.

Nd and Sr isotope compositions in modern and fossil bones - Proxies for vertebrate provenance and taphonomy

Rare earth elements (REE), while not essential for the physiologic functions of animals, are ingested and incorporated in ppb concentrations in bones and teeth. Nd isotope compositions of modern bones of animals from isotopically distinct habitats demonstrate that the (143)Nd/(144)Nd of the apatite can be used as a fingerprint for bedrock geology or ambient water mass. This potentially allows the provenance and migration of extant vertebrates to be traced, similar to the use of Sr isotopes. Although REE may be enriched by up to 5 orders of magnitude during diagenesis and recrystallization of bone apatite, in vivo (143)Nd/(144)Nd may be preserved in the inner cortex of fossil bones or enamel. However, tracking the provenance of ancient or extinct vertebrates is possible only for well-preserved archeological and paleontological skeletal remains with in vivo-like Nd contents at the ppb-level. Intra-bone and -tooth REE analysis can be used to screen for appropriate areas. Large intra-bone Nd concentration gradients of 10(1)-10(3) are often measured. Nd concentrations in the inner bone cortex increase over timescales of millions of years, while bone rims may be enriched over millenial timescales. Nevertheless, epsilon(Nd) values are often similar within one epsilon(Nd) unit within a single bone. Larger intra-bone differences in specimens may either reflect a partial preservation of in vivo values or changing epsilon(Nd) values of the diagenetic fluid during fossilization. However, most fossil specimens and the outer rims of bones will record taphonomic (143)Nd/(144)Nd incorporated post mortem during diagenesis. Unlike REE patterns, (143)Nd/(144)Nd are not biased by fractionation processes during REE-uptake into the apatite crystal lattice, hence the epsilon(Nd) value is an important tracer for taphonomy and reworking. Bones and teeth from autochthonous fossil assemblages have small variations of +/- 1 epsilon(Nd) unit only. In contrast, fossil bones and teeth from over 20 different marine and terrestrial fossil sites have a total range of epsilon(Nd) values from -13.0 to 4.9 (n = 80), often matching the composition of the embedding sediment. This implies that the surrounding sediment is the source of Nd in the fossil bones and that the specimens of this study seem not to have been reworked. Differences in epsilon(Nd) values between skeletal remains and embedding sediment may either indicate reworking of fossils and/or a REE-uptake from a diagenetic fluid with non-sediment derived epsilon(Nd) values. The latter often applies to fossil shark teeth, which may preserve paleo-seawater values. Complementary to epsilon(Nd) values, (87)Sr/(86)Sr can help to further constrain the fossil provenance and reworking. (C) 2011 Elsevier Ltd. All rights reserved.

The odd couple: contrasting phylogeographic patterns in two sympatric sibling species of woodlouse-­‐hunter spiders in the Canary Islands

Comparative phylogeography seeks for commonalities in the spatial demographic history of sympatric organisms to characterize the mechanisms that shaped such patterns. The unveiling of incongruent phylogeographic patterns in co-occurring species, on the other hand, may hint to overlooked differences in their life histories or microhabitat preferences. The woodlouse-hunter spiders of the genus Dysdera have undergone a major diversi cation on the Canary Islands. The species pair Dysdera alegranzaensis and Dysdera nesiotes are endemic to the island of Lanzarote and nearby islets, where they co-occur at most of their known localities. The two species stand in sharp contrast to other sympatric endemic Dysdera in showing no evidence of somatic (non-genitalic) differentiation. Phylogenetic and population genetic analyses of mitochondrial cox1 sequences from an exhaustive sample of D. alegranzaensis and D. nesiotes specimens, and additional mitochondrial (16S, L1, nad1) and nuclear genes (28S, H3) were analysed to reveal their phylogeographic patterns and clarify their phylogenetic relationships. Relaxed molecular clock models using ve calibration points were further used to estimate divergence times between species and populations. Striking differences in phylogeography and population structure between the two species were observed. Dysdera nesiotes displayed a metapopulation-like structure, while D. alegranzaensis was characterized by a weaker geographical structure but greater genetic divergences among its main haplotype lineages, suggesting more complex population dynamics. Our study con rms that co-distributed sibling species may exhibit contrasting phylogeographic patterns in the absence of somatic differentiation. Further ecological studies, however, will be necessary to clarify whether the contrasting phylogeographies may hint at an overlooked niche partitioning between the two species. In addition, further comparisons with available phylogeographic data of other eastern Canarian Dysdera endemics con rm the key role of lava ows in structuring local populations in oceanic islands and identify localities that acted as refugia during volcanic eruptions

The Imprint of Geologic History on Within‐Island Diversification of Woodlouse-­Hunter Spiders (Araneae, Dysderidae) in the Canary Islands.

Geological processes and ecological adaptation are major drivers of diversification on oceanic islands. Although diversification in these islands is often interpreted as resulting from dispersal or island hopping rather than vicariance, this may not be the case in islands with complex geological histories. The island of Tenerife, in the Canary Islands, emerged in the late Miocene as 3 precursor islands that were subsequently connected and reisolated by volcanic cycles. The spider Dysdera verneaui is endemic to the island of Tenerife, where it is widely distributed throughout most island habitats, providing an excellent model to investigate the role of physical barriers and ecological adaptation in shaping within-island diversity. Here, we present evidence that the phylogeographic patterns of this species trace back to the independent emergence of the protoislands. Molecular markers (mitochondrial genes cox1, 16S, and nad1 and the nuclear genes ITS-2 and 28S) analyzed from 100 specimens (including a thorough sampling of D. verneaui populations and additional outgroups) identify 2 distinct evolutionary lineages that correspond to 2 precursor islands, each with diagnostic genital characters indicative of separate species status. Episodic introgression events between these 2 main evolutionary lineages explain the observed incongruence between mitochondrial and nuclear markers, probably as a result of the homogenization of their ITS-2 sequence types. The most widespread lineage exhibits a complex population structure, which is compatible with either secondary contact, following connection of deeply divergent lineages, or alternatively, a back colonization from 1 precursor island to another.

The Imprint of Geologic History on Within‐Island Diversification of Woodlouse-­Hunter Spiders (Araneae, Dysderidae) in the Canary Islands.

Geological processes and ecological adaptation are major drivers of diversification on oceanic islands. Although diversification in these islands is often interpreted as resulting from dispersal or island hopping rather than vicariance, this may not be the case in islands with complex geological histories. The island of Tenerife, in the Canary Islands, emerged in the late Miocene as 3 precursor islands that were subsequently connected and reisolated by volcanic cycles. The spider Dysdera verneaui is endemic to the island of Tenerife, where it is widely distributed throughout most island habitats, providing an excellent model to investigate the role of physical barriers and ecological adaptation in shaping within-island diversity. Here, we present evidence that the phylogeographic patterns of this species trace back to the independent emergence of the protoislands. Molecular markers (mitochondrial genes cox1, 16S, and nad1 and the nuclear genes ITS-2 and 28S) analyzed from 100 specimens (including a thorough sampling of D. verneaui populations and additional outgroups) identify 2 distinct evolutionary lineages that correspond to 2 precursor islands, each with diagnostic genital characters indicative of separate species status. Episodic introgression events between these 2 main evolutionary lineages explain the observed incongruence between mitochondrial and nuclear markers, probably as a result of the homogenization of their ITS-2 sequence types. The most widespread lineage exhibits a complex population structure, which is compatible with either secondary contact, following connection of deeply divergent lineages, or alternatively, a back colonization from 1 precursor island to another.

The Imprint of Geologic History on Within‐Island Diversification of Woodlouse-­Hunter Spiders (Araneae, Dysderidae) in the Canary Islands.

Geological processes and ecological adaptation are major drivers of diversification on oceanic islands. Although diversification in these islands is often interpreted as resulting from dispersal or island hopping rather than vicariance, this may not be the case in islands with complex geological histories. The island of Tenerife, in the Canary Islands, emerged in the late Miocene as 3 precursor islands that were subsequently connected and reisolated by volcanic cycles. The spider Dysdera verneaui is endemic to the island of Tenerife, where it is widely distributed throughout most island habitats, providing an excellent model to investigate the role of physical barriers and ecological adaptation in shaping within-island diversity. Here, we present evidence that the phylogeographic patterns of this species trace back to the independent emergence of the protoislands. Molecular markers (mitochondrial genes cox1, 16S, and nad1 and the nuclear genes ITS-2 and 28S) analyzed from 100 specimens (including a thorough sampling of D. verneaui populations and additional outgroups) identify 2 distinct evolutionary lineages that correspond to 2 precursor islands, each with diagnostic genital characters indicative of separate species status. Episodic introgression events between these 2 main evolutionary lineages explain the observed incongruence between mitochondrial and nuclear markers, probably as a result of the homogenization of their ITS-2 sequence types. The most widespread lineage exhibits a complex population structure, which is compatible with either secondary contact, following connection of deeply divergent lineages, or alternatively, a back colonization from 1 precursor island to another.

The odd couple: contrasting phylogeographic patterns in two sympatric sibling species of woodlouse-­‐hunter spiders in the Canary Islands

Comparative phylogeography seeks for commonalities in the spatial demographic history of sympatric organisms to characterize the mechanisms that shaped such patterns. The unveiling of incongruent phylogeographic patterns in co-occurring species, on the other hand, may hint to overlooked differences in their life histories or microhabitat preferences. The woodlouse-hunter spiders of the genus Dysdera have undergone a major diversi cation on the Canary Islands. The species pair Dysdera alegranzaensis and Dysdera nesiotes are endemic to the island of Lanzarote and nearby islets, where they co-occur at most of their known localities. The two species stand in sharp contrast to other sympatric endemic Dysdera in showing no evidence of somatic (non-genitalic) differentiation. Phylogenetic and population genetic analyses of mitochondrial cox1 sequences from an exhaustive sample of D. alegranzaensis and D. nesiotes specimens, and additional mitochondrial (16S, L1, nad1) and nuclear genes (28S, H3) were analysed to reveal their phylogeographic patterns and clarify their phylogenetic relationships. Relaxed molecular clock models using ve calibration points were further used to estimate divergence times between species and populations. Striking differences in phylogeography and population structure between the two species were observed. Dysdera nesiotes displayed a metapopulation-like structure, while D. alegranzaensis was characterized by a weaker geographical structure but greater genetic divergences among its main haplotype lineages, suggesting more complex population dynamics. Our study con rms that co-distributed sibling species may exhibit contrasting phylogeographic patterns in the absence of somatic differentiation. Further ecological studies, however, will be necessary to clarify whether the contrasting phylogeographies may hint at an overlooked niche partitioning between the two species. In addition, further comparisons with available phylogeographic data of other eastern Canarian Dysdera endemics con rm the key role of lava ows in structuring local populations in oceanic islands and identify localities that acted as refugia during volcanic eruptions