Genomic and transcriptomic alterations that affect cellular processes, such as cell proliferation, differentiation, apoptosis and invasion, commonly occur in breast oncogenesis. Epidemiological evidence has proven that the risk of breast cancer predisposition varies among different ethnicities. This study aims to identify the transcriptome changes that commonly occur during the transition of normal breast epithelium to carcinoma in three local ethnic groups (Malays, Chinese and Indians). The gene expression patterns of 43 breast carcinomas with 43 patient-matched normal breast tissues were investigated using Affymetrix U133A GeneChip (containing 22,283 probe sets targeting approximately 18,400 different transcripts) and analyzed with GeneSpring GX10. Our findings revealed a total of 33 significantly differentially expressed genes, which showed>2-fold change at a 99.9% confidence interval level (p<0.001). The significantly differentially expressed genes included CD24, CD36, CD9, TACSTD1, TACSTD2, HBB, LEP, LPL, AKR1C1, AKR1C2 and AKR1C3. Our results indicate that the vast majority of gene expression changes, from normal breast epithelial to carcinoma, found in our three major ethnic populations are similar to those in the Caucasian population. Further study of the differentially expressed genes identified in our present study is needed to search for potential breast tumor biomarkers. This will eventually help to improve the therapeutic and treatment strategies for breast cancer patients in the future.
MicroRNAs (miRNAs) and long non-coding RNAs (lncRNAs) are capable of regulating gene expression post-transcriptionally. Since the past decade, a number of in vitro, in vivo, and clinical studies reported the roles of these non-coding RNAs (ncRNAs) in regulating angiogenesis, an important cancer hallmark that is associated with metastases and poor prognosis. The specific roles of various miRNAs and lncRNAs in regulating angiogenesis in breast cancer, with particular focus on the downstream targets and signalling pathways regulated by these ncRNAs will be discussed in this review. In light of the recent trend in exploiting ncRNAs as cancer therapeutics, the potential use of miRNAs and lncRNAs as biomarkers and novel therapeutic agent against angiogenesis was also discussed.
Activation of Akt signaling pathway has been documented in various human malignancies, including breast carcinoma. The objective of this study is to determine the incidence of Akt phosphorylation in breast tumours and its relationship with expression of ER-alpha, ER-beta, HER2, Ki-67 and phosphorylated Bcl-2 associated death domain (p-BAD). Immunohistochemical staining was performed to detect these molecules on 43 paraffin-embedded breast tumour tissues with commercially available antibodies. Eighteen (41.9%), 3 (7.0%), 23 (53.5%), 35 (81.4%), 21 (48.8%), 29 (67.4%), and 34 (81.0%) of breast tumours were positive for nuclear ER-alpha, nuclear ER-beta, membranous HER2, cytonuclear p-Akt (Thr308), p-Akt (Ser473), p-BAD and Ki-67, respectively. ER-alpha expression was inversely correlated with HER2 and Ki-67 (P = 0.041 and P = 0.040, respectively). The p-Akt (Ser473) was correlated with increased level of p-BAD (Ser136) (P = 0.012). No relationship of Akt phosphorylation with HER2, ER-alpha or ER-beta was found. The p-Akt (Ser473) immunoreactivity was significantly higher in stage IV than in stage I or II (P = 0.036 or P = 0.009). The higher Ki-67 and lower ER-alpha expression showed an association with patient age of <50 years (P = 0.004) and with positive nodal status (P = 0.033), respectively. Our data suggest that the Akt phosphorylation and inactivation of its downstream target, BAD may play a role in survival of breast cancer cell. This study does not support the simple model of linear HER2/PI3K/Akt pathway in breast cancer.
The aim of the present study was to evaluate the association between the Glyoxalase I (GLOI) Ala111Glu polymorphism and breast cancer risk among the major Malaysian ethnic groups, the Malays, Chinese and Indians, as well as clinico-pathological characteristics of these patients. Genotyping of GLOI gene was performed on blood samples obtained from 387 patients and 252 normal healthy women who had no history of any malignancy using the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method. The genotype and allele frequencies of GLOI polymorphism were not significantly different between the patients and normal individuals among the Malays (P= 0.721, 0.402), Chinese (P= 0.208, 0.079) and Indians (P= 0.612, 0.349), respectively. The Malay, Chinese and Indian women who were Glu/Glu homozygotes (P= 0.419, 0.093, 0.367), Ala/Glu heterozygotes (P= 0.648, 0.182, 0.402) and carriers of Glu allele (P= 0.402, 0.079, 0.349), respectively, were not associated with breast cancer risk. The Glu allele genotype was significantly associated with absence of progesterone receptor (P= 0.036). Thus, the polymorphic variant of the GLOI gene might not be a useful genetic marker to identify Malaysian Malay, Chinese or Indian women who could be at greater risk of developing breast cancer.
Study site: Universiti Malaya Medical Centre (UMMC)
The present study aimed to evaluate the association between the E-Selectin Ser128Arg polymorphism and breast cancer risk and clinicopathological characteristics of the patients.
Epithelial to mesenchymal transition (EMT) is a developmental program that has been implicated in progression, metastasis and therapeutic resistance of some carcinomas. To identify genes whose overexpression drives EMT, we screened a lentiviral expression library of 17000 human open reading frames (ORFs) using high-content imaging to quantitate cytoplasmic vimentin. Hits capable of increasing vimentin in the mammary carcinoma-derived cell line MDA-MB-468 were confirmed in the non-tumorigenic breast-epithelial cell line MCF10A. When overexpressed in this model, they increased the rate of cell invasion through Matrigel™, induced mesenchymal marker expression and reduced expression of the epithelial marker E-cadherin. In gene-expression datasets derived from breast cancer patients, the expression of several novel genes correlated with expression of known EMT marker genes, indicating their in vivo relevance. As EMT-associated properties are thought to contribute in several ways to cancer progression, genes identified in this study may represent novel targets for anti-cancer therapy.
The genotype analysis of the Gly and Arg allele at codon 388 of fibroblast growth factor receptor-4 (FGFR4) gene was evaluated using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method in a hospital-based Malaysian population. Peripheral blood samples were collected from 387 breast cancer patients and 252 normal and healthy women who had no history of any malignancy. The aim of the present study was to evaluate the association between the FGFR4 Gly388Arg polymorphism and breast cancer risk as well as clinicopathological parameters of the patients. The Gly/Gly, Gly/Arg, Arg/Arg, and Arg allele genotypes were detected in 46.3%, 44.4%, 9.3%, and 53.7% of breast cancer cases, respectively. The distribution of genotype (p = 0.204) and allele (p = 0.086) frequencies of FGFR4 polymorphism were not significantly different between the breast cancer cases and normal individuals. Women who were Arg/ Arg homozygotes (OR = 1.714, 95% CI 0.896-3.278), Gly/Arg heterozygotes (OR = 1.205, 95% CI 0.863-1.683), carriers of Arg allele genotype (OR = 1.269, 95% CI 0.921-1.750), or Arg allele (OR = 1.246, 95% CI 0.970-1.602) were not associated with breast cancer risk. The Arg allele genotype was significantly associated with lymph node metastases (p = 0.001) but not with other clinicopathological parameters. Our findings suggest that the polymorphic variant at codon 388 of FGFR4 gene does not confer increased risk to breast cancer development but it may be a potential genetic marker for tumor prognosis.
Overexpression and amplification of cyclin D1 were investigated by immunohistochemistry and differential polymerase chain reaction (dPCR) in 440 formalin-fixed primary breast carcinoma tissues. Overexpression of cyclin D1 was detected in 60% (263/440) and amplification of cyclin D1 was noted in 27% (119/440) of the primary breast carcinomas. Molecular analysis demonstrated that cyclin D1 was amplified in 30% (7/23) of the comedo DCIS, 22% (9/41) of the comedo DCIS and 32% (13/41) of the adjacent invasive ductal carcinomas, 30% (82/270) of the invasive ductal carcinomas, 27% (9/33) of the invasive lobular carcinomas, 19% (4/21) of the colloid carcinomas and 13% (2/15) of the medullary carcinomas. Cyclin D1 was amplified in 11% (2/19) of the invasive ductal carcinomas but not in the adjacent non-comedo DCIS lesions. Our observation showed that cyclin D1 was strongly positive in 61% (14/23) of the comedo subtype, 61% (11/18) of the non-comedo subtype, 59% (24/41) of the comedo DCIS and 63% (26/41) of the adjacent invasive ductal carcinomas, 53% (10/19) of the non-comedo DCIS and 58% (11/19) of the adjacent invasive lesions, 58% (157/270) of the invasive ductal carcinomas, 73% (24/33) of the invasive lobular carcinomas, 52% (11/21) of the colloid carcinomas and 27% (4/15) of the medullary carcinomas. A significant association was observed between in situ components and adjacent invasive lesions for cyclin D1 expression (p<0.05) and amplification (p<0.05). A significant relationship was noted between amplification of cyclin D1 and lymph node metastases (p<0.05) but not with histological grade (p>0.05), estrogen receptor status (p>0.05) and proliferation index (Ki-67 and PCNA) (p>0.05). However, overexpression of cyclin D1 was statistically associated with well differentiated tumors (p<0.05) and estrogen receptor positivity (p<0.05). No relationship was seen with nodal status (p>0.05) and proliferation index (Ki-67 and PCNA) (p>0.05). These observations suggest that tumors positive for cyclin D1 protein may have features of good prognosis but amplification of cyclin D1 gene could be an indicator of tumors with poor prognostic features. Although majority of the Malaysian patients belong to younger age group (<50 years old), amplification and expression of cyclin D1 was not statistically associated with patient age (p>0.05). These observations indicate that amplification and up-regulation of cyclin D1 may be independent of patient age. Moreover, overexpression and amplification of cyclin D1 in preinvasive, preinvasive and adjacent invasive lesions, and invasive carcinomas suggest that the gene may play an important role in early and late stages of breast carcinogenesis.
Kinase suppressor of Ras-1 (KSR1) facilitates signal transduction in Ras-dependent cancers, including pancreatic and lung carcinomas but its role in breast cancer has not been well studied. Here, we demonstrate for the first time it functions as a tumor suppressor in breast cancer in contrast to data in other tumors. Breast cancer patients (n>1000) with high KSR1 showed better disease-free and overall survival, results also supported by Oncomine analyses, microarray data (n=2878) and genomic data from paired tumor and cell-free DNA samples revealing loss of heterozygosity. KSR1 expression is associated with high breast cancer 1, early onset (BRCA1), high BRCA1-associated ring domain 1 (BARD1) and checkpoint kinase 1 (Chk1) levels. Phospho-profiling of major components of the canonical Ras-RAF-mitogen-activated protein kinases pathway showed no significant changes after KSR1 overexpression or silencing. Moreover, KSR1 stably transfected cells formed fewer and smaller size colonies compared to the parental ones, while in vivo mouse model also demonstrated that the growth of xenograft tumors overexpressing KSR1 was inhibited. The tumor suppressive action of KSR1 is BRCA1 dependent shown by 3D-matrigel and soft agar assays. KSR1 stabilizes BRCA1 protein levels by reducing BRCA1 ubiquitination through increasing BARD1 abundance. These data link these proteins in a continuum with clinical relevance and position KSR1 in the major oncoprotein pathways in breast tumorigenesis.
The C1772T, G1790A and C111A polymorphisms of Hypoxia-inducible factor-1alpha (HIF-1alpha) gene were analyzed in a hospital-based Malaysian population using PCR-RFLP method. Genomic DNA was extracted from the blood samples collected from 410 breast cancer patients and 275 normal and healthy women. We investigated the association between HIF-1alpha polymorphisms and breast cancer risk, and clinico-pathological parameters in the population. The genotype and allele frequencies of C1772T (P=0.0093 vs P=0.0024) polymorphism were significantly different between the breast cancer cases and normal subjects but similar association was not observed for G1790A (P>0.05) and C111A (P>0.05) polymorphisms, respectively. Women who were CT heterozygotes (OR=1.51; 95% CI, 1.01-2.25), TT homozygotes (OR=4.03; 95% CI, 1.09-17.60) and carriers of T allele genotype (OR=1.65; 95% CI, 1.13-2.43) were significantly associated with increased risk of breast cancer. Significant relationship was observed also between T allele and breast cancer risk (OR=1.69; 95% CI, 1.20-2.40). Clinico-pathological analysis showed that 1772T allele genotype was significantly associated with nodal metastases (P=0.0478) but independent of ER status, tumor grade and patients' age (P>0.05). Our observations suggest that the polymorphic allele of C1772T may be associated with increased risk of developing breast cancer, and presence of 1772T allele may be a useful genetic marker for tumor prognosis.
The HER2 codon Ile655Val and Cyclin D1 (CCND1) G870A polymorphisms were analyzed in a hospital-based Malaysian population using PCR-RFLP method. Peripheral blood samples were collected from 230 breast cancer patients, and 200 normal and healthy women who had no history of breast disease or breast cancer. We evaluated the association between HER2 or CCND1 polymorphisms and breast cancer risk, and clinico-pathological parameters in the population. The genotype and allele frequencies of HER2 (P=0.163 vs P=0.0622) and CCND1 (P=0.377 vs P=0.284) polymorphisms were not significantly different between the breast cancer cases and normal subjects, respectively. Women who were Ile/Val heterozygotes (OR=1.48; 95% CI, 0.91-2.43), Val/Val homozygotes (OR=1.93; 95% CI, 0.51-7.77) and carriers of Val allele genotype (OR=1.53; 95% CI, 0.95-2.45) were not significantly associated with increased breast cancer risk. Similarly, women who were homozygous (OR=1.34; 95% CI, 0.77-2.34) or heterozygous (OR=0.98; 95% CI, 0.60-1.60) for A allele, or carriers of A allele genotype (OR=1.10; 95% CI, 0.70-1.73) were not associated with breast cancer risk. Analysis on clinico-pathological parameters showed that Val allele genotype was significantly correlated with nodal metastases but A allele genotype was not associated with any of the variables. Our findings suggest that the polymorphic alleles of HER2 and CCND1 may not play an important role as genetic markers for breast cancer risk, but presence of Val allele may be useful for tumor prognosis.
Most common breast cancer susceptibility variants have been identified through genome-wide association studies (GWAS) of predominantly estrogen receptor (ER)-positive disease. We conducted a GWAS using 21,468 ER-negative cases and 100,594 controls combined with 18,908 BRCA1 mutation carriers (9,414 with breast cancer), all of European origin. We identified independent associations at P < 5 × 10-8 with ten variants at nine new loci. At P < 0.05, we replicated associations with 10 of 11 variants previously reported in ER-negative disease or BRCA1 mutation carrier GWAS and observed consistent associations with ER-negative disease for 105 susceptibility variants identified by other studies. These 125 variants explain approximately 16% of the familial risk of this breast cancer subtype. There was high genetic correlation (0.72) between risk of ER-negative breast cancer and breast cancer risk for BRCA1 mutation carriers. These findings may lead to improved risk prediction and inform further fine-mapping and functional work to better understand the biological basis of ER-negative breast cancer.
Genome-wide association studies have identified breast cancer risk variants in over 150 genomic regions, but the mechanisms underlying risk remain largely unknown. These regions were explored by combining association analysis with in silico genomic feature annotations. We defined 205 independent risk-associated signals with the set of credible causal variants in each one. In parallel, we used a Bayesian approach (PAINTOR) that combines genetic association, linkage disequilibrium and enriched genomic features to determine variants with high posterior probabilities of being causal. Potentially causal variants were significantly over-represented in active gene regulatory regions and transcription factor binding sites. We applied our INQUSIT pipeline for prioritizing genes as targets of those potentially causal variants, using gene expression (expression quantitative trait loci), chromatin interaction and functional annotations. Known cancer drivers, transcription factors and genes in the developmental, apoptosis, immune system and DNA integrity checkpoint gene ontology pathways were over-represented among the highest-confidence target genes.
Breast cancer susceptibility variants frequently show heterogeneity in associations by tumor subtype1-3. To identify novel loci, we performed a genome-wide association study including 133,384 breast cancer cases and 113,789 controls, plus 18,908 BRCA1 mutation carriers (9,414 with breast cancer) of European ancestry, using both standard and novel methodologies that account for underlying tumor heterogeneity by estrogen receptor, progesterone receptor and human epidermal growth factor receptor 2 status and tumor grade. We identified 32 novel susceptibility loci (P breast cancer subtypes and will inform the development of subtype-specific polygenic risk scores.
Matched MeSH terms: Breast Neoplasms/genetics*; Triple Negative Breast Neoplasms/genetics
In a three-stage genome-wide association study among East Asian women including 22,780 cases and 24,181 controls, we identified 3 genetic loci newly associated with breast cancer risk, including rs4951011 at 1q32.1 (in intron 2 of the ZC3H11A gene; P=8.82×10(-9)), rs10474352 at 5q14.3 (near the ARRDC3 gene; P=1.67×10(-9)) and rs2290203 at 15q26.1 (in intron 14 of the PRC1 gene; P=4.25×10(-8)). We replicated these associations in 16,003 cases and 41,335 controls of European ancestry (P=0.030, 0.004 and 0.010, respectively). Data from the ENCODE Project suggest that variants rs4951011 and rs10474352 might be located in an enhancer region and transcription factor binding sites, respectively. This study provides additional insights into the genetics and biology of breast cancer.
We analyzed 3,872 common genetic variants across the ESR1 locus (encoding estrogen receptor α) in 118,816 subjects from three international consortia. We found evidence for at least five independent causal variants, each associated with different phenotype sets, including estrogen receptor (ER(+) or ER(-)) and human ERBB2 (HER2(+) or HER2(-)) tumor subtypes, mammographic density and tumor grade. The best candidate causal variants for ER(-) tumors lie in four separate enhancer elements, and their risk alleles reduce expression of ESR1, RMND1 and CCDC170, whereas the risk alleles of the strongest candidates for the remaining independent causal variant disrupt a silencer element and putatively increase ESR1 and RMND1 expression.
Genome-wide association studies (GWAS) and large-scale replication studies have identified common variants in 79 loci associated with breast cancer, explaining ∼14% of the familial risk of the disease. To identify new susceptibility loci, we performed a meta-analysis of 11 GWAS, comprising 15,748 breast cancer cases and 18,084 controls together with 46,785 cases and 42,892 controls from 41 studies genotyped on a 211,155-marker custom array (iCOGS). Analyses were restricted to women of European ancestry. We generated genotypes for more than 11 million SNPs by imputation using the 1000 Genomes Project reference panel, and we identified 15 new loci associated with breast cancer at P < 5 × 10(-8). Combining association analysis with ChIP-seq chromatin binding data in mammary cell lines and ChIA-PET chromatin interaction data from ENCODE, we identified likely target genes in two regions: SETBP1 at 18q12.3 and RNF115 and PDZK1 at 1q21.1. One association appears to be driven by an amino acid substitution encoded in EXO1.
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.
Polygenic risk scores (PRS) have been shown to predict breast cancer risk in European women, but their utility in Asian women is unclear. Here we evaluate the best performing PRSs for European-ancestry women using data from 17,262 breast cancer cases and 17,695 controls of Asian ancestry from 13 case-control studies, and 10,255 Chinese women from a prospective cohort (413 incident breast cancers). Compared to women in the middle quintile of the risk distribution, women in the highest 1% of PRS distribution have a ~2.7-fold risk and women in the lowest 1% of PRS distribution has ~0.4-fold risk of developing breast cancer. There is no evidence of heterogeneity in PRS performance in Chinese, Malay and Indian women. A PRS developed for European-ancestry women is also predictive of breast cancer risk in Asian women and can help in developing risk-stratified screening programmes in Asia.
Molecular profiling of breast cancer has enabled the development of more robust molecular prognostic signatures and therapeutic options for breast cancer patients. However, non-Caucasian populations remain understudied. Here, we present the mutational, transcriptional, and copy number profiles of 560 Malaysian breast tumours and a comparative analysis of breast cancers arising in Asian and Caucasian women. Compared to breast tumours in Caucasian women, we show an increased prevalence of HER2-enriched molecular subtypes and higher prevalence of TP53 somatic mutations in ER+ Asian breast tumours. We also observe elevated immune scores in Asian breast tumours, suggesting potential clinical response to immune checkpoint inhibitors. Whilst HER2-subtype and enriched immune score are associated with improved survival, presence of TP53 somatic mutations is associated with poorer survival in ER+ tumours. Taken together, these population differences unveil opportunities to improve the understanding of this disease and lay the foundation for precision medicine in different populations.