3D genomic analysis reveals novel enhancer-hijacking caused by complex structural alterations that drive oncogene overexpression.

Cancer genomes are composed of many complex structural alterations on chromosomes and extrachromosomal DNA (ecDNA), making it difficult to identify non-coding enhancer regions that are hijacked to activate oncogene expression. Here, we describe a 3D genomics-based analysis called HAPI (Highly Active Promoter Interactions) to characterize enhancer hijacking. HAPI analysis of HiChIP data from 34 cancer cell lines identified enhancer hijacking events that activate both known and potentially novel oncogenes such as MYC, CCND1 , ETV1 , CRKL , and ID4 . Furthermore, we found enhancer hijacking among multiple oncogenes from different chromosomes, often including MYC , on the same complex amplicons such as ecDNA. We characterized a MYC - ERBB2 chimeric ecDNA, in which ERBB2 heavily hijacks MYC 's enhancers. Notably, CRISPRi of the MYC promoter led to increased interaction of ERBB2 with MYC enhancers and elevated ERBB2 expression. Our HAPI analysis tool provides a robust strategy to detect enhancer hijacking and reveals novel insights into oncogene activation.

TXNIP loss expands Myc-dependent transcriptional programs by increasing Myc genomic binding.

The c-Myc protooncogene places a demand on glucose uptake to drive glucose-dependent biosynthetic pathways. To meet this demand, c-Myc protein (Myc henceforth) drives the expression of glucose transporters, glycolytic enzymes, and represses the expression of thioredoxin interacting protein (TXNIP), which is a potent negative regulator of glucose uptake. A Mychigh/TXNIPlow gene signature is clinically significant as it correlates with poor clinical prognosis in triple-negative breast cancer (TNBC) but not in other subtypes of breast cancer, suggesting a functional relationship between Myc and TXNIP. To better understand how TXNIP contributes to the aggressive behavior of TNBC, we generated TXNIP null MDA-MB-231 (231:TKO) cells for our study. We show that TXNIP loss drives a transcriptional program that resembles those driven by Myc and increases global Myc genome occupancy. TXNIP loss allows Myc to invade the promoters and enhancers of target genes that are potentially relevant to cell transformation. Together, these findings suggest that TXNIP is a broad repressor of Myc genomic binding. The increase in Myc genomic binding in the 231:TKO cells expands the Myc-dependent transcriptome we identified in parental MDA-MB-231 cells. This expansion of Myc-dependent transcription following TXNIP loss occurs without an apparent increase in Myc's intrinsic capacity to activate transcription and without increasing Myc levels. Together, our findings suggest that TXNIP loss mimics Myc overexpression, connecting Myc genomic binding and transcriptional programs to the nutrient and progrowth signals that control TXNIP expression.

Loss of G0/G1 switch gene 2 (G0S2) promotes disease progression and drug resistance in chronic myeloid leukaemia (CML) by disrupting glycerophospholipid metabolism.

Tyrosine kinase inhibitors (TKIs) targeting BCR::ABL1 have turned chronic myeloid leukaemia (CML) from a fatal disease into a manageable condition for most patients. Despite improved survival, targeting drug-resistant leukaemia stem cells (LSCs) remains a challenge for curative CML therapy. Aberrant lipid metabolism can have a large impact on membrane dynamics, cell survival and therapeutic responses in cancer. While ceramide and sphingolipid levels were previously correlated with TKI response in CML, the role of lipid metabolism in TKI resistance is not well understood. We have identified downregulation of a critical regulator of lipid metabolism, G0/G1 switch gene 2 (G0S2), in multiple scenarios of TKI resistance, including (1) BCR::ABL1 kinase-independent TKI resistance, (2) progression of CML from the chronic to the blast phase of the disease, and (3) in CML versus normal myeloid progenitors. Accordingly, CML patients with low G0S2 expression levels had a worse overall survival. G0S2 downregulation in CML was not a result of promoter hypermethylation or BCR::ABL1 kinase activity, but was rather due to transcriptional repression by MYC. Using CML cell lines, patient samples and G0s2 knockout (G0s2-/- ) mice, we demonstrate a tumour suppressor role for G0S2 in CML and TKI resistance. Our data suggest that reduced G0S2 protein expression in CML disrupts glycerophospholipid metabolism, correlating with a block of differentiation that renders CML cells resistant to therapy. Altogether, our data unravel a new role for G0S2 in regulating myeloid differentiation and TKI response in CML, and suggest that restoring G0S2 may have clinical utility.

Comprehensive evaluation and efficient classification of BRCA1 RING domain missense substitutions.

BRCA1 is a high-risk susceptibility gene for breast and ovarian cancer. Pathogenic protein-truncating variants are scattered across the open reading frame, but all known missense substitutions that are pathogenic because of missense dysfunction are located in either the amino-terminal RING domain or the carboxy-terminal BRCT domain. Heterodimerization of the BRCA1 and BARD1 RING domains is a molecularly defined obligate activity. Hence, we tested every BRCA1 RING domain missense substitution that can be created by a single nucleotide change for heterodimerization with BARD1 in a mammalian two-hybrid assay. Downstream of the laboratory assay, we addressed three additional challenges: assay calibration, validation thereof, and integration of the calibrated results with other available data, such as computational evidence and patient/population observational data to achieve clinically applicable classification. Overall, we found that 15%-20% of BRCA1 RING domain missense substitutions are pathogenic. Using a Bayesian point system for data integration and variant classification, we achieved clinical classification of 89% of observed missense substitutions. Moreover, among missense substitutions not present in the human observational data used here, we find an additional 45 with concordant computational and functional assay evidence in favor of pathogenicity plus 223 with concordant evidence in favor of benignity; these are particularly likely to be classified as likely pathogenic and likely benign, respectively, once human observational data become available.

A human breast cancer-derived xenograft and organoid platform for drug discovery and precision oncology.

Models that recapitulate the complexity of human tumors are urgently needed to develop more effective cancer therapies. We report a bank of human patient-derived xenografts (PDXs) and matched organoid cultures from tumors that represent the greatest unmet need: endocrine-resistant, treatment-refractory and metastatic breast cancers. We leverage matched PDXs and PDX-derived organoids (PDxO) for drug screening that is feasible and cost-effective with in vivo validation. Moreover, we demonstrate the feasibility of using these models for precision oncology in real time with clinical care in a case of triple-negative breast cancer (TNBC) with early metastatic recurrence. Our results uncovered a Food and Drug Administration (FDA)-approved drug with high efficacy against the models. Treatment with this therapy resulted in a complete response for the individual and a progression-free survival (PFS) period more than three times longer than their previous therapies. This work provides valuable methods and resources for functional precision medicine and drug development for human breast cancer.

MS4A3 promotes differentiation in chronic myeloid leukemia by enhancing common ß-chain cytokine receptor endocytosis.

The chronic phase of chronic myeloid leukemia (CP-CML) is characterized by the excessive production of maturating myeloid cells. As CML stem/progenitor cells (LSPCs) are poised to cycle and differentiate, LSPCs must balance conservation and differentiation to avoid exhaustion, similar to normal hematopoiesis under stress. Since BCR-ABL1 tyrosine kinase inhibitors (TKIs) eliminate differentiating cells but spare BCR-ABL1-independent LSPCs, understanding the mechanisms that regulate LSPC differentiation may inform strategies to eliminate LSPCs. Upon performing a meta-analysis of published CML transcriptomes, we discovered that low expression of the MS4A3 transmembrane protein is a universal characteristic of LSPC quiescence, BCR-ABL1 independence, and transformation to blast phase (BP). Several mechanisms are involved in suppressing MS4A3, including aberrant methylation and a MECOM-C/EBPε axis. Contrary to previous reports, we find that MS4A3 does not function as a G1/S phase inhibitor but promotes endocytosis of common ß-chain (ßc) cytokine receptors upon GM-CSF/IL-3 stimulation, enhancing downstream signaling and cellular differentiation. This suggests that LSPCs downregulate MS4A3 to evade ßc cytokine-induced differentiation and maintain a more primitive, TKI-insensitive state. Accordingly, knockdown (KD) or deletion of MS4A3/Ms4a3 promotes TKI resistance and survival of CML cells ex vivo and enhances leukemogenesis in vivo, while targeted delivery of exogenous MS4A3 protein promotes differentiation. These data support a model in which MS4A3 governs response to differentiating myeloid cytokines, providing a unifying mechanism for the differentiation block characteristic of CML quiescence and BP-CML. Promoting MS4A3 reexpression or delivery of ectopic MS4A3 may help eliminate LSPCs in vivo.

A predominant enhancer co-amplified with the SOX2 oncogene is necessary and sufficient for its expression in squamous cancer.

Amplification and overexpression of the SOX2 oncogene represent a hallmark of squamous cancers originating from diverse tissue types. Here, we find that squamous cancers selectively amplify a 3' noncoding region together with SOX2, which harbors squamous cancer-specific chromatin accessible regions. We identify a single enhancer e1 that predominantly drives SOX2 expression. Repression of e1 in SOX2-high cells causes collapse of the surrounding enhancers, remarkable reduction in SOX2 expression, and a global transcriptional change reminiscent of SOX2 knockout. The e1 enhancer is driven by a combination of transcription factors including SOX2 itself and the AP-1 complex, which facilitates recruitment of the co-activator BRD4. CRISPR-mediated activation of e1 in SOX2-low cells is sufficient to rebuild the e1-SOX2 loop and activate SOX2 expression. Our study shows that squamous cancers selectively amplify a predominant enhancer to drive SOX2 overexpression, uncovering functional links among enhancer activation, chromatin looping, and lineage-specific copy number amplifications of oncogenes.

CaBagE: A Cas9-based Background Elimination strategy for targeted, long-read DNA sequencing.

A substantial fraction of the human genome is difficult to interrogate with short-read DNA sequencing technologies due to paralogy, complex haplotype structures, or tandem repeats. Long-read sequencing technologies, such as Oxford Nanopore's MinION, enable direct measurement of complex loci without introducing many of the biases inherent to short-read methods, though they suffer from relatively lower throughput. This limitation has motivated recent efforts to develop amplification-free strategies to target and enrich loci of interest for subsequent sequencing with long reads. Here, we present CaBagE, a method for target enrichment that is efficient and useful for sequencing large, structurally complex targets. The CaBagE method leverages the stable binding of Cas9 to its DNA target to protect desired fragments from digestion with exonuclease. Enriched DNA fragments are then sequenced with Oxford Nanopore's MinION long-read sequencing technology. Enrichment with CaBagE resulted in a median of 116X coverage (range 39-416) of target loci when tested on five genomic targets ranging from 4-20kb in length using healthy donor DNA. Four cancer gene targets were enriched in a single reaction and multiplexed on a single MinION flow cell. We further demonstrate the utility of CaBagE in two ALS patients with C9orf72 short tandem repeat expansions to produce genotype estimates commensurate with genotypes derived from repeat-primed PCR for each individual. With CaBagE there is a physical enrichment of on-target DNA in a given sample prior to sequencing. This feature allows adaptability across sequencing platforms and potential use as an enrichment strategy for applications beyond sequencing. CaBagE is a rapid enrichment method that can illuminate regions of the 'hidden genome' underlying human disease.

The lingering mysteries of metastatic recurrence in breast cancer.

Despite being the hallmark of cancer that is responsible for the highest number of deaths, very little is known about the biology of metastasis. Metastatic disease typically manifests after a protracted period of undetectable disease following surgery or systemic therapy, owing to relapse or recurrence. In the case of breast cancer, metastatic relapse can occur months to decades after initial diagnosis and treatment. In this review, we provide an overview of the known key factors that influence metastatic recurrence, with the goal of highlighting the critical unanswered questions that still need to be addressed to make a difference in the mortality of breast cancer patients.

STAT3 and GR Cooperate to Drive Gene Expression and Growth of Basal-Like Triple-Negative Breast Cancer.

Breast cancers are divided into subtypes with different prognoses and treatment responses based on global differences in gene expression. Luminal breast cancer gene expression and proliferation are driven by estrogen receptor alpha, and targeting this transcription factor is the most effective therapy for this subtype. By contrast, it remains unclear which transcription factors drive the gene expression signature that defines basal-like triple-negative breast cancer, and there are no targeted therapies approved to treat this aggressive subtype. In this study, we utilized integrated genomic analysis of DNA methylation, chromatin accessibility, transcription factor binding, and gene expression in large collections of breast cancer cell lines and patient tumors to identify transcription factors responsible for the basal-like gene expression program. Glucocorticoid receptor (GR) and STAT3 bind to the same genomic regulatory regions, which were specifically open and unmethylated in basal-like breast cancer. These transcription factors cooperated to regulate expression of hundreds of genes in the basal-like gene expression signature, which were associated with poor prognosis. Combination treatment with small-molecule inhibitors of both transcription factors resulted in synergistic decreases in cell growth in cell lines and patient-derived organoid models. This study demonstrates that GR and STAT3 cooperate to regulate the basal-like breast cancer gene expression program and provides the basis for improved therapy for basal-like triple-negative breast cancer through rational combination of STAT3 and GR inhibitors. SIGNIFICANCE: This study demonstrates that GR and STAT3 cooperate to activate the canonical gene expression signature of basal-like triple-negative breast cancer and that combination treatment with STAT3 and GR inhibitors could provide synergistic therapeutic efficacy.

Spatially-resolved quantification of proteins in triple negative breast cancers reveals differences in the immune microenvironment associated with prognosis.

Triple negative breast cancer (TNBC) is an aggressive breast cancer subtype. Recent studies have shown that MHC class II (MHCII) expression and tumor infiltrating lymphocytes are important prognostic factors in patients with TNBC, although the relative importance of lymphocyte subsets and associated protein expression is incompletely understood. NanoString Digital Spatial Profiling (DSP) allows for spatially resolved, highly multiplexed quantification of proteins in clinical samples. In this study, we sought to determine if DSP could be used to characterize expression of MHCII and other immune related proteins in tumor epithelial versus stromal compartments of patient-derived TNBCs (N = 10) using a panel of 39 markers. We confirmed that a subset of TNBCs have elevated expression of HLA-DR in tumor epithelial cells; HLA-DR expression was also significantly higher in the tumors of patients with long-term disease-free survival when compared to patients that relapsed. HLA-DR expression in the epithelial compartment was correlated with high expression of CD4 and ICOS in the stromal compartment of the same tumors. We also identified candidate protein biomarkers with significant differential expression between patients that relapsed versus those that did not. In conclusion, DSP is a powerful method that allows for quantification of proteins in the immune microenvironment of TNBCs.

TBCRC 002: a phase II, randomized, open-label trial of preoperative letrozole with or without bevacizumab in postmenopausal women with newly diagnosed stage 2/3 hormone receptor-positive and HER2-negative breast cancer.

BACKGROUND: In preclinical studies, the expression of vascular endothelial growth factor (VEGF) in hormone receptor-positive breast cancer is associated with estrogen-independent tumor growth and resistance to endocrine therapies. This study investigated whether the addition of bevacizumab, a monoclonal antibody against VEGF, to letrozole enhanced the antitumor activity of the letrozole in the preoperative setting. METHODS: Postmenopausal women with newly diagnosed stage 2 or 3 estrogen and/or progesterone receptor-positive, HER2-negative breast cancer were randomly assigned (2:1) between letrozole 2.5 mg PO daily plus bevacizumab 15 mg/kg IV every 3 weeks (Let/Bev) and letrozole 2.5 mg PO daily (Let) for 24 weeks prior to definitive surgery. Primary objective was within-arm pathologic complete remission (pCR) rate. Secondary objectives were safety, objective response, and downstaging rate. RESULTS: Seventy-five patients were randomized (Let/Bev n = 50, Let n = 25). Of the 45 patients evaluable for pathological response in the Let/Bev arm, 5 (11%; 95% CI, 3.7-24.1%) achieved pCR and 4 (9%; 95% CI, 2.5-21.2%) had microscopic residual disease; no pCRs or microscopic residual disease was seen in the Let arm (0%; 95% CI, 0-14.2%). The rates of downstaging were 44.4% (95% CI, 29.6-60.0%) and 37.5% (95% CI, 18.8-59.4%) in the Let/Bev and Let arms, respectively. Adverse events typically associated with letrozole (hot flashes, arthralgias, fatigue, myalgias) occurred in similar frequencies in the two arms. Hypertension, headache, and proteinuria were seen exclusively in the Let/Bev arm. The rates of grade 3 and 4 adverse events and discontinuation due to adverse events were 18% vs 8% and 16% vs none in the Let/Bev and Let arms, respectively. A small RNA-based classifier predictive of response to preoperative Let/Bev was developed and confirmed on an independent cohort. CONCLUSION: In the preoperative setting, the addition of bevacizumab to letrozole was associated with a pCR rate of 11%; no pCR was seen with letrozole alone. There was additive toxicity with the incorporation of bevacizumab. Responses to Let/Bev can be predicted from the levels of 5 small RNAs in a pretreatment biopsy. TRIAL REGISTRATION: This trial is registered with ClinicalTrials.gov (Identifier: NCT00161291), first posted on September 12, 2005, and is completed.

A Multigene Assay Determines Risk of Recurrence in Patients with Triple-Negative Breast Cancer.

Approximately 40% of patients with stage I-III triple-negative breast cancer (TNBC) recur after standard treatment, whereas the remaining 60% experience long-term disease-free survival (DFS). There are currently no clinical tests to assess the risk of recurrence in TNBC patients. We previously determined that TNBC patients with MHC class II (MHCII) pathway expression in their tumors experienced significantly longer DFS. To translate this discovery into a clinical test, we developed an MHCII Immune Activation assay, which measures expression of 36 genes using NanoString technology. Preanalytical testing confirmed that the assay is accurate and reproducible in formalin-fixed paraffin-embedded (FFPE) tumor specimens. The assay measurements were concordant with RNA-seq, MHCII protein expression, and tumor-infiltrating lymphocyte counts. In a training set of 44 primary TNBC tumors, the MHCII Immune Activation Score was significantly associated with longer DFS (HR = 0.17; P = 0.015). In an independent validation cohort of 56 primary FFPE TNBC tumors, the Immune Activation Score was significantly associated with longer DFS (HR = 0.19; P = 0.011) independent of clinical stage. An Immune Activation Score threshold for identifying patients with very low risk of relapse in the training set provided 100% specificity in the validation cohort. The assay format enables adoption as a standardized clinical prognostic test for identifying TNBC patients with a low risk of recurrence. Correlative data support future studies to determine if the assay can identify patients in whom chemotherapy can be safely deescalated and patients likely to respond to immunotherapy. SIGNIFICANCE: The MHCII Immune Activation assay identifies TNBC patients with a low risk of recurrence, addressing a critical need for prognostic biomarker tests that enable precision medicine for TNBC patients.

Single-Cell Transcriptomes Distinguish Stem Cell State Changes and Lineage Specification Programs in Early Mammary Gland Development.

The mammary gland consists of cells with gene expression patterns reflecting their cellular origins, function, and spatiotemporal context. However, knowledge of developmental kinetics and mechanisms of lineage specification is lacking. We address this significant knowledge gap by generating a single-cell transcriptome atlas encompassing embryonic, postnatal, and adult mouse mammary development. From these data, we map the chronology of transcriptionally and epigenetically distinct cell states and distinguish fetal mammary stem cells (fMaSCs) from their precursors and progeny. fMaSCs show balanced co-expression of factors associated with discrete adult lineages and a metabolic gene signature that subsides during maturation but reemerges in some human breast cancers and metastases. These data provide a useful resource for illuminating mammary cell heterogeneity, the kinetics of differentiation, and developmental correlates of tumorigenesis.

Genomic regulation of invasion by STAT3 in triple negative breast cancer.

Breast cancer is a heterogeneous disease comprised of four molecular subtypes defined by whether the tumor-originating cells are luminal or basal epithelial cells. Breast cancers arising from the luminal mammary duct often express estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth receptor 2 (HER2). Tumors expressing ER and/or PR are treated with anti-hormonal therapies, while tumors overexpressing HER2 are targeted with monoclonal antibodies. Immunohistochemical detection of ER, PR, and HER2 receptors/proteins is a critical step in breast cancer diagnosis and guided treatment. Breast tumors that do not express these proteins are known as "triple negative breast cancer" (TNBC) and are typically basal-like. TNBCs are the most aggressive subtype, with the highest mortality rates and no targeted therapy, so there is a pressing need to identify important TNBC tumor regulators. The signal transducer and activator of transcription 3 (STAT3) transcription factor has been previously implicated as a constitutively active oncogene in TNBC. However, its direct regulatory gene targets and tumorigenic properties have not been well characterized. By integrating RNA-seq and ChIP-seq data from 2 TNBC tumors and 5 cell lines, we discovered novel gene signatures directly regulated by STAT3 that were enriched for processes involving inflammation, immunity, and invasion in TNBC. Functional analysis revealed that STAT3 has a key role regulating invasion and metastasis, a characteristic often associated with TNBC. Our findings suggest therapies targeting STAT3 may be important for preventing TNBC metastasis.

Expression of the MHC Class II Pathway in Triple-Negative Breast Cancer Tumor Cells Is Associated with a Good Prognosis and Infiltrating Lymphocytes.

Triple-negative breast cancer (TNBC) is a subtype with heterogeneous patient outcomes. Approximately 40% of patients experience rapid relapse, while the remaining patients have long-term disease-free survival. To determine if there are molecular differences between primary tumors that predict prognosis, we performed RNA-seq on 47 macrodissected tumors from newly diagnosed patients with TNBC (n = 47; 22 relapse, 25 no relapse; follow-up median, 8 years; range, 2-11 years). We discovered that expression of the MHC class II (MHC II) antigen presentation pathway in tumor tissue was the most significant pathway associated with progression-free survival (HR, 0.36; log-rank P = 0.0098). The association between MHC II pathway expression and good prognosis was confirmed in a public gene expression database of 199 TNBC cases (HR, 0.28; log-rank P = 4.5 × 10(-8)). Further analysis of immunohistochemistry, laser-capture microdissected tumors, and TNBC cell lines demonstrated that tumor cells, in addition to immune cells, aberrantly express the MHC II pathway. MHC II pathway expression was also associated with B-cell and T-cell infiltration in the tumor. Together, these data support the model that aberrant expression of the MHC II pathway in TNBC tumor cells may trigger an antitumor immune response that reduces the rate of relapse and enhances progression-free survival. Cancer Immunol Res; 4(5); 390-9. ©2016 AACR.

TBCRC 019: A Phase II Trial of Nanoparticle Albumin-Bound Paclitaxel with or without the Anti-Death Receptor 5 Monoclonal Antibody Tigatuzumab in Patients with Triple-Negative Breast Cancer.

PURPOSE: Tigatuzumab (TIG), an agonistic anti-DR5 antibody, triggers apoptosis in DR5(+) human tumor cells without crosslinking. TIG has strong in vitro/in vivo activity against basal-like breast cancer cells enhanced by chemotherapy agents. This study evaluates activity of TIG and chemotherapy in patients with metastatic triple-negative breast cancer (TNBC). EXPERIMENTAL DESIGN: Randomized 2:1 phase II trial of albumin-bound paclitaxel (nab-PAC) ± TIG in patients with TNBC stratified by prior chemotherapy. Patients received nab-PAC weekly × 3 ± TIG every other week, every 28 days. Primary objective was within-arm objective response rate (ORR). Secondary objectives were safety, progression-free survival (PFS), clinical benefit, and TIG immunogenicity. Metastatic research biopsies were required. RESULTS: Among 64 patients (60 treated; TIG/nab-PAC n = 39 and nab-PAC n = 21), there were 3 complete remissions (CR), 8 partial remissions (PR; 1 almost CR), 11 stable diseases (SD), and 17 progressive diseases (PD) in the TIG/nab-PAC arm (ORR, 28%), and no CRs, 8 PRs, 4 SDs, and 9 PDs in the nab-PAC arm (ORR, 38%). There was a numerical increase in CRs and several patients had prolonged PFS (1,025+, 781, 672, 460, 334) in the TIG/nab-PAC arm. Grade 3 toxicities were 28% and 29%, respectively, with no grade 4-5. Exploratory analysis suggests an association of ROCK1 gene pathway activation with efficacy in the TIG/nab-PAC arm. CONCLUSIONS: ORR and PFS were similar in both. Preclinical activity of TIG in basal-like breast cancer and prolonged PFS in few patients in the combination arm support further investigation of anti-DR5 agents. ROCK pathway activation merits further evaluation.

Recurrent read-through fusion transcripts in breast cancer.

Read-through fusion transcripts that result from the splicing of two adjacent genes in the same coding orientation are a recently discovered type of chimeric RNA. We sought to determine if read-through fusion transcripts exist in breast cancer. We performed paired-end RNA-seq of 168 breast samples, including 28 breast cancer cell lines, 42 triple negative breast cancer primary tumors, 42 estrogen receptor positive (ER+) breast cancer primary tumors, and 56 non-malignant breast tissue samples. We analyzed the sequencing data to identify breast cancer associated read-through fusion transcripts. We discovered two recurrent read-through fusion transcripts that were identified in breast cancer cell lines, confirmed across breast cancer primary tumors, and were not detected in normal tissues (SCNN1A-TNFRSF1A and CTSD-IFITM10). Both fusion transcripts use canonical splice sites to join the last splice donor of the 5' gene to the first splice acceptor of the 3' gene, creating an in-frame fusion transcript. Western blots indicated that the fusion transcripts are translated into fusion proteins in breast cancer cells. Custom small interfering RNAs targeting the CTSD-IFITM10 fusion junction reduced expression of the fusion transcript and reduced breast cancer cell proliferation. Read-through fusion transcripts between adjacent genes with different biochemical functions represent a new type of recurrent molecular defect in breast cancer that warrant further investigation as potential biomarkers and therapeutic targets. Both breast cancer associated fusion transcripts identified in this study involve membrane proteins (SCNN1A-TNFRSF1A and CTSD-IFITM10), which raises the possibility that they could be breast cancer-specific cell surface markers.

Distinct properties of cell-type-specific and shared transcription factor binding sites.

Most human transcription factors bind a small subset of potential genomic sites and often use different subsets in different cell types. To identify mechanisms that govern cell-type-specific transcription factor binding, we used an integrative approach to study estrogen receptor α (ER). We found that ER exhibits two distinct modes of binding. Shared sites, bound in multiple cell types, are characterized by high-affinity estrogen response elements (EREs), inaccessible chromatin, and a lack of DNA methylation, while cell-specific sites are characterized by a lack of EREs, co-occurrence with other transcription factors, and cell-type-specific chromatin accessibility and DNA methylation. These observations enabled accurate quantitative models of ER binding that suggest tethering of ER to one-third of cell-specific sites. The distinct properties of cell-specific binding were also observed with glucocorticoid receptor and for ER in primary mouse tissues, representing an elegant genomic encoding scheme for generating cell-type-specific gene regulation.