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

The role of climatic tolerances and seed traits in reduced extinction rates of temperate polygonaceae.

The latitudinal diversity gradient (LDG) is one of the most striking and consistent biodiversity patterns across taxonomic groups. We investigate the species richness gradient in the buckwheat family, Polygonaceae, which exhibits a reverse LDG and is, thus, decoupled from dominant gradients of energy and environmental stability that increase toward the tropics and confound mechanistic interpretations. We test competing age and evolutionary diversification hypotheses, which may explain the diversification of this plant family over the past 70 million years. Our analyses show that the age hypothesis, which posits that clade richness is positively correlated with the ecological and evolutionary time since clade origin, fails to explain the richness gradient observed in Polygonaceae. However, an evolutionary diversification hypothesis is highly supported, with diversification rates being 3.5 times higher in temperate clades compared to tropical clades. We demonstrate that differences in rates of speciation, migration, and molecular evolution insufficiently explain the observed patterns of differential diversification rates. We suggest that reduced extinction rates in temperate clades may be associated with adaptive responses to selection, through which seed morphology and climatic tolerances potentially act to minimize risk in temporally variable environments. Further study is needed to understand causal pathways among these traits and factors correlated with latitude.

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.

A global marine sedimentary response to the end-Permian mass extinction: Examples from southern Turkey and the western United States

The end-Permian mass extinction greatly diminished marine diversity and brought about a whole-scale restructuring of marine ecosystems; these ecosystem changes also profoundly affected the sedimentary record. Data presented here, attained through facies analyses of strata deposited during the immediate aftermath of the end-Permian mass extinction (southern Turkey) and at the close of the Early Triassic (southwestern United States), in combination with a literature review, show that sedimentary systems were profoundly affected by: (1) a reduction in biotic diversity and abundance and (2) long-term environmental fluctuations that resulted from the end-Permian crisis. Lower Triassic strata display widespread microbialite and carbonate seafloor fan development and contain indicators of suppressed infaunal bioturbation such as flat-pebble conglomerates and wrinkle structures (facies considered unusual in post-Cambrian subtidal deposits). Our observations suggest that depositional systems, too, respond to biotic crises, and that certain facies may act as barometers of ecologic and environmental change independent of fossil assemblage analyses. Close investigation of facies changes during other critical times in Earth history may serve as an important tool in interpreting the ecology of metazoans and their environment.

Analyse des risques d'extinction des populations naturelles

A review of extinction risk analysis and viability methods is presented. The importance of environmental, demographic and genetic uncertainties, as well as the role of catastrophes are successively considered, and different approaches aiming at the integration of these risk factors in predictive population dynamic models are discussed.

Structure du temps de travail des infirmières de santé publique: effets de l'âge des patients, de la nature des prestations fournies et du lieu du contact. [Structure of work times of public health nurses: effect of the age of patients, nature of care given and location of contact].

In the framework of health services research sponsored by the Swiss National Science Foundation, a research was undertaken of the activity of the large majority of the public health nurses working in the Swiss cantons of Vaud and Fribourg (total population 700,000). During one week, 130 nurses gathered, with a specially devised instrument, data on 4165 patient visits. Studying the duration of the contacts, one has distinguished contact duration per se (DC), duration of the travel time preceding the contact (DD), and total duration in relation with the contact (DTC-addition of the first two). It was noted that the three durations increased significantly with patient age (as regard travel time, this is explained by the higher proportion of home visits in higher age groups, as compared with visits at a health center). Examined according to location of the visit, contact duration per se (without travel) is higher for visits at home and in nursing homes than for those taking place at a health center. Looked at in respect to the care given (technical care, or basic nursing care, or both simultaneously), our data show that the provision of basic nursing care (alone or with technical care) doubles contact duration (from 20 to 42-45'). The analyses according to patient age shows that, at an advanced age (beyond 80 years particularly), there is an important increase of the visits where both types of care are given. However, contact duration per se shows a significant raise with age only for the group "technical care only"; it can be demonstrated that this is due to the fact that older patients require more complex technical acts (e.g., bladder care, as compared with simpler acts such as injection). A model of the relationships between patient age and contact duration is proposed: it is because of the increase in the proportions of home visits, of visits including basic nursing care, and of more complex technical acts that older persons require more of the working time of public health nurses.

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.

MR spectroscopy of the human brain with enhanced signal intensity at ultrashort echo times on a clinical platform at 3T and 7T.

Recently, the spin-echo full-intensity acquired localized (SPECIAL) spectroscopy technique was proposed to unite the advantages of short TEs on the order of milliseconds (ms) with full sensitivity and applied to in vivo rat brain. In the present study, SPECIAL was adapted and optimized for use on a clinical platform at 3T and 7T by combining interleaved water suppression (WS) and outer volume saturation (OVS), optimized sequence timing, and improved shimming using FASTMAP. High-quality single voxel spectra of human brain were acquired at TEs below or equal to 6 ms on a clinical 3T and 7T system for six volunteers. Narrow linewidths (6.6 +/- 0.6 Hz at 3T and 12.1 +/- 1.0 Hz at 7T for water) and the high signal-to-noise ratio (SNR) of the artifact-free spectra enabled the quantification of a neurochemical profile consisting of 18 metabolites with Cramér-Rao lower bounds (CRLBs) below 20% at both field strengths. The enhanced sensitivity and increased spectral resolution at 7T compared to 3T allowed a two-fold reduction in scan time, an increased precision of quantification for 12 metabolites, and the additional quantification of lactate with CRLB below 20%. Improved sensitivity at 7T was also demonstrated by a 1.7-fold increase in average SNR (= peak height/root mean square [RMS]-of-noise) per unit-time.

Correlating the end-Triassic mass extinction and flood basalt volcanism at the 100 ka level

New high-precision U/Pb geochronology from volcanic ashes shows that the Triassic-Jurassic boundary and end-Triassic biological crisis from two independent marine stratigraphic sections correlate with the onset of terrestrial flood volcanism in the Central Atlantic Magmatic Province to <150 ka. This narrows the correlation between volcanism and mass extinction by an order of magnitude for any such catastrophe in Earth history. We also show that a concomitant drop and rise in sea level and negative delta C-13 spike in the very latest Triassic occurred locally in <290 ka. Such rapid sea-level fluctuations on a global scale require that global cooling and glaciation were closely associated with the end-Triassic extinction and potentially driven by Central Atlantic Magmatic Province volcanism.

Constraints on the upper Permian to upper Triassic marine carbon isotope curve : case studies from the Tethys

Secular variations of the seawater carbon isotopic composition provide evidence for paleoceanographic and paleoclimatic changes and may serve for chemiostratigraphic correlations. The present study aimed to improve the current knowledge on the Upper Permian and Triassic segment of the Phanerozoic marine carbon isotope curve, whose Triassic part was poorly constrained by previous studies. Profiles of inorganic carbon isotopes are provided for sections from Himalaya (Salt Range, Kashmir, Spiti and Nepal), Oman and North Dobrogea (Romania) on the basis of whole-rock carbonate analysis. The data acquired, together with a literature compilation confirmed that most of the Upper Permian is characterized by high δ13C values (averaging +40/00) but failed to detect a positive excursion as suggested by recent compilations. In the light of these observations, the large drop in δ13C values associated with the end-Permian mass extinction appears to be driven by a sudden transfer of previously stocked 13C depleted carbon, rather than by the overturn of a Late Permian stratified ocean. The Triassic data-set outlines significant secular variations. The best documented is a carbon isotope positive excursion just across the Lower-Middle Triassic boundary, globally developed since it was detected in various paleogeographic settings. It is interpreted to reflect variations in surface ocean chemistry, possibly related to increased primary productivity, at times when the biotic recovery after the end-Permian mass-extinction began to accelerate significantly and when a sharp rise in seawater δ34S values occurred globally. Strontium isotope data obtained from well preserved biogenic phosphates allow a refinement of the Middle Triassic segment of the seawater strontium isotope curve and show a major inflexion point of the seawater strontium isotope curve also near the Lower Triassic - Middle Triassic boundary. These facts suggest that the transition from the Early to the Middle Triassic was a time of revolutionary global change which represented an important step in the evolution of Mesozoic marine environments. A tentative carbon isotope curve for the Upper Permian to Upper Triassic time interval is proposed. Its major features are: ? high but constant δ13C values during the Late Permian ? a sharp drop in δ13C values in the latest Permian ? subsequent recovery of δ13C values ? a short-lived positive excursion across the Early-Middle Triassic boundary ? a gradual rise in δ13C values starting in the Late Ladinian or in the Early Carnian It is foreseen that these fluctuations of the carbon isotope curve may serve as chronostratigraphic markers and further assist in the correlation of Permian and Triassic carbonate deposits.

A Note on improving Contact Times in Panel Surveys

With increasing costs for conducting surveys, many survey agencies resort to implementing call strategies. Obtaining contact in panel surveys as early as possible, without annoying people by contacting them at undesired times and ultimately causing them to refuse, requires using efficient call time strategies. In this research, the author uses call data from the Swiss Household Panel (SHP), a centralized Computer Assisted Telephone Interview (CATI) survey with a randomized (experimental) call-household assignment. Using random effects models, the author analyzes the efficiency gains of obtaining initial contact by assigning optimal times to first calls, and times and spacing to second and later calls depending on household sociodemography and prior call patterns. The author concludes with some recommendations for making early and successful contact during fieldwork.

Prolonged survival for patients with newly diagnosed, inoperable glioblastoma with 3-times daily ultrafractionated radiation therapy.

Ultrafractionation of radiation therapy is a novel regimen consisting of irradiating tumors several times daily, delivering low doses (<0.75 Gy) at which hyperradiosensitivity occurs. We recently demonstrated the high efficiency of ultrafractionated radiotherapy (RT) on glioma xenografts and report here on a phase II clinical trial to determine the safety, tolerability, and efficacy of an ultrafractionation regimen in patients with newly and inoperable glioblastoma (GBM). Thirty-one patients with histologically proven, newly diagnosed, and unresectable supratentorial GBM (WHO grade IV) were enrolled. Three daily doses of 0.75 Gy were delivered at least 4 hours apart, 5 days per week over 6-7 consecutive weeks (90 fractions for a total of 67.5 Gy). Conformal irradiation included the tumor bulk with a margin of 2.5 cm. The primary end points were safety, toxicity, and tolerability, and the secondary end points were overall survival (OS) and progression-free survival (PFS). Multivariate analysis was used to compare the OS and PFS with the EORTC-NCIC trial 26981-22981/CE.3 of RT alone vs radiation therapy and temozolomide (TMZ). The ultrafractionation radiation regimen was safe and well tolerated. No acute Grade III and/or IV CNS toxicity was observed. Median PFS and OS from initial diagnosis were 5.1 and 9.5 months, respectively. When comparing with the EORTC/NCIC trial, in both PFS and OS multivariate analysis, ultrafractionation showed superiority over RT alone, but not over RT and TMZ. The ultrafractionation regimen is safe and may prolong the survival of patients with GBM. Further investigation is warranted and a trial associating ultra-fractionation and TMZ is ongoing.

Comparison of T1 relaxation times of the neurochemical profile in rat brain at 9.4 tesla and 14.1 tesla.

Knowledge of T(1) relaxation times can be important for accurate relative and absolute quantification of brain metabolites, for sensitivity optimizations, for characterizing molecular dynamics, and for studying changes induced by various pathological conditions. (1)H T(1) relaxation times of a series of brain metabolites, including J-coupled ones, were determined using a progressive saturation (PS) technique that was validated with an adiabatic inversion-recovery (IR) method. The (1)H T(1) relaxation times of 16 functional groups of the neurochemical profile were measured at 14.1T and 9.4T. Overall, the T(1) relaxation times found at 14.1T were, within the experimental error, identical to those at 9.4T. The T(1)s of some coupled spin resonances of the neurochemical profile were measured for the first time (e.g., those of gamma-aminobutyrate [GABA], aspartate [Asp], alanine [Ala], phosphoethanolamine [PE], glutathione [GSH], N-acetylaspartylglutamate [NAAG], and glutamine [Gln]). Our results suggest that T(1) does not increase substantially beyond 9.4T. Furthermore, the similarity of T(1) among the metabolites (approximately 1.5 s) suggests that T(1) relaxation time corrections for metabolite quantification are likely to be similar when using rapid pulsing conditions. We therefore conclude that the putative T(1) increase of metabolites has a minimal impact on sensitivity when increasing B(0) beyond 9.4T.