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.

Genetic Variant Associated With Survival of Patients With Stage II-III Colon Cancer.

BACKGROUND & AIMS: Many genetic variants have been associated with colorectal cancer risk, although few have been associated with survival times of patients. Identification of genetic variants associated with survival times might improve our understanding of disease progression and aid in outcome prediction. We performed a genome-wide association study to identify variants associated with colon cancer survival time. METHODS: We performed a post hoc analysis of data from NCCTG N0147 (Alliance), a randomized phase 3 trial of patients with resected stage III colon cancer, and from NSABP C-08 (NRG), a phase 3 trial that compared therapy regimens for patients with resected stage II or III colon cancer. Genotype analyses were performed on DNA from blood samples from 4974 patients. We used Cox proportional hazards regression to evaluate the association of each single nucleotide polymorphism with times of overall survival and disease-free survival, adjusting for age at diagnosis, sex, treatment group, and principal components of genetic ancestry. We performed the analysis for studies N0147 and C-08 separately, and results were combined in a fixed-effects meta-analysis. RESULTS: A locus on chromosome 7p15.2 was significantly associated with overall survival time (P ≤ 5x10-08). The most significant variant at this locus, rs76766811 (P = 1.6x10-08), is common among African Americans (minor allele frequency, approximately 18%) but rare in European Americans (minor allele frequency <0.1%). Within strata of self-reported ancestry, this variant was associated with times of overall survival and disease-free survival in only African Americans (hazard ratio for overall survival, 2.82; 95% CI, 1.88-4.23; P = 5.0x10-07 and hazard ratio for disease-free survival, 2.27; 95% CI, 1.62-3.18; P = 1.8x10-06). CONCLUSIONS: In an analysis of data from 2 trials of patients with stage II or III colon cancer, we identified rs76766811 as a potential prognostic variant in African American patients. This finding should be confirmed in additional study populations. ClinicalTrials.gov Identifiers: NCT00096278 (NSABP C-08) and NCT00079274 (NCCTG N0147).

Circulating bilirubin levels and risk of colorectal cancer: serological and Mendelian randomization analyses.

BACKGROUND: Bilirubin, a byproduct of hemoglobin breakdown and purported anti-oxidant, is thought to be cancer preventive. We conducted complementary serological and Mendelian randomization (MR) analyses to investigate whether alterations in circulating levels of bilirubin are associated with risk of colorectal cancer (CRC). We decided a priori to perform analyses separately in men and women based on suggestive evidence that associations may differ by sex. METHODS: In a case-control study nested in the European Prospective Investigation into Cancer and Nutrition (EPIC), pre-diagnostic unconjugated bilirubin (UCB, the main component of total bilirubin) concentrations were measured by high-performance liquid chromatography in plasma samples of 1386 CRC cases and their individually matched controls. Additionally, 115 single-nucleotide polymorphisms (SNPs) robustly associated (P < 5 × 10-8) with circulating total bilirubin were instrumented in a 2-sample MR to test for a potential causal effect of bilirubin on CRC risk in 52,775 CRC cases and 45,940 matched controls in the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO), the Colon Cancer Family Registry (CCFR), and the Colorectal Transdisciplinary (CORECT) study. RESULTS: The associations between circulating UCB levels and CRC risk differed by sex (Pheterogeneity = 0.008). Among men, higher levels of UCB were positively associated with CRC risk (odds ratio [OR] = 1.19, 95% confidence interval [CI] = 1.04-1.36; per 1-SD increment of log-UCB). In women, an inverse association was observed (OR = 0.86 (0.76-0.97)). In the MR analysis of the main UGT1A1 SNP (rs6431625), genetically predicted higher levels of total bilirubin were associated with a 7% increase in CRC risk in men (OR = 1.07 (1.02-1.12); P = 0.006; per 1-SD increment of total bilirubin), while there was no association in women (OR = 1.01 (0.96-1.06); P = 0.73). Raised bilirubin levels, predicted by instrumental variables excluding rs6431625, were suggestive of an inverse association with CRC in men, but not in women. These differences by sex did not reach formal statistical significance (Pheterogeneity ≥ 0.2). CONCLUSIONS: Additional insight into the relationship between circulating bilirubin and CRC is needed in order to conclude on a potential causal role of bilirubin in CRC development.

Adiposity, metabolites, and colorectal cancer risk: Mendelian randomization study.

BACKGROUND: Higher adiposity increases the risk of colorectal cancer (CRC), but whether this relationship varies by anatomical sub-site or by sex is unclear. Further, the metabolic alterations mediating the effects of adiposity on CRC are not fully understood. METHODS: We examined sex- and site-specific associations of adiposity with CRC risk and whether adiposity-associated metabolites explain the associations of adiposity with CRC. Genetic variants from genome-wide association studies of body mass index (BMI) and waist-to-hip ratio (WHR, unadjusted for BMI; N = 806,810), and 123 metabolites from targeted nuclear magnetic resonance metabolomics (N = 24,925), were used as instruments. Sex-combined and sex-specific Mendelian randomization (MR) was conducted for BMI and WHR with CRC risk (58,221 cases and 67,694 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium, Colorectal Cancer Transdisciplinary Study, and Colon Cancer Family Registry). Sex-combined MR was conducted for BMI and WHR with metabolites, for metabolites with CRC, and for BMI and WHR with CRC adjusted for metabolite classes in multivariable models. RESULTS: In sex-specific MR analyses, higher BMI (per 4.2 kg/m2) was associated with 1.23 (95% confidence interval (CI) = 1.08, 1.38) times higher CRC odds among men (inverse-variance-weighted (IVW) model); among women, higher BMI (per 5.2 kg/m2) was associated with 1.09 (95% CI = 0.97, 1.22) times higher CRC odds. WHR (per 0.07 higher) was more strongly associated with CRC risk among women (IVW OR = 1.25, 95% CI = 1.08, 1.43) than men (IVW OR = 1.05, 95% CI = 0.81, 1.36). BMI or WHR was associated with 104/123 metabolites at false discovery rate-corrected P ≤ 0.05; several metabolites were associated with CRC, but not in directions that were consistent with the mediation of positive adiposity-CRC relations. In multivariable MR analyses, associations of BMI and WHR with CRC were not attenuated following adjustment for representative metabolite classes, e.g., the univariable IVW OR for BMI with CRC was 1.12 (95% CI = 1.00, 1.26), and this became 1.11 (95% CI = 0.99, 1.26) when adjusting for cholesterol in low-density lipoprotein particles. CONCLUSIONS: Our results suggest that higher BMI more greatly raises CRC risk among men, whereas higher WHR more greatly raises CRC risk among women. Adiposity was associated with numerous metabolic alterations, but none of these explained associations between adiposity and CRC. More detailed metabolomic measures are likely needed to clarify the mechanistic pathways.

Genetic variant predictors of gene expression provide new insight into risk of colorectal cancer.

Genome-wide association studies have reported 56 independently associated colorectal cancer (CRC) risk variants, most of which are non-coding and believed to exert their effects by modulating gene expression. The computational method PrediXcan uses cis-regulatory variant predictors to impute expression and perform gene-level association tests in GWAS without directly measured transcriptomes. In this study, we used reference datasets from colon (n = 169) and whole blood (n = 922) transcriptomes to test CRC association with genetically determined expression levels in a genome-wide analysis of 12,186 cases and 14,718 controls. Three novel associations were discovered from colon transverse models at FDR ≤ 0.2 and further evaluated in an independent replication including 32,825 cases and 39,933 controls. After adjusting for multiple comparisons, we found statistically significant associations using colon transcriptome models with TRIM4 (discovery P = 2.2 × 10- 4, replication P = 0.01), and PYGL (discovery P = 2.3 × 10- 4, replication P = 6.7 × 10- 4). Interestingly, both genes encode proteins that influence redox homeostasis and are related to cellular metabolic reprogramming in tumors, implicating a novel CRC pathway linked to cell growth and proliferation. Defining CRC risk regions as one megabase up- and downstream of one of the 56 independent risk variants, we defined 44 non-overlapping CRC-risk regions. Among these risk regions, we identified genes associated with CRC (P < 0.05) in 34/44 CRC-risk regions. Importantly, CRC association was found for two genes in the previously reported 2q25 locus, CXCR1 and CXCR2, which are potential cancer therapeutic targets. These findings provide strong candidate genes to prioritize for subsequent laboratory follow-up of GWAS loci. This study is the first to implement PrediXcan in a large colorectal cancer study and findings highlight the utility of integrating transcriptome data in GWAS for discovery of, and biological insight into, risk loci.

Correction to: Genetic variant predictors of gene expression provide new insight into risk of colorectal cancer.

Every author has erroneously been assigned to the affiliation "62". The affiliation 62 belongs to the author Graham Casey.

Association Between Intake of Red and Processed Meat and Survival in Patients With Colorectal Cancer in a Pooled Analysis.

BACKGROUND & AIMS: Red and processed meat intake is associated with colorectal cancer (CRC) incidence, but it is not clear if intake is associated with patient survival after diagnosis. METHODS: We pooled data from 7627 patients with stage I-IV CRC from 10 studies in the International Survival Analysis in Colorectal Cancer Consortium. Cox proportional hazards regression models were used to evaluate the associations of intake of red and processed meat before diagnosis with overall and CRC-specific survival. RESULTS: Among 7627 patients with CRC, 2338 died, including 1576 from CRC, over a median follow-up time of 5.1 years. In multivariable-adjusted analyses, higher intake of red or processed meat was not associated with overall survival of patients with stage I-III CRC: Q4 vs Q1 red meat hazard ratio [HR], 1.08 (95% CI, 0.93-1.26) and Q4 vs Q1 processed meat HR, 1.10 (95% CI, 0.93-1.32) or with CRC-specific survival: Q4 vs Q1 red meat HR, 1.09 (95% CI, 0.89-1.33) and Q4 vs Q1 processed meat HR, 1.11 (95% CI, 0.87-1.42). Results were similar for patients with stage IV CRC. However, patients with stage I-III CRC who reported an intake of processed meat above the study-specific medians had a higher risk of death from any cause (HR, 1.12; 95% CI, 1.01-1.25) than patients who reported eating at or less than the median. CONCLUSION: In this large consortium of CRC patient cohorts, intake of red and processed meat before a diagnosis of CRC was not associated with shorter survival time after diagnosis, although a possible weak adverse association cannot be excluded. Studies that evaluate dietary data from several time points before and after cancer diagnosis are required to confirm these findings.

Mendelian randomisation study of age at menarche and age at menopause and the risk of colorectal cancer.

BACKGROUND: Substantial evidence supports an association between use of menopausal hormone therapy and decreased colorectal cancer (CRC) risk, indicating a role of exogenous sex hormones in CRC development. However, findings on endogenous oestrogen exposure and CRC are inconsistent. METHODS: We used a Mendelian randomisation approach to test for a causal effect of age at menarche and age at menopause as surrogates for endogenous oestrogen exposure on CRC risk. Weighted genetic risk scores based on 358 single-nucleotide polymorphisms associated with age at menarche and 51 single-nucleotide polymorphisms associated with age at menopause were used to estimate the association with CRC risk using logistic regression in 12,944 women diagnosed with CRC and 10,741 women without CRC from three consortia. Sensitivity analyses were conducted to address pleiotropy and possible confounding by body mass index. RESULTS: Genetic risk scores for age at menarche (odds ratio per year 0.98, 95% confidence interval: 0.95-1.02) and age at menopause (odds ratio 0.98, 95% confidence interval: 0.94-1.01) were not significantly associated with CRC risk. The sensitivity analyses yielded similar results. CONCLUSIONS: Our study does not support a causal relationship between genetic risk scores for age at menarche and age at menopause and CRC risk.

Phylogenomics of a rapid radiation: is chromosomal evolution linked to increased diversification in north american spiny lizards (Genus Sceloporus)?

BACKGROUND: Resolving the short phylogenetic branches that result from rapid evolutionary diversification often requires large numbers of loci. We collected targeted sequence capture data from 585 nuclear loci (541 ultraconserved elements and 44 protein-coding genes) to estimate the phylogenetic relationships among iguanian lizards in the North American genus Sceloporus. We tested for diversification rate shifts to determine if rapid radiation in the genus is correlated with chromosomal evolution. RESULTS: The phylogenomic trees that we obtained for Sceloporus using concatenation and coalescent-based species tree inference provide strong support for the monophyly and interrelationships among nearly all major groups. The diversification analysis supported one rate shift on the Sceloporus phylogeny approximately 20-25 million years ago that is associated with the doubling of the speciation rate from 0.06 species/million years (Ma) to 0.15 species/Ma. The posterior probability for this rate shift occurring on the branch leading to the Sceloporus species groups exhibiting increased chromosomal diversity is high (posterior probability = 0.997). CONCLUSIONS: Despite high levels of gene tree discordance, we were able to estimate a phylogenomic tree for Sceloporus that solves some of the taxonomic problems caused by previous analyses of fewer loci. The taxonomic changes that we propose using this new phylogenomic tree help clarify the number and composition of the major species groups in the genus. Our study provides new evidence for a putative link between chromosomal evolution and the rapid divergence and radiation of Sceloporus across North America.

Short Tree, Long Tree, Right Tree, Wrong Tree: New Acquisition Bias Corrections for Inferring SNP Phylogenies.

Single nucleotide polymorphisms (SNPs) are useful markers for phylogenetic studies owing in part to their ubiquity throughout the genome and ease of collection. Restriction site associated DNA sequencing (RADseq) methods are becoming increasingly popular for SNP data collection, but an assessment of the best practises for using these data in phylogenetics is lacking. We use computer simulations, and new double digest RADseq (ddRADseq) data for the lizard family Phrynosomatidae, to investigate the accuracy of RAD loci for phylogenetic inference. We compare the two primary ways RAD loci are used during phylogenetic analysis, including the analysis of full sequences (i.e., SNPs together with invariant sites), or the analysis of SNPs on their own after excluding invariant sites. We find that using full sequences rather than just SNPs is preferable from the perspectives of branch length and topological accuracy, but not of computational time. We introduce two new acquisition bias corrections for dealing with alignments composed exclusively of SNPs, a conditional likelihood method and a reconstituted DNA approach. The conditional likelihood method conditions on the presence of variable characters only (the number of invariant sites that are unsampled but known to exist is not considered), while the reconstituted DNA approach requires the user to specify the exact number of unsampled invariant sites prior to the analysis. Under simulation, branch length biases increase with the amount of missing data for both acquisition bias correction methods, but branch length accuracy is much improved in the reconstituted DNA approach compared to the conditional likelihood approach. Phylogenetic analyses of the empirical data using concatenation or a coalescent-based species tree approach provide strong support for many of the accepted relationships among phrynosomatid lizards, suggesting that RAD loci contain useful phylogenetic signal across a range of divergence times despite the presence of missing data. Phylogenetic analysis of RAD loci requires careful attention to model assumptions, especially if downstream analyses depend on branch lengths.

Longitudinal study of immunity to SARS-CoV2 in ocrelizumab-treated MS patients up to 2 years after COVID-19 vaccination.

OBJECTIVES: (1) To plot the trajectory of humoral and cellular immune responses to the primary (two-dose) COVID-19 mRNA series and the third/booster dose in B-cell-depleted multiple sclerosis (MS) patients up to 2 years post-vaccination; (2) to identify predictors of immune responses to vaccination; and (3) to assess the impact of intercurrent COVID-19 infections on SARS CoV-2-specific immunity. METHODS: Sixty ocrelizumab-treated MS patients were enrolled from NYU (New York) and University of Colorado (Anschutz) MS Centers. Samples were collected pre-vaccination, and then 4, 12, 24, and 48 weeks post-primary series, and 4, 12, 24, and 48 weeks post-booster. Binding anti-Spike antibody responses were assessed with multiplex bead-based immunoassay (MBI) and electrochemiluminescence (Elecsys®, Roche Diagnostics), and neutralizing antibody responses with live-virus immunofluorescence-based microneutralization assay. Spike-specific cellular responses were assessed with IFNγ/IL-2 ELISpot (Invitrogen) and, in a subset, by sequencing complementarity determining regions (CDR)-3 within T-cell receptors (Adaptive Biotechnologies). A linear mixed-effect model was used to compare antibody and cytokine levels across time points. Multivariate analyses identified predictors of immune responses. RESULTS: The primary vaccination induced an 11- to 208-fold increase in binding and neutralizing antibody levels and a 3- to 4-fold increase in IFNγ/IL-2 responses, followed by a modest decline in antibody but not cytokine responses. Booster dose induced a further 3- to 5-fold increase in binding antibodies and 4- to 5-fold increase in IFNγ/IL-2, which were maintained for up to 1 year. Infections had a variable impact on immunity. INTERPRETATION: Humoral and cellular benefits of COVID-19 vaccination in B-cell-depleted MS patients were sustained for up to 2 years when booster doses were administered.

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

BACKGROUND: Fusobacterium nucleatum (F. nucleatum) activates oncogenic signaling pathways and induces inflammation to promote colorectal carcinogenesis. METHODS: We characterized F. nucleatum and its subspecies in colorectal tumors and examined associations with tumor characteristics and colorectal cancer-specific survival. We conducted deep sequencing of nusA, nusG, and bacterial 16s rRNA genes in tumors from 1,994 patients with colorectal cancer and assessed associations between F. nucleatum presence and clinical characteristics, colorectal cancer-specific mortality, and somatic mutations. RESULTS: F. nucleatum, which was present in 10.3% of tumors, was detected in a higher proportion of right-sided and advanced-stage tumors, particularly subspecies animalis. Presence of F. nucleatum was associated with higher colorectal cancer-specific mortality (HR, 1.97; P = 0.0004). This association was restricted to nonhypermutated, microsatellite-stable tumors (HR, 2.13; P = 0.0002) and those who received chemotherapy [HR, 1.92; confidence interval (CI), 1.07-3.45; P = 0.029). Only F. nucleatum subspecies animalis, the main subspecies detected (65.8%), was associated with colorectal cancer-specific mortality (HR, 2.16; P = 0.0016), subspecies vincentii and nucleatum were not (HR, 1.07; P = 0.86). Additional adjustment for tumor stage suggests that the effect of F. nucleatum on mortality is partly driven by a stage shift. Presence of F. nucleatum was associated with microsatellite instable tumors, tumors with POLE exonuclease domain mutations, and ERBB3 mutations, and suggestively associated with TP53 mutations. CONCLUSIONS: F. nucleatum, and particularly subspecies animalis, was associated with a higher colorectal cancer-specific mortality and specific somatic mutated genes. IMPACT: Our findings identify the F. nucleatum subspecies animalis as negatively impacting colorectal cancer mortality, which may occur through a stage shift and its effect on chemoresistance.

Genetic Predictors of Severe Skin Toxicity in Patients with Stage III Colon Cancer Treated with Cetuximab: NCCTG N0147 (Alliance).

BACKGROUND: Cetuximab, an EGFR inhibitor used to treat multiple cancer types, including colon cancer, causes severe skin toxicity in 5%-20% of patients, leading to decreased quality of life and treatment delays. Our understanding of which patients have an increased risk of severe toxicities is limited. We conducted a genome-wide association study to identify germline variants predictive of cetuximab-induced severe skin toxicity. METHODS: Our study included 1,209 patients with stage III colon cancer randomized to receive cetuximab plus 5-fluorouracil and oxaliplatin as part of the NCCTG N0147 (Alliance) clinical trial. Skin toxicity outcomes were collected using the Common Toxicity Criteria for Adverse Events version 3.0. We performed genotyping, evaluating approximately 10 million genetic variants. We used logistic regression to evaluate the association of each genetic variant and severe (grade ≥ 3) skin toxicity, adjusting for age, sex, and genetic ancestry. Genome-wide significance was defined as P < 5 × 10-8. RESULTS: Participants were predominantly middle-aged white men; 20% (n = 243) experienced severe skin toxicity. Two genetic variants in the retinoic acid receptor alpha (RARA) gene were significantly associated with severe skin toxicity [OR, 3.93; 95% confidence interval (CI), 2.47-6.25; P < 7.8 × 10-9]. Functional annotations indicate these variants are in the RARA promoter. Additional significantly associated variants were identified in chromosome 2 intergenic regions. CONCLUSIONS: Identified variants could represent a potential target for risk stratification of patients with colon cancer receiving cetuximab. IMPACT: Retinoids have shown promise in the treatment of cetuximab-induced skin toxicity, so follow-up work could evaluate whether individuals with the RARA variant would benefit from retinoid therapy.

Association Between Smoking and Molecular Subtypes of Colorectal Cancer.

Background: Smoking is associated with colorectal cancer (CRC) risk. Previous studies suggested this association may be restricted to certain molecular subtypes of CRC, but large-scale comprehensive analysis is lacking. Methods: A total of 9789 CRC cases and 11 231 controls of European ancestry from 11 observational studies were included. We harmonized smoking variables across studies and derived sex study-specific quartiles of pack-years of smoking for analysis. Four somatic colorectal tumor markers were assessed individually and in combination, including BRAF mutation, KRAS mutation, CpG island methylator phenotype (CIMP), and microsatellite instability (MSI) status. A multinomial logistic regression analysis was used to assess the association between smoking and risk of CRC subtypes by molecular characteristics, adjusting for age, sex, and study. All statistical tests were 2-sided and adjusted for Bonferroni correction. Results: Heavier smoking was associated with higher risk of CRC overall and stratified by individual markers (P trend < .001). The associations differed statistically significantly between all molecular subtypes, which was the most statistically significant for CIMP and BRAF. Compared with never-smokers, smokers in the fourth quartile of pack-years had a 90% higher risk of CIMP-positive CRC (odds ratio = 1.90, 95% confidence interval = 1.60 to 2.26) but only 35% higher risk for CIMP-negative CRC (odds ratio = 1.35, 95% confidence interval = 1.22 to 1.49; P difference = 2.1 x 10-6). The association was also stronger in tumors that were CIMP positive, MSI high, or KRAS wild type when combined (P difference < .001). Conclusion: Smoking was associated with differential risk of CRC subtypes defined by molecular characteristics. Heavier smokers had particularly higher risk of CRC subtypes that were CIMP positive and MSI high in combination, suggesting that smoking may be involved in the development of colorectal tumors via the serrated pathway.

Telomere Maintenance Variants and Survival after Colorectal Cancer: Smoking- and Sex-Specific Associations.

BACKGROUND: Telomeres play an important role in colorectal cancer prognosis. Variation in telomere maintenance genes may be associated with survival after colorectal cancer diagnosis, but evidence is limited. In addition, possible interactions between telomere maintenance genes and prognostic factors, such as smoking and sex, also remain to be investigated. METHODS: We conducted gene-wide analyses of colorectal cancer prognosis in 4,896 invasive colorectal cancer cases from the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO); 1,871 common variants within 13 telomere maintenance genes were included. Cox models were fit to estimate associations of these variants individually with overall and colorectal cancer-specific survival. Likelihood ratio tests were used to test for interaction by smoking and sex. P values were adjusted using Bonferroni correction. RESULTS: The association between minor allele of rs7200950 (ACD) with colorectal cancer-specific survival varied significantly by smoking pack-years (corrected P = 0.049), but no significant trend was observed. By sex, minor alleles for rs2975843 (TERF1), rs75676021 (POT1), and rs74429678 (POT1) were associated with decreased overall and/or colorectal cancer-specific survival in women but not in men. CONCLUSIONS: Our study reported a gene-wide statistically significant interaction with sex (TERF1, POT1). Although significant interaction by smoking pack-years (ACD) was observed, there was no evidence of a dose response. Validation of these findings in other large studies and further functional annotation on these SNPs are warranted. IMPACT: Our study found a gene-smoking and gene-sex interaction on survival after colorectal cancer diagnosis, providing new insights into the role of genetic polymorphisms in telomere maintenance on colorectal cancer prognosis.

Genetic Predictors of Circulating 25-Hydroxyvitamin D and Prognosis after Colorectal Cancer.

BACKGROUND: Low serum 25-hydroxyvitamin D [25(OH)D] concentrations in patients with colorectal cancer have been consistently associated with higher mortality in observational studies. It is unclear whether low 25(OH)D levels directly influence colorectal cancer mortality. To minimize bias, we use genetic variants associated with vitamin D levels to evaluate the association with overall and colorectal cancer-specific survival. METHODS: Six genetic variants have been robustly identified to be associated with 25(OH)D levels in genome-wide association studies. On the basis of data from the International Survival Analysis in Colorectal Cancer Consortium, the individual genetic variants and a weighted genetic risk score were tested for association with overall and colorectal cancer-specific survival using Cox proportional hazards models in 7,657 patients with stage I to IV colorectal cancer, of whom 2,438 died from any cause and 1,648 died from colorectal cancer. RESULTS: The 25(OH)D decreasing allele of SNP rs2282679 (GC gene, encodes group-specific component/vitamin D-binding protein) was associated with poorer colorectal cancer-specific survival, although not significant after multiple-testing correction. None of the other five SNPs showed an association. The genetic risk score showed nonsignificant associations with increased overall [HR = 1.54; confidence interval (CI), 0.86-2.78] and colorectal cancer-specific mortality (HR = 1.76; 95% CI, 0.86-3.58). A significant increased risk of overall mortality was observed in women (HR = 3.26; 95% CI, 1.45-7.33; P heterogeneity = 0.01) and normal-weight individuals (HR = 4.14; 95% CI, 1.50-11.43, P heterogeneity = 0.02). CONCLUSIONS: Our results provided little evidence for an association of genetic predisposition of lower vitamin D levels with increased overall or colorectal cancer-specific survival, although power might have been an issue. IMPACT: Further studies are warranted to investigate the association in specific subgroups.

Mendelian Randomization of Circulating Polyunsaturated Fatty Acids and Colorectal Cancer Risk.

BACKGROUND: Results from epidemiologic studies examining polyunsaturated fatty acids (PUFA) and colorectal cancer risk are inconsistent. Mendelian randomization may strengthen causal inference from observational studies. Given their shared metabolic pathway, examining the combined effects of aspirin/NSAID use with PUFAs could help elucidate an association between PUFAs and colorectal cancer risk. METHODS: Information was leveraged from genome-wide association studies (GWAS) regarding PUFA-associated SNPs to create weighted genetic scores (wGS) representing genetically predicted circulating blood PUFAs for 11,016 non-Hispanic white colorectal cancer cases and 13,732 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). Associations per SD increase in the wGS were estimated using unconditional logistic regression. Interactions between PUFA wGSs and aspirin/NSAID use on colorectal cancer risk were also examined. RESULTS: Modest colorectal cancer risk reductions were observed per SD increase in circulating linoleic acid [ORLA = 0.96; 95% confidence interval (CI) = 0.93-0.98; P = 5.2 × 10-4] and α-linolenic acid (ORALA = 0.95; 95% CI = 0.92-0.97; P = 5.4 × 10-5), whereas modest increased risks were observed for arachidonic (ORAA = 1.06; 95% CI = 1.03-1.08; P = 3.3 × 10-5), eicosapentaenoic (OREPA = 1.04; 95% CI = 1.01-1.07; P = 2.5 × 10-3), and docosapentaenoic acids (ORDPA = 1.03; 95% CI = 1.01-1.06; P = 1.2 × 10-2). Each of these effects was stronger among aspirin/NSAID nonusers in the stratified analyses. CONCLUSIONS: Our study suggests that higher circulating shorter-chain PUFAs (i.e., LA and ALA) were associated with reduced colorectal cancer risk, whereas longer-chain PUFAs (i.e., AA, EPA, and DPA) were associated with an increased colorectal cancer risk. IMPACT: The interaction of PUFAs with aspirin/NSAID use indicates a shared colorectal cancer inflammatory pathway. Future research should continue to improve PUFA genetic instruments to elucidate the independent effects of PUFAs on colorectal cancer.

Landscape of somatic single nucleotide variants and indels in colorectal cancer and impact on survival.

Colorectal cancer (CRC) is a biologically heterogeneous disease. To characterize its mutational profile, we conduct targeted sequencing of 205 genes for 2,105 CRC cases with survival data. Our data shows several findings in addition to enhancing the existing knowledge of CRC. We identify PRKCI, SPZ1, MUTYH, MAP2K4, FETUB, and TGFBR2 as additional genes significantly mutated in CRC. We find that among hypermutated tumors, an increased mutation burden is associated with improved CRC-specific survival (HR = 0.42, 95% CI: 0.21-0.82). Mutations in TP53 are associated with poorer CRC-specific survival, which is most pronounced in cases carrying TP53 mutations with predicted 0% transcriptional activity (HR = 1.53, 95% CI: 1.21-1.94). Furthermore, we observe differences in mutational frequency of several genes and pathways by tumor location, stage, and sex. Overall, this large study provides deep insights into somatic mutations in CRC, and their potential relationships with survival and tumor features.

Does evolutionary innovation in pharyngeal jaws lead to rapid lineage diversification in labrid fishes?

BACKGROUND: Major modifications to the pharyngeal jaw apparatus are widely regarded as a recurring evolutionary key innovation that has enabled adaptive radiation in many species-rich clades of percomorph fishes. However one of the central predictions of this hypothesis, that the acquisition of a modified pharyngeal jaw apparatus will be positively correlated with explosive lineage diversification, has never been tested. We applied comparative methods to a new time-calibrated phylogeny of labrid fishes to test whether diversification rates shifted at two scales where major pharyngeal jaw innovations have evolved: across all of Labridae and within the subclade of parrotfishes. RESULTS: Diversification patterns within early labrids did not reflect rapid initial radiation. Much of modern labrid diversity stems from two recent rapid diversification events; one within julidine fishes and the other with the origin of the most species-rich clade of reef-associated parrotfishes. A secondary pharyngeal jaw innovation was correlated with rapid diversification within the parrotfishes. However diversification rate shifts within parrotfishes are more strongly correlated with the evolution of extreme dichromatism than with pharyngeal jaw modifications. CONCLUSION: The temporal lag between pharyngeal jaw modifications and changes in diversification rates casts doubt on the key innovation hypothesis as a simple explanation for much of the richness seen in labrids and scarines. Although the possession of a secondarily modified PJA was correlated with increased diversification rates, this pattern is better explained by the evolution of extreme dichromatism (and other social and behavioral characters relating to sexual selection) within Scarus and Chlorurus. The PJA-innovation hypothesis also fails to explain the most dominant aspect of labrid lineage diversification, the radiation of the julidines. We suggest that pharyngeal jaws might have played a more important role in enabling morphological evolution of the feeding apparatus in labrids and scarines rather than in accelerating lineage diversification.

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.