Probability of achieving optimal molecular response to imatinib in chronic myeloid leukemia (CML) patients: Pharmacokinetic/Pharmacodynamic (PK/PD) relationships observed under field-conditions

Objectives: Imatinib has been increasingly proposed for therapeutic drug monitoring (TDM), as trough concentrations (Cmin) correlate with response rates in CML patients. This analysis aimed to evaluate the impact of imatinib exposure on optimal molecular response rates in a large European cohort of patients followed by centralized TDM.¦Methods: Sequential PK/PD analysis was performed in NONMEM 7 on 2230 plasma (PK) samples obtained along with molecular response (PD) data from 1299 CML patients. Model-based individual Bayesian estimates of exposure, parameterized as to initial dose adjusted and log-normalized Cmin (log-Cmin) or clearance (CL), were investigated as potential predictors of optimal molecular response, while accounting for time under treatment (stratified at 3 years), gender, CML phase, age, potentially interacting comedication, and TDM frequency. PK/PD analysis used mixed-effect logistic regression (iterative two-stage method) to account for intra-patient correlation.¦Results: In univariate analyses, CL, log-Cmin, time under treatment, TDM frequency, gender (all p<0.01) and CML phase (p=0.02) were significant predictors of the outcome. In multivariate analyses, all but log-Cmin remained significant (p<0.05). Our model estimates a 54.1% probability of optimal molecular response in a female patient with a median CL of 14.4 L/h, increasing by 4.7% with a 35% decrease in CL (percentile 10 of CL distribution), and decreasing by 6% with a 45% increased CL (percentile 90), respectively. Male patients were less likely than female to be in optimal response (odds ratio: 0.62, p<0.001), with an estimated probability of 42.3%.¦Conclusions: Beyond CML phase and time on treatment, expectedly correlated to the outcome, an effect of initial imatinib exposure on the probability of achieving optimal molecular response was confirmed in field-conditions by this multivariate analysis. Interestingly, male patients had a higher risk of suboptimal response, which might not exclusively derive from their 18.5% higher CL, but also from reported lower adherence to the treatment. A prospective longitudinal study would be desirable to confirm the clinical importance of identified covariates and to exclude biases possibly affecting this observational survey.

Quantifying degradation rates of transmembrane receptor kinases.

Transmembrane receptor-kinases are widespread throughout eukaryotes and their activities are known to regulate all kinds of cellular responses in diverse organs and cell types. In order to guarantee the correct amplitude and duration of signals, receptor levels at the cellular surface need to be tightly controlled. The regulation of receptor degradation is the most direct way to achieve this and elaborate mechanisms are in place to control this process. Therefore, the rate of receptor degradation is a parameter of central importance for understanding the dynamics of a signal transduction cascade. Unfortunately, degradation of transmembrane receptors is a complicated multistep process that involves internalization from the plasma membrane, invagination into the lumen of endosomal compartments, and finally fusion with the vacuole for degradation by vacuolar proteases. Therefore, degradation should be measured in an as noninvasive way as possible, such as not to interfere with the complicated transport processes. Here, a method for minimally invasive, in vivo turn-over measurements in intact organs is provided. This technique was used for quantifying the turn-over rates of the Brassinosteroid receptor kinase BRI1 (BRASSINOSTEROID INSENSITIVE 1) in Arabidopsis thaliana root meristems. Pulse-chase expression of a fluorescently labeled BRI1 variant was used and its turn-over rate was determined by quantitative confocal microscopy. This method is well suited to measure turn-over of transmembrane kinases, but can evidently be extended to measure turn-over of any types of transmembrane proteins.

Molecular epidemiology of methicillin resistant Staphylococcus aureus isolated from newborns in a hospital in Rio de Janeiro, Brazil

Methicillin resistant Staphylococcus aureus (MRSA) is an organism that is frequently transmitted in hospitals and perinatal units. The MRSA is considered a public health problem in neonatology because of its strong potential for dissemination in the wards associated with high rates of morbidity and mortality. In this study we describe the bacteriological, epidemiological and molecular characteristics of MRSA isolated from anterior nares and blood cultures of newborns hospitalized in a public maternity hospital in the city of Rio de Janeiro, Brazil. The frequency of MRSA isolated from nasal swabs of newborns was 47.8% (43/90). The genetic analysis of MRSA strains from anterior nares, showed 8 different pulsed field gel electrophoresis patterns (PFGE). Upon analysis of PFGE patterns of the 12 MRSA strains isolated from blood cultures, 8 different patterns were observed, 9 (75%) strains were genetic related to nasal secretion isolates patterns. In conclusion, our data demonstrate the importance of screening of newborns for the presence of MRSA in Brazilian hospitals and the usefulness of genetic typing of these pathogen during epidemiologic studies. This should lead to a better knowledge on the significancy and spreading of MRSA in the hospitals.

Hepatitis B and C in the hemodialysis unit of Tocantins, Brazil: serological and molecular profiles

A survey was conducted in the hemodialysis population of the state of Tocantins, Brazil, aiming to assess the prevalence of hepatitis B virus (HBV) and hepatitis C virus (HCV) infections, to analyze associated risk factors, and also to investigate these viruses genotypes distribution. During January and March 2001, all patients (n = 100) were interviewed at the unique dialysis unit in Tocantins. Blood samples were collected and serum samples were screened for HBV serological markers. Hepatitis B surface antigen positive samples were tested for HBV DNA. All samples were also tested for anti-HCV antibodies and HCV RNA. An overall prevalence of 45% was found for HBV infection (4% were HBsAg/anti-HBc positive, 2% were anti-HBc only and 39% had anti-HBc/anti-HBs markers). Concerning HCV infection, anti-HCV and HCV RNA were detected in 13% and 14% of the subjects, respectively. Three patients were HCV RNA positive and anti-HCV negative, resulting in an overall HCV prevalence of 16%. Univariate analysis of risk factors showed that only shift and length of time on hemodialysis were associated with HBV and HCV positivity, respectively. Among the four HBsAg-positive samples, HBV DNA was detected in three of them, which were identified as genotype A by restriction fragment length polymorphism (RFLP) analysis. All 14 HCV RNA-positive samples were genotyped by INNO-LiPA. Genotypes 1a and 3a were found in 85% and 15%, respectively. The present data show low HBsAg and HCV prevalence rates. The risk factors associated with HBV and HCV positivity suggest that nosocomial transmission may influence in spreading these viruses in the dialysis unit studied.

Substitution rates at neutral genes depend on population size under fluctuating demography and overlapping generations.

It is widely accepted that the rate of evolution (substitution rate) at neutral genes is unaffected by population size fluctuations. This result has implications for the analysis of genetic data in population genetics and phylogenetics, and provides, in particular, a justification for the concept of the molecular clock. Here, we show that the substitution rate at neutral genes does depend on population size fluctuations in the presence of overlapping generations. As both population size fluctuations and overlapping generations are expected to be the norm rather than the exception in natural populations, this observation may be relevant for understanding variation in substitution rates within and between lineages.

Testing advances in molecular discrimination among Chinook salmon life histories: evidence from a blind test.

The application of DNA-based markers toward the task of discriminating among alternate salmon runs has evolved in accordance with ongoing genomic developments and increasingly has enabled resolution of which genetic markers associate with important life-history differences. Accurate and efficient identification of the most likely origin for salmon encountered during ocean fisheries, or at salvage from fresh water diversion and monitoring facilities, has far-reaching consequences for improving measures for management, restoration and conservation. Near-real-time provision of high-resolution identity information enables prompt response to changes in encounter rates. We thus continue to develop new tools to provide the greatest statistical power for run identification. As a proof of concept for genetic identification improvements, we conducted simulation and blind tests for 623 known-origin Chinook salmon (Oncorhynchus tshawytscha) to compare and contrast the accuracy of different population sampling baselines and microsatellite loci panels. This test included 35 microsatellite loci (1266 alleles), some known to be associated with specific coding regions of functional significance, such as the circadian rhythm cryptochrome genes, and others not known to be associated with any functional importance. The identification of fall run with unprecedented accuracy was demonstrated. Overall, the top performing panel and baseline (HMSC21) were predicted to have a success rate of 98%, but the blind-test success rate was 84%. Findings for bias or non-bias are discussed to target primary areas for further research and resolution.

A Novel Pulse-Chase SILAC Strategy Measures Changes in Protein Decay and Synthesis Rates Induced by Perturbation of Proteostasis with an Hsp90 Inhibitor.

Standard proteomics methods allow the relative quantitation of levels of thousands of proteins in two or more samples. While such methods are invaluable for defining the variations in protein concentrations which follow the perturbation of a biological system, they do not offer information on the mechanisms underlying such changes. Expanding on previous work [1], we developed a pulse-chase (pc) variant of SILAC (stable isotope labeling by amino acids in cell culture). pcSILAC can quantitate in one experiment and for two conditions the relative levels of proteins newly synthesized in a given time as well as the relative levels of remaining preexisting proteins. We validated the method studying the drug-mediated inhibition of the Hsp90 molecular chaperone, which is known to lead to increased synthesis of stress response proteins as well as the increased decay of Hsp90 "clients". We showed that pcSILAC can give information on changes in global cellular proteostasis induced by treatment with the inhibitor, which are normally not captured by standard relative quantitation techniques. Furthermore, we have developed a mathematical model and computational framework that uses pcSILAC data to determine degradation constants kd and synthesis rates Vs for proteins in both control and drug-treated cells. The results show that Hsp90 inhibition induced a generalized slowdown of protein synthesis and an increase in protein decay. Treatment with the inhibitor also resulted in widespread protein-specific changes in relative synthesis rates, together with variations in protein decay rates. The latter were more restricted to individual proteins or protein families than the variations in synthesis. Our results establish pcSILAC as a viable workflow for the mechanistic dissection of changes in the proteome which follow perturbations. Data are available via ProteomeXchange with identifier PXD000538.

Survival, prognostic factors and rates of leukemic transformation in 381 untreated patients with MDS and del(5q): a multicenter study.

Myelodysplastic syndromes (MDS) with del(5q) are considered to have a benign course of the disease. In order to address the issue of the propensity of those patients to progress to acute myeloid leukemia (AML), data on 381 untreated patients with MDS and del(5q) characterized by low or intermediate I International Prognostic Scoring System (IPSS) risk score were collected from nine centers and registries. Median survival of the entire group was 74 months. Transfusion-dependent patients had a median survival of 44 months vs 97 months for transfusion-independent patients (P<0.0001). Transfusion need at diagnosis was the most important patient characteristic for survival. Of the 381 patients, 48 (12.6%) progressed to AML. The cumulative progression rate calculated using the Kaplan-Meier method was 4.9% at 2 years and 17.6% at 5 years. Factors associated with the risk of AML transformation were high-risk World Health Organization adapted Prognostic Scoring System (WPSS) score, marrow blast count >5% and red-cell transfusion dependency at diagnosis. In conclusion, patients with MDS and del(5q) are facing a considerable risk of AML transformation. More detailed cytogenetic and molecular studies may help to identify the patients at risk of progression.

Molecular epidemiology of HIV-1 clades in Southern Brazil

Human immunodeficiency virus (HIV) clades B and C account for more than 60% of the HIV-1 infections worldwide. In this paper, we describe the profiles of patients infected with subtypes of HIV-1 from the state of Paraná, Southern Brazil, and correlate them with demographic and epidemiological findings. A retrospective analysis of HIV cases reported from 1999-2007 was also performed. Data from 293 patients were reviewed and 245 were older than 13 (58% female). The distribution of clades was as follows: B 140 (57%), C 67 (23%), F 24 (10%) and mosaic or unique recombinant forms (URFs) 24 (10%). Of the 48 patients younger than 13 years of age (62.5% male), vertical transmission occurred in 46 and the distribution of clades was as follows: B 14 (29%), C 24 (50%), F 7 (15%) and URFs 6 (13%). There was no significant difference in mortality between HIV-1 subtypes. In both groups, patients infected with clade C tended to have higher rates of injection drug use exposure risk.

Molecular characterisation of hepatitis B virus in the resident Chinese population in Panama City

Despite the effectiveness of current hepatitis B virus (HBV) vaccines, it is estimated that 350 million individuals suffer from chronic HBV infection and more than 50% of these affected individuals live on the Asian continent. Panama is a country with a great diversity of foreign groups; the Chinese community is a large example of this phenomenon. There is an urgent need to perform studies that evaluate the prevalence and the genetic diversity of HBV in this community. This study aimed to evaluate the prevalence of HBV and its genotypes and mutant variants in the Chinese population residing in Panama. In total, 320 subjects were enrolled in the study. Forty-two subjects (13.1%) were positive for HBsAg and HBV-DNA from 18 subjects revealed the presence of genotypes B2 and C1. Secondary mutations associated with drug resistance at positions rtV207L and rtN239T of the reverse transcriptase gene were identified. Additionally, the mutation pair A1762T/G1764A was found in three samples and the mutation G1896A was detected in an HBeAg-negative subject. In conclusion, to our knowledge, this is the first study to report high HBV prevalence rates in resident ethnic Chinese in Central America and the presence of genotypes B2 and C1 in this region.

Evaluation of parasitological and molecular techniques for the diagnosis and assessment of cure of schistosomiasis mansoni in a low transmission area

This study evaluated parasitological and molecular techniques for the diagnosis and assessment of cure of schistosomiasis mansoni. A population-based study was performed in 201 inhabitants from a low transmission locality named Pedra Preta, municipality of Montes Claros, state of Minas Gerais, Brazil. Four stool samples were analysed using two techniques, the Kato-Katz® (KK) technique (18 slides) and the TF-Test®, to establish the infection rate. The positivity rate of 18 KK slides of four stool samples was 28.9% (58/201) and the combined parasitological techniques (KK+TF-Test®) produced a 35.8% positivity rate (72/201). Furthermore, a polymerase chain reaction (PCR)-ELISA assay produced a positivity rate of 23.4% (47/201) using the first sample. All 72 patients with positive parasitological exams were treated with a single dose of Praziquantel® and these patients were followed-up 30, 90 and 180 days after treatment to establish the cure rate. Cure rates obtained by the analysis of 12 KK slides were 100%, 100% and 98.4% at 30, 90 and 180 days after treatment, respectively. PCR-ELISA revealed cure rates of 98.5%, 95.5% and 96.5%, respectively. The diagnostic and assessment of cure for schistosomiasis may require an increased number of KK slides or a test with higher sensitivity, such as PCR-ELISA, in situations of very low parasite load, such as after therapeutic interventions.

Functional and molecular evolutionary analysis of the "Brevix Radix" gene family in mono-and dicotyledonous plants

ABSTRACTIn contrast to animals, plants cannot move from their place of birth and, therefore, need to adapt to their particular habitat in order to survive. Thus, plant development is remarkably plastic, making plants an ideal system for the isolation of genes that account for intraspecific natural variation and possibly environmental adaptation. However, to date, this approach mostly identified null alleles and missed mutations with subtle effects. For instance, BREVIS RADIX (BRX) has been isolated as a key regulator of root growth through a naturally occurring loss-of-function allele in the Arabidopsis thaliana accession Uk-1 and is the founding member of a highly-conserved plant-specific gene family.In this work, we show that a strong selective pressure is acting on the BRX gene family and dates back before the monocot-dicot divergence. However, functional diversification is observed mainly in dicotyledon BRX family genes and is correlated with acceleration in the evolutionary rates in the N-terminal regions. Population genetic data revealed that BRX is highly conserved across Arabidopsis accessions and presents signatures of adaptation. Interestingly, a seven amino acid deletion polymorphism in BRX sequence was found in a few accessions, which seems to be responsible for their enhanced primary root growth. Nevertheless, BRX might not only be active in the root, as suggested by its expression in the shoot. Indeed, leaves and cotyledons of brx mutants are significantly smaller than wild- type. This phenotype is a direct consequence of the absence of BRX function in the shoot rather than an indirect effect of an altered root system growth. Interestingly, cotyledons of brx plants reflect the same physiological defects as the root. Moreover, phenotypes in BRX gain-of-function plants, such as epinastic leaves and increased epidermal cell size, could be associated with an increase in leaf brassinosteroid content.Collectively, these results indicate that BRX contributes to local adaptation by ubiquitously regulating plant growth, probably through the modulation of brassinosteroid biosynthesis.RÉSUMÉContrairement à la plupart des animaux, les plantes ne peuvent se mouvoir et doivent ainsi s'adapter à leur environnement pour survivre. Pour cette raison, elles représentent un système idéal pour l'identification de gènes contribuant à la variation naturelle intra- spécifique, ainsi qu'à l'adaptation. Cependant, cette approche a, jusqu'à présent, surtout permis d'isoler des allèles nuls et non des mutations conférant des effets plus subtiles. C'est le cas du gène Β REVIS RADIX (BRX), un régulateur clé de la croissance racinaire, qui a été identifié grâce à un allèle non-fonctionnel présent dans l'accession naturelle d'Arabidopsis thaliana Uk-1. BRX et ses homologues des plantes mono- et dicotylédones forment une famille très conservée et spécifique aux plantes.Dans ce travail, nous démontrons que la famille de gènes BRX est soumise à une forte pression de sélection qui remonte avant la divergence entre mono- et dicotylédones. Cependant, une diversification fonctionnelle a été observée chez les gènes des dicotylédones et corrèle avec une accélération de la vitesse d'évolution dans leur région N- terminale. Une analyse génétique de différentes accessions naturelles d'Arabidopsis a révélé que BRX est hautement conservé et présente des signatures d'adaptation. Remarquablement, un polymorphisme de délétion de sept acides aminés a été détecté dans quelques accessions et a pour conséquence une plus forte croissance de la racine primaire. Néanmoins, il semble que le rôle de BRX ne se limite pas qu'à la racine, comme indiqué par son expression dans les parties aériennes de la plante. En effet, les mutants brx présentent des cotylédons et des feuilles significativement plus petits que le type sauvage, une conséquence directe de l'absence d'activité de BRX dans ces organes. Nous avons aussi noté que les cotylédons des mutants brx, à l'instar des racines, ont une perception altérée de l'auxine et peuvent être complémentés par l'application exogène de brassinostéroïdes. De plus, dans des plantes présentant un gain de fonction BRX, les feuilles sont épinastiques et les cellules de leur épiderme plus grandes. Ces phénotypes sont accompagnés d'une augmentation de la concentration de brassinostéroïdes dans les feuilles. Conjointement, ces résultats démontrent que BRX contribue à une adaptation locale de la plante par la régulation générale de sa croissance, probablement en modulant la biosynthèse des brassinostéroïdes.

Can we continue to neglect genomic variation in introgression rates when inferring the history of speciation? A case study in a Mytilus hybrid zone.

The use of molecular data to reconstruct the history of divergence and gene flow between populations of closely related taxa represents a challenging problem. It has been proposed that the long-standing debate about the geography of speciation can be resolved by comparing the likelihoods of a model of isolation with migration and a model of secondary contact. However, data are commonly only fit to a model of isolation with migration and rarely tested against the secondary contact alternative. Furthermore, most demographic inference methods have neglected variation in introgression rates and assume that the gene flow parameter (Nm) is similar among loci. Here, we show that neglecting this source of variation can give misleading results. We analysed DNA sequences sampled from populations of the marine mussels, Mytilus edulis and M. galloprovincialis, across a well-studied mosaic hybrid zone in Europe and evaluated various scenarios of speciation, with or without variation in introgression rates, using an Approximate Bayesian Computation (ABC) approach. Models with heterogeneous gene flow across loci always outperformed models assuming equal migration rates irrespective of the history of gene flow being considered. By incorporating this heterogeneity, the best-supported scenario was a long period of allopatric isolation during the first three-quarters of the time since divergence followed by secondary contact and introgression during the last quarter. By contrast, constraining migration to be homogeneous failed to discriminate among any of the different models of gene flow tested. Our simulations thus provide statistical support for the secondary contact scenario in the European Mytilus hybrid zone that the standard coalescent approach failed to confirm. Our results demonstrate that genomic variation in introgression rates can have profound impacts on the biological conclusions drawn from inference methods and needs to be incorporated in future studies.

Molecular changes underlying mammalian phenotypic innovation

Mammals are characterized by specific phenotypic traits that include lactation, hair, and relatively large brains with unique structures. Individual mammalian lineages have, in turn, evolved characteristic traits that distinguish them from others. These include obvious anatom¬ical differences but also differences related to reproduction, life span, cognitive abilities, be¬havior. and disease susceptibility. However, the molecular basis of the diverse mammalian phenotypes and the selective pressures that shaped their evolution remain largely unknown. In the first part of my thesis, I analyzed the genetic factors associated with the origin of a unique mammalian phenotype lactation and I studied the selective pressures that forged the transition from oviparity to viviparity. Using a comparative genomics approach and evolutionary simulations, I showed that the emergence of lactation, as well as the appear¬ance of the casein gene family, significantly reduced selective pressure on the major egg-yolk proteins (the vitellogenin family). This led to a progressive loss of vitellogenins, which - in oviparous species - act as storage proteins for lipids, amino acids, phosphorous and calcium in the isolated egg. The passage to internal fertilization and placentation in therian mam¬mals rendered vitellogenins completely dispensable, which ended in the loss of the whole gene family in this lineage. As illustrated by the vitellogenin study, changes in gene content are one possible underlying factor for the evolution of mammalian-specific phenotypes. However, more subtle genomic changes, such as mutations in protein-coding sequences, can also greatly affect the phenotypes. In particular, it was proposed that changes at the level of gene reg¬ulation could underlie many (or even most) phenotypic differences between species. In the second part of my thesis, I participated in a major comparative study of mammalian tissue transcriptomes, with the goal of understanding how evolutionary forces affected expression patterns in the past 200 million years of mammalian evolution. I showed that, while com¬parisons of gene expressions are in agreement with the known species phylogeny, the rate of expression evolution varies greatly among lineages. Species with low effective population size, such as monotremes and hominoids, showed significantly accelerated rates of gene expression evolution. The most likely explanation for the high rate of gene expression evolution in these lineages is the accumulation of mildly deleterious mutations in regulatory regions, due to the low efficiency of purifying selection. Thus, our observations are in agreement with the nearly neutral theory of molecular evolution. I also describe substantial differences in evolutionary rates between tissues, with brain being the most constrained (especially in primates) and testis significantly accelerated. The rate of gene expression evolution also varies significantly between chromosomes. In particular, I observed an acceleration of gene expression changes on the X chromosome, probably as a result of adaptive processes associated with the origin of therian sex chromosomes. Lastly, I identified several individual genes as well as co-regulated expression modules that have undergone lineage specific expression changes and likely under¬lie various phenotypic innovations in mammals. The methods developed during my thesis, as well as the comprehensive gene content analyses and transcriptomics datasets made available by our group, will likely prove to be useful for further exploratory analyses of the diverse mammalian phenotypes.

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