Cost-effective DNA methylation profiling by FML-seq.

Current methods for profiling DNA methylation require costly reagents, sequencing, and labor time. We introduce fragmentation at methylated loci and sequencing (FML-seq), a sequencing library protocol that greatly reduces all these costs. Relative to other techniques tested on the same human cell lines, FML-seq produces similar measurements of absolute and differential cytosine methylation at a fraction of the price. FML-seq enables inexpensive, high-throughput experimental designs for large-scale epigenetics research projects.

Cost-effective DNA methylation profiling by FML-seq.

Current methods for profiling DNA methylation require costly reagents, sequencing, or labor time. We introduce FML-seq, a sequencing library protocol that greatly reduces all these costs. Relative to other techniques tested on the same human cell lines, FML-seq produces similar measurements of absolute and differential cytosine methylation at a fraction of the price. FML-seq enables inexpensive, high-throughput experimental designs for large-scale epigenetics research projects.

Molecular classification and biomarkers of clinical outcome in breast ductal carcinoma in situ: Analysis of TBCRC 038 and RAHBT cohorts.

Ductal carcinoma in situ (DCIS) is the most common precursor of invasive breast cancer (IBC), with variable propensity for progression. We perform multiscale, integrated molecular profiling of DCIS with clinical outcomes by analyzing 774 DCIS samples from 542 patients with 7.3 years median follow-up from the Translational Breast Cancer Research Consortium 038 study and the Resource of Archival Breast Tissue cohorts. We identify 812 genes associated with ipsilateral recurrence within 5 years from treatment and develop a classifier that predicts DCIS or IBC recurrence in both cohorts. Pathways associated with recurrence include proliferation, immune response, and metabolism. Distinct stromal expression patterns and immune cell compositions are identified. Our multiscale approach employed in situ methods to generate a spatially resolved atlas of breast precancers, where complementary modalities can be directly compared and correlated with conventional pathology findings, disease states, and clinical outcome.

The microdissected gene expression landscape of nasopharyngeal cancer reveals vulnerabilities in FGF and noncanonical NF-κB signaling.

Nasopharyngeal cancer (NPC) is an Epstein-Barr virus (EBV)-positive epithelial malignancy with an extensive inflammatory infiltrate. Traditional RNA-sequencing techniques uncovered only microenvironment signatures, while the gene expression of the tumor epithelial compartment has remained a mystery. Here, we use Smart-3SEQ to prepare transcriptome-wide gene expression profiles from microdissected NPC tumors, dysplasia, and normal controls. We describe changes in biological pathways across the normal to tumor spectrum and show that fibroblast growth factor (FGF) ligands are overexpressed in NPC tumors, while negative regulators of FGF signaling, including SPRY1, SPRY2, and LGALS3, are down-regulated early in carcinogenesis. Within the NF-κB signaling pathway, the critical noncanonical transcription factors, RELB and NFKB2, are enriched in the majority of NPC tumors. We confirm the responsiveness of EBV-positive NPC cell lines to targeted inhibition of these pathways, reflecting the heterogeneity in NPC patient tumors. Our data comprehensively describe the gene expression landscape of NPC and unravel the mysteries of receptor tyrosine kinase and NF-κB pathways in NPC.

Gene expression profiling of single cells from archival tissue with laser-capture microdissection and Smart-3SEQ.

RNA sequencing (RNA-seq) is a sensitive and accurate method for quantifying gene expression. Small samples or those whose RNA is degraded, such as formalin-fixed paraffin-embedded (FFPE) tissue, remain challenging to study with nonspecialized RNA-seq protocols. Here, we present a new method, Smart-3SEQ, that accurately quantifies transcript abundance even with small amounts of total RNA and effectively characterizes small samples extracted by laser-capture microdissection (LCM) from FFPE tissue. We also obtain distinct biological profiles from FFPE single cells, which have been impossible to study with previous RNA-seq protocols, and we use these data to identify possible new macrophage phenotypes associated with the tumor microenvironment. We propose Smart-3SEQ as a highly cost-effective method to enable large gene expression profiling experiments unconstrained by sample size and tissue availability. In particular, Smart-3SEQ's compatibility with FFPE tissue unlocks an enormous number of archived clinical samples; combined with LCM it allows unprecedented studies of small cell populations and single cells isolated by their in situ context.

Clinically Relevant Molecular Subtypes in Leiomyosarcoma.

PURPOSE: Leiomyosarcoma is a malignant neoplasm with smooth muscle differentiation. Little is known about its molecular heterogeneity and no targeted therapy currently exists for leiomyosarcoma. Recognition of different molecular subtypes is necessary to evaluate novel therapeutic options. In a previous study on 51 leiomyosarcomas, we identified three molecular subtypes in leiomyosarcoma. The current study was performed to determine whether the existence of these subtypes could be confirmed in independent cohorts. EXPERIMENTAL DESIGN: Ninety-nine cases of leiomyosarcoma were expression profiled with 3'end RNA-Sequencing (3SEQ). Consensus clustering was conducted to determine the optimal number of subtypes. RESULTS: We identified 3 leiomyosarcoma molecular subtypes and confirmed this finding by analyzing publically available data on 82 leiomyosarcoma from The Cancer Genome Atlas (TCGA). We identified two new formalin-fixed, paraffin-embedded tissue-compatible diagnostic immunohistochemical markers; LMOD1 for subtype I leiomyosarcoma and ARL4C for subtype II leiomyosarcoma. A leiomyosarcoma tissue microarray with known clinical outcome was used to show that subtype I leiomyosarcoma is associated with good outcome in extrauterine leiomyosarcoma while subtype II leiomyosarcoma is associated with poor prognosis in both uterine and extrauterine leiomyosarcoma. The leiomyosarcoma subtypes showed significant differences in expression levels for genes for which novel targeted therapies are being developed, suggesting that leiomyosarcoma subtypes may respond differentially to these targeted therapies. CONCLUSIONS: We confirm the existence of 3 molecular subtypes in leiomyosarcoma using two independent datasets and show that the different molecular subtypes are associated with distinct clinical outcomes. The findings offer an opportunity for treating leiomyosarcoma in a subtype-specific targeted approach.

Cell-lineage heterogeneity and driver mutation recurrence in pre-invasive breast neoplasia.

BACKGROUND: All cells in an individual are related to one another by a bifurcating lineage tree, in which each node is an ancestral cell that divided into two, each branch connects two nodes, and the root is the zygote. When a somatic mutation occurs in an ancestral cell, all its descendants carry the mutation, which can then serve as a lineage marker for the phylogenetic reconstruction of tumor progression. Using this concept, we investigate cell lineage relationships and genetic heterogeneity of pre-invasive neoplasias compared to invasive carcinomas. METHODS: We deeply sequenced over a thousand phylogenetically informative somatic variants in 66 morphologically independent samples from six patients that represent a spectrum of normal, early neoplasia, carcinoma in situ, and invasive carcinoma. For each patient, we obtained a highly resolved lineage tree that establishes the phylogenetic relationships among the pre-invasive lesions and with the invasive carcinoma. RESULTS: The trees reveal lineage heterogeneity of pre-invasive lesions, both within the same lesion, and between histologically similar ones. On the basis of the lineage trees, we identified a large number of independent recurrences of PIK3CA H1047 mutations in separate lesions in four of the six patients, often separate from the diagnostic carcinoma. CONCLUSIONS: Our analyses demonstrate that multi-sample phylogenetic inference provides insights on the origin of driver mutations, lineage heterogeneity of neoplastic proliferations, and the relationship of genomically aberrant neoplasias with the primary tumors. PIK3CA driver mutations may be comparatively benign inducers of cellular proliferation.

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.

Desktop transcriptome sequencing from archival tissue to identify clinically relevant translocations.

Somatic mutations, often translocations or single nucleotide variations, are pathognomonic for certain types of cancers and are increasingly of clinical importance for diagnosis and prediction of response to therapy. Conventional clinical assays only evaluate 1 mutation at a time, and targeted tests are often constrained to identify only the most common mutations. Genome-wide or transcriptome-wide high-throughput sequencing (HTS) of clinical samples offers an opportunity to evaluate for all clinically significant mutations with a single test. Recently a "desktop version" of HTS has become available, but most of the experience to date is based on data obtained from high-quality DNA from frozen specimens. In this study, we demonstrate, as a proof of principle, that translocations in sarcomas can be diagnosed from formalin-fixed paraffin-embedded (FFPE) tissue with desktop HTS. Using the first generation MiSeq platform, full transcriptome sequencing was performed on FFPE material from archival blocks of 3 synovial sarcomas, 3 myxoid liposarcomas, 2 Ewing sarcomas, and 1 clear cell sarcoma. Mapping the reads to the "sarcomatome" (all known 83 genes involved in translocations and mutations in sarcoma) and using a novel algorithm for ranking fusion candidates, the pathognomonic fusions and the exact breakpoints were identified in all cases of synovial sarcoma, myxoid liposarcoma, and clear cell sarcoma. The Ewing sarcoma fusion gene was detectable in FFPE material only with a sequencing platform that generates greater sequencing depth. The results show that a single transcriptome HTS assay, from FFPE, has the potential to replace conventional molecular diagnostic techniques for the evaluation of clinically relevant mutations in cancer.

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.

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.

Expression of subtype-specific group 1 leiomyosarcoma markers in a wide variety of sarcomas by gene expression analysis and immunohistochemistry.

Leiomyosarcomas (LMSs) constitute approximately one quarter of all sarcomas and are usually defined by morphologic criteria and/or immunoreactivity for actin or desmin. Among high-grade lesions, the distinction from undifferentiated pleomorphic sarcoma (UPS) can be problematic, and previous studies have shown that a significant number of LMS cases may be hiding under the diagnosis of UPS. We recently described 3 novel molecular LMS subtypes that are distributed similarly over LMSs of gyneocologic and non-gyneocologic origins. The group 1 subtype shows an improved disease-specific survival compared with the other 2 groups that is independent of histologic grade. Group 1 comprises approximately 25% of all LMSs, and is defined by a shared pattern of gene expression, a distinct pattern of genomic changes, and reactivity for at least 3 of 5 immunohistochemistry (IHC) markers (smooth muscle gamma actin, calsequestrin 2, human muscle cofilin2, myosin light chain kinase, and sarcolemmal membrane associated protein), as tested on 271 cases of LMS in tissue microarrays. These IHC markers have not been well characterized in non-LMS sarcomas. Here we provide a characterization of these 5 markers across normal tissues, an additional 59 cases of LMS, and a wide range of 565 non-LMS soft tissue tumors from 44 diagnostic categories, with a focus on UPS. When analyzed individually, the 5 markers were found to be expressed in many sarcomas other than LMSs. However, when analyzed by the same criteria used for the recognition of group 1 LMSs, in which a case is scored positive when at least 3 of 5 markers reacted, coordinate expression was seen in significant numbers of cases from only 3 diagnostic groups that included 22% of leiomyomas (n=22), 16% of gastrointestinal stromal tumors (n=43), and 18% of endometrial stromal sarcomas (n=11). In addition, 5% (n=57) of UPSs showed a staining pattern similar to that seen in group 1 LMSs. To further examine the possibility that group 1 LMS constitutes a small part of cases diagnosed as UPS, we examined the expression of the top 500 genes from the group 1 LMS expression signature in 29 UPSs by complementary DNA microarray. Unsupervised hierarchical clustering of 29 UPS expression showed that 2 (7%) had coordinated high levels of expression of genes from the group 1 LMS signature, a rate similar to that seen by IHC analysis. These findings show that group 1 LMS IHC markers smooth muscle gamma actin, calsequestrin 2, human muscle cofilin2, myosin light chain kinase, and sarcolemmal membrane associated protein when coordinately expressed have specificity for a subset of LMS when compared with other sarcomas, and may be useful for the recognition of group 1 LMS cases within cases diagnosed as UPS.

Endogenous versus tumor-specific host response to breast carcinoma: a study of stromal response in synchronous breast primaries and biopsy site changes.

PURPOSE: We recently described two types of stromal response in breast cancer derived from gene expression studies of tenosynovial giant cell tumors and fibromatosis. The purpose of this study is to elucidate the basis of this stromal response--whether they are elicited by individual tumors or whether they represent an endogenous host reaction produced by the patient. EXPERIMENTAL DESIGN: Stromal signatures from patients with synchronous dual primaries were analyzed by immunohistochemistry on a tissue microarray (n = 26 pairs) to evaluate the similarity of stromal responses in different tumors within the same patient. We also characterized the extent to which the stromal signatures were conserved between stromal response to injury compared to the stromal response to carcinoma using gene expression profiling and tissue microarray immunohistochemistry. RESULTS: The two stromal response signatures showed divergent associations in synchronous primaries: the DTF fibroblast response is more likely to be similar in a patient with multiple breast primaries (permutation analysis P = 0.0027), whereas CSF1 macrophage response shows no significant concordance in separate tumors within a given patient. The DTF fibroblast signature showed more concordance across normal, cancer, and biopsy site samples from within a patient, than across normal, cancer, and biopsy site samples from a random group of patients, whereas the CSF1 macrophage response did not. CONCLUSIONS: The results suggest that the DTF fibroblast response is host-specific, whereas the CSF1 response may be tumor-elicited. Our findings provide further insight into stromal response and may facilitate the development of therapeutic strategies to target particular stromal subtypes.

MicroRNA profiling of BRCA1/2 mutation-carrying and non-mutation-carrying high-grade serous carcinomas of ovary.

BACKGROUND: MicroRNAs (miRNA) are 20 approximately 25 nucleotide non-coding RNAs that inhibit the translation of targeted mRNA, and they have been implicated in the development of human malignancies. High grade serous ovarian carcinomas, the most common and lethal subtype of ovarian cancer, can occur sporadically or in the setting of BRCA1/2 syndromes. Little is known regarding the miRNA expression profiles of high grade serous carcinoma in relation to BRCA1/2 status, and compared to normal tubal epithelium, the putative tissue of origin for high grade serous carcinomas. METHODOLOGY/PRINCIPAL FINDINGS: Global miRNA expression profiling was performed on a series of 33 high grade serous carcinomas, characterized with respect to BRCA1/2 status (mutation, epigenetic silencing with loss of expression or normal), and with clinical follow-up, together with 2 low grade serous carcinomas, 2 serous borderline tumors, and 3 normal fallopian tube samples, using miRNA microarrays (328 human miRNA). Unsupervised hierarchical clustering based on miRNA expression profiles showed no clear separation between the groups of carcinomas with different BRCA1/2 status. There were relatively few miRNAs that were differentially expressed between the genotypic subgroups. Comparison of 33 high grade serous carcinomas to 3 normal fallopian tube samples identified several dysregulated miRNAs (false discovery rate <5%), including miR-422b and miR-34c. Quantitative RT-PCR analysis performed on selected miRNAs confirmed the pattern of differential expression shown by microarray analysis. Prognostically, lower level miR-422b and miR-34c in high grade serous carcinomas were both associated with decreased disease-specific survival by Kaplan-Meier analysis (p<0.05). CONCLUSIONS/SIGNIFICANCE: High grade serous ovarian carcinomas with and without BRCA1/2 abnormalities demonstrate very similar miRNA expression profiles. High grade serous carcinomas as a group exhibit significant miRNA dysregulation in comparison to tubal epithelium and the levels of miR-34c and miR-422b appear to be prognostically important.

Gene expression profiling identifies p63 as a diagnostic marker for giant cell tumor of the bone.

Giant cell tumor of the bone (GCTOB) is a primary bone tumor that occurs mainly in young adults and is capable of locally aggressive growth. Its histologic appearance can resemble a number of benign and malignant tumors but no useful diagnostic marker is known currently. To identify diagnostic markers for this tumor, global gene expression profiling using cDNA microarray was performed on 6 fresh-frozen GCTOB, 3 aneurysmal bone cysts, 4 fibrous dysplasias and 12 giant cell tumors of tendon sheath/diffuse-type giant cell tumors. Unsupervised hierarchical clustering separated the tumors based on their histopathologic types, and significance analysis of microarray identified several genes including TP73L (encoding the p63 protein) that are significantly highly expressed in GCTOB relative to these other tumors. The diagnostic utility of p63 was subsequently confirmed using anti-p63 antibody on a series of 26 GCTOB, 25 aneurysmal bone cysts, 15 chondroblastomas, 13 giant cell reparative granulomas, 13 chondromyxoid fibromas, 4 brown tumors, 4 fibrous dysplasias, 53 giant cell tumors of tendon sheath/diffuse-type giant cell tumors and 385 additional mesenchymal tumors in tissue microarrays. Strong p63 nuclear staining was present in 18 of 26 (69%) GCTOB, 3 of 15 (20%) chondroblastomas and in 1 of 25 (4%) aneurysmal bone cysts while none of the other tumors commonly considered in the differential diagnosis of GCTOB showed any detectable p63 staining. Strong p63 staining is rare in bone and soft-tissue tumors in general. In contrast to the pattern of p63 staining, the majority of the chondroblastomas (70%) demonstrated S-100 immunoreactivity while only a minority of the GCTOB (8%) was immunoreactive for S-100. These findings altogether show that p63 can be used as a diagnostic marker to aid the clinical diagnosis of GCTOB.

Placental S100 (S100P) and GATA3: markers for transitional epithelium and urothelial carcinoma discovered by complementary DNA microarray.

The morphologic distinction between prostate and urothelial carcinoma can be difficult. To identify novel diagnostic markers that may aid in the differential diagnosis of prostate versus urothelial carcinoma, we analyzed expression patterns in prostate and bladder cancer tissues using complementary DNA microarrays. Together with our prior studies on renal neoplasms and normal kidney, these studies suggested that the gene for placental S100 (S100P) is specifically expressed in benign and malignant urothelial cells. Using tissue microarrays, a polyclonal antiserum against S100P protein stained 86% of 295 urothelial carcinomas while only 3% of 260 prostatic adenocarcinomas and 1% of 133 renal cell carcinomas stained. A commercially available monoclonal antibody against S100P stained 78% of 300 urothelial carcinomas while only 2% of 256 prostatic adenocarcinomas and none of 137 renal cell carcinomas stained. A second gene, GATA3, also showed high level expression in urothelial tumors by cDNA array. A commercially available monoclonal antibody against GATA3 stained 67% of 308 urothelial carcinomas, but none of the prostate or renal carcinomas. For comparison, staining was also performed for p63 and cytokeratin 5/6. p63 stained 87% of urothelial carcinomas whereas CK5/6 stained 54%. Importantly, when S100P and p63 were combined 95% of urothelial carcinomas were labeled by one or both markers. We conclude that the detection of S100P and GATA3 protein expression may help distinguish urothelial carcinomas from other genitourinary neoplasms that enter into the differential diagnosis.