Magnus Fontes

Microarray-based classification of a consecutive series of 121 childhood acute leukemias: prediction of leukemic and genetic subtype as well as of minimal residual disease status.

Gene expression analyses were performed on 121 consecutive childhood leukemias (87 B-lineage acute lymphoblastic leukemias (ALLs), 11 T-cell ALLs and 23 acute myeloid leukemias (AMLs)), investigated during an 8-year period at a single center. The supervised learning algorithm k-nearest neighbor was utilized to build gene expression predictors that could classify the ALLs/AMLs according to clinically important subtypes with high accuracy. Validation experiments in an independent data set verified the high prediction accuracies of our classifiers. B-lineage ALLs with uncharacteristic cytogenetic aberrations or with a normal karyotype displayed heterogeneous gene expression profiles, resulting in low prediction accuracies. Minimal residual disease status (MRD) in T-cell ALLs with a high (>0.1%) MRD at day 29 could be classified with 100% accuracy already at the time of diagnosis. In pediatric leukemias with uncharacteristic cytogenetic aberrations or with a normal karyotype, unsupervised analysis identified two novel subgroups: one consisting mainly of cases remaining in complete remission (CR) and one containing a few patients in CR and all but one of the patients who relapsed. This study of a consecutive series of childhood leukemias confirms and extends further previous reports demonstrating that global gene expression profiling provides a valuable tool for genetic and clinical classification of childhood leukemias.

Transcriptional regulation of aquaporins in accessions of Arabidopsis in response to drought stress.

Aquaporins facilitate water transport over cellular membranes, and are therefore believed to play an important role in water homeostasis. In higher plants aquaporin-like proteins, also called major intrinsic proteins (MIPs), are divided into five subfamilies. We have previously shown that MIP transcription in Arabidopsis thaliana is generally downregulated in leaves upon drought stress, apart from two members of the plasma membrane intrinsic protein (PIP) subfamily, AtPIP1;4 and AtPIP2;5, which are upregulated. In order to assess whether this regulation is general or accession-specific we monitored the gene expression of all PIPs in five Arabidopsis accessions. The overall drought regulation of PIPs was well conserved for all five accessions tested, suggesting a general and fundamental physiological role of this drought response. In addition, significant differences among accessions were identified for transcripts of three PIP genes. Principal component analysis showed that most of the PIP transcriptional variation during drought stress could be explained by one variable linked to leaf water content. Promoter-GUS constructs of AtPIP1;4, AtPIP2;5 and also AtPIP2;6, which is unresponsive to drought stress, had distinct expression patterns concentrated in the base of the leaf petioles and parts of the flowers. The presence of drought stress response elements within the 1.6-kb promoter regions of AtPIP1;4 and AtPIP2;5 was demonstrated by comparing transcription of the promoter reporter construct and the endogenous gene upon drought stress. Analysis by ATTED-II and other web-based bioinformatical tools showed that several of the MIPs downregulated upon drought are strongly co-expressed, whereas AtPIP1;4, AtPIP2;5 and AtPIP2;6 are not co-expressed.

Gene expression profiling of leukemic cell lines reveals conserved molecular signatures among subtypes with specific genetic aberrations

Hematologic malignancies are characterized by fusion genes of biological/clinical importance. Immortalized cell lines with such aberrations are today widely used to model different aspects of leukemogenesis. Using cDNA microarrays, we determined the gene expression profiles of 40 cell lines as well as of primary leukemias harboring 11q23/MLL rearrangements, t(1;19)[TCF3/PBX1], t(12;21)[ETV6/RUNX1], t(8;21)[RUNX1/CBFA2T1], t(8;14)[IGH@/MYC], t(8;14)[TRA@/MYC], t(9;22)[BCR/ABL1], t(10;11)[PICALM/MLLT10], t(15;17)[PML/RARA], or inv(16)[CBFB/MYH11]. Unsupervised classification revealed that hematopoietic cell lines of diverse origin, but with the same primary genetic changes, segregated together, suggesting that pathogenetically important regulatory networks remain conserved despite numerous passages. Moreover, primary leukemias cosegregated with cell lines carrying identical genetic rearrangements, further supporting that critical regulatory pathways remain intact in hematopoietic cell lines. Transcriptional signatures correlating with clinical subtypes/primary genetic changes were identified and annotated based on their biological/molecular properties and chromosomal localization. Furthermore, the expression profile of tyrosine kinase-encoding genes was investigated, identifying several differentially expressed members, segregating with primary genetic changes, which may be targeted with tyrosine kinase inhibitors. The identified conserved signatures are likely to reflect regulatory networks of importance for the transforming abilities of the primary genetic changes and offer important pathogenetic insights as well as a number of targets for future rational drug design.

Early changes in the hypothalamic region in prodromal Huntington disease revealed by MRI analysis.

Huntington disease (HD) is a fatal neurodegenerative disorder caused by an expanded CAG repeat. Its length can be used to estimate the time of clinical diagnosis, which is defined by overt motor symptoms. Non-motor symptoms begin before motor onset, and involve changes in hypothalamus-regulated functions such as sleep, emotion and metabolism. Therefore we hypothesized that hypothalamic changes occur already prior to the clinical diagnosis. We performed voxel-based morphometry and logistic regression analyses of cross-sectional MR images from 220 HD gene carriers and 75 controls in the Predict-HD study. We show that changes in the hypothalamic region are detectable before clinical diagnosis and that its grey matter contents alone are sufficient to distinguish HD gene carriers from control cases. In conclusion, our study shows, for the first time, that alterations in grey matter contents in the hypothalamic region occur at least a decade before clinical diagnosis in HD using MRI.

Integrative analysis of gene expression and copy number alterations using canonical correlation analysis

BackgroundWith the rapid development of new genetic measurement methods, several types of genetic alterations can be quantified in a high-throughput manner. While the initial focus has been on investigating each data set separately, there is an increasing interest in studying the correlation structure between two or more data sets. Multivariate methods based on Canonical Correlation Analysis (CCA) have been proposed for integrating paired genetic data sets. The high dimensionality of microarray data imposes computational difficulties, which have been addressed for instance by studying the covariance structure of the data, or by reducing the number of variables prior to applying the CCA. In this work, we propose a new method for analyzing high-dimensional paired genetic data sets, which mainly emphasizes the correlation structure and still permits efficient application to very large data sets. The method is implemented by translating a regularized CCA to its dual form, where the computational complexity depends mainly on the number of samples instead of the number of variables. The optimal regularization parameters are chosen by cross-validation. We apply the regularized dual CCA, as well as a classical CCA preceded by a dimension-reducing Principal Components Analysis (PCA), to a paired data set of gene expression changes and copy number alterations in leukemia.ResultsUsing the correlation-maximizing methods, regularized dual CCA and PCA+CCA, we show that without pre-selection of known disease-relevant genes, and without using information about clinical class membership, an exploratory analysis singles out two patient groups, corresponding to well-known leukemia subtypes. Furthermore, the variables showing the highest relevance to the extracted features agree with previous biological knowledge concerning copy number alterations and gene expression changes in these subtypes. Finally, the correlation-maximizing methods are shown to yield results which are more biologically interpretable than those resulting from a covariance-maximizing method, and provide different insight compared to when each variable set is studied separately using PCA.ConclusionsWe conclude that regularized dual CCA as well as PCA+CCA are useful methods for exploratory analysis of paired genetic data sets, and can be efficiently implemented also when the number of variables is very large.

Group Selection and Contribution of Minority Variants during Virus Adaptation Determines Virus Fitness and Phenotype.

Understanding how a pathogen colonizes and adapts to a new host environment is a primary aim in studying emerging infectious diseases. Adaptive mutations arise among the thousands of variants generated during RNA virus infection, and identifying these variants will shed light onto how changes in tropism and species jumps can occur. Here, we adapted Coxsackie virus B3 to a highly permissive and less permissive environment. Using deep sequencing and bioinformatics, we identified a multi-step adaptive process to adaptation involving residues in the receptor footprints that correlated with receptor availability and with increase in virus fitness in an environment-specific manner. We show that adaptation occurs by selection of a dominant mutation followed by group selection of minority variants that together, confer the fitness increase observed in the population, rather than selection of a single dominant genotype.

Approximate geodesic distances reveal biologically relevant structures in microarray data

MOTIVATION Genome-wide gene expression measurements, as currently determined by the microarray technology, can be represented mathematically as points in a high-dimensional gene expression space. Genes interact with each other in regulatory networks, restricting the cellular gene expression profiles to a certain manifold, or surface, in gene expression space. To obtain knowledge about this manifold, various dimensionality reduction methods and distance metrics are used. For data points distributed on curved manifolds, a sensible distance measure would be the geodesic distance along the manifold. In this work, we examine whether an approximate geodesic distance measure captures biological similarities better than the traditionally used Euclidean distance. RESULTS We computed approximate geodesic distances, determined by the Isomap algorithm, for one set of lymphoma and one set of lung cancer microarray samples. Compared with the ordinary Euclidean distance metric, this distance measure produced more instructive, biologically relevant, visualizations when applying multidimensional scaling. This suggests the Isomap algorithm as a promising tool for the interpretation of microarray data. Furthermore, the results demonstrate the benefit and importance of taking nonlinearities in gene expression data into account.

Identification of ETV6-RUNX1 -like and DUX4 -rearranged subtypes in paediatric B-cell precursor acute lymphoblastic leukaemia

Fusion genes are potent driver mutations in cancer. In this study, we delineate the fusion gene landscape in a consecutive series of 195 paediatric B-cell precursor acute lymphoblastic leukaemia (BCP ALL). Using RNA sequencing, we find in-frame fusion genes in 127 (65%) cases, including 27 novel fusions. We describe a subtype characterized by recurrent IGH-DUX4 or ERG-DUX4 fusions, representing 4% of cases, leading to overexpression of DUX4 and frequently co-occurring with intragenic ERG deletions. Furthermore, we identify a subtype characterized by an ETV6-RUNX1-like gene-expression profile and coexisting ETV6 and IKZF1 alterations. Thus, this study provides a detailed overview of fusion genes in paediatric BCP ALL and adds new pathogenetic insights, which may improve risk stratification and provide therapeutic options for this disease.

The correlation pattern of acquired copy number changes in 164 ETV6/RUNX1-positive childhood acute lymphoblastic leukemias

The ETV6/RUNX1 fusion gene, present in 25% of B-lineage childhood acute lymphoblastic leukemia (ALL), is thought to represent an initiating event, which requires additional genetic changes for leukemia development. To identify additional genetic alterations, 24 ETV6/RUNX1-positive ALLs were analyzed using 500K single nucleotide polymorphism arrays. The results were combined with previously published data sets, allowing us to ascertain genomic copy number aberrations (CNAs) in 164 cases. In total, 45 recurrent CNAs were identified with an average number of 3.5 recurrent changes per case (range 0-13). Twenty-six percent of cases displayed a set of recurrent CNAs identical to that of other cases in the data set. The majority (74%), however, displayed a unique pattern of recurrent CNAs, indicating a large heterogeneity within this ALL subtype. As previously demonstrated, alterations targeting genes involved in B-cell development were common (present in 28% of cases). However, the combined analysis also identified alterations affecting nuclear hormone response (24%) to be a characteristic feature of ETV6/RUNX1-positive ALL. Studying the correlation pattern of the CNAs allowed us to highlight significant positive and negative correlations between specific aberrations. Furthermore, oncogenetic tree models identified ETV6, CDKN2A/B, PAX5, del(6q) and +16 as possible early events in the leukemogenic process.

The ABCs of viral hepatitis that define biomarker signatures of acute viral hepatitis

Viral hepatitis is the leading cause of liver disease worldwide and can be caused by several agents, including hepatitis A (HAV), B (HBV), and C (HCV) virus. We employed multiplexed protein immune assays to identify biomarker signatures of viral hepatitis in order to define unique and common responses for three different acute viral infections of the liver. We performed multianalyte profiling, measuring the concentrations of 182 serum proteins obtained from acute HAV‐ (18), HBV‐ (18), and HCV‐infected (28) individuals, recruited as part of a hospital‐based surveillance program in Cairo, Egypt. Virus‐specific biomarker signatures were identified and validation was performed using a unique patient population. A core signature of 46 plasma proteins was commonly modulated in all three infections, as compared to healthy controls. Principle component analysis (PCA) revealed a host response based upon 34 proteins, which could distinguish HCV patients from HAV‐ and HBV‐infected individuals or healthy controls. When HAV and HBV groups were compared directly, 34 differentially expressed serum proteins allowed the separation of these two patient groups. A validation study was performed on an additional 111 patients, confirming the relevance of our initial findings, and defining the 17 analytes that reproducibly segregated the patient populations. Conclusions: This combined discovery and biomarker validation approach revealed a previously unrecognized virus‐specific induction of host proteins. The identification of hepatitis virus specific signatures provides a foundation for functional studies and the identification of potential correlates of viral clearance. (Hepatology 2014;59:1273‐1282)