Integrative analyses identify modulators of response to neoadjuvant aromatase inhibitors in patients with early breast cancer

Introduction Aromatase inhibitors (AIs) are a vital component of estrogen receptor positive (ER+) breast cancer treatment. De novo and acquired resistance, however, is common. The aims of this study were to relate patterns of copy number aberrations to molecular and proliferative response to AIs, to study differences in the patterns of copy number aberrations between breast cancer samples pre- and post-AI neoadjuvant therapy, and to identify putative biomarkers for resistance to neoadjuvant AI therapy using an integrative analysis approach. Methods Samples from 84 patients derived from two neoadjuvant AI therapy trials were subjected to copy number profiling by microarray-based comparative genomic hybridisation (aCGH, n = 84), gene expression profiling (n = 47), matched pre- and post-AI aCGH (n = 19 pairs) and Ki67-based AI-response analysis (n = 39). Results Integrative analysis of these datasets identified a set of nine genes that, when amplified, were associated with a poor response to AIs, and were significantly overexpressed when amplified, including CHKA, LRP5 and SAPS3. Functional validation in vitro, using cell lines with and without amplification of these genes (SUM44, MDA-MB134-VI, T47D and MCF7) and a model of acquired AI-resistance (MCF7-LTED) identified CHKA as a gene that when amplified modulates estrogen receptor (ER)-driven proliferation, ER/estrogen response element (ERE) transactivation, expression of ER-regulated genes and phosphorylation of V-AKT murine thymoma viral oncogene homolog 1 (AKT1). Conclusions These data provide a rationale for investigation of the role of CHKA in further models of de novo and acquired resistance to AIs, and provide proof of concept that integrative genomic analyses can identify biologically relevant modulators of AI response

Identification of chemokine receptors as potential modulators of endocrine resistance in oestrogen receptor–positive breast cancers

Introduction Endocrine therapies target oestrogenic stimulation of breast cancer (BC) growth, but resistance remains problematic. Our aims in this study were (1) to identify genes most strongly associated with resistance to endocrine therapy by intersecting global gene transcription data from patients treated presurgically with the aromatase inhibitor anastrazole with those from MCF7 cells adapted to long-term oestrogen deprivation (LTED) (2) to assess the clinical value of selected genes in public clinical data sets and (3) to determine the impact of targeting these genes with novel agents. Methods Gene expression and Ki67 data were available from 69 postmenopausal women with oestrogen receptor–positive (ER+) early BC, at baseline and 2 weeks after anastrazole treatment, and from cell lines adapted to LTED. The functional consequences of target genes on proliferation, ER-mediated transcription and downstream cell signalling were assessed. Results By intersecting genes predictive of a poor change in Ki67 with those upregulated in LTED cells, we identified 32 genes strongly correlated with poor antiproliferative response that were associated with inflammation and/or immunity. In a panel of LTED cell lines, C-X-C chemokine receptor type 7 (CXCR7) and CXCR4 were upregulated compared to their wild types (wt), and CXCR7, but not CXCR4, was associated with reduced relapse-free survival in patients with ER+ BC. The CXCR4 small interfering RNA variant (siCXCR4) had no specific effect on the proliferation of wt-SUM44, wt-MCF7 and their LTED derivatives. In contrast, siCXCR7, as well as CCX733, a CXCR7 antagonist, specifically suppressed the proliferation of MCF7-LTED cells. siCXCR7 suppressed proteins associated with G1/S transition and inhibited ER transactivation in MCF7-LTED, but not wt-MCF7, by impeding association between ER and proline-, glutamic acid– and leucine-rich protein 1, an ER coactivator. Conclusions These data highlight CXCR7 as a potential therapeutic target warranting clinical investigation in endocrine-resistant BC

ESR1 Is Co-Expressed with Closely Adjacent Uncharacterised Genes Spanning a Breast Cancer Susceptibility Locus at 6q25.1

Approximately 80% of human breast carcinomas present as oestrogen receptor α-positive (ER+ve) disease, and ER status is a critical factor in treatment decision-making. Recently, single nucleotide polymorphisms (SNPs) in the region immediately upstream of the ER gene (ESR1) on 6q25.1 have been associated with breast cancer risk. Our investigation of factors associated with the level of expression of ESR1 in ER+ve tumours has revealed unexpected associations between genes in this region and ESR1 expression that are important to consider in studies of the genetic causes of breast cancer risk. RNA from tumour biopsies taken from 104 postmenopausal women before and after 2 weeks treatment with an aromatase (oestrogen synthase) inhibitor was analyzed on Illumina 48K microarrays. Multiple-testing corrected Spearman correlation revealed that three previously uncharacterized open reading frames (ORFs) located immediately upstream of ESR1, C6ORF96, C6ORF97, and C6ORF211 were highly correlated with ESR1 (Rs = 0.67, 0.64, and 0.55 respectively, FDR<1×10−7). Publicly available datasets confirmed this relationship in other groups of ER+ve tumours. DNA copy number changes did not account for the correlations. The correlations were maintained in cultured cells. An ERα antagonist did not affect the ORFs' expression or their correlation with ESR1, suggesting their transcriptional co-activation is not directly mediated by ERα. siRNA inhibition of C6ORF211 suppressed proliferation in MCF7 cells, and C6ORF211 positively correlated with a proliferation metagene in tumours. In contrast, C6ORF97 expression correlated negatively with the metagene and predicted for improved disease-free survival in a tamoxifen-treated published dataset, independently of ESR1. Our observations suggest that some of the biological effects previously attributed to ER could be mediated and/or modified by these co-expressed genes. The co-expression and function of these genes may be important influences on the recently identified relationship between SNPs in this region and breast cancer risk

The Orphan Nuclear Receptor LRH-1 and ERα Activate GREB1 Expression to Induce Breast Cancer Cell Proliferation

BACKGROUND: Liver Receptor Homolog 1 (LRH-1, NR5A2) is an orphan nuclear receptor that is over-expressed in cancers in tissues such as the breast, colon and pancreas. LRH-1 plays important roles in embryonic development, steroidogenesis and cholesterol homeostasis. In tumor cells, LRH-1 induces proliferation and cell cycle progression. High LRH-1 expression is demonstrated in breast cancers, positively correlating with ERα status and aromatase activity. LRH-1 dependent cellular mechanisms in breast cancer epithelial cells are poorly defined. Hence in the present study we investigated the actions of LRH-1 in estrogen receptor α (ERα) positive breast cancer cells. RESULTS: The study aimed to investigate LRH-1 dependent mechanisms that promote breast cancer proliferation. We identified that LRH-1 regulated the expression of Growth Regulation by Estrogen in Breast Cancer 1 (GREB1) in MCF-7 and MDA-MB-231 cells. Over-expression of LRH-1 increased GREB1 mRNA levels while knockdown of LRH-1 reduced its expression. GREB1 is a well characterised ERα target gene, with three estrogen response elements (ERE) located on its promoter. Chromatin immunoprecipitation studies provided evidence of the co-localisation of LRH-1 and ERα at all three EREs. With electrophoretic mobility shift assays, we demonstrated direct binding of LRH-1 to EREs located on GREB1 and Trefoil Factor 1 (TFF1, pS2) promoters. LRH-1 and ERα co-operatively activated transcription of ERE luciferase reporter constructs suggesting an overlap in regulation of target genes in breast cancer cells. Over-expression of LRH-1 resulted in an increase in cell proliferation. This effect was more pronounced with estradiol treatment. In the presence of ICI 182,780, an ERα antagonist, LRH-1 still induced proliferation. CONCLUSIONS: We conclude that in ER-positive breast cancer cells, LRH-1 promotes cell proliferation by enhancing ERα mediated transcription of target genes such as GREB-1. Collectively these findings indicate the importance of LRH-1 in the progression of hormone-dependent breast cancer and implicate LRH-1 as a potential avenue for drug development

Expression levels of uridine 5'-diphospho-glucuronosyltransferase genes in breast tissue from healthy women are associated with mammographic density

Introduction Mammographic density (MD), as assessed from film screen mammograms, is determined by the relative content of adipose, connective and epithelial tissue in the female breast. In epidemiological studies, a high percentage of MD confers a four to six fold risk elevation of developing breast cancer, even after adjustment for other known breast cancer risk factors. However, the biologic correlates of density are little known. Methods Gene expression analysis using whole genome arrays was performed on breast biopsies from 143 women; 79 women with no malignancy (healthy women) and 64 newly diagnosed breast cancer patients, both included from mammographic centres. Percent MD was determined using a previously validated, computerized method on scanned mammograms. Significance analysis of microarrays (SAM) was performed to identify genes influencing MD and a linear regression model was used to assess the independent contribution from different variables to MD. Results SAM-analysis identified 24 genes differentially expressed between samples from breasts with high and low MD. These genes included three uridine 5'-diphospho-glucuronosyltransferase (UGT) genes and the oestrogen receptor gene (ESR1). These genes were down-regulated in samples with high MD compared to those with low MD. The UGT gene products, which are known to inactivate oestrogen metabolites, were also down-regulated in tumour samples compared to samples from healthy individuals. Several single nucleotide polymorphisms (SNPs) in the UGT genes associated with the expression of UGT and other genes in their vicinity were identified. Conclusions Three UGT enzymes were lower expressed both in breast tissue biopsies from healthy women with high MD and in biopsies from newly diagnosed breast cancers. The association was strongest amongst young women and women using hormonal therapy. UGT2B10 predicts MD independently of age, hormone therapy and parity. Our results indicate that down-regulation of UGT genes in women exposed to female sex hormones is associated with high MD and might increase the risk of breast cancer

Tumour sampling method can significantly influence gene expression profiles derived from neoadjuvant window studies

Patient-matched transcriptomic studies using tumour samples before and after treatment allow inter-patient heterogeneity to be controlled, but tend not to include an untreated comparison. Here, Illumina BeadArray technology was used to measure dynamic changes in gene expression from thirty-seven paired diagnostic core and surgically excised breast cancer biopsies obtained from women receiving no treatment prior to surgery, to determine the impact of sampling method and tumour heterogeneity. Despite a lack of treatment and perhaps surprisingly, consistent changes in gene expression were identified during the diagnosis-surgery interval (48 up, 2 down; Siggenes FDR 0.05) in a manner independent of both subtype and sampling-interval length. Instead, tumour sampling method was seen to directly impact gene expression, with similar effects additionally identified in six published breast cancer datasets. In contrast with previous findings, our data does not support the concept of a significant wounding or immune response following biopsy in the absence of treatment and instead implicates a hypoxic response following the surgical biopsy. Whilst sampling-related gene expression changes are evident in treated samples, they are secondary to those associated with response to treatment. Nonetheless, sampling method remains a potential confounding factor for neoadjuvant study design

Genome-Wide Association Study in East Asians Identifies Novel Susceptibility Loci for Breast Cancer

Genetic factors play an important role in the etiology of both sporadic and familial breast cancer. We aimed to discover novel genetic susceptibility loci for breast cancer. We conducted a four-stage genome-wide association study (GWAS) in 19,091 cases and 20,606 controls of East-Asian descent including Chinese, Korean, and Japanese women. After analyzing 690,947 SNPs in 2,918 cases and 2,324 controls, we evaluated 5,365 SNPs for replication in 3,972 cases and 3,852 controls. Ninety-four SNPs were further evaluated in 5,203 cases and 5,138 controls, and finally the top 22 SNPs were investigated in up to 17,423 additional subjects (7,489 cases and 9,934 controls). SNP rs9485372, near the TGF-β activated kinase (TAB2) gene in chromosome 6q25.1, showed a consistent association with breast cancer risk across all four stages, with a P-value of 3.8×10−12 in the combined analysis of all samples. Adjusted odds ratios (95% confidence intervals) were 0.89 (0.85–0.94) and 0.80 (0.75–0.86) for the A/G and A/A genotypes, respectively, compared with the genotype G/G. SNP rs9383951 (P = 1.9×10−6 from the combined analysis of all samples), located in intron 5 of the ESR1 gene, and SNP rs7107217 (P = 4.6×10−7), located at 11q24.3, also showed a consistent association in each of the four stages. This study provides strong evidence for a novel breast cancer susceptibility locus represented by rs9485372, near the TAB2 gene (6q25.1), and identifies two possible susceptibility loci located in the ESR1 gene and 11q24.3, respectively