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.

Genome-wide interaction analysis of folate for colorectal cancer risk

BACKGROUND

Epidemiological and experimental evidence suggests that higher folate intake is associated with decreased colorectal cancer (CRC) risk; however, the mechanisms underlying this relationship are not fully understood. Genetic variation that may have a direct or indirect impact on folate metabolism can provide insights into folate's role in CRC.

OBJECTIVES

Our aim was to perform a genome-wide interaction analysis to identify genetic variants that may modify the association of folate on CRC risk.

METHODS

We applied traditional case-control logistic regression, joint 3-degree of freedom, and a 2-step weighted hypothesis approach to test the interactions of common variants (allele frequency >1%) across the genome and dietary folate, folic acid supplement use, and total folate in relation to risk of CRC in 30,550 cases and 42,336 controls from 51 studies from 3 genetic consortia (CCFR, CORECT, GECCO).

RESULTS

Inverse associations of dietary, total folate, and folic acid supplement with CRC were found (odds ratio [OR]: 0.93; 95% confidence interval [CI]: 0.90, 0.96; and 0.91; 95% CI: 0.89, 0.94 per quartile higher intake, and 0.82 (95% CI: 0.78, 0.88) for users compared with nonusers, respectively). Interactions (P-interaction < 5×10-8) of folic acid supplement and variants in the 3p25.2 locus (in the region of Synapsin II [SYN2]/tissue inhibitor of metalloproteinase 4 [TIMP4]) were found using traditional interaction analysis, with variant rs150924902 (located upstream to SYN2) showing the strongest interaction. In stratified analyses by rs150924902 genotypes, folate supplementation was associated with decreased CRC risk among those carrying the TT genotype (OR: 0.82; 95% CI: 0.79, 0.86) but increased CRC risk among those carrying the TA genotype (OR: 1.63; 95% CI: 1.29, 2.05), suggesting a qualitative interaction (P-interaction = 1.4×10-8). No interactions were observed for dietary and total folate.

CONCLUSIONS

Variation in 3p25.2 locus may modify the association of folate supplement with CRC risk. Experimental studies and studies incorporating other relevant omics data are warranted to validate this finding.

Genome-wide study of genetic polymorphisms predictive for outcome from first-line oxaliplatin-based chemotherapy in colorectal cancer patients

We conducted the first large genome-wide association study to identify novel genetic variants that predict better (or poorer) prognosis in colorectal cancer patients receiving standard first-line oxaliplatin-based chemotherapy vs chemotherapy without oxaliplatin. We used data from two phase III trials, NCCTG N0147 and NCCTG N9741 and a population-based patient cohort, DACHS. Multivariable Cox proportional hazards models were employed, including an interaction term between each SNP and type of treatment for overall survival (OS) and progression-free survival. The analysis was performed for studies individually, and the results were combined using fixed-effect meta-analyses separately for resected stage III colon cancer (3098 patients from NCCTG N0147 and 549 patients from DACHS) and mCRC (505 patients from NCCTG N9741 and 437 patients from DACHS). We further performed gene-based analysis as well as in silico bioinformatics analysis for CRC-relevant functional genomic annotation of identified loci. In stage III colon cancer patients, a locus on chr22 (rs11912167) was associated with significantly poorer OS after oxaliplatin-based chemotherapy vs chemotherapy without oxaliplatin (Pinteraction  < 5 × 10-8 ). For mCRC patients, three loci on chr1 (rs1234556), chr12 (rs11052270) and chr15 (rs11858406) were found to be associated with differential OS (P < 5 × 10-7 ). The locus on chr1 located in the intronic region of RCSD1 was replicated in an independent cohort of 586 mCRC patients from ALGB/SWOG 80405 (Pinteraction  = .04). The GWA gene-based analysis yielded for RCSD1 the most significant association with differential OS in mCRC (P = 6.6 × 10-6 ). With further investigation into its biological mechanisms, this finding could potentially be used to individualize first-line treatment and improve clinical outcomes.

Combining Asian and European genome-wide association studies of colorectal cancer improves risk prediction across racial and ethnic populations

Polygenic risk scores (PRS) have great potential to guide precision colorectal cancer (CRC) prevention by identifying those at higher risk to undertake targeted screening. However, current PRS using European ancestry data have sub-optimal performance in non-European ancestry populations, limiting their utility among these populations. Towards addressing this deficiency, we expand PRS development for CRC by incorporating Asian ancestry data (21,731 cases; 47,444 controls) into European ancestry training datasets (78,473 cases; 107,143 controls). The AUC estimates (95% CI) of PRS are 0.63(0.62-0.64), 0.59(0.57-0.61), 0.62(0.60-0.63), and 0.65(0.63-0.66) in independent datasets including 1681-3651 cases and 8696-115,105 controls of Asian, Black/African American, Latinx/Hispanic, and non-Hispanic White, respectively. They are significantly better than the European-centric PRS in all four major US racial and ethnic groups (p-values < 0.05). Further inclusion of non-European ancestry populations, especially Black/African American and Latinx/Hispanic, is needed to improve the risk prediction and enhance equity in applying PRS in clinical practice.

Assessing efficiency of fine-mapping obesity-associated variants through leveraging ancestry architecture and functional annotation using PAGE and UKBB cohorts

Inadequate representation of non-European ancestry populations in genome-wide association studies (GWAS) has limited opportunities to isolate functional variants. Fine-mapping in multi-ancestry populations should improve the efficiency of prioritizing variants for functional interrogation. To evaluate this hypothesis, we leveraged ancestry architecture to perform comparative GWAS and fine-mapping of obesity-related phenotypes in European ancestry populations from the UK Biobank (UKBB) and multi-ancestry samples from the Population Architecture for Genetic Epidemiology (PAGE) consortium with comparable sample sizes. In the investigated regions with genome-wide significant associations for obesity-related traits, fine-mapping in our ancestrally diverse sample led to 95% and 99% credible sets (CS) with fewer variants than in the European ancestry sample. Lead fine-mapped variants in PAGE regions had higher average coding scores, and higher average posterior probabilities for causality compared to UKBB. Importantly, 99% CS in PAGE loci contained strong expression quantitative trait loci (eQTLs) in adipose tissues or harbored more variants in tighter linkage disequilibrium (LD) with eQTLs. Leveraging ancestrally diverse populations with heterogeneous ancestry architectures, coupled with functional annotation, increased fine-mapping efficiency and performance, and reduced the set of candidate variants for consideration for future functional studies. Significant overlap in genetic causal variants across populations suggests generalizability of genetic mechanisms underpinning obesity-related traits across populations.

Effect of tetracycline hydrochloride application on dental pulp stem cell metabolism-booster or obstacle for tissue engineering?

Introduction: Stem cells and scaffolds are an important foundation and starting point for tissue engineering. Human dental pulp stem cells (DPSC) are mesenchymal stem cells with self-renewal and multi-directional differentiation potential, and are ideal candidates for tissue engineering due to their excellent biological properties and accessibility without causing major trauma at the donor site. Tetracycline hydrochloride (TCH), a broad-spectrum antibiotic, has been widely used in recent years for the synthesis of cellular scaffolds to reduce the incidence of postoperative infections. Methods: In order to evaluate the effects of TCH on DPSC, the metabolism of DPSC in different concentrations of TCH environment was tested. Moreover, cell morphology, survival rates, proliferation rates, cell migration rates and differentiation abilities of DPSC at TCH concentrations of 0-500 μg/ml were measured. Phalloidin staining, live-dead staining, MTS assay, cell scratch assay and real-time PCR techniques were used to detect the changes in DPSC under varies TCH concentrations. Results: At TCH concentrations higher than 250 μg/ml, DPSC cells were sequestered, the proportion of dead cells increased, and the cell proliferation capacity and cell migration capacity decreased. The osteogenic and adipogenic differentiation abilities of DPSC, however, were already inhibited at TCH con-centrations higher than 50 μg/ml. Here, the expression of the osteogenic genes, runt-related transcription factor 2 (RUNX2) and osteocalcin (OCN), the lipogenic genes lipase (LPL), as well as the peroxisome proliferator-activated receptor-γ (PPAR-γ) expression were found to be down-regulated. Discussion: The results of the study indicated that TCH in concentrations above 50 µg/ml negatively affects the differentiation capability of DPSC. In addition, TCH at concentrations above 250 µg/ml adversely affects the growth status, percentage of living cells, proliferation and migration ability of cells.

Whole Genome Sequencing Based Analysis of Inflammation Biomarkers in the Trans-Omics for Precision Medicine (TOPMed) Consortium

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38,465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program. We identified 22 distinct single-variant associations across 6 traits - E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin - that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

Genetic Susceptibility to Nonalcoholic Fatty Liver Disease and Risk for Pancreatic Cancer: Mendelian Randomization

BACKGROUND

There are conflicting data on whether nonalcoholic fatty liver disease (NAFLD) is associated with susceptibility to pancreatic cancer. Using Mendelian randomization (MR), we investigated the relationship between genetic predisposition to NAFLD and risk for pancreatic cancer.

METHODS

Data from genome-wide association studies (GWAS) within the Pancreatic Cancer Cohort Consortium (PanScan; cases n = 5,090, controls n = 8,733) and the Pancreatic Cancer Case Control Consortium (PanC4; cases n = 4,163, controls n = 3,792) were analyzed. We used data on 68 genetic variants with four different MR methods [inverse variance weighting (IVW), MR-Egger, simple median, and penalized weighted median] separately to predict genetic heritability of NAFLD. We then assessed the relationship between each of the four MR methods and pancreatic cancer risk, using logistic regression to calculate ORs and 95% confidence intervals (CI), adjusting for PC risk factors, including obesity and diabetes.

RESULTS

No association was found between genetically predicted NAFLD and pancreatic cancer risk in the PanScan or PanC4 samples [e.g., PanScan, IVW OR, 1.04; 95% confidence interval (CI), 0.88-1.22; MR-Egger OR, 0.89; 95% CI, 0.65-1.21; PanC4, IVW OR, 1.07; 95% CI, 0.90-1.27; MR-Egger OR, 0.93; 95% CI, 0.67-1.28]. None of the four MR methods indicated an association between genetically predicted NAFLD and pancreatic cancer risk in either sample.

CONCLUSIONS

Genetic predisposition to NAFLD is not associated with pancreatic cancer risk.

IMPACT

Given the close relationship between NAFLD and metabolic conditions, it is plausible that any association between NAFLD and pancreatic cancer might reflect host metabolic perturbations (e.g., obesity, diabetes, or metabolic syndrome) and does not necessarily reflect a causal relationship between NAFLD and pancreatic cancer.

Microbiological and Clinical Effects of a Proanthocyanidin-enriched Extract from Rumex acetosa in Periodontally Healthy Carriers of Porphyromonas gingivalis: a Randomized Controlled Pilot Study

Rumex acetosa significantly inhibits the adhesion of Porphyromonas gingivalis (P. g.) to eukaryotic host cells in vitro. The objective of this randomized placebo-controlled pilot-trial was to analyze effects of a mouth rinse containing 0.8% (w/w) of a quantified proanthocyanidin-enriched extract from Rumex acetosa (RA1) on microbiological, clinical, and cytological parameters in systemically healthy individuals without history of periodontitis, harboring P. g. intraorally. 35 subjects received a supragingival debridement (SD) followed by mouth rinsing (3 times daily) with either RA1 mouth rinse solution (test) or placebo (control) for 7 days as adjunct to routine oral hygiene. Supragingival biofilm samples were taken at screening visit, baseline (BL), 2, 4, 7 and 14 days after SD. P. g. and 11 other oral microorganisms were detected and quantified by rtPCR. Changes in the oral microbiota composition of one test and one control subject were assessed via high throughput 16S rRNS gene amplicon sequencing. Approximal Plaque Index (API) and the modified Sulcular Bleeding Index (SBI) were assessed at BL, 7- and 14-days following SD. Brush biopsies were taken at BL and 14 d following SD. Intergroup comparisons revealed no significant microbiological, cytological, and clinical differences at any timepoint. However, a significant reduction in SBI at day 14 (p = 0.003) and API at day 7 (p = 0.02) and day 14 (p = 0.009) was found in the test group by intragroup comparison. No severe adverse events were observed. The results indicate that RA1 mouth rinse is safe but does not seem to inhibit colonization of P. g. or improve periodontal health following SD.

Using DEPendency of Association on the Number of Top Hits (DEPTH) as a Complementary Tool to Identify Novel Colorectal Cancer Susceptibility Loci

BACKGROUND

DEPendency of association on the number of Top Hits (DEPTH) is an approach to identify candidate susceptibility regions by considering the risk signals from overlapping groups of sequential variants across the genome.

METHODS

We applied a DEPTH analysis using a sliding window of 200 SNPs to colorectal cancer data from the Colon Cancer Family Registry (CCFR; 5,735 cases and 3,688 controls), and Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO; 8,865 cases and 10,285 controls) studies. A DEPTH score > 1 was used to identify candidate susceptibility regions common to both analyses. We compared DEPTH results against those from conventional genome-wide association study (GWAS) analyses of these two studies as well as against 132 published susceptibility regions.

RESULTS

Initial DEPTH analysis revealed 2,622 (CCFR) and 3,686 (GECCO) candidate susceptibility regions, of which 569 were common to both studies. Bootstrapping revealed 40 and 49 candidate susceptibility regions in the CCFR and GECCO data sets, respectively. Notably, DEPTH identified at least 82 regions that would not be detected using conventional GWAS methods, nor had they been identified by previous colorectal cancer GWASs. We found four reproducible candidate susceptibility regions (2q22.2, 2q33.1, 6p21.32, 13q14.3). The highest DEPTH scores were in the human leukocyte antigen locus at 6p21 where the strongest associated SNPs were rs762216297, rs149490268, rs114741460, and rs199707618 for the CCFR data, and rs9270761 for the GECCO data.

CONCLUSIONS

DEPTH can identify candidate susceptibility regions for colorectal cancer not identified using conventional analyses of larger datasets.

IMPACT

DEPTH has potential as a powerful complementary tool to conventional GWAS analyses for discovering susceptibility regions within the genome.

WHOLE GENOME SEQUENCING ANALYSIS OF BODY MASS INDEX IDENTIFIES NOVEL AFRICAN ANCESTRY-SPECIFIC RISK ALLELE

Obesity is a major public health crisis associated with high mortality rates. Previous genome-wide association studies (GWAS) investigating body mass index (BMI) have largely relied on imputed data from European individuals. This study leveraged whole-genome sequencing (WGS) data from 88,873 participants from the Trans-Omics for Precision Medicine (TOPMed) Program, of which 51% were of non-European population groups. We discovered 18 BMI-associated signals (P < 5 × 10-9). Notably, we identified and replicated a novel low frequency single nucleotide polymorphism (SNP) in MTMR3 that was common in individuals of African descent. Using a diverse study population, we further identified two novel secondary signals in known BMI loci and pinpointed two likely causal variants in the POC5 and DMD loci. Our work demonstrates the benefits of combining WGS and diverse cohorts in expanding current catalog of variants and genes confer risk for obesity, bringing us one step closer to personalized medicine.

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

BACKGROUND

Diabetes is an established risk factor for colorectal cancer. However, the mechanisms underlying this relationship still require investigation and it is not known if the association is modified by genetic variants. To address these questions, we undertook a genome-wide gene-environment interaction analysis.

METHODS

We used data from 3 genetic consortia (CCFR, CORECT, GECCO; 31,318 colorectal cancer cases/41,499 controls) and undertook genome-wide gene-environment interaction analyses with colorectal cancer risk, including interaction tests of genetics(G)xdiabetes (1-degree of freedom; d.f.) and joint testing of Gxdiabetes, G-colorectal cancer association (2-d.f. joint test) and G-diabetes correlation (3-d.f. joint test).

RESULTS

Based on the joint tests, we found that the association of diabetes with colorectal cancer risk is modified by loci on chromosomes 8q24.11 (rs3802177, SLC30A8 - ORAA: 1.62, 95% CI: 1.34-1.96; ORAG: 1.41, 95% CI: 1.30-1.54; ORGG: 1.22, 95% CI: 1.13-1.31; p-value3-d.f.: 5.46 × 10-11) and 13q14.13 (rs9526201, LRCH1 - ORGG: 2.11, 95% CI: 1.56-2.83; ORGA: 1.52, 95% CI: 1.38-1.68; ORAA: 1.13, 95% CI: 1.06-1.21; p-value2-d.f.: 7.84 × 10-09).

DISCUSSION

These results suggest that variation in genes related to insulin signaling (SLC30A8) and immune function (LRCH1) may modify the association of diabetes with colorectal cancer risk and provide novel insights into the biology underlying the diabetes and colorectal cancer relationship.

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

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

SIGNIFICANCE

This gene-environment interaction analysis revealed a genetic locus in FMN1/GREM1 that interacts with body mass index in colorectal cancer risk, suggesting potential implications for precision prevention strategies.

Ancestral diversity in lipoprotein(a) studies helps address evidence gaps

INTRODUCTION

The independent and causal cardiovascular disease risk factor lipoprotein(a) (Lp(a)) is elevated in >1.5 billion individuals worldwide, but studies have prioritised European populations.

METHODS

Here, we examined how ancestrally diverse studies could clarify Lp(a)'s genetic architecture, inform efforts examining application of Lp(a) polygenic risk scores (PRS), enable causal inference and identify unexpected Lp(a) phenotypic effects using data from African (n=25 208), East Asian (n=2895), European (n=362 558), South Asian (n=8192) and Hispanic/Latino (n=8946) populations.

RESULTS

Fourteen genome-wide significant loci with numerous population specific signals of large effect were identified that enabled construction of Lp(a) PRS of moderate (R2=15% in East Asians) to high (R2=50% in Europeans) accuracy. For all populations, PRS showed promise as a 'rule out' for elevated Lp(a) because certainty of assignment to the low-risk threshold was high (88.0%-99.9%) across PRS thresholds (80th-99th percentile). Causal effects of increased Lp(a) with increased glycated haemoglobin were estimated for Europeans (p value =1.4×10-6), although inverse effects in Africans and East Asians suggested the potential for heterogeneous causal effects. Finally, Hispanic/Latinos were the only population in which known associations with coronary atherosclerosis and ischaemic heart disease were identified in external testing of Lp(a) PRS phenotypic effects.

CONCLUSIONS

Our results emphasise the merits of prioritising ancestral diversity when addressing Lp(a) evidence gaps.

Selection, optimization, and validation of ten chronic disease polygenic risk scores for clinical implementation in diverse populations

Polygenic risk scores (PRS) have improved in predictive performance supporting their use in clinical practice. Reduced predictive performance of PRS in diverse populations can exacerbate existing health disparities. The NHGRI-funded eMERGE Network is returning a PRS-based genome-informed risk assessment to 25,000 diverse adults and children. We assessed PRS performance, medical actionability, and potential clinical utility for 23 conditions. Standardized metrics were considered in the selection process with additional consideration given to strength of evidence in African and Hispanic populations. Ten conditions were selected with a range of high-risk thresholds: atrial fibrillation, breast cancer, chronic kidney disease, coronary heart disease, hypercholesterolemia, prostate cancer, asthma, type 1 diabetes, obesity, and type 2 diabetes. We developed a pipeline for clinical PRS implementation, used genetic ancestry to calibrate PRS mean and variance, created a framework for regulatory compliance, and developed a PRS clinical report. eMERGE's experience informs the infrastructure needed to implement PRS-based implementation in diverse clinical settings.

The association between genetically elevated polyunsaturated fatty acids and risk of cancer

BACKGROUND

The causal relevance of polyunsaturated fatty acids (PUFAs) for risk of site-specific cancers remains uncertain.

METHODS

Using a Mendelian randomization (MR) framework, we assessed the causal relevance of PUFAs for risk of cancer in European and East Asian ancestry individuals. We defined the primary exposure as PUFA desaturase activity, proxied by rs174546 at the FADS locus. Secondary exposures were defined as omega 3 and omega 6 PUFAs that could be proxied by genetic polymorphisms outside the FADS region. Our study used summary genetic data on 10 PUFAs and 67 cancers, corresponding to 562,871 cases and 1,619,465 controls, collected by the Fatty Acids in Cancer Mendelian Randomization Collaboration. We estimated odds ratios (ORs) for cancer per standard deviation increase in genetically proxied PUFA exposures.

FINDINGS

Genetically elevated PUFA desaturase activity was associated (P < 0.0007) with higher risk (OR [95% confidence interval]) of colorectal cancer (1.09 [1.07-1.11]), esophageal squamous cell carcinoma (1.16 [1.06-1.26]), lung cancer (1.06 [1.03-1.08]) and basal cell carcinoma (1.05 [1.02-1.07]). There was little evidence for associations with reproductive cancers (OR = 1.00 [95% CI: 0.99-1.01]; Pheterogeneity = 0.25), urinary system cancers (1.03 [0.99-1.06], Pheterogeneity = 0.51), nervous system cancers (0.99 [0.95-1.03], Pheterogeneity = 0.92) or blood cancers (1.01 [0.98-1.04], Pheterogeneity = 0.09). Findings for colorectal cancer and esophageal squamous cell carcinoma remained compatible with causality in sensitivity analyses for violations of assumptions. Secondary MR analyses highlighted higher omega 6 PUFAs (arachidonic acid, gamma-linolenic acid and dihomo-gamma-linolenic acid) as potential mediators. PUFA biosynthesis is known to interact with aspirin, which increases risk of bleeding and inflammatory bowel disease. In a phenome-wide MR study of non-neoplastic diseases, we found that genetic lowering of PUFA desaturase activity, mimicking a hypothetical intervention to reduce cancer risk, was associated (P < 0.0006) with increased risk of inflammatory bowel disease but not bleeding.

INTERPRETATION

The PUFA biosynthesis pathway may be an intervention target for prevention of colorectal cancer and esophageal squamous cell carcinoma but with potential for increased risk of inflammatory bowel disease.

FUNDING

Cancer Resesrch UK (C52724/A20138, C18281/A19169). UK Medical Research Council (MR/P014054/1). National Institute for Health Research (NIHR202411). UK Medical Research Council (MC_UU_00011/1, MC_UU_00011/3, MC_UU_00011/6, and MC_UU_00011/4). National Cancer Institute (R00 CA215360). National Institutes of Health (U01 CA164973, R01 CA60987, R01 CA72520, U01 CA74806, R01 CA55874, U01 CA164973 and U01 CA164973).

Application of Mendelian randomization to explore the causal role of the human gut microbiome in colorectal cancer

The role of the human gut microbiome in colorectal cancer (CRC) is unclear as most studies on the topic are unable to discern correlation from causation. We apply two-sample Mendelian randomization (MR) to estimate the causal relationship between the gut microbiome and CRC. We used summary-level data from independent genome-wide association studies to estimate the causal effect of 14 microbial traits (n = 3890 individuals) on overall CRC (55,168 cases, 65,160 controls) and site-specific CRC risk, conducting several sensitivity analyses to understand the nature of results. Initial MR analysis suggested that a higher abundance of Bifidobacterium and presence of an unclassified group of bacteria within the Bacteroidales order in the gut increased overall and site-specific CRC risk. However, sensitivity analyses suggested that instruments used to estimate relationships were likely complex and involved in many potential horizontal pleiotropic pathways, demonstrating that caution is needed when interpreting MR analyses with gut microbiome exposures. In assessing reverse causality, we did not find strong evidence that CRC causally affected these microbial traits. Whilst our study initially identified potential causal roles for two microbial traits in CRC, importantly, further exploration of these relationships highlighted that these were unlikely to reflect causality.

Abstract 5223: Circulating lipoprotein lipids and colorectal cancer risk: A Mendelian randomization analysis from the GECCO consortium

Background Conventional observational studies have reported conflicting results regarding the association between low density lipoprotein cholesterol (LDL-C) and risk of colorectal cancer (CRC). We conducted a Mendelian randomization analysis to address this association.

Methods Single-nucleotide polymorphisms (SNPs) associated with five blood lipids (total cholesterol, HDL-C, high-density lipoprotein cholesterol [HDL-C], non-HDL-C, and triglyceride) were obtained from a genome-wide association study (GWAS) meta-analysis of European ancestry in the Global Lipids Genetics Consortium (GLGC, N ≤ 1 319 982), and two lipids (apolipoprotein A1 and apolipoprotein B) from a GWAS in the UK Biobank (UKB, N ≤ 441 016). Summary statistics were obtained for these SNPs from a GWAS of CRC in the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO) including 34,869 cases and 29,051 controls. Associations with CRC risk per one standard deviation increase in the genetically predicted lipids level were generated using inverse-variance weighted random-effects models.

Results No overall association was observed between genetically predicted levels of blood lipids and CRC risk. However, increased risks were observed for all LDL-C related traits among women, including total cholesterol (OR = 1.10; 95%CI = 1.02, 1.18), LDL-C (1.07; 1.00, 1.14), non-HDL-C (1.09; 1.02, 1.16), and apolipoprotein B (1.11; 1.02, 1.21); whereas no significant association was found among men. Similar but ostensibly stronger associations of these traits

were seen with distal colon cancer cases, with no significant association showing on proximal colon or rectum cancer cases. We also observed similar positive associations of LDL-C related traits among those having CRC before their 50 years. Of note, risk reduction was found for apolipoprotein A1, a major component of HDL-C, in these early-onset cases (0.87; 0.78, 0.98).

Conclusion Results from this study suggest that high circulating LDL-C levels may increase the risk of CRC, particularly cancer of the distal colon, and the association may differ by sex and age at CRC onset.

Key words: Blood lipids; colorectal cancer; Mendelian randomization.

Citation Format: Lili Liu, Wanqing Wen, Jirong Long, Themistocles L Assimes, Luis Bujanda, Stephen B Gruber, Sébastien Küry, Brigid Lynch, Conghui Qu, Minta Thomas, Emily White, Michael O. Woods, Ulrike Peters, Wei Zheng. Circulating lipoprotein lipids and colorectal cancer risk: A Mendelian randomization analysis from the GECCO consortium. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 5223.

Abstract 5760: Molecular, immune, and microbial profiles of early-onset, intermediate-onset, and later-onset CRCs

Background: Despite heightened interest in early-onset CRC, little is known about the tumor molecular, immune, and microbial characteristics of early-onset CRC. It is also unclear whether CRCs diagnosed at or shortly after age 50 are similar to early-onset CRC. We hypothesized that tumor molecular, immune, and microbial characteristics in CRC tissue might show differential heterogeneity patterns between three age groups (&lt;50 "early-onset", 50-54 "intermediate-onset", ≥55 "later-onset").

Methods: We examined 3,395 CRC cases with available tissue data, including 660 early-onset and 243 intermediate-onset cases in the Nurses' Health Study, Health Professionals Follow-up Study, and Ontario Familial Colon Cancer Registry. We profiled the in situ T-cell landscape of 959 cases using digital imaging, machine learning, and a customized 9-plex multiplexed immunofluoresence panel with antibodies directed against CD3, CD4, CD8, CD45RA, CD45RO, FOXP3, and MKI67 (Ki-67). Using chi-squared test or Spearman’s correlation test, we assessed differences in tumor characteristics including microsatellite instability, CpG island methylator phenotype (CIMP), KRAS, and BRAF mutations, LINE-1 methylation levels, pks+ E. coli and Fusobacterium nucleatum positivity, histopathologic lymphocytic reaction, and T-cell densities between early-onset, intermediate-onset, and later-onset cases.

Results: Compared to later-onset CRC (21%), early-onset (4.8%) and intermediate-onset CRCs (5.1%) exhibited a lower prevalence of CIMP-high status (P&lt;0.001). The mean tumor LINE-1 methylation level increased with increasing age [59 (SD, 12) in early-onset, 61 (SD, 10) in intermediate-onset, and 64 (SD, 9.8) in later-onset CRC (P&lt;0.001)]. Early-onset CRC (4.1%) had fewer BRAF mutations than intermediate-onset (10.4%) and later-onset CRC (15%) (P&lt;0.001). Compared to later-onset CRC, early-onset CRC tended to show lower levels of tumor-infiltrating lymphocytes (P=0.013), and intratumoral periglandular reaction (P=0.025). Compared to later-onset CRC, early-onset CRC had a lower density of memory (both CD4+ and CD8+) T-cells (median 21 vs. 48 cells/mm2; P=0.002) and a higher density of MKI67+ immune cells (median 50 vs. 19 cells/mm2; P=0.003) in tumor epithelial areas but not in the surrounding stroma. No significant differences in other tested characteristics were identified.

Conclusions: Compared to later-onset CRC, early-onset and intermediate-onset CRCs tended to have aggressive tumor phenotypes such as LINE-1 hypomethylation, lower lymphocytic immune reaction, and a higher density of MKI67+ immune cells. These findings highlight the importance of the tumor microenvironment in the etiology of early-onset and intermediate-onset CRCs.

Citation Format: Tomotaka Ugai, Yasutoshi Takashima, Andressa Dias Costa, Daniel Buchanan, Jeroen Huyghe, Li Hsu, Conghui Qu, Claire Thomas, Steve Gallinger, Robert Grant, Ulrike Peters, Amanda I. Phipps, Jonathan Nowak, Shuji Ogino. Molecular, immune, and microbial profiles of early-onset, intermediate-onset, and later-onset CRCs. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 5760.

Molecular Characteristics of Early-Onset Colorectal Cancer According to Detailed Anatomical Locations: Comparison With Later-Onset Cases

INTRODUCTION

Early-onset colorectal cancer diagnosed before the age of 50 years has been increasing. Likely reflecting the pathogenic role of the intestinal microbiome, which gradually changes across the entire colorectal length, the prevalence of certain tumor molecular characteristics gradually changes along colorectal subsites. Understanding how colorectal tumor molecular features differ by age and tumor location is important in personalized patient management.

METHODS

Using 14,004 cases with colorectal cancer including 3,089 early-onset cases, we examined microsatellite instability (MSI), CpG island methylator phenotype (CIMP), and KRAS and BRAF mutations in carcinomas of the cecum, ascending colon, transverse colon, descending colon, sigmoid colon, and rectum and compared early-onset cases with later-onset cases.

RESULTS

The proportions of MSI-high, CIMP-high, and BRAF -mutated early-onset tumors were lowest in the rectum (8.8%, 3.4%, and 3.5%, respectively) and highest in the ascending colon (46% MSI-high; 15% CIMP-high) or transverse colon (8.6% BRAF -mutated) (all Ptrend <0.001 across the rectum to ascending colon). Compared with later-onset tumors, early-onset tumors showed a higher prevalence of MSI-high status and a lower prevalence of CIMP-high status and BRAF mutations in most subsites. KRAS mutation prevalence was higher in the cecum compared with that in the other subsites in both early-onset and later-onset tumors ( P < 0.001). Notably, later-onset MSI-high tumors showed a continuous decrease in KRAS mutation prevalence from the rectum (36%) to ascending colon (9%; Ptrend <0.001), followed by an increase in the cecum (14%), while early-onset MSI-high cancers showed no such trend.

DISCUSSION

Our findings support biogeographical and pathogenic heterogeneity of colorectal carcinomas in different colorectal subsites and age groups.

Causal effects on complex traits are similar for common variants across segments of different continental ancestries within admixed individuals

Individuals of admixed ancestries (for example, African Americans) inherit a mosaic of ancestry segments (local ancestry) originating from multiple continental ancestral populations. This offers the unique opportunity of investigating the similarity of genetic effects on traits across ancestries within the same population. Here we introduce an approach to estimate correlation of causal genetic effects (radmix) across local ancestries and analyze 38 complex traits in African-European admixed individuals (N = 53,001) to observe very high correlations (meta-analysis radmix = 0.95, 95% credible interval 0.93-0.97), much higher than correlation of causal effects across continental ancestries. We replicate our results using regression-based methods from marginal genome-wide association study summary statistics. We also report realistic scenarios where regression-based methods yield inflated heterogeneity-by-ancestry due to ancestry-specific tagging of causal effects, and/or polygenicity. Our results motivate genetic analyses that assume minimal heterogeneity in causal effects by ancestry, with implications for the inclusion of ancestry-diverse individuals in studies.

Leisure time sedentary behaviour and risks of breast, colorectal, and prostate cancer: A Mendelian randomization analysis

Sedentary behaviours have been associated with increased risks of some common cancers in epidemiological studies; however, it is unclear if these associations are causal. We examined potential causal associations between self-reported leisure television watching and computer use and risks of breast, colorectal, and prostate cancer using a two-sample Mendelian randomization framework. Genetic variants were identified from a recent genome-wide association study (GWAS). Cancer data were obtained from cancer GWAS consortia. Additional sensitivity analyses were applied to examine the robustness of the results. A 1-standard deviation increment in hours of television watching increased risk of breast (OR: 1.15, 95% confidence interval [CI]: 1.05,1.26) and colorectal cancer (OR: 1.32, 95%CI: 1.16,1.49) with little evidence of an association for prostate cancer risk. In multivariable models adjusted for years of education, the effect estimates for television watching were attenuated (breast cancer, OR: 1.08, 95%CI: 0.92,1.27; colorectal cancer, OR: 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 by sex (colorectal cancer), anatomical subsites, and cancer subtypes. There was little evidence of associations between computer use and cancer risk. We found evidence of positive associations between hours of television watching and risks of breast and colorectal cancer. However, these findings should be interpreted cautiously given the complex role of education. Future studies using objective measures of exposure can provide new insights into the possible role of sedentary behaviour in cancer development.

Novelty and impact

Evidence from observational studies that examined associations between sedentary behaviours and common cancers is mixed and causality is uncertain. In our Mendelian randomization analyses, higher levels of leisure television watching were found to increase the risks of breast and colorectal cancer, suggesting that the that the promotion of lowering sedentary behaviour time could be an effective strategy in the primary prevention of these commonly diagnosed cancers.

Article category

Cancer Epidemiology.

Integrated multi-omics approach to distinct molecular characterization and classification of early-onset colorectal cancer

Incidence of early-onset colorectal cancer (EOCRC), defined by a diagnosed age under 50 years, is increasing, but its heterogeneous etiologies that differ from general CRC remain undetermined. We initially characterize the genome, epigenome, transcriptome, and proteome of tumors from 79 patients in a Chinese CRC cohort. Data for an additional 126 EOCRC subjects are obtained from the International Cancer Genome Consortium Chinese cohort and The Cancer Genome Atlas European cohort. We observe that early-onset tumors have a high tumor mutation burden; increased DNA repair features by mutational signature 3 and multi-layer pathway enrichments; strong perturbations at effects of DNA methylation and somatic copy-number alteration on gene expression; and upregulated immune infiltration as hot tumors underlying immunophenotypes. Notably, LMTK3 exhibits ancestral mutation disparity, potentially being a functional modulator and biomarker that drives molecular alterations in EOCRC development and immunotherapies. This integrative omics study provides valuable knowledge for precision oncology of CRC.