The circular RNome of primary breast cancer.

Circular RNAs (circRNAs) are a class of RNAs that is under increasing scrutiny, although their functional roles are debated. We analyzed RNA-seq data of 348 primary breast cancers and developed a method to identify circRNAs that does not rely on unmapped reads or known splice junctions. We identified 95,843 circRNAs, of which 20,441 were found recurrently. Of the circRNAs that match exon boundaries of the same gene, 668 showed a poor or even negative (R < 0.2) correlation with the expression level of the linear gene. In silico analysis showed only a minority (8.5%) of circRNAs could be explained by known splicing events. Both these observations suggest that specific regulatory processes for circRNAs exist. We confirmed the presence of circRNAs of CNOT2, CREBBP, and RERE in an independent pool of primary breast cancers. We identified circRNA profiles associated with subgroups of breast cancers and with biological and clinical features, such as amount of tumor lymphocytic infiltrate and proliferation index. siRNA-mediated knockdown of circCNOT2 was shown to significantly reduce viability of the breast cancer cell lines MCF-7 and BT-474, further underlining the biological relevance of circRNAs. Furthermore, we found that circular, and not linear, CNOT2 levels are predictive for progression-free survival time to aromatase inhibitor (AI) therapy in advanced breast cancer patients, and found that circCNOT2 is detectable in cell-free RNA from plasma. We showed that circRNAs are abundantly present, show characteristics of being specifically regulated, are associated with clinical and biological properties, and thus are relevant in breast cancer.

Gene length corrected trimmed mean of M-values (GeTMM) processing of RNA-seq data performs similarly in intersample analyses while improving intrasample comparisons.

BACKGROUND: Current normalization methods for RNA-sequencing data allow either for intersample comparison to identify differentially expressed (DE) genes or for intrasample comparison for the discovery and validation of gene signatures. Most studies on optimization of normalization methods typically use simulated data to validate methodologies. We describe a new method, GeTMM, which allows for both inter- and intrasample analyses with the same normalized data set. We used actual (i.e. not simulated) RNA-seq data from 263 colon cancers (no biological replicates) and used the same read count data to compare GeTMM with the most commonly used normalization methods (i.e. TMM (used by edgeR), RLE (used by DESeq2) and TPM) with respect to distributions, effect of RNA quality, subtype-classification, recurrence score, recall of DE genes and correlation to RT-qPCR data. RESULTS: We observed a clear benefit for GeTMM and TPM with regard to intrasample comparison while GeTMM performed similar to TMM and RLE normalized data in intersample comparisons. Regarding DE genes, recall was found comparable among the normalization methods, while GeTMM showed the lowest number of false-positive DE genes. Remarkably, we observed limited detrimental effects in samples with low RNA quality. CONCLUSIONS: We show that GeTMM outperforms established methods with regard to intrasample comparison while performing equivalent with regard to intersample normalization using the same normalized data. These combined properties enhance the general usefulness of RNA-seq but also the comparability to the many array-based gene expression data in the public domain.

Confirmation of a metastasis-specific microRNA signature in primary colon cancer.

The identification of patients with high-risk stage II colon cancer who may benefit from adjuvant therapy may allow the clinical approach to be tailored for these patients based on an understanding of tumour biology. MicroRNAs have been proposed as markers of the prognosis or treatment response in colorectal cancer. Recently, a 2-microRNA signature (let-7i and miR-10b) was proposed to identify colorectal cancer patients at risk of developing distant metastasis. We assessed the prognostic value of this signature and additional candidate microRNAs in an independent, clinically well-defined, prospectively collected cohort of primary colon cancer patients including stage I-II colon cancer without and stage III colon cancer with adjuvant treatment. The 2-microRNA signature specifically predicted hepatic recurrence in the stage I-II group, but not the overall ability to develop distant metastasis. The addition of miR-30b to the 2-microRNA signature allowed the prediction of both distant metastasis and hepatic recurrence in patients with stage I-II colon cancer who did not receive adjuvant chemotherapy. Available gene expression data allowed us to associate miR-30b expression with axon guidance and let-7i expression with cell adhesion, migration, and motility.

High mRNA expression of splice variant SYK short correlates with hepatic disease progression in chemonaive lymph node negative colon cancer patients.

OBJECTIVE: Overall and splice specific expression of Spleen Tyrosine Kinase (SYK) has been posed as a marker predicting both poor and favorable outcome in various epithelial malignancies. However, its role in colorectal cancer is largely unknown. The aim of this study was to explore the prognostic role of SYK in three cohorts of colon cancer patients. METHODS: Total messenger RNA (mRNA) expression of SYK, SYK(T), and mRNA expression of its two splice variants SYK short (S) and SYK long (L) were measured using quantitative reverse transcriptase (RT-qPCR) in 240 primary colon cancer patients (n = 160 patients with chemonaive lymph node negative [LNN] and n = 80 patients with adjuvant treated lymph node positive [LNP] colon cancer) and related to microsatellite instability (MSI), known colorectal cancer mutations, and disease-free (DFS), hepatic metastasis-free (HFS) and overall survival (OS). Two independent cohorts of patients with respectively 48 and 118 chemonaive LNN colon cancer were used for validation. RESULTS: Expression of SYK and its splice variants was significantly lower in tumors with MSI, and in KRAS wild type, BRAF mutant and PTEN mutant tumors. In a multivariate Cox regression analysis, as a continuous variable, increasing SYK(S) mRNA expression was associated with worse HFS (Hazard Ratio[HR] = 1.83; 95% Confidence Interval[CI] = 1.08-3.12; p = 0.026) in the LNN group, indicating a prognostic role for SYK(S) mRNA in patients with chemonaive LNN colon cancer. However, only a non-significant trend between SYK(S) and HFS in one of the two validation cohorts was observed (HR = 4.68; 95%CI = 0.75-29.15; p = 0.098). CONCLUSION: In our cohort, we discovered SYK(S) as a significant prognostic marker for HFS for patients with untreated LNN colon cancer. This association could however not be confirmed in two independent smaller cohorts, suggesting that further extensive validation is needed to confirm the prognostic value of SYK(S) expression in chemonaive LNN colon cancer.

A Systematic Analysis of Oncogenic Gene Fusions in Primary Colon Cancer.

Genomic rearrangements that give rise to oncogenic gene fusions can offer actionable targets for cancer therapy. Here we present a systematic analysis of oncogenic gene fusions among a clinically well-characterized, prospectively collected set of 278 primary colon cancers spanning diverse tumor stages and clinical outcomes. Gene fusions and somatic genetic variations were identified in fresh frozen clinical specimens by Illumina RNA-sequencing, the STAR fusion gene detection pipeline, and GATK RNA-seq variant calling. We considered gene fusions to be pathogenically relevant when recurrent, producing divergent gene expression (outlier analysis), or as functionally important (e.g., kinase fusions). Overall, 2.5% of all specimens were defined as harboring a relevant gene fusion (kinase fusions 1.8%). Novel configurations of BRAF, NTRK3, and RET gene fusions resulting from chromosomal translocations were identified. An R-spondin fusion was found in only one tumor (0.35%), much less than an earlier reported frequency of 10% in colorectal cancers. We also found a novel fusion involving USP9X-ERAS formed by chromothripsis and leading to high expression of ERAS, a constitutively active RAS protein normally expressed only in embryonic stem cells. This USP9X-ERAS fusion appeared highly oncogenic on the basis of its ability to activate AKT signaling. Oncogenic fusions were identified only in lymph node-negative tumors that lacked BRAF or KRAS mutations. In summary, we identified several novel oncogenic gene fusions in colorectal cancer that may drive malignant development and offer new targets for personalized therapy. Cancer Res; 77(14); 3814-22. ©2017 AACR.

Breast cancer genome and transcriptome integration implicates specific mutational signatures with immune cell infiltration.

A recent comprehensive whole genome analysis of a large breast cancer cohort was used to link known and novel drivers and substitution signatures to the transcriptome of 266 cases. Here, we validate that subtype-specific aberrations show concordant expression changes for, for example, TP53, PIK3CA, PTEN, CCND1 and CDH1. We find that CCND3 expression levels do not correlate with amplification, while increased GATA3 expression in mutant GATA3 cancers suggests GATA3 is an oncogene. In luminal cases the total number of substitutions, irrespective of type, associates with cell cycle gene expression and adverse outcome, whereas the number of mutations of signatures 3 and 13 associates with immune-response specific gene expression, increased numbers of tumour-infiltrating lymphocytes and better outcome. Thus, while earlier reports imply that the sheer number of somatic aberrations could trigger an immune-response, our data suggests that substitutions of a particular type are more effective in doing so than others.

Molecular characteristics of circulating tumor cells resemble the liver metastasis more closely than the primary tumor in metastatic colorectal cancer.

BACKGROUND: CTCs are a promising alternative for metastatic tissue biopsies for use in precision medicine approaches. We investigated to what extent the molecular characteristics of circulating tumor cells (CTCs) resemble the liver metastasis and/or the primary tumor from patients with metastatic colorectal cancer (mCRC). RESULTS: The CTC profiles were concordant with the liver metastasis in 17/23 patients (74%) and with the primary tumor in 13 patients (57%). The CTCs better resembled the liver metastasis in 13 patients (57%), and the primary tumor in five patients (22%). The strength of the correlations was not associated with clinical parameters. Nine genes (CDH1, CDH17, CDX1, CEACAM5, FABP1, FCGBP, IGFBP3, IGFBP4, and MAPT) displayed significant differential expressions, all of which were downregulated, in CTCs compared to the tissues in the 23 patients. PATIENTS AND METHODS: Patients were retrospectively selected from a prospective study. Using the CellSearch System, CTCs were enumerated and isolated just prior to liver metastasectomy. A panel of 25 CTC-specific genes was measured by RT-qPCR in matching CTCs, primary tumors, and liver metastases. Spearman correlation coefficients were calculated and considered as continuous variables with r=1 representing absolute concordance and r=-1 representing absolute discordance. A cut-off of r>0.1 was applied in order to consider profiles to be concordant. CONCLUSIONS: In the majority of the patients, CTCs reflected the molecular characteristics of metastatic cells better than the primary tumors. Genes involved in cell adhesion and epithelial-to-mesenchymal transition were downregulated in the CTCs. Our results support the use of CTC characterization as a liquid biopsy for precision medicine.

An 8-gene mRNA expression profile in circulating tumor cells predicts response to aromatase inhibitors in metastatic breast cancer patients.

BACKGROUND: Molecular characterization of circulating tumor cells (CTC) is promising for personalized medicine. We aimed to identify a CTC gene expression profile predicting outcome to first-line aromatase inhibitors in metastatic breast cancer (MBC) patients. METHODS: CTCs were isolated from 78 MBC patients before treatment start. mRNA expression levels of 96 genes were measured by quantitative reverse transcriptase polymerase chain reaction. After applying predefined exclusion criteria based on lack of sufficient RNA quality and/or quantity, the data from 45 patients were used to construct a gene expression profile to predict poor responding patients, defined as disease progression or death <9 months, by a leave-one-out cross validation. RESULTS: Of the 45 patients, 19 were clinically classified as poor responders. To identify them, the 75% most variable genes were used to select genes differentially expressed between good and poor responders. An 8-gene CTC predictor was significantly associated with outcome (Hazard Ratio [HR] 4.40, 95% Confidence Interval [CI]: 2.17-8.92, P < 0.001). This predictor identified poor responding patients with a sensitivity of 63% and a positive predictive value of 75%, while good responding patients were correctly predicted in 85% of the cases. In multivariate Cox regression analysis, including CTC count at baseline, the 8-gene CTC predictor was the only factor independently associated with outcome (HR 4.59 [95% CI: 2.11-9.56], P < 0.001). This 8-gene signature was not associated with outcome in a group of 71 MBC patients treated with systemic treatments other than AI. CONCLUSIONS: An 8-gene CTC predictor was identified which discriminates good and poor outcome to first-line aromatase inhibitors in MBC patients. Although results need to be validated, this study underscores the potential of molecular characterization of CTCs.

Gene expression profiles of circulating tumor cells versus primary tumors in metastatic breast cancer.

Before using circulating tumor cells (CTCs) as liquid biopsy, insight into molecular discrepancies between CTCs and primary tumors is essential. We characterized CellSearch-enriched CTCs from 62 metastatic breast cancer (MBC) patients with ≥5 CTCs starting first-line systemic treatment. Expression levels of 35 tumor-associated, CTC-specific genes, including ESR1, coding for the estrogen receptor (ER), were measured by reverse transcription quantitative polymerase chain reaction and correlated to corresponding primary tumors. In 30 patients (48%), gene expression profiles of 35 genes were discrepant between CTCs and the primary tumor, but this had no prognostic consequences. In 15 patients (24%), the expression of ER was discrepant. Patients with ER-negative primary tumors and ER-positive CTCs had a longer median TTS compared to those with concordantly ER-negative CTCs (8.5 versus 2.1 months, P = 0.05). From seven patients, an axillary lymph node metastasis was available. In two patients, the CTC profiles better resembled the lymph node metastasis than the primary tumor. Our findings suggest that molecular discordances between CTCs and primary tumors frequently occur, but that this bears no prognostic consequences. Alterations in ER-status between primary tumors and CTCs might have prognostic implications.

mRNA expression profiles in circulating tumor cells of metastatic colorectal cancer patients.

INTRODUCTION: The molecular characterization of circulating tumor cells (CTCs) is a promising tool for the repeated and non-invasive evaluation of predictive and prognostic factors. Challenges associated with CTC characterization using the only FDA approved method for CTC enumeration, the CellSearch technique, include the presence of an excess of leukocytes in CTC-enriched blood fractions. Here we aimed to identify colorectal tumor-specific gene expression levels in the blood of patients with and without detectable CTCs according to CellSearch criteria. MATERIALS AND METHODS: Blood of 30 healthy donors (HDs) and 142 metastatic colorectal cancer (mCRC) patients was subjected to CellSearch CTC enumeration and isolation. In all samples, 95 mRNAs were measured by reverse transcriptase quantitative PCR (RT-qPCR). HD blood samples and patient samples with three or more CTCs were compared to identify CTC-specific mRNAs. Patient samples without detectable CTCs were separately analyzed. RESULTS: Thirty-four CTC-specific mRNAs were higher expressed in patients with ≥3 CTCs compared with HDs (Mann-Whitney U-test P < 0.05). Among patients without detectable CTCs, a HD-unlike subgroup was identified which could be distinguished from HDs by the expression of epithelial genes such as KRT19, KRT20 and AGR2. Also, in an independent patient set, a similar HD-unlike group could be identified among the patients without detectable CTCs according to the CellSearch system. CONCLUSION: Extensive molecular characterization of colorectal CTCs is feasible and a subgroup of patients without detectable CTCs according to CellSearch criteria bears circulating tumor load, which may have clinical consequences. This CTC-specific gene panel for mCRC patients may enable the exploration of CTC characterization as a novel means to further individualize cancer treatment.

Integrative analysis of genomics and proteomics data on clinical breast cancer tissue specimens extracted with acid guanidinium thiocyanate-phenol-chloroform.

Acid guanidinium thiocyanate, phenol, and chloroform extraction (AGPC) is a commonly used procedure to extract RNA from fresh frozen tissues and cell lines. In addition, DNA and proteins can be recovered, which makes AGPC an attractive source for integrative analysis on tissues of which little material is available, such as clinical specimens. Despite this potential, AGPC has only scarcely been used for proteomic analysis, mainly due to difficulties in extracting proteins. We have used a quantitative mass spectrometry method to show that proteins can readily be recovered from AGPC extracted tissues with high recovery and repeatability, which allows this method to be used for global proteomic profiling. Protein expression data for a selected number of clinically relevant markers, of which transcript and protein levels are known to be correlated, were in agreement with genomic and transcriptomic data obtained from the same AGPC lysate. Furthermore, global proteomic profiling successfully discriminated breast tumor tissues according to their clinical subtype. Lastly, a reference gene set of differentially expressed transcripts was strongly enriched in the differentially abundant proteins in our cohort. AGPC lysates are therefore well suited for comparative protein and integrative analyses.

KRAS and BRAF mutation status in circulating colorectal tumor cells and their correlation with primary and metastatic tumor tissue.

Although anti-EGFR therapy has established efficacy in metastatic colorectal cancer, only 10-20% of unselected patients respond. This is partly due to KRAS and BRAF mutations, which are currently assessed in the primary tumor. To improve patient selection, assessing mutation status in circulating tumor cells (CTCs), which possibly better represent metastases than the primary tumor, could be advantageous. We investigated the feasibility of KRAS and BRAF mutation detection in colorectal CTCs by comparing three sensitive methods and compared mutation status in matching primary tumor, liver metastasis and CTCs. CTCs were isolated from blood drawn from 49 patients before liver resection using CellSearch™. DNA and RNA was isolated from primary tumors, metastases and CTCs. Mutations were assessed by co-amplification at lower denaturation temperature-PCR (Transgenomic™), real-time PCR (EntroGen™) and nested Allele-Specific Blocker (ASB-)PCR and confirmed by Sanger sequencing. In 43 of the 49 patients, tissue RNA and DNA was of sufficient quantity and quality. In these 43 patients, discordance between primary and metastatic tumor was 23% for KRAS and 7% for BRAF mutations. RNA and DNA from CTCs was available from 42 of the 43 patients, in which ASB-PCR was able to detect the most mutations. Inconclusive results in patients with low CTC counts limited the interpretation of discrepancies between tissue and CTCs. Determination of KRAS and BRAF mutations in CTCs is challenging but feasible. Of the tested methods, nested ASB-PCR, enabling detection of KRAS and BRAF mutations in patients with as little as two CTCs, seems to be superior.

Hormone replacement therapy dependent changes in breast cancer-related gene expression in breast tissue of healthy postmenopausal women.

Risk assessment of future breast cancer risk through exposure to sex steroids currently relies on clinical scorings such as mammographic density. Knowledge about the gene expression patterns in existing breast cancer tumors may be used to identify risk factors in the breast tissue of women still free of cancer. The differential effects of estradiol, estradiol together with gestagens, or tibolone on breast cancer-related gene expression in normal breast tissue samples taken from postmenopausal women may be used to identify gene expression profiles associated with a higher breast cancer risk. Breast tissue samples were taken from 33 healthy postmenopausal women both before and after a six month treatment with either 2mg micronized estradiol [E2], 2mg micronized estradiol and 1mg norethisterone acetate [E2+NETA], 2.5mg tibolone [T] or [no HRT]. Except for [E2], which was only given to women after hysterectomy, the allocation to each of the three groups was randomized. The expression of 102 mRNAs and 46 microRNAs putatively involved in breast cancer was prospectively determined in the biopsies of 6 women receiving [no HRT], 5 women receiving [E2], 5 women receiving [E2+NETA], and 6 receiving [T]. Using epithelial and endothelial markers genes, non-representative biopsies from 11 women were eliminated. Treatment of postmenopausal women with [E2+NETA] resulted in the highest number of differentially (p<0.05) regulated genes (16.2%) compared to baseline, followed by [E2] (10.1%) and [T] (4.7%). Among genes that were significantly down-regulated by [E2+NETA] ranked estrogen-receptor-1 (ESR1, p=0.019) and androgen receptor (AR, p=0.019), whereas CYP1B1, a gene encoding an estrogen-metabolizing enzyme, was significantly up-regulated (p=0.016). Mammary cells triggered by [E2+NETA] and [E2] adjust for steroidogenic up-regulation through down-regulation of the estrogen-receptor pathway. In this prospective study, prolonged administration of [E2+NETA] and to a lesser extent of [E2] but not [T] were associated in otherwise healthy breast tissue with a change in the expression of genes putatively involved in breast cancer. Our data suggest that normal mammary cells triggered by [E2+NETA] adjust for steroidogenic up-regulation through down-regulation of the estrogen-receptor pathway. This feasibility study provides the basis for whole genome analyses to identify novel markers involved in increased breast cancer risk.

mRNA and microRNA expression profiles in circulating tumor cells and primary tumors of metastatic breast cancer patients.

PURPOSE: Molecular characterization of circulating tumor cells (CTC) holds great promise. Unfortunately, routinely isolated CTC fractions currently still contain contaminating leukocytes, which makes CTC-specific molecular characterization extremely challenging. In this study, we determined mRNA and microRNA (miRNA) expression of potentially CTC-specific genes that are considered to be clinically relevant in breast cancer. EXPERIMENTAL DESIGN: CTCs were isolated with the epithelial cell adhesion molecule-based CellSearch Profile Kit. Selected genes were measured by real-time reverse transcriptase PCR in CTCs of 50 metastatic breast cancer patients collected before starting first-line systemic therapy in blood from 53 healthy blood donors (HBD) and in primary tumors of 8 of the patients. The molecular profiles were associated with CTC counts and clinical parameters and compared with the profiles generated from the corresponding primary tumors. RESULTS: We identified 55 mRNAs and 10 miRNAs more abundantly expressed in samples from 32 patients with at least 5 CTCs in 7.5 mL of blood compared with samples from 9 patients without detectable CTCs and HBDs. Clustering analysis resulted in 4 different patient clusters characterized by 5 distinct gene clusters. Twice the number of patients from cluster 2 to 4 had developed both visceral and nonvisceral metastases. Comparing transcript levels in CTCs with those measured in corresponding primary tumors showed clinically relevant discrepancies in estrogen receptor and HER2 levels. CONCLUSIONS: Our study shows that molecular profiling of low numbers of CTCs in a high background of leukocytes is feasible and shows promise for further studies on the clinical relevance of molecular characterization of CTCs.

Clinical significance of the nuclear receptor co-regulator DC-SCRIPT in breast cancer: an independent retrospective validation study.

INTRODUCTION: In this study we aimed to validate the prognostic value of DC-SCRIPT mRNA expression in a large independent breast cancer cohort. In addition, since DC-SCRIPT is a transcriptional co-regulator of nuclear receptors, we explored its prognostic value in relation to estrogen-receptor-α (ESR1) and -ß (ESR2) and evaluated its predictive value for response to tamoxifen treatment. METHODS: DC-SCRIPT mRNA levels were measured by real-time PCR in 1,505 primary invasive breast cancers and associated with outcome (disease-free survival (DFS), metastasis-free survival (MFS) and overall survival (OS)) using univariate and multivariable Cox regression analysis. Logistic and Cox regressions were used to associate DC-SCRIPT levels with clinical benefit and progression-free survival (PFS) for 296 patients treated with first-line systemic tamoxifen for advanced disease. RESULTS: In univariate and multivariable analysis higher DC-SCRIPT levels were associated with a favorable outcome for both the entire cohort and patients with lymph node-negative (LNN) disease that did not receive adjuvant therapy (DFS, MFS and OS; all, P < 0.001). This association was most pronounced in small (pT1) tumors, in ESR1-positive tumors and in tumors with low ESR2 expression. For first-line endocrine therapy for advanced disease no predictive association was seen with clinical benefit or PFS. CONCLUSIONS: This study provides a higher level of evidence that DC-SCRIPT is indeed an independent, pure prognostic, factor for primary breast cancer and shows that DC-SCRIPT mRNA expression is most informative for either ESR1-positive and/or ESR2-low pT1 tumors.