Genetic architectures of proximal and distal colorectal cancer are partly distinct.

OBJECTIVE: An understanding of the etiologic heterogeneity of colorectal cancer (CRC) is critical for improving precision prevention, including individualized screening recommendations and the discovery of novel drug targets and repurposable drug candidates for chemoprevention. Known differences in molecular characteristics and environmental risk factors among tumors arising in different locations of the colorectum suggest partly distinct mechanisms of carcinogenesis. The extent to which the contribution of inherited genetic risk factors for CRC differs by anatomical subsite of the primary tumor has not been examined. DESIGN: To identify new anatomical subsite-specific risk loci, we performed genome-wide association study (GWAS) meta-analyses including data of 48 214 CRC cases and 64 159 controls of European ancestry. We characterised effect heterogeneity at CRC risk loci using multinomial modelling. RESULTS: We identified 13 loci that reached genome-wide significance (p<5×10-8) and that were not reported by previous GWASs for overall CRC risk. Multiple lines of evidence support candidate genes at several of these loci. We detected substantial heterogeneity between anatomical subsites. Just over half (61) of 109 known and new risk variants showed no evidence for heterogeneity. In contrast, 22 variants showed association with distal CRC (including rectal cancer), but no evidence for association or an attenuated association with proximal CRC. For two loci, there was strong evidence for effects confined to proximal colon cancer. CONCLUSION: Genetic architectures of proximal and distal CRC are partly distinct. Studies of risk factors and mechanisms of carcinogenesis, and precision prevention strategies should take into consideration the anatomical subsite of the tumour.

Lymphocytic infiltration in stage II microsatellite stable colorectal tumors: A retrospective prognosis biomarker analysis.

BACKGROUND: Identifying stage II patients with colorectal cancer (CRC) at higher risk of progression is a clinical priority in order to optimize the advantages of adjuvant chemotherapy while avoiding unnecessary toxicity. Recently, the intensity and the quality of the host immune response in the tumor microenvironment have been reported to have an important role in tumorigenesis and an inverse association with tumor progression. This association is well established in microsatellite instable CRC. In this work, we aim to assess the usefulness of measures of T-cell infiltration as prognostic biomarkers in 640 stage II, CRC tumors, 582 of them confirmed microsatellite stable. METHODS AND FINDINGS: We measured both the quantity and clonality index of T cells by means of T-cell receptor (TCR) immunosequencing in a discovery dataset (95 patients with colon cancer diagnosed at stage II and microsatellite stable, median age 67, 30% women) and replicated the results in 3 additional series of stage II patients from 2 countries. Series 1 and 2 were recruited in Barcelona, Spain and included 112 fresh frozen (FF, median age 69, 44% women) and 163 formalin-fixed paraffin-embedded (FFPE, median age 67, 39% women) samples, respectively. Series 3 included 270 FFPE samples from patients recruited in Haifa, Northern Israel, as part of a large case-control study of CRC (median age 73, 46% women). Median follow-up time was 81.1 months. Cox regression models were fitted to evaluate the prognostic value of T-cell abundance and Simpson clonality of TCR variants adjusting by sex, age, tumor location, and stage (IIA and IIB). In the discovery dataset, higher TCR abundance was associated with better prognosis (hazard ratio [HR] for ≥Q1 = 0.25, 95% CI 0.10-0.63, P = 0.003). A functional analysis of gene expression on these tumors revealed enrichment in pathways related to immune response. Higher values of clonality index (lower diversity) were not associated with worse disease-free survival, though the HR for ≥Q3 was 2.32 (95% CI 0.90-5.97, P = 0.08). These results were replicated in an independent FF dataset (TCR abundance: HR = 0.30, 95% CI 0.12-0.72, P = 0.007; clonality: HR = 3.32, 95% CI 1.38-7.94, P = 0.007). Also, the association with prognosis was tested in 2 independent FFPE datasets. The same association was observed with TCR abundance (HR = 0.41, 95% CI 0.18-0.93, P = 0.03 and HR = 0.56, 95% CI 0.31-1, P = 0.042, respectively, for each FFPE dataset). However, the clonality index was associated with prognosis only in the FFPE dataset from Israel (HR = 2.45, 95% CI 1.39-4.32, P = 0.002). Finally, a combined analysis combining all microsatellite stable (MSS) samples demonstrated a clear prognosis value both for TCR abundance (HR = 0.39, 95% CI 0.26-0.57, P = 1.3e-06) and the clonality index (HR = 2.13, 95% CI 1.44-3.15, P = 0.0002). These associations were also observed when variables were considered continuous in the models (HR per log2 of TCR abundance = 0.85, 95% CI 0.78-0.93, P = 0.0002; HR per log2 or clonality index = 1.16, 95% CI 1.03-1.31, P = 0.016). LIMITATIONS: This is a retrospective study, and samples had been preserved with different methods. Validation series lack complete information about microsatellite instability (MSI) status and pathology assessment. The Molecular Epidemiology of Colorectal Cancer (MECC) study had information about overall survival instead of progression-free survival. CONCLUSION: Results from this study demonstrate that tumor lymphocytes, assessed by TCR repertoire quantification based on a sequencing method, are an independent prognostic factor in microsatellite stable stage II CRC.

Tumour characteristics and survivorship in a cohort of breast cancer: the MCC-Spain study.

PURPOSE: The objective of this study is to analyse the relative survival with breast cancer in women diagnosed after new treatments were generalised and to ascertain the current effect that tumour characteristics such as grade, stage or subtype have on survival as well as the new AJCC-pathological prognostic score. METHODS: The breast cancer MCC-Spain follow-up study is a prospective cohort study of 1685 incident breast cancer cases. Women between 20 and 85 years old were recruited between the years 2008 and 2013 in 18 hospitals located in 10 Spanish provinces and they have been followed until 2017/2018. Relative survival was estimated after 3, 5 and 8 years of follow-up using Ederer II method. In addition, Weibull regression adjusted by age, hospital, grade and stage was used to investigate prognosis factors. RESULTS: Among components of TNM staging system, tumour size greater than 50 mm (i.e. T3 or T4) more than doubled the risk of dying, while N3 nodal involvement and presence of metastasis had a huge effect on mortality. The AJCC pathological prognostic score strongly correlated with survival; thus, hazard ratios increased as the score rose, being 2.31, 4.00, 4.94, 7.92, 2.26, 14.9 and 58.9 for scores IB, IIA, IIB, IIIA, IIIB, IIIC and IV, respectively. CONCLUSION: Both TNM staging and histological/molecular biomarkers are associated with overall survival in Spanish women with breast cancer; when both are combined in the AJCC pathological prognosis score, the prognostic value improved with risk indices that increased rapidly as the pathological prognosis score increased.

Colon-specific eQTL analysis to inform on functional SNPs.

BACKGROUND: Genome-wide association studies on colorectal cancer have identified more than 60 susceptibility loci, but for most of them there is no clear knowledge of functionality or the underlying gene responsible for the risk modification. Expression quantitative trail loci (eQTL) may provide functional information for such single nucleotide polymorphisms (SNPs). METHODS: We have performed detailed eQTL analysis specific for colon tissue on a series of 97 colon tumours, their paired adjacent normal mucosa and 47 colon mucosa samples donated by healthy individuals. R package MatrixEQTL was used to search for genome-wide cis-eQTL and trans-eQTL fitting linear models adjusted for age, gender and tissue type to rank transformed expression data. RESULTS: The cis-eQTL analyses has revealed 29,073 SNP-gene associations with permutation-adjusted P-values < 0.01. These correspond to 363 unique genes. The trans-eQTL analysis identified 10,665 significant SNP-gene associations, most of them in the same chromosome, further than 1 Mb of the gene. We provide a web tool to search for specific SNPs or genes. The tool calculates Pearson or Spearman correlation, and allows to select tissue type for analysis. Data and plots can be exported. CONCLUSIONS: This resource should be useful to prioritise SNPs for further functional studies and to identify relevant genes behind identified loci.

Telomere length alterations in microsatellite stable colorectal cancer and association with the immune response.

Telomeres are repetitive sequences (TTAGGG) located at the end of chromosomes. Telomeres progressively shorten with each cell replication cycle, ultimately leading to chromosomal instability and loss of cell viability. Telomere length anomaly appears to be one of the earliest and most prevalent genetic alterations in malignant transformation. Here we aim to estimate telomere length from whole-exome sequencing data in colon tumors and normal colonic mucosa, and to analyze the potential association of telomere length with clinical factors and gene expression in colon cancer. Reads containing at least five repetitions of the telomere sequence (TTAGGG) were extracted from the raw sequences of 42 adjacent normal-tumor paired samples. The number of reads from the tumor sample was normalized to build the Tumor Telomere Length Ratio (TTLR), considered an estimation of telomere length change in the tumor compared to the paired normal tissue. We evaluated the associations between TTLR and clinical factors, gene expression and copy number (CN) aberrations measured in the same tumor samples. Colon tumors showed significantly shorter telomeres than their paired normal samples. No significant association was observed between TTLR and gender, age, tumor location, prognosis, stromal infiltration or molecular subtypes. The functional gene set enrichment analysis showed pathways related to immune response significantly associated with TLLR. By extracting a relative measure of telomere length from whole-exome sequencing data, we have assessed that colon tumor cells predominantly shorten telomeres, and this alteration is associated with expression changes in genes related to immune response and inflammation in tumor cells.

Comprehensive analysis of copy number aberrations in microsatellite stable colon cancer in view of stromal component.

BACKGROUND: Somatic copy number aberrations (CNAs) are common acquired changes in cancer cells having an important role in the progression of colon cancer (colorectal cancer, CRC). This study aimed to perform a characterisation of CNA and their impact in gene expression. METHODS: Copy number aberrations were inferred from SNP array data in a series of 99 CRC. Copy number aberration events were calculated and used to assess the association between copy number dosage, clinical and molecular characteristics of the tumours, and gene expression changes. All analyses were adjusted for the quantity of stroma in each sample, which was inferred from gene expression data. RESULTS: High heterogeneity among samples was observed; the proportion of altered genome ranged between 0.04 and 26.6%. Recurrent CNA regions with gains were frequent in chromosomes 7p, 8q, 13q, and 20, whereas 8p, 17p, and 18 cumulated losses. A significant positive correlation was observed between the number of somatic mutations and total CNA (Spearman's r=0.42, P=0.006). Approximately 37% of genes located in CNA regions changed their level of expression and the average partial correlation (adjusted for stromal content) with copy number was 0.54 (interquartile range 0.20 to 0.81). Altered genes showed enrichment in pathways relevant for CRC. Tumours classified as CMS2 and CMS4 by the consensus molecular subtyping showed higher frequency of CNA. Losses of one small region in 1p36.33, with gene CDK11B, were associated with poor prognosis. More than 66% of the recurrent CNA were validated in the The Cancer Genome Atlas (TCGA) data when analysed with the same procedure. Furthermore, 79% of the genes with altered expression in our data were validated in the TCGA. CONCLUSIONS: Although CNA are frequent events in microsatellite stable CRC, few focal recurrent regions were found. These aberrations have strong effects on gene expression and contribute to deregulate relevant cancer pathways. Owing to the diploid nature of stromal cells, it is important to consider the purity of tumour samples to accurately calculate CNA events in CRC.

Use of non-steroidal anti-inflammatory drugs and risk of breast cancer: The Spanish Multi-Case-control (MCC) study.

BACKGROUND: The relationship between non-steroidal anti-inflammatory drug (NSAID) consumption and breast cancer has been repeatedly studied, although the results remain controversial. Most case-control studies reported that NSAID consumption protected against breast cancer, while most cohort studies did not find this effect. Most studies have dealt with NSAIDs as a whole group or with specific drugs, such aspirin, ibuprofen, or others, but not with NSAID subgroups according to the Anatomical Therapeutic Chemical Classification System; moreover, scarce attention has been paid to their effect on different tumor categories (i.e.: ductal/non-ductal, stage at diagnosis or presence of hormonal receptors). METHODS: In this case-control study, we report the NSAID - breast cancer relationship in 1736 breast cancer cases and 1895 healthy controls; results are reported stratifying by the women's characteristics (i.e.: menopausal status or body mass index category) and by tumor characteristics. RESULTS: In our study, NSAID use was associated with a 24 % reduction in breast cancer risk (Odds ratio [OR] = 0.76; 95 % Confidence Interval [CI]: 0.64-0.89), and similar results were found for acetic acid derivatives, propionic acid derivatives and COXIBs, but not for aspirin. Similar results were found in postmenopausal and premenopausal women. NSAID consumption also protected against hormone + or HER2+ cancers, but not against triple negative breast cancers. The COX-2 selectivity showed an inverse association with breast cancer (i.e. OR < 1), except in advanced clinical stage and triple negative cancers. CONCLUSION: Most NSAIDs, but not aspirin, showed an inverse association against breast cancer; this effect seems to be restricted to hormone + or HER2+ cancers.

Adoptive NK Cell Transfer as a Treatment in Colorectal Cancer Patients: Analyses of Tumour Cell Determinants Correlating With Efficacy In Vitro and In Vivo.

Background: Colorectal cancer (CRC) is a heterogeneous disease with variable mutational profile and tumour microenvironment composition that influence tumour progression and response to treatment. While chemoresistant and poorly immunogenic CRC remains a challenge, the development of new strategies guided by biomarkers could help stratify and treat patients. Allogeneic NK cell transfer emerges as an alternative against chemoresistant and poorly immunogenic CRC. Methods: NK cell-related immunological markers were analysed by transcriptomics and immunohistochemistry in human CRC samples and correlated with tumour progression and overall survival. The anti-tumour ability of expanded allogeneic NK cells using a protocol combining cytokines and feeder cells was analysed in vitro and in vivo and correlated with CRC mutational status and the expression of ligands for immune checkpoint (IC) receptors regulating NK cell activity. Results: HLA-I downmodulation and NK cell infiltration correlated with better overall survival in patients with a low-stage (II) microsatellite instability-high (MSI-H) CRC, suggesting a role of HLA-I as a prognosis biomarker and a potential benefit of NK cell immunotherapy. Activated allogeneic NK cells were able to eliminate CRC cultures without PD-1 and TIM-3 restriction but were affected by HLA-I expression. In vivo experiments confirmed the efficacy of the therapy against both HLA+ and HLA- CRC cell lines. Concomitant administration of pembrolizumab failed to improve tumour control. Conclusions: Our results reveal an immunological profile of CRC tumours in which immunogenicity (MSI-H) and immune evasion mechanisms (HLA downmodulation) favour NK cell immunosurveillance at early disease stages. Accordingly, we have shown that allogeneic NK cell therapy can target tumours expressing mutations conferring poor prognosis regardless of the expression of T cell-related inhibitory IC ligands. Overall, this study provides a rationale for a new potential basis for CRC stratification and NK cell-based therapy.

Tumor immune infiltration estimated from gene expression profiles predicts colorectal cancer relapse.

A substantial fraction of patients with stage I-III colorectal adenocarcinoma (CRC) experience disease relapse after surgery with curative intent. However, biomarkers for predicting the likelihood of CRC relapse have not been fully explored. Therefore, we assessed the association between tumor infiltration by a broad array of innate and adaptive immune cell types and CRC relapse risk. We implemented a discovery-validation design including a discovery dataset from Moffitt Cancer Center (MCC; Tampa, FL) and three independent validation datasets: (1) GSE41258 (2) the Molecular Epidemiology of Colorectal Cancer (MECC) study, and (3) GSE39582. Infiltration by 22 immune cell types was inferred from tumor gene expression data, and the association between immune infiltration by each cell type and relapse-free survival was assessed using Cox proportional hazards regression. Within each of the four independent cohorts, CD4+ memory activated T cell (HR: 0.93, 95% CI: 0.90-0.96; FDR = 0.0001) infiltration was associated with longer time to disease relapse, independent of stage, microsatellite instability, and adjuvant therapy. Based on our meta-analysis across the four datasets, 10 innate and adaptive immune cell types associated with disease relapse of which 2 were internally validated using multiplex immunofluorescence. Moreover, immune cell type infiltration was a better predictors of disease relapse than Consensus Molecular Subtype (CMS) and other expression-based biomarkers (Immune-AICMCC:238.1-238.9; CMS-AICMCC: 241.0). These data suggest that transcriptome-derived immune profiles are prognostic indicators of CRC relapse and quantification of both innate and adaptive immune cell types may serve as candidate biomarkers for predicting prognosis and guiding frequency and modality of disease surveillance.

Genetically determined telomere length and multiple myeloma risk and outcome.

Telomeres are involved in processes like cellular growth, chromosomal stability, and proper segregation to daughter cells. Telomere length measured in leukocytes (LTL) has been investigated in different cancer types, including multiple myeloma (MM). However, LTL measurement is prone to heterogeneity due to sample handling and study design (retrospective vs. prospective). LTL is genetically determined; genome-wide association studies identified 11 SNPs that, combined in a score, can be used as a genetic instrument to measure LTL and evaluate its association with MM risk. This approach has been already successfully attempted in various cancer types but never in MM. We tested the "teloscore" in 2407 MM patients and 1741 controls from the International Multiple Myeloma rESEarch (IMMeNSE) consortium. We observed an increased risk for longer genetically determined telomere length (gdTL) (OR = 1.69; 95% CI 1.36-2.11; P = 2.97 × 10-6 for highest vs. lowest quintile of the score). Furthermore, in a subset of 1376 MM patients we tested the relationship between the teloscore and MM patients survival, observing a better prognosis for longer gdTL compared with shorter gdTL (HR = 0.93; 95% CI 0.86-0.99; P = 0.049). In conclusion, we report convincing evidence that longer gdTL is a risk marker for MM risk, and that it is potentially involved in increasing MM survival.

Lung metastases share common immune features regardless of primary tumor origin.

BACKGROUND: Only certain disseminated cells are able to grow in secondary organs to create a metastatic tumor. Under the hypothesis that the immune microenvironment of the host tissue may play an important role in this process, we have categorized metastatic samples based on their immune features. METHODS: Gene expression data of metastatic samples (n=374) from four secondary sites (brain, bone, liver and lung) were used to characterize samples based on their immune and stromal infiltration using gene signatures and cell quantification tools. A clustering analysis was done that separated metastatic samples into three different immune categories: high, medium and low. RESULTS: Significant differences were found between the immune profiles of samples metastasizing in distinct organs. Metastases in lung showed a higher immunogenic score than metastases in brain, liver or bone, regardless of their primary site of origin. Also, they preferentially clustered in the high immune group. Samples in this cluster exhibited a clear inflammatory phenotype, higher levels of immune infiltrate, overexpression of programmed death-ligand 1 (PD-L1) and cytotoxic T-lymphocyte-associated protein 4 (CTLA4) pathways and upregulation of genes predicting clinical response to programmed cell death protein 1 (PD-1) blockade (T-cell inflammatory signature). A decision tree algorithm was used to select CD74 as a biomarker that identify samples belonging to this high-immune subtype of metastases, having specificity of 0.96 and sensitivity of 1. CONCLUSIONS: We have found a group of lung-enriched metastases showing an inflammatory phenotype susceptible to be treated with immunotherapy.

Changes of CD68, CD163, and PD-L1 tumor expression during high-dose-rate and pulsed-dose-rate brachytherapy for cervical cancer.

PURPOSE: We hypothesized that radiation doses delivered with high-dose-rate (HDR) and pulsed-dose-rate (PDR) brachytherapy in patients with cervical cancer could trigger immune stimulation by modulating immune cells in the tumor microenvironment. The objective was to determine CD68, CD163, and PD-L1 expression in biopsies from patients with cervical cancer and to compare the effects of HDR vs. PDR brachytherapy on the expression of these proteins. METHODS AND MATERIALS: Nineteen consecutive women (mean age, 55.9 years) with histologically proven cervical cancer scheduled for brachytherapy after standard external beam irradiation therapy combined with platinum-based chemotherapy were included in a prospective study. Core tissue biopsies were obtained before radiochemotherapy (biopsy #1), after completion of radiochemotherapy and before brachytherapy (biopsy #2), and 2 weeks after completion of brachytherapy (biopsy #3). HDR or PDR brachytherapy was delivered according to availability of equipment. CD68, CD163, and PD-L1 immunohistochemical expression was estimated using a quantitative method. Conditional logistic regression models were used to assess the relationship between gene expression and time of biopsy for each brachytherapy technique. RESULTS: In relation to CD68 and CD163, statistically significant relationships between gene expression and biopsy tissue samples were not found in any of the brachytherapy techniques, although there was trend toward downexpression of CD68 and CD163 in biopsies #2 and #3 in the HDR brachytherapy cohort only. There was a significant increase in PD-L1 expression in biopsy #3 also in the HDR cohort as compared with the PDR cohort. CONCLUSIONS: Decreased CD68 and CD163 expression did not reach statistical significance, but this trend may have clinical translational potential. Overexpression of PD-L1 in tissue biopsies taken at 14 days in the HDR brachytherapy cohort may tentatively suggest that this time period would be an adequate interval for blockade of the PD-1/PD-L1 axis.

DNA methylation events in transcription factors and gene expression changes in colon cancer.

Aim: Gain insight about the role of DNA methylation in the malignant growth of colon cancer. Patients & methods: Methylation and gene expression from 90 adjacent-tumor paired tissues and 48 healthy tissues were analyzed. Tumor genes whose change in expression was explained by changes in methylation were identified using linear models adjusted for tumor stromal content. Results: No differences in methylation were found between adjacent and healthy tissues, but clear differences were found between adjacent and tumor samples. We identified hypermethylated CpG islands located in promoter regions that drive differential gene expression of transcription factors and their target genes. Conclusion: Changes in methylation of a few genes provoke important changes in gene expression, by expanding the signal through transcription activation/repression.

NTHL1 biallelic mutations seldom cause colorectal cancer, serrated polyposis or a multi-tumor phenotype, in absence of colorectal adenomas.

The cancer-predisposing syndrome caused by biallelic mutations in NTHL1 may not be a solely colorectal cancer (CRC) and polyposis syndrome but rather a multi-tumor recessive disease. The presence of ≤10 adenomas in several mutation carriers suggests a possible causal role of NTHL1 in hereditary or early-onset nonpolyposis CRC. The involvement of NTHL1 in serrated/hyperplastic polyposis remains unexplored. The aim of our study is to elucidate the role of NTHL1 in the predisposition to personal or familial history of multiple tumor types, familial/early-onset nonpolyposis CRC, and serrated polyposis. NTHL1 mutational screening was performed in 312 cancer patients with personal or family history of multiple tumor types, 488 with hereditary nonpolyposis CRC, and 96 with serrated/hyperplastic polyposis. While no biallelic mutation carriers were identified in patients with personal and/or family history of multiple tumor types or with serrated polyposis, one was identified among the 488 nonpolyposis CRC patients. The carrier of c.268C>T (p.Q90*) and 550-1G>A was diagnosed with CRC and meningioma at ages 37 and 45 respectively, being reclassified as attenuated adenomatous polyposis after the cumulative detection of 26 adenomas. Our findings suggest that biallelic mutations in NTHL1 rarely cause CRC, a personal/familial multi-tumor history, or serrated polyposis, in absence of adenomas.

Discovery of common and rare genetic risk variants for colorectal cancer.

To further dissect the genetic architecture of colorectal cancer (CRC), we performed whole-genome sequencing of 1,439 cases and 720 controls, imputed discovered sequence variants and Haplotype Reference Consortium panel variants into genome-wide association study data, and tested for association in 34,869 cases and 29,051 controls. Findings were followed up in an additional 23,262 cases and 38,296 controls. We discovered a strongly protective 0.3% frequency variant signal at CHD1. In a combined meta-analysis of 125,478 individuals, we identified 40 new independent signals at P < 5 × 10-8, bringing the number of known independent signals for CRC to ~100. New signals implicate lower-frequency variants, Krüppel-like factors, Hedgehog signaling, Hippo-YAP signaling, long noncoding RNAs and somatic drivers, and support a role for immune function. Heritability analyses suggest that CRC risk is highly polygenic, and larger, more comprehensive studies enabling rare variant analysis will improve understanding of biology underlying this risk and influence personalized screening strategies and drug development.

Novel Common Genetic Susceptibility Loci for Colorectal Cancer.

BACKGROUND: Previous genome-wide association studies (GWAS) have identified 42 loci (P < 5 × 10-8) associated with risk of colorectal cancer (CRC). Expanded consortium efforts facilitating the discovery of additional susceptibility loci may capture unexplained familial risk. METHODS: We conducted a GWAS in European descent CRC cases and control subjects using a discovery-replication design, followed by examination of novel findings in a multiethnic sample (cumulative n = 163 315). In the discovery stage (36 948 case subjects/30 864 control subjects), we identified genetic variants with a minor allele frequency of 1% or greater associated with risk of CRC using logistic regression followed by a fixed-effects inverse variance weighted meta-analysis. All novel independent variants reaching genome-wide statistical significance (two-sided P < 5 × 10-8) were tested for replication in separate European ancestry samples (12 952 case subjects/48 383 control subjects). Next, we examined the generalizability of discovered variants in East Asians, African Americans, and Hispanics (12 085 case subjects/22 083 control subjects). Finally, we examined the contributions of novel risk variants to familial relative risk and examined the prediction capabilities of a polygenic risk score. All statistical tests were two-sided. RESULTS: The discovery GWAS identified 11 variants associated with CRC at P < 5 × 10-8, of which nine (at 4q22.2/5p15.33/5p13.1/6p21.31/6p12.1/10q11.23/12q24.21/16q24.1/20q13.13) independently replicated at a P value of less than .05. Multiethnic follow-up supported the generalizability of discovery findings. These results demonstrated a 14.7% increase in familial relative risk explained by common risk alleles from 10.3% (95% confidence interval [CI] = 7.9% to 13.7%; known variants) to 11.9% (95% CI = 9.2% to 15.5%; known and novel variants). A polygenic risk score identified 4.3% of the population at an odds ratio for developing CRC of at least 2.0. CONCLUSIONS: This study provides insight into the architecture of common genetic variation contributing to CRC etiology and improves risk prediction for individualized screening.

Colorectal cancer, sun exposure and dietary vitamin D and calcium intake in the MCC-Spain study.

OBJECTIVES: To explore the association of colorectal cancer with environmental solar radiation and sun exposure behavior, considering phenotypic variables (eye color, hair color and skin phenotype), dietary intake of vitamin D and calcium, and socio-demographic factors. STUDY DESIGN: Multicenter population-based frequency matched case-control study in Spain (MCC-Spain), with 2140 CRC cases and 3950 controls. METHODS: Data were obtained through personal interviews using a structured epidemiological questionnaire that included socio-demographic data, residential history, environmental exposures, behavior, phenotypic and dietary information. An environmental-lifetime sun exposure score was constructed combining residential history and average daily solar radiation, direct and diffuse. Logistic regression was used to explore the association between different variables. A structural equation model was used to verify the associations of the conceptual model. RESULTS: We found a lower risk of CRC in subjects frequently exposed to sunlight during the previous summer and skin burning due to sun exposure. No association was observed in relation to the residential solar radiation scores. Subjects with light eye or light hair colors had a lower risk of CRC that those with darker colors. Dietary calcium and vitamin D were also protective factors, but not in the multivariate model. The structural equation model analysis suggested that higher sun exposure was associated with a decreased risk of CRC, as well as dietary intake of calcium and vitamin D, and these factors are correlated among themselves and with environmental solar radiation and skin phenotypes. CONCLUSION: The results agree with previous observations that sun exposure, dietary vitamin D and calcium intake, and serum 25(OH)D concentration reduce the risk of CRC and indicate that these factors may be relevant for cancer prevention.

Risk Model for Colorectal Cancer in Spanish Population Using Environmental and Genetic Factors: Results from the MCC-Spain study.

Colorectal cancer (CRC) screening of the average risk population is only indicated according to age. We aim to elaborate a model to stratify the risk of CRC by incorporating environmental data and single nucleotide polymorphisms (SNP). The MCC-Spain case-control study included 1336 CRC cases and 2744 controls. Subjects were interviewed on lifestyle factors, family and medical history. Twenty-one CRC susceptibility SNPs were genotyped. The environmental risk model, which included alcohol consumption, obesity, physical activity, red meat and vegetable consumption, and nonsteroidal anti-inflammatory drug use, contributed to CRC with an average per factor OR of 1.36 (95% CI 1.27 to 1.45). Family history of CRC contributed an OR of 2.25 (95% CI 1.87 to 2.72), and each additional SNP contributed an OR of 1.07 (95% CI 1.04 to 1.10). The risk of subjects with more than 25 risk alleles (5th quintile) was 82% higher (OR 1.82, 95% CI 1.11 to 2.98) than subjects with less than 19 alleles (1st quintile). This risk model, with an AUROC curve of 0.63 (95% CI 0.60 to 0.66), could be useful to stratify individuals. Environmental factors had more weight than the genetic score, which should be considered to encourage patients to achieve a healthier lifestyle.

Mutanome and expression of immune response genes in microsatellite stable colon cancer.

The aim of this study was to analyze the impact of the mutanome in the prognosis of microsatellite stable stage II CRC tumors. The exome of 42 stage II, microsatellite stable, colon tumors (21 of them relapse) and their paired mucosa were sequenced and analyzed. Although some pathways accumulated more mutations in patients exhibiting good or poor prognosis, no single somatic mutation was associated with prognosis. Exome sequencing data is also valuable to infer tumor neoantigens able to elicit a host immune response. Hence, putative neoantigens were identified by combining information about missense mutations in each tumor and HLAs genotypes of the patients. Under the hypothesis that neoantigens should be correctly presented in order to activate the immune response, expression levels of genes involved in the antigen presentation machinery were also assessed. In addition, CD8A level (as a marker of T-cell infiltration) was measured. We found that tumors with better prognosis showed a tendency to generate a higher number of immunogenic epitopes, and up-regulated genes involved in the antigen processing machinery. Moreover, tumors with higher T-cell infiltration also showed better prognosis. Stratifying by consensus molecular subtype, CMS4 tumors showed the highest association of expression levels of genes involved in the antigen presentation machinery with prognosis. Thus, we hypothesize that a subset of stage II microsatellite stable CRC tumors are able to generate an immune response in the host via MHC class I antigen presentation, directly related with a better prognosis.

The Use of Antihypertensive Medication and the Risk of Breast Cancer in a Case-Control Study in a Spanish Population: The MCC-Spain Study.

INTRODUCTION: The evidence on the relationship between breast cancer and different types of antihypertensive drugs taken for at least 5 years is limited and inconsistent. Furthermore, the debate has recently been fueled again with new data reporting an increased risk of breast cancer among women with a long history of use of antihypertensive drugs compared with nonusers. METHODS: In this case-control study, we report the antihypertensive drugs-breast cancer relationship in 1,736 breast cancer cases and 1,895 healthy controls; results are reported stratifying by the women's characteristics (i.e., menopausal status or body mass index category) tumor characteristics and length of use of antihypertensive drugs. RESULTS: The relationship among breast cancer and use of calcium channel blockers (CCB) for 5 or more years had odds ratio (OR) = 1.77 (95% CI, 0.99 to 3.17). Stratifying by BMI, the OR increased significantly in the group with BMI ≥ 25 (OR 2.54, 95% CI, 1.24 to 5.22). CCBs were even more strongly associated with more aggressive tumors, (OR for invasive tumors = 1.96, 95% CI = 1.09 to 3.53; OR for non ductal cancers = 3.97, 95% CI = 1.73 to 9.05; OR for Erbb2+ cancer = 2.97, 95% CI: 1.20 to 7.32). On the other hand, premenopausal women were the only group in which angiotensin II receptor blockers may be associated with breast cancer (OR = 4.27, 95% CI = 1.32 to 13.84) but this could not be identified with any type or stage. Use of angiotensin-converting-enzyme inhibitors, beta blockers and diuretics were not associated with risk. CONCLUSIONS: In this large population-based study we found that long term use of calcium channel blockers is associated with some subtypes of breast cancer (and with breast cancer in overweight women).