A genome-wide association study to identify genetic susceptibility loci that modify ductal and lobular postmenopausal breast cancer risk associated with menopausal hormone therapy use: a two-stage design with replication

Potential Role of Zinc Finger 365 rs10822013 and rs10995190 in Mammographic Density, Sporadic Breast Cancer Risk, and Prognosis

Background

Despite suggesting many genetic risk markers as the outcome of Genome-wide association studies (GWAS) for breast cancer, replicating the results in different populations has remained the main issue. In this regard, this study assessed the association of two variations in Zinc Finger 365 (ZNF365) in an Iranian population.

Methods

In a case-control study conducted at Mashhad University of Medical Sciences, Mashhad, Iran, between 2017 and 2020, ZNF365-rs10822013 and rs10995190 were genotyped using Allele-Specific PCR (AS-PCR). Breast density was assessed using mammography images. PHASE software module version 2 and SPSS version 16.0 were used for haplotype and statistical analyses. Quantitative and qualitative variables were compared between groups using independent t tests and Chi square tests, respectively. Binary logistic regression analysis was performed to calculate odds ratios. Multivariate analysis was then undertaken for the baseline variables, with a P<0.05 in the univariate analysis. The survival analysis was performed using the Kaplan-Meier method and the log-rank test.

Results

In this survey, 732 females, including 342 breast cancer patients and 390 healthy subjects, were enrolled. rs10822013-T allele (P=0.014), rs10995190-G allele (P=0.003), and TG haplotype (P=0.002) were significantly associated with the increased risk of breast cancer. Moreover, rs10995190-GG genotype (P=0.042) and C-G haplotype (P=0.019) revealed a significant association with better overall survival. However, considered polymorphisms and their haplotypes indicated no association with breast density and clinical features of breast cancer.

Conclusion

ZNF365 variants might be a potential risk marker of breast cancer in the Iranian population. The interaction between alleles in haplotypes may modulate the amount of the risk conferred by these variants. Further studies on different ethnic groups can validate these results.

A genome-wide gene-environment interaction study of breast cancer risk for women of European ancestry

BACKGROUND

Genome-wide studies of gene-environment interactions (G×E) may identify variants associated with disease risk in conjunction with lifestyle/environmental exposures. We conducted a genome-wide G×E analysis of ~ 7.6 million common variants and seven lifestyle/environmental risk factors for breast cancer risk overall and for estrogen receptor positive (ER +) breast cancer.

METHODS

Analyses were conducted using 72,285 breast cancer cases and 80,354 controls of European ancestry from the Breast Cancer Association Consortium. Gene-environment interactions were evaluated using standard unconditional logistic regression models and likelihood ratio tests for breast cancer risk overall and for ER + breast cancer. Bayesian False Discovery Probability was employed to assess the noteworthiness of each SNP-risk factor pairs.

RESULTS

Assuming a 1 × 10-5 prior probability of a true association for each SNP-risk factor pairs and a Bayesian False Discovery Probability < 15%, we identified two independent SNP-risk factor pairs: rs80018847(9p13)-LINGO2 and adult height in association with overall breast cancer risk (ORint = 0.94, 95% CI 0.92-0.96), and rs4770552(13q12)-SPATA13 and age at menarche for ER + breast cancer risk (ORint = 0.91, 95% CI 0.88-0.94).

CONCLUSIONS

Overall, the contribution of G×E interactions to the heritability of breast cancer is very small. At the population level, multiplicative G×E interactions do not make an important contribution to risk prediction in breast cancer.

Genome-wide interaction analysis of menopausal hormone therapy use and breast cancer risk among 62,370 women

Use of menopausal hormone therapy (MHT) is associated with increased risk for breast cancer. However, the relevant mechanisms and its interaction with genetic variants are not fully understood. We conducted a genome-wide interaction analysis between MHT use and genetic variants for breast cancer risk in 27,585 cases and 34,785 controls from 26 observational studies. All women were post-menopausal and of European ancestry. Multivariable logistic regression models were used to test for multiplicative interactions between genetic variants and current MHT use. We considered interaction p-values < 5 × 10-8 as genome-wide significant, and p-values < 1 × 10-5 as suggestive. Linkage disequilibrium (LD)-based clumping was performed to identify independent candidate variants. None of the 9.7 million genetic variants tested for interactions with MHT use reached genome-wide significance. Only 213 variants, representing 18 independent loci, had p-values < 1 × 105. The strongest evidence was found for rs4674019 (p-value = 2.27 × 10-7), which showed genome-wide significant interaction (p-value = 3.8 × 10-8) with current MHT use when analysis was restricted to population-based studies only. Limiting the analyses to combined estrogen-progesterone MHT use only or to estrogen receptor (ER) positive cases did not identify any genome-wide significant evidence of interactions. In this large genome-wide SNP-MHT interaction study of breast cancer, we found no strong support for common genetic variants modifying the effect of MHT on breast cancer risk. These results suggest that common genetic variation has limited impact on the observed MHT-breast cancer risk association.

A genome-wide gene-based gene-environment interaction study of breast cancer in more than 90,000 women

Background

Genome-wide association studies (GWAS) have identified more than 200 susceptibility loci for breast cancer, but these variants explain less than a fifth of the disease risk. Although gene-environment interactions have been proposed to account for some of the remaining heritability, few studies have empirically assessed this.

Methods

We obtained genotype and risk factor data from 46,060 cases and 47,929 controls of European ancestry from population-based studies within the Breast Cancer Association Consortium (BCAC). We built gene expression prediction models for 4,864 genes with a significant (P<0.01) heritable component using the transcriptome and genotype data from the Genotype-Tissue Expression (GTEx) project. We leveraged predicted gene expression information to investigate the interactions between gene-centric genetic variation and 14 established risk factors in association with breast cancer risk, using a mixed-effects score test.

Results

After adjusting for number of tests using Bonferroni correction, no interaction remained statistically significant. The strongest interaction observed was between the predicted expression of the C13orf45 gene and age at first full-term pregnancy (PGXE=4.44×10-6).

Conclusion

In this transcriptome-informed genome-wide gene-environment interaction study of breast cancer, we found no strong support for the role of gene expression in modifying the associations between established risk factors and breast cancer risk.

Impact

Our study suggests a limited role of gene-environment interactions in breast cancer risk.

Gene-Environment Interactions Relevant to Estrogen and Risk of Breast Cancer: Can Gene-Environment Interactions Be Detected Only among Candidate SNPs from Genome-Wide Association Studies?

In this study we aim to examine gene-environment interactions (GxEs) between genes involved with estrogen metabolism and environmental factors related to estrogen exposure. GxE analyses were conducted with 1970 Korean breast cancer cases and 2052 controls in the case-control study, the Seoul Breast Cancer Study (SEBCS). A total of 11,555 SNPs from the 137 candidate genes were included in the GxE analyses with eight established environmental factors. A replication test was conducted by using an independent population from the Breast Cancer Association Consortium (BCAC), with 62,485 Europeans and 9047 Asians. The GxE tests were performed by using two-step methods in GxEScan software. Two interactions were found in the SEBCS. The first interaction was shown between rs13035764 of NCOA1 and age at menarche in the GE|2df model (p-2df = 1.2 × 10-3). The age at menarche before 14 years old was associated with the high risk of breast cancer, and the risk was higher when subjects had homozygous minor allele G. The second GxE was shown between rs851998 near ESR1 and height in the GE|2df model (p-2df = 1.1 × 10-4). Height taller than 160 cm was associated with a high risk of breast cancer, and the risk increased when the minor allele was added. The findings were not replicated in the BCAC. These results would suggest specificity in Koreans for breast cancer risk.

Genetic Variation and Hot Flashes: A Systematic Review

CONTEXT

Approximately 70% of women report experiencing vasomotor symptoms (VMS, hot flashes and/or night sweats). The etiology of VMS is not clearly understood but may include genetic factors.

EVIDENCE ACQUISITION

We searched PubMed and Embase in accordance with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidance. We included studies on associations between genetic variation and VMS. We excluded studies focused on medication interventions or prevention or treatment of breast cancer.

EVIDENCE SYNTHESIS

Of 202 unique citations, 18 citations met the inclusion criteria. Study sample sizes ranged from 51 to 17 695. Eleven of the 18 studies had fewer than 500 participants; 2 studies had 1000 or more. Overall, statistically significant associations with VMS were found for variants in 14 of the 26 genes assessed in candidate gene studies. The cytochrome P450 family 1 subfamily A member 1 (CYP1B1) gene was the focus of the largest number (n = 7) of studies, but strength and statistical significance of associations of CYP1B1 variants with VMS were inconsistent. A genome-wide association study reported statistically significant associations between 14 single-nucleotide variants in the tachykinin receptor 3 gene and VMS. Heterogeneity across trials regarding VMS measurement methods and effect measures precluded quantitative meta-analysis; there were few studies of each specific genetic variant.

CONCLUSIONS

Genetic variants are associated with VMS. The associations are not limited to variations in sex-steroid metabolism genes. However, studies were few and future studies are needed to confirm and extend these findings.

E3 Ubiquitin Ligase TRIP12: Regulation, Structure, and Physiopathological Functions

The Thyroid hormone Receptor Interacting Protein 12 (TRIP12) protein belongs to the 28-member Homologous to the E6-AP C-Terminus (HECT) E3 ubiquitin ligase family. First described as an interactor of the thyroid hormone receptor, TRIP12’s biological importance was revealed by the embryonic lethality of a murine model bearing an inactivating mutation in the TRIP12 gene. Further studies showed the participation of TRIP12 in the regulation of major biological processes such as cell cycle progression, DNA damage repair, chromatin remodeling, and cell differentiation by an ubiquitination-mediated degradation of key protein substrates. Moreover, alterations of TRIP12 expression have been reported in cancers that can serve as predictive markers of therapeutic response. The TRIP12 gene is also referenced as a causative gene associated to intellectual disorders such as Clark–Baraitser syndrome and is clearly implicated in Autism Spectrum Disorder. The aim of the review is to provide an exhaustive and integrated overview of the different aspects of TRIP12 ranging from its regulation, molecular functions and physio-pathological implications.

Interactions between a Polygenic Risk Score and Non-genetic Risk Factors in Young-Onset Breast Cancer

Most gene-environmental studies have focused on breast cancers generally, the preponderance of which occur after age 50. Young-onset breast cancers (YOBC) tend to be aggressive and may be etiologically different. The goal of this analysis was to assess interactions between an established 77-SNP polygenic risk score (PRS) and non-genetic risk factors for YOBC. We constructed the PRS using a family-based study of 1,291 women diagnosed with breast cancer before age 50 and their parents and unaffected sisters. We used conditional logistic regression to analyze interactions between the PRS and 14 established risk factors. In further analyses we assessed the same interactions, but for invasive cancer, estrogen receptor (ER) positive cancer and with broader inclusion of racial/ethnic groups. Results showed a decreased association between the PRS and YOBC risk for women who had ever used hormonal birth control (odds ratio [OR] = 2.20 versus 3.89) and a stronger association between the PRS and YOBC risk in pre-menopausal women (OR = 2.46 versus 1.23). Restricting the analysis to ER+ cancers or invasive cancers or using samples from all ethnic groups produced similar results. In conclusion, the PRS may interact with hormonal birth control use and with menopausal status on risk of YOBC.

Oncogene Delta/Notch-Like EGF-Related Receptor Promotes Cell Proliferation, Invasion, and Migration in Hepatocellular Carcinoma and Predicts a Poor Prognosis

PURPOSE

To determine the expression and function of Delta/Notch-like EGF-related receptor (DNER) in hepatocellular carcinoma (HCC).

METHODS

The expression of DNER in 84 HCC tissue samples and matched adjacent noncancerous specimens, as well as HCC cells, were evaluated by quantitative real-time polymerase chain reaction (qRT-PCR) and western blotting. Survival analysis was evaluated using Kaplan-Meier method. For experiments in vitro, cell viability was measured by Cell Counting Kit-8 Assay and Colony Formation Assay. Furthermore, cell invasion and migration assays were performed with Transwell Assay.

RESULTS

The results showed that DNER was overexpressed in the tissues and cell lines of HCC (all, p < 0.05), and the upregulated expression of DNER was significantly correlated with advanced pathologic stage (p = 0.013) and pathologic-M1 (p = 0.012) in HCC patients. Survival analysis revealed that patients with high DNER levels had worse overall survival (OS) than those with low DNER levels (p = 0.004). More importantly, DNER could be an independent predictor of prognosis for OS (HR = 2.582, 95% CI 1.239-5.380, p = 0.011). In vitro, knockdown of DNER significantly suppressed cell proliferation, colony formation, cell invasion, and cell migration in HepG2 cells. Moreover, inhibition of DNER inactivated PI3K/AKT signaling pathway by downregulating the expression of p-PI3K, p-AKT, and p-70s6k.

CONCLUSIONS

Taken together, DNER could promote proliferation, migration, and invasion of HCC cells by regulating the activation of PI3K/AKT pathway, and it might act as a potential prognostic biomarker for HCC.

Association between single nucleotide polymorphisms of delta/notch-like epidermal growth factor (EGF)-related receptor (DNER) and Delta-like 1 Ligand (DLL 1) with the risk of type 2 diabetes mellitus in a Chinese Han population

The aim of the study was to investigate the association between single nucleotide polymorphisms (SNPs), rs1861612 of delta/notch-like Epidermal Growth Factor (EGF)-related receptor (DNER) and rs1884190 in the Delta-like 1 ligand (DLL1) gene and type 2 diabetes mellitus (T2DM) susceptibility in a Chinese Han population. DNER rs1861612 and DLL1 rs1884190 polymorphisms were genotyped in patients with T2DM and age- and sex-matched T2DM-free controls from a Chinese Han population. A total of 298 patients with T2DM and 500 controls were enrolled in this study. We found that TC and TT genotypes of rs1861612 and variant T were associated with a significantly increased risk of T2DM. In contrast, the AG and AA genotypes of rs1884190 were not significantly associated with the risk of T2DM, even after further stratification analysis based on age or sex. Our results showed that DNER rs1861612 C to T change and variant T genotype may contribute to T2DM in a Chinese Han population.

Interactions between breast cancer susceptibility loci and menopausal hormone therapy in relationship to breast cancer in the Breast and Prostate Cancer Cohort Consortium

Current use of menopausal hormone therapy (MHT) has important implications for postmenopausal breast cancer risk, and observed associations might be modified by known breast cancer susceptibility loci. To provide the most comprehensive assessment of interactions of prospectively collected data on MHT and 17 confirmed susceptibility loci with invasive breast cancer risk, a nested case-control design among eight cohorts within the NCI Breast and Prostate Cancer Cohort Consortium was used. Based on data from 13,304 cases and 15,622 controls, multivariable-adjusted logistic regression analyses were used to estimate odds ratios (OR) and 95 % confidence intervals (CI). Effect modification of current and past use was evaluated on the multiplicative scale. P values <1.5 × 10(-3) were considered statistically significant. The strongest evidence of effect modification was observed for current MHT by 9q31-rs865686. Compared to never users of MHT with the rs865686 GG genotype, the association between current MHT use and breast cancer risk for the TT genotype (OR 1.79, 95 % CI 1.43-2.24; P interaction = 1.2 × 10(-4)) was less than expected on the multiplicative scale. There are no biological implications of the sub-multiplicative interaction between MHT and rs865686. Menopausal hormone therapy is unlikely to have a strong interaction with the common genetic variants associated with invasive breast cancer.

Gene-environment interaction and risk of breast cancer

Hereditary, genetic factors as well as lifestyle and environmental factors, for example, parity and body mass index, predict breast cancer development. Gene–environment interaction studies may help to identify subgroups of women at high-risk of breast cancer and can be leveraged to discover new genetic risk factors. A few interesting results in studies including over 30 000 breast cancer cases and healthy controls indicate that such interactions exist. Explorative gene–environment interaction studies aiming to identify new genetic or environmental factors are scarce and still underpowered. Gene–environment interactions might be stronger for rare genetic variants, but data are lacking. Ongoing initiatives to genotype larger sample sets in combination with comprehensive epidemiologic databases will provide further opportunities to study gene–environment interactions in breast cancer. However, based on the available evidence, we conclude that associations between the common genetic variants known today and breast cancer risk are only weakly modified by environmental factors, if at all.

Invasive lobular breast cancer: the prognostic impact of histopathological grade, E-cadherin and molecular subtypes

Aims

The aim of this study was to compare breast cancer specific survival (BCSS) for invasive lobular carcinoma (ILC) and invasive ductal carcinoma (IDC) and, further, to evaluate critically the prognostic value of histopathological grading of ILC and examine E-cadherin as a prognostic marker in ILC.

Methods and results

The study comprised 116 lobular and 611 ductal breast carcinomas occurring between 1961 and 2008. All cases had been classified previously according to histopathological type and grade, stained for oestrogen receptor (ER), progesterone receptor (PR), antigen Ki67 (Ki67), epithelial growth factor receptor (EGFR), cytokeratin 5 (CK5) and human epidermal growth factor receptor 2 (HER2) and classified into molecular subtypes. For the present study, immunohistochemical staining for E-cadherin was performed. The Kaplan–Meier method and Cox proportional hazards models were used in the analyses. Grade 2 tumours comprised 85.3% of the lobular tumours and 51.9% of the ductal tumours. BCSS in ILC grade 2 was comparable to that of IDC grade 3. E-cadherin-negative ILC had a poorer prognosis compared to E-cadherin positive ILC and to IDC regardless of E-cadherin status.

Conclusions

The implication of histopathological grading may differ in ILC compared to IDC. E-cadherin may be useful in prognostication in ILC and thereby influence the determination of treatment strategies for this group of women.

Identification of mutant genes with high-frequency, high-risk, and high-expression in lung adenocarcinoma

BACKGROUND

To identify mutant genes with high-frequency-risk-expression between lung adenocarcinoma and normal samples.

METHODS

The ribonucleic acid RNA-Seq data GSE34914 and GSE37765 were downloaded from the Gene Expression Omnibus database, including 12 lung adenocarcinoma samples and six controls. All RNA-Seq reads were processed and the gene-expression level was calculated. Single nucleotide variation (SNV) was analyzed and the locations of mutant sites were recorded. In addition, the frequency and risk-level of mutant genes were calculated. Gene Ontology (GO) functional analysis was performed. The reported cancer genes were searched in tumor suppressor genes, Cancer Genes, and the Catalogue of Somatic Mutations in Cancer (COSMIC) database.

RESULTS

The SNV annotations of somatic mutation sites showed that 70% of mutation sites in the exon region occurred in the coding sequence (CDS). Thyroid hormone receptor interactor (TRIP)12 was identified with the highest frequency. A total of 118 mutant genes with high frequency and high-risk were selected and significantly enriched into several GO terms. No base mutation of cyclin C (CCNC) or RAB11A was recorded. At fragments per kilobase per million reads (FPKM) ≥ 56.5, reported tumor suppressor genes catenin (cadherin-associated protein), delta (CTNND)1, dual specificity phosphatase (DUSP)6, malate dehydrogenase (MDH)1 and RNA binding motif protein (RBM)5, were identified. Notably, signal transducer and activator of transcription 2 (STAT2) was the only transcription factor (TF) with high-risk mutation and its expression was detected.

CONCLUSION

For the mutant genes with high-frequency-risk-expression, CTNND1, DUSP6, MDH1 and RBM5 were identified. TRIP12 might be a potential cancer-related gene, and expression of TF STAT2 with high-risk was detected. These mutant gene candidates might promote the development of lung adenocarcinoma and provide new diagnostic potential targets for treatment.

Genetic modifiers of menopausal hormone replacement therapy and breast cancer risk: a genome-wide interaction study

Women using menopausal hormone therapy (MHT) are at increased risk of developing breast cancer (BC). To detect genetic modifiers of the association between current use of MHT and BC risk, we conducted a meta-analysis of four genome-wide case-only studies followed by replication in 11 case-control studies. We used a case-only design to assess interactions between single-nucleotide polymorphisms (SNPs) and current MHT use on risk of overall and lobular BC. The discovery stage included 2920 cases (541 lobular) from four genome-wide association studies. The top 1391 SNPs showing P values for interaction (Pint) <3.0 × 10(-3) were selected for replication using pooled case-control data from 11 studies of the Breast Cancer Association Consortium, including 7689 cases (676 lobular) and 9266 controls. Fixed-effects meta-analysis was used to derive combined Pint. No SNP reached genome-wide significance in either the discovery or combined stage. We observed effect modification of current MHT use on overall BC risk by two SNPs on chr13 near POMP (combined Pint≤8.9 × 10(-6)), two SNPs in SLC25A21 (combined Pint≤4.8 × 10(-5)), and three SNPs in PLCG2 (combined Pint≤4.5 × 10(-5)). The association between lobular BC risk was potentially modified by one SNP in TMEFF2 (combined Pint≤2.7 × 10(-5)), one SNP in CD80 (combined Pint≤8.2 × 10(-6)), three SNPs on chr17 near TMEM132E (combined Pint≤2.2×10(-6)), and two SNPs on chr18 near SLC25A52 (combined Pint≤4.6 × 10(-5)). In conclusion, polymorphisms in genes related to solute transportation in mitochondria, transmembrane signaling, and immune cell activation are potentially modifying BC risk associated with current use of MHT. These findings warrant replication in independent studies.