Genome-wide interaction study of dietary intake of fibre, fruits, and vegetables with risk of colorectal cancer.

BACKGROUND: Consumption of fibre, fruits and vegetables have been linked with lower colorectal cancer (CRC) risk. A genome-wide gene-environment (G × E) analysis was performed to test whether genetic variants modify these associations. METHODS: A pooled sample of 45 studies including up to 69,734 participants (cases: 29,896; controls: 39,838) of European ancestry were included. To identify G × E interactions, we used the traditional 1--degree-of-freedom (DF) G × E test and to improve power a 2-step procedure and a 3DF joint test that investigates the association between a genetic variant and dietary exposure, CRC risk and G × E interaction simultaneously. FINDINGS: The 3-DF joint test revealed two significant loci with p-value <5 × 10-8. Rs4730274 close to the SLC26A3 gene showed an association with fibre (p-value: 2.4 × 10-3) and G × fibre interaction with CRC (OR per quartile of fibre increase = 0.87, 0.80, and 0.75 for CC, TC, and TT genotype, respectively; G × E p-value: 1.8 × 10-7). Rs1620977 in the NEGR1 gene showed an association with fruit intake (p-value: 1.0 × 10-8) and G × fruit interaction with CRC (OR per quartile of fruit increase = 0.75, 0.65, and 0.56 for AA, AG, and GG genotype, respectively; G × E -p-value: 0.029). INTERPRETATION: We identified 2 loci associated with fibre and fruit intake that also modify the association of these dietary factors with CRC risk. Potential mechanisms include chronic inflammatory intestinal disorders, and gut function. However, further studies are needed for mechanistic validation and replication of findings. FUNDING: National Institutes of Health, National Cancer Institute. Full funding details for the individual consortia are provided in acknowledgments.

Two genome-wide interaction loci modify the association of nonsteroidal anti-inflammatory drugs with colorectal cancer.

Regular, long-term aspirin use may act synergistically with genetic variants, particularly those in mechanistically relevant pathways, to confer a protective effect on colorectal cancer (CRC) risk. We leveraged pooled data from 52 clinical trial, cohort, and case-control studies that included 30,806 CRC cases and 41,861 controls of European ancestry to conduct a genome-wide interaction scan between regular aspirin/nonsteroidal anti-inflammatory drug (NSAID) use and imputed genetic variants. After adjusting for multiple comparisons, we identified statistically significant interactions between regular aspirin/NSAID use and variants in 6q24.1 (top hit rs72833769), which has evidence of influencing expression of TBC1D7 (a subunit of the TSC1-TSC2 complex, a key regulator of MTOR activity), and variants in 5p13.1 (top hit rs350047), which is associated with expression of PTGER4 (codes a cell surface receptor directly involved in the mode of action of aspirin). Genetic variants with functional impact may modulate the chemopreventive effect of regular aspirin use, and our study identifies putative previously unidentified targets for additional mechanistic interrogation.

Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes.

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

Novel insights into genetic susceptibility for colorectal cancer from transcriptome-wide association and functional investigation.

BACKGROUND: Transcriptome-wide association studies have been successful in identifying candidate susceptibility genes for colorectal cancer (CRC). To strengthen susceptibility gene discovery, we conducted a large transcriptome-wide association study and an alternative splicing transcriptome-wide association study in CRC using improved genetic prediction models and performed in-depth functional investigations. METHODS: We analyzed RNA-sequencing data from normal colon tissues and genotype data from 423 European descendants to build genetic prediction models of gene expression and alternative splicing and evaluated model performance using independent RNA-sequencing data from normal colon tissues of the Genotype-Tissue Expression Project. We applied the verified models to genome-wide association studies (GWAS) summary statistics among 58 131 CRC cases and 67 347 controls of European ancestry to evaluate associations of genetically predicted gene expression and alternative splicing with CRC risk. We performed in vitro functional assays for 3 selected genes in multiple CRC cell lines. RESULTS: We identified 57 putative CRC susceptibility genes, which included the 48 genes from transcriptome-wide association studies and 15 genes from splicing transcriptome-wide association studies, at a Bonferroni-corrected P value less than .05. Of these, 16 genes were not previously implicated in CRC susceptibility, including a gene PDE7B (6q23.3) at locus previously not reported by CRC GWAS. Gene knockdown experiments confirmed the oncogenic roles for 2 unreported genes, TRPS1 and METRNL, and a recently reported gene, C14orf166. CONCLUSION: This study discovered new putative susceptibility genes of CRC and provided novel insights into the biological mechanisms underlying CRC development.

Genome-wide interaction analysis of folate for colorectal cancer risk.

BACKGROUND: Epidemiological and experimental evidence suggests that higher folate intake is associated with decreased colorectal cancer (CRC) risk; however, the mechanisms underlying this relationship are not fully understood. Genetic variation that may have a direct or indirect impact on folate metabolism can provide insights into folate's role in CRC. OBJECTIVES: Our aim was to perform a genome-wide interaction analysis to identify genetic variants that may modify the association of folate on CRC risk. METHODS: We applied traditional case-control logistic regression, joint 3-degree of freedom, and a 2-step weighted hypothesis approach to test the interactions of common variants (allele frequency >1%) across the genome and dietary folate, folic acid supplement use, and total folate in relation to risk of CRC in 30,550 cases and 42,336 controls from 51 studies from 3 genetic consortia (CCFR, CORECT, GECCO). RESULTS: Inverse associations of dietary, total folate, and folic acid supplement with CRC were found (odds ratio [OR]: 0.93; 95% confidence interval [CI]: 0.90, 0.96; and 0.91; 95% CI: 0.89, 0.94 per quartile higher intake, and 0.82 (95% CI: 0.78, 0.88) for users compared with nonusers, respectively). Interactions (P-interaction < 5×10-8) of folic acid supplement and variants in the 3p25.2 locus (in the region of Synapsin II [SYN2]/tissue inhibitor of metalloproteinase 4 [TIMP4]) were found using traditional interaction analysis, with variant rs150924902 (located upstream to SYN2) showing the strongest interaction. In stratified analyses by rs150924902 genotypes, folate supplementation was associated with decreased CRC risk among those carrying the TT genotype (OR: 0.82; 95% CI: 0.79, 0.86) but increased CRC risk among those carrying the TA genotype (OR: 1.63; 95% CI: 1.29, 2.05), suggesting a qualitative interaction (P-interaction = 1.4×10-8). No interactions were observed for dietary and total folate. CONCLUSIONS: Variation in 3p25.2 locus may modify the association of folate supplement with CRC risk. Experimental studies and studies incorporating other relevant omics data are warranted to validate this finding.

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.

Assessment of Colorectal Cancer Risk Factors through the Application of Network-Based Approaches in a Racially Diverse Cohort of Colon Organoid Stem Cells.

Numerous demographic factors have been associated with colorectal cancer (CRC) risk. To better define biological mechanisms underlying these associations, we performed RNA sequencing of stem-cell-enriched organoids derived from the healthy colons of seven European Americans and eight African Americans. A weighted gene co-expression network analysis was performed following RNA sequencing. Module-trait relationships were determined through the association testing of each module and five CRC risk factors (age, body mass index, sex, smoking history, and race). Only modules that displayed a significantly positive correlation for gene significance and module membership were considered for further investigation. In total, 16 modules were associated with known CRC risk factors (p < 0.05). To contextualize the role of risk modules in CRC, publicly available RNA-sequencing data from TCGA-COAD were downloaded and re-analyzed. Differentially expressed genes identified between tumors and matched normal-adjacent tissue were overlaid across each module. Loci derived from CRC genome-wide association studies were additionally overlaid across modules to identify robust putative targets of risk. Among them, MYBL2 and RXRA represented strong plausible drivers through which cigarette smoking and BMI potentially modulated CRC risk, respectively. In summary, our findings highlight the potential of the colon organoid system in identifying novel CRC risk mechanisms in an ancestrally diverse and cellularly relevant population.

Probing the diabetes and colorectal cancer relationship using gene - environment interaction analyses.

BACKGROUND: Diabetes is an established risk factor for colorectal cancer. However, the mechanisms underlying this relationship still require investigation and it is not known if the association is modified by genetic variants. To address these questions, we undertook a genome-wide gene-environment interaction analysis. METHODS: We used data from 3 genetic consortia (CCFR, CORECT, GECCO; 31,318 colorectal cancer cases/41,499 controls) and undertook genome-wide gene-environment interaction analyses with colorectal cancer risk, including interaction tests of genetics(G)xdiabetes (1-degree of freedom; d.f.) and joint testing of Gxdiabetes, G-colorectal cancer association (2-d.f. joint test) and G-diabetes correlation (3-d.f. joint test). RESULTS: Based on the joint tests, we found that the association of diabetes with colorectal cancer risk is modified by loci on chromosomes 8q24.11 (rs3802177, SLC30A8 - ORAA: 1.62, 95% CI: 1.34-1.96; ORAG: 1.41, 95% CI: 1.30-1.54; ORGG: 1.22, 95% CI: 1.13-1.31; p-value3-d.f.: 5.46 × 10-11) and 13q14.13 (rs9526201, LRCH1 - ORGG: 2.11, 95% CI: 1.56-2.83; ORGA: 1.52, 95% CI: 1.38-1.68; ORAA: 1.13, 95% CI: 1.06-1.21; p-value2-d.f.: 7.84 × 10-09). DISCUSSION: These results suggest that variation in genes related to insulin signaling (SLC30A8) and immune function (LRCH1) may modify the association of diabetes with colorectal cancer risk and provide novel insights into the biology underlying the diabetes and colorectal cancer relationship.

Multi-omic analysis in normal colon organoids highlights MSH4 as a novel marker of defective mismatch repair in Lynch syndrome and microsatellite instability.

INTRODUCTION: Lynch syndrome (LS) is a hereditary condition that increases the risk of colorectal (CRC) and extracolonic cancers that exhibit microsatellite instability (MSI-H). MSI-H is driven by defective mismatch repair (dMMR), and approximately 15% of nonhereditary CRCs also exhibit MSI-H. Here, we aimed to better define mechanisms underlying tumor initiation in LS and MSI-H cancers through multi-omic analyses of LS normal colon organoids and MSI-H tumors. METHODS: Right (n = 35) and left (n = 23) colon organoids generated from normal colon biopsies at routine colonoscopy of LS and healthy individuals were subjected to Illumina EPIC array. Differentially methylated region (DMR) analysis was performed by DMRcate. RNA-sequencing (n = 16) and bisulfite-sequencing (n = 15) were performed on a subset of right colon organoids. CRISPR-cas9-mediated editing of MMR genes in colon organoids of healthy individuals was followed by quantitative PCR of MSH4. The relationship between MSH4 expression and tumor mutational burden was further explored in three independent tumor data sets. RESULTS: We identified a hypermethylated region of MSH4 in both the right and left colon organoids of LS versus healthy controls, which we validated using bisulfite-sequencing. DMR analysis in three gastrointestinal and one endometrial data set revealed that this region was also hypermethylated in MSI-H versus microsatellite stable (MSS) tumors. MSH4 expression was increased in colon organoids of LS versus healthy subjects and in publicly available MSI-H versus MSS tumors across four RNA-seq and four microarray data sets. CRISPR-cas9 editing of MLH1 and MSH2, but not MSH6, in normal colon organoids significantly increased MSH4 expression. MSH4 expression was significantly associated with tumor mutational burden in three publicly available data sets. CONCLUSIONS: Our findings implicate DNA methylation and gene expression differences of MSH4 as a marker of dMMR and as a potential novel biomarker of LS. Our study of LS colon organoids supports the hypothesis that dMMR exists in the colons of LS subjects prior to CRC.

A Genetic Locus within the FMN1/GREM1 Gene Region Interacts with Body Mass Index in Colorectal Cancer Risk.

Colorectal cancer risk can be impacted by genetic, environmental, and lifestyle factors, including diet and obesity. Gene-environment interactions (G × E) can provide biological insights into the effects of obesity on colorectal cancer risk. Here, we assessed potential genome-wide G × E interactions between body mass index (BMI) and common SNPs for colorectal cancer risk using data from 36,415 colorectal cancer cases and 48,451 controls from three international colorectal cancer consortia (CCFR, CORECT, and GECCO). The G × E tests included the conventional logistic regression using multiplicative terms (one degree of freedom, 1DF test), the two-step EDGE method, and the joint 3DF test, each of which is powerful for detecting G × E interactions under specific conditions. BMI was associated with higher colorectal cancer risk. The two-step approach revealed a statistically significant G×BMI interaction located within the Formin 1/Gremlin 1 (FMN1/GREM1) gene region (rs58349661). This SNP was also identified by the 3DF test, with a suggestive statistical significance in the 1DF test. Among participants with the CC genotype of rs58349661, overweight and obesity categories were associated with higher colorectal cancer risk, whereas null associations were observed across BMI categories in those with the TT genotype. Using data from three large international consortia, this study discovered a locus in the FMN1/GREM1 gene region that interacts with BMI on the association with colorectal cancer risk. Further studies should examine the potential mechanisms through which this locus modifies the etiologic link between obesity and colorectal cancer. SIGNIFICANCE: This gene-environment interaction analysis revealed a genetic locus in FMN1/GREM1 that interacts with body mass index in colorectal cancer risk, suggesting potential implications for precision prevention strategies.

Genome-wide Interaction Study with Smoking for Colorectal Cancer Risk Identifies Novel Genetic Loci Related to Tumor Suppression, Inflammation, and Immune Response.

BACKGROUND: Tobacco smoking is an established risk factor for colorectal cancer. However, genetically defined population subgroups may have increased susceptibility to smoking-related effects on colorectal cancer. METHODS: A genome-wide interaction scan was performed including 33,756 colorectal cancer cases and 44,346 controls from three genetic consortia. RESULTS: Evidence of an interaction was observed between smoking status (ever vs. never smokers) and a locus on 3p12.1 (rs9880919, P = 4.58 × 10-8), with higher associated risk in subjects carrying the GG genotype [OR, 1.25; 95% confidence interval (CI), 1.20-1.30] compared with the other genotypes (OR <1.17 for GA and AA). Among ever smokers, we observed interactions between smoking intensity (increase in 10 cigarettes smoked per day) and two loci on 6p21.33 (rs4151657, P = 1.72 × 10-8) and 8q24.23 (rs7005722, P = 2.88 × 10-8). Subjects carrying the rs4151657 TT genotype showed higher risk (OR, 1.12; 95% CI, 1.09-1.16) compared with the other genotypes (OR <1.06 for TC and CC). Similarly, higher risk was observed among subjects carrying the rs7005722 AA genotype (OR, 1.17; 95% CI, 1.07-1.28) compared with the other genotypes (OR <1.13 for AC and CC). Functional annotation revealed that SNPs in 3p12.1 and 6p21.33 loci were located in regulatory regions, and were associated with expression levels of nearby genes. Genetic models predicting gene expression revealed that smoking parameters were associated with lower colorectal cancer risk with higher expression levels of CADM2 (3p12.1) and ATF6B (6p21.33). CONCLUSIONS: Our study identified novel genetic loci that may modulate the risk for colorectal cancer of smoking status and intensity, linked to tumor suppression and immune response. IMPACT: These findings can guide potential prevention treatments.

Insights into Early Onset Colorectal Cancer through Analysis of Normal Colon Organoids of Familial Adenomatous Polyposis Patients.

Early onset colorectal cancer (EOCRC) rates have increased in recent decades. While lowering the recommended age for routine colonoscopies to 45 may reduce this burden, such measures do not address those who develop CRC before that age. Additional measures are needed to identify individuals at-risk for CRC. To better define transcriptomic events that precede the development of CRC, we performed RNA-sequencing analysis in colon organoids derived from seven healthy and six familial adenomatous polyposis (FAP) patients. This led to the identification of 2635 significant differentially expressed genes (FDR < 0.05). Through secondary analysis of publicly available datasets, we found that these genes were enriched for significant genes also present in FAP CRC and non-hereditary CRC datasets, including a subset that were unique to EOCRC. By exposing FAP colon organoids to a three-day ethanol treatment, we found that two EOCRC-relevant genes were also targets of CRC related lifestyle factors. Our data provides unique insight into the potential, early mechanisms of CRC development in colon epithelial cells, which may provide biomarkers for patient monitoring. We also show how modifiable lifestyle factors may further alter genes relevant to EOCRC, adding weight to the hypothesis that such factors represent an important contributor to increased EOCRC incidence.

Transcriptomic Response to Calcium in Normal Colon Organoids is Impacted by Colon Location and Sex.

Observational studies indicate that calcium supplementation may protect against colorectal cancer. Stratified analyses suggest that this protective effect may differ based on anatomic subsite and sex, but these hypotheses have been difficult to test experimentally. Here, we exposed 36 patient-derived organoid lines derived from normal colon biopsies (21 right colons, 15 left colons) of unrelated subjects (18 female, 18 male) to moderate (1.66 mmol/L) or high (5.0 mmol/L) concentrations of calcium for 72 hours. We performed bulk RNA-sequencing to measure gene expression, and cell composition was inferred using single-cell deconvolution in CIBERSORTx. We tested for significant differences in gene expression using generalized linear models in DESeq2. Exposure to higher levels of calcium was associated with changes in cell composition (P < 0.05), most notably increased goblet and reduced stem cell populations, and differential expression of 485 genes (FDR < 0.05). We found that 40 of these differentially expressed genes mapped to genomic loci identified through colorectal cancer genome-wide association studies, suggesting a potential biologic overlap between calcium supplementation and inherited colorectal cancer risk. Stratified analyses identified more differentially expressed genes in colon organoids derived from right sided colon and male subjects than those derived from left sided colon and female subjects. We confirmed the presence of a stronger right-sided effect for one of these genes, HSD17B2 using qPCR in a subset of matched right and left colon organoids (n = 4). By relating our findings to genetic data, we provide new insights into how nutritional and genetic factors may interact to influence colorectal cancer risk. PREVENTION RELEVANCE: A chemopreventive role for calcium in colorectal cancer is still unclear. Here, we identify mechanisms through which calcium supplementation may reduce risk. Calcium supplementation increased differentiation and altered expression of colorectal cancer-related genes in a large study of patient-derived colon organoids. These findings were influenced by colon location and sex.

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.

Transcriptome-Wide Association Study for Inflammatory Bowel Disease Reveals Novel Candidate Susceptibility Genes in Specific Colon Subsites and Tissue Categories.

BACKGROUND AND AIMS: Genome-wide association studies [GWAS] for inflammatory bowel disease [IBD] have identified 240 risk variants. However, the benefit of understanding the genetic architecture of IBD remains to be exploited. Transcriptome-wide association studies [TWAS] associate gene expression with genetic susceptibility to disease, providing functional insight into risk loci. In this study, we integrate relevant datasets for IBD and perform a TWAS to nominate novel genes implicated in IBD genetic susceptibility. METHODS: We applied elastic net regression to generate gene expression prediction models for the University of Barcelona and University of Virginia RNA sequencing project [BarcUVa-Seq] and correlated expression and disease association research [CEDAR] datasets. Together with Genotype-Tissue Expression project [GTEx] data, and GWAS results from about 60 000 individuals, we employed Summary-PrediXcan and Summary-MultiXcan for single and joint analyses of TWAS results, respectively. RESULTS: BarcUVa-Seq TWAS revealed 39 novel genes whose expression in the colon is associated with IBD genetic susceptibility. They included expression markers for specific colon cell types. TWAS meta-analysis including all tissues/cell types provided 186 novel candidate susceptibility genes. Additionally, we identified 78 novel susceptibility genes whose expression is associated with IBD exclusively in immune (N = 19), epithelial (N = 25), mesenchymal (N = 22) and neural (N = 12) tissue categories. Associated genes were involved in relevant molecular pathways, including pathways related to known IBD therapeutics, such as tumour necrosis factor signalling. CONCLUSION: These findings provide insight into tissue-specific molecular processes underlying IBD genetic susceptibility. Associated genes could be candidate targets for new therapeutics and should be prioritized in functional studies.

Genome-wide characterization of mitochondrial DNA methylation in human brain.

Background: There is growing interest in the role of DNA methylation in regulating the transcription of mitochondrial genes, particularly in brain disorders characterized by mitochondrial dysfunction. Here, we present a novel approach to interrogate the mitochondrial DNA methylome at single base resolution using targeted bisulfite sequencing. We applied this method to investigate mitochondrial DNA methylation patterns in post-mortem superior temporal gyrus and cerebellum brain tissue from seven human donors. Results: We show that mitochondrial DNA methylation patterns are relatively low but conserved, with peaks in DNA methylation at several sites, such as within the D-LOOP and the genes MT-ND2, MT-ATP6, MT-ND4, MT-ND5 and MT-ND6, predominantly in a non-CpG context. The elevated DNA methylation we observe in the D-LOOP we validate using pyrosequencing. We identify loci that show differential DNA methylation patterns associated with age, sex and brain region. Finally, we replicate previously reported differentially methylated regions between brain regions from a methylated DNA immunoprecipitation sequencing study. Conclusions: We have annotated patterns of DNA methylation at single base resolution across the mitochondrial genome in human brain samples. Looking to the future this approach could be utilized to investigate the role of mitochondrial epigenetic mechanisms in disorders that display mitochondrial dysfunction.

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.

Novel insights into the molecular mechanisms underlying risk of colorectal cancer from smoking and red/processed meat carcinogens by modeling exposure in normal colon organoids.

Tobacco smoke and red/processed meats are well-known risk factors for colorectal cancer (CRC). Most research has focused on studies of normal colon biopsies in epidemiologic studies or treatment of CRC cell lines in vitro. These studies are often constrained by challenges with accuracy of self-report data or, in the case of CRC cell lines, small sample sizes and lack of relationship to normal tissue at risk. In an attempt to address some of these limitations, we performed a 24-hour treatment of a representative carcinogens cocktail in 37 independent organoid lines derived from normal colon biopsies. Machine learning algorithms were applied to bulk RNA-sequencing and revealed cellular composition changes in colon organoids. We identified 738 differentially expressed genes in response to carcinogens exposure. Network analysis identified significantly different modules of co-expression, that included genes related to MSI-H tumor biology, and genes previously implicated in CRC through genome-wide association studies. Our study helps to better define the molecular effects of representative carcinogens from smoking and red/processed meat in normal colon epithelial cells and in the etiology of the MSI-H subtype of CRC, and suggests an overlap between molecular mechanisms involved in inherited and environmental CRC risk.

Transcriptome-wide In Vitro Effects of Aspirin on Patient-derived Normal Colon Organoids.

Mechanisms underlying aspirin chemoprevention of colorectal cancer remain unclear. Prior studies have been limited because of the inability of preclinical models to recapitulate human normal colon epithelium or cellular heterogeneity present in mucosal biopsies. To overcome some of these obstacles, we performed in vitro aspirin treatment of colon organoids derived from normal mucosal biopsies to reveal transcriptional networks relevant to aspirin chemoprevention. Colon organoids derived from 38 healthy individuals undergoing endoscopy were treated with 50 µmol/L aspirin or vehicle control for 72 hours and subjected to bulk RNA sequencing. Paired regression analysis using DESeq2 identified differentially expressed genes (DEG) associated with aspirin treatment. Cellular composition was determined using CIBERSORTx. Aspirin treatment was associated with 1,154 significant (q < 0.10) DEGs prior to deconvolution. We provide replication of these findings in an independent population-based RNA-sequencing dataset of mucosal biopsies (BarcUVa-Seq), where a significant enrichment for overlap of DEGs was observed (P < 2.2E-16). Single-cell deconvolution revealed changes in cell composition, including a decrease in transit-amplifying cells following aspirin treatment (P = 0.01). Following deconvolution, DEGs included novel putative targets for aspirin such as TRABD2A (q = 0.055), a negative regulator of Wnt signaling. Weighted gene co-expression network analysis identified 12 significant modules, including two that contained hubs for EGFR and PTGES2, the latter being previously implicated in aspirin chemoprevention. In summary, aspirin treatment of patient-derived colon organoids using physiologically relevant doses resulted in transcriptome-wide changes that reveal altered cell composition and improved understanding of transcriptional pathways, providing novel insight into its chemopreventive properties. PREVENTION RELEVANCE: Numerous studies have highlighted a role for aspirin in colorectal cancer chemoprevention, though the mechanisms driving this association remain unclear. We addressed this by showing that aspirin treatment of normal colon organoids diminished the transit-amplifying cell population, inhibited prostaglandin synthesis, and dysregulated expression of novel genes implicated in colon tumorigenesis.