Large-scale Integrated Analysis of Genetics and Metabolomic Data Reveals Potential Links Between Lipids and Colorectal Cancer Risk

BACKGROUND

The etiology of colorectal cancer is not fully understood.

METHODS

Using genetic variants and metabolomics data including 217 metabolites from the Framingham Heart Study (n = 1,357), we built genetic prediction models for circulating metabolites. Models with prediction R2 > 0.01 (Nmetabolite = 58) were applied to predict levels of metabolites in two large consortia with a combined sample size of approximately 46,300 cases and 59,200 controls of European and approximately 21,700 cases and 47,400 controls of East Asian (EA) descent. Genetically predicted levels of metabolites were evaluated for their associations with colorectal cancer risk in logistic regressions within each racial group, after which the results were combined by meta-analysis.

RESULTS

Of the 58 metabolites tested, 24 metabolites were significantly associated with colorectal cancer risk [Benjamini-Hochberg FDR (BH-FDR) < 0.05] in the European population (ORs ranged from 0.91 to 1.06; P values ranged from 0.02 to 6.4 × 10-8). Twenty one of the 24 associations were replicated in the EA population (ORs ranged from 0.26 to 1.69, BH-FDR < 0.05). In addition, the genetically predicted levels of C16:0 cholesteryl ester was significantly associated with colorectal cancer risk in the EA population only (OREA: 1.94, 95% CI, 1.60-2.36, P = 2.6 × 10-11; OREUR: 1.01, 95% CI, 0.99-1.04, P = 0.3). Nineteen of the 25 metabolites were glycerophospholipids and triacylglycerols (TAG). Eighteen associations exhibited significant heterogeneity between the two racial groups (PEUR-EA-Het < 0.005), which were more strongly associated in the EA population. This integrative study suggested a potential role of lipids, especially certain glycerophospholipids and TAGs, in the etiology of colorectal cancer.

CONCLUSIONS

This study identified potential novel risk biomarkers for colorectal cancer by integrating genetics and circulating metabolomics data.

IMPACT

The identified metabolites could be developed into new tools for risk assessment of colorectal cancer in both European and EA populations.

Associations Between Glycemic Traits and Colorectal Cancer: A Mendelian Randomization Analysis

BACKGROUND

Glycemic traits-such as hyperinsulinemia, hyperglycemia, and type 2 diabetes-have been associated with higher colorectal cancer risk in observational studies; however, causality of these associations is uncertain. We used Mendelian randomization (MR) to estimate the causal effects of fasting insulin, 2-hour glucose, fasting glucose, glycated hemoglobin (HbA1c), and type 2 diabetes with colorectal cancer.

METHODS

Genome-wide association study summary data were used to identify genetic variants associated with circulating levels of fasting insulin (n = 34), 2-hour glucose (n = 13), fasting glucose (n = 70), HbA1c (n = 221), and type 2 diabetes (n = 268). Using 2-sample MR, we examined these variants in relation to colorectal cancer risk (48 214 case patient and 64 159 control patients).

RESULTS

In inverse-variance models, higher fasting insulin levels increased colorectal cancer risk (odds ratio [OR] per 1-SD = 1.65, 95% confidence interval [CI] = 1.15 to 2.36). We found no evidence of any effect of 2-hour glucose (OR per 1-SD = 1.02, 95% CI = 0.86 to 1.21) or fasting glucose (OR per 1-SD = 1.04, 95% CI = 0.88 to 1.23) concentrations on colorectal cancer risk. Genetic liability to type 2 diabetes (OR per 1-unit increase in log odds = 1.04, 95% CI = 1.01 to 1.07) and higher HbA1c levels (OR per 1-SD = 1.09, 95% CI = 1.00 to 1.19) increased colorectal cancer risk, although these findings may have been biased by pleiotropy. Higher HbA1c concentrations increased rectal cancer risk in men (OR per 1-SD = 1.21, 95% CI = 1.05 to 1.40), but not in women.

CONCLUSIONS

Our results support a causal effect of higher fasting insulin, but not glucose traits or type 2 diabetes, on increased colorectal cancer risk. This suggests that pharmacological or lifestyle interventions that lower circulating insulin levels may be beneficial in preventing colorectal tumorigenesis.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

IMPACT

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

Genetic Regulation of DNA Methylation Yields Novel Discoveries in GWAS of Colorectal Cancer

BACKGROUND

Colorectal cancer has a strong epigenetic component that is accompanied by frequent DNA methylation (DNAm) alterations in addition to heritable genetic risk. It is of interest to understand the interrelationship of germline genetics, DNAm, and colorectal cancer risk.

METHODS

We performed a genome-wide methylation quantitative trait locus (meQTL) analysis in 1,355 people, assessing the pairwise associations between genetic variants and lymphocytes methylation data. In addition, we used penalized regression with cis-genetic variants ± 1 Mb of methylation to identify genome-wide heritable DNAm. We evaluated the association of genetically predicted methylation with colorectal cancer risk based on genome-wide association studies (GWAS) of over 125,000 cases and controls using the multivariate sMiST as well as univariately via examination of marginal association with colorectal cancer risk.

RESULTS

Of the 142 known colorectal cancer GWAS loci, 47 were identified as meQTLs. We identified four novel colorectal cancer-associated loci (NID2, ATXN10, KLHDC10, and CEP41) that reside over 1 Mb outside of known colorectal cancer loci and 10 secondary signals within 1 Mb of known loci.

CONCLUSIONS

Leveraging information of DNAm regulation into genetic association of colorectal cancer risk reveals novel pathways in colorectal cancer tumorigenesis. Our summary statistics-based framework sMiST provides a powerful approach by combining information from the effect through methylation and residual direct effects of the meQTLs on disease risk. Further validation and functional follow-up of these novel pathways are needed.

IMPACT

Using genotype, DNAm, and GWAS, we identified four new colorectal cancer risk loci. We studied the landscape of genetic regulation of DNAm via single-SNP and multi-SNP meQTL analyses.

Genetic variants associated with circulating C-reactive protein levels and colorectal cancer survival: Sex-specific and lifestyle factors specific associations

Elevated blood levels of C-reactive protein (CRP) have been linked to colorectal cancer (CRC) survival. We evaluated genetic variants associated with CRP levels and their interactions with sex and lifestyle factors in association with CRC-specific mortality. Our study included 16 142 CRC cases from the International Survival Analysis in Colorectal Cancer Consortium. We identified 618 common single nucleotide polymorphisms (SNPs) associated with CRP levels from the NHGRI-EBI GWAS Catalog. Cox proportional hazards regression was used to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for associations between SNPs and CRC-specific mortality adjusting for age, sex, genotyping platform/study and principal components. We investigated their interactions with sex and lifestyle factors using likelihood ratio tests. Of 5472 (33.9%) deaths accrued over up to 10 years of follow-up, 3547 (64.8%) were due to CRC. No variants were associated with CRC-specific mortality after multiple comparison correction. We observed strong evidence of interaction between variant rs1933736 at FRK gene and sex in relation to CRC-specific mortality (corrected Pinteraction  = .0004); women had higher CRC-specific mortality associated with the minor allele (HR = 1.11, 95% CI = 1.04-1.19) whereas an inverse association was observed for men (HR = 0.88, 95% CI = 0.82-0.94). There was no evidence of interactions between CRP-associated SNPs and alcohol, obesity or smoking. Our study observed a significant interaction between sex and a CRP-associated variant in relation to CRC-specific mortality. Future replication of this association and functional annotation of the variant are needed.

Risk Stratification for Early-Onset Colorectal Cancer Using a Combination of Genetic and Environmental Risk Scores: An International Multi-Center Study

BACKGROUND

The incidence of colorectal cancer (CRC) among individuals aged younger than 50 years has been increasing. As screening guidelines lower the recommended age of screening initiation, concerns including the burden on screening capacity and costs have been recognized, suggesting that an individualized approach may be warranted. We developed risk prediction models for early-onset CRC that incorporate an environmental risk score (ERS), including 16 lifestyle and environmental factors, and a polygenic risk score (PRS) of 141 variants.

METHODS

Relying on risk score weights for ERS and PRS derived from studies of CRC at all ages, we evaluated risks for early-onset CRC in 3486 cases and 3890 controls aged younger than 50 years. Relative and absolute risks for early-onset CRC were assessed according to values of the ERS and PRS. The discriminatory performance of these scores was estimated using the covariate-adjusted area under the receiver operating characteristic curve.

RESULTS

Increasing values of ERS and PRS were associated with increasing relative risks for early-onset CRC (odds ratio per SD of ERS = 1.14, 95% confidence interval [CI] = 1.08 to 1.20; odds ratio per SD of PRS = 1.59, 95% CI = 1.51 to 1.68), both contributing to case-control discrimination (area under the curve = 0.631, 95% CI = 0.615 to 0.647). Based on absolute risks, we can expect 26 excess cases per 10 000 men and 21 per 10 000 women among those scoring at the 90th percentile for both risk scores.

CONCLUSIONS

Personal risk scores have the potential to identify individuals at differential relative and absolute risk for early-onset CRC. Improved discrimination may aid in targeted CRC screening of younger, high-risk individuals, potentially improving outcomes.

Predicted gene expression in ancestrally diverse populations leads to discovery of susceptibility loci for lifestyle and cardiometabolic traits

One mechanism by which genetic factors influence complex traits and diseases is altering gene expression. Direct measurement of gene expression in relevant tissues is rarely tenable; however, genetically regulated gene expression (GReX) can be estimated using prediction models derived from large multi-omic datasets. These approaches have led to the discovery of many gene-trait associations, but whether models derived from predominantly European ancestry (EA) reference panels can map novel associations in ancestrally diverse populations remains unclear. We applied PrediXcan to impute GReX in 51,520 ancestrally diverse Population Architecture using Genomics and Epidemiology (PAGE) participants (35% African American, 45% Hispanic/Latino, 10% Asian, and 7% Hawaiian) across 25 key cardiometabolic traits and relevant tissues to identify 102 novel associations. We then compared associations in PAGE to those in a random subset of 50,000 White British participants from UK Biobank (UKBB50k) for height and body mass index (BMI). We identified 517 associations across 47 tissues in PAGE but not UKBB50k, demonstrating the importance of diverse samples in identifying trait-associated GReX. We observed that variants used in PrediXcan models were either more or less differentiated across continental-level populations than matched-control variants depending on the specific population reflecting sampling bias. Additionally, variants from identified genes specific to either PAGE or UKBB50k analyses were more ancestrally differentiated than those in genes detected in both analyses, underlining the value of population-specific discoveries. This suggests that while EA-derived transcriptome imputation models can identify new associations in non-EA populations, models derived from closely matched reference panels may yield further insights. Our findings call for more diversity in reference datasets of tissue-specific gene expression.

Antibacterial properties of functionalized silk fibroin and sericin membranes for wound healing applications in oral and maxillofacial surgery

Oral wounds are among the most troublesome injuries which easily affect the patients' quality of life. To date, the development of functional antibacterial dressings for oral wound healing remains a challenge. In this regard, we investigated antibacterial silk protein-based membranes for the application as wound dressings in oral and maxillofacial surgery. The present study includes five variants of casted membranes, i.e., i) membranes-silver nanoparticles (CM-Ag), ii) membranes-gentamicin (CM-G), iii) membranes-control (without functionalization) (CM-C), iv) membranes-silk sericin control (CM-SSC), and v) membranes-silk fibroin/silk sericin (CM-SF/SS), and three variants of nonwovens, i.e., i) silver nanoparticles (NW-Ag), ii) gentamicin (NW-G), iii) control (without functionalization) (NW-C). The surface structure of the samples was visualized with scanning electron microscopy. In addition, antibacterial testing was accomplished using agar diffusion assay, colony forming unit (CFU) analysis, and qrt-PCR. Following antibacterial assays, biocompatibility was evaluated by cell proliferation assay (XTT), cytotoxicity assay (LDH), and live-dead assay on L929 mouse fibroblasts. Findings indicated significantly lower bacterial colony growth and DNA counts for CM-Ag with a reduction of bacterial counts by 3log levels (99.9% reduction) in CFU and qrt-PCR assay compared to untreated control membranes (CM-C and CM-SSC) and membranes functionalized with gentamicin (CM-G and NW-G) (p < 0.001). Similarly, NW-G yielded significantly lower DNA and colony growth counts compared to NW-Ag and NW-C (p < 0.001). In conclusion, CM-Ag represented 1log level better antibacterial activity compared to NW-G, whereas NW-G showed better cytocompatibility for L929 cells. As data suggest, these two membranes have the potential of application in the field of bacteria-free oral wound healing. However, provided that loading strategy and cytocompatibility are adjusted according to the antibacterial agents' characteristic and fabrication technique of the membranes.

Sleep problems and risk of cancer incidence and mortality in an older cohort: The Cardiovascular Health Study (CHS)

BACKGROUND

Sleep problems (SP) can indicate underlying sleep disorders, such as obstructive sleep apnea, which may adversely impact cancer risk and mortality.

METHODS

We assessed the association of baseline and longitudinal sleep apnea and insomnia symptoms with incident cancer (N = 3930) and cancer mortality (N = 4580) in the Cardiovascular Health Study. We used Cox proportional hazards regression to calculate adjusted hazard ratios (HR) and 95% confidence intervals (CI) to evaluate the associations.

RESULTS

Overall, 885 incident cancers and 804 cancer deaths were identified over a median follow-up of 12 and 14 years, respectively. Compared to participants who reported no sleep apnea symptoms, the risk of incident cancer was inversely associated [(HR (95%CI)] with snoring [0.84 (0.71, 0.99)]. We noted an elevated prostate cancer incidence for apnea [2.34 (1.32, 4.15)] and snoring [1.69 (1.11, 2.57)]. We also noted an elevated HR for lymphatic or hematopoietic cancers [daytime sleepiness: 1.81 (1.06, 3.08)]. We found an inverse relationship for cancer mortality with respect to snoring [0.73 (0.62, 0.8)] and apnea [(0.69 (0.51, 0.94))]. We noted a significant inverse relationship between difficulty falling asleep and colorectal cancer death [0.32 (0.15, 0.69)] and snoring with lung cancer death [0.56 (0.35, 0.89)].

CONCLUSIONS

The relationship between SP and cancer risk and mortality was heterogeneous. Larger prospective studies addressing more cancer sites, molecular type-specific associations, and better longitudinal SP assessments are needed for improved delineation of SP-cancer risk dyad.

Genome-wide association study identifies tumor anatomical site-specific risk variants for colorectal cancer survival

Identification of new genetic markers may improve the prediction of colorectal cancer prognosis. Our objective was to examine genome-wide associations of germline genetic variants with disease-specific survival in an analysis of 16,964 cases of colorectal cancer. We analyzed genotype and colorectal cancer-specific survival data from a consortium of 15 studies. Approximately 7.5 million SNPs were examined under the log-additive model using Cox proportional hazards models, adjusting for clinical factors and principal components. Additionally, we ran secondary analyses stratifying by tumor site and disease stage. We used a genome-wide p-value threshold of 5 × 10-8 to assess statistical significance. No variants were statistically significantly associated with disease-specific survival in the full case analysis or in the stage-stratified analyses. Three SNPs were statistically significantly associated with disease-specific survival for cases with tumors located in the distal colon (rs698022, HR = 1.48, CI 1.30-1.69, p = 8.47 × 10-9) and the proximal colon (rs189655236, HR = 2.14, 95% CI 1.65-2.77, p = 9.19 × 10-9 and rs144717887, HR = 2.01, 95% CI 1.57-2.58, p = 3.14 × 10-8), whereas no associations were detected for rectal tumors. Findings from this large genome-wide association study highlight the potential for anatomical-site-stratified genome-wide studies to identify germline genetic risk variants associated with colorectal cancer-specific survival. Larger sample sizes and further replication efforts are needed to more fully interpret these findings.

Genome-wide and transcriptome-wide association studies of mammographic density phenotypes reveal novel loci

BACKGROUND

Mammographic density (MD) phenotypes, including percent density (PMD), area of dense tissue (DA), and area of non-dense tissue (NDA), are associated with breast cancer risk. Twin studies suggest that MD phenotypes are highly heritable. However, only a small proportion of their variance is explained by identified genetic variants.

METHODS

We conducted a genome-wide association study, as well as a transcriptome-wide association study (TWAS), of age- and BMI-adjusted DA, NDA, and PMD in up to 27,900 European-ancestry women from the MODE/BCAC consortia.

RESULTS

We identified 28 genome-wide significant loci for MD phenotypes, including nine novel signals (5q11.2, 5q14.1, 5q31.1, 5q33.3, 5q35.1, 7p11.2, 8q24.13, 12p11.2, 16q12.2). Further, 45% of all known breast cancer SNPs were associated with at least one MD phenotype at p < 0.05. TWAS further identified two novel genes (SHOX2 and CRISPLD2) whose genetically predicted expression was significantly associated with MD phenotypes.

CONCLUSIONS

Our findings provided novel insight into the genetic background of MD phenotypes, and further demonstrated their shared genetic basis with breast cancer.

Molecular and Pathology Features of Colorectal Tumors and Patient Outcomes Are Associated with Fusobacterium nucleatum and Its Subspecies animalis

BACKGROUND

Fusobacterium nucleatum (F. nucleatum) activates oncogenic signaling pathways and induces inflammation to promote colorectal carcinogenesis.

METHODS

We characterized F. nucleatum and its subspecies in colorectal tumors and examined associations with tumor characteristics and colorectal cancer-specific survival. We conducted deep sequencing of nusA, nusG, and bacterial 16s rRNA genes in tumors from 1,994 patients with colorectal cancer and assessed associations between F. nucleatum presence and clinical characteristics, colorectal cancer-specific mortality, and somatic mutations.

RESULTS

F. nucleatum, which was present in 10.3% of tumors, was detected in a higher proportion of right-sided and advanced-stage tumors, particularly subspecies animalis. Presence of F. nucleatum was associated with higher colorectal cancer-specific mortality (HR, 1.97; P = 0.0004). This association was restricted to nonhypermutated, microsatellite-stable tumors (HR, 2.13; P = 0.0002) and those who received chemotherapy [HR, 1.92; confidence interval (CI), 1.07-3.45; P = 0.029). Only F. nucleatum subspecies animalis, the main subspecies detected (65.8%), was associated with colorectal cancer-specific mortality (HR, 2.16; P = 0.0016), subspecies vincentii and nucleatum were not (HR, 1.07; P = 0.86). Additional adjustment for tumor stage suggests that the effect of F. nucleatum on mortality is partly driven by a stage shift. Presence of F. nucleatum was associated with microsatellite instable tumors, tumors with POLE exonuclease domain mutations, and ERBB3 mutations, and suggestively associated with TP53 mutations.

CONCLUSIONS

F. nucleatum, and particularly subspecies animalis, was associated with a higher colorectal cancer-specific mortality and specific somatic mutated genes.

IMPACT

Our findings identify the F. nucleatum subspecies animalis as negatively impacting colorectal cancer mortality, which may occur through a stage shift and its effect on chemoresistance.

A genome-wide search for determinants of survival in 1926 patients with advanced colorectal cancer with follow-up in over 22,000 patients

BACKGROUND

While genome-wide association studies (GWAS) have identified germline variants influencing the risk of developing colorectal cancer (CRC), there has been limited examination of the possible role of inherited variation as a determinant of patient outcome.

PATIENTS AND METHODS

We performed a GWAS for overall survival (OS) in 1926 patients with advanced CRC from the COIN and COIN-B clinical trials. For single nucleotide polymorphisms (SNPs) showing an association with OS (P < 1.0 × 10-5), we conducted sensitivity analyses based on the time from diagnosis to death and sought independent replications in 5675 patients from the Study of Colorectal Cancer in Scotland (SOCCS) and 16,964 patients from the International Survival Analysis in Colorectal cancer Consortium (ISACC). We analysed the Human Protein Atlas to determine if ERBB4 expression was associated with survival in 438 patients with colon adenocarcinomas.

RESULTS

The most significant SNP associated with OS was rs79612564 in ERBB4 (hazard ratio [HR] = 1.24, 95% confidence interval [CI] = 1.16-1.32, P = 1.9 × 10-7). SNPs at 17 loci had suggestive associations for OS and all had similar effects on the time from diagnosis to death. No lead SNPs were independently replicated in the meta-analysis of all patients from SOCCS and ISACC. However, rs79612564 was significant in stage-IV patients from SOCCS (P = 2.1 × 10-2) but not ISACC (P = 0.89) and SOCCS combined with COIN and COIN-B attained genome-wide significance (P = 1.7 × 10-8). Patients with high ERBB4 expression in their colon adenocarcinomas had worse survival (HR = 1.50, 95% CI = 1.1-1.9, P = 4.6 × 10-2).

CONCLUSIONS

Genetic and expression data support a potential role for rs79612564 in the receptor tyrosine kinase ERBB4 as a predictive biomarker of survival.

A Statistical Method for Association Analysis of Cell Type Compositions

Gene expression data are often collected from tissue samples that are composed of multiple cell types. Studies of cell type composition based on gene expression data from tissue samples have recently attracted increasing research interest and led to new method development for cell type composition estimation. This new information on cell type composition can be associated with individual characteristics (e.g., genetic variants) or clinical outcomes (e.g., survival time). Such association analysis can be conducted for each cell type separately followed by multiple testing correction. An alternative approach is to evaluate this association using the composition of all the cell types, thus aggregating association signals across cell types. A key challenge of this approach is to account for the dependence across cell types. We propose a new method to quantify the distances between cell types while accounting for their dependencies, and use this information for association analysis. We demonstrate our method in two applied examples: to assess the association between immune cell type composition in tumor samples of colorectal cancer patients versus survival time and SNP genotypes. We found immune cell composition has prognostic value, and our distance metric leads to more accurate survival time prediction than other distance metrics that ignore cell type dependencies. In addition, survival time-associated SNPs are enriched among the SNPs associated with immune cell composition.

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

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

Salicylic Acid and Risk of Colorectal Cancer: A Two-Sample Mendelian Randomization Study

Salicylic acid (SA) has observationally been shown to decrease colorectal cancer (CRC) risk. Aspirin (acetylsalicylic acid, that rapidly deacetylates to SA) is an effective primary and secondary chemopreventive agent. Through a Mendelian randomization (MR) approach, we aimed to address whether levels of SA affected CRC risk, stratifying by aspirin use. A two-sample MR analysis was performed using GWAS summary statistics of SA (INTERVAL and EPIC-Norfolk, N = 14,149) and CRC (CCFR, CORECT, GECCO and UK Biobank, 55,168 cases and 65,160 controls). The DACHS study (4410 cases and 3441 controls) was used for replication and stratification of aspirin-use. SNPs proxying SA were selected via three methods: (1) functional SNPs that influence the activity of aspirin-metabolising enzymes; (2) pathway SNPs present in enzymes' coding regions; and (3) genome-wide significant SNPs. We found no association between functional SNPs and SA levels. The pathway and genome-wide SNPs showed no association between SA and CRC risk (OR: 1.03, 95% CI: 0.84-1.27 and OR: 1.08, 95% CI: 0.86-1.34, respectively). Results remained unchanged upon aspirin use stratification. We found little evidence to suggest that an SD increase in genetically predicted SA protects against CRC risk in the general population and upon stratification by aspirin use.

Smoking Behavior and Prognosis After Colorectal Cancer Diagnosis: A Pooled Analysis of 11 Studies

Background

Smoking has been associated with colorectal cancer (CRC) incidence and mortality in previous studies, but current evidence on smoking in association with survival after CRC diagnosis is limited.

Methods

We pooled data from 12 345 patients with stage I-IV CRC from 11 epidemiologic studies in the International Survival Analysis in Colorectal Cancer Consortium. Cox proportional hazards regression models were used to evaluate the associations of prediagnostic smoking behavior with overall, CRC-specific, and non-CRC-specific survival.

Results

Among 12 345 patients with CRC, 4379 (35.5%) died (2515 from CRC) over a median follow-up time of 7.5 years. Smoking was strongly associated with worse survival in stage I-III patients, whereas no association was observed among stage IV patients. Among stage I-III patients, clear dose-response relationships with all survival outcomes were seen for current smokers. For example, current smokers with 40 or more pack-years had statistically significantly worse overall, CRC-specific, and non-CRC-specific survival compared with never smokers (hazard ratio [HR] =1.94, 95% confidence interval [CI] =1.68 to 2.25; HR = 1.41, 95% CI = 1.12 to 1.78; and HR = 2.67, 95% CI = 2.19 to 3.26, respectively). Similar associations with all survival outcomes were observed for former smokers who had quit for less than 10 years, but only a weak association with non-CRC-specific survival was seen among former smokers who had quit for more than 10 years.

Conclusions

This large consortium of CRC patient studies provides compelling evidence that smoking is strongly associated with worse survival of stage I-III CRC patients in a clear dose-response manner. The detrimental effect of smoking was primarily related to noncolorectal cancer events, but current heavy smoking also showed an association with CRC-specific survival.

Hepcidin-regulating iron metabolism genes and pancreatic ductal adenocarcinoma: a pathway analysis of genome-wide association studies

BACKGROUND

Epidemiological studies have suggested positive associations for iron and red meat intake with risk of pancreatic ductal adenocarcinoma (PDAC). Inherited pathogenic variants in genes involved in the hepcidin-regulating iron metabolism pathway are known to cause iron overload and hemochromatosis.

OBJECTIVES

The objective of this study was to determine whether common genetic variation in the hepcidin-regulating iron metabolism pathway is associated with PDAC.

METHODS

We conducted a pathway analysis of the hepcidin-regulating genes using single nucleotide polymorphism (SNP) summary statistics generated from 4 genome-wide association studies in 2 large consortium studies using the summary data-based adaptive rank truncated product method. Our population consisted of 9253 PDAC cases and 12,525 controls of European descent. Our analysis included 11 hepcidin-regulating genes [bone morphogenetic protein 2 (BMP2), bone morphogenetic protein 6 (BMP6), ferritin heavy chain 1 (FTH1), ferritin light chain (FTL), hepcidin (HAMP), homeostatic iron regulator (HFE), hemojuvelin (HJV), nuclear factor erythroid 2-related factor 2 (NRF2), ferroportin 1 (SLC40A1), transferrin receptor 1 (TFR1), and transferrin receptor 2 (TFR2)] and their surrounding genomic regions (±20 kb) for a total of 412 SNPs.

RESULTS

The hepcidin-regulating gene pathway was significantly associated with PDAC (P = 0.002), with the HJV, TFR2, TFR1, BMP6, and HAMP genes contributing the most to the association.

CONCLUSIONS

Our results support that genetic susceptibility related to the hepcidin-regulating gene pathway is associated with PDAC risk and suggest a potential role of iron metabolism in pancreatic carcinogenesis. Further studies are needed to evaluate effect modification by intake of iron-rich foods on this association.

Bacterial translocation and microgap formation at a novel conical indexed implant abutment system for single crowns

OBJECTIVES

A conometric concept was recently introduced in which conical implant abutments hold the matching crown copings by friction alone, eliminating the need for cement or screws. The aim of this in vitro study was to assess the presence of microgap formation and bacterial leakage at the Acuris conometric restorative interface of three different implant abutment systems.

MATERIAL AND METHODS

A total of 75 Acuris samples of three implant-abutment systems (Ankylos, Astra Tech EV, Xive) were subjected to microbiological (n = 60) and scanning electron microscopic (SEM) investigation (n = 15). Bacterial migration into and out of the conical coupling system were analyzed in an anaerobic workstation for 48, 96, 144, and 192 h. Bacterial DNA quantification using qrt-PCR was performed at each time point. The precision of the conometric coupling and internal fit of cemented CAD/CAM crowns on corresponding Acuris TiN copings were determined by means of SEM.

RESULTS

qrt-PCR results failed to demonstrate microbial leakage from or into the Acuris system. SEM analysis revealed minute punctate microgaps at the apical aspect of the conometric junction (2.04 to 2.64 µm), while mean cement gaps of 12 to 145 µm were observed at the crown-coping interface.

CONCLUSIONS

The prosthetic morse taper connection of all systems examined does not allow bacterial passage. Marginal integrity and internal luting gap between the ceramic crown and the coping remained within the clinically acceptable limits.

CLINICAL RELEVANCE

Conometrically seated single crowns provide sufficient sealing efficiency, relocating potential misfits from the crown-abutment interface to the crown-coping interface.

Association Between Smoking and Molecular Subtypes of Colorectal Cancer

Background

Smoking is associated with colorectal cancer (CRC) risk. Previous studies suggested this association may be restricted to certain molecular subtypes of CRC, but large-scale comprehensive analysis is lacking.

Methods

A total of 9789 CRC cases and 11 231 controls of European ancestry from 11 observational studies were included. We harmonized smoking variables across studies and derived sex study-specific quartiles of pack-years of smoking for analysis. Four somatic colorectal tumor markers were assessed individually and in combination, including BRAF mutation, KRAS mutation, CpG island methylator phenotype (CIMP), and microsatellite instability (MSI) status. A multinomial logistic regression analysis was used to assess the association between smoking and risk of CRC subtypes by molecular characteristics, adjusting for age, sex, and study. All statistical tests were 2-sided and adjusted for Bonferroni correction.

Results

Heavier smoking was associated with higher risk of CRC overall and stratified by individual markers (P trend < .001). The associations differed statistically significantly between all molecular subtypes, which was the most statistically significant for CIMP and BRAF. Compared with never-smokers, smokers in the fourth quartile of pack-years had a 90% higher risk of CIMP-positive CRC (odds ratio = 1.90, 95% confidence interval = 1.60 to 2.26) but only 35% higher risk for CIMP-negative CRC (odds ratio = 1.35, 95% confidence interval = 1.22 to 1.49; P difference = 2.1 x 10-6). The association was also stronger in tumors that were CIMP positive, MSI high, or KRAS wild type when combined (P difference < .001).

Conclusion

Smoking was associated with differential risk of CRC subtypes defined by molecular characteristics. Heavier smokers had particularly higher risk of CRC subtypes that were CIMP positive and MSI high in combination, suggesting that smoking may be involved in the development of colorectal tumors via the serrated pathway.

Genetic insights into biological mechanisms governing human ovarian ageing

Reproductive longevity is essential for fertility and influences healthy ageing in women1,2, but insights into its underlying biological mechanisms and treatments to preserve it are limited. Here we identify 290 genetic determinants of ovarian ageing, assessed using normal variation in age at natural menopause (ANM) in about 200,000 women of European ancestry. These common alleles were associated with clinical extremes of ANM; women in the top 1% of genetic susceptibility have an equivalent risk of premature ovarian insufficiency to those carrying monogenic FMR1 premutations3. The identified loci implicate a broad range of DNA damage response (DDR) processes and include loss-of-function variants in key DDR-associated genes. Integration with experimental models demonstrates that these DDR processes act across the life-course to shape the ovarian reserve and its rate of depletion. Furthermore, we demonstrate that experimental manipulation of DDR pathways highlighted by human genetics increases fertility and extends reproductive life in mice. Causal inference analyses using the identified genetic variants indicate that extending reproductive life in women improves bone health and reduces risk of type 2 diabetes, but increases the risk of hormone-sensitive cancers. These findings provide insight into the mechanisms that govern ovarian ageing, when they act, and how they might be targeted by therapeutic approaches to extend fertility and prevent disease.

Abstract LB090: Associations of somatically mutated genes and pathways with colorectal cancer specific survival in 4,500 colorectal cancer patients

Colorectal cancer (CRC) is a heterogenous disease that develops through somatic mutations in driver genes, leading to activation of diverse neoplastic pathways. To systematically examine if somatically mutated genes and pathways impact survival, we sequenced tumor and normal DNA samples for 4,512 CRC cases using a targeted panel. We performed Cox regression to estimate hazard ratios (HR) and 95% confidence intervals (CIs) for associations of disease-specific (DS)-survival with somatically mutated genes and pathways, adjusting for age and sex, and stratifying baseline hazards by study population. We assessed statistical significance using Bonferroni p-value thresholds to account for multiple testing of 214 genes (2.34x10-4) and 6 key CRC pathways (8.3x10-3). We limited analyses to non-silent mutations.

We observed that DS-survival was significantly more favorable among individuals with hypermutated (HM) tumors, primarily consisting of microsatellite unstable and POLE-mutated tumors, as compared to those with non-hypermutated (NHM) tumors (HR=0.4, 95% CI: 0.3-0.5, p=1.2x10-17). BRAF V600E mutations were associated with poorer survival (HR=2.0, 95% CI: 1.6-2.5, p=1.6x10-10). This association was more pronounced among NHM tumors (HR=2.3, 95% CI: 1.8-2.9, p=4.3x10-12). We identified suggestive associations (p-values &lt; 5.0x10-3 in overall or stratified analyses) between DS-survival and mutations in B2M, TP53, and SMAD4. Mutations in B2M may provide more favorable prognosis for survival (HR=0.5, 95% CI: 0.3-0.8, p=4.4x10-3), with similar effect sizes in HM and NHM tumors. Poorer survival may be associated with mutations in TP53 (HR=1.2, 95% CI:1.2-1.4, p=9.3x10-4) and SMAD4 (HR=1.3, 95% CI:1.1-1.6, p=3.2x10-3) in NHM tumors.

We further observed statistically significant associations between survival and three mutated pathways: TP53/ATM (HR=1.2, 95% CI:1.1-1.4, p=9.0x10-4), receptor tyrosine kinases (RTK) and RAS (HR=1.3, 95% CI:1.1-1.5, p=4.5x10-5), and TGF-beta (HR=1.3, 95% CI:1.1-1.5, p=2.9x10-4). Findings for TP53/ATM and RTK/RAS were primarily due to mutations in one gene within each pathway (TP53 and BRAF, respectively). However, the TGF-beta pathway finding, which was more pronounced in NHM tumors, included two genes with p-values less than 0.05 (SMAD4 p=3.2x10-3 and TGFBR2 p=0.01).

Despite our large sample size, relatively few somatically mutated genes were significantly associated with DS-survival. It may be that pathway-level testing affords more power for analyses and that larger sample sizes are needed. It is of interest that non-silent mutations in B2M, whose product is a component of the class I major histocompatibility complex, was associated with improved DS-survival, though our finding was not statistically significant after multiple testing correction so further exploration and replication is needed.

Citation Format: Tabitha A. Harrison, Syed H. Zaidi, Conghui Qu, Amanda I. Phipps, Robert S. Steinfelder, Quang M. Trinh, Sonja I. Berndt, Daniel D. Buchanan, Peter T. Campbell, Andrew T. Chan, Kim F. Doheny, David A. Drew, Jane C. Figueiredo, Steven J. Gallinger, Marios Giannakis, Andrea Gsur, Marc J. Gunter, Michael Hoffmeister, Wen-Yi Huang, Paul J. Limburg, Victor Moreno, Polly A. Newcomb, Shuji Ogino, Ross L. Prentice, Tameka Shelford, Wei Sun, Stephen N. Thibodeau, Li Hsu, Ulrike Peters. Associations of somatically mutated genes and pathways with colorectal cancer specific survival in 4,500 colorectal cancer patients [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr LB090.

Abstract 881: Benchmarking genome-wide polygenic risk score development techniques in colorectal cancer risk prediction

Colorectal cancer (CRC) is a leading cause of cancer death, yet many CRC deaths are preventable via CRC screening. Currently only age and family history are used to define screening eligibility. However, CRC risk varies substantially in the population. In recent years polygenic risk scores (PRS) have gained attention as powerful risk prediction tool to personalize interventions. PRS provides a quantitative measure of an individual's inherited risk based on the cumulative effect of many genetic risk variants. Here, we benchmark several genome wide PRS techniques to select the best performing models in CRC risk prediction.

We built CRC risk prediction models that incorporate genome-wide genotype data from large-scale research studies (55,105 cases and 65,079 controls, European ancestries) with the imputed genetic data on over 40 million variants. The risk prediction models were externally evaluated in the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort, including 101,987 genotyped individuals within the Kaiser Permanente Northern California (KPNC) integrated healthcare delivery system. We built genome-wide PRS using various methods including known CRC risk variants, thresholding and pruning followed by machine learning approaches (ML), LDpred, improved LDpred2, SBayesR, PRS-CS, Lassosum and empirical Bayes.

Among 55,033 individuals of European ancestry in the GERA cohort, we evaluated the performance of models in terms of the age and sex-adjusted AUC. We showed that LDpred, LDpred2, LDpred2-sparse, SBayesR and PRS-CS perform equally well in terms of discriminatory accuracy (AUC=0.65). In addition, the PRS developed using the above-mentioned techniques identified the top 30% of the GERA European population has a hazard ratio estimate of ~2.2 on CRC risk, which is comparable to that for having an affected first-degree relative. The developed CRC PRSs will provide way for risk-stratified CRC screening and other targeted interventions.

PRS derivation methodsNo. of variantsAUC(1,311 cases and 53,722 controls)Hazard ratio estimates (CI)Top 30% of population vs. remainingKnown variants1400.631.92 (1.75-2.23)PT Clumping + ML (Ridge)10,0000.631.94 (1.72-2.19)LDpred1.2M0.652.20 (1.94-2.47)LDPred21.2M0.652.20 (1.93-2.45)LDpred2 Sparse530K0.652.20 (1.90-2.41)SBayesR1.2M0.652.20 (1.88-2.38)PRS-CS1.2M0.652.20 (1.91-2.43)Lassosum1.2M0.621.76 (1.56-2.58)EBPRS1.2M0.621.81 (1.66-2.11)AUC based on family history in GERA cohort is 0.54

Citation Format: Minta Thomas, Lori C Sakoda, Jeffrey K Lee, Mark A Jenkins, Andrea Burnett-Hartman, Heather Hampel, Elisabeth A Rosenthal, Hermann Brenner, Jenny Chang-Claude, Marc J Gunter, Polly A Newcomb, Steven Gallinger, Tabitha A Harrison, Graham Casey, Victor Moreno, Gail P Jarvik, Stephen B Gruber, Robert E Schoen, Andrew T Chan, Richard B Hayes, Douglas A Corley, Ulrike Peters, Li Hsu. Benchmarking genome-wide polygenic risk score development techniques in colorectal cancer risk prediction [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 881.

Oral Microbiota Transplant in Dogs with Naturally Occurring Periodontitis

In periodontitis patients, dysbiosis of the oral microbiota is not only found at clinically diseased periodontal sites but also at clinically healthy periodontal sites, buccal mucosae, tongue, and saliva. The present study evaluated the safety and efficacy of an oral microbiota transplant (OMT) for the treatment of periodontitis in dogs. Eighteen systemically healthy beagle dogs with naturally occurring periodontitis were enrolled in the study and randomly assigned to a test or control group. A 4-y-old, periodontally healthy female beagle dog served as a universal OMT donor. To reduce periodontal inflammation, all dogs received full-mouth mechanical debridement of teeth and mucosae 2 wk before baseline. At baseline, full-mouth mechanical debridement was repeated and followed by adjunctive subgingival and oral irrigation with 0.1% NaOCl. Subsequently, test dogs were inoculated with an OMT from the healthy donor. No daily oral hygiene was performed after OMT transplantation. Adverse events were assessed throughout the observation period. Clinical examinations were performed and whole-mouth oral microbiota samples were collected at week 2, baseline, week 2, and week 12. The composition of oral microbiota samples was analyzed using high-throughput 16S ribosomal RNA gene amplicon sequencing followed by taxonomic assignment and downstream bioinformatic and statistical analyses. Results demonstrated that the intergroup difference in the primary outcome measure, probing pocket depth at week 12, was statistically insignificant. However, the single adjunctive OMT had an additional effect on the oral microbiota composition compared to the full-mouth mechanical and antimicrobial debridement alone. The OMT resulted in an "ecological shift" toward the composition of the donor microbiota, but this was transient in nature and was not observed at week 12. No local or systemic adverse events were observed throughout the study period. The results indicate that OMT may modulate the microbiota composition in dogs with naturally occurring periodontitis and can be applied safely.