Batch effect confounding leads to strong bias in performance estimates obtained by cross-validation.

BACKGROUND: With the large amount of biological data that is currently publicly available, many investigators combine multiple data sets to increase the sample size and potentially also the power of their analyses. However, technical differences ("batch effects") as well as differences in sample composition between the data sets may significantly affect the ability to draw generalizable conclusions from such studies. FOCUS: The current study focuses on the construction of classifiers, and the use of cross-validation to estimate their performance. In particular, we investigate the impact of batch effects and differences in sample composition between batches on the accuracy of the classification performance estimate obtained via cross-validation. The focus on estimation bias is a main difference compared to previous studies, which have mostly focused on the predictive performance and how it relates to the presence of batch effects. DATA: We work on simulated data sets. To have realistic intensity distributions, we use real gene expression data as the basis for our simulation. Random samples from this expression matrix are selected and assigned to group 1 (e.g., 'control') or group 2 (e.g., 'treated'). We introduce batch effects and select some features to be differentially expressed between the two groups. We consider several scenarios for our study, most importantly different levels of confounding between groups and batch effects. METHODS: We focus on well-known classifiers: logistic regression, Support Vector Machines (SVM), k-nearest neighbors (kNN) and Random Forests (RF). Feature selection is performed with the Wilcoxon test or the lasso. Parameter tuning and feature selection, as well as the estimation of the prediction performance of each classifier, is performed within a nested cross-validation scheme. The estimated classification performance is then compared to what is obtained when applying the classifier to independent data.

Strong time dependence of the 76-gene prognostic signature for node-negative breast cancer patients in the TRANSBIG multicenter independent validation series.

PURPOSE: Recently, a 76-gene prognostic signature able to predict distant metastases in lymph node-negative (N(-)) breast cancer patients was reported. The aims of this study conducted by TRANSBIG were to independently validate these results and to compare the outcome with clinical risk assessment. EXPERIMENTAL DESIGN: Gene expression profiling of frozen samples from 198 N(-) systemically untreated patients was done at the Bordet Institute, blinded to clinical data and independent of Veridex. Genomic risk was defined by Veridex, blinded to clinical data. Survival analyses, done by an independent statistician, were done with the genomic risk and adjusted for the clinical risk, defined by Adjuvant! Online. RESULTS: The actual 5- and 10-year time to distant metastasis were 98% (88-100%) and 94% (83-98%), respectively, for the good profile group and 76% (68-82%) and 73% (65-79%), respectively, for the poor profile group. The actual 5- and 10-year overall survival were 98% (88-100%) and 87% (73-94%), respectively, for the good profile group and 84% (77-89%) and 72% (63-78%), respectively, for the poor profile group. We observed a strong time dependence of this signature, leading to an adjusted hazard ratio of 13.58 (1.85-99.63) and 8.20 (1.10-60.90) at 5 years and 5.11 (1.57-16.67) and 2.55 (1.07-6.10) at 10 years for time to distant metastasis and overall survival, respectively. CONCLUSION: This independent validation confirmed the performance of the 76-gene signature and adds to the growing evidence that gene expression signatures are of clinical relevance, especially for identifying patients at high risk of early distant metastases.

Integrated Analysis of High-Resolution Gene Expression and Copy Number Profiling Identified Biallelic Deletion of CDKN2A/2B Tumor Suppressor Locus As the Most Frequent and Unique Genomic Abnormality in Diffuse Large B-Cell Lymphoma (DLBCL) with Strong Prognostic Value in Both GCB and ABC Subtypes and Not Overcome by a Dose-Intensive Immunochemotherapy Regimen Plus Rituximab. Results of a Prospective GELA Clinical Trial Program

Background and aim of the study: Genomic gains and losses play a crucial role in the development and progression of DLBCL and are closely related to gene expression profiles (GEP), including the germinal center B-cell like (GCB) and activated B-cell like (ABC) cell of origin (COO) molecular signatures. To identify new oncogenes or tumor suppressor genes (TSG) involved in DLBCL pathogenesis and to determine their prognostic values, an integrated analysis of high-resolution gene expression and copy number profiling was performed. Patients and methods: Two hundred and eight adult patients with de novo CD20+ DLBCL enrolled in the prospective multicentric randomized LNH-03 GELA trials (LNH03-1B, -2B, -3B, 39B, -5B, -6B, -7B) with available frozen tumour samples, centralized reviewing and adequate DNA/RNA quality were selected. 116 patients were treated by Rituximab(R)-CHOP/R-miniCHOP and 92 patients were treated by the high dose (R)-ACVBP regimen dedicated to patients younger than 60 years (y) in frontline. Tumour samples were simultaneously analysed by high resolution comparative genomic hybridization (CGH, Agilent, 144K) and gene expression arrays (Affymetrix, U133+2). Minimal common regions (MCR), as defined by segments that affect the same chromosomal region in different cases, were delineated. Gene expression and MCR data sets were merged using Gene expression and dosage integrator algorithm (GEDI, Lenz et al. PNAS 2008) to identify new potential driver genes. Results: A total of 1363 recurrent (defined by a penetrance > 5%) MCRs within the DLBCL data set, ranging in size from 386 bp, affecting a single gene, to more than 24 Mb were identified by CGH. Of these MCRs, 756 (55%) showed a significant association with gene expression: 396 (59%) gains, 354 (52%) single-copy deletions, and 6 (67%) homozygous deletions. By this integrated approach, in addition to previously reported genes (CDKN2A/2B, PTEN, DLEU2, TNFAIP3, B2M, CD58, TNFRSF14, FOXP1, REL...), several genes targeted by gene copy abnormalities with a dosage effect and potential physiopathological impact were identified, including genes with TSG activity involved in cell cycle (HACE1, CDKN2C) immune response (CD68, CD177, CD70, TNFSF9, IRAK2), DNA integrity (XRCC2, BRCA1, NCOR1, NF1, FHIT) or oncogenic functions (CD79b, PTPRT, MALT1, AUTS2, MCL1, PTTG1...) with distinct distribution according to COO signature. The CDKN2A/2B tumor suppressor locus (9p21) was deleted homozygously in 27% of cases and hemizygously in 9% of cases. Biallelic loss was observed in 49% of ABC DLBCL and in 10% of GCB DLBCL. This deletion was strongly correlated to age and associated to a limited number of additional genetic abnormalities including trisomy 3, 18 and short gains/losses of Chr. 1, 2, 19 regions (FDR < 0.01), allowing to identify genes that may have synergistic effects with CDKN2A/2B inactivation. With a median follow-up of 42.9 months, only CDKN2A/2B biallelic deletion strongly correlates (FDR p.value < 0.01) to a poor outcome in the entire cohort (4y PFS = 44% [32-61] respectively vs. 74% [66-82] for patients in germline configuration; 4y OS = 53% [39-72] vs 83% [76-90]). In a Cox proportional hazard prediction of the PFS, CDKN2A/2B deletion remains predictive (HR = 1.9 [1.1-3.2], p = 0.02) when combined with IPI (HR = 2.4 [1.4-4.1], p = 0.001) and GCB status (HR = 1.3 [0.8-2.3], p = 0.31). This difference remains predictive in the subgroup of patients treated by R-CHOP (4y PFS = 43% [29-63] vs. 66% [55-78], p=0.02), in patients treated by R-ACVBP (4y PFS = 49% [28-84] vs. 83% [74-92], p=0.003), and in GCB (4y PFS = 50% [27-93] vs. 81% [73-90], p=0.02), or ABC/unclassified (5y PFS = 42% [28-61] vs. 67% [55-82] p = 0.009) molecular subtypes (Figure 1). Conclusion: We report for the first time an integrated genetic analysis of a large cohort of DLBCL patients included in a prospective multicentric clinical trial program allowing identifying new potential driver genes with pathogenic impact. However CDKN2A/2B deletion constitutes the strongest and unique prognostic factor of chemoresistance to R-CHOP, regardless the COO signature, which is not overcome by a more intensified immunochemotherapy. Patients displaying this frequent genomic abnormality warrant new and dedicated therapeutic approaches.

Brain Tissue Hypoxia is a Strong Predictor of Outcome after Severe Traumatic Brain Injury Independent from Elevated Intracranial Pressure

Introduction: Low brain tissue oxygen pressure (PbtO2) is associated with worse outcome in patients with severe traumatic brain injury (TBI). However, it is unclear whether brain tissue hypoxia is merely a marker of injury severity or a predictor of prognosis, independent from intracranial pressure (ICP) and injury severity. Hypothesis: We hypothesized that brain tissue hypoxia was an independent predictor of outcome in patients wih severe TBI, irrespective of elevated ICP and of the severity of cerebral and systemic injury. Methods: This observational study was conducted at the Neurological ICU, Hospital of the University of Pennsylvania, an academic level I trauma center. Patients admitted with severe TBI who had PbtO2 and ICP monitoring were included in the study. PbtO2, ICP, mean arterial pressure (MAP) and cerebral perfusion pressure (CPP = MAP-ICP) were monitored continuously and recorded prospectively every 30 min. Using linear interpolation, duration and cumulative dose (area under the curve, AUC) of brain tissue hypoxia (PbtO2 < 15 mm Hg), elevated ICP >20 mm Hg and low CPP <60 mm Hg were calculated, and the association with outcome at hospital discharge, dichotomized as good (Glasgow Outcome Score [GOS] 4-5) vs. poor (GOS 1-3), was analyzed. Results: A total of 103 consecutive patients, monitored for an average of 5 days, was studied. Brain tissue hypoxia was observed in 66 (64%) patients despite ICP was < 20 mm Hg and CPP > 60 mm Hg (72 +/- 39% and 49 +/- 41% of brain hypoxic time, respectively). Compared with patients with good outcome, those with poor outcome had a longer duration of brain hypoxia (1.7 +/- 3.7 vs. 8.3 +/- 15.9 hrs, P<0.01), as well as a longer duration (11.5 +/- 16.5 vs. 21.6 +/- 29.6 hrs, P=0.03) and a greater cumulative dose (56 +/- 93 vs. 143 +/- 218 mm Hg*hrs, P<0.01) of elevated ICP. By multivariable logistic regression, admission Glasgow Coma Scale (OR, 0.83, 95% CI: 0.70-0.99, P=0.04), Marshall CT score (OR 2.42, 95% CI: 1.42-4.11, P<0.01), APACHE II (OR 1.20, 95% CI: 1.03-1.43, P=0.03), and the duration of brain tissue hypoxia (OR 1.13; 95% CI: 1.01-1.27; P=0.04) were all significantly associated with poor outcome. No independent association was found between the AUC for elevated ICP and outcome (OR 1.01, 95% CI 0.97-1.02, P=0.11) in our prospective cohort. Conclusions: In patients with severe TBI, brain tissue hypoxia is frequent, despite normal ICP and CPP, and is associated with poor outcome, independent of intracranial hypertension and the severity of cerebral and systemic injury. Our findings indicate that PbtO2 is a strong physiologic prognostic marker after TBI. Further study is warranted to examine whether PbtO2-directed therapy improves outcome in severely head-injured patients .

Mutational and computational analysis of the alpha(1b)-adrenergic receptor. Involvement of basic and hydrophobic residues in receptor activation and G protein coupling.

To investigate their role in receptor coupling to G(q), we mutated all basic amino acids and some conserved hydrophobic residues of the cytosolic surface of the alpha(1b)-adrenergic receptor (AR). The wild type and mutated receptors were expressed in COS-7 cells and characterized for their ligand binding properties and ability to increase inositol phosphate accumulation. The experimental results have been interpreted in the context of both an ab initio model of the alpha(1b)-AR and of a new homology model built on the recently solved crystal structure of rhodopsin. Among the twenty-three basic amino acids mutated only mutations of three, Arg(254) and Lys(258) in the third intracellular loop and Lys(291) at the cytosolic extension of helix 6, markedly impaired the receptor-mediated inositol phosphate production. Additionally, mutations of two conserved hydrophobic residues, Val(147) and Leu(151) in the second intracellular loop had significant effects on receptor function. The functional analysis of the receptor mutants in conjunction with the predictions of molecular modeling supports the hypothesis that Arg(254), Lys(258), as well as Leu(151) are directly involved in receptor-G protein interaction and/or receptor-mediated activation of the G protein. In contrast, the residues belonging to the cytosolic extensions of helices 3 and 6 play a predominant role in the activation process of the alpha(1b)-AR. These findings contribute to the delineation of the molecular determinants of the alpha(1b)-AR/G(q) interface.

Regulation of glucose transporter expression in cardiac myocytes: p38 MAPK is a strong inducer of GLUT4.

OBJECTIVE: In vivo differentiation of cardiac myocytes is associated with downregulation of the glucose transporter isoform GLUT1 and upregulation of the isoform GLUT4. Adult rat cardiomyocytes in primary culture undergo spontaneous dedifferentiation, followed by spreading and partial redifferentiation, which can be influenced by growth factors. We used this model to study the signaling mechanisms modifying the expression of GLUT4 in cardiac myocytes. RESULTS: Adult rat cardiomyocytes in primary culture exhibited spontaneous upregulation of GLUT1 and downregulation of GLUT4, suggesting resumption of a fetal program of GLUT gene expression. Treatment with IGF-1 and, to a minor extent, FGF-2 resulted in restored expression of GLUT4 protein and mRNA. Activation of p38 MAPK mediated the increased expression of GLUT4 in response to IGF-1. Transient transfection experiments in neonatal cardiac myocytes confirmed that p38 MAPK could activate the glut4 promoter. Electrophoretic mobility shift assay in adult rat cardiomyocytes and transient transfection experiments in neonatal cardiac myocytes indicated that MEF2 was the main transcription factor transducing the effect of p38 MAPK activation on the glut4 promoter. CONCLUSION: Spontaneous dedifferentiation of adult rat cardiomyocytes in vitro is associated with downregulation of GLUT4, which can be reversed by treatment with IGF-1. The effect of IGF-1 is mediated by the p38 MAPK/MEF2 axis, which is a strong inducer of GLUT4 expression.

The impact of cyclin D1 mRNA isoforms, morphology and p53 in mantle cell lymphoma: p53 alterations and blastoid morphology are strong predictors of a high proliferation index.

BACKGROUND: Mantle cell lymphoma is a clinically heterogeneous disease characterized by overexpression of cyclin D1 protein. Blastoid morphology, high proliferation, and secondary genetic aberrations are markers of aggressive behavior. Expression profiling of mantle cell lymphoma revealed that predominance of the 3'UTR-deficient, short cyclin D1 mRNA isoform was associated with high cyclin D1 levels, a high "proliferation signature" and poor prognosis. DESIGN AND METHODS: Sixty-two cases of mantle cell lymphoma were analyzed for cyclin D1 mRNA isoforms and total cyclin D1 levels by real-time reverse transcriptase polymerase chain reaction, and TP53 alterations were assessed by immunohistochemistry and molecular analysis. Results were correlated with proliferation index and clinical outcome. RESULTS: Predominance of the short cyclin D1 mRNA was found in 14 (23%) samples, including four with complete loss of the standard transcript. TP53 alterations were found in 15 (24%) cases. Predominance of 3'UTR-deficient mRNA was significantly associated with high cyclin D1 mRNA levels (P=0.009) and more commonly found in blastoid mantle cell lymphoma (5/11, P=0.060) and cases with a proliferation index of >20% (P=0.026). Both blastoid morphology (11/11, P<0.001) and TP53 alterations (15/15, P<0.001) were significantly correlated with a high proliferation index. A proliferation index of 10% was determined to be a significant threshold for survival in multivariate analysis (P=0.01). CONCLUSIONS: TP53 alterations are strongly associated with a high proliferation index and aggressive behavior in mantle cell lymphoma. Predominance of the 3'UTR-deficient transcript correlates with higher cyclin D1 levels and may be a secondary contributing factor to high proliferation, but failed to reach prognostic significance in this study.

Strong smooth muscle differentiation is dependent on myocardin gene amplification in most human retroperitoneal leiomyosarcomas.

Myocardin (MYOCD), a serum response factor (SRF) transcriptional cofactor, is essential for cardiac and smooth muscle development and differentiation. We show here by array-based comparative genomic hybridization, fluorescence in situ hybridization, and expression analysis approaches that MYOCD gene is highly amplified and overexpressed in human retroperitoneal leiomyosarcomas (LMS), a very aggressive well-differentiated tumor. MYOCD inactivation by shRNA in a human LMS cell line with MYOCD locus amplification leads to a dramatic decrease of smooth muscle differentiation and strongly reduces cell migration. Moreover, forced MYOCD expression in three undifferentiated sarcoma cell lines and in one liposarcoma cell line confers a strong smooth muscle differentiation phenotype and increased migration abilities. Collectively, these results show that human retroperitoneal LMS differentiation is dependent on MYOCD amplification/overexpression, suggesting that in these well-differentiated LMS, differentiation could be a consequence of an acquired genomic alteration. In this hypothesis, these tumors would not necessarily derive from cells initially committed to smooth muscle differentiation. These data also provide new insights on the cellular origin of these sarcomas and on the complex connections between oncogenesis and differentiation in mesenchymal tumors.

Numerical simulation of gel electrophoresis of DNA knots in weak and strong electric fields.

Gel electrophoresis allows one to separate knotted DNA (nicked circular) of equal length according to the knot type. At low electric fields, complex knots, being more compact, drift faster than simpler knots. Recent experiments have shown that the drift velocity dependence on the knot type is inverted when changing from low to high electric fields. We present a computer simulation on a lattice of a closed, knotted, charged DNA chain drifting in an external electric field in a topologically restricted medium. Using a Monte Carlo algorithm, the dependence of the electrophoretic migration of the DNA molecules on the knot type and on the electric field intensity is investigated. The results are in qualitative and quantitative agreement with electrophoretic experiments done under conditions of low and high electric fields.

Parvalbumin interneurons mediate neuronal circuitry-neurogenesis coupling in the adult hippocampus.

Using immunohistology, electron microscopy, electrophysiology and optogenetics, we found that proliferating adult mouse hippocampal neural precursors received immature GABAergic synaptic inputs from parvalbumin-expressing interneurons. Recently shown to suppress adult quiescent neural stem cell activation, parvalbumin interneuron activation promoted newborn neuronal progeny survival and development. Our results suggest a niche mechanism involving parvalbumin interneurons that couples local circuit activity to the diametric regulation of two critical early phases of adult hippocampal neurogenesis.

Low E2F1 transcript levels are a strong determinant of favorable breast cancer outcome.

INTRODUCTION: We investigated whether mRNA levels of E2F1, a key transcription factor involved in proliferation, differentiation and apoptosis, could be used as a surrogate marker for the determination of breast cancer outcome. METHODS: E2F1 and other proliferation markers were measured by quantitative RT-PCR in 317 primary breast cancer patients from the Stiftung Tumorbank Basel. Correlations to one another as well as to the estrogen receptor and ERBB2 status and clinical outcome were investigated. Results were validated and further compared with expression-based prognostic profiles using The Netherlands Cancer Institute microarray data set reported by Fan and colleagues. RESULTS: E2F1 mRNA expression levels correlated strongly with the expression of other proliferation markers, and low values were mainly found in estrogen receptor-positive and ERBB2-negative phenotypes. Patients with low E2F1-expressing tumors were associated with favorable outcome (hazard ratio = 4.3 (95% confidence interval = 1.8-9.9), P = 0.001). These results were consistent in univariate and multivariate Cox analyses, and were successfully validated in The Netherlands Cancer Institute data set. Furthermore, E2F1 expression levels correlated well with the 70-gene signature displaying the ability of selecting a common subset of patients at good prognosis. Breast cancer patients' outcome was comparably predictable by E2F1 levels, by the 70-gene signature, by the intrinsic subtype gene classification, by the wound response signature and by the recurrence score. CONCLUSION: Assessment of E2F1 at the mRNA level in primary breast cancer is a strong determinant of breast cancer patient outcome. E2F1 expression identified patients at low risk of metastasis irrespective of the estrogen receptor and ERBB2 status, and demonstrated similar prognostic performance to different gene expression-based predictors.

Strong EBV-specific CD8+ T-cell response in patients with early multiple sclerosis

Epstein-Barr virus (EBV) has been associated with multiple sclerosis (MS), however, most studies examining the relationship between the virus and the disease have been based on serologies, and if EBV is linked to MS, CD8+ T cells are likely to be involved as they are important both in MS pathogenesis and in controlling viruses. We hypothesized that valuable information on the link between MS and EBV would be ascertained from the study of frequency and activation levels of EBV-specific CD8+ T cells in different categories of MS patients and control subjects. We investigated EBV-specific cellular immune responses using proliferation and enzyme linked immunospot assays, and humoral immune responses by analysis of anti-EBV antibodies, in a cohort of 164 subjects, including 108 patients with different stages of MS, 35 with other neurological diseases and 21 healthy control subjects. Additionally, the cohort were all tested against cytomegalovirus (CMV), another neurotropic herpes virus not convincingly associated with MS, nor thought to be deleterious to the disease. We corrected all data for age using linear regression analysis over the total cohorts of EBV- and CMV-infected subjects. In the whole cohort, the rate of EBV and CMV infections were 99% and 51%, respectively. The frequency of IFN-gamma secreting EBV-specific CD8+ T cells in patients with clinically isolated syndrome (CIS) was significantly higher than that found in patients with relapsing-remitting MS (RR-MS), secondary-progressive MS, primary-progressive MS, patients with other neurological diseases and healthy controls. The shorter the interval between MS onset and our assays, the more intense was the EBV-specific CD8+ T-cell response. Confirming the above results, we found that EBV-specific CD8+ T-cell responses decreased in 12/13 patients with CIS followed prospectively for 1.0 +/- 0.2 years. In contrast, there was no difference between categories for EBV-specific CD4+ T cell, or for CMV-specific CD4+ and CD8+ T-cell responses. Anti-EBV-encoded nuclear antigen-1 (EBNA-1)-specific antibodies correlated with EBV-specific CD8+ T cells in patients with CIS and RR-MS. However, whereas EBV-specific CD8+ T cells were increased the most in early MS, EBNA-1-specific antibodies were increased in early as well as in progressive forms of MS. Our data show high levels of CD8+ T-cell activation against EBV--but not CMV--early in the course of MS, which support the hypothesis that EBV might be associated with the onset of this disease.