Roland Eils

Signatures of mutational processes in human cancer

All cancers are caused by somatic mutations; however, understanding of the biological processes generating these mutations is limited. The catalogue of somatic mutations from a cancer genome bears the signatures of the mutational processes that have been operative. Here we analysed 4,938,362 mutations from 7,042 cancers and extracted more than 20 distinct mutational signatures. Some are present in many cancer types, notably a signature attributed to the APOBEC family of cytidine deaminases, whereas others are confined to a single cancer class. Certain signatures are associated with age of the patient at cancer diagnosis, known mutagenic exposures or defects in DNA maintenance, but many are of cryptic origin. In addition to these genome-wide mutational signatures, hypermutation localized to small genomic regions, ‘kataegis’, is found in many cancer types. The results reveal the diversity of mutational processes underlying the development of cancer, with potential implications for understanding of cancer aetiology, prevention and therapy.

Systemic Spread Is an Early Step in Breast Cancer

It is widely accepted that metastasis is a late event in cancer progression. Here, however, we show that tumor cells can disseminate systemically from earliest epithelial alterations in HER-2 and PyMT transgenic mice and from ductal carcinoma in situ in women. Wild-type mice transplanted with single premalignant HER-2 transgenic glands displayed disseminated tumor cells and micrometastasis in bone marrow and lungs. The number of disseminated cancer cells and their karyotypic abnormalities were similar for small and large tumors in patients and mouse models. When activated by bone marrow transplantation into wild-type recipients, 80 early-disseminated cancer cells sufficed to induce lethal carcinosis. Therefore, release from dormancy of early-disseminated cancer cells may frequently account for metachronous metastasis.

Phenotypic profiling of the human genome by time-lapse microscopy reveals cell division genes.

Despite our rapidly growing knowledge about the human genome, we do not know all of the genes required for some of the most basic functions of life. To start to fill this gap we developed a high-throughput phenotypic screening platform combining potent gene silencing by RNA interference, time-lapse microscopy and computational image processing. We carried out a genome-wide phenotypic profiling of each of the ∼21,000 human protein-coding genes by two-day live imaging of fluorescently labelled chromosomes. Phenotypes were scored quantitatively by computational image processing, which allowed us to identify hundreds of human genes involved in diverse biological functions including cell division, migration and survival. As part of the Mitocheck consortium, this study provides an in-depth analysis of cell division phenotypes and makes the entire high-content data set available as a resource to the community.

Three-Dimensional Maps of All Chromosomes in Human Male Fibroblast Nuclei and Prometaphase Rosettes

Studies of higher-order chromatin arrangements are an essential part of ongoing attempts to explore changes in epigenome structure and their functional implications during development and cell differentiation. However, the extent and cell-type-specificity of three-dimensional (3D) chromosome arrangements has remained controversial. In order to overcome technical limitations of previous studies, we have developed tools that allow the quantitative 3D positional mapping of all chromosomes simultaneously. We present unequivocal evidence for a probabilistic 3D order of prometaphase chromosomes, as well as of chromosome territories (CTs) in nuclei of quiescent (G0) and cycling (early S-phase) human diploid fibroblasts (46, XY). Radial distance measurements showed a probabilistic, highly nonrandom correlation with chromosome size: small chromosomes—independently of their gene density—were distributed significantly closer to the center of the nucleus or prometaphase rosette, while large chromosomes were located closer to the nuclear or rosette rim. This arrangement was independently confirmed in both human fibroblast and amniotic fluid cell nuclei. Notably, these cell types exhibit flat-ellipsoidal cell nuclei, in contrast to the spherical nuclei of lymphocytes and several other human cell types, for which we and others previously demonstrated gene-density-correlated radial 3D CT arrangements. Modeling of 3D CT arrangements suggests that cell-type-specific differences in radial CT arrangements are not solely due to geometrical constraints that result from nuclear shape differences. We also found gene-density-correlated arrangements of higher-order chromatin shared by all human cell types studied so far. Chromatin domains, which are gene-poor, form a layer beneath the nuclear envelope, while gene-dense chromatin is enriched in the nuclear interior. We discuss the possible functional implications of this finding.

Dissecting the genomic complexity underlying medulloblastoma

Medulloblastoma is an aggressively growing tumour, arising in the cerebellum or medulla/brain stem. It is the most common malignant brain tumour in children, and shows tremendous biological and clinical heterogeneity. Despite recent treatment advances, approximately 40% of children experience tumour recurrence, and 30% will die from their disease. Those who survive often have a significantly reduced quality of life. Four tumour subgroups with distinct clinical, biological and genetic profiles are currently identified. WNT tumours, showing activated wingless pathway signalling, carry a favourable prognosis under current treatment regimens. SHH tumours show hedgehog pathway activation, and have an intermediate prognosis. Group 3 and 4 tumours are molecularly less well characterized, and also present the greatest clinical challenges. The full repertoire of genetic events driving this distinction, however, remains unclear. Here we describe an integrative deep-sequencing analysis of 125 tumour–normal pairs, conducted as part of the International Cancer Genome Consortium (ICGC) PedBrain Tumor Project. Tetraploidy was identified as a frequent early event in Group 3 and 4 tumours, and a positive correlation between patient age and mutation rate was observed. Several recurrent mutations were identified, both in known medulloblastoma-related genes (CTNNB1, PTCH1, MLL2, SMARCA4) and in genes not previously linked to this tumour (DDX3X, CTDNEP1, KDM6A, TBR1), often in subgroup-specific patterns. RNA sequencing confirmed these alterations, and revealed the expression of what are, to our knowledge, the first medulloblastoma fusion genes identified. Chromatin modifiers were frequently altered across all subgroups. These findings enhance our understanding of the genomic complexity and heterogeneity underlying medulloblastoma, and provide several potential targets for new therapeutics, especially for Group 3 and 4 patients.

From latent disseminated cells to overt metastasis: Genetic analysis of systemic breast cancer progression

According to the present view, metastasis marks the end in a sequence of genomic changes underlying the progression of an epithelial cell to a lethal cancer. Here, we aimed to find out at what stage of tumor development transformed cells leave the primary tumor and whether a defined genotype corresponds to metastatic disease. To this end, we isolated single disseminated cancer cells from bone marrow of breast cancer patients and performed single-cell comparative genomic hybridization. We analyzed disseminated tumor cells from patients after curative resection of the primary tumor (stage M0), as presumptive progenitors of manifest metastasis, and from patients with manifest metastasis (stage M1). Their genomic data were compared with those from microdissected areas of matched primary tumors. Disseminated cells from M0-stage patients displayed significantly fewer chromosomal aberrations than primary tumors or cells from M1-stage patients (P < 0.008 and P < 0.0001, respectively), and their aberrations appeared to be randomly generated. In contrast, primary tumors and M1 cells harbored different and characteristic chromosomal imbalances. Moreover, applying machine-learning methods for the classification of the genotypes, we could correctly identify the presence or absence of metastatic disease in a patient on the basis of a single-cell genome. We suggest that in breast cancer, tumor cells may disseminate in a far less progressed genomic state than previously thought, and that they acquire genomic aberrations typical of metastatic cells thereafter. Thus, our data challenge the widely held view that the precursors of metastasis are derived from the most advanced clone within the primary tumor.

Recruitment and activation of a lipid kinase by hepatitis C virus NS5A is essential for integrity of the membranous replication compartment

Hepatitis C virus (HCV) is a major causative agent of chronic liver disease in humans. To gain insight into host factor requirements for HCV replication, we performed a siRNA screen of the human kinome and identified 13 different kinases, including phosphatidylinositol-4 kinase III alpha (PI4KIIIα), as being required for HCV replication. Consistent with elevated levels of the PI4KIIIα product phosphatidylinositol-4-phosphate (PI4P) detected in HCV-infected cultured hepatocytes and liver tissue from chronic hepatitis C patients, the enzymatic activity of PI4KIIIα was critical for HCV replication. Viral nonstructural protein 5A (NS5A) was found to interact with PI4KIIIα and stimulate its kinase activity. The absence of PI4KIIIα activity induced a dramatic change in the ultrastructural morphology of the membranous HCV replication complex. Our analysis suggests that the direct activation of a lipid kinase by HCV NS5A contributes critically to the integrity of the membranous viral replication complex.

Recurrent somatic alterations of FGFR1 and NTRK2 in pilocytic astrocytoma

Pilocytic astrocytoma, the most common childhood brain tumor, is typically associated with mitogen-activated protein kinase (MAPK) pathway alterations. Surgically inaccessible midline tumors are therapeutically challenging, showing sustained tendency for progression and often becoming a chronic disease with substantial morbidities. Here we describe whole-genome sequencing of 96 pilocytic astrocytomas, with matched RNA sequencing (n = 73), conducted by the International Cancer Genome Consortium (ICGC) PedBrain Tumor Project. We identified recurrent activating mutations in FGFR1 and PTPN11 and new NTRK2 fusion genes in non-cerebellar tumors. New BRAF-activating changes were also observed. MAPK pathway alterations affected all tumors analyzed, with no other significant mutations identified, indicating that pilocytic astrocytoma is predominantly a single-pathway disease. Notably, we identified the same FGFR1 mutations in a subset of H3F3A-mutated pediatric glioblastoma with additional alterations in the NF1 gene. Our findings thus identify new potential therapeutic targets in distinct subsets of pilocytic astrocytoma and childhood glioblastoma.

Acute myeloid leukemias with reciprocal rearrangements can be distinguished by specific gene expression profiles

Acute myeloid leukemia (AML) is a heterogeneous group of genetically defined diseases. Their classification is important with regard to prognosis and treatment. We performed microarray analyses for gene expression profiling on bone marrow samples of 37 patients with newly diagnosed AML. All cases had either of the distinct subtypes AML M2 with t(8;21), AML M3 or M3v with t(15;17), or AML M4eo with inv(16). Diagnosis was established by cytomorphology, cytogenetics, fluorescence in situ hybridization, and reverse transcriptase–PCR in every sample. By using two different strategies for microarray data analyses, this study revealed a unique correlation between AML-specific cytogenetic aberrations and gene expression profiles.

Omic data from evolved E. coli are consistent with computed optimal growth from genome‐scale models

After hundreds of generations of adaptive evolution at exponential growth, Escherichia coli grows as predicted using flux balance analysis (FBA) on genome‐scale metabolic models (GEMs). However, it is not known whether the predicted pathway usage in FBA solutions is consistent with gene and protein expression in the wild‐type and evolved strains. Here, we report that >98% of active reactions from FBA optimal growth solutions are supported by transcriptomic and proteomic data. Moreover, when E. coli adapts to growth rate selective pressure, the evolved strains upregulate genes within the optimal growth predictions, and downregulate genes outside of the optimal growth solutions. In addition, bottlenecks from dosage limitations of computationally predicted essential genes are overcome in the evolved strains. We also identify regulatory processes that may contribute to the development of the optimal growth phenotype in the evolved strains, such as the downregulation of known regulons and stringent response suppression. Thus, differential gene and protein expression from wild‐type and adaptively evolved strains supports observed growth phenotype changes, and is consistent with GEM‐computed optimal growth states.