RSPO3 antagonism inhibits growth and tumorigenicity in colorectal tumors harboring common Wnt pathway mutations.

Activating mutations in the Wnt pathway are a characteristic feature of colorectal cancer (CRC). The R-spondin (RSPO) family is a group of secreted proteins that enhance Wnt signaling and RSPO2 and RSPO3 gene fusions have been reported in CRC. We have previously shown that Wnt pathway blockers exhibit potent combinatorial activity with taxanes to inhibit tumor growth. Here we show that RSPO3 antagonism synergizes with paclitaxel based chemotherapies in patient-derived xenograft models (PDX) with RSPO3 fusions and in tumors with common CRC mutations such as APC, ß-catenin, or RNF43. In these latter types of tumors that represent over 90% of CRC, RSPO3 is produced by stromal cells in the tumor microenvironment and the activating mutations appear to sensitize the tumors to Wnt-Rspo synergy. The combination of RSPO3 inhibition and taxane treatment provides an approach to effectively target oncogenic WNT signaling in a significant number of patients with colorectal and other intestinal cancers.

Long noncoding RNA EWSAT1-mediated gene repression facilitates Ewing sarcoma oncogenesis.

Chromosomal translocation that results in fusion of the genes encoding RNA-binding protein EWS and transcription factor FLI1 (EWS-FLI1) is pathognomonic for Ewing sarcoma. EWS-FLI1 alters gene expression through mechanisms that are not completely understood. We performed RNA sequencing (RNAseq) analysis on primary pediatric human mesenchymal progenitor cells (pMPCs) expressing EWS-FLI1 in order to identify gene targets of this oncoprotein. We determined that long noncoding RNA-277 (Ewing sarcoma-associated transcript 1 [EWSAT1]) is upregulated by EWS-FLI1 in pMPCs. Inhibition of EWSAT1 expression diminished the ability of Ewing sarcoma cell lines to proliferate and form colonies in soft agar, whereas EWSAT1 inhibition had no effect on other cell types tested. Expression of EWS-FLI1 and EWSAT1 repressed gene expression, and a substantial fraction of targets that were repressed by EWS-FLI1 were also repressed by EWSAT1. Analysis of RNAseq data from primary human Ewing sarcoma further supported a role for EWSAT1 in mediating gene repression. We identified heterogeneous nuclear ribonucleoprotein (HNRNPK) as an RNA-binding protein that interacts with EWSAT1 and found a marked overlap in HNRNPK-repressed genes and those repressed by EWS-FLI1 and EWSAT1, suggesting that HNRNPK participates in EWSAT1-mediated gene repression. Together, our data reveal that EWSAT1 is a downstream target of EWS-FLI1 that facilitates the development of Ewing sarcoma via the repression of target genes.

Co-regulation of invected and engrailed by a complex array of regulatory sequences in Drosophila.

invected (inv) and engrailed (en) form a gene complex that extends about 115 kb. These two genes encode highly related homeodomain proteins that are co-regulated in a complex manner throughout development. Our dissection of inv/en regulatory DNA shows that most enhancers are spread throughout a 62 kb region. We used two types of constructs to analyze the function of this DNA: P-element based reporter constructs with small pieces of DNA fused to the en promoter driving lacZ expression and large constructs with HA-tagged en and inv inserted in the genome with the phiC31 system. In addition, we generated deletions of inv and en DNA in situ and assayed their effects on inv/en expression. Our results support and extend our knowledge of inv/en regulation. First, inv and en share regulatory DNA, most of which is flanking the en transcription unit. In support of this, a 79-kb HA-en transgene can rescue inv en double mutants to viable, fertile adults. In contrast, an 84-kb HA-inv transgene lacks most of the enhancers for inv/en expression. Second, there are multiple enhancers for inv/en stripes in embryos; some of these may be redundant but others play discrete roles at different stages of embryonic development. Finally, no small reporter construct gave expression in the posterior compartment of imaginal discs, a hallmark of inv/en expression. Robust expression of HA-en in the posterior compartment of imaginal discs is evident from the 79-kb HA-en transgene, while a 45-kb HA-en transgene gives weaker, variable imaginal disc expression. We suggest that the activity of the imaginal disc enhancer(s) is dependent on the chromatin structure of the inv/en domain.

A shared transcriptional program in early breast neoplasias despite genetic and clinical distinctions.

BACKGROUND: The earliest recognizable stages of breast neoplasia are lesions that represent a heterogeneous collection of epithelial proliferations currently classified based on morphology. Their role in the development of breast cancer is not well understood but insight into the critical events at this early stage will improve efforts in breast cancer detection and prevention. These microscopic lesions are technically difficult to study so very little is known about their molecular alterations. RESULTS: To characterize the transcriptional changes of early breast neoplasia, we sequenced 3'- end enriched RNAseq libraries from formalin-fixed paraffin-embedded tissue of early neoplasia samples and matched normal breast and carcinoma samples from 25 patients. We find that gene expression patterns within early neoplasias are distinct from both normal and breast cancer patterns and identify a pattern of pro-oncogenic changes, including elevated transcription of ERBB2, FOXA1, and GATA3 at this early stage. We validate these findings on a second independent gene expression profile data set generated by whole transcriptome sequencing. Measurements of protein expression by immunohistochemistry on an independent set of early neoplasias confirms that ER pathway regulators FOXA1 and GATA3, as well as ER itself, are consistently upregulated at this early stage. The early neoplasia samples also demonstrate coordinated changes in long non-coding RNA expression and microenvironment stromal gene expression patterns. CONCLUSIONS: This study is the first examination of global gene expression in early breast neoplasia, and the genes identified here represent candidate participants in the earliest molecular events in the development of breast cancer.

Next generation sequencing-based expression profiling identifies signatures from benign stromal proliferations that define stromal components of breast cancer.

INTRODUCTION: Multiple studies have shown that the tumor microenvironment (TME) of carcinomas can play an important role in the initiation, progression, and metastasis of cancer. Here we test the hypothesis that specific benign fibrous soft tissue tumor gene expression profiles may represent distinct stromal fibroblastic reaction types that occur in different breast cancers. The discovered stromal profiles could classify breast cancer based on the type of stromal reaction patterns in the TME. METHODS: Next generation sequencing-based gene expression profiling (3SEQ) was performed on formalin fixed, paraffin embedded (FFPE) samples of 10 types of fibrous soft tissue tumors. We determined the extent to which these signatures could identify distinct subsets of breast cancers in four publicly available breast cancer datasets. RESULTS: A total of 53 fibrous tumors were sequenced by 3SEQ with an average of 29 million reads per sample. Both the gene signatures derived from elastofibroma (EF) and fibroma of tendon sheath (FOTS) demonstrated robust outcome results for survival in the four breast cancer datasets. The breast cancers positive for the EF signature (20-33% of the cohort) demonstrated significantly better outcome for survival. In contrast, the FOTS signature-positive breast cancers (11-35% of the cohort) had a worse outcome. CONCLUSIONS: We defined and validated two new stromal signatures in breast cancer (EF and FOTS), which are significantly associated with prognosis. Our group has previously identified novel cancer stromal gene expression signatures associated with outcome differences in breast cancer by gene expression profiling of three soft tissue tumors, desmoid-type fibromatosis (DTF), solitary fibrous tumor (SFT), and tenosynovial giant cell tumor (TGCT/CSF1), as surrogates for stromal expression patterns. By combining the stromal signatures of EF and FOTS, with our previously identified DTF and TGCT/CSF1 signatures we can now characterize clinically relevant stromal expression profiles in the TME for between 74% to 90% of all breast cancers.

Genome evolution during progression to breast cancer.

Cancer evolution involves cycles of genomic damage, epigenetic deregulation, and increased cellular proliferation that eventually culminate in the carcinoma phenotype. Early neoplasias, which are often found concurrently with carcinomas and are histologically distinguishable from normal breast tissue, are less advanced in phenotype than carcinomas and are thought to represent precursor stages. To elucidate their role in cancer evolution we performed comparative whole-genome sequencing of early neoplasias, matched normal tissue, and carcinomas from six patients, for a total of 31 samples. By using somatic mutations as lineage markers we built trees that relate the tissue samples within each patient. On the basis of these lineage trees we inferred the order, timing, and rates of genomic events. In four out of six cases, an early neoplasia and the carcinoma share a mutated common ancestor with recurring aneuploidies, and in all six cases evolution accelerated in the carcinoma lineage. Transition spectra of somatic mutations are stable and consistent across cases, suggesting that accumulation of somatic mutations is a result of increased ancestral cell division rather than specific mutational mechanisms. In contrast to highly advanced tumors that are the focus of much of the current cancer genome sequencing, neither the early neoplasia genomes nor the carcinomas are enriched with potentially functional somatic point mutations. Aneuploidies that occur in common ancestors of neoplastic and tumor cells are the earliest events that affect a large number of genes and may predispose breast tissue to eventual development of invasive carcinoma.

Cyclin D1 as a diagnostic immunomarker for endometrial stromal sarcoma with YWHAE-FAM22 rearrangement.

Endometrial stromal sarcoma (ESS) characterized by YWHAE-FAM22 genetic fusion is histologically higher grade and clinically more aggressive than ESS with JAZF1-SUZ12 or equivalent genetic rearrangements, hence it is clinically important to recognize this subset of ESS. To identify diagnostic immunomarkers for this biologically defined ESS subset, we compared gene expression profiles between YWHAE-FAM22 ESS and JAZF1-rearranged ESS. These studies showed consistent upregulation of cyclin D1 in YWHAE-FAM22 ESS compared with JAZF1-SUZ12 ESS. Immunohistochemically, the high-grade round cell component of all 12 YWHAE-FAM22 ESS demonstrated diffuse (≥70%) moderate to strong nuclear cyclin D1 staining, and this diffuse positivity was not seen in 34 ESSs with JAZF1 and equivalent genetic rearrangements or in 21 low-grade ESS with no demonstrable genetic rearrangements. In a series of 243 non-ESS pure uterine mesenchymal and mixed epithelial-mesenchymal tumors, only 2 of 8 undifferentiated endometrial sarcomas with nuclear uniformity and 1 of 80 uterine leiomyosarcomas demonstrate diffuse cyclin D1 immunoreactivity. Both cyclin D1-positive undifferentiated endometrial sarcomas showed diffuse strong CD10 staining, which is consistently absent in the high-grade round cell component of YWHAE-FAM22 ESS. The low-grade spindle cell component of YWHAE-FAM22 ESS showed a spatially heterogenous cyclin D1 staining pattern that was weaker and less diffuse overall. Our findings indicate that cyclin D1 is a sensitive and specific diagnostic immunomarker for YWHAE-FAM22 ESS. When evaluating high-grade uterine sarcomas, cyclin D1 can be included in the immunohistochemical panel as an indicator of YWHAE-FAM22 ESS.

Transcriptional profiling of long non-coding RNAs and novel transcribed regions across a diverse panel of archived human cancers.

BACKGROUND: Molecular characterization of tumors has been critical for identifying important genes in cancer biology and for improving tumor classification and diagnosis. Long non-coding RNAs, as a new, relatively unstudied class of transcripts, provide a rich opportunity to identify both functional drivers and cancer-type-specific biomarkers. However, despite the potential importance of long non-coding RNAs to the cancer field, no comprehensive survey of long non-coding RNA expression across various cancers has been reported. RESULTS: We performed a sequencing-based transcriptional survey of both known long non-coding RNAs and novel intergenic transcripts across a panel of 64 archival tumor samples comprising 17 diagnostic subtypes of adenocarcinomas, squamous cell carcinomas and sarcomas. We identified hundreds of transcripts from among the known 1,065 long non-coding RNAs surveyed that showed variability in transcript levels between the tumor types and are therefore potential biomarker candidates. We discovered 1,071 novel intergenic transcribed regions and demonstrate that these show similar patterns of variability between tumor types. We found that many of these differentially expressed cancer transcripts are also expressed in normal tissues. One such novel transcript specifically expressed in breast tissue was further evaluated using RNA in situ hybridization on a panel of breast tumors. It was shown to correlate with low tumor grade and estrogen receptor expression, thereby representing a potentially important new breast cancer biomarker. CONCLUSIONS: This study provides the first large survey of long non-coding RNA expression within a panel of solid cancers and also identifies a number of novel transcribed regions differentially expressed across distinct cancer types that represent candidate biomarkers for future research.

Phosphatidylinositol-3-kinase pathway mutations are common in breast columnar cell lesions.

The phosphatidylinositol-3-kinase pathway is one of the most commonly mutated pathways in invasive breast carcinoma, with PIK3CA mutations in ∼25% of invasive carcinomas, and AKT1 mutations in up to 5%. Ductal carcinoma in situ, and benign papillomas harbor similar mutations. However, activating point mutations in breast columnar cell lesions have been infrequently studied. Twenty-three breast resection specimens containing columnar cell lesions were identified; 14 with associated invasive carcinoma or carcinoma in situ. DNA extracts were prepared from formalin-fixed paraffin-embedded tissue and screened for a panel of point mutations (321 mutations in 30 genes) using a multiplex PCR panel with mass-spectroscopy readout. PIK3CA mutations were identified in 13/24 columnar cell lesions (54%) and 3/8 invasive carcinomas (37%). The mutation status of columnar cell lesions and associated carcinoma was frequently discordant. Of the 14 cases, only 5 demonstrated the same genotype in matched samples of columnar cell lesions and carcinoma (4 wild type, 1 PIK3CA H1047R). Interestingly, five patients had mutations in columnar cell lesions with wild-type carcinoma; two patients had different point mutations in columnar cell lesions and carcinoma. Only three cases had wild-type columnar cell lesion and mutated carcinoma. The 50% PIK3CA mutation prevalence in columnar cell lesions is greater than reported in most studies of invasive breast cancer. Further, columnar cell lesions and carcinoma were frequently discordant for PIK3CA/AKT1 mutation status. These findings raise interesting questions about the role of PIK3CA/AKT pathway in breast carcinogenesis, and the biologic/precursor potential of columnar cell lesions.

14-3-3 fusion oncogenes in high-grade endometrial stromal sarcoma.

14-3-3 proteins are ubiquitously expressed regulators of various cellular functions, including proliferation, metabolism, and differentiation, and altered 14-3-3 expression is associated with development and progression of cancer. We report a transforming 14-3-3 oncoprotein, which we identified through conventional cytogenetics and whole-transcriptome sequencing analysis as a highly recurrent genetic mechanism in a clinically aggressive form of uterine sarcoma: high-grade endometrial stromal sarcoma (ESS). The 14-3-3 oncoprotein results from a t(10;17) genomic rearrangement, leading to fusion between 14-3-3ε (YWHAE) and either of two nearly identical FAM22 family members (FAM22A or FAM22B). Expression of YWHAE-FAM22 fusion oncoproteins was demonstrated by immunoblot in t(10;17)-bearing frozen tumor and cell line samples. YWHAE-FAM22 fusion gene knockdowns were performed with shRNAs and siRNAs targeting various FAM22A exons in an t(10;17)-bearing ESS cell line (ESS1): Fusion protein expression was inhibited, with corresponding reduction in cell growth and migration. YWHAE-FAM22 maintains a structurally and functionally intact 14-3-3ε (YWHAE) protein-binding domain, which is directed to the nucleus by a FAM22 nuclear localization sequence. In contrast to classic ESS, harboring JAZF1 genetic fusions, YWHAE-FAM22 ESS display high-grade histologic features, a distinct gene-expression profile, and a more aggressive clinical course. Fluorescence in situ hybridization analysis demonstrated absolute specificity of YWHAE-FAM22A/B genetic rearrangement for high-grade ESS, with no fusions detected in other uterine and nonuterine mesenchymal tumors (55 tumor types, n = 827). These discoveries reveal diagnostically and therapeutically relevant models for characterizing aberrant 14-3-3 oncogenic functions.

Distinct DNA methylation patterns characterize differentiated human embryonic stem cells and developing human fetal liver.

To investigate the role of DNA methylation during human development, we developed Methyl-seq, a method that assays DNA methylation at more than 90,000 regions throughout the genome. Performing Methyl-seq on human embryonic stem cells (hESCs), their derivatives, and human tissues allowed us to identify several trends during hESC and in vivo liver differentiation. First, differentiation results in DNA methylation changes at a minimal number of assayed regions, both in vitro and in vivo (2%-11%). Second, in vitro hESC differentiation is characterized by both de novo methylation and demethylation, whereas in vivo fetal liver development is characterized predominantly by demethylation. Third, hESC differentiation is uniquely characterized by methylation changes specifically at H3K27me3-occupied regions, bivalent domains, and low density CpG promoters (LCPs), suggesting that these regions are more likely to be involved in transcriptional regulation during hESC differentiation. Although both H3K27me3-occupied domains and LCPs are also regions of high variability in DNA methylation state during human liver development, these regions become highly unmethylated, which is a distinct trend from that observed in hESCs. Taken together, our results indicate that hESC differentiation has a unique DNA methylation signature that may not be indicative of in vivo differentiation.