Gene fusions created by somatic genomic rearrangements are known to play an important role in the onset and development of some cancers, such as lymphomas and sarcomas. RNA-Seq (whole transcriptome ...shotgun sequencing) is proving to be a useful tool for the discovery of novel gene fusions in cancer transcriptomes. However, algorithmic methods for the discovery of gene fusions using RNA-Seq data remain underdeveloped. We have developed deFuse, a novel computational method for fusion discovery in tumor RNA-Seq data. Unlike existing methods that use only unique best-hit alignments and consider only fusion boundaries at the ends of known exons, deFuse considers all alignments and all possible locations for fusion boundaries. As a result, deFuse is able to identify fusion sequences with demonstrably better sensitivity than previous approaches. To increase the specificity of our approach, we curated a list of 60 true positive and 61 true negative fusion sequences (as confirmed by RT-PCR), and have trained an adaboost classifier on 11 novel features of the sequence data. The resulting classifier has an estimated value of 0.91 for the area under the ROC curve. We have used deFuse to discover gene fusions in 40 ovarian tumor samples, one ovarian cancer cell line, and three sarcoma samples. We report herein the first gene fusions discovered in ovarian cancer. We conclude that gene fusions are not infrequent events in ovarian cancer and that these events have the potential to substantially alter the expression patterns of the genes involved; gene fusions should therefore be considered in efforts to comprehensively characterize the mutational profiles of ovarian cancer transcriptomes.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
We introduce PyClone, a statistical model for inference of clonal population structures in cancers. PyClone is a Bayesian clustering method for grouping sets of deeply sequenced somatic mutations ...into putative clonal clusters while estimating their cellular prevalences and accounting for allelic imbalances introduced by segmental copy-number changes and normal-cell contamination. Single-cell sequencing validation demonstrates PyClone's accuracy.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK
Human cancer cell lines are the workhorse of cancer research. Although cell lines are known to evolve in culture, the extent of the resultant genetic and transcriptional heterogeneity and its ...functional consequences remain understudied. Here we use genomic analyses of 106 human cell lines grown in two laboratories to show extensive clonal diversity. Further comprehensive genomic characterization of 27 strains of the common breast cancer cell line MCF7 uncovered rapid genetic diversification. Similar results were obtained with multiple strains of 13 additional cell lines. Notably, genetic changes were associated with differential activation of gene expression programs and marked differences in cell morphology and proliferation. Barcoding experiments showed that cell line evolution occurs as a result of positive clonal selection that is highly sensitive to culture conditions. Analyses of single-cell-derived clones demonstrated that continuous instability quickly translates into heterogeneity of the cell line. When the 27 MCF7 strains were tested against 321 anti-cancer compounds, we uncovered considerably different drug responses: at least 75% of compounds that strongly inhibited some strains were completely inactive in others. This study documents the extent, origins and consequences of genetic variation within cell lines, and provides a framework for researchers to measure such variation in efforts to support maximally reproducible cancer research.
Full text
Available for:
KISLJ, NUK, SBMB, UL, UM, UPUK
Prostate cancers (PCs) with loss of the potent tumor suppressors TP53 and RB1 exhibit poor outcomes. TP53 and RB1 also influence cell plasticity and are frequently lost in PCs with neuroendocrine ...(NE) differentiation. Therapeutic strategies that address these aggressive variant PCs are urgently needed. Using deep genomic profiling of 410 metastatic biopsies, we determine the relationships between combined TP53 and RB1 loss and PC phenotypes. Notably, 40% of TP53/RB1-deficient tumors are classified as AR-active adenocarcinomas, indicating that NE differentiation is not an obligate consequence of TP53/RB1 inactivation. A gene expression signature reflecting TP53/RB1 loss is associated with diminished responses to AR antagonists and reduced survival. These tumors exhibit high proliferation rates and evidence of elevated DNA repair processes. While tumor cells lacking TP53/RB1 are highly resistant to all single-agent therapeutics tested, the combination of PARP and ATR inhibition is found to produce significant responses, reflecting a clinically exploitable vulnerability resulting from replication stress.
Display omitted
•Prostate cancers (PCs) with TP53 and RB1 loss exhibit very poor clinical outcomes•PCs with TP53/RB1 loss exhibit stem cell-like features and loss of AR activity•Loss of TP53/RB1 does not uniformly promote neuroendocrine transdifferentiation•TP53/RB1-null PCs exhibit replication stress and respond to PARP and ATR inhibition
Nyquist et al. demonstrate that TP53 and RB1 loss in prostate carcinoma (PC) attenuates AR signaling and enhances cell proliferation but does not uniformly induce neuroendocrine phenotypes. PCs with TP53/RB1 loss resist a wide range of cancer therapeutics but respond to PARP and ATR inhibition, likely reflecting enhanced replication stress.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
Adenocarcinoma in situ and minimally invasive adenocarcinoma are the pre-invasive forms of lung adenocarcinoma. The genomic and immune profiles of these lesions are poorly understood. Here we report ...exome and transcriptome sequencing of 98 lung adenocarcinoma precursor lesions and 99 invasive adenocarcinomas. We have identified EGFR, RBM10, BRAF, ERBB2, TP53, KRAS, MAP2K1 and MET as significantly mutated genes in the pre/minimally invasive group. Classes of genome alterations that increase in frequency during the progression to malignancy are revealed. These include mutations in TP53, arm-level copy number alterations, and HLA loss of heterozygosity. Immune infiltration is correlated with copy number alterations of chromosome arm 6p, suggesting a link between arm-level events and the tumor immune environment.
We performed phylogenetic analysis of high-grade serous ovarian cancers (68 samples from seven patients), identifying constituent clones and quantifying their relative abundances at multiple ...intraperitoneal sites. Through whole-genome and single-nucleus sequencing, we identified evolutionary features including mutation loss, convergence of the structural genome and temporal activation of mutational processes that patterned clonal progression. We then determined the precise clonal mixtures comprising each tumor sample. The majority of sites were clonally pure or composed of clones from a single phylogenetic clade. However, each patient contained at least one site composed of polyphyletic clones. Five patients exhibited monoclonal and unidirectional seeding from the ovary to intraperitoneal sites, and two patients demonstrated polyclonal spread and reseeding. Our findings indicate that at least two distinct modes of intraperitoneal spread operate in clonal dissemination and highlight the distribution of migratory potential over clonal populations comprising high-grade serous ovarian cancers.
Full text
Available for:
IJS, NUK, SBMB, UL, UM, UPUK
Multiple myeloma (MM) remains an incurable disease, with a treatment-refractory state eventually developing in all patients. Constant clonal evolution and genetic heterogeneity of MM are a likely ...explanation for the emergence of drug-resistant disease. Monitoring of MM genomic evolution on therapy by serial bone marrow biopsy is unfortunately impractical because it involves an invasive and painful procedure. We describe how noninvasive and highly sensitive isolation and characterization of circulating tumor cells (CTCs) from peripheral blood at single-cell resolution recapitulate MM in the bone marrow. We demonstrate that CTCs provide the same genetic information as bone marrow MM cells and even reveal mutations with greater sensitivity than bone marrow biopsies in some cases. Single CTC RNA sequencing enables classification of MM and quantitative assessment of genes that are relevant for prognosis. We propose that the genomic characterization of CTCs should be included in clinical trials to follow the emergence of resistant subclones after MM therapy.
Loss of heterozygosity (LOH) and copy number alteration (CNA) feature prominently in the somatic genomic landscape of tumors. As such, karyotypic aberrations in cancer genomes have been studied ...extensively to discover novel oncogenes and tumor-suppressor genes. Advances in sequencing technology have enabled the cost-effective detection of tumor genome and transcriptome mutation events at single-base-pair resolution; however, computational methods for predicting segmental regions of LOH in this context are not yet fully explored. Consequently, whole transcriptome, nucleotide-level resolution analysis of monoallelic expression patterns associated with LOH has not yet been undertaken in cancer. We developed a novel approach for inference of LOH from paired tumor/normal sequence data and applied it to a cohort of 23 triple-negative breast cancer (TNBC) genomes. Following extensive benchmarking experiments, we describe the nucleotide-resolution landscape of LOH in TNBC and assess the consequent effect of LOH on the transcriptomes of these tumors using RNA-seq-derived measurements of allele-specific expression. We show that the majority of monoallelic expression in the transcriptomes of triple-negative breast cancer can be explained by genomic regions of LOH and establish an upper bound for monoallelic expression that may be explained by other tumor-specific modifications such as epigenetics or mutations. Monoallelically expressed genes associated with LOH reveal that cell cycle, homologous recombination and actin-cytoskeletal functions are putatively disrupted by LOH in TNBC. Finally, we show how inference of LOH can be used to interpret allele frequencies of somatic mutations and postulate on temporal ordering of mutations in the evolutionary history of these tumors.