Assessment of variation in immunosuppressive pathway genes reveals TGFBR2 to be associated with prognosis of estrogen receptor-negative breast cancer after chemotherapy.

INTRODUCTION: Tumor lymphocyte infiltration is associated with clinical response to chemotherapy in estrogen receptor (ER) negative breast cancer. To identify variants in immunosuppressive pathway genes associated with prognosis after adjuvant chemotherapy for ER-negative patients, we studied stage I-III invasive breast cancer patients of European ancestry, including 9,334 ER-positive (3,151 treated with chemotherapy) and 2,334 ER-negative patients (1,499 treated with chemotherapy). METHODS: We pooled data from sixteen studies from the Breast Cancer Association Consortium (BCAC), and employed two independent studies for replications. Overall 3,610 single nucleotide polymorphisms (SNPs) in 133 genes were genotyped as part of the Collaborative Oncological Gene-environment Study, in which phenotype and clinical data were collected and harmonized. Multivariable Cox proportional hazard regression was used to assess genetic associations with overall survival (OS) and breast cancer-specific survival (BCSS). Heterogeneity according to chemotherapy or ER status was evaluated with the log-likelihood ratio test. RESULTS: Three independent SNPs in TGFBR2 and IL12B were associated with OS (P <10⁻³) solely in ER-negative patients after chemotherapy (267 events). Poorer OS associated with TGFBR2 rs1367610 (G > C) (per allele hazard ratio (HR) 1.54 (95% confidence interval (CI) 1.22 to 1.95), P = 3.08 × 10⁻4) was not found in ER-negative patients without chemotherapy or ER-positive patients with chemotherapy (P for interaction <10-3). Two SNPs in IL12B (r² = 0.20) showed different associations with ER-negative disease after chemotherapy: rs2546892 (G > A) with poorer OS (HR 1.50 (95% CI 1.21 to 1.86), P = 1.81 × 10⁻4), and rs2853694 (A > C) with improved OS (HR 0.73 (95% CI 0.61 to 0.87), P = 3.67 × 10⁻4). Similar associations were observed with BCSS. Association with TGFBR2 rs1367610 but not IL12B variants replicated using BCAC Asian samples and the independent Prospective Study of Outcomes in Sporadic versus Hereditary Breast Cancer Study and yielded a combined HR of 1.57 ((95% CI 1.28 to 1.94), P = 2.05 × 10⁻5) without study heterogeneity. CONCLUSIONS: TGFBR2 variants may have prognostic and predictive value in ER-negative breast cancer patients treated with adjuvant chemotherapy. Our findings provide further insights into the development of immunotherapeutic targets for ER-negative breast cancer.

Expression of an estrogen-regulated variant transcript of the peroxisomal branched chain fatty acid oxidase ACOX2 in breast carcinomas.

BACKGROUND: Alternate transcripts from a single gene locus greatly enhance the combinatorial flexibility of the human transcriptome. Different patterns of exon usage have been observed when comparing normal tissue to cancers, suggesting that variant transcripts may play a role in the tumor phenotype. METHODS: Ribonucleic acid-sequencing (RNA-seq) data from breast cancer samples was used to identify an intronic start variant transcript of Acyl-CoA oxidase 2, ACOX2 (ACOX2-i9). Difference in expression between Estrogen Receptor (ER) positive and ER negative patients was assessed by the Wilcoxon rank sum test, and the findings validated in The Cancer Genome Atlas (TCGA) breast cancer dataset (BRCA). ACOX2-i9 expression was also assessed in cell lines using both quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) and Western blot analysis. Knock down by short hairpin RNA (shRNA) and colony formation assays were used to determine whether ACOX2-i9 expression would influence cellular fitness. The effect of ACOX2-i9 expression on patient survival was assessed by the Kaplan-Meier survival function, and association to clinical parameters was analyzed using a Fisher exact test. RESULTS: The expression and translation of ACOX2-i9 into a 25 kDa protein was demonstrated in HepG2 cells as well as in several breast cancer cell lines. shRNA knock down of the ACOX2-i9 variant resulted in decreased cell viability of T47D and MDA-MB 436 cells. Moreover, expression of ACOX2-i9 was shown to be estrogen regulated, being induced by propyl pyrazoletriol and inhibited by tamoxifen and fulvestrant in ER+ T47D and Mcf-7 cells, but not in the ER- MDA-MB 436 cell line. This variant transcript showed expression predominantly in ER-positive breast tumors as assessed in our initial set of 53 breast cancers and further validated in 87 tumor/normal pairs from the TCGA breast cancer dataset, and expression was associated with better outcome in ER positive patients. CONCLUSIONS: ACOX2-i9 is specifically enriched in ER+ breast cancers where expression of the variant is associated with improved outcome. These data identify variant ACOX2 as a potential novel therapeutic biomarker in ER+ breast tumors.

Widespread alternative exon usage in clinically distinct subtypes of Invasive Ductal Carcinoma.

Cancer cells can have different patterns of exon usage of individual genes when compared to normal tissue, suggesting that alternative splicing may play a role in shaping the tumor phenotype. The discovery and identification of gene variants has increased dramatically with the introduction of RNA-sequencing technology, which enables whole transcriptome analysis of known, as well as novel isoforms. Here we report alternative splicing and transcriptional events among subtypes of invasive ductal carcinoma in The Cancer Genome Atlas (TCGA) Breast Invasive Carcinoma (BRCA) cohort. Alternative exon usage was widespread, and although common events were shared among three subtypes, ER+ HER2-, ER- HER2-, and HER2+, many events on the exon level were subtype specific. Additional RNA-seq analysis was carried out in an independent cohort of 43 ER+ HER2- and ER- HER2- primary breast tumors, confirming many of the exon events identified in the TCGA cohort. Alternative splicing and transcriptional events detected in five genes, MYO6, EPB41L1, TPD52, IQCG, and ACOX2 were validated by qRT-PCR in a third cohort of 40 ER+ HER2- and ER- HER2- patients, showing that these events were truly subtype specific.

Heterogeneous DNA Methylation Patterns in the GSTP1 Promoter Lead to Discordant Results between Assay Technologies and Impede Its Implementation as Epigenetic Biomarkers in Breast Cancer.

Altered DNA methylation patterns are found in many diseases, particularly in cancer, where the analysis of DNA methylation holds the promise to provide diagnostic, prognostic and predictive information of great clinical value. Methylation of the promoter-associated CpG island of GSTP1 occurs in many hormone-sensitive cancers, has been shown to be a biomarker for the early detection of cancerous lesions and has been associated with important clinical parameters, such as survival and response to treatment. In the current manuscript, we assessed the performance of several widely-used sodium bisulfite conversion-dependent methods (methylation-specific PCR, MethyLight, pyrosequencing and MALDI mass-spectrometry) for the analysis of DNA methylation patterns in the GSTP1 promoter. We observed large discordances between the results obtained by the different technologies. Cloning and sequencing of the investigated region resolved single-molecule DNA methylation patterns and identified heterogeneous DNA methylation patterns as the underlying cause of the differences. Heterogeneous DNA methylation patterns in the GSTP1 promoter constitute a major obstacle to the implementation of DNA methylation-based analysis of GSTP1 and might explain some of the contradictory findings in the analysis of the significance of GSTP1 promoter methylation in breast cancer.

Novel Associations between Common Breast Cancer Susceptibility Variants and Risk-Predicting Mammographic Density Measures.

Mammographic density measures adjusted for age and body mass index (BMI) are heritable predictors of breast cancer risk, but few mammographic density-associated genetic variants have been identified. Using data for 10,727 women from two international consortia, we estimated associations between 77 common breast cancer susceptibility variants and absolute dense area, percent dense area and absolute nondense area adjusted for study, age, and BMI using mixed linear modeling. We found strong support for established associations between rs10995190 (in the region of ZNF365), rs2046210 (ESR1), and rs3817198 (LSP1) and adjusted absolute and percent dense areas (all P < 10(-5)). Of 41 recently discovered breast cancer susceptibility variants, associations were found between rs1432679 (EBF1), rs17817449 (MIR1972-2: FTO), rs12710696 (2p24.1), and rs3757318 (ESR1) and adjusted absolute and percent dense areas, respectively. There were associations between rs6001930 (MKL1) and both adjusted absolute dense and nondense areas, and between rs17356907 (NTN4) and adjusted absolute nondense area. Trends in all but two associations were consistent with those for breast cancer risk. Results suggested that 18% of breast cancer susceptibility variants were associated with at least one mammographic density measure. Genetic variants at multiple loci were associated with both breast cancer risk and the mammographic density measures. Further understanding of the underlying mechanisms at these loci could help identify etiologic pathways implicated in how mammographic density predicts breast cancer risk.

2q36.3 is associated with prognosis for oestrogen receptor-negative breast cancer patients treated with chemotherapy.

Large population-based registry studies have shown that breast cancer prognosis is inherited. Here we analyse single-nucleotide polymorphisms (SNPs) of genes implicated in human immunology and inflammation as candidates for prognostic markers of breast cancer survival involving 1,804 oestrogen receptor (ER)-negative patients treated with chemotherapy (279 events) from 14 European studies in a prior large-scale genotyping experiment, which is part of the Collaborative Oncological Gene-environment Study (COGS) initiative. We carry out replication using Asian COGS samples (n=522, 53 events) and the Prospective Study of Outcomes in Sporadic versus Hereditary breast cancer (POSH) study (n=315, 108 events). Rs4458204_A near CCL20 (2q36.3) is found to be associated with breast cancer-specific death at a genome-wide significant level (n=2,641, 440 events, combined allelic hazard ratio (HR)=1.81 (1.49-2.19); P for trend=1.90 × 10(-9)). Such survival-associated variants can represent ideal targets for tailored therapeutics, and may also enhance our current prognostic prediction capabilities.

Individual and combined effects of DNA methylation and copy number alterations on miRNA expression in breast tumors.

BACKGROUND: The global effect of copy number and epigenetic alterations on miRNA expression in cancer is poorly understood. In the present study, we integrate genome-wide DNA methylation, copy number and miRNA expression and identify genetic mechanisms underlying miRNA dysregulation in breast cancer. RESULTS: We identify 70 miRNAs whose expression was associated with alterations in copy number or methylation, or both. Among these, five miRNA families are represented. Interestingly, the members of these families are encoded on different chromosomes and are complementarily altered by gain or hypomethylation across the patients. In an independent breast cancer cohort of 123 patients, 41 of the 70 miRNAs were confirmed with respect to aberration pattern and association to expression. In vitro functional experiments were performed in breast cancer cell lines with miRNA mimics to evaluate the phenotype of the replicated miRNAs. let-7e-3p, which in tumors is found associated with hypermethylation, is shown to induce apoptosis and reduce cell viability, and low let-7e-3p expression is associated with poorer prognosis. The overexpression of three other miRNAs associated with copy number gain, miR-21-3p, miR-148b-3p and miR-151a-5p, increases proliferation of breast cancer cell lines. In addition, miR-151a-5p enhances the levels of phosphorylated AKT protein. CONCLUSIONS: Our data provide novel evidence of the mechanisms behind miRNA dysregulation in breast cancer. The study contributes to the understanding of how methylation and copy number alterations influence miRNA expression, emphasizing miRNA functionality through redundant encoding, and suggests novel miRNAs important in breast cancer.