Quantitation of Ten Urinary Nicotine Metabolites, Including 4-Hydroxy-4-(3-pyridyl) Butanoic Acid, a Product of Nicotine 2'-Oxidation, and CYP2A6 Activity in Japanese Americans, Native Hawaiians, and Whites.

Smoking intensity varies across smokers and is influenced by individual variability in the metabolism of nicotine, the major addictive agent in tobacco. Therefore, lung cancer risk, which varies by racial ethnic group, is influenced by the primary catalyst of nicotine metabolism, cytochrome P450 2A6 (CYP2A6). In smokers, CYP2A6 catalyzes nicotine 5'-oxidation. In vitro, CYP2A6 also catalyzes, to a much lower extent, 2'-oxidation, which leads to the formation of 4-hydroxy-4-(3-pyridyl) butanoic acid (hydroxy acid). The urinary concentration of hydroxy acid has been quantified in only a few small studies of White smokers. To quantitatively assess the importance of nicotine 2'-oxidation in smokers, an LC-MS/MS-based method was developed for the analysis of nicotine and ten metabolites in urine. The concentrations of nicotine and these metabolites were measured in 303 smokers (99 Whites, 99 Native Hawaiians, and 105 Japanese Americans), and the relative metabolism of nicotine by four pathways was determined. Metabolism by these pathways was also compared across quartiles of CYP2A6 activity (measured as the plasma ratio of 3-hydroxycotinine to cotinine). As reported previously and consistent with their average CYP2A6 activity, nicotine 5'-oxidation was highest in Whites and lowest in Japanese Americans. Nicotine N-glucuronidation and N-oxidation increased with decreasing CYP2A6 activity. However, the relative urinary concentration of hydroxy acid (mean, 2.3%; 95% CI, 2.2-2.4%) did not vary by ethnic group or by CYP2A6 activity. In summary, CYP2A6 is not an important catalyst of nicotine 2'-oxidation in smokers, nor does nicotine 2'-oxidation compensate for decreased CYP2A6 activity.

Identification of novel epithelial ovarian cancer loci in women of African ancestry.

Women of African ancestry have lower incidence of epithelial ovarian cancer (EOC) yet worse survival compared to women of European ancestry. We conducted a genome-wide association study in African ancestry women with 755 EOC cases, including 537 high-grade serous ovarian carcinomas (HGSOC) and 1,235 controls. We identified four novel loci with suggestive evidence of association with EOC (p < 1 × 10-6 ), including rs4525119 (intronic to AKR1C3), rs7643459 (intronic to LOC101927394), rs4286604 (12 kb 3' of UGT2A2) and rs142091544 (5 kb 5' of WWC1). For HGSOC, we identified six loci with suggestive evidence of association including rs37792 (132 kb 5' of follistatin [FST]), rs57403204 (81 kb 3' of MAGEC1), rs79079890 (LOC105376360 intronic), rs66459581 (5 kb 5' of PRPSAP1), rs116046250 (GABRG3 intronic) and rs192876988 (32 kb 3' of GK2). Among the identified variants, two are near genes known to regulate hormones and diseases of the ovary (AKR1C3 and FST), and two are linked to cancer (AKR1C3 and MAGEC1). In follow-up studies of the 10 identified variants, the GK2 region SNP, rs192876988, showed an inverse association with EOC in European ancestry women (p = 0.002), increased risk of ER positive breast cancer in African ancestry women (p = 0.027) and decreased expression of GK2 in HGSOC tissue from African ancestry women (p = 0.004). A European ancestry-derived polygenic risk score showed positive associations with EOC and HGSOC in women of African ancestry suggesting shared genetic architecture. Our investigation presents evidence of variants for EOC shared among European and African ancestry women and identifies novel EOC risk loci in women of African ancestry.

Influence of UGT2B10 Genotype on Urinary Excretion of 4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanol- N-glucuronide by African American Smokers.

At similar smoking levels, African American's lung cancer risk is as much as twice that of whites. We hypothesized that racial/ethnic differences in UDP-glucuronosyltransferase (UGT)-catalyzed glucuronidation of 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanol (NNAL), a detoxication pathway for the tobacco-specific lung carcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) may contribute to this variable risk. UGT2B10 catalyzes NNAL- N-glucuronidation, and a UGT2B10 splice variant is common among African Americans. Smokers from two independent studies were genotyped for this variant (rs116294140) and an Asp67Tyr variant (rs61750900), and urinary NNAL and NNAL-glucuronide concentrations were quantified. In the first, no significant differences in NNAL- N-glucuronidation between African Americans ( n = 257) and whites ( n = 354) or between homozygous carriers of UGT2B10 variants (genetic score 2) and noncarriers (score 0) were detected. However, total NNAL glucuronidation by score 2 compared to score 0 smokers was lower (68.9 vs 71.2%, p < 0.0001). For NNAL- N-glucuronide to be more precisely quantified in a second study, a sensitive high-resolution LC-MS/MS-based method, which separated NNAL, NNAL- O-glucuronide, and NNAL- N-glucuronide prior to analysis, was developed. In this study, the excretion of total NNAL (free plus glucuronides) by African American ( n = 52) and white ( n = 54) smokers was not different; however, total NNAL glucuronidation by African Americans (64.0%) was slightly less than by whites (68.3%, p = 0.05). The mean NNAL- N-glucuronidation by African Americans was much lower than for whites (14 vs 24.9%, p < 0.00001), but the NNAL- O-glucuronidation was greater (50.0 vs 43.3%, p = 0.013). UGT2B10 genotype influenced NNAL- N-glucuronidation; the geometric mean percentage N-glucuronidation was 22.5% for smokers with genetic score 0 ( n = 57) and 11.2% for score 2 ( n = 11). In summary, the high prevalence of a UGT2B10 splice variant among African Americans results in lower NNAL- N-glucuronidation but only a small decrease in total NNAL glucuronidation. Therefore, despite the significant contribution of UGT2B10 to NNAL- N-glucuronidation, the UGT2B10 genotype does not play a large role in NNAL detoxication. Any decrease in N-glucuronidation was accompanied by a parallel increase in O-glucuronidation.

Comparison of the clinical prediction model PREMM(1,2,6) and molecular testing for the systematic identification of Lynch syndrome in colorectal cancer.

BACKGROUND: Lynch syndrome is caused by germline mismatch repair (MMR) gene mutations. The PREMM(1,2,6) model predicts the likelihood of a MMR gene mutation based on personal and family cancer history. OBJECTIVE: To compare strategies using PREMM(1,2,6) and tumour testing (microsatellite instability (MSI) and/or immunohistochemistry (IHC) staining) to identify mutation carriers. DESIGN: Data from population-based or clinic-based patients with colorectal cancers enrolled through the Colon Cancer Family Registry were analysed. Evaluation included MSI, IHC and germline mutation analysis for MLH1, MSH2, MSH6 and PMS2. Personal and family cancer histories were used to calculate PREMM(1,2,6) predictions. Discriminative ability to identify carriers from non-carriers using the area under the receiver operating characteristic curve (AUC) was assessed. Predictions were based on logistic regression models for (1) cancer assessment using PREMM(1,2,6), (2) MSI, (3) IHC for loss of any MMR protein expression, (4) MSI+IHC, (5) PREMM(1,2,6)+MSI, (6) PREMM(1,2,6)+IHC, (7) PREMM(1,2,6)+IHC+MSI. RESULTS: Among 1651 subjects, 239 (14%) had mutations (90 MLH1, 125 MSH2, 24 MSH6). PREMM(1,2,6) discriminated well with AUC 0.90 (95% CI 0.88 to 0.92). MSI alone, IHC alone, or MSI+IHC each had lower AUCs: 0.77, 0.82 and 0.82, respectively. The added value of IHC+PREMM(1,2,6) was slightly greater than PREMM(1,2,6)+MSI (AUC 0.94 vs 0.93). Adding MSI to PREMM(1,2,6)+IHC did not improve discrimination. CONCLUSION: PREMM(1,2,6) and IHC showed excellent performance in distinguishing mutation carriers from non-carriers and performed best when combined. MSI may have a greater role in distinguishing Lynch syndrome from other familial colorectal cancer subtypes among cases with high PREMM(1,2,6) scores where genetic evaluation does not disclose a MMR mutation.

Use of Deep Learning to Evaluate Tumor Microenvironmental Features for Prediction of Colon Cancer Recurrence.

Deep learning may detect biologically important signals embedded in tumor morphologic features that confer distinct prognoses. Tumor morphologic features were quantified to enhance patient risk stratification within DNA mismatch repair (MMR) groups using deep learning. Using a quantitative segmentation algorithm (QuantCRC) that identifies 15 distinct morphologic features, we analyzed 402 resected stage III colon carcinomas [191 deficient (d)-MMR; 189 proficient (p)-MMR] from participants in a phase III trial of FOLFOX-based adjuvant chemotherapy. Results were validated in an independent cohort (176 d-MMR; 1,094 p-MMR). Association of morphologic features with clinicopathologic variables, MMR, KRAS, BRAFV600E, and time-to-recurrence (TTR) was determined. Multivariable Cox proportional hazards models were developed to predict TTR. Tumor morphologic features differed significantly by MMR status. Cancers with p-MMR had more immature desmoplastic stroma. Tumors with d-MMR had increased inflammatory stroma, epithelial tumor-infiltrating lymphocytes (TIL), high-grade histology, mucin, and signet ring cells. Stromal subtype did not differ by BRAFV600E or KRAS status. In p-MMR tumors, multivariable analysis identified tumor-stroma ratio (TSR) as the strongest feature associated with TTR [HRadj 2.02; 95% confidence interval (CI), 1.14-3.57; P = 0.018; 3-year recurrence: 40.2% vs. 20.4%; Q1 vs. Q2-4]. Among d-MMR tumors, extent of inflammatory stroma (continuous HRadj 0.98; 95% CI, 0.96-0.99; P = 0.028; 3-year recurrence: 13.3% vs. 33.4%, Q4 vs. Q1) and N stage were the most robust prognostically. Association of TSR with TTR was independently validated. In conclusion, QuantCRC can quantify morphologic differences within MMR groups in routine tumor sections to determine their relative contributions to patient prognosis, and may elucidate relevant pathophysiologic mechanisms driving prognosis. SIGNIFICANCE: A deep learning algorithm can quantify tumor morphologic features that may reflect underlying mechanisms driving prognosis within MMR groups. TSR was the most robust morphologic feature associated with TTR in p-MMR colon cancers. Extent of inflammatory stroma and N stage were the strongest prognostic features in d-MMR tumors. TIL density was not independently prognostic in either MMR group.

A multifactorial likelihood model for MMR gene variant classification incorporating probabilities based on sequence bioinformatics and tumor characteristics: a report from the Colon Cancer Family Registry.

Mismatch repair (MMR) gene sequence variants of uncertain clinical significance are often identified in suspected Lynch syndrome families, and this constitutes a challenge for both researchers and clinicians. Multifactorial likelihood model approaches provide a quantitative measure of MMR variant pathogenicity, but first require input of likelihood ratios (LRs) for different MMR variation-associated characteristics from appropriate, well-characterized reference datasets. Microsatellite instability (MSI) and somatic BRAF tumor data for unselected colorectal cancer probands of known pathogenic variant status were used to derive LRs for tumor characteristics using the Colon Cancer Family Registry (CFR) resource. These tumor LRs were combined with variant segregation within families, and estimates of prior probability of pathogenicity based on sequence conservation and position, to analyze 44 unclassified variants identified initially in Australasian Colon CFR families. In addition, in vitro splicing analyses were conducted on the subset of variants based on bioinformatic splicing predictions. The LR in favor of pathogenicity was estimated to be ~12-fold for a colorectal tumor with a BRAF mutation-negative MSI-H phenotype. For 31 of the 44 variants, the posterior probabilities of pathogenicity were such that altered clinical management would be indicated. Our findings provide a working multifactorial likelihood model for classification that carefully considers mode of ascertainment for gene testing.

Genetic Susceptibility to Nonalcoholic Fatty Liver Disease and Risk for Pancreatic Cancer: Mendelian Randomization.

BACKGROUND: There are conflicting data on whether nonalcoholic fatty liver disease (NAFLD) is associated with susceptibility to pancreatic cancer. Using Mendelian randomization (MR), we investigated the relationship between genetic predisposition to NAFLD and risk for pancreatic cancer. METHODS: Data from genome-wide association studies (GWAS) within the Pancreatic Cancer Cohort Consortium (PanScan; cases n = 5,090, controls n = 8,733) and the Pancreatic Cancer Case Control Consortium (PanC4; cases n = 4,163, controls n = 3,792) were analyzed. We used data on 68 genetic variants with four different MR methods [inverse variance weighting (IVW), MR-Egger, simple median, and penalized weighted median] separately to predict genetic heritability of NAFLD. We then assessed the relationship between each of the four MR methods and pancreatic cancer risk, using logistic regression to calculate ORs and 95% confidence intervals (CI), adjusting for PC risk factors, including obesity and diabetes. RESULTS: No association was found between genetically predicted NAFLD and pancreatic cancer risk in the PanScan or PanC4 samples [e.g., PanScan, IVW OR, 1.04; 95% confidence interval (CI), 0.88-1.22; MR-Egger OR, 0.89; 95% CI, 0.65-1.21; PanC4, IVW OR, 1.07; 95% CI, 0.90-1.27; MR-Egger OR, 0.93; 95% CI, 0.67-1.28]. None of the four MR methods indicated an association between genetically predicted NAFLD and pancreatic cancer risk in either sample. CONCLUSIONS: Genetic predisposition to NAFLD is not associated with pancreatic cancer risk. IMPACT: Given the close relationship between NAFLD and metabolic conditions, it is plausible that any association between NAFLD and pancreatic cancer might reflect host metabolic perturbations (e.g., obesity, diabetes, or metabolic syndrome) and does not necessarily reflect a causal relationship between NAFLD and pancreatic cancer.

Asthma and lung cancer risk: a systematic investigation by the International Lung Cancer Consortium.

Asthma has been hypothesized to be associated with lung cancer (LC) risk. We conducted a pooled analysis of 16 studies in the International Lung Cancer Consortium (ILCCO) to quantitatively assess this association and compared the results with 36 previously published studies. In total, information from 585 444 individuals was used. Study-specific measures were combined using random effects models. A meta-regression and subgroup meta-analyses were performed to identify sources of heterogeneity. The overall LC relative risk (RR) associated with asthma was 1.28 [95% confidence intervals (CIs) = 1.16-1.41] but with large heterogeneity (I(2) = 73%, P < 0.001) between studies. Among ILCCO studies, an increased risk was found for squamous cell (RR = 1.69, 95%, CI = 1.26-2.26) and for small-cell carcinoma (RR = 1.71, 95% CI = 0.99-2.95) but was weaker for adenocarcinoma (RR = 1.09, 95% CI = 0.88-1.36). The increased LC risk was strongest in the 2 years after asthma diagnosis (RR = 2.13, 95% CI = 1.09-4.17) but subjects diagnosed with asthma over 10 years prior had no or little increased LC risk (RR = 1.10, 95% CI = 0.94-1.30). Because the increased incidence of LC was chiefly observed in small cell and squamous cell lung carcinomas, primarily within 2 years of asthma diagnosis and because the association was weak among never smokers, we conclude that the association may not reflect a causal effect of asthma on the risk of LC.

Parent of origin effects on age at colorectal cancer diagnosis.

Genomic imprinting refers to a parent-of-origin specific effect on gene expression. At least 1% of genes in the human genome are modulated in this manner. We sought evidence for genomic imprinting in colorectal cancer by studying the ages at diagnosis in the offspring of 2,061 parent-child pairs in which both parent and child were affected by nonsyndromic colorectal cancer. Families were ascertained through the colon Cancer Family Registry [http://epi.grants.cancer.gov/CFR/] from both population-based and clinic-based sources. We found that the affected offspring of affected fathers were on average younger than offspring of affected mothers (55.8 vs. 53.7 years; p = 0.0003), but when divided into sons and daughters, this difference was driven entirely by younger age at diagnosis in daughters of affected fathers compared to sons (52.3 years vs. 55.1 years; p = 0.0004). A younger age at diagnosis in affected daughters of affected fathers was also observable in various subsets including families that met Amsterdam II Criteria, families that did not meet Amsterdam Criteria, and in families with documented normal DNA mismatch repair in tumors. Imprinting effects are not expected to be affected by the sex of the offspring. Possible explanations for these unexpected findings include: (i) an imprinted gene on the pseudoautosomal regions of the X chromosome; (ii) an imprinted autosomal gene that affects a sex-specific pathway; or (iii) an X-linked gene unmasked because of colonic tissue-specific preferential inactivation of the maternal X chromosome.

Germline MutY human homologue mutations and colorectal cancer: a multisite case-control study.

BACKGROUND & AIMS: The MutY human homologue (MYH) gene is a member of the base-excision repair pathway involved in the repair of oxidative DNA damage. The objective of this study was to determine colorectal cancer (CRC) risk associated with mutations in the MYH gene. METHODS: A total of 3811 CRC cases and 2802 controls collected from a multisite CRC registry were screened for 9 germline MYH mutations; subjects with any mutation underwent screening of the entire MYH gene. Logistic regression was used to estimate age- and sex-adjusted odds ratios (AOR). Clinicopathologic and epidemiologic data were reviewed to describe the phenotype associated with MYH mutation status and assess for potential confounding and effect modification. RESULTS: Twenty-seven cases and 1 control subject carried homozygous or compound heterozygous MYH mutations (AOR, 18.1; 95% confidence interval, 2.5-132.7). CRC cases with homozygous/compound heterozygous mutations were younger at diagnosis (P=.01), had a higher proportion of right-sided (P=.01), synchronous cancers (P<.01), and personal history of adenomatous polyps (P=.003). Heterozygous MYH mutations were identified in 87 CRC cases and 43 controls; carriers were at increased risk of CRC (AOR, 1.48; 95% confidence interval, 1.02-2.16). There was a higher prevalence of low-frequency microsatellite instability (MSI) in tumors from heterozygous and homozygous/compound heterozygous MYH mutation carriers (P=.02); MSI status modified the CRC risk associated with heterozygous MYH mutations (P interaction<.001). CONCLUSIONS: Homozygous/compound heterozygous MYH mutations account for less than 1% of CRC cases. Heterozygous carriers are at increased risk of CRC. Further studies are needed to understand the possible interaction between the base excision repair and low-frequency MSI pathways.

Alpha-1-antitrypsin deficiency and smoking as risk factors for mismatch repair deficient colorectal cancer: a study from the colon cancer family registry.

In a previous study, alpha-1-antitrypsin (A1AT) deficiency alleles were found to be over represented among individuals with microsatellite unstable (MSI-high) colorectal cancers, and this was most significant in former or current smokers. We evaluated this association in a larger case-control study, stratified by microsatellite instability phenotypes. Concordant with prior observations, gender (female) and smoking history were positively associated with colorectal cancers having an MSI-high phenotype. No difference in frequency of A1AT deficiency alleles was found between cases and controls, irrespective of the MSI subtype.

Contribution of Germline Predisposition Gene Mutations to Breast Cancer Risk in African American Women.

BACKGROUND: The risks of breast cancer in African American (AA) women associated with inherited mutations in breast cancer predisposition genes are not well defined. Thus, whether multigene germline hereditary cancer testing panels are applicable to this population is unknown. We assessed associations between mutations in panel-based genes and breast cancer risk in 5054 AA women with breast cancer and 4993 unaffected AA women drawn from 10 epidemiologic studies. METHODS: Germline DNA samples were sequenced for mutations in 23 cancer predisposition genes using a QIAseq multiplex amplicon panel. Prevalence of mutations and odds ratios (ORs) for associations with breast cancer risk were estimated with adjustment for study design, age, and family history of breast cancer. RESULTS: Pathogenic mutations were identified in 10.3% of women with estrogen receptor (ER)-negative breast cancer, 5.2% of women with ER-positive breast cancer, and 2.3% of unaffected women. Mutations in BRCA1, BRCA2, and PALB2 were associated with high risks of breast cancer (OR = 47.55, 95% confidence interval [CI] = 10.43 to >100; OR = 7.25, 95% CI = 4.07 to 14.12; OR = 8.54, 95% CI = 3.67 to 24.95, respectively). RAD51D mutations were associated with high risk of ER-negative disease (OR = 7.82, 95% CI = 1.61 to 57.42). Moderate risks were observed for CHEK2, ATM, ERCC3, and FANCC mutations with ER-positive cancer, and RECQL mutations with all breast cancer. CONCLUSIONS: The study identifies genes that predispose to breast cancer in the AA population, demonstrates the validity of current breast cancer testing panels for use in AA women, and provides a basis for increased referral of AA patients for cancer genetic testing.

A Transcriptome-Wide Association Study Identifies Novel Candidate Susceptibility Genes for Pancreatic Cancer.

BACKGROUND: Although 20 pancreatic cancer susceptibility loci have been identified through genome-wide association studies in individuals of European ancestry, much of its heritability remains unexplained and the genes responsible largely unknown. METHODS: To discover novel pancreatic cancer risk loci and possible causal genes, we performed a pancreatic cancer transcriptome-wide association study in Europeans using three approaches: FUSION, MetaXcan, and Summary-MulTiXcan. We integrated genome-wide association studies summary statistics from 9040 pancreatic cancer cases and 12 496 controls, with gene expression prediction models built using transcriptome data from histologically normal pancreatic tissue samples (NCI Laboratory of Translational Genomics [n = 95] and Genotype-Tissue Expression v7 [n = 174] datasets) and data from 48 different tissues (Genotype-Tissue Expression v7, n = 74-421 samples). RESULTS: We identified 25 genes whose genetically predicted expression was statistically significantly associated with pancreatic cancer risk (false discovery rate < .05), including 14 candidate genes at 11 novel loci (1p36.12: CELA3B; 9q31.1: SMC2, SMC2-AS1; 10q23.31: RP11-80H5.9; 12q13.13: SMUG1; 14q32.33: BTBD6; 15q23: HEXA; 15q26.1: RCCD1; 17q12: PNMT, CDK12, PGAP3; 17q22: SUPT4H1; 18q11.22: RP11-888D10.3; and 19p13.11: PGPEP1) and 11 at six known risk loci (5p15.33: TERT, CLPTM1L, ZDHHC11B; 7p14.1: INHBA; 9q34.2: ABO; 13q12.2: PDX1; 13q22.1: KLF5; and 16q23.1: WDR59, CFDP1, BCAR1, TMEM170A). The association for 12 of these genes (CELA3B, SMC2, and PNMT at novel risk loci and TERT, CLPTM1L, INHBA, ABO, PDX1, KLF5, WDR59, CFDP1, and BCAR1 at known loci) remained statistically significant after Bonferroni correction. CONCLUSIONS: By integrating gene expression and genotype data, we identified novel pancreatic cancer risk loci and candidate functional genes that warrant further investigation.

Underlying liver disease and advanced stage liver cancer are associated with elevated neutrophil-lymphocyte ratio.

BACKGROUND/AIMS: Inflammation-based scores, such as the neutrophil-to-lymphocyte ratio (NLR), have been associated with prognosis in hepatocellular carcinoma (HCC); but variable cut-off values and potential lack of specificity have limited the utility of NLR. This study evaluates NLR in a large cohort of HCC patients. METHODS: We retrospectively reviewed 789 HCC cases (1993-2017) for demographics, tumor characteristics, treatment, and survival. NLR was stratified into NLR ≥1.5 and NLR ≥3 and analyzed for correlation with American Joint Committee on Cancer (AJCC) and Barcelona Clinic Liver Cancer (BCLC) stages. In 235 patients who underwent liver resection, survival and recurrence were evaluated by NLR. RESULTS: In 789 HCC cases, mean NLR was increased with advanced AJCC and BCLC stages. Hepatitis C patients were less likely to have NLR ≥1.5 and ≥3. Non-alcoholic steatohepatitis patients were more likely to have NLR ≥3. Patients with tumor size >5 cm, rupture, or macrovascular invasion were more likely to have NLR ≥3. In patients treated with resection, NLR ≥3 predicted early recurrence (odds ratio [OR] 4.14, P<0.01) and overall recurrence (OR 4.05, P<0.01). Mean NLR was 4.30 in those with recurrence and 2.75 in those without recurrence. Patients with NLR ≥3 showed significantly worse survival compared to those with NLR <3 (P<0.01 by log-rank test). CONCLUSION: Elevated NLR is associated with advanced cancer stage and aggressive tumor characteristics, such as large size, rupture, and invasion. NLR ≥3 was associated with early and overall recurrence after resection but varied with etiology. NLR may be a useful biomarker in predicting recurrence for HCC patients undergoing curative resection, but further studies are required to elucidate the effect of disease etiology.

Agnostic Pathway/Gene Set Analysis of Genome-Wide Association Data Identifies Associations for Pancreatic Cancer.

BACKGROUND: Genome-wide association studies (GWAS) identify associations of individual single-nucleotide polymorphisms (SNPs) with cancer risk but usually only explain a fraction of the inherited variability. Pathway analysis of genetic variants is a powerful tool to identify networks of susceptibility genes. METHODS: We conducted a large agnostic pathway-based meta-analysis of GWAS data using the summary-based adaptive rank truncated product method to identify gene sets and pathways associated with pancreatic ductal adenocarcinoma (PDAC) in 9040 cases and 12 496 controls. We performed expression quantitative trait loci (eQTL) analysis and functional annotation of the top SNPs in genes contributing to the top associated pathways and gene sets. All statistical tests were two-sided. RESULTS: We identified 14 pathways and gene sets associated with PDAC at a false discovery rate of less than 0.05. After Bonferroni correction (P ≤ 1.3 × 10-5), the strongest associations were detected in five pathways and gene sets, including maturity-onset diabetes of the young, regulation of beta-cell development, role of epidermal growth factor (EGF) receptor transactivation by G protein-coupled receptors in cardiac hypertrophy pathways, and the Nikolsky breast cancer chr17q11-q21 amplicon and Pujana ATM Pearson correlation coefficient (PCC) network gene sets. We identified and validated rs876493 and three correlating SNPs (PGAP3) and rs3124737 (CASP7) from the Pujana ATM PCC gene set as eQTLs in two normal derived pancreas tissue datasets. CONCLUSION: Our agnostic pathway and gene set analysis integrated with functional annotation and eQTL analysis provides insight into genes and pathways that may be biologically relevant for risk of PDAC, including those not previously identified.

Evaluation of vitamin D biosynthesis and pathway target genes reveals UGT2A1/2 and EGFR polymorphisms associated with epithelial ovarian cancer in African American Women.

An association between genetic variants in the vitamin D receptor (VDR) gene and epithelial ovarian cancer (EOC) was previously reported in women of African ancestry (AA). We sought to examine associations between genetic variants in VDR and additional genes from vitamin D biosynthesis and pathway targets (EGFR, UGT1A, UGT2A1/2, UGT2B, CYP3A4/5, CYP2R1, CYP27B1, CYP24A1, CYP11A1, and GC). Genotyping was performed using the custom-designed 533,631 SNP Illumina OncoArray with imputation to the 1,000 Genomes Phase 3 v5 reference set in 755 EOC cases, including 537 high-grade serous (HGSOC), and 1,235 controls. All subjects are of African ancestry (AA). Logistic regression was performed to estimate odds ratios (OR) and 95% confidence intervals (CI). We further evaluated statistical significance of selected SNPs using the Bayesian False Discovery Probability (BFDP). A significant association with EOC was identified in the UGT2A1/2 region for the SNP rs10017134 (per allele OR = 1.4, 95% CI = 1.2-1.7, P = 1.2 × 10-6 , BFDP = 0.02); and an association with HGSOC was identified in the EGFR region for the SNP rs114972508 (per allele OR = 2.3, 95% CI = 1.6-3.4, P = 1.6 × 10-5 , BFDP = 0.29) and in the UGT2A1/2 region again for rs1017134 (per allele OR = 1.4, 95% CI = 1.2-1.7, P = 2.3 × 10-5 , BFDP = 0.23). Genetic variants in the EGFR and UGT2A1/2 may increase susceptibility of EOC in AA women. Future studies to validate these findings are warranted. Alterations in EGFR and UGT2A1/2 could perturb enzyme efficacy, proliferation in ovaries, impact and mark susceptibility to EOC.

Haplotypes of the estrogen receptor beta gene and breast cancer risk.

Exposure to exogenous (oral contraceptives, postmenopausal hormone therapy) and endogenous (number of ovulatory cycles, adiposity) steroid hormones is associated with breast cancer risk. Breast cancer risk associated with these exposures could hypothetically be modified by genes in the steroid hormone synthesis, metabolism and signaling pathways. Estrogen receptors are the first step along the path of signaling cell growth and development upon stimulation with estrogens. The National Cancer Institute Breast and Prostate Cancer Cohort Consortium has systematically selected haplotype tagging SNPs in genes along the steroid hormone synthesis, metabolism and binding pathways, including the estrogen receptor beta (ESR2) gene. Four htSNPs tag the 6 major (>5% frequency) haplotypes of the ESR2 gene. These polymorphisms have been genotyped in 5,789 breast cancer cases and 7,761 controls nested within the American Cancer Society Cancer Prevention Study II, European Prospective Investigation into Cancer and Nutrition, Multiethnic Cohort, Nurses' Health Study and Women's Health Study cohorts. None of the SNPs were independently associated with breast cancer risk. One haplotype of the ESR2 gene was associated with breast cancer risk before correction for multiple testing (OR 1.17, 95% CI 1.07-1.28, p = 0.0007). This haplotype remained associated with breast cancer risk after adjustment for multiple testing using a permutation procedure. There was no statistically significant heterogeneity in SNP or haplotype odds ratios across cohorts. These data suggest that inherited variants in ESR2 (while possibly conferring a small increased risk of breast cancer) are not associated with appreciable (OR > 1.2) changes in breast cancer risk among Caucasian women.

International Lung Cancer Consortium: pooled analysis of sequence variants in DNA repair and cell cycle pathways.

BACKGROUND: The International Lung Cancer Consortium was established in 2004. To clarify the role of DNA repair genes in lung cancer susceptibility, we conducted a pooled analysis of genetic variants in DNA repair pathways, whose associations have been investigated by at least 3 individual studies. METHODS: Data from 14 studies were pooled for 18 sequence variants in 12 DNA repair genes, including APEX1, OGG1, XRCC1, XRCC2, XRCC3, ERCC1, XPD, XPF, XPG, XPA, MGMT, and TP53. The total number of subjects included in the analysis for each variant ranged from 2,073 to 13,955 subjects. RESULTS: Four of the variants were found to be weakly associated with lung cancer risk with borderline significance: these were XRCC3 T241M [heterozygote odds ratio (OR), 0.89; 95% confidence interval (95% CI), 0.79-0.99 and homozygote OR, 0.84; 95% CI, 0.71-1.00] based on 3,467 cases and 5,021 controls from 8 studies, XPD K751Q (heterozygote OR, 0.99; 95% CI, 0.89-1.10 and homozygote OR, 1.19; 95% CI, 1.02-1.39) based on 6,463 cases and 6,603 controls from 9 studies, and TP53 R72P (heterozygote OR, 1.14; 95% CI, 1.00-1.29 and homozygote OR, 1.20; 95% CI, 1.02-1.42) based on 3,610 cases and 5,293 controls from 6 studies. OGG1 S326C homozygote was suggested to be associated with lung cancer risk in Caucasians (homozygote OR, 1.34; 95% CI, 1.01-1.79) based on 2,569 cases and 4,178 controls from 4 studies but not in Asians. The other 14 variants did not exhibit main effects on lung cancer risk. DISCUSSION: In addition to data pooling, future priorities of International Lung Cancer Consortium include coordinated genotyping and multistage validation for ongoing genome-wide association studies.

Haplotype analysis of the HSD17B1 gene and risk of breast cancer: a comprehensive approach to multicenter analyses of prospective cohort studies.

The 17beta-hydroxysteroid dehydrogenase 1 gene (HSD17B1) encodes 17HSD1, which catalyzes the final step of estradiol biosynthesis. Despite the important role of HSD17B1 in hormone metabolism, few epidemiologic studies of HSD17B1 and breast cancer have been conducted. This study includes 5,370 breast cancer cases and 7,480 matched controls from five large cohorts in the Breast and Prostate Cancer Cohort Consortium. We characterized variation in HSD17B1 by resequencing and dense genotyping a multiethnic sample and identified haplotype-tagging single nucleotide polymorphisms (htSNP) that capture common variation within a 33.3-kb region around HSD17B1. Four htSNPs, including the previously studied SNP rs605059 (S312G), were genotyped to tag five common haplotypes in all cases and controls. Conditional logistic regression was used to estimate odds ratios (OR) for disease. We found no evidence of association between common HSD17B1 haplotypes or htSNPs and overall risk of breast cancer. The OR for each haplotype relative to the most common haplotype ranged from 0.98 to 1.07 (omnibus test for association: X2 = 3.77, P = 0.58, 5 degrees of freedom). When cases were subdivided by estrogen receptor (ER) status, two common haplotypes were associated with ER-negative tumors (test for trend, Ps = 0.0009 and 0.0076; n = 353 cases). HSD17B1 variants that are common in Caucasians are not associated with overall risk of breast cancer; however, there was an association among the subset of ER-negative tumors. Although the probability that these ER-negative findings are false-positive results is high, these findings were consistent across each cohort examined and warrant further study.

Sequence variants of estrogen receptor beta and risk of prostate cancer in the National Cancer Institute Breast and Prostate Cancer Cohort Consortium.

BACKGROUND: Estrogen receptor beta (ESR2) may play a role in modulating prostate carcinogenesis through the regulation of genes related to cell proliferation and apoptosis. METHODS: We conducted nested case-control studies in the Breast and Prostate Cancer Cohort Consortium (BPC3) that pooled 8,323 prostate cancer cases and 9,412 controls from seven cohorts. Whites were the predominant ethnic group. We characterized genetic variation in ESR2 by resequencing exons in 190 breast and prostate cancer cases and genotyping a dense set of single nucleotide polymorphisms (SNP) spanning the locus in a multiethnic panel of 349 cancer-free subjects. We selected four haplotype-tagging SNPs (htSNP) to capture common ESR2 variation in Whites; these htSNPs were then genotyped in all cohorts. Conditional logistic regression models were used to assess the association between sequence variants of ESR2 and the risk of prostate cancer. We also investigated the effect modification by age, body mass index, and family history, as well as the association between sequence variants of ESR2 and advanced-stage (>or=T3b, N1, or M1) and high-grade (Gleason sum >or=8) prostate cancer, respectively. RESULTS: The four tag SNPs in ESR2 were not significantly associated with prostate cancer risk, individually. The global test for the influence of any haplotype on the risk of prostate cancer was not significant (P = 0.31). However, we observed that men carrying two copies of one of the variant haplotypes (TACC) had a 1.46-fold increased risk of prostate cancer (99% confidence interval, 1.06-2.01) compared with men carrying zero copies of this variant haplotype. No SNPs or haplotypes were associated with advanced stage or high grade of prostate cancer. CONCLUSION: In our analysis focused on genetic variation common in Whites, we observed little evidence for any substantial association of inherited variation in ESR2 with risk of prostate cancer. A nominally significant (P < 0.01) association between the TACC haplotype and prostate cancer risk under the recessive model could be a chance finding and, in any event, would seem to contribute only slightly to the overall burden of prostate cancer.