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/.

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.

Using a Combination of Phylogenetic Comparative and Paleobotanical Methods to Elucidate Patterns of Lineage Selection in Rosales (Plantae: Anthophyta)

An increasing focus in evolutionary biology is on the interplay between mesoscale ecological and evolutionary processes such as population demographics, habitat tolerance, and especially geographic distribution, as potential drivers responsible for patterns of diversification and extinction over geologic time. However, few studies to date connect organismal processes such as survival and reproduction through mesoscale patterns to long-term macroevolutionary trends. In my dissertation, I investigate how mechanism of seed dispersal, mediated through geographic range size, influences diversification rates in the Rosales (Plantae: Anthophyta). In my first chapter, I validate the phylogenetic comparative methods that I use in my second and third chapters. Available state speciation and extinction (SSE) models assumptions about evolution known to be false through fossil data. I show, however, that as long as net diversification rates remain positive – a condition likely true for the Rosales – these violations of SSE’s assumptions do not cause significantly biased results. With SSE methods validated, my second chapter reconstructs three associations that appear to increase diversification rate for Rosalean genera: (1) herbaceous habit; (2) a three-way interaction combining animal dispersal, high within-genus species richness, and geographic range on multiple continents; (3) a four-way interaction combining woody habit with the other three characteristics of (2). I suggest that the three- and four-way interactions represent colonization ability and resulting extinction resistance in the face of late Cenozoic climate change; however, there are other possibilities as well that I hope to investigate in future research. My third chapter reconstructs the phylogeographic history of the Rosales using both non-fossil-assisted SSE methods as well as fossil-informed traditional phylogeographic analysis. Ancestral state reconstructions indicate that the Rosaceae diversified in North America while the other Rosalean families diversified elsewhere, possibly in Eurasia. SSE is able to successfully identify groups of genera that were likely to have been ancestrally widespread, but has poorer taxonomic resolution than methods that use fossil data. In conclusion, these chapters together suggest several potential causal links between organismal, mesoscale, and geologic scale processes, but further work will be needed to test the hypotheses that I raise here.

Disentangling the effects of key innovations on the diversification of Bromelioideae (bromeliaceae).

The evolution of key innovations, novel traits that promote diversification, is often seen as major driver for the unequal distribution of species richness within the tree of life. In this study, we aim to determine the factors underlying the extraordinary radiation of the subfamily Bromelioideae, one of the most diverse clades among the neotropical plant family Bromeliaceae. Based on an extended molecular phylogenetic data set, we examine the effect of two putative key innovations, that is, the Crassulacean acid metabolism (CAM) and the water-impounding tank, on speciation and extinction rates. To this aim, we develop a novel Bayesian implementation of the phylogenetic comparative method, binary state speciation and extinction, which enables hypotheses testing by Bayes factors and accommodates the uncertainty on model selection by Bayesian model averaging. Both CAM and tank habit were found to correlate with increased net diversification, thus fulfilling the criteria for key innovations. Our analyses further revealed that CAM photosynthesis is correlated with a twofold increase in speciation rate, whereas the evolution of the tank had primarily an effect on extinction rates that were found five times lower in tank-forming lineages compared to tank-less clades. These differences are discussed in the light of biogeography, ecology, and past climate change.

Immigration, species radiation and extinction in a highly diverse songbird lineage: white-eyes on Indian Ocean islands

Molecular phylogenetic hypotheses of species-rich lineages in regions where geological history can be reliably inferred may provide insights into the scale of processes driving diversification. Here we sample all extant or recently extinct white-eye (Zosterops) taxa of the southwest Indian Ocean, combined with samples from all principal continental lineages. Results support a high dispersal capability, with at least two independent continental sources for white-eyes of the region. An early (within 1.8 million years ago) expansion into the Indian Ocean may have originated either from Asia or Africa; the three resulting lineages show a disparate distribution consistent with considerable extinction following their arrival. Africa is supported as the origin of a later expansion into the region (within 1.2 million years ago). On two islands, a pair of Zosterops species derived from independent immigrations into the Indian Ocean co-occur or may have formerly co-occurred, providing strong support for their origin by double-island colonization rather than within-island (sympatric or microallopatric) speciation. On Mauritius and La Reunion, phylogenetic placement of sympatric white-eyes allow us to rule out a scenario in which independent within-island speciation occurred on both islands; one of the species pairs must have arisen by double colonization, while the other pair is likely to have arisen by the same mechanism. Long-distance immigration therefore appears to be responsible for much of the region's white-eye diversity. Independent immigrations into the region have resulted in lineages with mutually exclusive distributions and it seems likely that competition with congeneric species, rather than arrival frequency, may limit present-day diversity.

A Mismatch-Based Model for Memory Reconsolidation and Extinction in Attractor Networks

OSAN, R. , TORT, A. B. L. , AMARAL, O. B. . A mismatch-based model for memory reconsolidation and extinction in attractor networks. Plos One, v. 6, p. e23113, 2011.

Iron speciation and magnetic susceptibility in bottom sediments of the Black Sea, Core SMG91-903

12 cores of Late Pleistocene - Holocene deposits were studied. They were collected by gravity cores on the continental slope and in the deep-water part of the Black Sea within the Adler-Tuapse polygon. In four of them in New Euxinian deposits at the base of a packet of hydrotroilite laminae paleomagnetic anomalies likely resulting from the Gothenburg magnetic excursion occur. Comparison with results of similar studies in the western Black Sea, where the Gothenburg magnetic excursion was previously found, let to validate stratigraphic synchronism of the hydrotroilite horizon in the eastern and western parts of the Black Sea and to confirm the authors' views about peculiarities of paleogeographical development of the Black Sea basin in the Late Pleistocene - Holocene.

A Mismatch-Based Model for Memory Reconsolidation and Extinction in Attractor Networks

OSAN, R. , TORT, A. B. L. , AMARAL, O. B. . A mismatch-based model for memory reconsolidation and extinction in attractor networks. Plos One, v. 6, p. e23113, 2011.

A Mismatch-Based Model for Memory Reconsolidation and Extinction in Attractor Networks

OSAN, R. , TORT, A. B. L. , AMARAL, O. B. . A mismatch-based model for memory reconsolidation and extinction in attractor networks. Plos One, v. 6, p. e23113, 2011.

Phylogenetics, Ancestral State Reconstruction, And A New Infrafamilial Classification Of The Pantropical Ochnaceae (medusagynaceae, Ochnaceae S.str., Quiinaceae) Based On Five Dna Regions.

Ochnaceae s.str. (Malpighiales) are a pantropical family of about 500 species and 27 genera of almost exclusively woody plants. Infrafamilial classification and relationships have been controversial partially due to the lack of a robust phylogenetic framework. Including all genera except Indosinia and Perissocarpa and DNA sequence data for five DNA regions (ITS, matK, ndhF, rbcL, trnL-F), we provide for the first time a nearly complete molecular phylogenetic analysis of Ochnaceae s.l. resolving most of the phylogenetic backbone of the family. Based on this, we present a new classification of Ochnaceae s.l., with Medusagynoideae and Quiinoideae included as subfamilies and the former subfamilies Ochnoideae and Sauvagesioideae recognized at the rank of tribe. Our data support a monophyletic Ochneae, but Sauvagesieae in the traditional circumscription is paraphyletic because Testulea emerges as sister to the rest of Ochnoideae, and the next clade shows Luxemburgia+Philacra as sister group to the remaining Ochnoideae. To avoid paraphyly, we classify Luxemburgieae and Testuleeae as new tribes. The African genus Lophira, which has switched between subfamilies (here tribes) in past classifications, emerges as sister to all other Ochneae. Thus, endosperm-free seeds and ovules with partly to completely united integuments (resulting in an apparently single integument) are characters that unite all members of that tribe. The relationships within its largest clade, Ochnineae (former Ochneae), are poorly resolved, but former Ochninae (Brackenridgea, Ochna) are polyphyletic. Within Sauvagesieae, the genus Sauvagesia in its broad circumscription is polyphyletic as Sauvagesia serrata is sister to a clade of Adenarake, Sauvagesia spp., and three other genera. Within Quiinoideae, in contrast to former phylogenetic hypotheses, Lacunaria and Touroulia form a clade that is sister to Quiina. Bayesian ancestral state reconstructions showed that zygomorphic flowers with adaptations to buzz-pollination (poricidal anthers), a syncarpous gynoecium (a near-apocarpous gynoecium evolved independently in Quiinoideae and Ochninae), numerous ovules, septicidal capsules, and winged seeds with endosperm are the ancestral condition in Ochnoideae. Although in some lineages poricidal anthers were lost secondarily, the evolution of poricidal superstructures secured the maintenance of buzz-pollination in some of these genera, indicating a strong selective pressure on keeping that specialized pollination system.

Colonization and extinction of land planarians (Platyhelminthes, Tricladida) in a Brazilian Atlantic Forest regrowth remnant

Long-term assessments of species assemblages are valuable tools for detecting species ecological preferences and their dispersal tracks, as well as for assessing the possible effects of alien species on native communities. Here we report a 50-year-long study on population dynamics of the four species of land flatworms (Platyhelminthes, Tricladida, Terricola) that have colonized or become extinct in a 70-year-old Atlantic Forest regrowth remnant through the period 1955-2006. On the one hand, the two initially most abundant species, which are native to the study site, Notogynaphallia ernesti and Geoplana multicolor have declined over decades and at present do not exist in the forest remnant. The extinction of these species is most likely related with their preference for open vegetation areas, which presently do not exist in the forest remnant. On the other hand, the neotropical Geoplaninae 1 and the exotic Endeavouria septemlineata were detected in the forest only very recently. The long-term study allowed us to conclude that Geoplaninae 1 was introduced into the study area, although it is only known from the study site. Endeavouria septemlineata, an active predator of the exotic giant African snail, is originally known from Hawaii. This land flatworm species was observed repeatedly in Brazilian anthropogenic areas, and this is the first report of the species in relatively well preserved native forest, which may be evidence of an ongoing adaptive process. Monitoring of its geographic spread and its ecological role would be a good practice for preventing potential damaging effects, since it also feeds on native mollusk fauna, as we observed in lab conditions.

Effects of geological inhomogeneity on high Rayleigh number steady state heat and mass transfer in fluid-saturated porous media heated from below

A parametric study is carried out to investigate how geological inhomogeneity affects the pore-fluid convective flow field, the temperature distribution, and the mass concentration distribution in a fluid-saturated porous medium. The related numerical results have demonstrated that (1) the effects of both medium permeability inhomogeneity and medium thermal conductivity inhomogeneity are significant on the pore-fluid convective flow and the species concentration distribution in the porous medium; (2) the effect of medium thermal conductivity inhomogeneity is dramatic on the temperature distribution in the porous medium, but the effect of medium permeability inhomogeneity on the temperature distribution may be considerable, depending on the Rayleigh number involved in the analysis; (3) if the coupling effect between pore-fluid flow and mass transport is weak, the effect of the Lewis number is negligible on the pore-fluid convective flow and temperature distribution, hut it is significant on the species concentration distribution in the medium.

Speciation and absolute bioavailability: risk assessment of arsenic-contaminated sites in a residential suburb in Canberra

Watson is a fully developed suburb of some 30 years in Canberra (the capital city of Australia), A plunge dip using arsenical pesticides for tick control was operated there between 1946 and 1960, Chemical investigations revealed that many soil samples obtained from the study area contained levels of arsenic exceeding the current health-based investigation levels of 100 mg kg(-1) set by the National Health and Medical Research Council in Australia, For the speciation study, nine composite samples of surface and sub-surface soils and a composite sample of rocks were selected. ICP-MS analysis showed that arsenic levels in these samples ranged from 32 to 1597 mg kg(-1), Chemical speciation of arsenic showed that the arsenite (trivalent) components were 0.32-56% in the soil and 44.8% in the rock composite samples. Using a rat model, the absolute bioavailability of these contaminated soils relative to As3+ or As5+ ranged from 1.02 to 9.87% and 0.26 to 2.98%, respectively, An attempt was made to develop a suitable leachate test as an index of bioavailability. However, the results indicated that there was no significant correlation between the bioavailability and leachates using neutral pH water or 1 M HCl. Our results indicate that speciation is highly significant for the interpretation of bioavailability and risk assessment data; the bioavailable fractions of arsenic in soils from Watson are small and therefore the health impact upon the environment and humans due to this element is limited.

Mesoionic pyridopyrimidinylium and pyridooxazinylium olates and non-mesoionic pyridopyrimidinones. Structures in the solid state, solution, and matrices

X-Ray crystal structures, C-13 NMR spectra and theoretical calculations (B3LYP/6-31G*) are reported for the mesoionic (zwitterionic) pyridopyrimidinylium- and pyridooxazinyliumolates 2a, 3a and 5a,b as well as the enol ether 11b and the enamine 11c. The 1-NH compounds like 1a, 2a and 3a exist in the mesoionic form in the crystal and in solution, but the OH tautomers such as 1b and 2b dominate in the gas phase as revealed by the Ar matrix IR spectra in conjunction with DFT calculations. All data indicate that the mesoionic compounds can be regarded as intramolecular pyridine-ketene zwitterions (cf. 16 --> 17) with a high degree of positive charge on the pyridinium nitrogen, a long pyridinium N-CO bond (ca. 1.44-1.49 Angstrom), and normal C=O double bonds (ca. 1.22 Angstrom). All mesoionic compounds exhibit a pronounced tilting of the olate C=O groups (the C=O groups formally derived from a ketene) towards the pyridinium nitrogen, giving NCO angles of 110-118 degrees. Calculations reveal a hydrogen bond with 6-CH, analogous to what is found in ketene-pyridine zwitterions and the C3O2-pyridine complex. The 2-OH tautomers of type 1b, 2b, and 11 also show a high degree of zwitterionic character as indicated by the canonical structures 11 12.