Global copy number profiling of cancer genomes

Summary: In this article, we introduce a robust and efficient strategy for deriving global and allele-specific copy number alternations (CNA) from cancer whole exome sequencing data based on Log R ratios and B-allele frequencies. Applying the approach to the analysis of over 200 skin cancer samples, we demonstrate its utility for discovering distinct CNA events and for deriving ancillary information such as tumor purity

Medoidshift clustering applied to genomic bulk tumor data.

Despite the enormous medical impact of cancers and intensive study of their biology, detailed characterization of tumor growth and development remains elusive. This difficulty occurs in large part because of enormous heterogeneity in the molecular mechanisms of cancer progression, both tumor-to-tumor and cell-to-cell in single tumors. Advances in genomic technologies, especially at the single-cell level, are improving the situation, but these approaches are held back by limitations of the biotechnologies for gathering genomic data from heterogeneous cell populations and the computational methods for making sense of those data. One popular way to gain the advantages of whole-genome methods without the cost of single-cell genomics has been the use of computational deconvolution (unmixing) methods to reconstruct clonal heterogeneity from bulk genomic data. These methods, too, are limited by the difficulty of inferring genomic profiles of rare or subtly varying clonal subpopulations from bulk data, a problem that can be computationally reduced to that of reconstructing the geometry of point clouds of tumor samples in a genome space. Here, we present a new method to improve that reconstruction by better identifying subspaces corresponding to tumors produced from mixtures of distinct combinations of clonal subpopulations. We develop a nonparametric clustering method based on medoidshift clustering for identifying subgroups of tumors expected to correspond to distinct trajectories of evolutionary progression. We show on synthetic and real tumor copy-number data that this new method substantially improves our ability to resolve discrete tumor subgroups, a key step in the process of accurately deconvolving tumor genomic data and inferring clonal heterogeneity from bulk data

Molecular Portraits of Cancer Evolution and Ecology

Research on the molecular lesions that drive cancers holds the translational promise of unmasking distinct disease subtypes in otherwise pathologically identical patients. Yet clinical adoption is hindered by the reproducibility crisis for cancer biomarkers. In this thesis, a novel metric uncovered transcriptional diversity within individual non-small cell lung cancers, driven by chromosomal instability. Existing prognostic biomarkers were confounded by tumour sampling bias, arising from this diversity, in ~50% of patients assessed. An atlas of consistently expressed genes was derived to address this diagnostic challenge, yielding a clonal biomarker robust to sampling bias. This diagnostic based on cancer evolutionary principles maintained prognostic value in a metaanalysis of >900 patients, and over known risk factors in stage I disease, motivating further development as a clinical assay. Next, in situ RNA profiles of immune, fibroblast and endothelial cell subsets were generated from cancerous and adjacent non-malignant lung tissue. The phenotypic adaptation of stromal cells in the tumour microenvironment undermined the performance of existing molecular signatures for cell-type enumeration. Transcriptome-wide analysis delineated ~10% of genes displaying cell-type-specific expression, paving the way for high-fidelity signatures for the accurate digital dissection of tumour ecology. Lastly, the impact of branching, Darwinian evolution on the detection of epistatic interactions was evaluated in a pan-cancer analysis. The clonal status of driver genes was associated with the proportion of significant epistatic findings in 44-78% of the cancer-types assessed. Integrating the clonal architecture of tumours in future analyses could help decipher evolutionary dependencies. This work provides pragmatic solutions for refining molecular portraits of cancer in the light of their evolutionary and ecological features, moving the needle for precision cancer diagnostics

Methods and practice of detecting selection in human cancers

Cancer development and progression is an evolutionary process, understanding these evolutionary dynamics is important for treatment and diagnosis as how a cancer evolves determines its future prognosis. This thesis focuses on elucidating selective evolutionary pressures in cancers and somatic tissues using population genetics models and cancer genomics data. First a model for the expected diversity in the absence of selection was developed. This neutral model of evolution predicts that under neutrality the frequency of subclonal mutations is expected to follow a power law distribution. Surprisingly more than 30% of cancer across multiple cohorts fitted this model. The next part of the thesis develops models to explore the effects of selection given these should be observable as deviations from the neutral prediction. For this I developed two approaches. The first approach investigated selection at the level of individual samples and showed that a characteristic pattern of clusters of mutations is observed in deep sequencing experiments. Using a mathematical model, information encoded within these clusters can be used to measure the relative fitness of subclones and the time they emerge during tumour evolution. With this I observed strikingly high fitness advantages for subclones of above 20%. The second approach enables measuring recurrent patterns of selection in cohorts of sequenced cancers using dN/dS, the ratio of non-synonymous to synonymous mutations, a method originally developed for molecular species evolution. This approach demonstrates how selection coefficients can be extracted by combining measurements of dN/dS with the size of mutational lineages. With this approach selection coefficients were again observed to be strikingly high. Finally I looked at population dynamics in normal colonic tissue given that many mutations accumulate in physiologically normal tissue. I found that the current view of stem cell dynamics was unable to explain sequencing data from individual colonic crypts. Some new models were proposed that introduce a longer time scale evolution that suppresses the accumulation of mutations which appear consistent with the data

A telescope for the RNA universe : novel bioinformatic approaches to analyze RNA sequencing data

In this thesis I focus on the application of bioinformatics to analyze RNA. The type of experimental data of interest is sequencing data generated with various Next Generation Sequencing technique: nuclear RNA, cytoplasmic RNA, captured polyadenylated RNA fragments, etc. I highlight the necessity in developing new tools (e.g., to analyze nuclear RNA) and give a showcase example of implementing such tool and showing its usability on a real biological experiment. The thesis also covers existing tools to perform various types of RNA analysis and shows how these tools can be twigged and expanded to answer certain biological questions (e.g., studying changes in RNA specific to human aging). I also show how current bioinformatic approaches can be used in a particularly complex study such as investigating cancer (in this thesis, breast cancer) pathogenesis.UBL - phd migration 201