Search

Although next-generation sequencing has demonstrated great potential for novel gene discovery, confirming disease-causing genes after initial discovery remains challenging. Here, we applied a network analysis approach to prioritize candidate genes identified from whole-exome sequencing analysis of 98 cutaneous melanoma patients from 27 families. Using a network propagation method, we ranked candidate genes by their similarity to known disease genes in protein–protein interaction networks and identified gene clusters with functional connectivity. Using this approach, we identified several new candidate susceptibility genes that warrant future investigations such as NGLY1, IL1RN, FABP2, PRKDC, and PROSER2. The propagated network analysis also allowed us to link families that did not have common underlying genes but that carried variants in genes that interact on protein–protein interaction networks. In conclusion, our study provided an analysis perspective for gene prioritization in the context of genetic heterogeneity across families and prioritized top potential candidate susceptibility genes in our dataset.

Mounting evidence suggests that copy number variations (CNVs) can contribute to cancer susceptibility. The main goal of this study was to evaluate the role of germline CNVs in melanoma predisposition in high-risk melanoma families. We used genome-wide tiling comparative genomic hybridization and single nucleotide polymorphism arrays to characterize CNVs in 335 individuals (240 melanoma cases) from American melanoma-prone families (22 with germline CDKN2A or CDK4 mutations). We found that the global burden of overall CNVs (or deletions or duplications separately) was not significantly associated with case-control or CDKN2A/CDK4 mutation status after accounting for the familial dependence. However, we identified several rare CNVs that either involved known melanoma genes (e.g., PARP1, CDKN2A) or cosegregated with melanoma (duplication on 10q23.23, 3p12.2 and deletions on 8q424.3, 2q22.1) in families without mutations in known melanoma high-risk genes. Some of these CNVs were correlated with expression changes in disrupted genes based on RNASeq data from a subset of melanoma cases included in the CNV study. These results suggest that rare cosegregating CNVs may influence melanoma susceptibility in some melanoma-prone families and genes found in our study warrant further evaluation in future genetic analyses of melanoma.

For unknown reasons, there is huge variability in risk conferred by different HPV types and, remarkably, strong differences even between closely related variant lineages within each type. HPV16 is a uniquely powerful carcinogenic type, causing approximately half of cervical cancer and most other HPV-related cancers. To permit the large-scale study of HPV genome variability and precancer/cancer, starting with HPV16 and cervical cancer, we developed a high-throughput next-generation sequencing (NGS) whole-genome method. We designed a custom HPV16 AmpliSeq™ panel that generated 47 overlapping amplicons covering 99% of the genome sequenced on the Ion Torrent Proton platform. After validating with Sanger, the current “gold standard” of sequencing, in 89 specimens with concordance of 99.9%, we used our NGS method and custom annotation pipeline to sequence 796 HPV16-positive exfoliated cervical cell specimens. The median completion rate per sample was 98.0%.Our method enabled us to discover novel SNPs, large contiguous deletions suggestive of viral integration (OR of 27.3, 95% CI 3.3–222, P=0.002), and the sensitive detection of variant lineage coinfections. This method represents an innovative high-throughput, ultra-deep coverage technique for HPV genomic sequencing, which, in turn, enables the investigation of the role of genetic variation in HPV epidemiology and carcinogenesis. Keywords: HPV16, HPV epidemiology, HPV genomics

Waldenström macroglobulinemia (WM)/lymphoplasmacytic lymphoma (LPL) is a rare, chronic B-cell lymphoma with high heritability. We conduct a two-stage genome-wide association study of WM/LPL in 530 unrelated cases and 4362 controls of European ancestry and identify two high-risk loci associated with WM/LPL at 6p25.3 (rs116446171, near EXOC2 and IRF4; OR = 21.14, 95% CI: 14.40–31.03, P = 1.36 × 10−54) and 14q32.13 (rs117410836, near TCL1; OR = 4.90, 95% CI: 3.45–6.96, P = 8.75 × 10−19). Both risk alleles are observed at a low frequency among controls (~2–3%) and occur in excess in affected cases within families. In silico data suggest that rs116446171 may have functional importance, and in functional studies, we demonstrate increased reporter transcription and proliferation in cells transduced with the 6p25.3 risk allele. Although further studies are needed to fully elucidate underlying biological mechanisms, together these loci explain 4% of the familial risk and provide insights into genetic susceptibility to this malignancy., Waldenström macroglobulinemia (WM)/lymphoplasmacytic lymphoma (LPL) is a non-Hodgkin-type B cell lymphoma. Here, the authors identify two risk loci for WM/LPL in a two-stage GWAS involving a family-oversampling approach and provide evidence for a functional role of the non-coding SNP rs116446171.

We previously carried out a multi-stage genome-wide association study (GWAS) on lung cancer among never smokers in the Female Lung Cancer Consortium in Asia (FLCCA) (6,609 cases, 7,457 controls) that identified novel susceptibility loci at 10q25.2, 6q22.2, and 6p21.32, and confirmed two previously identified loci at 5p15.33 and 3q28. Household air pollution (HAP) attributed to solid fuel burning for heating and cooking, is the leading cause of the overall disease burden in Southeast Asia, and is known to contain lung carcinogens. To evaluate the gene–HAP interactions associated with lung cancer in loci independent of smoking, we analyzed data from studies participating in FLCCA with fuel use information available (n = 3; 1,731 cases; 1,349 controls). Coal use was associated with a 30 % increased risk of lung cancer (OR 1.3, 95 % CI 1.0–1.6). Among the five a priori SNPs identified by our GWAS, two showed a significant interaction with coal use (HLA Class II rs2395185, p = 0.02; TP63 rs4488809 (rs4600802), p = 0.04). The risk of lung cancer associated with coal exposure varied with the respective alleles for these two SNPs. Our observations provide evidence that genetic variation in HLA Class II and TP63 may modify the association between HAP and lung cancer risk. The roles played in the cell cycle and inflammation pathways by the proteins encoded by these two genes provide biological plausibility for these interactions; however, additional replication studies are needed in other non-smoking populations.

Known high-risk cutaneous malignant melanoma (CMM) genes account for melanoma risk in

A genome-wide association study (GWAS) of bladder cancer identified a genetic marker rs8102137 within the 19q12 region as a novel susceptibility variant. This marker is located upstream of the CCNE1 gene, which encodes cyclin E, a cell-cycle protein. We performed genetic fine-mapping analysis of the CCNE1 region using data from two bladder cancer GWAS (5,942 cases and 10,857 controls). We found that the original GWAS marker rs8102137 represents a group of 47 linked SNPs (with r2 ≥ 0.7) associated with increased bladder cancer risk. From this group, we selected a functional promoter variant rs7257330, which showed strong allele-specific binding of nuclear proteins in several cell lines. In both GWASs, rs7257330 was associated only with aggressive bladder cancer, with a combined per-allele OR = 1.18 [95% confidence interval (CI), 1.09–1.27, P = 4.67 × 10−5] versus OR = 1.01 (95% CI, 0.93–1.10, P = 0.79) for nonaggressive disease, with P = 0.0015 for case-only analysis. Cyclin E protein expression analyzed in 265 bladder tumors was increased in aggressive tumors (P = 0.013) and, independently, with each rs7257330-A risk allele (Ptrend = 0.024). Overexpression of recombinant cyclin E in cell lines caused significant acceleration of cell cycle. In conclusion, we defined the 19q12 signal as the first GWAS signal specific for aggressive bladder cancer. Molecular mechanisms of this genetic association may be related to cyclin E overexpression and alteration of cell cycle in carriers of CCNE1 risk variants. In combination with established bladder cancer risk factors and other somatic and germline genetic markers, the CCNE1 variants could be useful for inclusion into bladder cancer risk prediction models. Cancer Res; 74(20); 5808–18. ©2014 AACR.

We conducted a joint (pooled) analysis of three genome-wide association studies (GWAS) of esophageal squamous cell carcinoma (ESCC) in individuals of Chinese ancestry (5,337 ESCC cases and 5,787 controls) with 9,654 ESCC cases and 10,058 controls for follow-up. In a logistic regression model adjusted for age, sex, study and two eigenvectors, two new loci achieved genome-wide significance, marked by rs7447927 at 5q31.2 (per-allele odds ratio (OR) = 0.85, 95% confidence interval (CI) = 0.82-0.88; P = 7.72 × 10(-20)) and rs1642764 at 17p13.1 (per-allele OR = 0.88, 95% CI = 0.85-0.91; P = 3.10 × 10(-13)). rs7447927 is a synonymous SNP in TMEM173, and rs1642764 is an intronic SNP in ATP1B2, near TP53. Furthermore, a locus in the HLA class II region at 6p21.32 (rs35597309) achieved genome-wide significance in the two populations at highest risk for ESSC (OR = 1.33, 95% CI = 1.22-1.46; P = 1.99 × 10(-10)). Our joint analysis identifies new ESCC susceptibility loci overall as well as a new locus unique to the population in the Taihang Mountain region at high risk of ESCC.

Genome-wide association studies (GWAS) have mapped risk alleles for at least 10 distinct cancers to a small region of 63 000 bp on chromosome 5p15.33. This region harbors the TERT and CLPTM1L genes; the former encodes the catalytic subunit of telomerase reverse transcriptase and the latter may play a role in apoptosis. To investigate further the genetic architecture of common susceptibility alleles in this region, we conducted an agnostic subset-based meta-analysis (association analysis based on subsets) across six distinct cancers in 34 248 cases and 45 036 controls. Based on sequential conditional analysis, we identified as many as six independent risk loci marked by common single-nucleotide polymorphisms: five in the TERT gene (Region 1: rs7726159, P = 2.10 x 10(-39); Region 3: rs2853677, P = 3.30 x 10(-36) and PConditional = 2.36 x 10(-8); Region 4: rs2736098, P = 3.87 x 10(-12) and PConditional = 5.19 x 10(-6), Region 5: rs13172201, P = 0.041 and PConditional = 2.04 x 10(-6); and Region 6: rs10069690, P = 7.49 x 10(-15) and PConditional = 5.35 x 10(-7)) and one in the neighboring CLPTM1L gene (Region 2: rs451360; P = 1.90 x 10(-18) and PConditional = 7.06 x 10(-16)). Between three and five cancers mapped to each independent locus with both risk-enhancing and protective effects. Allele-specific effects on DNA methylation were seen for a subset of risk loci, indicating that methylation and subsequent effects on gene expression may contribute to the biology of risk variants on 5p15.33. Our results provide strong support for extensive pleiotropy across this region of 5p15.33, to an extent not previously observed in other cancer susceptibility loci.

Dysplastic nevi (DN) is a strong risk factor for cutaneous malignant melanoma (CMM), and it frequently occurs in melanoma-prone families. To identify genetic variants for DN, we genotyped 677 tagSNPs in 38 melanoma candidate genes that are involved in pigmentation, DNA repair, cell cycle control, and melanocyte proliferation pathways in a total of 504 individuals (310 with DN, 194 without DN) from 53 melanoma-prone families (23 CDKN2A mutation positive and 30 negative). Conditional logistic regression, conditioning on families, was used to estimate the association between DN and each single-nucleotide polymorphism (SNP) separately, adjusted for age, sex, CMM, and CDKN2A status. P-values for SNPs in the same gene were combined to yield gene-specific P-values. Two genes, CDK6 (cyclin-dependent kinase 6) and XRCC1, were significantly associated with DN after Bonferroni correction for multiple testing (P=0.0001 and 0.00025, respectively), whereas neither gene was significantly associated with CMM. Associations for CDK6 SNPs were stronger in CDKN2A mutation–positive families (rs2079147, Pinteraction=0.0033), whereas XRCC1 SNPs had similar effects in mutation-positive and -negative families. The association for one of the associated SNPs in XRCC1 (rs25487) was replicated in two independent data sets (random-effect meta-analysis: P

Genome-wide association studies (GWAS) have identified 76 variants associated with prostate cancer risk predominantly in populations of European ancestry. To identify additional susceptibility loci for this common cancer, we conducted a meta-analysis of > 10 million SNPs in 43,303 prostate cancer cases and 43,737 controls from studies in populations of European, African, Japanese and Latino ancestry. Twenty-three new susceptibility loci were identified at association P < 5 × 10(-8); 15 variants were identified among men of European ancestry, 7 were identified in multi-ancestry analyses and 1 was associated with early-onset prostate cancer. These 23 variants, in combination with known prostate cancer risk variants, explain 33% of the familial risk for this disease in European-ancestry populations. These findings provide new regions for investigation into the pathogenesis of prostate cancer and demonstrate the usefulness of combining ancestrally diverse populations to discover risk loci for disease.

Chronic inflammation has been implicated in the etiology of colorectal adenoma and cancer; however, few key inflammatory genes mediating this relationship have been identified. In this study, we investigated the association of germline variation in innate immunity genes in relation to the risk of colorectal neoplasia. Our study was based on the analysis of samples collected from the prostate, lung, colorectal and ovarian (PLCO) Cancer Screening Trial. We investigated the association between 196 tag single nucleotide polymorphisms (SNPs) in 20 key innate immunity genes with risk of advanced colorectal adenoma and cancer in 719 adenoma cases, 481 cancer cases and 719 controls. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CIs). After Bonferroni correction, the AG/GG genotype of rs5995355, which is upstream of NCF4, was associated with an increased risk of colorectal cancer (OR = 2.43, 95% CI = 1.73-3.39; p < 0.0001). NCF4 is part of the NAPDH complex, a key factor in biochemical pathways and the innate immune response. While not definitive, our analyses suggest that the variant allele does not affect expression of NCF4, but rather modulates activity of the NADPH complex. Additional studies on the functional consequences of rs5995355 in NCF4 may help to clarify the mechanistic link between inflammation and colorectal cancer.

We performed a multistage genome-wide association study including 7,683 individuals with pancreatic cancer and 14,397 controls of European descent. Four new loci reached genome-wide significance: rs6971499 at 7q32.3 (LINC-PINT, per-allele odds ratio (OR) = 0.79, 95% confidence interval (CI) 0.74-0.84, P = 3.0 x 10(-12)), rs7190458 at 16q23.1 (BCAR1/CTRB1/CTRB2, OR = 1.46, 95% CI 1.30-1.65, P = 1.1 x 10(-10)), rs9581943 at 13q12.2 (PDX1, OR = 1.15, 95% CI 1.10-1.20, P = 2.4 x 10(-9)) and rs16986825 at 22q12.1 (ZNRF3, OR = 1.18, 95% CI 1.12-1.25, P = 1.2 x 10(-8)). We identified an independent signal in exon 2 of TERT at the established region 5p15.33 (rs2736098, OR = 0.80, 95% CI 0.76-0.85, P = 9.8 x 10(-14)). We also identified a locus at 8q24.21 (rs1561927, P = 1.3 x 10(-7)) that approached genome-wide significance located 455 kb telomeric of PVT1. Our study identified multiple new susceptibility alleles for pancreatic cancer that are worthy of follow-up studies.

Cutaneous malignant melanoma (CMM) is an etiologically heterogeneous disease with genetic, environmental (sun exposure) and host (pigmentation/nevi) factors and their interactions contributing to risk. Genetic variants in DNA repair genes may be particularly important since their altered function in response to sun exposure-related DNA damage maybe related to risk for CMM. However, systematic evaluations of genetic variants in DNA repair genes are limited, particularly in high-risk families. We comprehensively analyzed DNA repair gene polymorphisms and CMM risk in melanoma-prone families with/without CDKN2A mutations. A total of 586 individuals (183 CMM) from 53 families (23 CDKN2A (+), 30 CDKN2A (-)) were genotyped for 2964 tagSNPs in 131 DNA repair genes. Conditional logistic regression, conditioning on families, was used to estimate trend p-values, odds ratios and 95% confidence intervals for the association between CMM and each SNP separately, adjusted for age and sex. p-Values for SNPs in the same gene were combined to yield gene specific p-values. Two genes, POLN and PRKDC, were significantly associated with melanoma after Bonferroni correction for multiple testing (p = 0.0003 and 0.00035, respectively). DCLRE1B showed suggestive association (p = 0.0006). 28 ∼ 56% of genotyped SNPs in these genes had single SNP p < 0.05. The most significant SNPs in POLN and PRKDC had similar effects in CDKN2A (+) and CDKN2A (-) families. Our finding suggests that polymorphisms in DNA repair genes, POLN and PRKDC, were associated with increased melanoma risk in melanoma families with and without CDKN2A mutations.

Few studies have examined the relationship between human leukocyte antigen (HLA) polymorphisms and adult glioma, particularly at class II loci. We evaluated the association between selected HLA class II polymorphisms and adult glioma in a large, hospital-based case-control study, using unconditional logistic regression. DQB1*06 (OR = 1.67, 95% CI = 1.17–2.39) and DRB1*13 (OR = 1.69, 95% CI = 1.08–2.64) alleles were associated with an increased risk of glioma, while the DQB1*05 allele showed an inverse association (OR = 0.63, 95% CI = 0.43–0.93). These results, which were of borderline significance once controlled for the false discovery rate, suggest a potential role for the DQB1*06 , DQB1*05 , and DRB1*13 alleles in glioma susceptibility.

As genome-wide association studies of breast cancer are replicating findings and refinement studies are narrowing the signal location, additional efforts are necessary to elucidate the underlying functional relationships. One approach is to evaluate variation in risk by genotype based on known breast carcinogens, such as ionizing radiation. Given the public health concerns associated with recent increases in medical radiation exposure, this approach may also identify potentially susceptible subpopulations. We examined interaction between 27 newly identified breast cancer risk alleles (identified within the NCI Cancer Genetic Markers of Susceptibility and the Breast Cancer Association Consortium genome-wide association studies) and occupational and medical diagnostic radiation exposure among 859 cases and 1083 controls nested within the United States Radiologic Technologists cohort. We did not find significant variation in the radiation-related breast cancer risk for the variant in RAD51L1 (rs10483813) on 14q24.1 as we had hypothesized. In exploratory analyses, we found that the radiation-associated breast cancer risk varied significantly by linked markers in 5p12 (rs930395, rs10941679, rs2067980 and rs4415084) in the mitochondrial ribosomal protein S30 (MRPS30) gene (P(interaction) = 0.04). Chance, however, may explain these findings, and as such, these results need to be confirmed in other populations with low to moderate levels of radiation exposure. Even though a complete understanding of the way(s) in which these variants may increase breast cancer risk remains elusive, this approach may yield clues for further investigation.

Epidemiological and mechanistic evidence on the association of quercetin-rich food intake with lung cancer risk and carcinogenesis are inconclusive. We investigated the role of dietary quercetin and the interaction between quercetin and P450 and glutathione S-transferase (GST) polymorphisms on lung cancer risk in 1822 incident lung cancer cases and 1991 frequency-matched controls from the Environment And Genetics in Lung cancer Etiology study. In non-tumor lung tissue from 38 adenocarcinoma patients, we assessed the correlation between quercetin intake and messenger RNA expression of the same P450 and GST metabolic genes. Multivariate odds ratios (ORs) and 95% confidence intervals (CIs) for sex-specific quintiles of intake were calculated using unconditional logistic regression adjusting for putative risk factors. Frequent intake of quercetin-rich foods was inversely associated with lung cancer risk (OR = 0.49; 95% CI: 0.37–0.67; P-trend < 0.001) and did not differ by P450 or GST genotypes, gender or histological subtypes. The association was stronger in subjects who smoked >20 cigarettes per day (OR = 0.35; 95% CI: 0.19–0.66; P-trend = 0.003). Based on a two-sample t-test, we compared gene expression and high versus low consumption of quercetin-rich foods and observed an overall upregulation of GSTM1, GSTM2, GSTT2, and GSTP1 as well as a downregulation of specific P450 genes (P-values < 0.05, adjusted for age and smoking status). In conclusion, we observed an inverse association of quercetin-rich food with lung cancer risk and identified a possible mechanism of quercetin-related changes in the expression of genes involved in the metabolism of tobacco carcinogens in humans. Our findings suggest an interplay between quercetin intake, tobacco smoking, and lung cancer risk. Further research on this relationship is warranted.

Three genetic loci for lung cancer risk have been identified by genome-wide association studies (GWAS), but inherited susceptibility to specific histologic types of lung cancer is not well established. We conducted a GWAS of lung cancer and its major histologic types, genotyping 515,922 single-nucleotide polymorphisms (SNPs) in 5739 lung cancer cases and 5848 controls from one population-based case-control study and three cohort studies. Results were combined with summary data from ten additional studies, for a total of 13,300 cases and 19,666 controls of European descent. Four studies also provided histology data for replication, resulting in 3333 adenocarcinomas (AD), 2589 squamous cell carcinomas (SQ), and 1418 small cell carcinomas (SC). In analyses by histology, rs2736100 (TERT), on chromosome 5p15.33, was associated with risk of adenocarcinoma (odds ratio [OR] = 1.23, 95% confidence interval [CI] = 1.13–1.33, p = 3.02 × 10−7), but not with other histologic types (OR = 1.01, p = 0.84 and OR = 1.00, p = 0.93 for SQ and SC, respectively). This finding was confirmed in each replication study and overall meta-analysis (OR = 1.24, 95% CI = 1.17–1.31, p = 3.74 × 10−14 for AD; OR = 0.99, p = 0.69 and OR = 0.97, p = 0.48 for SQ and SC, respectively). Other previously reported association signals on 15q25 and 6p21 were also refined, but no additional loci reached genome-wide significance. In conclusion, a lung cancer GWAS identified a distinct hereditary contribution to adenocarcinoma.

© 2017, BMJ Publishing Group. All rights reserved. Objective: Although several genome-wide association studies (GWAS) of non-cardia gastric cancer have been published, more novel association signals could be exploited by combining individual studies together, which will further elucidate the genetic susceptibility of non-cardia gastric cancer. Design: We conducted a meta-analysis of two published Chinese GWAS studies (2031 non-cardia gastric cancer cases and 4970 cancer-free controls) and followed by genotyping of additional 3564 cases and 4637 controls in two stages. Results: The overall meta-analysis revealed two new association signals. The first was a novel locus at 5q14.3 and marked by rs7712641 (per-allele OR=0.84, 95% CI 0.80 to 0.88; p=1.21×10-11). This single-nucleotide polymorphism (SNP) marker maps to the intron of the long non-coding RNA, lnc-POLR3G-4 (XLOC-004464), which we observed has lower expression in non-cardia gastric tumour compared with matched normal tissue (Pwilcoxon signed-rank=7.20×10-4). We also identified a new signal at the 1q22 locus, rs80142782 (per-allele OR=0.62; 95% CI 0.56 to 0.69; p=1.71×10-19), which was independent of the previously reported SNP at the same locus, rs4072037 (per-allele OR=0.74; 95% CI 0.69 to 0.79; p=6.28×10-17). Analysis of the new SNP conditioned on the known SNP showed that the new SNP remained genome-wide significant (Pconditional=3.47×10-8). Interestingly, rs80142782 has a minor allele frequency of 0.05 in East Asians but is monomorphic in both European and African populations. Conclusion: These findings add new evidence for inherited genetic susceptibility to non-cardia gastric cancer and provide further clues to its aetiology in the Han Chinese population.

Interpretation of biological mechanisms underlying genetic risk associations for prostate cancer is complicated by the relatively large number of risk variants (n = 100) and the thousands of surrogate SNPs in linkage disequilibrium. Here, we combined three distinct approaches: multiethnic fine-mapping, putative functional annotation (based upon epigenetic data and genomeencoded features), and expression quantitative trait loci (eQTL) analyses, in an attempt to reduce this complexity. We examined 67 risk regions using genotyping and imputation-based fine-mapping in populations of European (cases/controls: 8600/6946), African (cases/controls: 5327/5136), Japanese (cases/controls: 2563/4391) and Latino (cases/controls: 1034/1046) ancestry. Markers at 55 regions passed a region-specific significance threshold (P-value cutoff range: 3.9 × 10−4 –5.6 × 10−3 ) and in 30 regions 5604 | Human Molecular Genetics, 2015, Vol. 24, No. 19 we identified markers that were more significantly associated with risk than the previously reported variants in the multiethnic sample. Novel secondary signals (P < 5.0 × 10−6 ) were also detected in two regions (rs13062436/3q21 and rs17181170/3p12). Among 666 variants in the 55 regions with P-values within one order of magnitude of the most-associated marker, 193 variants (29%) in 48 regions overlapped with epigenetic or other putative functional marks. In 11 of the 55 regions, cis-eQTLs were detected with nearby genes. For 12 of the 55 regions (22%), the most significant region-specific, prostate-cancer associated variant represented the strongest candidate functional variant based on our annotations; the number of regions increased to 20 (36%) and 27 (49%) when examining the 2 and 3 most significantly associated variants in each region, respectively. These results have prioritized subsets of candidate variants for downstream functional evaluation.

Waldenström macroglobulinemia (WM) is a unique subset of lymphoplasmacytic lymphoma (LPL) that is defined by the presence of an LPL infiltrate in the bone marrow together with a monoclonal IgM protein in the serum. A somatic activating mutation, MYD88 L265P, occurs in 85+% of WM and in 25%-50% of patients with the precursor condition, IgM monoclonal gammopathy of undetermined significance (MGUS); however, germline MYD88 mutations have not been observed in WM patients, and the genetic basis for WM predisposition remains undefined. To identify novel WM susceptibility loci we conducted a two-stage genome-wide association study (GWAS) in over 450 WM cases and 4300 controls of European ancestry. Discovery (stage 1) included 217 WM cases (40% familial) and 3798 controls genotyped on the Illumina Omni Express or Illumina Omni2.5 platforms following standard quality control procedures. The genotyped data were imputed using the Haplotype Reference Consortium panel as a reference and analyzed using logistic regression. In stage 1, we identified three loci on chromosomes 6, 14 and 3 significantly associated (P Citation Format: Mary L. McMaster, Sonja I. Berndt, Shengchao A. Li, Susan Slager, Joseph Vijai, Charles C. Chung, Bin Zhu, Laurie Burdette, Brenda Birmann, Elizabeth E. Brown, James R. Cerhan, Karin Ekstrom-Smedby, Henrik Hjalgrim, Geffen Kleinstern, Brian K. Link, James McKay, Alain Monnereau, Lindsay M. Morton, Alexandra Nieters, Nathaniel Rothman, Christine F. Skibola, Alex Smith, Lauren R. Teras, Claire M. Vajdic, Roel Vermeulen, Belynda Hicks, Lynn R. Goldin, Neil E. Caporaso. A genome-wide association study of Waldenström macroglobulinemia/lymphoplasmacytic lymphoma demonstrates association with chromosome 6 [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1318. doi:10.1158/1538-7445.AM2017-1318

Human papillomavirus (HPV) is a very common sexually transmitted infection, however only a small proportion of women progress to cervical precancer or cancer. HPV16 is the most carcinogenic type, causing more than half of the cervical cancer globally. The HPV16 genome is 7,906base-pairs coding for 8 genes (E6, E7, E1, E2, E4, E5, L2, L1) and one upstream regulatory region (URR). Within HPV16 there are 4 main lineages (A, B, C, D) that are strongly associated with disease risk. Human APOBEC3A (hA3A) cytidine deaminases have been shown to have antiviral effects. The APOBEC mutational process results in a C to T base change at specific motifs (5’ [C/T]•C>T•W 3’). Previous studies established that there was evidence of APOBEC3 editing on a small number of HPV16 samples. It’s unknown how these mutations are related to infection clearance or the long-term accumulation of genomic mutations that contribute to HPV-associated cancers. We conducted detailed analyses to comprehensively evaluate APOBEC3 editing on HPV16 genomes using HPV16 whole-genome sequencing data from 3,215 HPV16-infected women in the NCI-HPV Persistence and Progression (PaP) cohort. Cases were defined as women with cervical precancer (CIN3, N = 1,093) or cancer (N = 109) and controls were women with no histologic evidence of precancer or cancer (T•W 3’) out of 96 potential motifs of 3 base-pairs. Using logistic regression, we compared the number of APOBEC3-associated variants in cases and controls, among HPV16 variant lineages, and among genome regions of the virus. We discovered that there is evidence of APOBEC3 editing throughout the HPV16 genome. Specifically, we observed that women with precancer or cancer had less APOBEC3-associated variants compared to the controls (OR = 0.84, p-value = 0.06). We further showed that women with an HPV16 A lineage infection have more APOBEC3-associated variants compared to those with a non-A lineage infection (OR = 1.35, p-value = 0.02). After controlling for the number of APOBEC3 vulnerable loci, we observed that the L1 (OR = 0.23, p-value = 0.04) and E7 (OR = 0.29, p-value = 0.07) genes have less APOBEC3 footprints overall, and particularly in the cases compared to the controls in these regions, compared to the viral non-coding upstream regulatory region (URR). Overall, we found that APOBEC3 is not affecting the HPV16 genome in a uniform way, and instead, it appears to be targeting specific regions which could suggest antiviral activity. Importantly, we determined that APOBEC3-associated variants are less prevalent in cases which could be related to disease progression in these individuals. Further evaluation is underway. Citation Format: Yanzi Xiao, Bin Zhu, Meredith Yeager, Michael Cullen, Joseph Boland, Nicolas Wentzensen, Tina Raine-Bennett, Zigui Chen, Kai Yu, Qi Yang, Mia Steinberg, David Roberson, Sara Bass, Laurie Burdette, Thomas Lorey, Philip Castle, Robert Burk, Mark Schiffman, Lisa Mirabello. Evidence of APOBEC3 editing in the HPV16 genome [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 4249. doi:10.1158/1538-7445.AM2017-4249

Candidate gene and genome-wide association studies (GWAS) have identified 11 independent susceptibility loci associated with bladder cancer risk. To discover additional risk variants, we conducted a new GWAS of 2422 bladder cancer cases and 5751 controls, followed by a meta-analysis with two independently published bladder cancer GWAS, resulting in a combined analysis of 6911 cases and 11 814 controls of European descent. TaqMan genotyping of 13 promising single nucleotide polymorphisms with P < 1 × 10(-5) was pursued in a follow-up set of 801 cases and 1307 controls. Two new loci achieved genome-wide statistical significance: rs10936599 on 3q26.2 (P = 4.53 × 10(-9)) and rs907611 on 11p15.5 (P = 4.11 × 10(-8)). Two notable loci were also identified that approached genome-wide statistical significance: rs6104690 on 20p12.2 (P = 7.13 × 10(-7)) and rs4510656 on 6p22.3 (P = 6.98 × 10(-7)); these require further studies for confirmation. In conclusion, our study has identified new susceptibility alleles for bladder cancer risk that require fine-mapping and laboratory investigation, which could further understanding into the biological underpinnings of bladder carcinogenesis.

Background: Vitamin D and its metabolites are believed to impede carcinogenesis by stimulating cell differentiation, inhibiting cell proliferation, and inducing apoptosis. Certain pesticides have been shown to deregulate vitamin D's anticarcinogenic properties. We hypothesize that certain pesticides may be linked to prostate cancer via an interaction with vitamin D genetic variants. Methods: We evaluated interactions between 41 pesticides and 152 single-nucleotide polymorphisms (SNP) in nine vitamin D pathway genes among 776 prostate cancer cases and 1,444 male controls in a nested case–control study of Caucasian pesticide applicators within the Agricultural Health Study. We assessed Pinteraction values using likelihood ratio tests from unconditional logistic regression and a false discovery rate (FDR) to account for multiple comparisons. Results: Five significant interactions (P < 0.01) displayed a monotonic increase in prostate cancer risk with individual pesticide use in one genotype and no association in the other. These interactions involved parathion and terbufos use and three vitamin D genes (VDR, RXRB, and GC). The exposure–response pattern among participants with increasing parathion use with the homozygous CC genotype for GC rs7041 compared with unexposed participants was noteworthy [low vs. no exposure: OR, 2.58, 95% confidence interval (CI), 1.07–6.25; high vs. no exposure: OR, 3.09, 95% CI, 1.10–8.68; Pinteraction = 3.8 × 10−3]. Conclusions: In this study, genetic variations in vitamin D pathway genes, particularly GC rs7041, an SNP previously linked to lower circulating vitamin D levels, modified pesticide associations with prostate cancer risk. Impact: Because our study is the first to examine this relationship, additional studies are needed to rule out chance findings. Cancer Epidemiol Biomarkers Prev; 22(9); 1557–66. ©2013 AACR.

INTRODUCTION: Recent evidence suggests a link between constitutional telomere length (TL) and cancer risk. Previous studies have suggested that longer telomeres were associated with an increased risk of melanoma and larger size and number of nevi. The goal of this study was to examine whether TL modified the risk of melanoma in melanoma-prone families with and without CDKN2A germline mutations. MATERIALS AND METHODS: We measured TL in blood DNA in 119 cutaneous malignant melanoma (CMM) cases and 208 unaffected individuals. We also genotyped 13 tagging SNPs in TERT. RESULTS: We found that longer telomeres were associated with an increased risk of CMM (adjusted OR = 2.81, 95% CI = 1.02-7.72, P = 0.04). The association of longer TL with CMM risk was seen in CDKN2A- cases but not in CDKN2A+ cases. Among CMM cases, the presence of solar injury was associated with shorter telomeres (P = 0.002). One SNP in TERT, rs2735940, was significantly associated with TL (P = 0.002) after Bonferroni correction. DISCUSSION: Our findings suggest that TL regulation could be variable by CDKN2A mutation status, sun exposure, and pigmentation phenotype. Therefore, TL measurement alone may not be a good marker for predicting CMM risk.

Background: The chromosome 9p21.3 region has been implicated in the pathogenesis of multiple cancers. Methods: We systematically examined up to 203 tagging SNPs of 22 genes on 9p21.3 (19.9–32.8 Mb) in eight case–control studies: thyroid cancer, endometrial cancer (EC), renal cell carcinoma, colorectal cancer (CRC), colorectal adenoma (CA), oesophageal squamous cell carcinoma (ESCC), gastric cardia adenocarcinoma and osteosarcoma (OS). We used logistic regression to perform single SNP analyses for each study separately, adjusting for study-specific covariates. We combined SNP results across studies by fixed-effect meta-analyses and a newly developed subset-based statistical approach (ASSET). Gene-based P-values were obtained by the minP method using the Adaptive Rank Truncated Product program. We adjusted for multiple comparisons by Bonferroni correction. Results: Rs3731239 in cyclin-dependent kinase inhibitors 2A (CDKN2A) was significantly associated with ESCC (P=7 × 10−6). The CDKN2A-ESCC association was further supported by gene-based analyses (Pgene=0.0001). In the meta-analyses by ASSET, four SNPs (rs3731239 in CDKN2A, rs615552 and rs573687 in CDKN2B and rs564398 in CDKN2BAS) showed significant associations with ESCC and EC (P

Genome-wide association studies (GWAS) have previously identified 13 loci associated with risk of chronic lymphocytic leukemia or small lymphocytic lymphoma (CLL). To identify additional CLL susceptibility loci, we conducted the largest meta-analysis for CLL thus far, including four GWAS with a total of 3,100 individuals with CLL (cases) and 7,667 controls. In the meta-analysis, we identified ten independent associated SNPs in nine new loci at 10q23.31 (ACTA2 or FAS (ACTA2/FAS), P = 1.22 × 10-14), 18q21.33 (BCL2, P = 7.76 × 10-11), 11p15.5 (C11orf21, P = 2.15 × 10 -10), 4q25 (LEF1, P = 4.24 × 10-10), 2q33.1 (CASP10 or CASP8 (CASP10/CASP8), P = 2.50 × 10-9), 9p21.3 (CDKN2B-AS1, P = 1.27 × 10-8), 18q21.32 (PMAIP1, P = 2.51 × 10 -8), 15q15.1 (BMF, P = 2.71 × 10-10) and 2p22.2 (QPCT, P = 1.68 × 10-8), as well as an independent signal at an established locus (2q13, ACOXL, P = 2.08 × 10-18). We also found evidence for two additional promising loci below genome-wide significance at 8q22.3 (ODF1, P = 5.40 × 10-8) and 5p15.33 (TERT, P = 1.92 × 10-7). Although further studies are required, the proximity of several of these loci to genes involved in apoptosis suggests a plausible underlying biological mechanism. © 2013 Nature America, Inc. All rights reserved.

Background: Folate is essential for one-carbon metabolism, a pathway required by DNA synthesis, methylation, and repair. Low dietary and circulating folate and polymorphic variation in this pathway are associated with increased risk of colorectal adenoma and cancer. Methods: We genotyped 882 single nucleotide polymorphisms (SNP) in 82 one-carbon metabolism genes for 1,331 cases of advanced colorectal adenoma, identified by sigmoidoscopy at baseline, and 1,501 controls from the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial (PLCO). We evaluated associations between one-carbon genes and adenoma risk in all subjects and stratified by folate intake. We applied the Adaptive Rank Truncated Product (ARTP) method to assess statistical significance at the gene and pathway levels. Results: Folate intake was inversely associated with advanced colorectal adenoma risk [odds ratio (OR) by quartile = 0.85, P = 1.9 × 10−5]. We found no statistically significant associations between one-carbon genes and adenoma risk in all subjects. As hypothesized, we observed a statistically significant pathway-level association (P = 0.038) in the lowest quartile of folate; no significant associations were found in higher quartiles. Several genes including adenosine deaminase (ADA) and cysteine dioxygenase (CDO1) contributed to this signal (gene-level P = 0.001 and 0.0073, respectively). The most statistically significant SNP was rs244072 in ADA (P = 2.37 × 10−5). Conclusions and Impact: Stratification by dietary folate and application of the ARTP method revealed statistically significant pathway- and gene-level associations between one-carbon metabolism genes and risk of advanced colorectal adenoma, which were not apparent in analysis of the entire population. Folate intake may interact with associations between common variants in one-carbon metabolism genes and colorectal adenoma risk. Cancer Epidemiol Biomarkers Prev; 21(3); 417–27. ©2012 AACR.

Previous research demonstrates increased prostate cancer risk for pesticide applicators and pesticide manufacturing workers. Although underlying mechanisms are unknown, human biomonitoring studies indicate increased genetic damage (e.g. chromosomal aberrations) with pesticide exposure. Given that the nucleotide excision repair (NER) pathway repairs a broad range of DNA damage, we evaluated interactions between pesticide exposure and 324 single-nucleotide polymorphisms (SNPs) tagging 27 NER genes among 776 prostate cancer cases and 1444 male controls in a nested case-control study of white Agricultural Health Study pesticide applicators. We determined interaction P values using likelihood ratio tests from logistic regression models and three-level pesticide variables (none/low/high) based on lifetime days of use weighted to an intensity score. We adjusted for multiple comparisons using the false discovery rate (FDR) method. Of the 17 interactions that met FDR

Background: Host genetic factors might affect the risk of progression from infection with carcinogenic human papillomavirus (HPV), the etiologic agent for cervical cancer, to persistent HPV infection, and hence to cervical precancer and cancer. Methodology/principal findings: We assessed 18,310 tag single nucleotide polymorphisms (SNPs) from 1113 genes in 416 cervical intraepithelial neoplasia 3 (CIN3)/cancer cases, 356 women with persistent carcinogenic HPV infection (median persistence of 25 months) and 425 randomly selected women (non-cases and non-HPV persistent) from the 10,049 women from the Guanacaste, Costa Rica HPV natural history cohort. For gene and SNP associations, we computed age-adjusted odds ratio and p-trend. Three comparisons were made: 1) association with CIN3/cancer (compared CIN3/cancer cases to random controls), 2) association with persistence (compared HPV persistence to random controls), and 3) progression (compared CIN3/cancers with HPV-persistent group). Regions statistically significantly associated with CIN3/cancer included genes for peroxiredoxin 3 PRDX3, and ribosomal protein S19 RPS19. The single most significant SNPs from each gene associated with CIN3/cancer were PRDX3 rs7082598 (P(trend)

To explore associations with prostate cancer and farming, it is important to investigate the relationship between pesticide use and single nucleotide polymorphisms (SNPs) in xenobiotic metabolic enzyme (XME) genes.[corrected] We evaluated pesticide-SNP interactions between 45 pesticides and 1913 XME SNPs with respect to prostrate cancer among 776 cases and 1444 controls in the Agricultural Health Study.We used unconditional logistic regression to estimate odds ratios (ORs) and 95% confidence intervals (CIs). Multiplicative SNP-pesticide interactions were calculated using a likelihood ratio test.A positive monotonic interaction was observed between petroleum oil/petroleum distillate use and rs1883633 in the oxidative stress gene glutamate cysteine ligase (GCLC; P interaction=1.0×10(-4)); men carrying at least one variant allele (minor allele) experienced an increased prostate cancer risk (OR=3.7, 95% CI: 1.9-7.3). Among men carrying the variant allele for thioredoxin reductase 2 (TXNRD2) rs4485648, microsomal epoxide hydrolase 1 (EPHX1) rs17309872, or myeloperoxidase (MPO) rs11079344, an increased prostate cancer risk was observed with high, compared with no, petroleum oil/petroleum distillate (OR=1.9, 95% CI: 1.1-3.2, P interaction=0.01; OR=2.1, 95% CI: 1.1-4.0, P interaction=0.01), or terbufos (OR=3.0, 95% CI: 1.5-6.0, P interaction=2.0×10(-3)) use, respectively. No interactions were deemed noteworthy at the false discovery rate=0.20 level; the number of observed interactions in XMEs was comparable with the number expected by chance alone.We observed several pesticide-SNP interactions in oxidative stress and phase I/II enzyme genes and risk of prostate cancer. Additional work is needed to explain the joint contribution of genetic variation in XMEs, pesticide use, and prostate cancer risk.

Chromosome 9p21 has been implicated in the pathogenesis of cutaneous malignant melanoma (CMM). In addition to CDKN2A, the major known high-risk susceptibility gene for CMM, recent studies suggest that other 9p21 genes may be involved in melanoma/nevi development. To identify 9p21 variants that influence susceptibility to CMM and number of nevi in CMM-prone families with and without CDKN2A mutations, we analyzed 562 individuals (183 CMM) from 53 families (23 CDKN2A+, 30 CDKN2A-) for 233 tagging SNPs in 21 genes at 9p21. Single SNP- and gene-based regression analyses were used to assess the risk of CMM, nevi count, skin complexion, and tanning ability associated with these SNPs and genes. We found that SNP rs7023329 in the MTAP gene was associated with number of nevi (P (trend) = 0.003) confirming a recent finding by a genome-wide association study. In addition, three SNPs in the ACO1 gene, rs7855483 (P (trend) = 0.002), rs17288067 (P (trend) = 0.0009), and rs10813813 (P (trend) = 0.005), showed the strongest associations with CMM risk. None of the examined 9p21 SNPs was associated with skin complexion, whereas two SNPs, rs10964862 in IFNW1 (P (trend) = 0.003), and rs13290968 in TUSC1 (P (trend) = 0.0006), were associated with tanning ability. Gene-based analyses suggested that the ACO1 gene was significantly associated with CMM (P = 0.0004); genes IFNW1 (P = 0.002) and ACO1 (P = 0.0002) were significantly associated with tanning ability. Our findings are consistent with recent proposals that additional 9p21 genes may contribute to CMM susceptibility in CMM-prone families. These genetic variants may, at least partially, exert their effects through nevi and tanning ability.

Background. Estrogen plays a major role in endometrial carcinogenesis, suggesting that common variants of genes in the sex hormone metabolic pathway may be related to endometrial cancer risk. In support of this view, variants in CYP19A1 [cytochrome P450 (CYP), family 19, subfamily A, polypeptide 1] have been associated with both circulating estrogen levels and endometrial cancer risk. Associations with variants in other genes have been suggested, but findings have been inconsistent. Methods. We examined 36 sex hormone-related genes using a tagging approach in a population-based case–control study of 417 endometrial cancer cases and 407 controls conducted in Poland. We evaluated common variation in these genes in relation to endometrial cancer risk using sequential haplotype scan, variable-sized sliding window and adaptive rank-truncated product (ARTP) methods. Results. In our case–control study, the strongest association with endometrial cancer risk was for AR (androgen receptor; ARTP P = 0.006). Multilocus analyses also identified boundaries for a region of interest in AR and in CYP19A1 around a previously identified susceptibility loci. We did not find evidence for consistent associations between previously reported candidate single-nucleotide polymorphisms in this pathway and endometrial cancer risk. Discussion. In summary, we identified regions in AR and CYP19A1 that are of interest for further evaluation in relation to endometrial cancer risk in future haplotype and subsequent fine mapping studies in larger study populations.

Cutaneous malignant melanoma (CMM) is an etiologically heterogeneous disease with genetic, host, environmental factors, and their interactions contributing to its development. CDKN2A and CDK4 are the two established major susceptibility genes for melanoma identified so far. Recent evidence suggests that copy number variations (CNVs) may contribute to disease susceptibility in several inherited diseases including cancer. The goals of this study were 1) to assess whether the frequency of CNVs varied by CMM or CDKN2A/CDK4 mutation status; and 2) to identify rare CNVs that were related to melanoma predisposition in these high-risk families. We used genome-wide tiling CGH arrays (Nimblegen 720K exon-focused) to investigate characteristics of CNVs in 174 CMM cases, 44 high-risk unaffected family members (with dysplastic nevi/large number of moles or germline CDKN2A mutations), and 48 unrelated spouses from 50 American melanoma-prone families (21 CDKN2A+, 2 CDK4+, 27 mutation negative). We used the Nexus Copy Number™ built-in FASST2 algorithm to identify significant CNVs (significant threshold = 0.000001; minimal number of probes per segment = 5; log2 ratio>0.3 for gains and -0.3 for losses). We found that the median number of total CNVs, or gains or losses separately, did not show significant differences in CMM cases, high-risk unaffected family members, and unrelated controls. Among CMM cases, CNV frequencies were not significantly associated with germline CDKN2A/CDK4 mutation status, age at melanoma diagnosis, or number of melanomas. Similar results were obtained when number of genes and lengths of DNA segments affected by CNVs were analyzed. Restriction to large CNVs (>10 kb or >100 kb) or rare CNVs (not reported in the Toronto CNV database) did not change results significantly. On the other hand, we identified several rare large CNVs (>10 kb, not reported in unrelated controls) that either involved known melanoma genes or co-segregated with melanoma (observed in multiple CMM cases) within families. These included a 1.3Mb deletion in PARP1 in a single CMM case, a 10 kb deletion in CDKN2A in 4 of 5 CMM cases and two obligate gene carriers in a large family that was negative for CDKN2A mutations by sequencing, a 175 kb deletion in LINGO2 in 2 of 3 CMM cases in one family, a 10 kb deletion in 8q24 in 3 of 4 CMM cases in one family, a 110 kb deletion in 2q22.1 in all three cases in a family, and a 57 kb duplication in 4q32.2 in all five cases in a CDKN2A mutation positive family. The role of these CNVs, particularly those involving genes that have unknown function related to CMM development, in CMM susceptibility remain to be investigated. Citation Format: Xiaohong (Rose) Yang, Jianxin Shi, Hunter Bennett, Laura Burke, Casey Dagnall, Laurie Burdette, Belynda Hicks, Margaret Tucker, Alisa Goldstein. Germline copy number variations in melanoma families with/without CDKN2A/CDK4 mutations. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2755. doi:10.1158/1538-7445.AM2015-2755

In spite of progress in the identification of the somatic events associated with Waldenström macroglobulinemia (WM), the genetic determinants of WM susceptibility have not been identified. We have conducted exome sequencing in 32 individuals from 9 well-characterized families at high risk for WM. We sequenced 3 or more patients or obligate carriers from each family. We used Nimblegen v2.0 and v3.0 for exome capture followed by sequencing on the Illumina HiSeq2000. Among quality control measures, we required 80% of coding sequences to achieve 15x coverage. To eliminate other possible sequencing artifacts, variants found in more than 1% of samples sequenced in our laboratory from other studies were excluded. Novoalign v.2.07.14 was used for alignment and GATK was used for local re-alignment and variant calling. A large number of rare ( Citation Format: Mary L. McMaster, Lynn R. Goldin, Melissa Rotunno, Ji He, Laurie Burdette, Amy Hutchinson, Joseph Boland, Meredith Yeager, Margaret A. Tucker, Stephen J. Chanock, Neil E. Caporaso. Exploration of rare variants from exome sequencing in families with Waldenstrom macroglobulinemia (WM). [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 1300. doi:10.1158/1538-7445.AM2014-1300

Genome wide association studies (GWAS) are an important tool for identifying complex human disease loci. A key variable for GWAS is the degree of admixture or recent mixing of previously isolated populations. Undetected admixture can lead to spurious associations, although it may also be used as a strategy for mapping by admixture linkage disequilibrium (MALD). Development of reagents to identify admixture and to perform MALD may help to identify genetic modifiers of sickle cell disease (SCD). We have examined ancestry and genetic admixture in SCD patients living in North America. We first examined ancestry from family histories for 649 adults recruited to the Bethesda Sickle Cell Cohort Study in the eastern US, where 61% of subjects are African American, 22% African, 12% Caribbean or South American and 1.8% of other origins. We hypothesized that the high proportion of subjects from the Africa ancestry group would diminish the expected degree of admixture in the Bethesda cohort compared to other SCD cohorts. To assess admixture at the genetic level, we identified ancestry informative markers (AIMs) from 3,804,602 SNPs genotyped by HapMap Phase III in CEU (European) and YRI (Nigeria) populations. These SNPs were combined with a published admixture mapping panel, yielding a new MALD panel of 2251 AIMs. We assessed the performance of our AIM panel by genotyping 221 HapMap individuals and 489 adults from the Bethesda cohort with an overall completion rate of 98.51% using Illumina iSelect arrays. HapMap samples had 98% concordance with publicly available data. Out of 2251 candidate markers, only 1806 successfully passed all quality controls (289 SNPs could not be genotyped; 82 were not AIMs; 74 had low concordance with HapMap). Using our HapMap genotype data, this admixture panel has a mean δ (difference in allele frequencies between populations) of 0.712 (SD 0.118) with an average inter-marker distance of 1.862 Mb. The utility of this 1806 marker panel for genome wide MALD in SCD was demonstrated with a proof of concept scan to map the known genetic locus underlying SCD in the Bethesda cohort. To compare admixture in the Bethesda cohort to other SCD populations, we defined a subset of AIMs (n=360) that distinguish 6 major geo-ethnic groups from the Human Genome Diversity Panel (n=952 individuals). These results suggest that these 360 markers can distinguish contributions from geo-ethnic ancestry groups besides African and European populations. Using these markers in the Bethesda cohort, principle components analysis (PCA) showed African American or Caribbean/South American subjects had a wider range of admixture compared to Africans or HapMap populations (CEU, YRI and CHB). We then compared admixture proportions in the Bethesda cohort to 2 additional SCD populations, including 469 anonymous samples from a newborn screen SCD cohort in the western US and 439 SCD adults from a clinical trial in North America and Europe (WalkPhasst). Comparison of these SCD cohorts to 3 HapMap populations using PCA showed nearly 50% of the variation is explained by 2 major vectors that distinguish European/African and Asian/Native American ancestry. Finally, we compared admixture proportions across all 3 SCD cohorts using Structure. Here, we observed significant differences in ancestry proportions arising from Africa, Europe and Asia/Americas (P Disclosures: No relevant conflicts of interest to declare.

Chromosome 9p21 has been implicated in the pathogenesis of cutaneous malignant melanoma (CMM), glioma, non-small cell lung cancer, and leukemia. This region includes cyclin-dependent kinase inhibitor 2A (CDKN2A), the major known high-risk susceptibility gene for CMM; CDKN2B; antisense noncoding RNA in the INK4 locus (ANRIL), a GWAS hotspot for multiple phenotypes including cancers; methylthioadenosine phosphorylase (MTAP), a gene with tumor suppressor function and a cluster of type I interferon (IFN) genes. The goal of this study was to examine the association between genetic variation in the 9p21 region and multiple cancer outcomes using genotyping data from studies participating in the National Cancer Institute's (NCI) iSelect project. We evaluated up to 202 tagSNPs in 22 genes in the 9p21 region (19.9-32.8 Mb), in case-control studies of thyroid cancer (THC), endometrial cancer (EC), testicular cancer (TC), renal cell carcinoma (RCC), colorectal cancer (CRC), colorectal adenoma (CA), esophageal squamous cell carcinoma (ESCC), gastric cardia cancer (GCC) and osteosarcoma (OS). The number of cases in each of the nine studies ranged from 96 to 1234. We first performed single SNP-based analyses for each study separately, using logistic regression adjusted for age and study-specific covariates. To combine SNP based analyses, we used a newly developed subset-based statistical approach (ASSET), which allow for heterogeneity of SNP effect on different outcomes. This approach exhaustively explores subsets of studies for the presence of association signals, and then evaluates the significance from subsets with the strongest association signal, meanwhile accounting for multiple tests required by the subset search. Gene based P-values were computed using an adaptive rank truncated product (ARTP) statistic and a permutation-based sampling procedure (10,000 permutations) that took into account the numbers of SNPs in each gene and their linkage disequilibrium structure. When analyzing each study separately, we found that a SNP (rs3731257) in CDKN2A was significantly associated with the risk of ESCC after Bonferroni correction for number of SNPs and studies (P=7 x 10-6). Meta-analyses by ASSET found that this SNP was significantly associated with both ESCC and EC (PASSET Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr LB-330. doi:1538-7445.AM2012-LB-330

The population of Costa Rica (CR) represents an admixture of major continental populations. An investigation of the CR population structure would provide an important foundation for mapping genetic variants underlying common diseases and traits. We conducted an analysis of 1,301 women from the Guanacaste region of CR using 27,904 single nucleotide polymorphisms (SNPs) genotyped on a custom Illumina InfiniumII iSelect chip. The program STRUCTURE was used to compare the CR Guanacaste sample with four continental reference samples, including HapMap Europeans (CEU), East Asians (JPT+CHB), West African Yoruba (YRI), as well as Native Americans (NA) from the Illumina iControl database. Our results show that the CR Guanacaste sample comprises a three-way admixture estimated to be 43% European, 38% Native American and 15% West African. An estimated 4% residual Asian ancestry may be within the error range. Results from principal components analysis reveal a correlation between genetic and geographic distance. The magnitude of linkage disequilibrium (LD) measured by the number of tagging SNPs required to cover the same region in the genome in the CR Guanacaste sample appeared to be weaker than that observed in CEU, JPT+CHB and NA reference samples but stronger than that of the HapMap YRI sample. Based on the clustering pattern observed in both STRUCTURE and principal components analysis, two subpopulations were identified that differ by approximately 20% in LD block size averaged over all LD blocks identified by Haploview. We also show in a simulated association study conducted within the two subpopulations, that the failure to account for population stratification (PS) could lead to a noticeable inflation in the false positive rate. However, we further demonstrate that existing PS adjustment approaches can reduce the inflation to an acceptable level for gene discovery.

Background: Quercetin is a flavonol found in fruits (apple, grape) and vegetables (onion, broccoli). Quercetin possesses anticarcinogenic properties that protect against cancer. However, the epidemiological evidence on the association of quercetin-rich food intake with lung cancer risk is inconclusive, and how quercetin affects lung carcinogenesis in humans is largely unknown. Laboratory analyses have shown that quercetin may inhibit Phase I cytochrome P450s enzymes responsible for the activation of tobacco-related procarcinogens, and induce Phase II glutathione S-transferase (GST) enzymes responsible for the removal of activated carcinogens. Methods: We investigated the role of dietary quercetin and the interaction between quercetin and P450 and GST polymorphisms on lung cancer risk in incident lung cancer cases and matched controls from Environment and Genetics in Lung cancer Etiology (EAGLE), a population-based case-control study including 2100 primary lung cancer cases and 2120 controls. Dietary information and genotype data were available for 1740 cases and 1919 controls. In fresh-frozen non-tumor lung tissue from 38 adenocarcinoma patients, we also assessed the correlation between dietary quercetin intake and mRNA expression of tobacco-related Phase I and Phase II metabolic genes. Multivariable odds ratios (ORs) and 95% confidence intervals (CIs) for sex-specific quintile of intake were calculated using unconditional logistic regression, adjusted for putative risk factors. Two sample t-test was used to assess difference in expression by quercetin-rich consumption status. Results: Frequent intake of quercetin-rich foods was inversely associated with lung cancer risk (ORQ5vsQ1=0.49;95% CI=0.37–0.67; p-trend: Conclusions: Quercetin-rich diet was associated with lower risk of lung cancer and with down-regulation of Phase I P450 genes and up-regulation of Phase II GST genes in lung tissue. These findings suggest a possible mechanism for quercetin-related protection against lung cancer risk through quercetin-induced changes in the expression of genes involved in the metabolism of tobacco-carcinogens. Further studies exploring the interplay between quercetin intake, tobacco smoking, and lung cancer risk are warranted. Citation Information: Cancer Prev Res 2010;3(1 Suppl):A105.

Polymorphisms in genes coding for enzymes that activate tobacco lung carcinogens may generate inter-individual differences in lung cancer risk. Previous studies had limited sample sizes, poor exposure characterization, and a few single nucleotide polymorphisms (SNPs) tested in candidate genes. We analyzed 25 SNPs (some previously untested) in 2101 primary lung cancer cases and 2120 population controls from the Environment And Genetics in Lung cancer Etiology (EAGLE) study from six phase I metabolic genes, including cytochrome P450s, microsomal epoxide hydrolase, and myeloperoxidase. We evaluated the main genotype effects and genotype-smoking interactions in lung cancer risk overall and in the major histology subtypes. We tested the combined effect of multiple SNPs on lung cancer risk and on gene expression. Findings were prioritized based on significance thresholds and consistency across different analyses, and accounted for multiple testing and prior knowledge. Two haplotypes in EPHX1 were significantly associated with lung cancer risk in the overall population. In addition, CYP1B1 and CYP2A6 polymorphisms were inversely associated with adenocarcinoma and squamous cell carcinoma risk, respectively. Moreover, the association between CYP1A1 rs2606345 genotype and lung cancer was significantly modified by intensity of cigarette smoking, suggesting an underlying dose-response mechanism. Finally, increasing number of variants at CYP1A1/A2 genes revealed significant protection in never smokers and risk in ever smokers. Results were supported by differential gene expression in non-tumor lung tissue samples with down-regulation of CYP1A1 in never smokers and up-regulation in smokers from CYP1A1/A2 SNPs. The significant haplotype associations emphasize that the effect of multiple SNPs may be important despite null single SNP-associations, and warrants consideration in genome-wide association studies (GWAS). Our findings emphasize the necessity of post-GWAS fine mapping and SNP functional assessment to further elucidate cancer risk associations.

The tumor suppressor p53 plays an essential role in the prevention of tumor development by its ability to respond to stress signals, resulting in the inhibition of cell growth, via cell-cycle arrest or induction of apoptosis. p53 also contributes to the repair of genotoxic DNA damages and is involved in cellular senescence. Previous studies have shown that both germ-line and somatic mutations in the TP53 gene, which encodes p53, can occur in a wide spectrum of cancers, including acute lymphoblastic leukemia, ALL. Recently, we characterized the common germ-line genetic variation across TP53 by extensive re-sequence analysis (http://snp500cancer.nci.nih.gov). Patients and Methods: We conducted a genetic association study in childhood ALL to investigate the possible contribution of common single nucleotide polymorphisms, SNPs, and the common haplotypes on which they reside. For TP53, we analyzed haplotype-tagging SNPs, which are estimated to capture greater than 95% of common haplotypes across the gene (IVS2+38C, Ex4+119 (R72P), IVS4-91, Ex6-34 (R213R), IVS6+62, IVS6-36 and 1846bp 3′ of STP). Cases of ALL (n = 513) were drawn from anonymized samples from the German co-operative therapy study (COALL-06-97) and 751 healthy anonymized blood donor controls collected from regional blood banks in Germany. Genotype analysis was performed by validated Taqman assays posted on the SNP500cancer website. Haplotypes were deduced by maximization-estimation analysis of unphased genotypes using PHASEv2.0.1 (http://www.stat.washington.edu/stephens/software.html). Results: The single locus analysis revealed an increased risk for ALL for four of the seven SNPs individually: C allele at IVS2+38C (304/982 in patients vs. 399/1476 in controls, OR 1.21 95% CI 1.01 – 1.45), G allele at IVS4-91 (150/988 vs. 168/1456, OR 1.37 95% CI 1.08 – 1.75), A allele at Ex6-34 (R213R) (164/996 vs. 178/1450, OR 1.41 95% CI 1.11 – 1.78) and T allele at 1846bp 3′ of STP (185/990 vs. 203/1420, OR 1.38 95% CI 1.10 – 1.72). Individuals homozygous AA at IVS6+62 (16/350 vs. 9/556, OR 2.82, 95% CI 1.16 – 6.99) or TT at 1846bp 3′ of STP (24/334 vs. 15/522 OR 2.50, 955 CI 1.24 – 5.09) showed further risk for disease. The non-synonymous coding SNP at Ex4+119 (R72P) did not show any association with ALL. The distribution of the inferred haplotypes differed significantly between cases and controls (global haplotype test p = 0.02). Children carrying the haplotype CCGAAGT have an increased risk for ALL (138/611 vs. 155/972, OR 1.42, 95% CI 1.09 −1.83) in reference to the most common haplotype GGAAGGC. Conclusions: Our data suggest that susceptibility to pediatric ALL could be associated with at least one copy of an “at risk” haplotype. The study results imply that germ-line genetic variation in TP53, a key tumor suppressor gene, could contribute to susceptibility to childhood ALL. Additional studies are required to confirm our findings and to characterize the functional elements in the at-risk haplotypes.

Pathway analysis of genome-wide association studies (GWAS) offer a unique opportunity to collectively evaluate genetic variants with effects that are too small to be detected individually. We applied a pathway analysis to a bladder cancer GWAS containing data from 3,532 cases and 5,120 controls of European background (n = 5 studies). Thirteen hundred and ninety-nine pathways were drawn from five publicly available resources (Biocarta, Kegg, NCI-PID, HumanCyc, and Reactome), and we constructed 22 additional candidate pathways previously hypothesized to be related to bladder cancer. In total, 1421 pathways, 5647 genes and 90,000 SNPs were included in our study. Logistic regression model adjusting for age, sex, study, DNA source, and smoking status was used to assess the marginal trend effect of SNPs on bladder cancer risk. Two complementary pathway-based methods (gene-set enrichment analysis [GSEA], and adapted rank-truncated product [ARTP]) were used to assess the enrichment of association signals within each pathway. Eighteen pathways were detected by either GSEA or ARTP at P≤0.01. To minimize false positives, we used the I2 statistic to identify SNPs displaying heterogeneous effects across the five studies. After removing these SNPs, seven pathways (‘Aromatic amine metabolism’ [PGSEA = 0.0100, PARTP = 0.0020], ‘NAD biosynthesis’ [PGSEA = 0.0018, PARTP = 0.0086], ‘NAD salvage’ [PARTP = 0.0068], ‘Clathrin derived vesicle budding’ [PARTP = 0.0018], ‘Lysosome vesicle biogenesis’ [PGSEA = 0.0023, PARTP