Potential role of fructose on human colon DNA methylation in racial disparities observed for colorectal cancer risk

Background and aims

An increasing body of observational studies has linked fructose intake to colorectal cancer (CRC). African Americans (AAs) are significantly more likely than European Americans to consume greater quantities of fructose and to develop right-side colon cancer. Yet, a mechanistic link between these two associations remains poorly defined. We aimed to identify differentially methylated regions (DMRs) associated with dietary fructose consumption measures obtained from food frequency questionnaires in a cohort of normal colon biopsies derived from AA men and women (n=79).

Methods

DNA methylation data from this study was obtained using the Illumina Infinium MethylationEPIC kit and is housed under accession GSE151732. DMR analysis was carried out using DMRcate in right and matched left colon, separately. Secondary analysis of CRC tumors was carried out using data derived from TCGA-COAD, GSE101764 and GSE193535. Differential expression analysis was carried out on CRC tumors from TCGA-COAD using DESeq2 .

Results

We identified 4,263 right-side fructose-DMRs. In contrast, only 24 DMRs survived multiple testing corrections (FDR<0.05) in matched, left colon. To identify targets by which dietary fructose drives CRC risk, we overlaid these findings with data from three CRC tumor datasets. Remarkably, almost 50% of right-side fructose-DMRs overlapped regions associated with CRC in at least one of three datasets. TNXB and CDX2 ranked among the most significant fructose risk DMRs in right and left colon respectively that also displayed altered gene expression in CRC tumors.

Conclusions

Our mechanistic data support the notion that fructose has a greater CRC-related effect in right than left AA colon, alluding to a potential role for fructose in contributing to racial disparities in CRC.

A Race-Specific, DNA Methylation Analysis of Aging in Normal Rectum: Implications for the Biology of Aging and Its Relationship to Rectal Cancer

Approximately 90% of colorectal cancer (CRC) develop over the age of 50, highlighting the important role of aging in CRC risk. African Americans (AAs) shoulder a greater CRC burden than European Americans (EA) and are more likely to develop CRC at a younger age. The effects of aging in AA and EA normal rectal tissue have yet to be defined. Here, we performed epigenome-wide DNA methylation analysis in the first, large-scale biracial cohort of normal rectum (n = 140 samples). We identified increased epigenetic age acceleration in EA than AA rectum (p = 3.91 × 10^-4) using linear regression. We also identified differentially methylated regions (DMRs) associated with chronological aging in AA and EA, separately using DMRcate. Next, a consensus set of regions associated with cancer was identified through DMR analysis of two rectal cancer cohorts. The vast majority of AA DMRs were present in our analysis of aging in rectum of EA subjects, though rates of epigenetic drift were significantly greater in AA (p = 1.94 × 10^-45). However, 3.66-fold more DMRs were associated with aging in rectum of EA subjects, many of which were also associated with rectal cancer. Our findings reveal a novel relationship between race, age, DNA methylation and rectal cancer risk that warrants further investigation.

DNA methylation analysis of normal colon organoids from familial adenomatous polyposis patients reveals novel insight into colon cancer development

BACKGROUND

Familial adenomatous polyposis (FAP) is an inherited colorectal cancer (CRC) syndrome resulting from germ line mutations in the adenomatous polyposis coli (APC) gene. While FAP accounts for less than 1% of all CRC cases, loss of APC expression is seen in > 80% of non-hereditary CRCs. To better understand molecular mechanisms underlying APC-driven CRC, we performed an epigenome-wide analysis of colon organoids derived from normal-appearing colons of FAP patients versus healthy subjects to identify differentially methylated regions (DMRs) that may precede the onset of CRC.

RESULTS

We identified 358 DMRs when comparing colon organoids of FAP patients to those of healthy subjects (FDR < 0.05, |mean beta difference| = 5%). Of these, nearly 50% of DMRs were also differentially methylated in at least one of three CRC tumor and normal adjacent tissue (NAT) cohorts (TCGA-COAD, GSE193535 and ColoCare). Moreover, 27 of the DMRs mapped to CRC genome-wide association study (GWAS) loci. We provide evidence suggesting that some of these DMRs led to significant differences in gene expression of adjacent genes using quantitative PCR. For example, we identified significantly greater expression of five genes: Kazal-type serine peptidase inhibitor domain 1 (KAZALD1, P = 0.032), F-Box and leucine-rich repeat protein 8 (FBXL8, P = 0.036), TRIM31 antisense RNA 1 (TRIM31-AS1, P = 0.036), Fas apoptotic inhibitory molecule 2 (FAIM2, P = 0.049) and (Collagen beta (1-0)galactosyltransferase 2 (COLGALT2, P = 0.049). Importantly, both FBXL8 and TRIM31-AS1 were also significantly differentially expressed in TCGA-COAD tumor versus matched NAT, supporting a role for these genes in CRC tumor development.

CONCLUSIONS

We performed the first DNA methylome-wide analysis of normal colon organoids derived from FAP patients compared to those of healthy subjects. Our results reveal that normal colon organoids from FAP patients exhibit extensive epigenetic differences compared to those of healthy subjects that appear similar to those exhibited in CRC tumor. Our analyses therefore identify DMRs and candidate target genes that are potentially important in CRC tumor development in FAP, with potential implications for non-hereditary CRC.

Beyond GWAS of Colorectal Cancer: Evidence of Interaction with Alcohol Consumption and Putative Causal Variant for the 10q24.2 Region

BACKGROUND

Currently known associations between common genetic variants and colorectal cancer explain less than half of its heritability of 25%. As alcohol consumption has a J-shape association with colorectal cancer risk, nondrinking and heavy drinking are both risk factors for colorectal cancer.

METHODS

Individual-level data was pooled from the Colon Cancer Family Registry, Colorectal Transdisciplinary Study, and Genetics and Epidemiology of Colorectal Cancer Consortium to compare nondrinkers (≤1 g/day) and heavy drinkers (>28 g/day) with light-to-moderate drinkers (1-28 g/day) in GxE analyses. To improve power, we implemented joint 2df and 3df tests and a novel two-step method that modifies the weighted hypothesis testing framework. We prioritized putative causal variants by predicting allelic effects using support vector machine models.

RESULTS

For nondrinking as compared with light-to-moderate drinking, the hybrid two-step approach identified 13 significant SNPs with pairwise r2 > 0.9 in the 10q24.2/COX15 region. When stratified by alcohol intake, the A allele of lead SNP rs2300985 has a dose-response increase in risk of colorectal cancer as compared with the G allele in light-to-moderate drinkers [OR for GA genotype = 1.11; 95% confidence interval (CI), 1.06-1.17; OR for AA genotype = 1.22; 95% CI, 1.14-1.31], but not in nondrinkers or heavy drinkers. Among the correlated candidate SNPs in the 10q24.2/COX15 region, rs1318920 was predicted to disrupt an HNF4 transcription factor binding motif.

CONCLUSIONS

Our study suggests that the association with colorectal cancer in 10q24.2/COX15 observed in genome-wide association study is strongest in nondrinkers. We also identified rs1318920 as the putative causal regulatory variant for the region.

IMPACT

The study identifies multifaceted evidence of a possible functional effect for rs1318920.

SNPs associated with colorectal cancer at 15q13.3 affect risk enhancers that modulate GREM1 gene expression

Several genome wide association studies of colorectal cancer (CRC) have identified single nucleotide polymorphisms (SNPs) on chromosome 15q13.3 associated with CRC risk. To identify functional variant(s) underlying this association, we investigated SNPs in linkage disequilibrium with the risk‐associated SNP rs4779584 that overlapped regulatory regions/enhancer elements characterized in colon‐related tissues and cells. We identified several SNP‐containing regulatory regions that exhibited enhancer activity in vitro, including one SNP (rs1406389) that correlated with allele‐specific effects on enhancer activity. Deletion of either this enhancer or another enhancer that had previously been reported in this region correlated with decreased expression of GREM1 following CRISPR/Cas9 genome editing. That GREM1 is one target of these enhancers was further supported by an expression quantitative trait loci correlation between rs1406389 and GREM1 expression in the transverse but not sigmoid colon in the Genotype‐Tissue Expression dataset. Taken together, we conclude that the 15q13.3 region contains at least two functional variants that map to distinct enhancers and impact CRC risk through modulation of GREM1 expression.