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

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

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.

NAD pool as an antitumor target against cancer stem cells in head and neck cancer.

Head and neck squamous cell carcinoma (HNSCC) is a heterogeneous group of tumors that affect different anatomical locations. Despite this heterogeneity, HNSCC treatment depends on the anatomical location, TNM stage and resectability of the tumor. Classical chemotherapy is based on platinum-derived drugs (cisplatin, carboplatin and oxaliplatin), taxanes (docetaxel, paclitaxel) and 5-fluorouracil1. Despite advances in HNSCC treatment, the rate of tumor recurrence and patient mortality remain high. Therefore, the search for new prognostic identifiers and treatments targeting therapy-resistant tumor cells is vital. Our work demonstrates that there are different subgroups with high phenotypic plasticity within the CSC population in HNSCC. CD10, CD184, and CD166 may identify some of these CSC subpopulations with NAMPT as a common metabolic gene for the resilient cells of these subpopulations. We observed that NAMPT reduction causes a decrease in tumorigenic and stemness properties, migration capacity and CSC phenotype through NAD pool depletion. However, NAMPT-inhibited cells can acquire resistance by activating the NAPRT enzyme of the Preiss-Handler pathway. We observed that coadministration of the NAMPT inhibitor with the NAPRT inhibitor cooperated inhibiting tumor growth. The use of an NAPRT inhibitor as an adjuvant improved NAMPT inhibitor efficacy and reduced the dose and toxicity of these inhibitors. Therefore, it seems that the reduction in the NAD pool could have efficacy in tumor therapy. This was confirmed by in vitro assays supplying the cells with products of inhibited enzymes (NA, NMN or NAD) and restoring their tumorigenic and stemness properties. In conclusion, the coinhibition of NAMPT and NAPRT improved the efficacy of antitumor treatment, indicating that the reduction in the NAD pool is important to prevent tumor growth.

BMPR2 as a Novel Predisposition Gene for Hereditary Colorectal Polyposis.

BACKGROUND & AIMS: Colorectal cancer (CRC) is one of the most prevalent tumors worldwide, with incidence quickly increasing (particularly in the context of early-onset cases), despite important prevention efforts, mainly in the form of population-wide screening programs. Although many cases present a clear familial component, the current list of hereditary CRC genes leaves a considerable proportion of the cases unexplained. METHODS: In this work, we used whole-exome sequencing approaches on 19 unrelated patients with unexplained colonic polyposis to identify candidate CRC predisposition genes. The candidate genes were then validated in an additional series of 365 patients. CRISPR-Cas9 models were used to validate BMPR2 as a potential candidate for CRC risk. RESULTS: We found 8 individuals carrying 6 different variants in the BMPR2 gene (approximately 2% of our cohort of patients with unexplained colonic polyposis). CRISPR-Cas9 models of 3 of these variants showed that the p.(Asn442Thrfs∗32) truncating variant completely abrogated BMP pathway function in a similar way to the BMPR2 knockout. Missense variants p.(Asn565Ser), p.(Ser967Pro) had varying effects on cell proliferation levels, with the former impairing cell control inhibition via noncanonical pathways. CONCLUSIONS: Collectively, these results support loss-of-function BMPR2 variants as candidates to be involved in CRC germline predisposition.

Combining Asian-European Genome-Wide Association Studies of Colorectal Cancer Improves Risk Prediction Across Race and Ethnicity.

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 expanded 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 were 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 1,681-3,651 cases and 8,696-115,105 controls of Asian, Black/African American, Latinx/Hispanic, and non-Hispanic White, respectively. They were 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.

Deciphering colorectal cancer genetics through multi-omic analysis of 100,204 cases and 154,587 controls of European and east Asian ancestries.

Colorectal cancer (CRC) is a leading cause of mortality worldwide. We conducted a genome-wide association study meta-analysis of 100,204 CRC cases and 154,587 controls of European and east Asian ancestry, identifying 205 independent risk associations, of which 50 were unreported. We performed integrative genomic, transcriptomic and methylomic analyses across large bowel mucosa and other tissues. Transcriptome- and methylome-wide association studies revealed an additional 53 risk associations. We identified 155 high-confidence effector genes functionally linked to CRC risk, many of which had no previously established role in CRC. These have multiple different functions and specifically indicate that variation in normal colorectal homeostasis, proliferation, cell adhesion, migration, immunity and microbial interactions determines CRC risk. Crosstissue analyses indicated that over a third of effector genes most probably act outside the colonic mucosa. Our findings provide insights into colorectal oncogenesis and highlight potential targets across tissues for new CRC treatment and chemoprevention strategies.

Transcriptomic and Proteomic Profiles for Elucidating Cisplatin Resistance in Head-and-Neck Squamous Cell Carcinoma.

To identify the novel genes involved in chemoresistance in head and neck squamous cell carcinoma (HNSCC), we explored the expression profiles of the following cisplatin (CDDP) resistant (R) versus parental (sensitive) cell lines by RNA-sequencing (RNA-seq): JHU029, HTB-43 and CCL-138. Using the parental condition as a control, 30 upregulated and 85 downregulated genes were identified for JHU029-R cells; 263 upregulated and 392 downregulated genes for HTB-43-R cells, and 154 upregulated and 68 downregulated genes for CCL-138-R cells. Moreover, we crossed-checked the RNA-seq results with the proteomic profiles of HTB-43-R (versus HTB-43) and CCL-138-R (versus CCL-138) cell lines. For the HTB-43-R cells, 21 upregulated and 72 downregulated targets overlapped between the proteomic and transcriptomic data; whereas in CCL-138-R cells, four upregulated and three downregulated targets matched. Following an extensive literature search, six genes from the RNA-seq (CLDN1, MAGEB2, CD24, CEACAM6, IL1B and ISG15) and six genes from the RNA-seq and proteomics crossover (AKR1C3, TNFAIP2, RAB7A, LGALS3BP, PSCA and SSRP1) were selected to be studied by qRT-PCR in 11 HNSCC patients: six resistant and five sensitive to conventional therapy. Interestingly, the high MAGEB2 expression was associated with resistant tumours and is revealed as a novel target to sensitise resistant cells to therapy in HNSCC patients.

Liquid Biopsy Biomarkers for Immunotherapy in Non-Small Cell Lung Carcinoma: Lessons Learned and the Road Ahead.

Over the recent years, advances in the development of anti-cancer treatments, particularly the implementation of ICIs (immune checkpoint inhibitors), have resulted in increased survival rates in NSCLC (non-small cell lung cancer) patients. However, a significant proportion of patients does not seem respond to immunotherapy, and some individuals even develop secondary resistance to treatment. Therefore, it is imperative to correctly identify the patients that will benefit from ICI therapy in order to tailor therapeutic options in an individualised setting, ultimately benefitting both the patient and the health system. Many different biomarkers have been explored to correctly stratify patients and predict response to immunotherapy, but liquid biopsy approaches have recently arisen as an interesting opportunity to predict and monitor treatment response due to their logistic accessibility. This review summarises the current data and efforts in the field of ICI response biomarkers in NSCLC patients and highlights advantages and limitations as we discuss the road to clinical implementation.

The Road so Far in Colorectal Cancer Pharmacogenomics: Are We Closer to Individualised Treatment?

In recent decades, survival rates in colorectal cancer have improved greatly due to pharmacological treatment. However, many patients end up developing adverse drug reactions that can be severe or even life threatening, and that affect their quality of life. These remain a limitation, as they may force dose reduction or treatment discontinuation, diminishing treatment efficacy. From candidate gene approaches to genome-wide analysis, pharmacogenomic knowledge has advanced greatly, yet there is still huge and unexploited potential in the use of novel technologies such as next-generation sequencing strategies. This review summarises the road of colorectal cancer pharmacogenomics so far, presents considerations and directions to be taken for further works and discusses the path towards implementation into clinical practice.

Tumor Profiling at the Service of Cancer Therapy.

Cancer treatment options have evolved significantly in the past few years. From the initial surgical procedures, to the latest next-generation technologies, we are now in the position to analyze and understand tumors in a one-by-one basis and use that to our advantage to provide with individualized treatment options that may increase patient survival. In this review, we will focus on how tumor profiling has evolved over the past decades to deliver more efficient and personalized treatment options, and how novel technologies can help us envisage the future of precision oncology toward a better management and, ultimately, increased survival.

Systematic meta-analyses, field synopsis and global assessment of the evidence of genetic association studies in colorectal cancer.

OBJECTIVE: To provide an understanding of the role of common genetic variations in colorectal cancer (CRC) risk, we report an updated field synopsis and comprehensive assessment of evidence to catalogue all genetic markers for CRC (CRCgene2). DESIGN: We included 869 publications after parallel literature review and extracted data for 1063 polymorphisms in 303 different genes. Meta-analyses were performed for 308 single nucleotide polymorphisms (SNPs) in 158 different genes with at least three independent studies available for analysis. Scottish, Canadian and Spanish data from genome-wide association studies (GWASs) were incorporated for the meta-analyses of 132 SNPs. To assess and classify the credibility of the associations, we applied the Venice criteria and Bayesian False-Discovery Probability (BFDP). Genetic associations classified as 'positive' and 'less-credible positive' were further validated in three large GWAS consortia conducted in populations of European origin. RESULTS: We initially identified 18 independent variants at 16 loci that were classified as 'positive' polymorphisms for their highly credible associations with CRC risk and 59 variants at 49 loci that were classified as 'less-credible positive' SNPs; 72.2% of the 'positive' SNPs were successfully replicated in three large GWASs and the ones that were not replicated were downgraded to 'less-credible' positive (reducing the 'positive' variants to 14 at 11 loci). For the remaining 231 variants, which were previously reported, our meta-analyses found no evidence to support their associations with CRC risk. CONCLUSION: The CRCgene2 database provides an updated list of genetic variants related to CRC risk by using harmonised methods to assess their credibility.

The polymorphic variant rs1800734 influences methylation acquisition and allele-specific TFAP4 binding in the MLH1 promoter leading to differential mRNA expression.

Expression of the mismatch repair gene MutL homolog 1 (MLH1) is silenced in a clinically important subgroup of sporadic colorectal cancers. These cancers exhibit hypermutability with microsatellite instability (MSI) and differ from microsatellite-stable (MSS) colorectal cancers in both prognosis and response to therapies. Loss of MLH1 is usually due to epigenetic silencing with associated promoter methylation; coding somatic mutations rarely occur. Here we use the presence of a colorectal cancer (CRC) risk variant (rs1800734) within the MLH1 promoter to investigate the poorly understood mechanisms of MLH1 promoter methylation and loss of expression. We confirm the association of rs1800734 with MSI+ but not MSS cancer risk in our own data and by meta-analysis. Using sensitive allele-specific detection methods, we demonstrate that MLH1 is the target gene for rs1800734 mediated cancer risk. In normal colon tissue, small allele-specific differences exist only in MLH1 promoter methylation, but not gene expression. In contrast, allele-specific differences in both MLH1 methylation and expression are present in MSI+ cancers. We show that MLH1 transcriptional repression is dependent on DNA methylation and can be reversed by a methylation inhibitor. The rs1800734 allele influences the rate of methylation loss and amount of re-expression. The transcription factor TFAP4 binds to the rs1800734 region but with much weaker binding to the risk than the protective allele. TFAP4 binding is absent on both alleles when promoter methylation is present. Thus we propose that TFAP4 binding shields the protective rs1800734 allele of the MLH1 promoter from BRAF induced DNA methylation more effectively than the risk allele.

Early Colorectal Cancers Provide New Evidence for a Lynch Syndrome-to-CMMRD Phenotypic Continuum.

Lynch syndrome (LS) is the most common hereditary colorectal cancer (CRC) syndrome, caused by heterozygous mutations in the mismatch repair (MMR) genes. Biallelic mutations in these genes lead however, to constitutive mismatch repair deficiency (CMMRD). In this study, we follow the diagnostic journey of a 12-year old patient with CRC, with a clinical phenotype overlapping CMMRD. We perform molecular and functional assays to discard a CMMRD diagnosis then identify by exome sequencing and validation in a cohort of 134 LS patients, a candidate variant in the MLH1 UTR region in homozygosis. We propose that this variant, together with other candidates, could be responsible for age-of-onset modulation. Our data support the idea that low-risk modifier alleles may influence early development of cancer in LS leading to a LS-to-CMMRD phenotypic continuum. Therefore, it is essential that larger efforts are directed to the identification and study of these genetic modifiers, in order to provide optimal cancer prevention strategies to these patients.

Association analyses identify 31 new risk loci for colorectal cancer susceptibility.

Colorectal cancer (CRC) is a leading cause of cancer-related death worldwide, and has a strong heritable basis. We report a genome-wide association analysis of 34,627 CRC cases and 71,379 controls of European ancestry that identifies SNPs at 31 new CRC risk loci. We also identify eight independent risk SNPs at the new and previously reported European CRC loci, and a further nine CRC SNPs at loci previously only identified in Asian populations. We use in situ promoter capture Hi-C (CHi-C), gene expression, and in silico annotation methods to identify likely target genes of CRC SNPs. Whilst these new SNP associations implicate target genes that are enriched for known CRC pathways such as Wnt and BMP, they also highlight novel pathways with no prior links to colorectal tumourigenesis. These findings provide further insight into CRC susceptibility and enhance the prospects of applying genetic risk scores to personalised screening and prevention.

Patterns of genetic differentiation and the footprints of historical migrations in the Iberian Peninsula.

The Iberian Peninsula is linguistically diverse and has a complex demographic history, including a centuries-long period of Muslim rule. Here, we study the fine-scale genetic structure of its population, and the genetic impacts of historical events, leveraging powerful, haplotype-based statistical methods to analyse 1413 individuals from across Spain. We detect extensive fine-scale population structure at extremely fine scales (below 10 Km) in some regions, including Galicia. We identify a major east-west axis of genetic differentiation, and evidence of historical north to south population movement. We find regionally varying fractions of north-west African ancestry (0-11%) in modern-day Iberians, related to an admixture event involving European-like and north-west African-like source populations. We date this event to 860-1120 CE, implying greater genetic impacts in the early half of Muslim rule in Iberia. Together, our results indicate clear genetic impacts of population movements associated with both the Muslim conquest and the subsequent Reconquista.

Author Correction: Telomere length and genetics are independent colorectal tumour risk factors in an evaluation of biomarkers in normal bowel.

Since the publication of this paper, the authors noticed that James E. East was assigned to the incorrect affiliation. The affiliation information is provided correctly, above.

Telomere length and genetics are independent colorectal tumour risk factors in an evaluation of biomarkers in normal bowel.

BACKGROUND: Colorectal cancer (CRC) screening might be improved by using a measure of prior risk to modulate screening intensity or the faecal immunochemical test threshold. Intermediate molecular biomarkers could aid risk prediction by capturing both known and unknown risk factors. METHODS: We sampled normal bowel mucosa from the proximal colon, distal colon and rectum of 317 individuals undergoing colonoscopy. We defined cases as having a personal history of colorectal polyp(s)/cancer, and controls as having no history of colorectal neoplasia. Molecular analyses were performed for: telomere length (TL); global methylation; and the expression of genes in molecular pathways associated with colorectal tumourigenesis. We also calculated a polygenic risk score (PRS) based on CRC susceptibility polymorphisms. RESULTS: Bowel TL was significantly longer in cases than controls, but was not associated with blood TL. PRS was significantly and independently higher in cases. Hypermethylation showed a suggestive association with case:control status. No gene or pathway was differentially expressed between cases and controls. Gene expression often varied considerably between bowel locations. CONCLUSIONS: PRS and bowel TL (but not blood TL) may be clinically-useful predictors of CRC risk. Sample collection to assess these biomarkers is feasible in clinical practice, especially where population screening uses flexible sigmoidoscopy or colonoscopy.

Thyroid cancer GWAS identifies 10q26.12 and 6q14.1 as novel susceptibility loci and reveals genetic heterogeneity among populations.

Thyroid cancer is the most heritable cancer of all those not displaying typical Mendelian inheritance. However, most of the genetic factors that would explain the high heritability remain unknown. Our aim was to identify additional common genetic variants associated with susceptibility to this disease. In order to do so, we performed a genome-wide association study in a series of 398 cases and 502 controls from Spain, followed by a replication in four well-defined Southern European case-control collections contributing a total of 1,422 cases and 1,908 controls. The association between the variation at the 9q22 locus near FOXE1 and thyroid cancer risk was consistent across all series, with several SNPs identified (rs7028661: OR = 1.64, p = 1.0 × 10(-22) , rs7037324: OR = 1.54, p = 1.2 × 10(-17) ). Moreover, the rare alleles of three SNPs (rs2997312, rs10788123 and rs1254167) at 10q26.12 showed suggestive evidence of association with higher risk of the disease (OR = 1.35, p = 1.2 × 10(-04) , OR = 1.26, p = 5.2 × 10(-04) and OR = 1.38, p = 5.9 × 10(-05) , respectively). Finally, the rare allele of rs4075570 at 6q14.1 conferred protection in the series studied (OR = 0.82, p = 2.0 × 10(-04) ). This study suggests that heterogeneity in genetic susceptibility between populations is a key feature to take into account when exploring genetic risk factors related to this disease.

A candidate gene study of capecitabine-related toxicity in colorectal cancer identifies new toxicity variants at DPYD and a putative role for ENOSF1 rather than TYMS.

OBJECTIVE: Capecitabine is an oral 5-fluorouracil (5-FU) pro-drug commonly used to treat colorectal carcinoma and other tumours. About 35% of patients experience dose-limiting toxicity. The few proven genetic biomarkers of 5-FU toxicity are rare variants and polymorphisms, respectively, at candidate loci dihydropyrimidine dehydrogenase (DPYD) and thymidylate synthase (TYMS). DESIGN: We investigated 1456 polymorphisms and rare coding variants near 25 candidate 5-FU pathway genes in 968 UK patients from the QUASAR2 clinical trial. RESULTS: We identified the first common DPYD polymorphisms to be consistently associated with capecitabine toxicity, rs12132152 (toxicity allele frequency (TAF)=0.031, OR=3.83, p=4.31×10(-6)) and rs12022243 (TAF=0.196, OR=1.69, p=2.55×10(-5)). rs12132152 was particularly strongly associated with hand-foot syndrome (OR=6.1, p=3.6×10(-8)). The rs12132152 and rs12022243 associations were independent of each other and of previously reported DPYD toxicity variants. Next-generation sequencing additionally identified rare DPYD variant p.Ala551Thr in one patient with severe toxicity. Using functional predictions and published data, we assigned p.Ala551Thr as causal for toxicity. We found that polymorphism rs2612091, which lies within an intron of ENOSF1, was also associated with capecitabine toxicity (TAF=0.532, OR=1.59, p=5.28×10(-6)). ENSOF1 is adjacent to TYMS and there is a poorly characterised regulatory interaction between the two genes/proteins. Unexpectedly, rs2612091 fully explained the previously reported associations between capecitabine toxicity and the supposedly functional TYMS variants, 5'VNTR 2R/3R and 3'UTR 6 bp ins-del. rs2612091 genotypes were, moreover, consistently associated with ENOSF1 mRNA levels, but not with TYMS expression. CONCLUSIONS: DPYD harbours rare and common capecitabine toxicity variants. The toxicity polymorphism in the TYMS region may actually act through ENOSF1.