Genetic variation in sex-steroid receptors and synthesizing enzymes and colorectal cancer risk in women

Prognostic value of fatty acid metabolism-related genes in colorectal cancer and their potential implications for immunotherapy

Introduction

Colorectal cancer is one of the most common gastrointestinal cancers and the second leading cause of cancer-related death. Although colonoscopy screening has greatly improved the early diagnosis of colorectal cancer, its recurrence and metastasis are still significant problems. Tumour cells usually have the hallmark of metabolic reprogramming, while fatty acids play important roles in energy storage, cell membrane synthesis, and signal transduction. Many pathways of fatty acid metabolism (FAM) are involved in the occurrence and development of colon cancer, and the complex molecular interaction network contains a variety of genes encoding key enzymes and related products.

Methods

Clinical information and RNA sequencing data were collected from TCGA and GEO databases. The prognosis model of colon cancer was constructed by LASSO-Cox regression analysis among the selected fatty acid metabolism genes with differential expression. Nomogram for the prognosis model was also constructed in order to analyze its value in evaluating the survival and clinical stage of the colon cancer patients. The differential expression of the selected genes was verified by qPCR and immunohistochemistry. GSEA and GSVA were used to analyze the enrichment pathways for high- and low-risk groups. CIBERSORT was used to analyze the immune microenvironment of colon cancer and to compare the infiltration of immune cells in the high- and low-risk groups. The "circlize" package was used to explore the correlation between the risk score signature and immunotherapy for colon cancer.

Results

We analysed the differential expression of 704 FAM-related genes between colon tumour and normal tissue and screened 10 genes with prognostic value. Subsequently, we constructed a prognostic model for colon cancer based on eight optimal FAM genes through LASSO Cox regression analysis in the TCGA-COAD dataset, and its practicality was validated in the GSE39582 dataset. Moreover, the risk score calculated based on the prognostic model was validated as an independent prognostic factor for colon cancer patients. We further constructed a nomogram composed of the risk score signature, age and American Joint Committee on Cancer (AJCC) stage for clinical application. The colon cancer cohort was divided into high- and low-risk groups according to the optimal cut-off value, and different enrichment pathways and immune microenvironments were depicted in the groups.

Discussion

Since the risk score signature was significantly correlated with the expression of immune checkpoint molecules, the prognostic model might be able to predict the immunotherapy response of colon cancer patients. In summary, our findings expand the prognostic value of FAM-related genes in colon cancer and provide evidence for their application in guiding immunotherapy.

Common variability in oestrogen-related genes and pancreatic ductal adenocarcinoma risk in women

The incidence of pancreatic ductal adenocarcinoma (PDAC) is different among males and females. This disparity cannot be fully explained by the difference in terms of exposure to known risk factors; therefore, the lower incidence in women could be attributed to sex-specific hormones. A two-phase association study was conducted in 12,387 female subjects (5436 PDAC cases and 6951 controls) to assess the effect on risk of developing PDAC of single nucleotide polymorphisms (SNPs) in 208 genes involved in oestrogen and pregnenolone biosynthesis and oestrogen-mediated signalling. In the discovery phase 14 polymorphisms showed a statistically significant association (P < 0.05). In the replication none of the findings were validated. In addition, a gene-based analysis was performed on the 208 selected genes. Four genes (NR5A2, MED1, NCOA2 and RUNX1) were associated with PDAC risk, but only NR5A2 showed an association (P = 4.08 × 10-5) below the Bonferroni-corrected threshold of statistical significance. In conclusion, despite differences in incidence between males and females, our study did not identify an effect of common polymorphisms in the oestrogen and pregnenolone pathways in relation to PDAC susceptibility. However, we validated the previously reported association between NR5A2 gene variants and PDAC risk.

Genome-Wide Interaction Analysis of Genetic Variants With Menopausal Hormone Therapy for Colorectal Cancer Risk

BACKGROUND

The use of menopausal hormone therapy (MHT) may interact with genetic variants to influence colorectal cancer (CRC) risk.

METHODS

We conducted a genome-wide, gene-environment interaction between single nucleotide polymorphisms and the use of any MHT, estrogen only, and combined estrogen-progestogen therapy with CRC risk, among 28 486 postmenopausal women (11 519 CRC patients and 16 967 participants without CRC) from 38 studies, using logistic regression, 2-step method, and 2- or 3-degree-of-freedom joint test. A set-based score test was applied for rare genetic variants.

RESULTS

The use of any MHT, estrogen only and estrogen-progestogen were associated with a reduced CRC risk (odds ratio [OR] = 0.71, 95% confidence interval [CI] = 0.64 to 0.78; OR = 0.65, 95% CI = 0.53 to 0.79; and OR = 0.73, 95% CI = 0.59 to 0.90, respectively). The 2-step method identified a statistically significant interaction between a GRIN2B variant rs117868593 and MHT use, whereby MHT-associated CRC risk was statistically significantly reduced in women with the GG genotype (OR = 0.68, 95% CI = 0.64 to 0.72) but not within strata of GC or CC genotypes. A statistically significant interaction between a DCBLD1 intronic variant at 6q22.1 (rs10782186) and MHT use was identified by the 2-degree-of-freedom joint test. The MHT-associated CRC risk was reduced with increasing number of rs10782186-C alleles, showing odds ratios of 0.78 (95% CI = 0.70 to 0.87) for TT, 0.68 (95% CI = 0.63 to 0.73) for TC, and 0.66 (95% CI = 0.60 to 0.74) for CC genotypes. In addition, 5 genes in rare variant analysis showed suggestive interactions with MHT (2-sided P < 1.2 × 10-4).

CONCLUSION

Genetic variants that modify the association between MHT and CRC risk were identified, offering new insights into pathways of CRC carcinogenesis and potential mechanisms involved.

Evaluation of the epidemiological and prognosis significance of ESR2 rs3020450 polymorphism in ovarian cancer

AIM

To investigate the correlation between the polymorphism of estrogen receptor β gene (ESR2) rs3020450 and cancer susceptibility, and explore the epidemiological significance and the effect of ESR2 expression levels on the prognosis of ovarian cancer.

METHODS

Based on meta-analysis the association between ESR2 rs3020450 polymorphism and cancer susceptibility was estimated and a case-control design was used to verify this result in ovarian cancer. The epidemiological effect of ESR2 rs3020450 polymorphism was assessed by attributable risk percentage (ARP) and population attributable risk percentage (PARP). Kaplan Meier plotters were used to evaluate overall survival (OS) and progression-free survival (PFS) in ovarian cancer patients and GEPIA for the differential expression of ESR2 levels in ovarian cancer and adjacent normal tissues.

RESULTS

The pooled analysis indicated no significant correlation between the ESR2 rs3020450 polymorphism and the cancer susceptibility. In the stratified analysis by cancer types, significantly decreased risk was found in ovarian cancer (AG vs GG: OR = 0.73, 95%CI: 0.53-0.97, P = 0.03). Unconditional logistic regression results of case-control study in ovarian cancer observed significant differences in all comparisons (AG vs GG: OR = 0.81, 95%CI: 0.62-0.98, P = 0.04; AA vs GG: OR = 0.63, 95%CI: 0.42-0.92, P = 0.01 and AG + AA vs GG: OR = 0.73, 95%CI: 0.53-0.96, P < 0.001). Based on meta-analysis and case-control pooled results, ARP and PARP were evaluated respectively in allele (21.95% and7.97%), heterozygote (36.99% and 12.11%) and dominant model (36.84% and 12.97%) of rs3020450 polymorphism in ovarian cancer. The expression levels of ESR2 in normal tissues was significantly higher than that in cancer tissues (OV, Median, 4.7:0.21), and significant correlations were observed between high ESR2 expression levels and long OS (HR = 0.80, 95%CI: 0.70-0.92, P = 0.002) and PFS (HR = 0.767, 95%Cl: 0.67-0.88, P < 0.001).

CONCLUSION

Our results indicated that ESR2 rs3020450 polymorphism was associated with ovarian cancer risk from epidemiological perspective, and high ESR2 expression levels was associated with long survival in patients with ovarian cancer.

Association between estrogen receptor 1 (ESR1) and leukemia inhibitory factor (LIF) polymorphisms can help in the prediction of recurrent implantation failure

OBJECTIVE

To investigate an association between polymorphisms related to the implantation process that together could help in the prediction of recurrent implantation failure (RIF).

DESIGN

Cohort study.

SETTING

Private fertility center and reproductive genetics laboratory.

PATIENT(S)

Forty-four women presenting RIF, who were included in study group (RIF group), and two control groups, one with 63 women who were attended at our service and became pregnant after the first IVF/intracytoplasmic sperm injection attempt (control group I) and other with 65 fertile women who had at least two children without any treatment and no history of miscarriage (control group II).

INTERVENTION(S)

None.

MAIN OUTCOME MEASURE(S)

Genotyping was performed in the intron region of TP63, VEGFA, MMP2, ESR1, and ESR2 genes and in the 3' untranslated region of the LIF gene on genomic DNA using real-time polymerase chain reaction.

RESULT(S)

The presence of ESR1/AA (rs12199722) and LIF/GT (rs929271) genotypes was more frequent in the RIF group, leading to a 7.9-fold increase in the chance of women presenting with RIF when compared with women who became pregnant on their first cycle of IVF/intracytoplasmic sperm injection and a 2.8-fold increase when compared with women who became pregnant without treatment.

CONCLUSION(S)

The association between ESR1 and LIF polymorphisms can help in the prediction of RIF.

Association of PvuII and XbaI polymorphisms on estrogen receptor alpha (ESR1) gene to changes into serum lipid profile of post-menopausal women: Effects of aging, body mass index and breast cancer incidence

Estrogen is a steroidal hormone involved in several physiological functions in the female body including regulation of serum lipid metabolism and breast cancer (BC). Estrogen actions on serum lipids mostly occur through its binding to intracellular Estrogen Receptor alpha (ERalpha) isoform, expressed in most of tissues. This gene (ESR1) exhibit many polymorphic sites (SNPs) located either on translated and non-translated regions that regulate ERalpha protein expression and function. This paper aimed to investigate the association of two intronic SNPs of ESR1 gene, namely c454-397T>C (PvuII) and c454-351A>G (XbaI) to alterations in serum levels of total cholesterol (T-chol), total lipid (TL), low density lipoprotein cholesterol (LDL), high density lipoprotein (HDL), and triglycerides (TG) in a cohort of post-menopausal women. In addition, we aimed to associate presence of these SNPs to development of BC along 5 years period. To do so, a group of healthy 499, highly miscigenated, post-menopausal Brazilian women, were carried using PCR-FRLP technique and further confirmed by automatic sequence analysis as well followed through 5 years for BC incidence. Measurements of serum lipid profile by standard commercial methods were carried individually whereas Dual Energy X-ray Absorciometry (DXA) measured Body Mass Indexes (BMI), Fat Mass (FM), Lean Body Mass (LBM), and Body Water Content (BWC). No effects of PvuII SNP on ESR1 gene were observed on patient´s serum T-chol, TL, LDL, HDL, and TG. However, c454-397T>C PvuII SNP is associated to lower body fat and higher levels of lean mass and body water and lower incidence of BC. On the other hand, statistically significant effect of XbaI c454-351A>G SNP on serum TG and TL levels. Patients homozygous for X allele were followed up from 2010–2015. They showed higher incidence of breast cancer (BC) when compared to either heterozygous and any P allele combination. Moreover, the increasing of TG and TL serum concentrations associated to SNP XbaI c454-351A>G on ESR1 gene is enhanced in both obese (higher BMI) and elder women. Taken together, these results suggested that XbaI c454-351A>G SNP is associated to changes in lipid profile, particularly in serum TG and TL, in this cohort of post-menopausal woman. Also, this paper shows another link between obesity and BC corroborating the current thesis that both diseases are interlinked.

CYP24A1 variant modifies the association between use of oestrogen plus progestogen therapy and colorectal cancer risk

BACKGROUND

Menopausal hormone therapy (MHT) use has been consistently associated with a decreased risk of colorectal cancer (CRC) in women. Our aim was to use a genome-wide gene-environment interaction analysis to identify genetic modifiers of CRC risk associated with use of MHT.

METHODS

We included 10 835 postmenopausal women (5419 cases and 5416 controls) from 10 studies. We evaluated use of any MHT, oestrogen-only (E-only) and combined oestrogen-progestogen (E+P) hormone preparations. To test for multiplicative interactions, we applied the empirical Bayes (EB) test as well as the Wald test in conventional case-control logistic regression as primary tests. The Cocktail test was used as secondary test.

RESULTS

The EB test identified a significant interaction between rs964293 at 20q13.2/CYP24A1 and E+P (interaction OR (95% CIs)=0.61 (0.52-0.72), P=4.8 × 10(-9)). The secondary analysis also identified this interaction (Cocktail test OR=0.64 (0.52-0.78), P=1.2 × 10(-5) (alpha threshold=3.1 × 10(-4)). The ORs for association between E+P and CRC risk by rs964293 genotype were as follows: C/C, 0.96 (0.61-1.50); A/C, 0.61 (0.39-0.95) and A/A, 0.40 (0.22-0.73), respectively.

CONCLUSIONS

Our results indicate that rs964293 modifies the association between E+P and CRC risk. The variant is located near CYP24A1, which encodes an enzyme involved in vitamin D metabolism. This novel finding offers additional insight into downstream pathways of CRC etiopathogenesis.

Estrogen receptor 1 (ESR1) genetic variations in cancer risk: a systematic review and meta-analysis

BACKGROUND

Emerging published data on the association between single nucleotide polymorphisms (SNPs) in the estrogen receptor 1 (ESR1) gene and cancer susceptibility are inconsistent. This review and meta-analysis is performed to derive a more precise evaluation of this relationship.

METHODS

The literature search of PubMed, Embase, Web of Science and CNKI databases was conducted from their inception through June 2014. Crude odds ratios (ORs) with 95% confidence intervals (CIs) were calculated to assess the association.

RESULTS

Twenty-two literatures were enrolled in this meta-analysis. The results indicated that ESR1 rs1801132 (C>G) was associated with cancer risk in Caucasian populations. However, the results of stratified analysis by cancer type and source of controls indicated that no significant association was found. Furthermore, rs2077647 (A>G) was only associated with an increased risk of hepatocellular carcinoma, but was an adverse effect on cancer risk in Caucasian populations.

CONCLUSIONS

This present meta-analysis indicated that rs1801132 (C>G) and rs2077647 (A>G) may be protective factors in Caucasian populations. Meanwhile, rs2077647 (A>G) may be closely related with hepatocellular carcinoma.

No association between germline variation in catechol-O-methyltransferase and colorectal cancer survival in postmenopausal women

OBJECTIVE

Sex steroid hormones play a role in colorectal cancer (CRC) development, but little is known about their influence on tumor progression and metastasis. Because catechol-O-methyltransferase (COMT; 22q11.21) activity is an important component of estrogen-mediated carcinogenesis, we hypothesized that germline variation in COMT may be associated with CRC survival.

METHODS

We identified 10 single nucleotide polymorphisms that tagged variation across two isoforms of COMT in 2,458 women with CRC from the Nurses' Health Study, Postmenopausal Hormones Supplementary Study to the Colon Cancer Family Registry, VITamins And Lifestyle Study, and Women's Health Initiative. All four studies participated in the Genetics and Epidemiology of Colorectal Cancer Consortium.

RESULTS

During a median follow-up of 7 years across all studies, there were 799 deaths, including 566 deaths from CRC. Based on multiple comparisons, no associations between single nucleotide polymorphisms and CRC-specific or overall survival reached statistical significance, including the well-characterized Val108/158Met polymorphism (rs4680; CRC-specific survival: hazard ratio per minor allele, 1.04; 95% CI, 0.92-1.17; overall survival: hazard ratio per minor allele, 1.01; 95% CI, 0.90-1.14).

CONCLUSIONS

In this large study of women with CRC, we find no evidence that common inherited variation in COMT is associated with survival time after diagnosis.

Oestrogen and colorectal cancer: mechanisms and controversies

AIM

The role of oestrogen metabolism and action in colorectal cancer (CRC) is controversial. An extensive review of the current literature, encompassing epidemiological evidence, systemic and peripheral oestrogen concentrations, 17β-hydroxysteroid dehydrogenase (17β-HSD) and aromatase in CRC, steroid sulphatase (STS)/oestrone sulphotransferase (EST) and in vitro and in vivo genomic effects was therefore undertaken.

METHODS

A literature search (key words: colorectal cancer, oestrogen, oestrogen receptor, 17β-HSD, STS, organic anion transporter) was performed using Embase, Medline, and Pubmed and papers were evaluated on scientific relevance on an individual basis.

RESULTS

Epidemiological data highlights that premenopausal women, or postmenopausal women taking hormone replacement therapy, are significantly less likely than males to develop CRC. This implies that oestrogen signalling is most likely involved in CRC physiology and aetiology. Little is known about oestrogen metabolism in the colon. However, the expression of 17β-HSD, STS, and EST, enzymes involved in oestrogen metabolism, have shown prognostic significance. Evidence also suggests that protective effects are modulated through oestrogen receptor beta, although which metabolite of oestrogen, oestradiol (E2) or oestrone (E1), is more active remains undefined. To complicate matters, the changes in the peripheral ratios of these enzymes, oestrogens and receptors most likely influences CRC progression.

CONCLUSION

Epidemiological evidence, now supported by in vitro and in vivo studies, strongly associates oestrogen action and metabolism with CRC. Initially protective against CRC, once developed, results suggests that oestrogens increase proliferation. Consequently, hormone-ablation therapy, already successful against breast and prostate cancer, may be effective against CRC.

Common single-nucleotide polymorphisms in the estrogen receptor β promoter are associated with colorectal cancer survival in postmenopausal women

Loss of estrogen receptor β (ERβ) expression in the gut is associated with colorectal cancer (CRC) initiation and progression. Germline single-nucleotide polymorphisms (SNP) in genes for the sex-steroid hormone receptors are not strongly associated with CRC risk; however, these SNPs have not previously been evaluated in relation to survival after diagnosis. We enrolled 729 women, ages 50 to 74, diagnosed with invasive CRC between 1997 and 2002 in 13 counties covered by the Seattle-Puget Sound Surveillance Epidemiology and End Results cancer registry. Participants provided germline DNA. We selected 99 tag-SNPs for the androgen receptor (AR), ERα (ESR1), ERβ (ESR2), and progesterone receptor (PGR) genes. Mortality outcomes were ascertained from the National Death Index. During a median of 6.6 years of follow-up, 244 deaths occurred (161 from CRC). We identified 20 SNPs (12 of ESR2 and 8 of PGR) for replication in 1,729 women diagnosed with incident invasive CRC (555 deaths; 405 from CRC) from three prospective cohort studies that participate in the Genetics and Epidemiology of Colorectal Cancer Consortium. Three correlated SNPs in the promoter of ESR2 (rs2987983, rs3020443, and rs2978381) were statistically significant predictors of CRC-specific and overall survival. Minor alleles of each were associated with improved survival [for rs2987983, CRC-specific HR, 0.77; 95% confidence interval (CI), 0.60-0.99 in the initial study, and HR, 0.79; CI, 0.64-0.98 in replication]. No associations were noted for SNPs of AR, ESR1, or PGR. SNPs in the promoter of ESR2 may be important to pathways related to the association between ERβ and tumor progression and metastasis.

Genetic variation in insulin pathway genes and distal colorectal adenoma risk

BACKGROUND

Insulin, glucose, and other insulin-related proteins that mediate insulin signaling are associated with colorectal neoplasia risk, but associations with common genetic variation in insulin axis genes are less clear. In this study, we used a comprehensive tag single-nucleotide polymorphisms (SNPs) approach to define genetic variation in six insulin axis genes (IGF1, IGF2, IGFBP1, IGFBP3, IRS1, and IRS2) and three genes associated with estrogen signaling (ESR1, ESR2, and PGR).

METHODS

We assessed associations between SNPs and distal colorectal adenoma (CRA) risk in a case-control study of 1,351 subjects. Cases were individuals with one or more adenomas diagnosed during sigmoidoscopy, and controls were individuals with no adenomas at the sigmoidoscopy exam. We used unconditional logistic regression assuming an additive model to assess SNP-specific risks adjusting for multiple comparisons with P (act).

RESULTS

Distal adenoma risk was significantly increased for one SNP in IGF2 [per minor allele OR = 1.41; 95 % confidence interval (CI) = 1.16, 1.67; P (act) = 0.005] and decreased for an ESR2 SNP (per minor allele OR = 0.78; 95 % CI = 0.66, 0.91; P (act) = 0.041). There was no statistically significant heterogeneity of these associations by race, sex, BMI, physical activity, or, in women, hormone replacement therapy use. Risk estimates did not differ in the colon versus rectum or for smaller (<1 cm) versus larger (>1 cm) adenomas.

CONCLUSIONS

These data suggest that selected genetic variability in IGF2 and ESR2 may be modifiers of CRA risk.

Estrogen and progesterone-related gene variants and colorectal cancer risk in women

Background

Observational studies and randomized trials have suggested that estrogens and/or progesterone may lower the risk for colorectal cancer. Inherited variation in the sex-hormone genes may be one mechanism by which sex hormones affect colorectal cancer, although data are limited.

Method

We conducted a comprehensive evaluation of single nucleotide polymorphisms (SNPs) in genes encoding 3 hormone receptors (ESR1, ESR2, PGR) and 5 hormone synthesizers (CYP19A1 and CYP17A1, HSD17B1, HSD17B2, HSD17B4) among 427 women with incident colorectal cancer and 871 matched controls who were Caucasians of European ancestry from 93676 postmenopausal women enrolled in the Women's Health Initiative Observational cohort. A total of 242 haplotype-tagging and functional SNPs in the 8 genes were included for analysis. Unconditional logistic regression with adjustment for age and hysterectomy status was used to estimate odds ratios (ORs) and 95% confidence intervals (CIs).

Results

We observed a weak association between the CYP17A1 rs17724534 SNP and colorectal cancer risk (OR per risk allele (A) = 1.39, 95% CI = 1.09-1.78, corrected p-value = 0.07). In addition, a suggestive interaction between rs17724534 and rs10883782 in 2 discrete LD blocks of CYP17A1 was observed in relation to colorectal cancer (empirical p value = 0.04). Moreover, one haplotype block of CYP19A1 was associated with colorectal cancer (corrected global p value = 0.02), which likely reflected the association with the tagging SNP, rs1902584, in the block.

Conclusion

Our findings offer some support for a suggestive association of CYP17A1 and CYP19A1 variants with colorectal cancer risk.

Variation in the CYP19A1 gene and risk of colon and rectal cancer

CYP19A1, or aromatase, influences estrogen-metabolizing enzymes and may influence cancer risk. We examine variation in the CYP19A1 gene and risk of colorectal cancer using data from population-based case-control studies (colon n = 1,574 cases, 1,970 controls; rectal n = 791 cases, 999 controls). Four SNPs were statistically significantly associated with colon cancer and four were associated with rectal cancer. After adjustment for multiple comparisons, the AA genotype of rs12591359 was associated with an increased risk of colon cancer (OR 1.44 95% CI 1.16-1.80) and the AA genotype of rs2470144 was associated with a reduced risk of rectal cancer (OR 0.65 95% CI 0.50-0.84). Variants of CYP19A1 were associated with CIMP+ and CIMP+/KRAS2-mutated tumors. CT/TT genotypes of rs1961177 were significantly associated with an increased likelihood of a MSI+ colon tumor (OR 1.77 95% CI 1.26-2.37). We observed statistically significant interactions between genetic variation in NFκB1 and CYP19A1 for both colon and rectal cancer. Our data suggest the importance of CYP19A1 in the development of colon and rectal cancer and that estrogen may influence risk through an inflammation-related mechanism.