Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data.

Despite the central role of quantitative PCR (qPCR) in the quantification of mRNA transcripts, most analyses of qPCR data are still delegated to the software that comes with the qPCR apparatus. This is especially true for the handling of the fluorescence baseline. This article shows that baseline estimation errors are directly reflected in the observed PCR efficiency values and are thus propagated exponentially in the estimated starting concentrations as well as 'fold-difference' results. Because of the unknown origin and kinetics of the baseline fluorescence, the fluorescence values monitored in the initial cycles of the PCR reaction cannot be used to estimate a useful baseline value. An algorithm that estimates the baseline by reconstructing the log-linear phase downward from the early plateau phase of the PCR reaction was developed and shown to lead to very reproducible PCR efficiency values. PCR efficiency values were determined per sample by fitting a regression line to a subset of data points in the log-linear phase. The variability, as well as the bias, in qPCR results was significantly reduced when the mean of these PCR efficiencies per amplicon was used in the calculation of an estimate of the starting concentration per sample.

Revisió morfológica de smd amb del(5q) aïllada i amb del (5q) + 1

Van den Berghe et al. va descriure el “Sdr. 5q-” com un SMD amb del(5q) aïllada, curs indolent, predomini femení, anèmia macrocítica, amb o sense trombocitosi i megacariocits monolobulats. El SMD amb del(5q) aïllada constitueix un subtipus específic de SMD en la classificació de la OMS. L’objectiu del estudi era revisar mostres de SP i MO de SMD amb del(5q) aïllada o amb del(5q) més una alteració citogenètica. Els resultats mostren la dificultat de classificació dels SMD amb les classificacions actuals, i suggereix que no hi ha diferencies entre els SMD amb del(5q) aïllada i els SMD amb del(5q) +1.

Toxicity at three years with and without irradiation of the internal mammary and medial supraclavicular lymph node chain in stage I to III breast cancer (EORTC trial 22922/10925).

INTRODUCTION: The EORTC 22922/10925 trial investigated the potential survival benefit and toxicity of elective irradiation of the internal mammary and medial supraclavicular (IM-MS) nodes Accrual completed in January 2004 and first results are expected in 2012. We present the toxicity reported until year 3 after treatment. PATIENTS AND METHODS: At each visit, toxicity was reported but severity was not graded routinely. Toxicity rates and performance status (PS) changes at three years were compared by chi(2) tests and logistic regression models in all the 3,866 of 4,004 patients eligible to the trial who received the allocated treatment. RESULTS: Only lung (fibrosis; dyspnoea; pneumonitis; any lung toxicities) (4.3% vs. 1.3%; p < 0.0001) but not cardiac toxicity (0.3% vs. 0.4%; p = 0.55) significantly increased with IM-MS treatment. No significant worsening of the PS was observed (p = 0.79), suggesting that treatment-related toxicity does not impair patient's daily activities. CONCLUSIONS: IM-MS irradiation seems well tolerated and does not significantly impair WHO PS at three years. A follow-up period of at least 10 years is needed to determine whether cardiac toxicity is increased after radiotherapy.

Does lamotrigine use in pregnancy increase orofacial cleft risk relative to other malformations?

OBJECTIVE: To investigate whether first trimester exposure to lamotrigine (LTG) monotherapy is specifically associated with an increased risk of orofacial clefts (OCs) relative to other malformations, in response to a signal regarding increased OC risk. METHODS: Population-based case-control study with malformed controls based on EUROCAT congenital anomaly registers. The study population covered 3.9 million births from 19 registries 1995-2005. Registrations included congenital anomaly among livebirths, stillbirths, and terminations of pregnancy following prenatal diagnosis. Cases were 5,511 nonsyndromic OC registrations, of whom 4,571 were isolated, 1,969 were cleft palate (CP), and 1,532 were isolated CP. Controls were 80,052 nonchromosomal, non-OC registrations. We compared first trimester LTG and antiepileptic drug (AED) use vs nonepileptic non-AED use, for mono and polytherapy, adjusting for maternal age. An additional exploratory analysis compared the observed and expected distribution of malformation types associated with LTG use. RESULTS: There were 72 LTG exposed (40 mono- and 32 polytherapy) registrations. The ORs for LTG monotherapy vs no AED use were 0.67 (95% CI 0.10-2.34) for OC relative to other malformations, 0.80 (95% CI 0.11-2.85) for isolated OC, 0.79 (95% CI 0.03-4.35) for CP, and 1.01 (95% CI 0.03-5.57) for isolated CP. ORs for any AED use vs no AED use were 1.43 (95% CI 1.03-1.93) for OC, 1.21 (95% CI 0.82-1.72) for isolated OC, 2.37 (95% CI 1.54-3.43) for CP, and 1.86 (95% CI 1.07-2.94) for isolated CP. The distribution of other nonchromosomal malformation types with LTG exposure was similar to non-AED exposed. CONCLUSION: We find no evidence of a specific increased risk of isolated orofacial clefts relative to other malformations due to lamotrigine (LTG) monotherapy. Our study is not designed to assess whether there is a generalized increased risk of malformations with LTG exposure.

Personality and melanin-based colouration traits in barn owls, tawny owl and kestrels

Individuals need to adapt to their local environment in order to survive. When selection pressures differ in local populations, polymorphism can evolve. Colour polymorphism is one of the most obvious polymorphisms since it is readily observable. Different sources of colouration exist, but melanin-based colouration is one of the most common in birds. The melanocortin system produces this colouration and because the melanocortin system has pleiotropic effects on behavioural and physiological traits, it is a good candidate to be an underlying mechanism to explain the maintenance of colour polymorphism. In this thesis I studied three different raptors which all display melanin-based colouration; barn owls (Tyto alba), tawny owls (Strix aluco) and Eurasian kestrels (Falco tinnunculus). The main question was if there was a relationship between melanin-based colouration and individual behavioural differences. The underlying hypothesis is that colour could be a signal of certain adaptive traits. Our goal was to find evolutionary explanations for the persistence of colour polymorphism. I found that nestling kestrels and barn owls differ in anti-predatory behaviour, with respect to their melanic colouration (chapters 1 and 2). Darker individuals show less reaction to human handling, but in kestrels aggression and colouration are related in opposite ways than in barn owls. More reddish barn owls travel greater distances in natal dispersal and this behaviour is repeatable between parents and same sex offspring (chapter 3). Dark reddish tawny owls defend their nests more intensely against intruders and appear to suffer less from nest predation (chapter 4). Finally I show that polymorphism in the Melanocortin 1 receptor gene (MC1R), which is strongly correlated with reddish colouration in the barn owl, is related to natal dispersal distance, providing a first indication for a genetic basis of the relation between this behaviour and colouration (chapter 5). My results demonstrate a clear link between melanin-based colouration and animal personality traits. I demonstrated this relation in three different species, which shows there is most likely a general underlying mechanism responsible. Different predation pressures might have shaped the reactions to predation, but also differences in sex-related colouration. Male-like and female-like colouration might signal more or less aggressive behaviour. Fluctuating environmental conditions might cause different individual strategies to produce equal reproductive success. The melanocortin system with its pleiotropic effects might be an underlying mechanism, as suggested by the results from the genetic polymorphism, the similar results found in these three species and by the similar relations reported in other species. This thesis demonstrates that colouration and individual differences are correlated and it provides the first glimpse of an underlying system. We can now conduct a more directed search for underlying mechanisms and evolutionary explanations with the use of quantitative genetic methods.

Isolated congenital atrioventricular block diagnosed in utero: natural history and outcome.

BACKGROUND: Isolated congenital atrioventricular block (CAVB) diagnosed in utero is associated with a high morbidity and mortality. Prognosis is especially poor when heart rate drops below 55 beats per minute (bpm) and when fetal hydrops develops. We describe the natural history and outcome of 24 infants with isolated CAVB diagnosed in utero, review the literature, and assess the risk factors that could predict outcome. METHODS: This was a retrospective multicenter study of 24 patients with isolated CAVB diagnosed in utero. RESULTS: CAVB was detected at a mean gestational age (GA) of 24.7 +/- 5.1 weeks. Ten fetuses initially presented with complete heart block. Low heart rate or incomplete heart block was the first documentation of bradyarrhythmia in the other 14 fetuses. In 11 of them, CAVB developed during pregnancy after a median time of 3 (range 1-16) weeks. Fetal hydrops developed in 10 of 24 (42%) fetuses at a mean GA of 27.6 +/- 5.1 weeks. Hydropic fetuses showed lower heart rates during pregnancy (47 +/- 10 bpm) than non-hydropic fetuses (57 +/- 10 bpm). There were three intrauterine deaths; all were hydropic and female. Nine viable females and 12 males were born at a mean GA of 37.1 +/- 6.1 weeks with an average birth weight of 3097 +/- 852 g. Fifteen CAVB patients required pacemaker (PM) intervention, 10 of them immediately after birth. Dilated cardiomyopathy (DCM) developed in three infants of whom two died of congestive heart failure, shortly after the diagnosis was made; one is still alive. Mortality before or after birth was 21%, and was associated with heart rates below 50 bpm and development of fetal hydrops. Poor outcome, defined as death, PM implantation, or development of DCM, occurred in 83% of cases and was associated with heart rates below 60 bpm during pregnancy. CONCLUSIONS: Isolated CAVB diagnosed in utero is associated with high morbidity and mortality. Patients who develop fetal hydrops show lower heart rates during pregnancy than patients who do not. A fetal heart rate below 50 bpm and development of fetal hydrops is associated with increased mortality. Rates below 60 bpm are associated with PM requirement and/or DCM.

Mixing compositions and other scales

Theory of compositional data analysis is often focused on the composition only. However in practical applications we often treat a composition together with covariableswith some other scale. This contribution systematically gathers and develop statistical tools for this situation. For instance, for the graphical display of the dependenceof a composition with a categorical variable, a colored set of ternary diagrams mightbe a good idea for a first look at the data, but it will fast hide important aspects ifthe composition has many parts, or it takes extreme values. On the other hand colored scatterplots of ilr components could not be very instructive for the analyst, if theconventional, black-box ilr is used.Thinking on terms of the Euclidean structure of the simplex, we suggest to set upappropriate projections, which on one side show the compositional geometry and on theother side are still comprehensible by a non-expert analyst, readable for all locations andscales of the data. This is e.g. done by defining special balance displays with carefully-selected axes. Following this idea, we need to systematically ask how to display, explore,describe, and test the relation to complementary or explanatory data of categorical, real,ratio or again compositional scales.This contribution shows that it is sufficient to use some basic concepts and very fewadvanced tools from multivariate statistics (principal covariances, multivariate linearmodels, trellis or parallel plots, etc.) to build appropriate procedures for all these combinations of scales. This has some fundamental implications in their software implementation, and how might they be taught to analysts not already experts in multivariateanalysis

Statistical treatment of grain-size curves and empirical distributions: densities as compositions?

The preceding two editions of CoDaWork included talks on the possible considerationof densities as infinite compositions: Egozcue and D´ıaz-Barrero (2003) extended theEuclidean structure of the simplex to a Hilbert space structure of the set of densitieswithin a bounded interval, and van den Boogaart (2005) generalized this to the setof densities bounded by an arbitrary reference density. From the many variations ofthe Hilbert structures available, we work with three cases. For bounded variables, abasis derived from Legendre polynomials is used. For variables with a lower bound, westandardize them with respect to an exponential distribution and express their densitiesas coordinates in a basis derived from Laguerre polynomials. Finally, for unboundedvariables, a normal distribution is used as reference, and coordinates are obtained withrespect to a Hermite-polynomials-based basis.To get the coordinates, several approaches can be considered. A numerical accuracyproblem occurs if one estimates the coordinates directly by using discretized scalarproducts. Thus we propose to use a weighted linear regression approach, where all k-order polynomials are used as predictand variables and weights are proportional to thereference density. Finally, for the case of 2-order Hermite polinomials (normal reference)and 1-order Laguerre polinomials (exponential), one can also derive the coordinatesfrom their relationships to the classical mean and variance.Apart of these theoretical issues, this contribution focuses on the application of thistheory to two main problems in sedimentary geology: the comparison of several grainsize distributions, and the comparison among different rocks of the empirical distribution of a property measured on a batch of individual grains from the same rock orsediment, like their composition

EORTC study 26041-22041: phase I/II study on concomitant and adjuvant temozolomide (TMZ) and radiotherapy (RT) with PTK787/ZK222584 (PTK/ZK) in newly diagnosed glioblastoma.

BACKGROUND: Glioblastoma is a highly vascularised tumour with a high expression of both vascular endothelial growth factor (VEGF) and VEGFR. PTK787/ZK222584 (PTK/ZK, vatalanib), a multiple VEGF receptor inhibitor, blocks the intracellular tyrosine kinase activity of all known VEGF receptors and is therefore suitable for long-term therapy of pathologic tumour neovascularisation. PATIENTS AND METHODS: The study was designed as an open-label, phase I/II study. A classic 3+3 design was selected. PTK/ZK was added to standard concomitant and adjuvant treatment, beginning in the morning of day 1 of radiotherapy (RT), and given continuously until disease progression or toxicity. PTK/ZK doses started from 500 mg with subsequent escalations to 1000 and 1250 mg/d. Adjuvant or maintenance PTK after the end of radiochemotherapy was given at a previously established dose of 750 mg twice daily continuously with TMZ at the standard adjuvant dose. RESULTS: Twenty patients were enrolled. Dose-limiting toxicities at a once daily dose of 1250 mg were grade 3 diarrhoea (n=1), grade 3 ALT increase (n=2), and myelosuppression with grade 4 thrombocytopenia and neutropenia (n=1). The recommended dose of PTK/ZK in combination with radiotherapy and temozolomide (TMZ) is 1000 mg once a day. This treatment is safe and well tolerated. CONCLUSION: In our phase I study once daily administration of up to 1000 mg of PTK/ZK in conjunction with concomitant temozolomide and radiotherapy was feasible and safe. Prolonged administration of this oral agent is manageable. The planned randomised phase II trial was discontinued right at its onset due to industry decision not to further develop this agent.

A genome-wide significant linkage for severe depression on chromosome 3: the depression network study.

Objective: The purpose of this study was to find loci for major depression via linkage analysis of a large sibling pair sample. Method: The authors conducted a genome-wide linkage analysis of 839 families consisting of 971 affected sibling pairs with severe recurrent major depression, comprising waves I and II of the Depression Network Study cohort. In addition to examining affected status, linkage analyses in the full data set were performed using diagnoses restricted by impairment severity, and association mapping of hits in a large case-control data set was attempted. Results: The authors identified genome-wide significant linkage to chromosome 3p25-26 when the diagnoses were restricted by severity, which was a maximum LOD score of 4.0 centered at the linkage marker D3S1515. The linkage signal identified was genome-wide significant after correction for the multiple phenotypes tested, although subsequent association mapping of the region in a genome-wide association study of a U.K. depression sample did not provide significant results. Conclusions: The authors report a genome-wide significant locus for depression that implicates genes that are highly plausible for involvement in the etiology of recurrent depression. Despite the fact that association mapping in the region was negative, the linkage finding was replicated by another group who found genome-wide-significant linkage for depression in the same region. This suggests that 3p25-26 is a new locus for severe recurrent depression. This represents the first report of a genome-wide significant locus for depression that also has an independent genome-wide significant replication.

Appropriate end-points for right results in the age of antiangiogenic agents: Future options for phase II trials in patients with recurrent glioblastoma.

The progression-free survival rate at 6months (PFS-6) has long been considered the best end-point for assessing the efficacy of new agents in phase II trials in patients with recurrent glioblastoma. However, due to the introduction of antiangiogenic agents in this setting, and their intrinsic propensity to alter neuroradiological disease assessment by producing pseudoregression, any end-point based on neuroradiological modifications should be reconsidered. Further, statistically significant effects on progression-free survival (PFS) only should not automatically be considered reliable evidence of meaningful clinical benefit. In this context, because of its direct and unquestionable clinical relevance, overall survival (OS) represents the gold standard end-point for measuring clinical efficacy, despite the disadvantage that it is influenced by subsequent therapies and usually takes longer time to be evaluated. Therefore, while awaiting novel imaging criteria for response evaluation and/or new imaging tools to distinguish between 'true' and 'pseudo'-responses to antiangiogenic agents, the measurement of OS or OS rates should be considered primary end-points, also in phase II trials with these agents.

Analyzing shapes as compositions of distances

We propose to analyze shapes as “compositions” of distances in Aitchison geometry asan alternate and complementary tool to classical shape analysis, especially when sizeis non-informative.Shapes are typically described by the location of user-chosen landmarks. Howeverthe shape – considered as invariant under scaling, translation, mirroring and rotation– does not uniquely define the location of landmarks. A simple approach is to usedistances of landmarks instead of the locations of landmarks them self. Distances arepositive numbers defined up to joint scaling, a mathematical structure quite similar tocompositions. The shape fixes only ratios of distances. Perturbations correspond torelative changes of the size of subshapes and of aspect ratios. The power transformincreases the expression of the shape by increasing distance ratios. In analogy to thesubcompositional consistency, results should not depend too much on the choice ofdistances, because different subsets of the pairwise distances of landmarks uniquelydefine the shape.Various compositional analysis tools can be applied to sets of distances directly or afterminor modifications concerning the singularity of the covariance matrix and yield resultswith direct interpretations in terms of shape changes. The remaining problem isthat not all sets of distances correspond to a valid shape. Nevertheless interpolated orpredicted shapes can be backtransformated by multidimensional scaling (when all pairwisedistances are used) or free geodetic adjustment (when sufficiently many distancesare used)

A compositional data analysis package for R providing multiple approaches

”compositions” is a new R-package for the analysis of compositional and positive data.It contains four classes corresponding to the four different types of compositional andpositive geometry (including the Aitchison geometry). It provides means for computation,plotting and high-level multivariate statistical analysis in all four geometries.These geometries are treated in an fully analogous way, based on the principle of workingin coordinates, and the object-oriented programming paradigm of R. In this way,called functions automatically select the most appropriate type of analysis as a functionof the geometry. The graphical capabilities include ternary diagrams and tetrahedrons,various compositional plots (boxplots, barplots, piecharts) and extensive graphical toolsfor principal components. Afterwards, ortion and proportion lines, straight lines andellipses in all geometries can be added to plots. The package is accompanied by ahands-on-introduction, documentation for every function, demos of the graphical capabilitiesand plenty of usage examples. It allows direct and parallel computation inall four vector spaces and provides the beginner with a copy-and-paste style of dataanalysis, while letting advanced users keep the functionality and customizability theydemand of R, as well as all necessary tools to add own analysis routines. A completeexample is included in the appendix

News from compositions, the R package

The R-package “compositions”is a tool for advanced compositional analysis. Its basicfunctionality has seen some conceptual improvement, containing now some facilitiesto work with and represent ilr bases built from balances, and an elaborated subsys-tem for dealing with several kinds of irregular data: (rounded or structural) zeroes,incomplete observations and outliers. The general approach to these irregularities isbased on subcompositions: for an irregular datum, one can distinguish a “regular” sub-composition (where all parts are actually observed and the datum behaves typically)and a “problematic” subcomposition (with those unobserved, zero or rounded parts, orelse where the datum shows an erratic or atypical behaviour). Systematic classificationschemes are proposed for both outliers and missing values (including zeros) focusing onthe nature of irregularities in the datum subcomposition(s).To compute statistics with values missing at random and structural zeros, a projectionapproach is implemented: a given datum contributes to the estimation of the desiredparameters only on the subcompositon where it was observed. For data sets withvalues below the detection limit, two different approaches are provided: the well-knownimputation technique, and also the projection approach.To compute statistics in the presence of outliers, robust statistics are adapted to thecharacteristics of compositional data, based on the minimum covariance determinantapproach. The outlier classification is based on four different models of outlier occur-rence and Monte-Carlo-based tests for their characterization. Furthermore the packageprovides special plots helping to understand the nature of outliers in the dataset.Keywords: coda-dendrogram, lost values, MAR, missing data, MCD estimator,robustness, rounded zeros