Germline genetic regulation of the colorectal tumor immune microenvironment.

OBJECTIVE: To evaluate the contribution of germline genetics to regulating the briskness and diversity of T cell responses in CRC, we conducted a genome-wide association study to examine the associations between germline genetic variation and quantitative measures of T cell landscapes in 2,876 colorectal tumors from participants in the Molecular Epidemiology of Colorectal Cancer Study (MECC). METHODS: Germline DNA samples were genotyped and imputed using genome-wide arrays. Tumor DNA samples were extracted from paraffin blocks, and T cell receptor clonality and abundance were quantified by immunoSEQ (Adaptive Biotechnologies, Seattle, WA). Tumor infiltrating lymphocytes per high powered field (TILs/hpf) were scored by a gastrointestinal pathologist. Regression models were used to evaluate the associations between each variant and the three T-cell features, adjusting for sex, age, genotyping platform, and global ancestry. Three independent datasets were used for replication. RESULTS: We identified a SNP (rs4918567) near RBM20 associated with clonality at a genome-wide significant threshold of 5 × 10- 8, with a consistent direction of association in both discovery and replication datasets. Expression quantitative trait (eQTL) analyses and in silico functional annotation for these loci provided insights into potential functional roles, including a statistically significant eQTL between the T allele at rs4918567 and higher expression of ADRA2A (P = 0.012) in healthy colon mucosa. CONCLUSIONS: Our study suggests that germline genetic variation is associated with the quantity and diversity of adaptive immune responses in CRC. Further studies are warranted to replicate these findings in additional samples and to investigate functional genomic mechanisms.

Genetic risk impacts the association of menopausal hormone therapy with colorectal cancer risk.

BACKGROUND: Menopausal hormone therapy (MHT), a common treatment to relieve symptoms of menopause, is associated with a lower risk of colorectal cancer (CRC). To inform CRC risk prediction and MHT risk-benefit assessment, we aimed to evaluate the joint association of a polygenic risk score (PRS) for CRC and MHT on CRC risk. METHODS: We used data from 28,486 postmenopausal women (11,519 cases and 16,967 controls) of European descent. A PRS based on 141 CRC-associated genetic variants was modeled as a categorical variable in quartiles. Multiplicative interaction between PRS and MHT use was evaluated using logistic regression. Additive interaction was measured using the relative excess risk due to interaction (RERI). 30-year cumulative risks of CRC for 50-year-old women according to MHT use and PRS were calculated. RESULTS: The reduction in odds ratios by MHT use was larger in women within the highest quartile of PRS compared to that in women within the lowest quartile of PRS (p-value = 2.7 × 10-8). At the highest quartile of PRS, the 30-year CRC risk was statistically significantly lower for women taking any MHT than for women not taking any MHT, 3.7% (3.3%-4.0%) vs 6.1% (5.7%-6.5%) (difference 2.4%, P-value = 1.83 × 10-14); these differences were also statistically significant but smaller in magnitude in the lowest PRS quartile, 1.6% (1.4%-1.8%) vs 2.2% (1.9%-2.4%) (difference 0.6%, P-value = 1.01 × 10-3), indicating 4 times greater reduction in absolute risk associated with any MHT use in the highest compared to the lowest quartile of genetic CRC risk. CONCLUSIONS: MHT use has a greater impact on the reduction of CRC risk for women at higher genetic risk. These findings have implications for the development of risk prediction models for CRC and potentially for the consideration of genetic information in the risk-benefit assessment of MHT use.

Characterization of additive gene-environment interactions for colorectal cancer risk.

BACKGROUND: Colorectal cancer (CRC) is a common, fatal cancer. Identifying subgroups who may benefit more from intervention is of critical public health importance. Previous studies have assessed multiplicative interaction between genetic risk scores and environmental factors, but few have assessed additive interaction, the relevant public health measure. METHODS: Using resources from colorectal cancer consortia including 45,247 CRC cases and 52,671 controls, we assessed multiplicative and additive interaction (relative excess risk due to interaction, RERI) using logistic regression between 13 harmonized environmental factors and genetic risk score including 141 variants associated with CRC risk. RESULTS: There was no evidence of multiplicative interaction between environmental factors and genetic risk score. There was additive interaction where, for individuals with high genetic susceptibility, either heavy drinking [RERI = 0.24, 95% confidence interval, CI, (0.13, 0.36)], ever smoking [0.11 (0.05, 0.16)], high BMI [female 0.09 (0.05, 0.13), male 0.10 (0.05, 0.14)], or high red meat intake [highest versus lowest quartile 0.18 (0.09, 0.27)] was associated with excess CRC risk greater than that for individuals with average genetic susceptibility. Conversely, we estimate those with high genetic susceptibility may benefit more from reducing CRC risk with aspirin/NSAID use [-0.16 (-0.20, -0.11)] or higher intake of fruit, fiber, or calcium [highest quartile versus lowest quartile -0.12 (-0.18, -0.050); -0.16 (-0.23, -0.09); -0.11 (-0.18, -0.05), respectively] than those with average genetic susceptibility. CONCLUSIONS: Additive interaction is important to assess for identifying subgroups who may benefit from intervention. The subgroups identified in this study may help inform precision CRC prevention.

Statistical methods for Mendelian models with multiple genes and cancers.

Risk evaluation to identify individuals who are at greater risk of cancer as a result of heritable pathogenic variants is a valuable component of individualized clinical management. Using principles of Mendelian genetics, Bayesian probability theory, and variant-specific knowledge, Mendelian models derive the probability of carrying a pathogenic variant and developing cancer in the future, based on family history. Existing Mendelian models are widely employed, but are generally limited to specific genes and syndromes. However, the upsurge of multigene panel germline testing has spurred the discovery of many new gene-cancer associations that are not presently accounted for in these models. We have developed PanelPRO, a flexible, efficient Mendelian risk prediction framework that can incorporate an arbitrary number of genes and cancers, overcoming the computational challenges that arise because of the increased model complexity. We implement an 11-gene, 11-cancer model, the largest Mendelian model created thus far, based on this framework. Using simulations and a clinical cohort with germline panel testing data, we evaluate model performance, validate the reverse-compatibility of our approach with existing Mendelian models, and illustrate its usage. Our implementation is freely available for research use in the PanelPRO R package.

Probing the diabetes and colorectal cancer relationship using gene - environment interaction analyses.

BACKGROUND: Diabetes is an established risk factor for colorectal cancer. However, the mechanisms underlying this relationship still require investigation and it is not known if the association is modified by genetic variants. To address these questions, we undertook a genome-wide gene-environment interaction analysis. METHODS: We used data from 3 genetic consortia (CCFR, CORECT, GECCO; 31,318 colorectal cancer cases/41,499 controls) and undertook genome-wide gene-environment interaction analyses with colorectal cancer risk, including interaction tests of genetics(G)xdiabetes (1-degree of freedom; d.f.) and joint testing of Gxdiabetes, G-colorectal cancer association (2-d.f. joint test) and G-diabetes correlation (3-d.f. joint test). RESULTS: Based on the joint tests, we found that the association of diabetes with colorectal cancer risk is modified by loci on chromosomes 8q24.11 (rs3802177, SLC30A8 - ORAA: 1.62, 95% CI: 1.34-1.96; ORAG: 1.41, 95% CI: 1.30-1.54; ORGG: 1.22, 95% CI: 1.13-1.31; p-value3-d.f.: 5.46 × 10-11) and 13q14.13 (rs9526201, LRCH1 - ORGG: 2.11, 95% CI: 1.56-2.83; ORGA: 1.52, 95% CI: 1.38-1.68; ORAA: 1.13, 95% CI: 1.06-1.21; p-value2-d.f.: 7.84 × 10-09). DISCUSSION: These results suggest that variation in genes related to insulin signaling (SLC30A8) and immune function (LRCH1) may modify the association of diabetes with colorectal cancer risk and provide novel insights into the biology underlying the diabetes and colorectal cancer relationship.

Clinical implications of conflicting variant interpretations in the cancer genetics clinic.

PURPOSE: The aim of this study was to describe the clinical impact of commercial laboratories issuing conflicting classifications of genetic variants. METHODS: Results from 2000 patients undergoing a multigene hereditary cancer panel by a single laboratory were analyzed. Clinically significant discrepancies between the laboratory-provided test reports and other major commercial laboratories were identified, including differences between pathogenic/likely pathogenic and variant of uncertain significance (VUS) classifications, via review of ClinVar archives. For patients carrying a VUS, clinical documentation was assessed for evidence of provider awareness of the conflict. RESULTS: Fifty of 975 (5.1%) patients with non-negative results carried a variant with a clinically significant conflict, 19 with a pathogenic/likely pathogenic variant reported in APC or MUTYH, and 31 with a VUS reported in CDKN2A, CHEK2, MLH1, MSH2, MUTYH, RAD51C, or TP53. Only 10 of 28 (36%) patients with a VUS with a clinically significant conflict had a documented discussion by a provider about the conflict. Discrepant counseling strategies were used for different patients with the same variant. Among patients with a CDKN2A variant or a monoallelic MUTYH variant, providers were significantly more likely to make recommendations based on the laboratory-reported classification. CONCLUSION: Our findings highlight the frequency of variant interpretation discrepancies and importance of clinician awareness. Guidance is needed on managing patients with discrepant variants to support accurate risk assessment.

A general framework for functionally informed set-based analysis: Application to a large-scale colorectal cancer study.

Genome-wide association studies (GWAS) have successfully identified tens of thousands of genetic variants associated with various phenotypes, but together they explain only a fraction of heritability, suggesting many variants have yet to be discovered. Recently it has been recognized that incorporating functional information of genetic variants can improve power for identifying novel loci. For example, S-PrediXcan and TWAS tested the association of predicted gene expression with phenotypes based on GWAS summary statistics by leveraging the information on genetic regulation of gene expression and found many novel loci. However, as genetic variants may have effects on more than one gene and through different mechanisms, these methods likely only capture part of the total effects of these variants. In this paper, we propose a summary statistics-based mixed effects score test (sMiST) that tests for the total effect of both the effect of the mediator by imputing genetically predicted gene expression, like S-PrediXcan and TWAS, and the direct effects of individual variants. It allows for multiple functional annotations and multiple genetically predicted mediators. It can also perform conditional association analysis while adjusting for other genetic variants (e.g., known loci for the phenotype). Extensive simulation and real data analyses demonstrate that sMiST yields p-values that agree well with those obtained from individual level data but with substantively improved computational speed. Importantly, a broad application of sMiST to GWAS is possible, as only summary statistics of genetic variant associations are required. We apply sMiST to a large-scale GWAS of colorectal cancer using summary statistics from ∼120, 000 study participants and gene expression data from the Genotype-Tissue Expression (GTEx) project. We identify several novel and secondary independent genetic loci.

Exploring various polygenic risk scores for skin cancer in the phenomes of the Michigan genomics initiative and the UK Biobank with a visual catalog: PRSWeb.

Polygenic risk scores (PRS) are designed to serve as single summary measures that are easy to construct, condensing information from a large number of genetic variants associated with a disease. They have been used for stratification and prediction of disease risk. The primary focus of this paper is to demonstrate how we can combine PRS and electronic health records data to better understand the shared and unique genetic architecture and etiology of disease subtypes that may be both related and heterogeneous. PRS construction strategies often depend on the purpose of the study, the available data/summary estimates, and the underlying genetic architecture of a disease. We consider several choices for constructing a PRS using data obtained from various publicly-available sources including the UK Biobank and evaluate their abilities to predict not just the primary phenotype but also secondary phenotypes derived from electronic health records (EHR). This study was conducted using data from 30,702 unrelated, genotyped patients of recent European descent from the Michigan Genomics Initiative (MGI), a longitudinal biorepository effort within Michigan Medicine. We examine the three most common skin cancer subtypes in the USA: basal cell carcinoma, cutaneous squamous cell carcinoma, and melanoma. Using these PRS for various skin cancer subtypes, we conduct a phenome-wide association study (PheWAS) within the MGI data to evaluate PRS associations with secondary traits. PheWAS results are then replicated using population-based UK Biobank data and compared across various PRS construction methods. We develop an accompanying visual catalog called PRSweb that provides detailed PheWAS results and allows users to directly compare different PRS construction methods.

Psychosocial outcomes following germline multigene panel testing in an ethnically and economically diverse cohort of patients.

BACKGROUND: Little is known about the psychological outcomes of germline multigene panel testing, particularly among diverse patients and those with moderate-risk pathogenic variants (PVs). METHODS: Study participants (N = 1264) were counseled and tested with a 25- or 28-gene panel and completed a 3-month postresult survey including the Multidimensional Impact of Cancer Risk Assessment (MICRA). RESULTS: The mean age was 52 years, 80% were female, and 70% had cancer; 45% were non-Hispanic White, 37% were Hispanic, 10% were Asian, 3% were Black, and 5% had another race/ethnicity. Approximately 28% had a high school education or less, and 23% were non-English-speaking. The genetic test results were as follows: 7% had a high-risk PV, 6% had a moderate-risk PV, 35% had a variant of uncertain significance (VUS), and 52% were negative. Most participants (92%) had a total MICRA score ≤ 38, which corresponded to a mean response of "never," "rarely," or only "sometimes" reacting negatively to results. A multivariate analysis found that mean total MICRA scores were significantly higher (more uncertainty/distress) among high- and moderate-risk PV carriers (29.7 and 24.8, respectively) than those with a VUS or negative results (17.4 and 16.1, respectively). Having cancer or less education was associated with a significantly higher total MICRA score; race/ethnicity was not associated with the total MICRA score. High- and moderate-risk PV carriers did not differ significantly from one another in the total MICRA score, uncertainty, distress, or positive experiences. CONCLUSIONS: In a diverse population undergoing genetic counseling and multigene panel testing for hereditary cancer risk, the psychological response corresponded to test results and showed low distress and uncertainty. Further studies are needed to assess patient understanding and subsequent cancer screening among patients from diverse backgrounds. LAY SUMMARY: Multigene panel tests for hereditary cancer have become widespread despite concerns about adverse psychological reactions among carriers of moderate-risk pathogenic variants (mutations) and among carriers of variants of uncertain significance. This large study of an ethnically and economically diverse cohort of patients undergoing panel testing found that 92% "never," "rarely," or only "sometimes" reacted negatively to results. Somewhat higher uncertainty and distress were identified among carriers of high- and moderate-risk pathogenic variants, and lower levels were identified among those with a variant of uncertain significance or a negative result. Although the psychological response corresponded to risk, reactions to testing were favorable, regardless of results.

Association of Polygenic Risk Scores for Multiple Cancers in a Phenome-wide Study: Results from The Michigan Genomics Initiative.

Health systems are stewards of patient electronic health record (EHR) data with extraordinarily rich depth and breadth, reflecting thousands of diagnoses and exposures. Measures of genomic variation integrated with EHRs offer a potential strategy to accurately stratify patients for risk profiling and discover new relationships between diagnoses and genomes. The objective of this study was to evaluate whether polygenic risk scores (PRS) for common cancers are associated with multiple phenotypes in a phenome-wide association study (PheWAS) conducted in 28,260 unrelated, genotyped patients of recent European ancestry who consented to participate in the Michigan Genomics Initiative, a longitudinal biorepository effort within Michigan Medicine. PRS for 12 cancer traits were calculated using summary statistics from the NHGRI-EBI catalog. A total of 1,711 synthetic case-control studies was used for PheWAS analyses. There were 13,490 (47.7%) patients with at least one cancer diagnosis in this study sample. PRS exhibited strong association for several cancer traits they were designed for, including female breast cancer, prostate cancer, melanoma, basal cell carcinoma, squamous cell carcinoma, and thyroid cancer. Phenome-wide significant associations were observed between PRS and many non-cancer diagnoses. To differentiate PRS associations driven by the primary trait from associations arising through shared genetic risk profiles, the idea of "exclusion PRS PheWAS" was introduced. Further analysis of temporal order of the diagnoses improved our understanding of these secondary associations. This comprehensive PheWAS used PRS instead of a single variant.

Circulating Levels of Insulin-like Growth Factor 1 and Insulin-like Growth Factor Binding Protein 3 Associate With Risk of Colorectal Cancer Based on Serologic and Mendelian Randomization Analyses.

BACKGROUND & AIMS: Human studies examining associations between circulating levels of insulin-like growth factor 1 (IGF1) and insulin-like growth factor binding protein 3 (IGFBP3) and colorectal cancer risk have reported inconsistent results. We conducted complementary serologic and Mendelian randomization (MR) analyses to determine whether alterations in circulating levels of IGF1 or IGFBP3 are associated with colorectal cancer development. METHODS: Serum levels of IGF1 were measured in blood samples collected from 397,380 participants from the UK Biobank, from 2006 through 2010. Incident cancer cases and cancer cases recorded first in death certificates were identified through linkage to national cancer and death registries. Complete follow-up was available through March 31, 2016. For the MR analyses, we identified genetic variants associated with circulating levels of IGF1 and IGFBP3. The association of these genetic variants with colorectal cancer was examined with 2-sample MR methods using genome-wide association study consortia data (52,865 cases with colorectal cancer and 46,287 individuals without [controls]) RESULTS: After a median follow-up period of 7.1 years, 2665 cases of colorectal cancer were recorded. In a multivariable-adjusted model, circulating level of IGF1 associated with colorectal cancer risk (hazard ratio per 1 standard deviation increment of IGF1, 1.11; 95% confidence interval [CI] 1.05-1.17). Similar associations were found by sex, follow-up time, and tumor subsite. In the MR analyses, a 1 standard deviation increment in IGF1 level, predicted based on genetic factors, was associated with a higher risk of colorectal cancer risk (odds ratio 1.08; 95% CI 1.03-1.12; P = 3.3 × 10-4). Level of IGFBP3, predicted based on genetic factors, was associated with colorectal cancer risk (odds ratio per 1 standard deviation increment, 1.12; 95% CI 1.06-1.18; P = 4.2 × 10-5). Colorectal cancer risk was associated with only 1 variant in the IGFBP3 gene region (rs11977526), which also associated with anthropometric traits and circulating level of IGF2. CONCLUSIONS: In an analysis of blood samples from almost 400,000 participants in the UK Biobank, we found an association between circulating level of IGF1 and colorectal cancer. Using genetic data from 52,865 cases with colorectal cancer and 46,287 controls, a higher level of IGF1, determined by genetic factors, was associated with colorectal cancer. Further studies are needed to determine how this signaling pathway might contribute to colorectal carcinogenesis.

A likelihood-based approach to assessing frequency of pathogenicity among variants of unknown significance in susceptibility genes.

Commercialized multigene panel testing brings unprecedented opportunities to understand germline genetic contributions to hereditary cancers. Most genetic testing companies classify the pathogenicity of variants as pathogenic, benign, or variants of unknown significance (VUSs). The unknown pathogenicity of VUSs poses serious challenges to clinical decision-making. This study aims to assess the frequency of VUSs that are likely pathogenic in disease-susceptibility genes. Using estimates of probands' probability of having a pathogenic mutation (ie, the carrier score) based on a family history probabilistic risk prediction model, we assume the carrier score distribution for probands with VUSs is a mixture of the carrier score distribution for probands with positive results and the carrier score distribution for probands with negative results. Under this mixture model, we propose a likelihood-based approach to assess the frequency of pathogenicity among probands with VUSs, while accounting for the existence of possible pathogenic mutations on genes not tested. We conducted simulations to assess the performance of the approach and show that under various settings, the approach performs well with very little bias in the estimated proportion of VUSs that are likely pathogenic. We also estimate the positive predictive value across the entire range of carrier scores. We apply our approach to the USC-Stanford Hereditary Cancer Panel Testing cohort, and estimate the proportion of probands that have VUSs in BRCA1/2 that are likely pathogenic to be 10.12% [95%CI: 0%, 43.04%]. This approach will enable clinicians to target high-risk patients who have VUSs, allowing for early prevention interventions.

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.

Large-Scale Genome-Wide Association Study of East Asians Identifies Loci Associated With Risk for Colorectal Cancer.

BACKGROUND & AIMS: Genome-wide association studies (GWASs) have associated approximately 50 loci with risk of colorectal cancer (CRC)-nearly one third of these loci were initially associated with CRC in studies conducted in East Asian populations. We conducted a GWAS of East Asians to identify CRC risk loci and evaluate the generalizability of findings from GWASs of European populations to Asian populations. METHODS: We analyzed genetic data from 22,775 patients with CRC (cases) and 47,731 individuals without cancer (controls) from 14 studies in the Asia Colorectal Cancer Consortium. First, we performed a meta-analysis of 7 GWASs (10,625 cases and 34,595 controls) and identified 46,554 promising risk variants for replication by adding them to the Multi-Ethnic Global Array (MEGA) for genotype analysis in 6445 cases and 7175 controls. These data were analyzed, along with data from an additional 5705 cases and 5961 controls genotyped using the OncoArray. We also obtained data from 57,976 cases and 67,242 controls of European descent. Variants at identified risk loci were functionally annotated and evaluated in correlation with gene expression levels. RESULTS: A meta-analyses of all samples from people of Asian descent identified 13 loci and 1 new variant at a known locus (10q24.2) associated with risk of CRC at the genome-wide significance level of P < 5 × 10-8. We did not perform experiments to replicate these associations in additional individuals of Asian ancestry. However, the lead risk variant in 6 of these loci was also significantly associated with risk of CRC in European descendants. A strong association (44%-75% increase in risk per allele) was found for 2 low-frequency variants: rs201395236 at 1q44 (minor allele frequency, 1.34%) and rs77969132 at 12p11.21 (minor allele frequency, 1.53%). For 8 of the 13 associated loci, the variants with the highest levels of significant association were located inside or near the protein-coding genes L1TD1, EFCAB2, PPP1R21, SLCO2A1, HLA-G, NOTCH4, DENND5B, and GNAS. For other intergenic loci, we provided evidence for the possible involvement of the genes ALDH7A1, PRICKLE1, KLF5, WWOX, and GLP2R. We replicated findings for 41 of 52 previously reported risk loci. CONCLUSIONS: We showed that most of the risk loci previously associated with CRC risk in individuals of European descent were also associated with CRC risk in East Asians. Furthermore, we identified 13 loci significantly associated with risk for CRC in Asians. Many of these loci contained genes that regulate the immune response, Wnt signaling to ß-catenin, prostaglandin E2 catabolism, and cell pluripotency and proliferation. Further analyses of these genes and their variants is warranted, particularly for the 8 loci for which the lead CRC risk variants were not replicated in persons of European descent.

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.

Inhibition of poly(ADP-ribose) polymerase induces synthetic lethality in BRIP1 deficient ovarian epithelial cells.

OBJECTIVE: Pathogenic variations in the homologous recombination (HR) gene, BRCA1 interacting protein C-terminal helicase 1 (BRIP1) increase the risk for ovarian cancer. PARP inhibitors (PARPi) exert a synthetic lethal effect in BRCA-mutated ovarian cancers. Effective HR requires cooperation between BRCA1 and BRIP1; therefore, BRIP1-incompetancy may predict vulnerability to synthetic lethality. Here we investigated the response of ovarian epithelial cells with defective BRIP1 function to PARPi, and compared these cells to those lacking BRCA1 activity. METHODS: We engineered Chinese Hamster ovarian (CHO) epithelial cells to express deficient BRIP1 or BRCA1, and exposed them to olaparib with or without carboplatin or cisplatin. We assessed cellular proliferation and survival; we calculated inhibitory concentrations and combination and reduction drug indices. RESULTS: BRIP1 and BRCA1 inactivation impedes HR activity, decreases cellular proliferation and compromises DNA damage recovery. Platinum agent exposure impairs cellular survival. Olaparib exposure alone decreases cell viability in BRCA1-deficient cells, although has no effect on BRIP1-deficient cells. Combining carboplatin or cisplatin with olaparib synergistically attenuates cellular survival, consistent with synthetic lethality. CONCLUSIONS: BRIP1-deficient ovarian epithelial cells exhibit defective HR, resulting in synthetic lethality when exposed to a platinum agent/PARPi combination. PARPi alone had no effect; this lack of effect may result from distinguishing molecular properties of BRIP1and/or consequences of genomic background. Our study identifies altered BRIP1 as a target for precision medicine-based therapies for ovarian cancers. This investigation supports consideration of the use of a platinum agent/PARPi combination in ovarian cancers depending upon genetic profile and genomic background.

Germline Genetic Features of Young Individuals With Colorectal Cancer.

BACKGROUND & AIMS: The incidence of colorectal cancer (CRC) in individuals younger than 50 years is increasing. We sought to ascertain the proportion of young CRC cases associated with genetic predisposition. METHODS: We performed a retrospective study of individuals diagnosed with CRC at an age younger than 50 years, evaluated by the clinical genetics service at a single tertiary care cancer center from 1998 through 2015. We collected data on patient histories, tumor phenotypes, and results of germline DNA sequencing. For subjects with uninformative clinical evaluations, germline DNA samples were (re)sequenced using a research-based next-generation sequencing multigene panel. The primary outcome was identification of a pathogenic germline mutation associated with cancer predisposition. RESULTS: Of 430 young CRC cases, 111 (26%) had a first-degree relative with CRC. Forty-one of the subjects with CRC (10%) had tumors with histologic evidence for mismatch repair deficiency. Of 315 subjects who underwent clinical germline sequencing, 79 had mutations associated with a hereditary cancer syndrome and 21 had variants of uncertain significance. Fifty-six subjects had pathogenic variants associated with Lynch syndrome (25 with mutations in MSH2, 24 with mutations in MLH1, 5 with mutations in MSH6, and 2 with mutations in PMS2) and 10 subjects had pathogenic variants associated with familial adenomatous polyposis. Thirteen subjects had mutations in other cancer-associated genes (8 in MUTYH, 2 in SMAD4, 1 in BRCA1, 1 in TP53, and 1 in CHEK2), all identified through multigene panel tests. Among 117 patients with uninformative clinical evaluations, next-generation sequence analysis using a multigene panel detected actionable germline variants in 6 patients (5%). Only 43 of the 85 subjects with germline mutations associated with a hereditary cancer syndrome (51%) reported a CRC diagnosis in a first-degree relative. CONCLUSIONS: Approximately 1 in 5 individuals diagnosed with CRC at age younger than 50 years carries a germline mutation associated with cancer; nearly half of these do not have clinical histories typically associated with the identified syndrome. Germline testing with multigene cancer panels should be considered for all young patients with CRC.

Determining Risk of Colorectal Cancer and Starting Age of Screening Based on Lifestyle, Environmental, and Genetic Factors.

BACKGROUND & AIMS: Guidelines for initiating colorectal cancer (CRC) screening are based on family history but do not consider lifestyle, environmental, or genetic risk factors. We developed models to determine risk of CRC, based on lifestyle and environmental factors and genetic variants, and to identify an optimal age to begin screening. METHODS: We collected data from 9748 CRC cases and 10,590 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium and the Colorectal Transdisciplinary study, from 1992 through 2005. Half of the participants were used to develop the risk determination model and the other half were used to evaluate the discriminatory accuracy (validation set). Models of CRC risk were created based on family history, 19 lifestyle and environmental factors (E-score), and 63 CRC-associated single-nucleotide polymorphisms identified in genome-wide association studies (G-score). We evaluated the discriminatory accuracy of the models by calculating area under the receiver operating characteristic curve values, adjusting for study, age, and endoscopy history for the validation set. We used the models to project the 10-year absolute risk of CRC for a given risk profile and recommend ages to begin screening in comparison to CRC risk for an average individual at 50 years of age, using external population incidence rates for non-Hispanic whites from the Surveillance, Epidemiology, and End Results program registry. RESULTS: In our models, E-score and G-score each determined risk of CRC with greater accuracy than family history. A model that combined both scores and family history estimated CRC risk with an area under the receiver operating characteristic curve value of 0.63 (95% confidence interval, 0.62-0.64) for men and 0.62 (95% confidence interval, 0.61-0.63) for women; area under the receiver operating characteristic curve values based on only family history ranged from 0.53 to 0.54 and those based only E-score or G-score ranged from 0.59 to 0.60. Although screening is recommended to begin at age 50 years for individuals with no family history of CRC, starting ages calculated based on combined E-score and G-score differed by 12 years for men and 14 for women, for individuals with the highest vs the lowest 10% of risk. CONCLUSIONS: We used data from 2 large international consortia to develop CRC risk calculation models that included genetic and environmental factors along with family history. These determine risk of CRC and starting ages for screening with greater accuracy than the family history only model, which is based on the current screening guideline. These scoring systems might serve as a first step toward developing individualized CRC prevention strategies.

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.

Inherited variation in circadian rhythm genes and risks of prostate cancer and three other cancer sites in combined cancer consortia.

Circadian disruption has been linked to carcinogenesis in animal models, but the evidence in humans is inconclusive. Genetic variation in circadian rhythm genes provides a tool to investigate such associations. We examined associations of genetic variation in nine core circadian rhythm genes and six melatonin pathway genes with risk of colorectal, lung, ovarian and prostate cancers using data from the Genetic Associations and Mechanisms in Oncology (GAME-ON) network. The major results for prostate cancer were replicated in the Prostate, Lung, Colorectal and Ovarian (PLCO) cancer screening trial, and for colorectal cancer in the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). The total number of cancer cases and controls was 15,838/18,159 for colorectal, 14,818/14,227 for prostate, 12,537/17,285 for lung and 4,369/9,123 for ovary. For each cancer site, we conducted gene-based and pathway-based analyses by applying the summary-based Adaptive Rank Truncated Product method (sARTP) on the summary association statistics for each SNP within the candidate gene regions. Aggregate genetic variation in circadian rhythm and melatonin pathways were significantly associated with the risk of prostate cancer in data combining GAME-ON and PLCO, after Bonferroni correction (ppathway < 0.00625). The two most significant genes were NPAS2 (pgene = 0.0062) and AANAT (pgene = 0.00078); the latter being significant after Bonferroni correction. For colorectal cancer, we observed a suggestive association with the circadian rhythm pathway in GAME-ON (ppathway = 0.021); this association was not confirmed in GECCO (ppathway = 0.76) or the combined data (ppathway = 0.17). No significant association was observed for ovarian and lung cancer. These findings support a potential role for circadian rhythm and melatonin pathways in prostate carcinogenesis. Further functional studies are needed to better understand the underlying biologic mechanisms.