Probabilistic base calling of Solexa sequencing data.

BACKGROUND: Solexa/Illumina short-read ultra-high throughput DNA sequencing technology produces millions of short tags (up to 36 bases) by parallel sequencing-by-synthesis of DNA colonies. The processing and statistical analysis of such high-throughput data poses new challenges; currently a fair proportion of the tags are routinely discarded due to an inability to match them to a reference sequence, thereby reducing the effective throughput of the technology. RESULTS: We propose a novel base calling algorithm using model-based clustering and probability theory to identify ambiguous bases and code them with IUPAC symbols. We also select optimal sub-tags using a score based on information content to remove uncertain bases towards the ends of the reads. CONCLUSION: We show that the method improves genome coverage and number of usable tags as compared with Solexa's data processing pipeline by an average of 15%. An R package is provided which allows fast and accurate base calling of Solexa's fluorescence intensity files and the production of informative diagnostic plots.

EUROCAT website data on prenatal detection rates of congenital anomalies.

The EUROCAT website www.eurocat-network.eu publishes prenatal detection rates for major congenital anomalies using data from European population-based congenital anomaly registers, covering 28% of the EU population as well as non-EU countries. Data are updated annually. This information can be useful for comparative purposes to clinicians and public health service managers involved in the antenatal care of pregnant women as well as those interested in perinatal epidemiology.

Estimation of electronic coupling in π -stacked donor-bridge-acceptor systems: correction of the two-state model

Comparison of donor-acceptor electronic couplings calculated within two-state and three-state models suggests that the two-state treatment can provide unreliable estimates of Vda because of neglecting the multistate effects. We show that in most cases accurate values of the electronic coupling in a π stack, where donor and acceptor are separated by a bridging unit, can be obtained as Ṽ da = (E2 - E1) μ12 Rda + (2 E3 - E1 - E2) 2 μ13 μ23 Rda2, where E1, E2, and E3 are adiabatic energies of the ground, charge-transfer, and bridge states, respectively, μij is the transition dipole moments between the states i and j, and Rda is the distance between the planes of donor and acceptor. In this expression based on the generalized Mulliken-Hush approach, the first term corresponds to the coupling derived within a two-state model, whereas the second term is the superexchange correction accounting for the bridge effect. The formula is extended to bridges consisting of several subunits. The influence of the donor-acceptor energy mismatch on the excess charge distribution, adiabatic dipole and transition moments, and electronic couplings is examined. A diagnostic is developed to determine whether the two-state approach can be applied. Based on numerical results, we showed that the superexchange correction considerably improves estimates of the donor-acceptor coupling derived within a two-state approach. In most cases when the two-state scheme fails, the formula gives reliable results which are in good agreement (within 5%) with the data of the three-state generalized Mulliken-Hush model

Anti-Nogo-A Antibody Treatment Promotes Recovery of Manual Dexterity after Unilateral Cervical Lesion in Adult Primates--re-examination and Extension of Behavioral Data.

In rodents and nonhuman primates subjected to spinal cord lesion, neutralizing the neurite growth inhibitor Nogo-A has been shown to promote regenerative axonal sprouting and functional recovery. The goal of the present report was to re-examine the data on the recovery of the primate manual dexterity using refined behavioral analyses and further statistical assessments, representing secondary outcome measures from the same manual dexterity test. Thirteen adult monkeys were studied; seven received an anti-Nogo-A antibody whereas a control antibody was infused into the other monkeys. Monkeys were trained to perform the modified Brinkman board task requiring opposition of index finger and thumb to grasp food pellets placed in vertically and horizontally oriented slots. Two parameters were quantified before and following spinal cord injury: (i) the standard 'score' as defined by the number of pellets retrieved within 30 s from the two types of slots; (ii) the newly introduced 'contact time' as defined by the duration of digit contact with the food pellet before successful retrieval. After lesion the hand was severely impaired in all monkeys; this was followed by progressive functional recovery. Remarkably, anti-Nogo-A antibody-treated monkeys recovered faster and significantly better than control antibody-treated monkeys, considering both the score for vertical and horizontal slots (Mann-Whitney test: P = 0.05 and 0.035, respectively) and the contact time (P = 0.008 and 0.005, respectively). Detailed analysis of the lesions excluded the possibility that this conclusion may have been caused by differences in lesion properties between the two groups of monkeys.

What is "clinical data"? Why and how can they be collected during field surveys on medicinal plants?

ETHNOPHARMACOLOGICAL RELEVANCE: "Reverse pharmacology", also called "bedside-to-bench" or "field to pharmacy" approach, is a research process starting with documentation of clinical outcome as observed by patients with different therapeutic regimens. The treatment most significantly associated with cure is selected for future studies: first, clinical safety and efficacy; then in vivo and vitro studies. Some clinical data, i.e. details on patient status and progress, can be collected during ethnobotanical surveys; they will help clinical researchers and, once effectiveness and safety are established, will also help users of traditional medicine make safer and more effective choices. To gather clinical data successfully, ethnopharmacologists need to be backed by an appropriate team of specialists in medicine and epidemiology. Ethnopharmacologists can also gather important data on traditional medicine safety. MATERIALS AND METHODS: The first step is to create a consensus on the meaning of "clinical data", their interest and importance. An understanding of why "a cure is not a proof of effectiveness" is a starting point to avoid faulty interpretation of the clinical observations. RESULTS: Experience showed that, with the "bedside-to-bench" approach, a treatment derived from traditional recipe can be scientifically validated (in terms of safety and effectiveness) with a cost of less than a million euros, thus providing an end-product that is affordable, available and sustainable. CONCLUSIONS: With rigorous clinical study results, medicinal plant users gain the possibility to refine heath strategies. The field surveyor may gain a better relationship with the population, once she/he is seen as bringing information useful for the quality of care in the community.

Estimating the basic reproductive number from viral sequence data.

Epidemiological processes leave a fingerprint in the pattern of genetic structure of virus populations. Here, we provide a new method to infer epidemiological parameters directly from viral sequence data. The method is based on phylogenetic analysis using a birth-death model (BDM) rather than the commonly used coalescent as the model for the epidemiological transmission of the pathogen. Using the BDM has the advantage that transmission and death rates are estimated independently and therefore enables for the first time the estimation of the basic reproductive number of the pathogen using only sequence data, without further assumptions like the average duration of infection. We apply the method to genetic data of the HIV-1 epidemic in Switzerland.

Intraoral open reduction and internal fixation of displaced mandibular angle fractures using a specific ad hoc reduction-compression forceps: a preliminary study

OBJECTIVES: To preliminarily evaluate prospectively the accuracy and reliability of a specific ad hoc reduction-compression forceps in intraoral open reduction of transverse and displaced mandibular angle fractures. STUDY DESIGN: We analyzed the clinical and radiologic data of 7 patients with 7 single transverse and displaced angle fractures. An intraoral approach was performed in all of the patients without using perioperative intermaxillary fixation. A single Arbeitsgemeinschaft Osteosynthese (AO) unilock reconstruction plate was fixed to each stable fragment with 3 locking screws (2.0 mm in 5 patients and 2.4 mm in 2 patients) at the basilar border of the mandible, according to AO/American Society of Internal Fixation (ASIF) principles. Follow-up was at 1, 3, 6, and 12 months, and we noted the status of healing and complications, if any. RESULTS: All of the patients had satisfactory fracture reduction as well as a successful treatment outcome without complications. CONCLUSION: This preliminary study demonstrated that the intraoral reduction of transverse and displaced angle fractures using a specific ad hoc reduction-forceps results in a high rate of success.

Impaired musculoskeletal response to age and exercise in PPARβ(-/-) diabetic mice.

Fragility fractures are recognized complication of diabetes, but yet the underlying mechanisms remain poorly understood. This is particularly pronounced in type 2 diabetes in which the propensity to fall is increased but bone mass is not necessarily low. Thus, whether factors implicated in the development of insulin resistance and diabetes directly impact on the musculoskeletal system remains to be investigated. PPARβ(-/-) mice have reduced metabolic activity and are glucose intolerant. We examined changes in bone and muscle in PPARβ(-/-) mice and investigated both the mechanism behind those changes with age as well as their response to exercise. Compared with their wild type, PPARβ(-/-) mice had an accelerated and parallel decline in both muscle and bone strength with age. These changes were accompanied by increased myostatin expression, low bone formation, and increased resorption. In addition, mesenchymal cells from PPARβ(-/-) had a reduced proliferation capacity and appeared to differentiate into more of an adipogenic phenotype. Concomitantly we observed an increased expression of PPARγ, characteristic of adipocytes. The anabolic responses of muscle and bone to exercise were also diminished in PPARβ(-/-) mice. The periosteal bone formation response to direct bone compression was, however, maintained, indicating that PPARβ controls periosteal bone formation through muscle contraction and/or metabolism. Taken together, these data indicate that PPARβ deficiency leads to glucose intolerance, decreased muscle function, and reduced bone strength. On a molecular level, PPARβ appears to regulate myostatin and PPARγ expression in muscle and bone, thereby providing potential new targets to reverse bone fragility in patients with metabolic disturbances.

The Exotic Palaeo-Thethys Terrane in Central Iran: New geological data from Anarak, Jandaq and Posht-e-Badam areas

Résumé Le « terrane » d'Anarak-Jandak occupe une position géologique clé au nord-ouest du Microcontinent Centre-East Iranien (CE1M), connecté avec le Bloc du Grand Kavir et la ceinture métamorphique de Sanandaj-Sirjan. Nous discutons ici l'origine de ces différentes unités, reliées jusqu'à présent à des épisodes orogéniques d'âge Précambrien à Paléozoïque inférieur, pour conclure finalement de leur affinité paléotéthysienne. Leur histoire commence par un épisode de rifting d'âge Ordovicien supérieur-Dévonien inférieur, pour se terminer au Trias par la collision des blocs Cimmériens dérivé du Gondwana avec le Bloc du Turan d'affinité asiatique (événement Eocimmérien). La plus importante unité métamorphique affleurant au sud-ouest de la région de Jandak-Anarak-Kaboudan est une épaisse séquence silicoclastique à grains fins contenant des blocs ophiolitiques (marginal-sea-type), et des associations basalte-gabbro à signatures géochimiques de type supra-subduction. Dans la région de Nakhlak, nous avons daté ces gabbros par la méthode U-Pb à 387f0.11 Ma ; les roches métamorphiques pélitiques ont donné des âges de refroidissement Ar-Ar pour la muscovite de 320 à 333 Ma. Ce complexe d'accrétion "varisque" a été métamorphisé dans le faciès schiste vert-amphibolite au cours de l'accrétion de la ceinture granitique d'Airekan, d'âge Cambrien inférieur (549±15 Ma par la méthode U/Pb), qui affleure aujourd'hui à l'extrémité nord-ouest du terrane d'Anarak-Jandak . La subduction vers le nord de l'océan Paléotéthys depuis le Paléazoïque supérieur jusqu'au Trias, a permis l'accumulation de grandes quantités de matériel océanique dans la zone de subduction. Par exemple, une succession de guyots (Anarak, Kaboudan, et Meraji Seamounts) et de hauts sous-marins, entrés en collision oblique avec le prisme d'accrétion, est à l'origine d'un léger métamorphisme de type HP qui affecte ces séries {âges Ar-Ar de 280 à 230 Ma). De plus, le magmatisme bimodal de Chah Gorbeh est caractérisé d'une part par des roches de type trondjémite-gabbros (262 Ma), d'autre part par des laves en coussin de type basaltes alcalins-rhyolites; ces roches magmatiques ont recoupé l'ophiolite d'Anarak lors de la mise en place de cette dernière dans la fosse interne de subduction. Quant au prisme d'accrétion de Doshakh, d'âge essentiellement Permien supérieur, i1 a été accrété le long de la marge continentale et métamorphisé dans le faciès schiste vert. La fermeture de la Paléotéthys s'enregistre finalement par la sédimentation dans le bassin d'avant pays du flysch de Bayazeh, d'âge probable Triasique. Le matériel issu de l'arc magmatique de la Paléotéthys est très bien préservé dans les dépôts infra-arc Dévonien supérieur-Carbonifère de Godar-e-Siah, ainsi que dans la succession d'avant-arc de Nakhlak. Pendant l'intervalle Paléozoïque supérieur-Trias, la région de Jandak a été soumise à un régime extensif de type bassin d'arrière-arc, dont un témoin pourrait être la ceinture ophiolitique d'Arusan, elle-même comparable aux écailles ophiolitiques d'Aghdarband au nord-est de l'Iran. Cet ensemble métamorphique est recoupé par des granites d'arc à collisionnel datés à 215±15 Ma. Dans la région de Yazd, témoin de la marge passive Cimmérienne, la sédimentation syn-rift Silurienne à Dévonienne inférieure a été interrompue pendant l'intervalle Trias moyen-Trias supérieur; il en a été de même pour les dépôts de plate-forme Paléozoïque supérieur. L'érosion, qui dans ce dernier cas a atteint le Permien, pourrait être liée au bombement flexural de la marge passive. La collision finale n'a pas induit de déformations trop importantes, et se caractérise par la mise en place de nappes sur la marge passive. Cet événement est scellé par des dépôts molassique du Lias. D'un point de vue régional, la zone s'étendant actuellement de la Mer Noire au Pamir a été soumise à six épisodes d'extension-compression du Jurassique inférieur (début du l'ouverture en position arrière-arc de la Néotéthys) à l'Eocène moyen. Par exemple, le terrane d'AnarakJandak, probablement situé entre le Kopeh Dagh et la plate-forme nord Afghane, s'est complètement détaché de sa patrie d'origine au début du Crétacé supérieur. Des preuves de cet événement se retrouvent dans les séries de plate-forme de Khur (préservation de séries syn-rift puis de marge passive). Les ophiolites de Nain et de Sabzevar sont de plus interprétée comme un témoin de l'existence de ce bassin d'arrière-arc. Dans l'intervalle Eocène-Oligocène, l'indentation par la plaque indienne de l'Eurasie a été contemporaine de la rotation horaire de fragments de l'ancien microcontinent Iranien et de la formation du CEIM. Cette rotation est responsable du transport du terrane d'Anarak-Jandak vers sa position actuelle en Iran Central, et de la dislocation de Terranes de moindre importance, comme le bloc de Posht-e Badam. Depuis le Miocène supérieur, et à la suite de la collision entre l'Arabie et l'Iran, le ternane d'Anarak-Jandak a subi des déformations liées à l'activité d'une zone de cisaillement dextre parallèle à la suture du Zagros, à l'arrière de l'arc magmatique d'Uromieh-Dokhtar. Résumé large public Le Microcontinent Centre-Est Iranien occupe une position géologique clé au centre de l'Iran. Les différentes unités qui le composent, reliées jusqu'à présent à des épisodes orogéniques d'âge Précambrien à Paléozoïque inférieur, sont maintenant rajeunies et liés à la fermeture de l'océean Paléotéthys. Leur histoire commence par un épisode de rifting d'âge Ordovicien supérieur à Dévonien inférieur, pour se terminer au Trias par la collision des- blocs Cimmériens, dérivés du Gondwana, avec le Bloc du Turan d'affinité asiatique. Dans la marge active asiatique de la Paléotéthys, nous avons daté les restes d'un océan marginal à 387±0.11 Ma. Ce complexe d'accrétion a été métamorphisé au cours de la réaccrétion de la ceinture granitique d'Airekan, d'âge Cambrien inférieur (549±15 Ma), qui affleure aujourd'hui à l'extrémité nord-ouest du « terrane » d'Anarak-Jandak correspondant à la plus grande partie de la région étudiée. Le matériel issu de l'arc magmatique de la Paléotéthys est très bien préservé et daté du Dévonien supérieur-Carbonifère. Pendant l'intervalle Paléozoïque supérieur-Trias, la région a été soumise à un régime extensif de type bassin d'arrière-arc, dont un témoin pourrait être la ceinture ophiolitique d'Arusan, comparable aux écailles ophiolitiques d'Aghdarband au nord-est de l'Iran. Cet ensemble métamorphique est recoupé par des granites datés à 215±15 Ma. La subduction vers le nord de l'océan Paléotéthys depuis le Paléozoïque supérieur jusqu'au Trias, a permis l'accumulation de grandes quantités de matériel océanique dans la zone de subduction. Par exemple, une succession de volcans sous-marins, entrés en collision avec le prisme d'accrétion, est à l'origine d'un léger métamorphisme de type HP qui affecte ces séries (280 à 230 Ma). Quant au prisme d'accrétion de Doshakh, d'âge essentiellement Permien supérieur, il a été mis en place le long de la marge continentale et métamorphisé dans le faciès schiste vert. La fermeture de la Paléotéthys s'enregistre finalement par la sédimentation dans le bassin d'avant pays du flysch de Bayazeh, d'âge Triasique. Dans la région de Yazd, on trouve les témoins de la marge passive Cimmérienne, la sédimentation syn-rift Silurienne à Dévonienne inférieure a été interrompue pendant l'intervalle Trias moyen-Trias supérieur, marqué par la flexuration de la marge passive lorsqu'elle rentra en collision avec la marge active asiatique. Cet événement est scellé par des dépôts molassique à charbon du Lias. Le «terrane» d'Anarak-Jandak, probablement situé à l'origine entre le Kopeh Dagh et la plate-forme nord Afghane, s'est complètement détaché de cette région au début du Crétacé supérieur lors de l'ouverture d'un bassin d'arrière-arc, engendré, cette fois, par la subduction de l'océan Néotéthys situé au sud des blocs cimmériens. Des preuves de cet événement se retrouvent dans les séries syn-rift, puis de marge passive de Khour. Les ophiolites de Nain et de Sabzevar sont interprétées comme un témoin de l'existence de ce bassin d'arrière-arc. Dans l'intervalle Eocène-Oligocène, l'indentation de l'Eurasie par la plaque indienne a été contemporaine de la rotation horaire de fragments de l'ancien microcontinent centre-Iranien. Cette rotation de près de 90° est responsable du transport du « terrane » d'Anarak-Jandak vers sa position actuelle. Abstract The Anarak-Jandaq terrane occupies a strategic geological situation at the north-western part of the Central-East Iranian Microcontinent (CEIM) and in connection with the Great Kavir Block and Sanandaj-Sirjan metamorphic belt. Our recent findings redefine the origin of these mentioned areas so far attributed to the Precambrian-Early Palaeozoic orogenic episodes, to be now directly related to the tectonic evolution of the Palaeo-Tethys Ocean, commenced by Late Ordovician-Early Devonian rifting events and terminated in the Triassic by the Eocimmerian tectonic event due to the collision of the Cimmerian blocks with the Asiatic Turan block. The most distributed metamorphic unit that is exposed from the south-west of Jandaq to the Anarak and Kaboudan areas is a thick and fine grain siliciclastic sequence accompanied by marginal-sea-basin ophiolitic blocks including basalt-gabbro association with supra-subduction-geochemical signature. These gabbros in the Nakhlak area were dated by U/Pb method at 387.6 ± 0.11 Ma and the metamorphic pelitic rocks yielded a range of 320 to 333 Ma muscovite-cooling ages based on 40Ar/39 Ar method. This "Variscan" accretionary complex was metamorphosed in greenschist-amphibolite facies during accretion to the Lower Cambrian Airekan granitic belt (549 ± 15 Ma by U/Pb method) that crops out at the northwestern edge of the Anarak-Jandaq terrane. Continued northward subduction of the Palaeo-Tethys Ocean during the entire Late Palaeozoic-Middle Triassic brought huge amount of oceanic material to the subduction zone. One chain of Carboniferous-Triassic oceanic rises and seamounts (the Anarak, Kaboudan, and Meraji Seamounts) obliquely collided with the accretionary wedge and created a mild HP metamorphic event (280-230 Ma based on 40Ar/39Ar results). Bimodal magmatism of the Chah Gorbeh area is characterized by a 262 Ma trondjemite-gabbro as well as pillow alkalibasalts-rhyolites which intruded the Anarak ophiolite when it was being emplaced within the inner-wall trench. The mainly Late Permian-Triassic Doshakh wedge was accreted along the continent and metamorphosed under lower greenschist facies and the probable Triassic Bayazeh flysch filled the foreland basin during the final closure. The Palaeo-Tethys magmatic arc products have been well preserved in the Late Devonian-Carboniferous Godar-e-Siah intra-arc deposits and the Triassic Nakhlak fore-arc succession. During the Late Palaeozoic-Triassic times, the Jandaq area has been affected by back-arc extension and probably the Arusan ophiolitic belt is the remnant of this narrow basin comparable to the Aqdarband ophiolitic remnant in north-east Iran. This metamorphic belt was intruded by 215 ± 15 Ma arc to collisional granites. In the passive margin of the Cimmerian block, on the Yazd region, the Silurian-Early Devonian syn-rift succession as well as the nearly continuous Upper Palaeozoic platform-type deposition was interrupted during the Middle to Late Triassic time, local erosion down to Devonian levels may be related to flexural bulge erosion. The collision event was not so strong to generate intensive deformation but was accompanied by some nappe thrusting onto the passive margin. It is finally unconformably covered by Liassic continental molassic deposits. Related to the onset of Neo-Tethyan back-arc opening in Early Jurassic to Mid-Eocene times, six periods of extensional-compressional events have differently influenced an elongated area, extending from the West Black Sea to Pamir. The Anarak-Jandaq terrane which was situated somewhere in this affected area, probably between the Kopeh Dagh and North Afghan platform, was completely detached from its source at the beginning of the Late Cretaceous

Drug Utilization Studies and Data Registries in Primary Care.

This article reviews the methodology of the studies on drug utilization with particular emphasis on primary care. Population based studies of drug inappropriateness can be done with microdata from Health Electronic Records and e-prescriptions. Multilevel models estimate the influence of factors affecting the appropriateness of drug prescription at different hierarchical levels: patient, doctor, health care organization and regulatory environment. Work by the GIUMAP suggest that patient characteristics are the most important factor in the appropriateness of prescriptions with significant effects at the general practicioner level.

High-resolution 3D whole-heart coronary MRA: a study on the combination of data acquisition in multiple breath-holds and 1D residual respiratory motion compensation.

OBJECT: To study a scan protocol for coronary magnetic resonance angiography based on multiple breath-holds featuring 1D motion compensation and to compare the resulting image quality to a navigator-gated free-breathing acquisition. Image reconstruction was performed using L1 regularized iterative SENSE. MATERIALS AND METHODS: The effects of respiratory motion on the Cartesian sampling scheme were minimized by performing data acquisition in multiple breath-holds. During the scan, repetitive readouts through a k-space center were used to detect and correct the respiratory displacement of the heart by exploiting the self-navigation principle in image reconstruction. In vivo experiments were performed in nine healthy volunteers and the resulting image quality was compared to a navigator-gated reference in terms of vessel length and sharpness. RESULTS: Acquisition in breath-hold is an effective method to reduce the scan time by more than 30 % compared to the navigator-gated reference. Although an equivalent mean image quality with respect to the reference was achieved with the proposed method, the 1D motion compensation did not work equally well in all cases. CONCLUSION: In general, the image quality scaled with the robustness of the motion compensation. Nevertheless, the featured setup provides a positive basis for future extension with more advanced motion compensation methods.

Mutation identification by direct comparison of whole-genome sequencing data from mutant and wild-type individuals using k-mers.

Genes underlying mutant phenotypes can be isolated by combining marker discovery, genetic mapping and resequencing, but a more straightforward strategy for mapping mutations would be the direct comparison of mutant and wild-type genomes. Applying such an approach, however, is hampered by the need for reference sequences and by mutational loads that confound the unambiguous identification of causal mutations. Here we introduce NIKS (needle in the k-stack), a reference-free algorithm based on comparing k-mers in whole-genome sequencing data for precise discovery of homozygous mutations. We applied NIKS to eight mutants induced in nonreference rice cultivars and to two mutants of the nonmodel species Arabis alpina. In both species, comparing pooled F2 individuals selected for mutant phenotypes revealed small sets of mutations including the causal changes. Moreover, comparing M3 seedlings of two allelic mutants unambiguously identified the causal gene. Thus, for any species amenable to mutagenesis, NIKS enables forward genetics without requiring segregating populations, genetic maps and reference sequences.

Rapid fatal outcome from pulmonary arteries compression in transitional cell carcinoma.

Transitional cell carcinoma of the urinary bladder is a malignancy that metastasizes frequently to lymph nodes including the mediastinal lymph nodes. This occurrence may produce symptoms due to compression of adjacent structures such as the superior vena cava syndrome or dysphagia from esophageal compression. We report the case of a 59-year-old man with metastatic transitional cell carcinoma for whom mediastinal lymphadenopathy led to pulmonary artery compression and a rapidly fatal outcome. This rare occurrence has to be distinguished from pulmonary embolism, a much more frequent event in cancer patients, in order that proper and prompt treatment be initiated.

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