Defining novel causal SNPs and linked phenotypes at melanoma-associated loci.

A number of genomic regions have been associated with melanoma risk through genome-wide association studies; however, the causal variants underlying the majority of these associations remain unknown. Here, we sequenced either the full locus or the functional regions including exons of 19 melanoma-associated loci in 1959 British melanoma cases and 737 controls. Variant filtering followed by Fisher's exact test analyses identified 66 variants associated with melanoma risk. Sequential conditional logistic regression identified the distinct haplotypes on which variants reside, and massively parallel reporter assays provided biological insights into how these variants influence gene function. We performed further analyses to link variants to melanoma risk phenotypes and assessed their association with melanoma-specific survival. Our analyses replicate previously known associations in the melanocortin 1 receptor (MC1R) and tyrosinase (TYR) loci, while identifying novel potentially causal variants at the MTAP/CDKN2A and CASP8 loci. These results improve our understanding of the architecture of melanoma risk and outcome.

Mutually exclusive genetic interactions and gene essentiality shape the genomic landscape of primary melanoma.

Melanoma is a heterogenous malignancy with an unpredictable clinical course. Most patients who present in the clinic are diagnosed with primary melanoma, yet large-scale sequencing efforts have focused primarily on metastatic disease. In this study we sequence-profiled 524 American Joint Committee on Cancer Stage I-III primary tumours. Our analysis of these data reveals recurrent driver mutations, mutually exclusive genetic interactions, where two genes were never or rarely co-mutated, and an absence of co-occurring genetic events. Further, we intersected copy number calls from our primary melanoma data with whole-genome CRISPR screening data to identify the transcription factor interferon regulatory factor 4 (IRF4) as a melanoma-associated dependency. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Population-based analysis of POT1 variants in a cutaneous melanoma case-control cohort.

Pathogenic germline variants in the protection of telomeres 1 gene (POT1) have been associated with predisposition to a range of tumour types, including melanoma, glioma, leukaemia and cardiac angiosarcoma. We sequenced all coding exons of the POT1 gene in 2928 European-descent melanoma cases and 3298 controls, identifying 43 protein-changing genetic variants. We performed POT1-telomere binding assays for all missense and stop-gained variants, finding nine variants that impair or disrupt protein-telomere complex formation, and we further define the role of variants in the regulation of telomere length and complex formation through molecular dynamics simulations. We determine that POT1 coding variants are a minor contributor to melanoma burden in the general population, with only about 0.5% of melanoma cases carrying germline pathogenic variants in this gene, but should be screened in individuals with a strong family history of melanoma and/or multiple malignancies.

Total serum N-glycans associate with response to immune checkpoint inhibition therapy and survival in patients with advanced melanoma.

BACKGROUND: Immune checkpoint inhibitors (ICIs) have revolutionized the treatment of melanoma and other cancers. However, no reliable biomarker of survival or response has entered the clinic to identify those patients with melanoma who are most likely to benefit from ICIs. Glycosylation affects proteins and lipids' structure and functions. Tumours are characterized by aberrant glycosylation which may contribute to their progression and hinder an effective antitumour immune response. METHODS: We aim at identifying novel glyco-markers of response and survival by leveraging the N-glycome of total serum proteins collected in 88 ICI-naive patients with advanced melanoma from two European countries. Samples were collected before and during ICI treatment. RESULTS: We observe that responders to ICIs present with a pre-treatment N-glycome profile significantly shifted towards higher abundancy of low-branched structures containing lower abundances of antennary fucose, and that this profile is positively associated with survival and a better predictor of response than clinical variables alone. CONCLUSION: While changes in serum protein glycosylation have been previously implicated in a pro-metastatic melanoma behaviour, we show here that they are also associated with response to ICI, opening new avenues for the stratification of patients and the design of adjunct therapies aiming at improving immune response.

Inherited duplications of PPP2R3B predispose to nevi and melanoma via a C21orf91-driven proliferative phenotype.

PURPOSE: Much of the heredity of melanoma remains unexplained. We sought predisposing germline copy-number variants using a rare disease approach. METHODS: Whole-genome copy-number findings in patients with melanoma predisposition syndrome congenital melanocytic nevus were extrapolated to a sporadic melanoma cohort. Functional effects of duplications in PPP2R3B were investigated using immunohistochemistry, transcriptomics, and stable inducible cellular models, themselves characterized using RNAseq, quantitative real-time polymerase chain reaction (qRT-PCR), reverse phase protein arrays, immunoblotting, RNA interference, immunocytochemistry, proliferation, and migration assays. RESULTS: We identify here a previously unreported genetic susceptibility to melanoma and melanocytic nevi, familial duplications of gene PPP2R3B. This encodes PR70, a regulatory unit of critical phosphatase PP2A. Duplications increase expression of PR70 in human nevus, and increased expression in melanoma tissue correlates with survival via a nonimmunological mechanism. PPP2R3B overexpression induces pigment cell switching toward proliferation and away from migration. Importantly, this is independent of the known microphthalmia-associated transcription factor (MITF)-controlled switch, instead driven by C21orf91. Finally, C21orf91 is demonstrated to be downstream of MITF as well as PR70. CONCLUSION: This work confirms the power of a rare disease approach, identifying a previously unreported copy-number change predisposing to melanocytic neoplasia, and discovers C21orf91 as a potentially targetable hub in the control of phenotype switching.

Melanoma Genomics.

The incidence of cutaneous melanoma continues to increase in pale skinned peoples in Europe and elsewhere. Epidemiological studies identified genetically determined phenotypes such as pale skin, freckles and red hair, and sunburn as risk factors for this cancer. The development of many melanocytic naevi is also genetically determined and a strong melanoma risk phenotype. Not surprisingly then, genome wide association studies have identified pigmentation genes as common risk genes, and to a lesser extent, genes associated with melanocytic naevi. More unexpectedly, genes associated with telomere length have also been identified as risk genes. Higher risk susceptibility genes have been identified, particularly CDKN2A as the most common cause, and very rarely genes such as CDK4, POT1, TERT and other genes in coding for proteins in the shelterin complex are found to be mutated. Familial melanoma genes are associated with an increased number of melanocytic naevi but not invariably and the atypical naevus phenotype is therefore an imperfect marker of gene carrier status. At a somatic level, the most common driver mutation is BRAF, second most common NRAS, third NF1 and increasing numbers of additional rarer mutations are being identified such as in TP53. It is of note that the BRAF and NRAS mutations are not C>T accepted as characteristic of ultraviolet light induced mutations.

A population-based analysis of germline BAP1 mutations in melanoma.

Germline mutation of the BRCA1 associated protein-1 (BAP1) gene has been linked to uveal melanoma, mesothelioma, meningioma, renal cell carcinoma and basal cell carcinoma. Germline variants have also been found in familial cutaneous melanoma pedigrees, but their contribution to sporadic melanoma has not been fully assessed. We sequenced BAP1 in 1,977 melanoma cases and 754 controls and used deubiquitinase assays, a pedigree analysis, and a histopathological review to assess the consequences of the mutations found. Sequencing revealed 30 BAP1 variants in total, of which 27 were rare (ExAc allele frequency <0.002). Of the 27 rare variants, 22 were present in cases (18 missense, one splice acceptor, one frameshift and two near splice regions) and five in controls (all missense). A missense change (S98R) in a case that completely abolished BAP1 deubiquitinase activity was identified. Analysis of cancers in the pedigree of the proband carrying the S98R variant and in two other pedigrees carrying clear loss-of-function alleles showed the presence of BAP1-associated cancers such as renal cell carcinoma, mesothelioma and meningioma, but not uveal melanoma. Two of these three probands carrying BAP1 loss-of-function variants also had melanomas with histopathological features suggestive of a germline BAP1 mutation. The remaining cases with germline mutations, which were predominantly missense mutations, were associated with less typical pedigrees and tumours lacking a characteristic BAP1-associated histopathological appearances, but may still represent less penetrant variants. Germline BAP1 alleles defined as loss-of-function or predicted to be deleterious/damaging are rare in cutaneous melanoma.

Estimating CDKN2A mutation carrier probability among global familial melanoma cases using GenoMELPREDICT.

BACKGROUND: Although rare in the general population, highly penetrant germline mutations in CDKN2A are responsible for 5%-40% of melanoma cases reported in melanoma-prone families. We sought to determine whether MELPREDICT was generalizable to a global series of families with melanoma and whether performance improvements can be achieved. METHODS: In total, 2116 familial melanoma cases were ascertained by the international GenoMEL Consortium. We recapitulated the MELPREDICT model within our data (GenoMELPREDICT) to assess performance improvements by adding phenotypic risk factors and history of pancreatic cancer. We report areas under the curve (AUC) with 95% confidence intervals (CIs) along with net reclassification indices (NRIs) as performance metrics. RESULTS: MELPREDICT performed well (AUC 0.752, 95% CI 0.730-0.775), and GenoMELPREDICT performance was similar (AUC 0.748, 95% CI 0.726-0.771). Adding a reported history of pancreatic cancer yielded discriminatory improvement (P < .0001) in GenoMELPREDICT (AUC 0.772, 95% CI 0.750-0.793, NRI 0.40). Including phenotypic risk factors did not improve performance. CONCLUSION: The MELPREDICT model functioned well in a global data set of familial melanoma cases. Adding pancreatic cancer history improved model prediction. GenoMELPREDICT is a simple tool for predicting CDKN2A mutational status among melanoma patients from melanoma-prone families and can aid in directing these patients to receive genetic testing or cancer risk counseling.

Identification of a novel susceptibility locus at 13q34 and refinement of the 20p12.2 region as a multi-signal locus associated with bladder cancer risk in individuals of European ancestry.

Candidate gene and genome-wide association studies (GWAS) have identified 15 independent genomic regions associated with bladder cancer risk. In search for additional susceptibility variants, we followed up on four promising single-nucleotide polymorphisms (SNPs) that had not achieved genome-wide significance in 6911 cases and 11 814 controls (rs6104690, rs4510656, rs5003154 and rs4907479, P < 1 × 10(-6)), using additional data from existing GWAS datasets and targeted genotyping for studies that did not have GWAS data. In a combined analysis, which included data on up to 15 058 cases and 286 270 controls, two SNPs achieved genome-wide statistical significance: rs6104690 in a gene desert at 20p12.2 (P = 2.19 × 10(-11)) and rs4907479 within the MCF2L gene at 13q34 (P = 3.3 × 10(-10)). Imputation and fine-mapping analyses were performed in these two regions for a subset of 5551 bladder cancer cases and 10 242 controls. Analyses at the 13q34 region suggest a single signal marked by rs4907479. In contrast, we detected two signals in the 20p12.2 region-the first signal is marked by rs6104690, and the second signal is marked by two moderately correlated SNPs (r(2) = 0.53), rs6108803 and the previously reported rs62185668. The second 20p12.2 signal is more strongly associated with the risk of muscle-invasive (T2-T4 stage) compared with non-muscle-invasive (Ta, T1 stage) bladder cancer (case-case P ≤ 0.02 for both rs62185668 and rs6108803). Functional analyses are needed to explore the biological mechanisms underlying these novel genetic associations with risk for bladder cancer.

Germline CDKN2A/P16INK4A mutations contribute to genetic determinism of sarcoma.

BACKGROUND: Sarcomas are rare mesenchymal malignancies whose pathogenesis is poorly understood; both environmental and genetic risk factors could contribute to their aetiology. METHODS AND RESULTS: We performed whole-exome sequencing (WES) in a familial aggregation of three individuals affected with soft-tissue sarcoma (STS) without TP53 mutation (Li-Fraumeni-like, LFL) and found a shared pathogenic mutation in CDKN2A tumour suppressor gene. We searched for individuals with sarcoma among 474 melanoma-prone families with a CDKN2A-/+ genotype and for CDKN2A mutations in 190 TP53-negative LFL families where the index case was a sarcoma. Including the initial family, eight independent sarcoma cases carried a germline mutation in the CDKN2A/p16INK4A gene. In five out of seven formalin-fixed paraffin-embedded sarcomas, heterozygosity was lost at germline CDKN2A mutations sites demonstrating complete loss of function. As sarcomas are rare in CDKN2A/p16INK4A carriers, we searched in constitutional WES of nine carriers for potential modifying rare variants and identified three in platelet-derived growth factor receptor (PDGFRA) gene. Molecular modelling showed that two never-described variants could impact the PDGFRA extracellular domain structure. CONCLUSION: Germline mutations in CDKN2A/P16INK4A, a gene known to predispose to hereditary melanoma, pancreatic cancer and tobacco-related cancers, account also for a subset of hereditary sarcoma. In addition, we identified PDGFRA as a candidate modifier gene.

A novel recurrent mutation in MITF predisposes to familial and sporadic melanoma.

So far, two genes associated with familial melanoma have been identified, accounting for a minority of genetic risk in families. Mutations in CDKN2A account for approximately 40% of familial cases, and predisposing mutations in CDK4 have been reported in a very small number of melanoma kindreds. Here we report the whole-genome sequencing of probands from several melanoma families, which we performed in order to identify other genes associated with familial melanoma. We identify one individual carrying a novel germline variant (coding DNA sequence c.G1075A; protein sequence p.E318K; rs149617956) in the melanoma-lineage-specific oncogene microphthalmia-associated transcription factor (MITF). Although the variant co-segregated with melanoma in some but not all cases in the family, linkage analysis of 31 families subsequently identified to carry the variant generated a log of odds (lod) score of 2.7 under a dominant model, indicating E318K as a possible intermediate risk variant. Consistent with this, the E318K variant was significantly associated with melanoma in a large Australian case-control sample. Likewise, it was similarly associated in an independent case-control sample from the United Kingdom. In the Australian sample, the variant allele was significantly over-represented in cases with a family history of melanoma, multiple primary melanomas, or both. The variant allele was also associated with increased naevus count and non-blue eye colour. Functional analysis of E318K showed that MITF encoded by the variant allele had impaired sumoylation and differentially regulated several MITF targets. These data indicate that MITF is a melanoma-predisposition gene and highlight the utility of whole-genome sequencing to identify novel rare variants associated with disease susceptibility.

Fine mapping of genetic susceptibility loci for melanoma reveals a mixture of single variant and multiple variant regions.

At least 17 genomic regions are established as harboring melanoma susceptibility variants, in most instances with genome-wide levels of significance and replication in independent samples. Based on genome-wide single nucleotide polymorphism (SNP) data augmented by imputation to the 1,000 Genomes reference panel, we have fine mapped these regions in over 5,000 individuals with melanoma (mainly from the GenoMEL consortium) and over 7,000 ethnically matched controls. A penalized regression approach was used to discover those SNP markers that most parsimoniously explain the observed association in each genomic region. For the majority of the regions, the signal is best explained by a single SNP, which sometimes, as in the tyrosinase region, is a known functional variant. However in five regions the explanation is more complex. At the CDKN2A locus, for example, there is strong evidence that not only multiple SNPs but also multiple genes are involved. Our results illustrate the variability in the biology underlying genome-wide susceptibility loci and make steps toward accounting for some of the "missing heritability."

Inherited variation in the PARP1 gene and survival from melanoma.

We report the association of an inherited variant located upstream of the poly(adenosine diphosphate-ribose) polymerase 1 (PARP1) gene (rs2249844), with survival in 11 BioGenoMEL melanoma cohorts. The gene encodes a protein involved in a number of cellular processes including single-strand DNA repair. Survival analysis was conducted for each cohort using proportional hazards regression adjusting for factors known to be associated with survival. Survival was measured as overall survival (OS) and, where available, melanoma-specific survival (MSS). Results were combined using random effects meta-analysis. Evidence for a role of the PARP1 protein in melanoma ulceration and survival was investigated by testing gene expression levels taken from formalin-fixed paraffin-embedded tumors. A significant association was seen for inheritance of the rarer variant of PARP1, rs2249844 with OS (hazard ratio (HR) = 1.16 per allele, 95% confidence interval (CI) 1.04-1.28, p = 0.005, eleven cohorts) and MSS (HR = 1.20 per allele, 95% CI 1.01-1.39, p = 0.03, eight cohorts). We report bioinformatic data supportive of a functional effect for rs2249844. Higher levels of PARP1 gene expression in tumors were shown to be associated with tumor ulceration and poorer OS.

Melanoma prone families with CDK4 germline mutation: phenotypic profile and associations with MC1R variants.

BACKGROUND: CDKN2A and CDK4 are high risk susceptibility genes for cutaneous malignant melanoma. Melanoma families with CDKN2A germline mutations have been extensively characterised, whereas CDK4 families are rare and lack a systematic investigation of their phenotype. METHODS: All known families with CDK4 germline mutations (n=17) were recruited for the study by contacting the authors of published papers or by requests via the Melanoma Genetics Consortium (GenoMEL). Phenotypic data related to primary melanoma and pigmentation characteristics were collected. The CDK4 exon 2 and the complete coding region of the MC1R gene were sequenced. RESULTS: Eleven families carried the CDK4 R24H mutation whereas six families had the R24C mutation. The total number of subjects with verified melanoma was 103, with a median age at first melanoma diagnosis of 39 years. Forty-three (41.7%) subjects had developed multiple primary melanomas (MPM). A CDK4 mutation was found in 89 (including 62 melanoma cases) of 209 tested subjects. CDK4 positive family members (both melanoma cases and unaffected subjects) were more likely to have clinically atypical nevi than CDK4 negative family members (p<0.001). MPM subjects had a higher frequency of MC1R red hair colour variants compared with subjects with one tumour (p=0.010). CONCLUSION: Our study shows that families with CDK4 germline mutations cannot be distinguished phenotypically from CDKN2A melanoma families, which are characterised by early onset of disease, increased occurrence of clinically atypical nevi, and development of MPM. In a clinical setting, the CDK4 gene should therefore always be examined when a melanoma family tests negative for CDKN2A mutation.

Prevalence and predictors of germline CDKN2A mutations for melanoma cases from Australia, Spain and the United Kingdom.

BACKGROUND: Mutations in the CDKN2A and CDK4 genes predispose to melanoma. From three case-control studies of cutaneous melanoma, we estimated the prevalence and predictors of these mutations for people from regions with widely differing latitudes and melanoma incidence. METHODS: Population-based cases and controls from the United Kingdom (1586 cases, 499 controls) and Australia (596 early-onset cases, 476 controls), and a hospital-based series from Spain (747 cases, 109 controls), were screened for variants in all exons of CDKN2A and the p16INK4A binding domain of CDK4. RESULTS: The prevalence of mutations for people with melanoma was similar across regions: 2.3%, 2.5% and 2.0% for Australia, Spain and the United Kingdom respectively. The strongest predictors of carrying a mutation were having multiple primaries (odds ratio (OR) = 5.4, 95% confidence interval (CI: 2.5, 11.6) for 2 primaries and OR = 32.4 (95% CI: 14.7, 71.2) for 3 or more compared with 1 primary only); and family history (OR = 3.8; 95% CI:1.89, 7.5) for 1 affected first- or second-degree relative and OR = 23.2 (95% CI: 11.3, 47.6) for 2 or more compared with no affected relatives). Only 1.1% of melanoma cases with neither a family history nor multiple primaries had mutations. CONCLUSIONS: There is a low probability (<2%) of detecting a germline CDKN2A mutation in people with melanoma except for those with a strong family history of melanoma (≥2 affected relatives, 25%), three or more primary melanomas (29%), or more than one primary melanoma who also have other affected relatives (27%).

Systemic Inflammation, the Peripheral Blood Transcriptome, and Primary Melanoma.

Peripheral blood transcriptomes from 383 patients with newly diagnosed melanoma were subjected to differential gene expression analysis. The hypotheses were that impaired systemic immunity is associated with poorer prognosis (thicker tumors and fewer tumor-infiltrating lymphocytes) and evidence of systemic inflammation (high-sensitivity CRP and fibrinogen levels). Higher fibrinogen levels were associated with thicker primary tumors. In single-gene analysis, high-sensitivity CRP levels were significantly associated with higher blood CD274 expression (coding for PD-L1), but each was independently prognostic, with high-sensitivity CRP associated with increased mortality and higher CD274 protective, independent of age. Pathway analysis identified downregulation of immune cell signaling pathways in the blood of people with thicker tumors and notable upregulation of signal transducer and activator of transcription 1 gene STAT1 in people with brisk tumor-infiltrating lymphocytes. Transcriptomic data provided evidence for increased NF-kB signaling with higher inflammatory markers but with reduction in expression of HLA class II molecules and higher CD274, suggesting that aberrant systemic inflammation is a significant mediator of reduced immune function in melanoma. In summary, transcriptomic data revealed evidence of reduced immune function in patients with thicker tumors and fewer tumor-infiltrating lymphocytes at diagnosis. Inflammatory markers were associated with thicker primaries and independently with death from melanoma, suggesting that systemic inflammation contributes to that reduced immune function.

APOE genotype, eicosapentaenoic acid (EPA) supplementation and n-3 highly unsaturated fatty acid (HUFA) levels in patients with multiple colorectal polyps: A secondary analysis of the seAFOod polyp prevention trial.

INTRODUCTION: We examined the relationship between Apolipoprotein E (APOE) genotype and n-3 highly unsaturated fatty acid (HUFA) levels in participants of the seAFOod trial, who were undergoing colonoscopy surveillance after removal of colorectal polyps. METHODS: Baseline and on-treatment (eicosapentaenoic acid [EPA] 2 g daily or placebo for 6 months) levels of n-3 HUFAs, and plasma 18-hydroxyeicosapentaenoic acid (HEPE), were analysed according to APOE genotype (based on polymorphisms rs429358 and rs7412) in 584 participants. RESULTS: Before treatment, APOE2/2 individuals had lower levels, and APOE4/4 participants had higher levels, of n-3 HUFAs, including EPA, than APOE3/3 counterparts (P < 0.01 for the APOE2/2 versus APOE4/4 comparison). After EPA supplementation, n-3 HUFA levels were not significantly different when stratified by APOE genotype, although APOE4 carriers displayed lower plasma 18-HEPE levels than individuals without an APOE4 allele (P = 0.002). CONCLUSIONS: APOE genotype is associated with differential n-3 HUFA and 18-HEPE levels in individuals with multiple colorectal polyps.

Longitudinal gut microbiome changes in immune checkpoint blockade-treated advanced melanoma.

Multiple clinical trials targeting the gut microbiome are being conducted to optimize treatment outcomes for immune checkpoint blockade (ICB). To improve the success of these interventions, understanding gut microbiome changes during ICB is urgently needed. Here through longitudinal microbiome profiling of 175 patients treated with ICB for advanced melanoma, we show that several microbial species-level genome bins (SGBs) and pathways exhibit distinct patterns from baseline in patients achieving progression-free survival (PFS) of 12 months or longer (PFS ≥12) versus patients with PFS shorter than 12 months (PFS <12). Out of 99 SGBs that could discriminate between these two groups, 20 were differentially abundant only at baseline, while 42 were differentially abundant only after treatment initiation. We identify five and four SGBs that had consistently higher abundances in patients with PFS ≥12 and <12 months, respectively. Constructing a log ratio of these SGBs, we find an association with overall survival. Finally, we find different microbial dynamics in different clinical contexts including the type of ICB regimen, development of immune-related adverse events and concomitant medication use. Insights into the longitudinal dynamics of the gut microbiome in association with host factors and treatment regimens will be critical for guiding rational microbiome-targeted therapies aimed at enhancing ICB efficacy.

Revisiting p38 Mitogen-Activated Protein Kinases (MAPK) in Inflammatory Arthritis: A Narrative of the Emergence of MAPK-Activated Protein Kinase Inhibitors (MK2i).

The p38 mitogen-activated protein kinase (p38-MAPK) is a crucial signaling pathway closely involved in several physiological and cellular functions, including cell cycle, apoptosis, gene expression, and responses to stress stimuli. It also plays a central role in inflammation and immunity. Owing to disparate p38-MAPK functions, it has thus far formed an elusive drug target with failed clinical trials in inflammatory diseases due to challenges including hepatotoxicity, cardiac toxicity, lack of efficacy, and tachyphylaxis, which is a brief initial improvement with rapid disease rebound. To overcome these limitations, downstream antagonism of the p38 pathway with a MAPK-activated protein kinase (MAPKAPK, also known as MK2) blockade has demonstrated the potential to abrogate inflammation without the prior recognized toxicities. Such MK2 inhibition (MK2i) is associated with robust suppression of key pro-inflammatory cytokines, including TNFα and IL-6 and others in experimental systems and in vitro. Considering this recent evidence regarding MK2i in inflammatory arthritis, we revisit the p38-MAPK pathway and discuss the literature encompassing the challenges of p38 inhibitors with a focus on this pathway. We then highlight how novel MK2i strategies, although encouraging in the pre-clinical arena, may either show evidence for efficacy or the lack of efficacy in emergent human trials data from different disease settings.

Polymorphisms in Cyclooxygenase, Lipoxygenase, and TP53 Genes Predict Colorectal Polyp Risk Reduction by Aspirin in the seAFOod Polyp Prevention Trial.

Aspirin and eicosapentaenoic acid (EPA) reduce colorectal adenomatous polyp risk and affect synthesis of oxylipins including prostaglandin E2. We investigated whether 35 SNPs in oxylipin metabolism genes such as cyclooxygenase (PTGS) and lipoxygenase (ALOX), as well as 7 SNPs already associated with colorectal cancer risk reduction by aspirin (e.g., TP53; rs104522), modified the effects of aspirin and EPA on colorectal polyp recurrence in the randomized 2 × 2 factorial seAFOod trial. Treatment effects were reported as the incidence rate ratio (IRR) and 95% confidence interval (CI) by stratifying negative binomial and Poisson regression analyses of colorectal polyp risk on SNP genotype. Statistical significance was reported with adjustment for the false discovery rate as the P and q value. 542 (of 707) trial participants had both genotype and colonoscopy outcome data. Reduction in colorectal polyp risk in aspirin users compared with nonaspirin users was restricted to rs4837960 (PTGS1) common homozygotes [IRR, 0.69; 95% confidence interval (CI), 0.53-0.90); q = 0.06], rs2745557 (PTGS2) compound heterozygote-rare homozygotes [IRR, 0.60 (0.41-0.88); q = 0.06], rs7090328 (ALOX5) rare homozygotes [IRR 0.27 (0.11-0.64); q = 0.05], rs2073438 (ALOX12) common homozygotes [IRR, 0.57 (0.41-0.80); q = 0.05], and rs104522 (TP53) rare homozygotes [IRR, 0.37 (0.17-0.79); q = 0.06]. No modification of colorectal polyp risk in EPA users was observed. In conclusion, genetic variants relevant to the proposed mechanism of action on oxylipins are associated with differential colorectal polyp risk reduction by aspirin in individuals who develop multiple colorectal polyps. SNP genotypes should be considered during development of personalized, predictive models of colorectal cancer chemoprevention by aspirin. PREVENTION RELEVANCE: Single-nucleotide polymorphisms in genes controlling lipid mediator signaling may modify the colorectal polyp prevention activity of aspirin. Further investigation is required to determine whether testing for genetic variants can be used to target cancer chemoprevention by aspirin to those who will benefit most.