Discovery of common and rare genetic risk variants for colorectal cancer

Casual associations between blood metabolites and colon cancer

BACKGROUND

Limited knowledge exists regarding the casual associations linking blood metabolites and the risk of developing colorectal cancer.

AIM

To investigate causal associations between blood metabolites and colon cancer.

METHODS

The study utilized a two-sample Mendelian randomization (MR) analysis to investigate the causal impact of 486 blood metabolites on colorectal cancer. The primary method of analysis used was the inverse variance weighted model. To further validate the results several sensitivity analyses were performed, including Cochran's Q test, MR-Egger intercept test, and MR robust adjusted profile score. These additional analyses were conducted to ensure the reliability and robustness of the findings.

RESULTS

After rigorous selection for genetic variation, 486 blood metabolites were included in the MR analysis. We found Mannose [odds ratio (OR) = 2.09 (1.10-3.97), P = 0.024], N-acetylglycine [OR = 3.14 (1.78-5.53), P = 7.54 × 10^-8], X-11593-O-methylascorbate [OR = 1.68 (1.04-2.72), P = 0.034], 1-arachidonoylglycerophosphocholine [OR = 4.23 (2.51-7.12), P = 6.35 × 10^-8] and 1-arachidonoylglycerophosphoethanolamine 4 [OR = 3.99 (1.17-13.54), P = 0.027] were positively causally associated with colorectal cancer, and we also found a negative causal relationship between Tyrosine [OR = 0.08 (0.01-0.63), P = 0.014], Urate [OR = 0.25 (0.10-0.62), P = 0.003], N-acetylglycine [0.73 (0.54-0.98), P = 0.033], X-12092 [OR = 0.89 (0.81-0.99), P = 0.028], Succinylcarnitine [OR = 0.48 (0.27-0.84), P = 0.09] with colorectal cancer. A series of sensitivity analyses were performed to confirm the rigidity of the results.

CONCLUSION

This study showed a causal relationship between 10 blood metabolites and colorectal cancer, of which 5 blood metabolites were found to be causal for the development of colorectal cancer and were confirmed as risk factors. The other five blood metabolites are protective factors.

Colorectal cancer initiation: Understanding early-stage disease for intervention

How tumors arise or the cause of precancerous lesions is a fundamental question in cancer biology. It is generally accepted that tumors originate from normal cells that undergo uncontrolled proliferation owing to genetic alterations. At the onset of adenoma formation, cancer driver mutations confer clonal growth advantage, enabling mutant cells to outcompete and eliminate the surrounding healthy cells. Hence, the development of precancerous lesions is not only attributed to the expansion of pre-malignant clones, but also relies on the relative fitness of mutated cells compared to the neighboring cells. Colorectal cancer (CRC) is an excellent model to investigate cancer origin as it follows a stereotypical process from mutant cell hyperplasia to adenoma formation and progression. Here, we review the evolving understanding of colonic tumor development, focusing on how cell intrinsic and extrinsic factors impact cell competition and the "clone war" between cancer-initiating cells and normal stem cells. We also discuss the promises and limitations of targeting cell competitiveness in cancer prevention and early intervention. The field of tumor initiation is currently in its infancy, elucidating the adenoma origin is crucial for designing effective prevention strategies and early treatments before cancer becomes incurable.

Cellular senescence gene TACC3 associated with colorectal cancer risk via genetic and DNA methylated alteration

Cell senescence genes play a vital role in the pathogenesis of colorectal cancer, a process that may involve the triggering of genetic variations and reversible phenotypes caused by epigenetic modifications. However, the specific regulatory mechanisms remain unclear. Using CellAge and The Cancer Genome Atlas databases and in-house RNA-seq data, DNA methylation-modified cellular senescence genes (DMCSGs) were validated by Support Vector Machine and correlation analyses. In 1150 cases and 1342 controls, we identified colorectal cancer risk variants in DMCSGs. The regulatory effects of gene, variant, and DNA methylation were explored through dual-luciferase and 5-azacytidine treatment experiments, complemented by multiple database analyses. Biological functions of key gene were evaluated via cell proliferation assays, SA-β-gal staining, senescence marker detection, and immune infiltration analyses. The genetic variant rs4558926 in the downstream of TACC3 was significantly associated with colorectal cancer risk (OR = 1.35, P = 3.22 × 10-4). TACC3 mRNA expression increased due to rs4558926 C > G and decreased DNA methylation levels. The CpG sites in the TACC3 promoter region were regulated by rs4558926. TACC3 knockdown decreased proliferation and senescence in colorectal cancer cells. In addition, subjects with high-TACC3 expression presented an immunosuppressive microenvironment. These findings provide insights into the involvement of genetic variants of cellular senescence genes in the development and progression of colorectal cancer.

Colorectal Cancer Risk between Mendelian and Non-Mendelian Inheritance

Hereditary colorectal cancer has been an area of focus for research and public health practitioners due to our ability to quantify risk and then act based on such results by enrolling patients in surveillance programs. The wide access to genetic testing and whole-genome sequencing has resulted in identifying many low/moderate penetrance genes. Above all, our understanding of the family component of colorectal cancer has been improving. Polygenic scores are becoming part of the risk assessment for many cancers, and the data about polygenic risk scores for colorectal cancer is promising. The challenge is determining how we incorporate this data in clinical care.

Allometric versus traditional body-shape indices and risk of colorectal cancer: a Mendelian randomization analysis

BACKGROUND

Traditional body-shape indices such as Waist Circumference (WC), Hip Circumference (HC), and Waist-to-Hip Ratio (WHR) are associated with colorectal cancer (CRC) risk, but are correlated with Body Mass Index (BMI), and adjustment for BMI introduces a strong correlation with height. Thus, new allometric indices have been developed, namely A Body Shape Index (ABSI), Hip Index (HI), and Waist-to-Hip Index (WHI), which are uncorrelated with weight and height; these have also been associated with CRC risk in observational studies, but information from Mendelian randomization (MR) studies is missing.

METHODS

We used two-sample MR to examine potential causal cancer site- and sex-specific associations of the genetically-predicted allometric body-shape indices with CRC risk, and compared them with BMI-adjusted traditional body-shape indices, and BMI. Data were obtained from UK Biobank and the GIANT consortium, and from GECCO, CORECT and CCFR consortia.

RESULTS

WHI was positively associated with CRC in men (OR per SD: 1.20, 95% CI: 1.03-1.39) and in women (1.15, 1.06-1.24), and similarly for colon and rectal cancer. ABSI was positively associated with colon and rectal cancer in men (1.27, 1.03-1.57; and 1.40, 1.10-1.77, respectively), and with colon cancer in women (1.20, 1.07-1.35). There was little evidence for association between HI and colon or rectal cancer. The BMI-adjusted WHR and HC showed similar associations to WHI and HI, whereas WC showed similar associations to ABSI only in women.

CONCLUSIONS

This large MR study provides strong evidence for a potential causal positive association of the allometric indices ABSI and WHI with CRC in both sexes, thus establishing the association between abdominal fat and CRC without the limitations of the traditional waist size indices and independently of BMI. Among the BMI-adjusted traditional indices, WHR and HC provided equivalent associations with WHI and HI, while differences were observed between WC and ABSI.

Colon Polyp Surveillance: Separating the Wheat From the Chaff

One goal of colorectal cancer (CRC) screening is to prevent CRC incidence by removing precancerous colonic polyps, which are detected in up to 50% of screening examinations. Yet, the lifetime risk of CRC is 3.9%-4.3%, so it is clear that most of these individuals with polyps would not develop CRC in their lifetime. It is, therefore, a challenge to determine which individuals with polyps will benefit from follow-up, and at what intervals. There is some evidence that individuals with advanced polyps, based on size and histology, benefit from intensive surveillance. However, a large proportion of individuals will have small polyps without advanced histologic features (ie, "nonadvanced"), where the benefits of surveillance are uncertain and controversial. Demand for surveillance will further increase as more polyps are detected due to increased screening uptake, recent United States recommendations to expand screening to younger individuals, and emergence of polyp detection technology. We review the current understanding and clinical implications of the natural history, biology, and outcomes associated with various categories of colon polyps based on size, histology, and number. Our aims are to highlight key knowledge gaps, specifically focusing on certain categories of polyps that may not be associated with future CRC risk, and to provide insights to inform research priorities and potential management strategies. Optimization of CRC prevention programs based on updated knowledge about the future risks associated with various colon polyps is essential to ensure cost-effective screening and surveillance, wise use of resources, and inform efforts to personalize recommendations.

Characterization of cell cycle, inflammation, and oxidative stress signaling role in non-communicable diseases: Insights into genetic variants, microRNAs and pathways

Non-Communicable Diseases (NCDs) significantly impact global health, contributing to over 70% of premature deaths, as reported by the World Health Organization (WHO). These diseases have complex and multifactorial origins, involving genetic, epigenetic, environmental and lifestyle factors. While Genome-Wide Association Study (GWAS) is widely recognized as a valuable tool for identifying variants associated with complex phenotypes; the multifactorial nature of NCDs necessitates a more comprehensive exploration, encompassing not only the genetic but also the epigenetic aspect. For this purpose, we employed a bioinformatics-multiomics approach to examine the genetic and epigenetic characteristics of NCDs (i.e. colorectal cancer, coronary atherosclerosis, squamous cell lung cancer, psoriasis, type 2 diabetes, and multiple sclerosis), aiming to identify novel biomarkers for diagnosis and prognosis. Leveraging GWAS summary statistics, we pinpointed Single Nucleotide Polymorphisms (SNPs) independently associated with each NCD. Subsequently, we identified genes linked to cell cycle, inflammation and oxidative stress mechanisms, revealing shared genes across multiple diseases, suggesting common functional pathways. From an epigenetic perspective, we identified microRNAs (miRNAs) with regulatory functions targeting these genes of interest. Our findings underscore critical genetic pathways implicated in these diseases. In colorectal cancer, the dysregulation of the "Cytokine Signaling in Immune System" pathway, involving LAMA5 and SMAD7, regulated by Hsa-miR-21-5p, Hsa-miR-103a-3p, and Hsa-miR-195-5p, emerged as pivotal. In coronary atherosclerosis, the pathway associated with "binding of TCF/LEF:CTNNB1 to target gene promoters" displayed noteworthy implications, with the MYC factor controlled by Hsa-miR-16-5p as a potential regulatory factor. Squamous cell lung carcinoma analysis revealed significant pathways such as "PTK6 promotes HIF1A stabilization," regulated by Hsa-let-7b-5p. In psoriasis, the "Endosomal/Vacuolar pathway," involving HLA-C and Hsa-miR-148a-3p and Hsa-miR-148b-3p, was identified as crucial. Type 2 Diabetes implicated the "Regulation of TP53 Expression" pathway, controlled by Hsa-miR-106a-5p and Hsa-miR-106b-5p. In conclusion, our study elucidates the genetic framework and molecular mechanisms underlying NCDs, offering crucial insights into potential genetic/epigenetic biomarkers for diagnosis and prognosis. The specificity of pathways and related miRNAs in different pathologies highlights promising candidates for further clinical validation, with the potential to advance personalized treatments and alleviate the global burden of NCDs.

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.

Germline genetic regulation of the colorectal tumor immune microenvironment

OBJECTIVE

To evaluate the contribution of germline genetics to regulating the briskness and diversity of T cell responses in CRC, we conducted a genome-wide association study to examine the associations between germline genetic variation and quantitative measures of T cell landscapes in 2,876 colorectal tumors from participants in the Molecular Epidemiology of Colorectal Cancer Study (MECC).

METHODS

Germline DNA samples were genotyped and imputed using genome-wide arrays. Tumor DNA samples were extracted from paraffin blocks, and T cell receptor clonality and abundance were quantified by immunoSEQ (Adaptive Biotechnologies, Seattle, WA). Tumor infiltrating lymphocytes per high powered field (TILs/hpf) were scored by a gastrointestinal pathologist. Regression models were used to evaluate the associations between each variant and the three T-cell features, adjusting for sex, age, genotyping platform, and global ancestry. Three independent datasets were used for replication.

RESULTS

We identified a SNP (rs4918567) near RBM20 associated with clonality at a genome-wide significant threshold of 5 × 10- 8, with a consistent direction of association in both discovery and replication datasets. Expression quantitative trait (eQTL) analyses and in silico functional annotation for these loci provided insights into potential functional roles, including a statistically significant eQTL between the T allele at rs4918567 and higher expression of ADRA2A (P = 0.012) in healthy colon mucosa.

CONCLUSIONS

Our study suggests that germline genetic variation is associated with the quantity and diversity of adaptive immune responses in CRC. Further studies are warranted to replicate these findings in additional samples and to investigate functional genomic mechanisms.

Loss of Y in regulatory T lymphocytes in the tumor micro-environment of primary colorectal cancers and liver metastases

Male sex is a risk factor for colorectal cancer (CRC) with higher illness burden and earlier onset. Thus, we hypothesized that loss of chromosome Y (LOY) in the tumor micro-environment (TME) might be involved in oncogenesis. Previous studies show that LOY in circulating leukocytes of aging men was associated with shorter survival and non-hematological cancer, as well as higher LOY in CD4 + T-lymphocytes in men with prostate cancer vs. controls. However, nothing is known about LOY in leukocytes infiltrating TME and we address this aspect here. We studied frequency and functional effects of LOY in blood, TME and non-tumorous tissue. Regulatory T-lymphocytes (Tregs) in TME had the highest frequency of LOY (22%) in comparison to CD4 + T-lymphocytes and cytotoxic CD8 + T-lymphocytes. LOY score using scRNA-seq was also linked to higher expression of PDCD1, TIGIT and IKZF2 in Tregs. PDCD1 and TIGIT encode immune checkpoint receptors involved in the regulation of Tregs function. Our study sets the direction for further functional research regarding a probable role of LOY in intensifying features related to the suppressive phenotype of Tregs in TME and consequently a possible influence on immunotherapy response in CRC patients.

Tissue-specific genetic variation suggests distinct molecular pathways between body shape phenotypes and colorectal cancer

It remains unknown whether adiposity subtypes are differentially associated with colorectal cancer (CRC). To move beyond single-trait anthropometric indicators, we derived four multi-trait body shape phenotypes reflecting adiposity subtypes from principal components analysis on body mass index, height, weight, waist-to-hip ratio, and waist and hip circumference. A generally obese (PC1) and a tall, centrally obese (PC3) body shape were both positively associated with CRC risk in observational analyses in 329,828 UK Biobank participants (3728 cases). In genome-wide association studies in 460,198 UK Biobank participants, we identified 3414 genetic variants across four body shapes and Mendelian randomization analyses confirmed positive associations of PC1 and PC3 with CRC risk (52,775 cases/45,940 controls from GECCO/CORECT/CCFR). Brain tissue-specific genetic instruments, mapped to PC1 through enrichment analysis, were responsible for the relationship between PC1 and CRC, while the relationship between PC3 and CRC was predominantly driven by adipose tissue-specific genetic instruments. This study suggests distinct putative causal pathways between adiposity subtypes and CRC.

A genome-wide association study provides insights into the genetic etiology of 57 essential and non-essential trace elements in humans

Trace elements are important for human health but may exert toxic or adverse effects. Mechanisms of uptake, distribution, metabolism, and excretion are partly under genetic control but have not yet been extensively mapped. Here we report a comprehensive multi-element genome-wide association study of 57 essential and non-essential trace elements. We perform genome-wide association meta-analyses of 14 trace elements in up to 6564 Scandinavian whole blood samples, and genome-wide association studies of 43 trace elements in up to 2819 samples measured only in the Trøndelag Health Study (HUNT). We identify 11 novel genetic loci associated with blood concentrations of arsenic, cadmium, manganese, selenium, and zinc in genome-wide association meta-analyses. In HUNT, several genome-wide significant loci are also indicated for other trace elements. Using two-sample Mendelian randomization, we find several indications of weak to moderate effects on health outcomes, the most precise being a weak harmful effect of increased zinc on prostate cancer. However, independent validation is needed. Our current understanding of trace element-associated genetic variants may help establish consequences of trace elements on human health.

Genetic Risk, Health-Associated Lifestyle, and Risk of Early-onset Total Cancer and Breast Cancer

Importance

Early-onset cancer (diagnosed under 50 years of age) is associated with aggressive disease characteristics and its rising incidence is a global concern. The association between healthy lifestyle and early-onset cancer and whether it varies by common genetic variants is unknown.

Objective

To examine the associations between genetic risk, lifestyle, and risk of early-onset cancers.

Design Setting and Participants

We analyzed a prospective cohort of 66,308 white British participants who were under age 50 and free of cancer at baseline in the UK Biobank.

Exposures

Sex-specific composite total cancer polygenic risk scores (PRSs), a breast cancer-specific PRS, and sex-specific health-associated lifestyle scores (HLSs, which summarize smoking status, body mass index [males only], physical activity, alcohol consumption, and diet).

Main Outcomes and Measures

Hazard ratios (HRs) and 95% confidence intervals (CIs) for early-onset total and breast cancer.

Results

A total of 1,247 incident invasive early-onset cancer cases (female: 820, male: 427, breast: 386) were documented. In multivariable-adjusted analyses with 2-year latency, higher genetic risk (highest vs. lowest tertile of PRS) was associated with significantly increased risks of early-onset total cancer in females (HR, 95% CI: 1.85, 1.50-2.29) and males (1.94, 1.45-2.59) as well as early-onset breast cancer in females (3.06, 2.20-4.25). An unfavorable lifestyle (highest vs. lowest category of HLS) was associated with higher risk of total cancer and breast cancer in females across genetic risk categories; the association with total cancer was stronger in the highest genetic risk category than the lowest: HRs in females and men were 1.85 (1.02, 3.36), 3.27 (0.78, 13.72) in the highest genetic risk category and 1.15 (0.44, 2.98), 1.16 (0.39, 3.40) in the lowest.

Conclusions and Relevance

Both genetic and lifestyle factors were independently associated with early-onset total and breast cancer risk. Compared to those with low genetic risk, individuals with a high genetic risk may benefit more from adopting a healthy lifestyle in preventing early-onset cancer.

Epidemiologic Factors in Relation to Colorectal Cancer Risk and Survival by Genotoxic Colibactin Mutational Signature

BACKGROUND

The genotoxin colibactin causes a tumor single-base substitution (SBS) mutational signature, SBS88. It is unknown whether epidemiologic factors' association with colorectal cancer risk and survival differs by SBS88.

METHODS

Within the Genetic Epidemiology of Colorectal Cancer Consortium and Colon Cancer Family Registry, we measured SBS88 in 4,308 microsatellite stable/microsatellite instability low tumors. Associations of epidemiologic factors with colorectal cancer risk by SBS88 were assessed using multinomial regression (N = 4,308 cases, 14,192 controls; cohort-only cases N = 1,911), and with colorectal cancer-specific survival using Cox proportional hazards regression (N = 3,465 cases).

RESULTS

392 (9%) tumors were SBS88 positive. Among all cases, the highest quartile of fruit intake was associated with lower risk of SBS88-positive colorectal cancer than SBS88-negative colorectal cancer [odds ratio (OR) = 0.53, 95% confidence interval (CI) 0.37-0.76; OR = 0.75, 95% CI 0.66-0.85, respectively, Pheterogeneity = 0.047]. Among cohort studies, associations of body mass index (BMI), alcohol, and fruit intake with colorectal cancer risk differed by SBS88. BMI ≥30 kg/m2 was associated with worse colorectal cancer-specific survival among those SBS88-positive [hazard ratio (HR) = 3.40, 95% CI 1.47-7.84], but not among those SBS88-negative (HR = 0.97, 95% CI 0.78-1.21, Pheterogeneity = 0.066).

CONCLUSIONS

Most epidemiologic factors did not differ by SBS88 for colorectal cancer risk or survival. Higher BMI may be associated with worse colorectal cancer-specific survival among those SBS88-positive; however, validation is needed in samples with whole-genome or whole-exome sequencing available.

IMPACT

This study highlights the importance of identification of tumor phenotypes related to colorectal cancer and understanding potential heterogeneity for risk and survival.

Genetic risk impacts the association of menopausal hormone therapy with colorectal cancer risk

BACKGROUND

Menopausal hormone therapy (MHT), a common treatment to relieve symptoms of menopause, is associated with a lower risk of colorectal cancer (CRC). To inform CRC risk prediction and MHT risk-benefit assessment, we aimed to evaluate the joint association of a polygenic risk score (PRS) for CRC and MHT on CRC risk.

METHODS

We used data from 28,486 postmenopausal women (11,519 cases and 16,967 controls) of European descent. A PRS based on 141 CRC-associated genetic variants was modeled as a categorical variable in quartiles. Multiplicative interaction between PRS and MHT use was evaluated using logistic regression. Additive interaction was measured using the relative excess risk due to interaction (RERI). 30-year cumulative risks of CRC for 50-year-old women according to MHT use and PRS were calculated.

RESULTS

The reduction in odds ratios by MHT use was larger in women within the highest quartile of PRS compared to that in women within the lowest quartile of PRS (p-value = 2.7 × 10-8). At the highest quartile of PRS, the 30-year CRC risk was statistically significantly lower for women taking any MHT than for women not taking any MHT, 3.7% (3.3%-4.0%) vs 6.1% (5.7%-6.5%) (difference 2.4%, P-value = 1.83 × 10-14); these differences were also statistically significant but smaller in magnitude in the lowest PRS quartile, 1.6% (1.4%-1.8%) vs 2.2% (1.9%-2.4%) (difference 0.6%, P-value = 1.01 × 10-3), indicating 4 times greater reduction in absolute risk associated with any MHT use in the highest compared to the lowest quartile of genetic CRC risk.

CONCLUSIONS

MHT use has a greater impact on the reduction of CRC risk for women at higher genetic risk. These findings have implications for the development of risk prediction models for CRC and potentially for the consideration of genetic information in the risk-benefit assessment of MHT use.

Association between 19 medication use and risk of common cancers: A cross-sectional and Mendelian randomisation study

Background

Previous studies have yielded inconsistent results concerning drug use and the risk of cancers. We conducted a large-scale cross-sectional study and a two-sample Mendelian randomisation (MR) study to reveal the causal effect between the use of 19 medications and the risk of four common cancers (breast, lung, colorectal, and prostate).

Methods

We obtained information on medication use and cancer diagnosis from National Health and Nutrition Examination Survey participants. After propensity score matching, we conducted survey-weighted multivariate logistic regression and restricted cubic spline analysis to assess the observed correlation between medication use and cancer while adjusting for multiple covariates. We also performed MR analysis to investigate causality based on summary data from genome-wide association studies on medication use and cancers. We performed sensitivity analyses, replication analysis, genetic correlation analysis, and reverse MR analysis to improve the reliability of MR findings.

Results

We found that the use of agents acting on the renin-angiotensin system was associated with reduced risk of prostate cancer (odds ratio (OR) = 0.42; 95% confidence interval (CI) = 0.27-0.63, P < 0.001), and there was a nonlinear association of 'decrease-to-increase-to-decrease' (P < 0.0001). The random-effects inverse variance weighted (IVW) model-based primary MR analysis (OR = 0.94, 95% CI = 0.91-0.97, P = 0.0007) and replication MR analysis (OR = 0.90, 95% CI = 0.85-0.96, P = 0.0006) both provided robust evidence of the causality of genetic liability for the use of agents acting on the renin-angiotensin system on a decreased risk of prostate cancer.

Conclusions

Our study provides robust evidence that the use of drugs acting on the renin-angiotensin system can reduce prostate cancer risk. Given the high prevalence of prostate cancer, these findings have important implications for drug selection and prostate cancer prevention in patients with cardiovascular disease.

Artificially sweetened beverages consumption and risk of obesity-related cancers: a wide-angled Mendelian randomization study

Background

The impact of artificially sweetened beverages (ASBs) consumption on obesity-related cancers (ORCs) risk remains controversial. To address this challenging issue, this study employed wide-angle mendelian randomization (MR) analyses to explore the genetic causality between ASB consumption and the risk of ORCs, thereby effectively minimizing the impact of external confounders.

Methods

We conducted a suite of analyses encompassing univariable, multivariable, and two-step MR to evaluate causal associations between ASB consumption (samples = 85,852) and risk of ORCs (total samples = 2,974,770) using summary statistics from genome-wide association studies (GWAS). Total, direct, and intermediary effects were derived by performing inverse-variance weighted (IVW), MR-Egger, weighted mode, weighted median, and lasso method. Additionally, we performed an extensive range of sensitivity analyses to counteract the potential effects of confounders, heterogeneity, and pleiotropy, enhancing the robustness and reliability of the findings.

Results

Genetically predicted ASB consumption was positively associated with the risk of colorectal cancer (CRC, p = 0.011; OR: 6.879; 95% CI: 1.551, 30.512 by IVW) and breast cancer (p = 0.022; OR: 3.881; 95% CI: 2.023, 9.776 by IVW). Multivariable analysis yielded similar results. The results of the two-step MR unveiled that body mass index (BMI) assumes a pivotal role in mediating the association between ASB consumption and CRC risk (intermediary effect = 0.068, p = 0.024).

Conclusion

No causal connection exists between ASB consumption and the majority of ORCs, in addition to CRC and breast cancer. Additionally, our findings suggest that BMI might be a potential mediator in the association between ASB consumption and CRC.

Genome-wide polygenic risk scores for colorectal cancer have implications for risk-based screening

BACKGROUND

Hereditary factors, including single genetic variants and family history, can be used for targeting colorectal cancer (CRC) screening, but limited data exist on the impact of polygenic risk scores (PRS) on risk-based CRC screening.

METHODS

Using longitudinal health and genomics data on 453,733 Finnish individuals including 8801 CRC cases, we estimated the impact of a genome-wide CRC PRS on CRC screening initiation age through population-calibrated incidence estimation over the life course in men and women.

RESULTS

Compared to the cumulative incidence of CRC at age 60 in Finland (the current age for starting screening in Finland), a comparable cumulative incidence was reached 5 and 11 years earlier in persons with high PRS (80-99% and >99%, respectively), while those with a low PRS (< 20%) reached comparable incidence 7 years later. The PRS was associated with increased risk of post-colonoscopy CRC after negative colonoscopy (hazard ratio 1.76 per PRS SD, 95% CI 1.54-2.01). Moreover, the PRS predicted colorectal adenoma incidence and improved incident CRC risk prediction over non-genetic risk factors.

CONCLUSIONS

Our findings demonstrate that a CRC PRS can be used for risk stratification of CRC, with further research needed to optimally integrate the PRS into risk-based screening.

Body size and risk of colorectal cancer molecular defined subtypes and pathways: Mendelian randomization analyses

BACKGROUND

Obesity has been positively associated with most molecular subtypes of colorectal cancer (CRC); however, the magnitude and the causality of these associations is uncertain.

METHODS

We used Mendelian randomization (MR) to examine potential causal relationships between body size traits (body mass index [BMI], waist circumference, and body fat percentage) with risks of Jass classification types and individual subtypes of CRC (microsatellite instability [MSI] status, CpG island methylator phenotype [CIMP] status, BRAF and KRAS mutations). Summary data on tumour markers were obtained from two genetic consortia (CCFR, GECCO).

FINDINGS

A 1-standard deviation (SD:5.1 kg/m2) increment in BMI levels was found to increase risks of Jass type 1MSI-high,CIMP-high,BRAF-mutated,KRAS-wildtype (odds ratio [OR]: 2.14, 95% confidence interval [CI]: 1.46, 3.13; p-value = 9 × 10-5) and Jass type 2non-MSI-high,CIMP-high,BRAF-mutated,KRAS-wildtype CRC (OR: 2.20, 95% CI: 1.26, 3.86; p-value = 0.005). The magnitude of these associations was stronger compared with Jass type 4non-MSI-high,CIMP-low/negative,BRAF-wildtype,KRAS-wildtype CRC (p-differences: 0.03 and 0.04, respectively). A 1-SD (SD:13.4 cm) increment in waist circumference increased risk of Jass type 3non-MSI-high,CIMP-low/negative,BRAF-wildtype,KRAS-mutated (OR 1.73, 95% CI: 1.34, 2.25; p-value = 9 × 10-5) that was stronger compared with Jass type 4 CRC (p-difference: 0.03). A higher body fat percentage (SD:8.5%) increased risk of Jass type 1 CRC (OR: 2.59, 95% CI: 1.49, 4.48; p-value = 0.001), which was greater than Jass type 4 CRC (p-difference: 0.03).

INTERPRETATION

Body size was more strongly linked to the serrated (Jass types 1 and 2) and alternate (Jass type 3) pathways of colorectal carcinogenesis in comparison to the traditional pathway (Jass type 4).

FUNDING

Cancer Research UK, National Institute for Health Research, Medical Research Council, National Institutes of Health, National Cancer Institute, American Institute for Cancer Research, Brigham and Women's Hospital, Prevent Cancer Foundation, Victorian Cancer Agency, Swedish Research Council, Swedish Cancer Society, Region Västerbotten, Knut and Alice Wallenberg Foundation, Lion's Cancer Research Foundation, Insamlingsstiftelsen, Umeå University. Full funding details are provided in acknowledgements.

Polygenic risk scores, radiation treatment exposures and subsequent cancer risk in childhood cancer survivors

Survivors of childhood cancer are at increased risk for subsequent cancers attributable to the late effects of radiotherapy and other treatment exposures; thus, further understanding of the impact of genetic predisposition on risk is needed. Combining genotype data for 11,220 5-year survivors from the Childhood Cancer Survivor Study and the St Jude Lifetime Cohort, we found that cancer-specific polygenic risk scores (PRSs) derived from general population, genome-wide association study, cancer loci identified survivors of European ancestry at increased risk of subsequent basal cell carcinoma (odds ratio per s.d. of the PRS: OR = 1.37, 95% confidence interval (CI) = 1.29-1.46), female breast cancer (OR = 1.42, 95% CI = 1.27-1.58), thyroid cancer (OR = 1.48, 95% CI = 1.31-1.67), squamous cell carcinoma (OR = 1.20, 95% CI = 1.00-1.44) and melanoma (OR = 1.60, 95% CI = 1.31-1.96); however, the association for colorectal cancer was not significant (OR = 1.19, 95% CI = 0.94-1.52). An investigation of joint associations between PRSs and radiotherapy found more than additive increased risks of basal cell carcinoma, and breast and thyroid cancers. For survivors with radiotherapy exposure, the cumulative incidence of subsequent cancer by age 50 years was increased for those with high versus low PRS. These findings suggest a degree of shared genetic etiology for these malignancy types in the general population and survivors, which remains evident in the context of strong radiotherapy-related risk.

Genome-Wide Gene-Environment Interaction Analyses to Understand the Relationship between Red Meat and Processed Meat Intake and Colorectal Cancer Risk

BACKGROUND

High red meat and/or processed meat consumption are established colorectal cancer risk factors. We conducted a genome-wide gene-environment (GxE) interaction analysis to identify genetic variants that may modify these associations.

METHODS

A pooled sample of 29,842 colorectal cancer cases and 39,635 controls of European ancestry from 27 studies were included. Quantiles for red meat and processed meat intake were constructed from harmonized questionnaire data. Genotyping arrays were imputed to the Haplotype Reference Consortium. Two-step EDGE and joint tests of GxE interaction were utilized in our genome-wide scan.

RESULTS

Meta-analyses confirmed positive associations between increased consumption of red meat and processed meat with colorectal cancer risk [per quartile red meat OR = 1.30; 95% confidence interval (CI) = 1.21-1.41; processed meat OR = 1.40; 95% CI = 1.20-1.63]. Two significant genome-wide GxE interactions for red meat consumption were found. Joint GxE tests revealed the rs4871179 SNP in chromosome 8 (downstream of HAS2); greater than median of consumption ORs = 1.38 (95% CI = 1.29-1.46), 1.20 (95% CI = 1.12-1.27), and 1.07 (95% CI = 0.95-1.19) for CC, CG, and GG, respectively. The two-step EDGE method identified the rs35352860 SNP in chromosome 18 (SMAD7 intron); greater than median of consumption ORs = 1.18 (95% CI = 1.11-1.24), 1.35 (95% CI = 1.26-1.44), and 1.46 (95% CI = 1.26-1.69) for CC, CT, and TT, respectively.

CONCLUSIONS

We propose two novel biomarkers that support the role of meat consumption with an increased risk of colorectal cancer.

IMPACT

The reported GxE interactions may explain the increased risk of colorectal cancer in certain population subgroups.

Unveiling challenges in Mendelian randomization for gene-environment interaction

Gene-environment (GxE) interactions play a crucial role in understanding the complex etiology of various traits, but assessing them using observational data can be challenging due to unmeasured confounders for lifestyle and environmental risk factors. Mendelian randomization (MR) has emerged as a valuable method for assessing causal relationships based on observational data. This approach utilizes genetic variants as instrumental variables (IVs) with the aim of providing a valid statistical test and estimation of causal effects in the presence of unmeasured confounders. MR has gained substantial popularity in recent years largely due to the success of genome-wide association studies. Many methods have been developed for MR; however, limited work has been done on evaluating GxE interaction. In this paper, we focus on two primary IV approaches: the two-stage predictor substitution and the two-stage residual inclusion, and extend them to accommodate GxE interaction under both the linear and logistic regression models for continuous and binary outcomes, respectively. Comprehensive simulation study and analytical derivations reveal that resolving the linear regression model is relatively straightforward. In contrast, the logistic regression model presents a considerably more intricate challenge, which demands additional effort.

Genome-wide association studies and Mendelian randomization analyses provide insights into the causes of early-onset colorectal cancer

BACKGROUND

The incidence of early-onset colorectal cancer (EOCRC; diagnosed <50 years of age) is rising globally; however, the causes underlying this trend are largely unknown. CRC has strong genetic and environmental determinants, yet common genetic variants and causal modifiable risk factors underlying EOCRC are unknown. We conducted the first EOCRC-specific genome-wide association study (GWAS) and Mendelian randomization (MR) analyses to explore germline genetic and causal modifiable risk factors associated with EOCRC.

PATIENTS AND METHODS

We conducted a GWAS meta-analysis of 6176 EOCRC cases and 65 829 controls from the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO), the Colorectal Transdisciplinary Study (CORECT), the Colon Cancer Family Registry (CCFR), and the UK Biobank. We then used the EOCRC GWAS to investigate 28 modifiable risk factors using two-sample MR.

RESULTS

We found two novel risk loci for EOCRC at 1p34.1 and 4p15.33, which were not previously associated with CRC risk. We identified a deleterious coding variant (rs36053993, G396D) at polyposis-associated DNA repair gene MUTYH (odds ratio 1.80, 95% confidence interval 1.47-2.22) but show that most of the common genetic susceptibility was from noncoding signals enriched in epigenetic markers present in gastrointestinal tract cells. We identified new EOCRC-susceptibility genes, and in addition to pathways such as transforming growth factor (TGF) β, suppressor of Mothers Against Decapentaplegic (SMAD), bone morphogenetic protein (BMP) and phosphatidylinositol kinase (PI3K) signaling, our study highlights a role for insulin signaling and immune/infection-related pathways in EOCRC. In our MR analyses, we found novel evidence of probable causal associations for higher levels of body size and metabolic factors-such as body fat percentage, waist circumference, waist-to-hip ratio, basal metabolic rate, and fasting insulin-higher alcohol drinking, and lower education attainment with increased EOCRC risk.

CONCLUSIONS

Our novel findings indicate inherited susceptibility to EOCRC and suggest modifiable lifestyle and metabolic targets that could also be used to risk-stratify individuals for personalized screening strategies or other interventions.

Five genes identified as prognostic markers for colorectal cancer through the integration of genome-wide association study and expression quantitative trait loci data

Background: Colorectal cancer (CRC) is a prominent form of cancer globally, ranking second in terms of prevalence and serving as a leading cause of cancer-related deaths, but the underlying biological interpretation remains largely unknown. Methods: We used the summary data-based Mendelian randomization method to integrate CRC genome-wide association studies (ncase = 7062; ncontrol = 195,745) and expression quantitative trait loci summary data in peripheral whole blood (Consortium for Architecture of Gene Expression: n = 2765; Genotype-Tissue Expression [v8]: n = 755) and colon tissue (colon-transverse: n = 406; colon-sigmoid: n = 373) and identified related genes. Results: Genes ABTB1, CYP21A2, NLRP1, PHKG1 and PIP5K1C have emerged as significant prognostic markers for CRC patient survival. Functional analysis revealed their involvement in cancer cell migration and invasion mechanisms, providing valuable insights for the development of future anti-CRC drugs. Conclusion: We successfully identified five CRC risk genes, providing new insights and research directions for the effective mechanisms of CRC.

Immunotherapy for colorectal cancer: insight from inherited genetics

Immunotherapy shows efficacy for multiple cancer types and potential for expanded use. However, current immune checkpoint inhibitors (ICIs) are ineffective against microsatellite-stable colorectal cancer (CRC), which is more commonly diagnosed. Immunotherapy strategies for non-responsive CRC, including new targets and new combination therapies, are being tested to address this need. Importantly, a subset of inherited germline genetic variants associated with CRC risk are predicted to regulate genes with immune functions, including genes related to existing ICIs, as well as new potential targets in the major histocompatibility complex (MHC) region and immunoregulatory cytokines. We review discoveries in the inherited genetics of CRC related to the immune system and draw connections with ongoing developments and emerging immunotherapy targets.

Healthy dietary patterns, genetic risk, and gastrointestinal cancer incident risk: a large-scale prospective cohort study

BACKGROUND

Dietary patterns have been associated with several cancers, especially gastrointestinal cancer (GIC). However, whether a healthy dietary pattern could modify the risk of GIC among people with different genetic backgrounds is not clear.

OBJECTIVE

The objective of the study was to investigate how dietary patterns and genetic susceptibility contribute to the risk of GIC independently and jointly.

METHODS

This large-scale prospective cohort study included 105,463 participants in UK Biobank who were aged 40-72 y and cancer-free at baseline. Dietary intake (Oxford WebQ) was used to calculate dietary pattern scores including dietary approach to stop hypertension (DASH) score and healthful plant-based diet index (hPDI). Genetic risk was quantified by a polygenic risk score (PRS) comprising 129 known GIC-associated loci. Cox proportional hazards regression was performed to estimate the associations of dietary patterns and PRS with GIC incidence after adjusting for potential confounders.

RESULTS

Over a median follow-up of 11.70 y, 1,661 participants were diagnosed with GIC. DASH and hPDI were associated with 20% and 36% reductions, respectively, in GIC risk. Low PRS was associated with a 30 % decrease in GIC risk (HR: 0.70; 95% CI: 0.62, 0.79). Participants with healthy dietary scores at high-genetic risk had a lower GIC risk with HR of 0.77 (95% CI: 0.60, 0.98) for DASH and 0.66 (95% CI: 0.52, 0.84) for hPDI than those with unhealthy dietary score. Participants with both high-dietary score and low-genetic risk showed the lowest risk of GIC, with HR of 0.58 (95% CI: 0.45, 0.75) for DASH and 0.45 (95% CI: 0.34, 0.58) for hPDI.

CONCLUSIONS

Adherence to DASH and hPDI were associated with a lower risk of some gastrointestinal cancers, and these 2 dietary patterns may partly compensate for genetic predispositions to cancer. Our results advance the development of precision medicine strategies that consider both dietary patterns and genetics to improve gastrointestinal health.

Causal relationships between risk of venous thromboembolism and 18 cancers: a bidirectional Mendelian randomization analysis

BACKGROUND

People with cancer experience high rates of venous thromboembolism (VTE). Risk of subsequent cancer is also increased in people experiencing their first VTE. The causal mechanisms underlying this association are not completely understood, and it is unknown whether VTE is itself a risk factor for cancer.

METHODS

We used data from large genome-wide association study meta-analyses to perform bidirectional Mendelian randomization analyses to estimate causal associations between genetic liability to VTE and risk of 18 different cancers.

RESULTS

We found no conclusive evidence that genetic liability to VTE was causally associated with an increased incidence of cancer, or vice versa. We observed an association between liability to VTE and pancreatic cancer risk [odds ratio for pancreatic cancer: 1.23 (95% confidence interval: 1.08-1.40) per log-odds increase in VTE risk, P = 0.002]. However, sensitivity analyses revealed this association was predominantly driven by a variant proxying non-O blood group, with inadequate evidence to suggest a causal relationship.

CONCLUSIONS

These findings do not support the hypothesis that genetic liability to VTE is a cause of cancer. Existing observational epidemiological associations between VTE and cancer are therefore more likely to be driven by pathophysiological changes which occur in the setting of active cancer and anti-cancer treatments. Further work is required to explore and synthesize evidence for these mechanisms.

Impact of weight loss on cancer-related proteins in serum: results from a cluster randomised controlled trial of individuals with type 2 diabetes

BACKGROUND

Type 2 diabetes is associated with higher risk of several cancer types. However, the biological intermediates driving this relationship are not fully understood. As novel interventions for treating and managing type 2 diabetes become increasingly available, whether they also disrupt the pathways leading to increased cancer risk is currently unknown. We investigated the effect of a type 2 diabetes intervention, in the form of intentional weight loss, on circulating proteins associated with cancer risk to gain insight into potential mechanisms linking type 2 diabetes and adiposity with cancer development.

METHODS

Fasting serum samples from participants with diabetes enrolled in the Diabetes Remission Clinical Trial (DiRECT) receiving the Counterweight-Plus weight-loss programme (intervention, N = 117, mean weight-loss 10 kg, 46% diabetes remission) or best-practice care by guidelines (control, N = 143, mean weight-loss 1 kg, 4% diabetes remission) were subject to proteomic analysis using the Olink Oncology-II platform (48% of participants were female; 52% male). To identify proteins which may be altered by the weight-loss intervention, the difference in protein levels between groups at baseline and 1 year was examined using linear regression. Mendelian randomization (MR) was performed to extend these results to evaluate cancer risk and elucidate possible biological mechanisms linking type 2 diabetes and cancer development. MR analyses were conducted using independent datasets, including large cancer meta-analyses, UK Biobank, and FinnGen, to estimate potential causal relationships between proteins modified during intentional weight loss and the risk of colorectal, breast, endometrial, gallbladder, liver, and pancreatic cancers.

FINDINGS

Nine proteins were modified by the intervention: glycoprotein Nmb; furin; Wnt inhibitory factor 1; toll-like receptor 3; pancreatic prohormone; erb-b2 receptor tyrosine kinase 2; hepatocyte growth factor; endothelial cell specific molecule 1 and Ret proto-oncogene (Holm corrected P-value <0.05). Mendelian randomization analyses indicated a causal relationship between predicted circulating furin and glycoprotein Nmb on breast cancer risk (odds ratio (OR) = 0.81, 95% confidence interval (CI) = 0.67-0.99, P-value = 0.03; and OR = 0.88, 95% CI = 0.78-0.99, P-value = 0.04 respectively), though these results were not supported in sensitivity analyses examining violations of MR assumptions.

INTERPRETATION

Intentional weight loss among individuals with recently diagnosed diabetes may modify levels of cancer-related proteins in serum. Further evaluation of the proteins identified in this analysis could reveal molecular pathways that mediate the effect of adiposity and type 2 diabetes on cancer risk.

FUNDING

The main sources of funding for this work were Diabetes UK, Cancer Research UK, World Cancer Research Fund, and Wellcome.

A large-scale microRNA transcriptome-wide association study identifies two susceptibility microRNAs, miR-1307-5p and miR-192-3p, for colorectal cancer risk

Transcriptome-wide association studies (TWAS) have identified many putative susceptibility genes for colorectal cancer (CRC) risk. However, susceptibility miRNAs, critical dysregulators of gene expression, remain unexplored. We genotyped DNA samples from 313 CRC East Asian patients and performed small RNA sequencing in their normal colon tissues distant from tumors to build genetic models for predicting miRNA expression. We applied these models and data from genome-wide association studies (GWAS) including 23 942 cases and 217 267 controls of East Asian ancestry to investigate associations of predicted miRNA expression with CRC risk. Perturbation experiments separately by promoting and inhibiting miRNAs expressions and further in vitro assays in both SW480 and HCT116 cells were conducted. At a Bonferroni-corrected threshold of P < 4.5 × 10-4, we identified two putative susceptibility miRNAs, miR-1307-5p and miR-192-3p, located in regions more than 500 kb away from any GWAS-identified risk variants in CRC. We observed that a high predicted expression of miR-1307-5p was associated with increased CRC risk, while a low predicted expression of miR-192-3p was associated with increased CRC risk. Our experimental results further provide strong evidence of their susceptible roles by showing that miR-1307-5p and miR-192-3p play a regulatory role, respectively, in promoting and inhibiting CRC cell proliferation, migration, and invasion, which was consistently observed in both SW480 and HCT116 cells. Our study provides additional insights into the biological mechanisms underlying CRC development.

Association between circulating inflammatory markers and adult cancer risk: a Mendelian randomization analysis

BACKGROUND

Tumour-promoting inflammation is a "hallmark" of cancer and conventional epidemiological studies have reported links between various inflammatory markers and cancer risk. The causal nature of these relationships and, thus, the suitability of these markers as intervention targets for cancer prevention is unclear.

METHODS

We meta-analysed 6 genome-wide association studies of circulating inflammatory markers comprising 59,969 participants of European ancestry. We then used combined cis-Mendelian randomization and colocalisation analysis to evaluate the causal role of 66 circulating inflammatory markers in risk of 30 adult cancers in 338,294 cancer cases and up to 1,238,345 controls. Genetic instruments for inflammatory markers were constructed using genome-wide significant (P < 5.0 × 10-8) cis-acting SNPs (i.e., in or ±250 kb from the gene encoding the relevant protein) in weak linkage disequilibrium (LD, r2 < 0.10). Effect estimates were generated using inverse-variance weighted random-effects models and standard errors were inflated to account for weak LD between variants with reference to the 1000 Genomes Phase 3 CEU panel. A false discovery rate (FDR)-corrected P-value ("q-value") <0.05 was used as a threshold to define "strong evidence" to support associations and 0.05 ≤ q-value < 0.20 to define "suggestive evidence". A colocalisation posterior probability (PPH4) >70% was employed to indicate support for shared causal variants across inflammatory markers and cancer outcomes. Findings were replicated in the FinnGen study and then pooled using meta-analysis.

FINDINGS

We found strong evidence to support an association of genetically-proxied circulating pro-adrenomedullin concentrations with increased breast cancer risk (OR: 1.19, 95% CI: 1.10-1.29, q-value = 0.033, PPH4 = 84.3%) and suggestive evidence to support associations of interleukin-23 receptor concentrations with increased pancreatic cancer risk (OR: 1.42, 95% CI: 1.20-1.69, q-value = 0.055, PPH4 = 73.9%), prothrombin concentrations with decreased basal cell carcinoma risk (OR: 0.66, 95% CI: 0.53-0.81, q-value = 0.067, PPH4 = 81.8%), and interleukin-1 receptor-like 1 concentrations with decreased triple-negative breast cancer risk (OR: 0.92, 95% CI: 0.88-0.97, q-value = 0.15, PPH4 = 85.6%). These findings were replicated in pooled analyses with the FinnGen study. Though suggestive evidence was found to support an association of macrophage migration inhibitory factor concentrations with increased bladder cancer risk (OR: 2.46, 95% CI: 1.48-4.10, q-value = 0.072, PPH4 = 76.1%), this finding was not replicated when pooled with the FinnGen study. For 22 of 30 cancer outcomes examined, there was little evidence (q-value ≥0.20) that any of the 66 circulating inflammatory markers examined were associated with cancer risk.

INTERPRETATION

Our comprehensive joint Mendelian randomization and colocalisation analysis of the role of circulating inflammatory markers in cancer risk identified potential roles for 4 circulating inflammatory markers in risk of 4 site-specific cancers. Contrary to reports from some prior conventional epidemiological studies, we found little evidence of association of circulating inflammatory markers with the majority of site-specific cancers evaluated.

FUNDING

Cancer Research UK (C68933/A28534, C18281/A29019, PPRCPJT∖100005), World Cancer Research Fund (IIG_FULL_2020_022), National Institute for Health Research (NIHR202411, BRC-1215-20011), Medical Research Council (MC_UU_00011/1, MC_UU_00011/3, MC_UU_00011/6, and MC_UU_00011/4), Academy of Finland Project 326291, European Union's Horizon 2020 grant agreement no. 848158 (EarlyCause), French National Cancer Institute (INCa SHSESP20, 2020-076), Versus Arthritis (21173, 21754, 21755), National Institutes of Health (U19 CA203654), National Cancer Institute (U19CA203654).

Benzo[a]pyrene exposure affects colorectal cancer susceptibility by regulating ERβ-mediated LINC02977 transcription

Environmental pollutants known as polycyclic aromatic hydrocarbons (PAHs) are produced through the incomplete combustion of organic material. While PAHs have been investigated as genotoxicants, they can also operate through nongenotoxic pathways in estrogen-dependent malignancies, such as breast, cervical and ovarian cancer. However, whether PAHs induce colorectal cancer (CRC) risk through estrogenic effects is still illusive. Here, we systematically investigated the abnormal expression and activation of estrogen receptor beta (ERβ) regulated by PAHs in CRC as well as the underlying mechanisms of ERβ-mediated CRC risk. Based on the 300 plasma samples from CRC patients and healthy controls detected by GC-MS/MS, we found that the plasma concentrations of benzo[a]pyrene (BaP) were significantly higher in CRC cases than in healthy controls, with significant estrogenic effects. Moreover, histone deacetylase 2 (HDAC2)-induced deacetylation of the promoter decreases ERβ expression, which is associated with poor overall survival and advanced tumor stage. The study also revealed that BaP and estradiol (E2) had different carcinogenic effects, with BaP promoting cell proliferation and inhibiting apoptosis, while E2 had the opposite effects. Additionally, this study mapped ERβ genomic binding regions by performing ChIP-seq and ATAC-seq and identified genetic variants of rs1411680 and its high linkage disequilibrium SNP rs6477937, which were significantly associated with CRC risk through meta-analysis of two independent Chinese population genome-wide association studies comprising 2,248 cases and 3,173 controls and then validation in a large-scale European population. By integrating data from functional genomics, we validated the regulatory effect of rs6477937 as an ERβ binding-disrupting SNP that mediated allele-specific expression of LINC02977 in a long-range chromosomal interaction manner, which was found to be highly expressed in CRC tissues. Overall, this study suggests that the different active effects on ERβ by PAHs and endogenous E2 may play a crucial role in the development and progression of CRC and highlights the potential of targeting ERβ and its downstream targets for CRC prevention and treatment.

Association between serum urea concentrations and the risk of colorectal cancer, particularly in individuals with type 2 diabetes: A cohort study

Dysregulation of the urea cycle (UC) has been detected in colorectal cancer (CRC). However, the impact of the UC's end product, urea, on CRC development remains unclear. We investigated the association between serum urea and CRC risk based on the data of 348 872 participants cancer-free at recruitment from the UK Biobank. Multivariable Cox proportional hazards models were fitted to conduct risk estimates. Stratification analyses based on sex, diet pattern, metabolic factors (including body mass index [BMI], the estimated glomerular filtration rate [eGFR] and type 2 diabetes [T2D]) and genetic profiles (the polygenic risk score [PRS] of CRC) were conducted to find potential modifiers. During an average of 9.0 years of follow-up, we identified 3408 (1.0%) CRC incident cases. Serum urea showed a nonlinear relationship with CRC risk (P-nonlinear: .035). Lower serum urea levels were associated with a higher CRC risk, with a fully-adjusted hazard ratio (HR) of 1.26 (95% confidence interval [CI]: 1.13-1.41) in the first quartile (Q1) of urea, compared to the Q4. This association was largely consistent across subgroups of sex, protein diet, BMI, eGFR and CRC-PRSs (P-interaction >.05); however, it was stronger in the T2D, with an interaction between urea and T2D on both additive (synergy index: 3.32, [95% CI: 1.24-8.88]) and multiplicative scales (P-interaction: .019). Lower serum urea concentrations were associated with an increased risk of CRC, with a more pronounced effect observed in individuals with T2D. Maintaining stable levels of serum urea has important implications for CRC prevention, particularly in individuals with T2D.

Bidirectional Mendelian randomization analysis investigating the genetic association between primary breast cancer and colorectal cancer

Purpose

With the advancement in early diagnosis and treatment, the prognosis for individuals diagnosed with breast cancer (BC) has improved significantly. The prognosis of primary breast cancer (PBC) survivors can be significantly influenced by the occurrence of colorectal cancer (CRC) as a secondary primary cancer (SPC). The objective of this study is to explore the possible genetic association between PBC and CRC, aiming to lay a groundwork for the development of preventive strategies against SPC-CRC following BC surgery.

Methods

We employed a bidirectional two-sample Mendelian randomization (MR) approach to thoroughly examine genetic instrumental variables (IVs) derived from genome-wide association studies (GWAS) conducted on PBC and CRC. And applied inverse variance weighted (IVW) and multiple other MR methods (weighted median, simple median, MR-PRESSO and MR-RAPS) to evaluate the association between the two cancers (PBC and CRC) at genetic level. Furthermore, the robustness of the findings was further confirmed through the utilization of the genetic risk score (GRS) method in a secondary analysis.

Results

Forward MR analysis, a total of 179 BC genetic IVs, 25 estrogen receptor-negative (ER-) genetic IVs and 135 ER-positive (ER+) genetic IVs were screened. Reverse MR analysis, 179 genetic IVs of CRC, 25 genetic IVs of colon cancer, 135 genetic IVs of rectal cancer, 25 genetic IVs of left colon cancer and 135 genetic IVs of right colon cancer were screened. IVW and other MR methods found no significant genetic association between PBC and CRC (P > 0.05). Subgroup analysis also showed that ER- BC and ER+ BC were not correlated with the occurrence of CRC (P > 0.05). The findings of the secondary analysis using GRS were consistent with those obtained from the primary analysis, thereby confirming the robustness and reliability of this study.

Conclusions

Our findings do not provide any evidence supporting the association between PBC and CRC at the genetic level. Further large-scale prospective studies are warranted to replicate our findings.

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.

Observational and genetic associations between cardiorespiratory fitness and cancer: a UK Biobank and international consortia study

BACKGROUND

The association of fitness with cancer risk is not clear.

METHODS

We used Cox proportional hazards models to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for risk of lung, colorectal, endometrial, breast, and prostate cancer in a subset of UK Biobank participants who completed a submaximal fitness test in 2009-12 (N = 72,572). We also investigated relationships using two-sample Mendelian randomisation (MR), odds ratios (ORs) were estimated using the inverse-variance weighted method.

RESULTS

After a median of 11 years of follow-up, 4290 cancers of interest were diagnosed. A 3.5 ml O2⋅min-1⋅kg-1 total-body mass increase in fitness (equivalent to 1 metabolic equivalent of task (MET), approximately 0.5 standard deviation (SD)) was associated with lower risks of endometrial (HR = 0.81, 95% CI: 0.73-0.89), colorectal (0.94, 0.90-0.99), and breast cancer (0.96, 0.92-0.99). In MR analyses, a 0.5 SD increase in genetically predicted O2⋅min-1⋅kg-1 fat-free mass was associated with a lower risk of breast cancer (OR = 0.92, 95% CI: 0.86-0.98). After adjusting for adiposity, both the observational and genetic associations were attenuated.

DISCUSSION

Higher fitness levels may reduce risks of endometrial, colorectal, and breast cancer, though relationships with adiposity are complex and may mediate these relationships. Increasing fitness, including via changes in body composition, may be an effective strategy for cancer prevention.

Genetic and microbial determinants of azoxymethane-induced colorectal tumor susceptibility in Collaborative Cross mice and their implication in human cancer

The insights into interactions between host genetics and gut microbiome (GM) in colorectal tumor susceptibility (CTS) remains lacking. We used Collaborative Cross mouse population model to identify genetic and microbial determinants of Azoxymethane-induced CTS. We identified 4417 CTS-associated single nucleotide polymorphisms (SNPs) containing 334 genes that were transcriptionally altered in human colorectal cancers (CRCs) and consistently clustered independent human CRC cohorts into two subgroups with different prognosis. We discovered a set of genera in early-life associated with CTS and defined a 16-genus signature that accurately predicted CTS, the majority of which were correlated with human CRCs. We identified 547 SNPs associated with abundances of these genera. Mediation analysis revealed GM as mediators partially exerting the effect of SNP UNC3869242 within Duox2 on CTS. Intestine cell-specific depletion of Duox2 altered GM composition and contribution of Duox2 depletion to CTS was significantly influenced by GM. Our findings provide potential novel targets for personalized CRC prevention and treatment.

Prioritization of risk genes in colorectal cancer by integrative analysis of multi-omics data and gene networks

Genome-wide association studies (GWASs) have identified over 140 colorectal cancer (CRC)-associated loci; however, target genes at the majority of loci and underlying molecular mechanisms are poorly understood. Here, we utilized a Bayesian approach, integrative risk gene selector (iRIGS), to prioritize risk genes at CRC GWAS loci by integrating multi-omics data. As a result, a total of 105 high-confidence risk genes (HRGs) were identified, which exhibited strong gene dependencies for CRC and enrichment in the biological processes implicated in CRC. Among the 105 HRGs, CEBPB, located at the 20q13.13 locus, acted as a transcription factor playing critical roles in cancer. Our subsequent assays indicated the tumor promoter function of CEBPB that facilitated CRC cell proliferation by regulating multiple oncogenic pathways such as MAPK, PI3K-Akt, and Ras signaling. Next, by integrating a fine-mapping analysis and three independent case-control studies in Chinese populations consisting of 8,039 cases and 12,775 controls, we elucidated that rs1810503, a putative functional variant regulating CEBPB, was associated with CRC risk (OR=0.90, 95%CI=0.86-0.93, P=1.07×10-7). The association between rs1810503 and CRC risk was further validated in three additional multi-ancestry populations consisting of 24,254 cases and 58,741 controls. Mechanistically, the rs1810503 A to T allele change weakened the enhancer activity in an allele-specific manner to decrease CEBPB expression via long-range promoter-enhancer interactions, mediated by the transcription factor, REST, and thus decreased CRC risk. In summary, our study provides a genetic resource and a generalizable strategy for CRC etiology investigation, and highlights the biological implications of CEBPB in CRC tumorigenesis, shedding new light on the etiology of CRC.

Leisure time television watching, computer use and risks of breast, colorectal and prostate cancer: A Mendelian randomisation analysis

BACKGROUND

Sedentary behaviours have been associated with increased risks of some common cancers in epidemiological studies; however, it is unclear if these associations are causal.

METHODS

We used univariable and multivariable two-sample Mendelian randomisation (MR) to examine potential causal relationships between sedentary behaviours and risks of breast, colorectal and prostate cancer. Genetic variants associated with self-reported leisure television watching and computer use were identified from a recent genome-wide association study (GWAS). Data related to cancer risk were obtained from cancer GWAS consortia. A series of sensitivity analyses were applied to examine the robustness of the results to the presence of confounding.

RESULTS

A 1-standard deviation (SD: 1.5 h/day) increment in hours of television watching increased risk of breast cancer (OR per 1-SD: 1.15, 95% confidence interval [CI]: 1.05-1.26) and colorectal cancer (OR per 1-SD: 1.32, 95% CI: 1.16-1.49) while there was little evidence of an association for prostate cancer risk (OR per 1-SD: 0.94, 95% CI: 0.84-1.06). After adjusting for years of education, the effect estimates for television watching were attenuated (breast cancer, OR per 1-SD: 1.08, 95% CI: 0.92-1.27; colorectal cancer, OR per 1-SD: 1.08, 95% CI: 0.90-1.31). Post hoc analyses showed that years of education might have a possible confounding and mediating role in the association between television watching with breast and colorectal cancer. Consistent results were observed for each cancer site according to sex (colorectal cancer), anatomical subsites and cancer subtypes. There was little evidence of associations between genetically predicted computer use and cancer risk.

CONCLUSIONS

Our univariable analysis identified some positive associations between hours of television watching and risks of breast and colorectal cancer. However, further adjustment for additional lifestyle factors especially years of education attenuated these results. Future studies using objective measures of exposure can provide new insights into the possible role of sedentary behaviour in cancer development.

Identifying metabolic features of colorectal cancer liability using Mendelian randomization

Background

Recognizing the early signs of cancer risk is vital for informing prevention, early detection, and survival.

Methods

To investigate whether changes in circulating metabolites characterize the early stages of colorectal cancer (CRC) development, we examined the associations between a genetic risk score (GRS) associated with CRC liability (72 single-nucleotide polymorphisms) and 231 circulating metabolites measured by nuclear magnetic resonance spectroscopy in the Avon Longitudinal Study of Parents and Children (N = 6221). Linear regression models were applied to examine the associations between genetic liability to CRC and circulating metabolites measured in the same individuals at age 8 y, 16 y, 18 y, and 25 y.

Results

The GRS for CRC was associated with up to 28% of the circulating metabolites at FDR-P < 0.05 across all time points, particularly with higher fatty acids and very-low- and low-density lipoprotein subclass lipids. Two-sample reverse Mendelian randomization (MR) analyses investigating CRC liability (52,775 cases, 45,940 controls) and metabolites measured in a random subset of UK Biobank participants (N = 118,466, median age 58 y) revealed broadly consistent effect estimates with the GRS analysis. In conventional (forward) MR analyses, genetically predicted polyunsaturated fatty acid concentrations were most strongly associated with higher CRC risk.

Conclusions

These analyses suggest that higher genetic liability to CRC can cause early alterations in systemic metabolism and suggest that fatty acids may play an important role in CRC development.

Funding

This work was supported by the Elizabeth Blackwell Institute for Health Research, University of Bristol, the Wellcome Trust, the Medical Research Council, Diabetes UK, the University of Bristol NIHR Biomedical Research Centre, and Cancer Research UK. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. This work used the computational facilities of the Advanced Computing Research Centre, University of Bristol - http://www.bristol.ac.uk/acrc/.

[Associations of circulating leptin levels with colorectal adenoma and colorectal cancer: a case-control and Mendelian randomization study]

OBJECTIVE

To explore the causal association between circulating leptin levels and the risk of colorectal adenoma and colorectal cancer.

METHODS

We collected demographic and clinical data and serum samples from 497 patients with colorectal adenoma, 955 patients with colorectal cancer, and 911 healthy individuals from the First Affiliated Hospital of Zhejiang University School of Medicine, Zhejiang Cancer Hospital, Zhuji People's Hospital, and Lin'an District First People's Hospital. Instrumental variables of leptin were selected and genotyping tests were performed. A logistic regression model and stratified analysis were used to evaluate the association of serum leptin levels with colorectal adenoma, colorectal cancer, and the progression of colorectal adenoma to colorectal cancer. Genetic risk score (GRS) and single nucleotide polymorphisms (SNPs) were further used as instrumental variables in one-sample and two-sample Mendelian randomization analyses leveraging two-stage least squares and inverse-variance weighted methods to estimate the causal association of leptin levels with the risk of colorectal adenoma, colorectal cancer, and progression of colorectal adenoma to colorectal cancer.

RESULTS

High levels of leptin, compared with its lowest quartile, were positively correlated with colorectal adenoma (P=0.005) and negatively with colorectal cancer (P < 0.001) and the risk of progression of colorectal adenoma to colorectal cancer (P < 0.001). Mendelian randomization analysis showed that GRS of leptin, either weighted or not, was not significantly correlated with the risk of colorectal adenoma, colorectal cancer, or the progression of colorectal adenoma to colorectal cancer, nor did the two-sample Mendelian randomization study support an association between leptin and the risk of colorectal cancer (P>0.05).

CONCLUSION

Although the case-control study suggests probable correlations of leptin with the risk of colorectal adenoma, colorectal cancer, and colorectal adenoma progression to colorectal cancer, Mendelian randomization studies did not support a causal association of leptin with the risks of colorectal adenoma, colorectal cancer, or colorectal adenoma progression to colorectal cancer.

Causal genetic regulation of DNA replication on immune microenvironment in colorectal tumorigenesis: Evidenced by an integrated approach of trans-omics and GWAS

The interplay between DNA replication stress and immune microenvironment alterations is known to play a crucial role in colorectal tumorigenesis, but a comprehensive understanding of their association with and relevant biomarkers involved in colorectal tumorigenesis is lacking. To address this gap, we conducted a study aiming to investigate this association and identify relevant biomarkers. We analyzed transcriptomic and proteomic profiles of 904 colorectal tumor tissues and 342 normal tissues to examine pathway enrichment, biological activity, and the immune microenvironment. Additionally, we evaluated genetic effects of single variants and genes on colorectal cancer susceptibility using data from genome-wide association studies (GWASs) involving both East Asian (7062 cases and 195745 controls) and European (24476 cases and 23073 controls) populations. We employed mediation analysis to infer the causal pathway, and applied multiplex immunofluorescence to visualize colocalized biomarkers in colorectal tumors and immune cells. Our findings revealed that both DNA replication activity and the flap structure-specific endonuclease 1 ( FEN1) gene were significantly enriched in colorectal tumor tissues, compared with normal tissues. Moreover, a genetic variant rs4246215 G>T in FEN1 was associated with a decreased risk of colorectal cancer (odds ratio = 0.94, 95% confidence interval: 0.90-0.97, P meta = 4.70 × 10 -9). Importantly, we identified basophils and eosinophils that both exhibited a significantly decreased infiltration in colorectal tumors, and were regulated by rs4246215 through causal pathways involving both FEN1 and DNA replication. In conclusion, this trans-omics incorporating GWAS data provides insights into a plausible pathway connecting DNA replication and immunity, expanding biological knowledge of colorectal tumorigenesis and therapeutic targets.

A Rare Variant in MDH2 (rs111879470) Is Associated with Predisposition to Recurrent Breast Cancer in an Extended High-Risk Pedigree

A significant fraction of breast cancer recurs, with lethal outcome, but specific genetic variants responsible have yet to be identified. Five cousin pairs with recurrent breast cancer from pedigrees with a statistical excess of recurrent breast cancer were sequenced to identify rare, shared candidate predisposition variants. The candidates were tested for association with breast cancer risk with UKBiobank data. Additional breast cancer cases were assayed for a subset of candidate variants to test for co-segregation. Three-dimensional protein structure prediction methods were used to investigate how the mutation under consideration is predicted to change structural and electrostatic properties in the mutated protein. One hundred and eighty-one rare candidate predisposition variants were shared in at least one cousin pair from a high-risk pedigree. A rare variant in MDH2 was found to segregate with breast-cancer-affected relatives in one extended pedigree. MDH2 is an estrogen-stimulated gene encoding the protein malate dehydrogenase, which catalyzes the reversible oxidation of malate to oxaloacetate. The molecular simulation results strongly suggest that the mutation changes the NAD+ binding pocket electrostatics of MDH2. This small sequencing study, using a powerful approach based on recurrent breast cancer cases from high-risk pedigrees, identified a set of strong candidate variants for inherited predisposition for breast cancer recurrence, including MDH2, which should be pursued in other resources.

Polygenic risk alters the penetrance of monogenic kidney disease

Chronic kidney disease (CKD) is determined by an interplay of monogenic, polygenic, and environmental risks. Autosomal dominant polycystic kidney disease (ADPKD) and COL4A-associated nephropathy (COL4A-AN) represent the most common forms of monogenic kidney diseases. These disorders have incomplete penetrance and variable expressivity, and we hypothesize that polygenic factors explain some of this variability. By combining SNP array, exome/genome sequence, and electronic health record data from the UK Biobank and All-of-Us cohorts, we demonstrate that the genome-wide polygenic score (GPS) significantly predicts CKD among ADPKD monogenic variant carriers. Compared to the middle tertile of the GPS for noncarriers, ADPKD variant carriers in the top tertile have a 54-fold increased risk of CKD, while ADPKD variant carriers in the bottom tertile have only a 3-fold increased risk of CKD. Similarly, the GPS significantly predicts CKD in COL4A-AN carriers. The carriers in the top tertile of the GPS have a 2.5-fold higher risk of CKD, while the risk for carriers in the bottom tertile is not different from the average population risk. These results suggest that accounting for polygenic risk improves risk stratification in monogenic kidney disease.

Deleterious and ethnic-related BRCA1/2 mutations in tissue and blood of Egyptian colorectal cancer patients and its correlation with human papillomavirus

This study aimed to identify BRCA1/2 mutational patterns in the tissue and blood of Egyptian colorectal cancer (CRC) patients and to study the possible correlation of this mutational pattern with Human papillomavirus (HPV) infection. Eighty-two colonoscopic biopsies and forty-six blood samples were collected from Egyptian CRC patients, as well as blood samples of age and sex-matched healthy controls (n = 43) were enrolled. The libraries were performed using Qiaseq Human BRCA1 and BRCA2 targeted DNA panel and sequenced via Ion proton sequencer. Also, the CRC tissues were subjected to conventional PCR targeting the HPV Late 1 (L1) region. Our analysis revealed that the BRCA-DNA damage pathway had been altered in more than 65% of the CRC patients. Comparing tissue and blood samples from CRC patients, 25 somatic mutations were found exclusively in tissue, while 41 germline mutations were found exclusively in blood. Additionally, we identified 23 shared BRCA1/2 pathogenic (PVs) mutations in both blood and tissue samples, with a significantly higher frequency in blood samples compared to tissue samples. The most affected exon in BRCA1 was exon 10, while the most affected exons in BRCA2 were 11, 14, 18, 24, and 27 exons. Notably, we revealed an ethnic-related cluster of polymorphism variants in our population closely related to South Asian and African ethnicities. Novel PVs were identified and submitted to the ClinVar database. HPV was found in 23.8% of the CRC tissues, and 54% of HPV-positive cases had somatic BRCA1/2 PVs. The results of this research point to a possible connection between infection with HPV and BRCA1/2 mutations in the occurrence of colorectal cancer in the Egyptian population, which has a mixed ethnic background. Our data also indicate that liquid biopsy (blood samples) may be more representative than tissue samples for detecting BRCA1/2 mutations. These findings may have implications for cancer screening and the development of personalized, targeted therapies, such as PARP inhibitors, which can effectively target BRCA1/2 mutations.

Diet-Wide Association, Genetic Susceptibility and Colorectal Cancer Risk: A Prospective Cohort Study

BACKGROUND

Both genetic and dietary factors play significant roles in the etiology of colorectal cancer (CRC). To evaluate the relationship between certain food exposures and the risk of CRC, we carried out a large-scale association analysis in the UK Biobank.

METHODS

The associations of 139 foods and nutrients' intake with CRC risk were assessed among 118,210 participants. A polygenic risk score (PRS) of CRC was created to explore any interaction between dietary factors and genetic susceptibility in CRC risk. The hazard ratio (HR) and 95% confidence interval (CI) of CRC risk linked to dietary variables and PRS were estimated using Cox regression models. Multiple comparisons were corrected using the error discovery rate (FDR).

RESULTS

During a mean follow-up of 12.8 years, 1466 incidents of CRC were identified. In the UK Biobank, alcohol and white bread were associated with increased CRC risk, and their HRs were 1.08 (95% CI: 1.03-1.14; FDRP = 0.028) and 1.10 (95% CI: 1.05-1.16; FDRP = 0.003), whereas dietary fiber, calcium, magnesium, phosphorus, and manganese intakes were inversely associated. We found no evidence of any PRS-nutrient interaction relationship in relation to CRC risk.

CONCLUSIONS

Our results show that higher intakes of alcohol and white bread are associated with increased CRC risk, whilst dietary fiber, calcium, magnesium, phosphorus, and manganese are inversely associated.

Clustering of colon, lung, and other cancer susceptibility genes with protein tyrosine phosphatases and protein kinases in multiple short genomic regions

Interactions of large gene families are poorly understood. We found that human, mouse, and rat colon and lung cancer susceptibility genes, presently considered as separate gene families, were frequently pairwise linked. The orthologous mouse map positions of 142 of 159 early discovered colon and lung cancer susceptibility genes formed 41 genomic clusters conserved >70 million years. These linked gene pairs concordantly affected both tumors and their majority was linked with two other gene families - protein tyrosine phosphatases and cancer driver protein kinases. 25% of both protein tyrosine phosphatases and protein kinases mapped <1 cM from a colon or lung cancer susceptibility gene, and 50% in <3 cM. Similar linkage was detected with most other human susceptibility genes that controlled 29 different cancer types. This concentration of tumor susceptibility genes with protein tyrosine phosphatases and driver protein kinases in multiple relatively short genomic regions suggests their possible functional diversity.

Identifying metabolic features of colorectal cancer liability using Mendelian randomization

Background

Recognizing the early signs of cancer risk is vital for informing prevention, early detection, and survival.

Methods

To investigate whether changes in circulating metabolites characterise the early stages of colorectal cancer (CRC) development, we examined associations between a genetic risk score (GRS) associated with CRC liability (72 single nucleotide polymorphisms) and 231 circulating metabolites measured by nuclear magnetic resonance spectroscopy in the Avon Longitudinal Study of Parents and Children (N=6,221). Linear regression models were applied to examine associations between genetic liability to colorectal cancer and circulating metabolites measured in the same individuals at age 8, 16, 18 and 25 years.

Results

The GRS for CRC was associated with up to 28% of the circulating metabolites at FDR-P<0.05 across all time points, particularly with higher fatty acids and very-low- and low-density lipoprotein subclass lipids. Two-sample reverse Mendelian randomization (MR) analyses investigating CRC liability (52,775 cases, 45,940 controls) and metabolites measured in a random subset of UK Biobank participants (N=118,466, median age 58y) revealed broadly consistent effect estimates with the GRS analysis. In conventional (forward) MR analyses, genetically predicted polyunsaturated fatty acid concentrations were most strongly associated with higher CRC risk.

Conclusions

These analyses suggest that higher genetic liability to CRC can cause early alterations in systemic metabolism, and suggest that fatty acids may play an important role in CRC development.

Funding

This work was supported by the Elizabeth Blackwell Institute for Health Research, University of Bristol, the Wellcome Trust, the Medical Research Council, Diabetes UK, the University of Bristol NIHR Biomedical Research Centre, and Cancer Research UK. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. This work used the computational facilities of the Advanced Computing Research Centre, University of Bristol - http://www.bristol.ac.uk/acrc/.

Large-scale alternative polyadenylation (APA)-wide association studies to identify putative susceptibility genes in human common cancers

Alternative polyadenylation (APA) modulates mRNA processing in the 3' untranslated regions (3'UTR), which affect mRNA stability and translation efficiency. Here, we build genetic models to predict APA levels in multiple tissues using sequencing data of 1,337 samples from the Genotype-Tissue Expression, and apply these models to assess associations between genetically predicted APA levels and cancer risk with data from large genome-wide association studies of six common cancers, including breast, ovary, prostate, colorectum, lung, and pancreas among European-ancestry populations. At a Bonferroni-corrected P □<□0.05, we identify 58 risk genes, including seven in newly identified loci. Using luciferase reporter assays, we demonstrate that risk alleles of 3'UTR variants, rs324015 ( STAT6 ), rs2280503 ( DIP2B ), rs1128450 ( FBXO38 ) and rs145220637 ( LDAH ), could significantly increase post-transcriptional activities of their target genes compared to reference alleles. Further gene knockdown experiments confirm their oncogenic roles. Our study provides additional insight into the genetic susceptibility of these common cancers.

Genetic Variants That Impact Alternative Polyadenylation in Cancer Represent Candidate Causal Risk Loci

Alternative polyadenylation (APA) is emerging as a major mechanism of posttranscriptional regulation. APA can impact the development and progression of cancer, suggesting that the genetic determinants of APA might play an important role in regulating cancer risk. Here, we depicted a pan-cancer atlas of human APA quantitative trait loci (apaQTL), containing approximately 0.7 million apaQTLs across 32 cancer types. Systematic multiomics analyses indicated that cancer apaQTLs could contribute to APA regulation by altering poly(A) motifs, RNA-binding proteins (RBP), and chromatin regulatory elements and were preferentially enriched in genome-wide association studies (GWAS)-identified cancer susceptibility loci. Moreover, apaQTL-related genes (aGene) were broadly related to cancer signaling pathways, high mutational burden, immune infiltration, and drug response, implicating their potential as therapeutic targets. Furthermore, apaQTLs were mapped in Chinese colorectal cancer tumor tissues and then screened for functional apaQTLs associated with colorectal cancer risk in 17,789 cases and 19,951 controls using GWAS-ChIP data, with independent validation in a large-scale population consisting of 6,024 cases and 10,022 controls. A multi-ancestry-associated apaQTL variant rs1020670 with a C>G change in DNM1L was identified, and the G allele contributed to an increased risk of colorectal cancer. Mechanistically, the risk variant promoted aberrant APA and facilitated higher usage of DNM1L proximal poly(A) sites mediated by the RBP CSTF2T, which led to higher expression of DNM1L with a short 3'UTR. This stabilized DNM1L to upregulate its expression, provoking colorectal cancer cell proliferation. Collectively, these findings generate a resource for understanding APA regulation and the genetic basis of human cancers, providing insights into cancer etiology.

SIGNIFICANCE

Cancer risk is mediated by alternative polyadenylation quantitative trait loci, including the rs1020670-G variant that promotes alternative polyadenylation of DNM1L and increases colorectal cancer risk.

Ability of a polygenic risk score to refine colorectal cancer risk in Lynch syndrome

BACKGROUND

Polygenic risk scores (PRSs) have been used to stratify colorectal cancer (CRC) risk in the general population, whereas its role in Lynch syndrome (LS), the most common type of hereditary CRC, is still conflicting. We aimed to assess the ability of PRS to refine CRC risk prediction in European-descendant individuals with LS.

METHODS

1465 individuals with LS (557 MLH1, 517 MSH2/EPCAM, 299 MSH6 and 92 PMS2) and 5656 CRC-free population-based controls from two independent cohorts were included. A 91-SNP PRS was applied. A Cox proportional hazard regression model with 'family' as a random effect and a logistic regression analysis, followed by a meta-analysis combining both cohorts were conducted.

RESULTS

Overall, we did not observe a statistically significant association between PRS and CRC risk in the entire cohort. Nevertheless, PRS was significantly associated with a slightly increased risk of CRC or advanced adenoma (AA), in those with CRC diagnosed <50 years and in individuals with multiple CRCs or AAs diagnosed <60 years.

CONCLUSION

The PRS may slightly influence CRC risk in individuals with LS in particular in more extreme phenotypes such as early-onset disease. However, the study design and recruitment strategy strongly influence the results of PRS studies. A separate analysis by genes and its combination with other genetic and non-genetic risk factors will help refine its role as a risk modifier in LS.