Identifying DNase I hypersensitive sites as driver distal regulatory elements in breast cancer

Efforts to identify driver mutations in cancer have largely focused on genes, whereas non-coding sequences remain relatively unexplored. Here we develop a statistical method based on characteristics known to influence local mutation rate and a series of enrichment filters in order to identify distal regulatory elements harboring putative driver mutations in breast cancer. We identify ten DNase I hypersensitive sites that are significantly mutated in breast cancers and associated with the aberrant expression of neighboring genes. A pan-cancer analysis shows that three of these elements are significantly mutated across multiple cancer types and have mutation densities similar to protein-coding driver genes. Functional characterization of the most highly mutated DNase I hypersensitive sites in breast cancer (using in silico and experimental approaches) confirms that they are regulatory elements and affect the expression of cancer genes. Our study suggests that mutations of regulatory elements in tumors likely play an important role in cancer development.

Cancer driver mutations can occur within noncoding genomic sequences. Here, the authors develop a statistical approach to identify candidate noncoding driver mutations in DNase I hypersensitive sites in breast cancer and experimentally demonstrate they are regulatory elements of known cancer genes.

Whole Genome Sequencing Revealed Mutations in Two Independent Genes as the Underlying Cause of Retinal Degeneration in an Ashkenazi Jewish Pedigree

Retinitis pigmentosa (RP) causes progressive photoreceptor loss resulting from mutations in over 80 genes. This study identified the genetic cause of RP in three members of a non-consanguineous pedigree. Detailed ophthalmic evaluation was performed in the three affected family members. Whole exome sequencing (WES) and whole genome sequencing (WGS) were performed in the three affected and the two unaffected family members and variants were filtered to detect rare, potentially deleterious variants segregating with disease. WES and WGS did not identify potentially pathogenic variants shared by all three affected members. However, WES identified a previously reported homozygous nonsense mutation in KIZ (c.226C>T, p.Arg76*) in two affected sisters, but not in their affected second cousin. WGS revealed a novel 1.135 kb homozygous deletion in a retina transcript of C21orf2 and a novel 30.651 kb heterozygous deletion in CACNA2D4 in the affected second cousin. The sisters with the KIZ mutation carried no copies of the C21orf2 or CACNA2D4 deletions, while the second cousin with the C21orf2 and CACNA2D4 deletions carried no copies of the KIZ mutation. This study identified two independent, homozygous mutations in genes previously reported in autosomal recessive RP in a non-consanguineous family, and demonstrated the value of WGS when WES fails to identify likely disease-causing mutations.

Glioblastoma cellular cross-talk converges on NF-κB to attenuate EGFR inhibitor sensitivity

Zanca et al. show that heterogeneous expression of the wild-type EGFR receptor and its constitutively active mutant form, EGFRvIII, limits sensitivity to EGFR-directed therapies through an interclonal communication mechanism mediated by IL-6 cytokine secreted from EGFRvIII-positive tumor cells.

In glioblastoma (GBM), heterogeneous expression of amplified and mutated epidermal growth factor receptor (EGFR) presents a substantial challenge for the effective use of EGFR-directed therapeutics. Here we demonstrate that heterogeneous expression of the wild-type receptor and its constitutively active mutant form, EGFRvIII, limits sensitivity to these therapies through an interclonal communication mechanism mediated by interleukin-6 (IL-6) cytokine secreted from EGFRvIII-positive tumor cells. IL-6 activates a NF-κB signaling axis in a paracrine and autocrine manner, leading to bromodomain protein 4 (BRD4)-dependent expression of the prosurvival protein survivin (BIRC5) and attenuation of sensitivity to EGFR tyrosine kinase inhibitors (TKIs). NF-κB and survivin are coordinately up-regulated in GBM patient tumors, and functional inhibition of either protein or BRD4 in in vitro and in vivo models restores sensitivity to EGFR TKIs. These results provide a rationale for improving anti-EGFR therapeutic efficacy through pharmacological uncoupling of a convergence point of NF-κB-mediated survival that is leveraged by an interclonal circuitry mechanism established by intratumoral mutational heterogeneity.

iPSCORE: A Resource of 222 iPSC Lines Enabling Functional Characterization of Genetic Variation across a Variety of Cell Types

Large-scale collections of induced pluripotent stem cells (iPSCs) could serve as powerful model systems for examining how genetic variation affects biology and disease. Here we describe the iPSCORE resource: a collection of systematically derived and characterized iPSC lines from 222 ethnically diverse individuals that allows for both familial and association-based genetic studies. iPSCORE lines are pluripotent with high genomic integrity (no or low numbers of somatic copy-number variants) as determined using high-throughput RNA-sequencing and genotyping arrays, respectively. Using iPSCs from a family of individuals, we show that iPSC-derived cardiomyocytes demonstrate gene expression patterns that cluster by genetic background, and can be used to examine variants associated with physiological and disease phenotypes. The iPSCORE collection contains representative individuals for risk and non-risk alleles for 95% of SNPs associated with human phenotypes through genome-wide association studies. Our study demonstrates the utility of iPSCORE for examining how genetic variants influence molecular and physiological traits in iPSCs and derived cell lines.

High-Throughput and Cost-Effective Characterization of Induced Pluripotent Stem Cells

Reprogramming somatic cells to induced pluripotent stem cells (iPSCs) offers the possibility of studying the molecular mechanisms underlying human diseases in cell types difficult to extract from living patients, such as neurons and cardiomyocytes. To date, studies have been published that use small panels of iPSC-derived cell lines to study monogenic diseases. However, to study complex diseases, where the genetic variation underlying the disorder is unknown, a sizable number of patient-specific iPSC lines and controls need to be generated. Currently the methods for deriving and characterizing iPSCs are time consuming, expensive, and, in some cases, descriptive but not quantitative. Here we set out to develop a set of simple methods that reduce cost and increase throughput in the characterization of iPSC lines. Specifically, we outline methods for high-throughput quantification of surface markers, gene expression analysis of in vitro differentiation potential, and evaluation of karyotype with markedly reduced cost.

Large-Scale Profiling Reveals the Influence of Genetic Variation on Gene Expression in Human Induced Pluripotent Stem Cells

In this study, we used whole-genome sequencing and gene expression profiling of 215 human induced pluripotent stem cell (iPSC) lines from different donors to identify genetic variants associated with RNA expression for 5,746 genes. We were able to predict causal variants for these expression quantitative trait loci (eQTLs) that disrupt transcription factor binding and validated a subset of them experimentally. We also identified copy-number variant (CNV) eQTLs, including some that appear to affect gene expression by altering the copy number of intergenic regulatory regions. In addition, we were able to identify effects on gene expression of rare genic CNVs and regulatory single-nucleotide variants and found that reactivation of gene expression on the X chromosome depends on gene chromosomal position. Our work highlights the value of iPSCs for genetic association analyses and provides a unique resource for investigating the genetic regulation of gene expression in pluripotent cells.

Aberrant DNA Methylation in Human iPSCs Associates with MYC-Binding Motifs in a Clone-Specific Manner Independent of Genetics

Induced pluripotent stem cells (iPSCs) show variable methylation patterns between lines, some of which reflect aberrant differences relative to embryonic stem cells (ESCs). To examine whether this aberrant methylation results from genetic variation or non-genetic mechanisms, we generated human iPSCs from monozygotic twins to investigate how genetic background, clone, and passage number contribute. We found that aberrantly methylated CpGs are enriched in regulatory regions associated with MYC protein motifs and affect gene expression. We classified differentially methylated CpGs as being associated with genetic and/or non-genetic factors (clone and passage), and we found that aberrant methylation preferentially occurs at CpGs associated with clone-specific effects. We further found that clone-specific effects play a strong role in recurrent aberrant methylation at specific CpG sites across different studies. Our results argue that a non-genetic biological mechanism underlies aberrant methylation in iPSCs and that it is likely based on a probabilistic process involving MYC that takes place during or shortly after reprogramming.

PI3Kγ is a molecular switch that controls immune suppression

Macrophages play critical, but opposite, roles in acute and chronic inflammation and cancer. In response to pathogens or injury, inflammatory macrophages express cytokines that stimulate cytotoxic T cells, whereas macrophages in neoplastic and parasitic diseases express anti-inflammatory cytokines that induce immune suppression and may promote resistance to T cell checkpoint inhibitors. Here we show that macrophage PI 3-kinase γ controls a critical switch between immune stimulation and suppression during inflammation and cancer. PI3Kγ signalling through Akt and mTor inhibits NFκB activation while stimulating C/EBPβ activation, thereby inducing a transcriptional program that promotes immune suppression during inflammation and tumour growth. By contrast, selective inactivation of macrophage PI3Kγ stimulates and prolongs NFκB activation and inhibits C/EBPβ activation, thus promoting an immunostimulatory transcriptional program that restores CD8+ T cell activation and cytotoxicity. PI3Kγ synergizes with checkpoint inhibitor therapy to promote tumour regression and increased survival in mouse models of cancer. In addition, PI3Kγ-directed, anti-inflammatory gene expression can predict survival probability in cancer patients. Our work thus demonstrates that therapeutic targeting of intracellular signalling pathways that regulate the switch between macrophage polarization states can control immune suppression in cancer and other disorders.

Establishing the involvement of the novel gene AGBL5 in retinitis pigmentosa by whole genome sequencing

While more than 250 genes are known to cause inherited retinal degenerations (IRD), nearly 40-50% of families have the genetic basis for their disease unknown. In this study we sought to identify the underlying cause of IRD in a family by whole genome sequence (WGS) analysis. Clinical characterization including standard ophthalmic examination, fundus photography, visual field testing, electroretinography, and review of medical and family history was performed. WGS was performed on affected and unaffected family members using Illumina HiSeq X10. Sequence reads were aligned to hg19 using BWA-MEM and variant calling was performed with Genome Analysis Toolkit. The called variants were annotated with SnpEff v4.11, PolyPhen v2.2.2, and CADD v1.3. Copy number variations were called using Genome STRiP (svtoolkit 2.00.1611) and SpeedSeq software. Variants were filtered to detect rare potentially deleterious variants segregating with disease. Candidate variants were validated by dideoxy sequencing. Clinical evaluation revealed typical adolescent-onset recessive retinitis pigmentosa (arRP) in affected members. WGS identified about 4 million variants in each individual. Two rare and potentially deleterious compound heterozygous variants p.Arg281Cys and p.Arg487* were identified in the gene ATP/GTP binding protein like 5 (AGBL5) as likely causal variants. No additional variants in IRD genes that segregated with disease were identified. Mutation analysis confirmed the segregation of these variants with the IRD in the pedigree. Homology models indicated destabilization of AGBL5 due to the p.Arg281Cys change. Our findings establish the involvement of mutations in AGBL5 in RP and validate the WGS variant filtering pipeline we designed.

ADAR1 Activation Drives Leukemia Stem Cell Self-Renewal by Impairing Let-7 Biogenesis

Post-transcriptional adenosine-to-inosine RNA editing mediated by adenosine deaminase acting on RNA1 (ADAR1) promotes cancer progression and therapeutic resistance. However, ADAR1 editase-dependent mechanisms governing leukemia stem cell (LSC) generation have not been elucidated. In blast crisis chronic myeloid leukemia (BC CML), we show that increased JAK2 signaling and BCR-ABL1 amplification activate ADAR1. In a humanized BC CML mouse model, combined JAK2 and BCR-ABL1 inhibition prevents LSC self-renewal commensurate with ADAR1 downregulation. Lentiviral ADAR1 wild-type, but not an editing-defective ADAR1(E912A) mutant, induces self-renewal gene expression and impairs biogenesis of stem cell regulatory let-7 microRNAs. Combined RNA sequencing, qRT-PCR, CLIP-ADAR1, and pri-let-7 mutagenesis data suggest that ADAR1 promotes LSC generation via let-7 pri-microRNA editing and LIN28B upregulation. A small-molecule tool compound antagonizes ADAR1's effect on LSC self-renewal in stromal co-cultures and restores let-7 biogenesis. Thus, ADAR1 activation represents a unique therapeutic vulnerability in LSCs with active JAK2 signaling.

Associations Between Common and Rare Exonic Genetic Variants and Serum Levels of 20 Cardiovascular-Related Proteins: The Tromsø Study

Supplemental Digital Content is available in the text.

Background—

Genetic variation can be used to study causal relationships between biomarkers and diseases. Here, we identify new common and rare genetic variants associated with cardiovascular-related protein levels (protein quantitative trait loci [pQTLs]). We functionally annotate these pQTLs, predict and experimentally confirm a novel molecular interaction, and determine which pQTLs are associated with diseases and physiological phenotypes.

Methods and Results—

As part of a larger case–control study of venous thromboembolism, serum levels of 51 proteins implicated in cardiovascular diseases were measured in 330 individuals from the Tromsø Study. Exonic genetic variation near each protein’s respective gene (cis) was identified using sequencing and arrays. Using single site and gene-based tests, we identified 27 genetic associations between pQTLs and the serum levels of 20 proteins: 14 associated with common variation in cis, of which 6 are novel (ie, not previously reported); 7 associations with rare variants in cis, of which 4 are novel; and 6 associations in trans. Of the 20 proteins, 15 were associated with single sites and 7 with rare variants. cis-pQTLs for kallikrein and F12 also show trans associations for proteins (uPAR, kininogen) known to be cleaved by kallikrein and with NTproBNP. We experimentally demonstrate that kallikrein can cleave proBNP (NTproBNP precursor) in vitro. Nine of the pQTLs have previously identified associations with 17 disease and physiological phenotypes.

Conclusions—

We have identified cis and trans genetic variation associated with the serum levels of 20 proteins and utilized these pQTLs to study molecular mechanisms underlying disease and physiological phenotypes.

Joint effects of cancer and variants in the factor 5 gene on the risk of venous thromboembolism

Venous thromboembolism occurs frequently in cancer patients. Two variants in the factor 5 gene (F5), rs6025 encoding for the factor V Leiden mutation R506Q, and rs4524 encoding K858R, have been found to be associated with venous thromboembolism. We assessed the joint effect of active cancer and these two F5 variants on venous thromboembolism risk in a case-cohort study. Cases with a first venous thromboembolism (n=609) and a randomly selected age-weighted cohort (n=1,691) were sampled from the general population in Tromsø, Norway. Venous thromboembolism was classified as cancer-related if it occurred in the period 6 months before to 2 years after a diagnosis of cancer. Active cancer was associated with an 8.9-fold higher risk of venous thromboembolism (95% CI 7.2-10.9). The risk of cancer-related venous thromboembolism was 16.7-fold (95% CI 9.9-28.0) higher in subjects heterozygous for rs6025 compared with non-carriers of this variant without active cancer. In subjects with active cancer the risk of venous thromboembolism was 15.9-fold higher (95% CI 9.1-27.9) in those with one risk allele at rs4524, and 21.1-fold (95% CI 12.4-35.8) higher in those with two risk alleles compared with non-carriers without active cancer. A synergistic interaction was observed between active cancer and factor V Leiden (relative excess risk due to interaction 7.0; 95% CI 0.5-14.4) and rs4524 (relative excess risk due to interaction 15.0; 95% CI 7.5-29.2). The incidence of venous thromboembolism during the initial 6 months following a diagnosis of cancer was particularly high in subjects with risk alleles at these loci. This implies that the combination of cancer and F5 variants synergistically increases venous thromboembolism risk.

Mutational Profiling Can Establish Clonal or Independent Origin in Synchronous Bilateral Breast and Other Tumors

Background

Synchronous tumors can be independent primary tumors or a primary-metastatic (clonal) pair, which may have clinical implications. Mutational profiling of tumor DNA is increasingly common in the clinic. We investigated whether mutational profiling can distinguish independent from clonal tumors in breast and other cancers, using a carefully defined test based on the Clonal Likelihood Score (CLS = 100 x # shared high confidence (HC) mutations/ # total HC mutations).

Methods

Statistical properties of a formal test using the CLS were investigated. A high CLS is evidence in favor of clonality; the test is implemented as a one-sided binomial test of proportions. Test parameters were empirically determined using 16,422 independent breast tumor pairs and 15 primary-metastatic tumor pairs from 10 cancer types using The Cancer Genome Atlas.

Results

We validated performance of the test with its established parameters, using five published data sets comprising 15,758 known independent tumor pairs (maximum CLS = 4.1%, minimum p-value = 0.48) and 283 known tumor clonal pairs (minimum CLS 13%, maximum p-value <0.01), across renal cell, testicular, and colorectal cancer. The CLS test correctly classified all validation samples but one, which it appears may have been incorrectly classified in the published data. As proof-of-concept we then applied the CLS test to two new cases of invasive synchronous bilateral breast cancer at our institution, each with one hormone receptor positive (ER+/PR+/HER2-) lobular and one triple negative ductal carcinoma. High confidence mutations were identified by exome sequencing and results were validated using deep targeted sequencing. The first tumor pair had CLS of 81% (p-value < 10–15), supporting clonality. In the second pair, no common mutations of 184 variants were validated (p-value >0.99), supporting independence. A plausible molecular mechanism for the shift from hormone receptor positive to triple negative was identified in the clonal pair.

Conclusion

We have developed the statistical properties of a carefully defined Clonal Likelihood Score test from mutational profiling of tumor DNA. Under identified conditions, the test appears to reliably distinguish between synchronous tumors of clonal and of independent origin in several cancer types. This approach may have scientific and clinical utility.

Abstract 4803: BAP1 mutation is a frequent somatic event in peritoneal malignant mesothelioma

Malignant mesothelioma (MM) arises from mesothelial cells that line the pleural, peritoneal and pericardial surfaces. The majority of MMs are pleural and have been associated with asbestos exposure. Previously, pleural MMs have been genetically characterized by the loss of BAP1 (40-60%) as well as loss of NF2 (75%) and CDKN2A (60%). The rare peritoneal form of MM occurs in ∼10% cases. With only ∼300 cases diagnosed in the US per year, its link to asbestos exposure is not clear and its mutational landscape unknown. We analyzed the somatic mutational landscape of 12 peritoneal MM of epitheloid subtype using copy number analysis (N = 9), whole exome sequencing (N = 7) and targeted sequencing (N = 12). Peritoneal MM display few copy number alterations, with most samples having less than 10 Mbp total changes, mostly through deletions and no high copy number amplification. Chromosome band 3p21 encoding BAP1 is the most recurrently deleted region (5/9), while, in contrast to pleural MM, NF2 and CDKN2A are not affected. We further identified 87 non-silent mutations across 7 sequenced tumors, with a median of 8 mutated genes per tumor, resulting in a very low mutation rate (median 1.3 10-6). BAP1 was the only recurrently mutated gene (N = 3/7). In one additional case, loss of the entire chromosome 3 leaves a non-functional copy of BAP1 carrying a rare nonsense germline variant, thus suggesting a potential genetic predisposition in this patient. Finally, with targeted sequencing of BAP1 in 3 additional cases, we conclude that BAP1 is frequently altered through copy number losses (N = 3/12), mutations (N = 3/12) or both (N = 2/12) sometimes at a sub-clonal level. Our findings suggest a major role for BAP1 in peritoneal MM susceptibility and oncogenesis and indicate important molecular differences to pleural MM as well as potential strategies for therapy and prevention.

Citation Format: Olivier Harismendy, Hakan Alakus, Shawn Yost, Brian Woo, Randall French, Grace Y. Lin, Kristen Jepsen, Kelly A. Frazer, Andrew M. Lowy. BAP1 mutation is a frequent somatic event in peritoneal malignant mesothelioma. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4803. doi:10.1158/1538-7445.AM2015-4803

Brief Report: Oxidative Stress Mediates Cardiomyocyte Apoptosis in a Human Model of Danon Disease and Heart Failure

Danon disease is a familial cardiomyopathy associated with impaired autophagy due to mutations in the gene encoding lysosomal-associated membrane protein type 2 (LAMP-2). Emerging evidence has highlighted the importance of autophagy in regulating cardiomyocyte bioenergetics, function, and survival. However, the mechanisms responsible for cellular dysfunction and death in cardiomyocytes with impaired autophagic flux remain unclear. To investigate the molecular mechanisms responsible for Danon disease, we created induced pluripotent stem cells (iPSCs) from two patients with different LAMP-2 mutations. Danon iPSC-derived cardiomyocytes (iPSC-CMs) exhibited impaired autophagic flux and key features of heart failure such as increased cell size, increased expression of natriuretic peptides, and abnormal calcium handling compared to control iPSC-CMs. Additionally, Danon iPSC-CMs demonstrated excessive amounts of mitochondrial oxidative stress and apoptosis. Using the sulfhydryl antioxidant N-acetylcysteine to scavenge free radicals resulted in a significant reduction in apoptotic cell death in Danon iPSC-CMs. In summary, we have modeled Danon disease using human iPSC-CMs from patients with mutations in LAMP-2, allowing us to gain mechanistic insight into the pathogenesis of this disease. We demonstrate that LAMP-2 deficiency leads to an impairment in autophagic flux, which results in excessive oxidative stress, and subsequent cardiomyocyte apoptosis. Scavenging excessive free radicals with antioxidants may be beneficial for patients with Danon disease. In vivo studies will be necessary to validate this new treatment strategy.

BAP1 mutation is a frequent somatic event in peritoneal malignant mesothelioma

BACKGROUND

Malignant mesothelioma (MM) arises from mesothelial cells that line the pleural, peritoneal and pericardial surfaces. The majority of MMs are pleural and have been associated with asbestos exposure. Previously, pleural MMs have been genetically characterized by the loss of BAP1 (40-60%) as well as loss of NF2 (75%) and CDKN2A (60%). The rare peritoneal form of MM occurs in ~10% cases. With only ~300 cases diagnosed in the US per year, its link to asbestos exposure is not clear and its mutational landscape unknown.

METHODS

We analyzed the somatic mutational landscape of 12 peritoneal MM of epitheloid subtype using copy number analysis (N = 9), whole exome sequencing (N = 7) and targeted sequencing (N = 12).

RESULTS

Peritoneal MM display few copy number alterations, with most samples having less than 10 Mbp total changes, mostly through deletions and no high copy number amplification. Chromosome band 3p21 encoding BAP1 is the most recurrently deleted region (5/9), while, in contrast to pleural MM, NF2 and CDKN2A are not affected. We further identified 87 non-silent mutations across 7 sequenced tumors, with a median of 8 mutated genes per tumor, resulting in a very low mutation rate (median 1.3 10(-6)). BAP1 was the only recurrently mutated gene (N = 3/7). In one additional case, loss of the entire chromosome 3 leaves a non-functional copy of BAP1 carrying a rare nonsense germline variant, thus suggesting a potential genetic predisposition in this patient. Finally, with targeted sequencing of BAP1 in 3 additional cases, we conclude that BAP1 is frequently altered through copy number losses (N = 3/12), mutations (N = 3/12) or both (N = 2/12) sometimes at a sub-clonal level.

CONCLUSION

Our findings suggest a major role for BAP1 in peritoneal MM susceptibility and oncogenesis and indicate important molecular differences to pleural MM as well as potential strategies for therapy and prevention.

Genetic and epigenetic profiling of CLL disease progression reveals limited somatic evolution and suggests a relationship to memory-cell development

We examined genetic and epigenetic changes that occur during disease progression from indolent to aggressive forms of chronic lymphocytic leukemia (CLL) using serial samples from 27 patients. Analysis of DNA mutations grouped the leukemia cases into three categories: evolving (26%), expanding (26%) and static (47%). Thus, approximately three-quarters of the CLL cases had little to no genetic subclonal evolution. However, we identified significant recurrent DNA methylation changes during progression at 4752 CpGs enriched for regions near Polycomb 2 repressive complex (PRC2) targets. Progression-associated CpGs near the PRC2 targets undergo methylation changes in the same direction during disease progression as during normal development from naive to memory B cells. Our study shows that CLL progression does not typically occur via subclonal evolution, but that certain CpG sites undergo recurrent methylation changes. Our results suggest CLL progression may involve developmental processes shared in common with the generation of normal memory B cells.

Transcriptome sequencing reveals potential mechanism of cryptic 3' splice site selection in SF3B1-mutated cancers

Mutations in the splicing factor SF3B1 are found in several cancer types and have been associated with various splicing defects. Using transcriptome sequencing data from chronic lymphocytic leukemia, breast cancer and uveal melanoma tumor samples, we show that hundreds of cryptic 3’ splice sites (3’SSs) are used in cancers with SF3B1 mutations. We define the necessary sequence context for the observed cryptic 3’ SSs and propose that cryptic 3’SS selection is a result of SF3B1 mutations causing a shift in the sterically protected region downstream of the branch point. While most cryptic 3’SSs are present at low frequency (<10%) relative to nearby canonical 3’SSs, we identified ten genes that preferred out-of-frame cryptic 3’SSs. We show that cancers with mutations in the SF3B1 HEAT 5-9 repeats use cryptic 3’SSs downstream of the branch point and provide both a mechanistic model consistent with published experimental data and affected targets that will guide further research into the oncogenic effects of SF3B1 mutation.

A key goal of cancer genomics studies is to identify genes that are recurrently mutated at a rate above background and likely contribute to cancer development. Many such recurrently mutated genes have been identified over the last few years, but we often do not know the underlying mechanisms by which they contribute to cancer growth. Unexpectedly, several genes in the spliceosome, the collection of RNAs and proteins that remove introns from transcribed RNAs, are recurrently mutated in different cancers. Here, we have examined mutations in the splicing factor SF3B1, a key component of the spliceosome, and identified a global splicing defect present in different cancers with SF3B1 mutations by comparing the expression of splice junctions using generalized linear models. While prior studies have reported a limited number of aberrant splicing events in SF3B1-mutated cancers, we have established that SF3B1 mutations are associated with usage of hundreds of atypical splice sites at the 3’ end of the intron. We have identified nucleotide sequence requirements for these cryptic splice sites that are consistent with a proposed mechanistic model. These findings greatly expand our understanding of the effect of SF3B1 mutations on splicing and provide new targets for determining the oncogenic effect of SF3B1 mutations.

An RNA editing fingerprint of cancer stem cell reprogramming

BACKGROUND

Deregulation of RNA editing by adenosine deaminases acting on dsRNA (ADARs) has been implicated in the progression of diverse human cancers including hematopoietic malignancies such as chronic myeloid leukemia (CML). Inflammation-associated activation of ADAR1 occurs in leukemia stem cells specifically in the advanced, often drug-resistant stage of CML known as blast crisis. However, detection of cancer stem cell-associated RNA editing by RNA sequencing in these rare cell populations can be technically challenging, costly and requires PCR validation. The objectives of this study were to validate RNA editing of a subset of cancer stem cell-associated transcripts, and to develop a quantitative RNA editing fingerprint assay for rapid detection of aberrant RNA editing in human malignancies.

METHODS

To facilitate quantification of cancer stem cell-associated RNA editing in exons and intronic or 3'UTR primate-specific Alu sequences using a sensitive, cost-effective method, we established an in vitro RNA editing model and developed a sensitive RNA editing fingerprint assay that employs a site-specific quantitative PCR (RESSq-PCR) strategy. This assay was validated in a stably-transduced human leukemia cell line, lentiviral-ADAR1 transduced primary hematopoietic stem and progenitor cells, and in primary human chronic myeloid leukemia stem cells.

RESULTS

In lentiviral ADAR1-expressing cells, increased RNA editing of MDM2, APOBEC3D, GLI1 and AZIN1 transcripts was detected by RESSq-PCR with improved sensitivity over sequencing chromatogram analysis. This method accurately detected cancer stem cell-associated RNA editing in primary chronic myeloid leukemia samples, establishing a cancer stem cell-specific RNA editing fingerprint of leukemic transformation that will support clinical development of novel diagnostic tools to predict and prevent cancer progression.

CONCLUSIONS

RNA editing quantification enables rapid detection of malignant progenitors signifying cancer progression and therapeutic resistance, and will aid future RNA editing inhibitor development efforts.

Abstract 4798: The niche specific role of CD44 splice isoform expression in blast crisis leukemia stem cell generation

Introduction

Chronic myeloid leukemia (CML) is a myeloproliferative neoplasm initiated in hematopoietic stem cells by expression of the BCR-ABL1 fusion oncogene and its protein product, which enhances ABL1 kinase activity. Targeted BCR-ABL tyrosine kinase inhibitors (TKIs), such as Imatinib, and the second generation TKIs like Nilotinib as well as the dual specific SCR and ABL inhibitor Dasatinib have significantly slowed disease progression by eradicating the bulk of BCR-ABL1 expressing cells in the circulation. However, progression to a therapeutically resistant blast crisis (BC) phase is driven by CD34+CD38+Lin- progenitors that co-opt stem cell properties, such as enhanced self-renewal and survival, albeit in a deregulated manner. These BC leukemia stem cells (BC LSC) harbor enhanced BCL2 and beta-catenin expression. Seminal research shows that a key Wnt/beta-catenin target gene, CD44, plays a vital role in cancer stem cell survival and retention in the malignant microenvironment. Because alternative mRNA splicing in humans generates many CD44 isoforms, we investigated CD44 variant expression and retention of human BC LSC in specific hematopoietic niches.

Methods and Results

We performed whole transcriptome RNA sequencing (RNA Seq) of FACS sorted progenitors (Lin-CD34+CD38+) from chronic phase (CP)(n=8) and blast crisis (BC)(n=8) CML patients as well as from normal cord blood (CB) (n=7) and adult peripheral blood (NPB)(n=4). A number of CD44 transcript variants were detected: 3, 4 (CD44s), 5, 6, 7, 8 plus additional variants. RNA seq analysis uncovered a higher overall CD44 expression in BC compared to CP progenitors. Notably, one isoform of CD44, variant 3, was highly upregulated in BC compared to CP progenitors and minimally expressed by their normal counterparts. Moreover, BC progenitors that engrafted in RAG2-/-γc-/- harbored CD44v3 expression in the bone marrow and the splenic niche and this differential BC LSC isoform expression was reduced after dasatinib treatment. In addition, BC LSC in the spleen showed a reduction of migration specific markers, such as RHAMM, ICAM-1 and Osteopontin, compared with the bone marrow resident BC LSC. Furthermore, a comparative splice isoform specific qRT-PCR analysis of CD44 variants expression in young versus old bone marrow showed no correlation with aging. In fact, human embryonic stem cells harbored variant 3 expression.

Conclusions

These data suggest that CD44 transcript variant 3 upregulation occurs following malignant reprogramming of hematopoietic progenitors that enables them to adopt features of an embryonic transcriptional program in select microenvironments. Detection of CD44 variant 3 has the potential to more precisely identify patients at risk for relapse and progression, suggest that combined therapeutic strategies involving a BCR-ABL specific TKI and a CD44 monoclonal antibody may decrease relapse and BC transformation rates.

Citation Format: Frida L. Holm, Eva Hellqvist, Cayla Mason, Christian Barrett, Kelly A. Frazer, Anil Sadarangani, Catriona HM Jamieson. The niche specific role of CD44 splice isoform expression in blast crisis leukemia stem cell generation. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 4798. doi:10.1158/1538-7445.AM2014-4798

Biased estimates of clonal evolution and subclonal heterogeneity can arise from PCR duplicates in deep sequencing experiments

Accurate allele frequencies are important for measuring subclonal heterogeneity and clonal evolution. Deep-targeted sequencing data can contain PCR duplicates, inflating perceived read depth. Here we adapted the Illumina TruSeq Custom Amplicon kit to include single molecule tagging (SMT) and show that SMT-identified duplicates arise from PCR. We demonstrate that retention of PCR duplicate reads can imply clonal evolution when none exists, while their removal effectively controls the false positive rate. Additionally, PCR duplicates alter estimates of subclonal heterogeneity in tumor samples. Our method simplifies PCR duplicate identification and emphasizes their removal in studies of tumor heterogeneity and clonal evolution.

Network-based analysis identifies epigenetic biomarkers of esophageal squamous cell carcinoma progression

MOTIVATION

A rapid progression of esophageal squamous cell carcinoma (ESCC) causes a high mortality rate because of the propensity for metastasis driven by genetic and epigenetic alterations. The identification of prognostic biomarkers would help prevent or control metastatic progression. Expression analyses have been used to find such markers, but do not always validate in separate cohorts. Epigenetic marks, such as DNA methylation, are a potential source of more reliable and stable biomarkers. Importantly, the integration of both expression and epigenetic alterations is more likely to identify relevant biomarkers.

RESULTS

We present a new analysis framework, using ESCC progression-associated gene regulatory network (GRN escc), to identify differentially methylated CpG sites prognostic of ESCC progression. From the CpG loci differentially methylated in 50 tumor-normal pairs, we selected 44 CpG loci most highly associated with survival and located in the promoters of genes more likely to belong to GRN escc. Using an independent ESCC cohort, we confirmed that 8/10 of CpG loci in the promoter of GRN escc genes significantly correlated with patient survival. In contrast, 0/10 CpG loci in the promoter genes outside the GRN escc were correlated with patient survival. We further characterized the GRN escc network topology and observed that the genes with methylated CpG loci associated with survival deviated from the center of mass and were less likely to be hubs in the GRN escc. We postulate that our analysis framework improves the identification of bona fide prognostic biomarkers from DNA methylation studies, especially with partial genome coverage.