Search

Top genotyped SNP associations with metastasis at diagnosis in the discovery and all replication studies.

Top 30 loci associated with metastasis at diagnosis in the discovery set of 541 European osteosarcoma cases.

Relationship of variables to metastasis at diagnosis in the discovery set.

SNPs in the NFIB locus associated with metastasis at diagnosis with P{less than or equal to}1x10-5 in the discovery stage, and functionally annotated with information from the ENCODE and 1000 Genomes Project data.

Associations of the top GWAS and imputed SNP with metastasis at diagnosis in the discovery stage (European ancestry) by inheritance model.

Supplementary Figure 1. Manhattan plot of the Discovery stage. SNPs in the NFIB risk locus are highlighted in green. Supplementary Figure 2. LD structure of the NFIB risk locus. Plot was made using the CEU population of the 1000 Genomes Project data with LDlink (ref 54). Supplementary Figure 3. Significant relationship between the genotype of the metastasis associated SNP rs7034162, and expression of NFIB and IGFBP5 in osteosarcoma cell lines and tumors. Significance is based on linear regression comparing the distribution of expression of NFIBbetween the TT homozygous non-risk genotype (combined N=34) of rs7034162 genotypes to the heterozygous risk AT (combined N=10) and the homozygous risk AA genotypes (combined N=2). The green triangle represents the expression levels in the OSA cell line, the blue circle is expression levels in the HOS cell line, and the red square the U2OS cell line expression. These analyses were based on the publically available genotype and expression data from 17 osteosarcoma cell lines and 29 osteosarcoma tumors (ref 20). Supplementary Figure 4. No significant relationship between the genotype of the metastasis associated SNP, rs7034162, and expression of protein-encoding genes FREM1, ZDHHC1, and MPDZ in the neighborhood of NFIB in osteosarcoma cell lines and tumors. Significance is based on linear regression comparing the distribution of expression of FREM1, ZDHHC1, and MPDZ between the TT homozygous non-risk genotype (combined N=34) of rs7034162 genotypes to the heterozygous risk AT (combined N=10) and the homozygous risk AA genotypes (combined N=2). The green triangle represents the expression levels in the OSA cell line, the blue circle is expression levels in the HOS cell line, and the red square the U2OS cell line expression. These analyses were based on the publically available genotype and expression data from 17 osteosarcoma cell lines and 29 osteosarcoma tumors (ref 20). Supplementary Figure 5. The human osteosarcoma cell lines, U2OS and HOS, had higher NFIB protein levels then OSA cells. NFIB levels were determined using immunoblot analysis in the OSA, HOS and U2OS cells both basal (A) and after treatment with siRNA against NFIB (B). Supplementary Figure 6. NFIB expression significantly correlates with IGFBP5 expression in osteosarcoma cell lines and tumors. Significance is based on linear regression comparing the expression levels of NFIB with the corresponding IGFBP5 expression levels. Supplementary Figure 7. Expression levels of IGFBP5 decreases after NFIB suppression of human osteosarcoma cell lines. IGFBP5 expression of OSA, HOS and U2OS cells was determined 48 hours after transfection with control siRNA (si-NEG) or siRNA targeting NFIB (si-NFIB). Graph showing relative expression compared to control treated U2OS cells. Supplementary Figure 8. A model of how the NFIB risk locus leads to changes in NFIB expression and potentially may affect osteosarcoma metastatic potential. The risk allele (right panel) leads to lowered expression of NFIB and NFIB-mediated lower expression of IGFPB5, which then leads to less IGFBP5-mediated inhibition of IGF-1 (right panel, grey arrows), leading to an increase in proliferation, survival and metastasis of osteosarcoma cells (right panel, black arrows). The reference allele (left panel) leads to both higher expression of NFIB and higher expression of NFIB-mediated IGFPB5. In turn, higher IGFBP5 levels leads to increased inhibition of IGF-1 (left panel, black arrows), leading to a decrease in proliferation, survival and metastasis of osteosarcoma cells (left panel, grey arrows). Supplementary Figure 9. Plot of the PCA Eigenvectors of the Discovery stage. Supplementary Figure 10. QQ-plot of the Discovery stage.

Description of participating case studies.

We investigated whether genetic susceptibility to tuberculosis (TB) influences lung adenocarcinoma development among never-smokers using TB genome-wide association study (GWAS) results within the Female Lung Cancer Consortium in Asia. Pathway analysis with the adaptive rank truncated product method was used to assess the association between a TB-related gene-set and lung adenocarcinoma using GWAS data from 5512 lung adenocarcinoma cases and 6277 controls. The gene-set consisted of 31 genes containing known/suggestive associations with genetic variants from previous TB-GWAS. Subsequently, we followed-up with Mendelian Randomization to evaluate the association between TB and lung adenocarcinoma using three genome-wide significant variants from previous TB-GWAS in East Asians. The TB-related gene-set was associated with lung adenocarcinoma (p = 0.016). Additionally, the Mendelian Randomization showed an association between TB and lung adenocarcinoma (OR = 1.31, 95% CI: 1.03, 1.66, p = 0.027). Our findings support TB as a causal risk factor for lung cancer development among never-smoking Asian women.

HPV16 causes half of cervical cancers worldwide; for unknown reasons, most infections resolve within two years. Here, we analyze the viral genomes of 5,328 HPV16-positive case-control samples to investigate mutational signatures and the role of human APOBEC3-induced mutations in viral clearance and cervical carcinogenesis. We identify four de novo mutational signatures, one of which matches the COSMIC APOBEC-associated signature 2. The viral genomes of the precancer/cancer cases are less likely to contain within-host somatic HPV16 APOBEC3-induced mutations (Fisher’s exact test, P = 6.2 x 10−14), and have a 30% lower nonsynonymous APOBEC3 mutation burden compared to controls. We replicate the low prevalence of HPV16 APOBEC3-induced mutations in 1,749 additional cases. APOBEC3 mutations also historically contribute to the evolution of HPV16 lineages. We demonstrate that cervical infections with a greater burden of somatic HPV16 APOBEC3-induced mutations are more likely to be benign or subsequently clear, suggesting they may reduce persistence, and thus progression, within the host., The APOBEC mutational signature is prevalent in different tumour types. Here, using HPV16- positive cervical samples, the authors show that the signature is more prevalent in the viral genome of benign or clearing HPV16 infections compared to the viral genomes of the more advanced precancerous lesions or cervical cancer.

Human papillomavirus (HPV) type 31 (HPV31) is closely related to the most carcinogenic type, HPV16, but only accounts for 4% of cervical cancer cases worldwide. Viral genetic and epigenetic variations have been associated with carcinogenesis for other high-risk HPV types, but little is known about HPV31. We sequenced 2093 HPV31 viral whole genomes from two large studies, one from the U.S. and one international. In addition, we investigated CpG methylation in a subset of 175 samples. We evaluated the association of HPV31 lineages/sublineages, single nucleotide polymorphisms (SNPs) and viral methylation with cervical carcinogenesis. HPV31 A/B clade was &gt, 1.8-fold more associated with cervical intraepithelial neoplasia grade 3 and cancer (CIN3+) compared to the most common C lineage. Lineage/sublineage distribution varied by race/ethnicity and geographic region. A viral genome-wide association analysis identified SNPs within the A/B clade associated with CIN3+, including H23Y (C626T) (odds ratio = 1.60, confidence intervals = 1.17–2.19) located in the pRb CR2 binding-site within the E7 oncogene. Viral CpG methylation was higher in lineage B, compared to the other lineages, and was most elevated in CIN3+. In conclusion, these data support the increased oncogenicity of the A/B lineages and suggest variation of E7 as a contributing risk factor.

Based on our initial genome-wide association study (GWAS) on esophageal squamous cell carcinoma (ESCC) in Han Chinese, we conducted a follow-up study to examine the single nucleotide polymorphisms (SNPs) associated with family history (FH) of upper gastrointestinal cancer (UGI) cancer in cases with ESCC. We evaluated the association between SNPs and FH of UGI cancer among ESCC cases in a stage-1 case-only analysis of the National Cancer Institute (NCI, 541 cases with FH and 1399 without FH) and Henan GWAS (493 cases with FH and 869 without FH) data (discovery phase). The top SNPs (or their surrogates) from discovery were advanced to a stage-2 evaluation in additional Henan subjects (2801 cases with FH and 3136 without FH, replication phase). A total of 19 SNPs were associated with FH of UGI cancer in ESCC cases with P −5 in the stage-1 meta-analysis of NCI and Henan GWAS data. In stage-2, the association for rs79747906 (located at 18p11.31, P = 5.79 × 10−6 in discovery) was replicated (P = 0.006), with a pooled-OR of 1.59 (95%CI: 1.11-2.28). We identified potential genetic variants associated with FH of UGI cancer. Our findings may provide important insights into new low-penetrance susceptibility regions involved in the susceptibility of families with multiple UGI cancer cases.

International audience; Severe aplastic anemia (SAA) is a rare disorder characterized by hypoplastic bone marrow and progressive pancytopenia. The etiology of acquired SAA is not understood but is likely related to abnormal immune responses and environmental exposures. We conducted a genome-wide association study of individuals with SAA genetically matched to healthy controls in discovery (359 cases, 1,396 controls) and validation sets (175 cases, 1,059 controls). Combined analyses identified linked SNPs in distinct blocks within the major histocompatibility complex on 6p21. The top SNP encodes p.Met76Val in the P4 binding pocket of the HLA class II gene HLA-DPB1 (rs1042151A>G, odds ratio [OR] 1.75, 95% confidence interval [CI] 1.50–2.03, p = 1.94 × 10−13) and was associated with HLA-DP cell surface expression in healthy individuals (p = 2.04 × 10−6). Phylogenetic analyses indicate that Val76 is not monophyletic and likely occurs in conjunction with different HLA-DP binding groove conformations. Imputation of HLA-DPB1 alleles revealed increased risk of SAA associated with Val76-encoding alleles DPB1∗03:01, (OR 1.66, p = 1.52 × 10−7), DPB1∗10:01 (OR 2.12, p = 0.0003), and DPB1∗01:01 (OR 1.60, p = 0.0008). A second SNP near HLA-B, rs28367832G>A, reached genome-wide significance (OR 1.49, 95% CI 1.22–1.78, p = 7.27 × 10−9) in combined analyses; the association remained significant after excluding cases with clonal copy-neutral loss-of-heterozygosity affecting class I HLA genes (8.6% of cases and 0% of controls). SNPs in the HLA class II gene HLA-DPB1 and possibly class I (HLA-B) are associated with SAA. The replacement of Met76 to Val76 in certain HLA-DPB1 alleles might influence risk of SAA through mechanisms involving DP peptide binding specificity, expression, and/or other factors affecting DP function.

PLOS ONE 15(9), e0237792 (2020). doi:10.1371/journal.pone.0237792, Published by PLOS, San Francisco, California, US

HPV35 has been found in only ∼2% of invasive cervical cancers (ICC) worldwide but up to 10% in Sub-Saharan Africa, warranting further investigation and consideration of impact on preventive strategies. We studied HPV35 and ethnicity, in relation to the known steps in cervical carcinogenesis, using multiple large epidemiologic studies in the U.S. and internationally. Combining five U.S. studies, we measured HPV35 positivity and, in Northern California, observed HPV35 type-specific population prevalence and estimated 5-year risk of developing precancer when HPV35-positive. HPV35 genetic variation was examined for differences in carcinogenicity in 1053 HPV35+ cervical specimens from a U.S. cohort and an international collection. African-American women had more HPV35 (12.1% vs 5.1%, P < .001) and more HPV35-associated precancers (7.4% vs 2.1%, P < .001) compared to other ethnicities. Precancer risks after HPV35 infection did not vary by ethnicity (global P = .52). The HPV35 A2 sublineage showed an increased association with precancer/cancer in African-Americans (OR = 5.6 vs A1, 95% CI = 1.3-24.8) and A2 was more prevalent among ICC in Africa than other world regions (41.9% vs 10.4%, P < .01). Our analyses support a strong link between HPV35 and cervical carcinogenesis in women of African ancestry. Current HPV vaccines cover the majority of cervical precancer/cancer across all ethnic groups; additional analyses are required to determine whether the addition of HPV35 to the already highly effective nine-valent HPV vaccine would provide better protection for women in Africa or of African ancestry.

Renal cell carcinoma (RCC) has an undisputed genetic component and a stable 2:1 male to female sex ratio in its incidence across populations, suggesting possible sexual dimorphism in its genetic susceptibility. We conducted the first sex-specific genome-wide association analysis of RCC for men (3227 cases, 4916 controls) and women (1992 cases, 3095 controls) of European ancestry from two RCC genome-wide scans and replicated the top findings using an additional series of men (2261 cases, 5852 controls) and women (1399 cases, 1575 controls) from two independent cohorts of European origin. Our study confirmed sex-specific associations for two known RCC risk loci at 14q24.2 (DPF3) and 2p21(EPAS1). We also identified two additional suggestive male-specific loci at 6q24.3 (SAMD5, male odds ratio (OR(male)) = 0.83 [95% CI = 0.78-0.89], P(male) = 1.71 × 10(−8) compared with female odds ratio (OR(female)) = 0.98 [95% CI = 0.90–1.07], P(female) = 0.68) and 12q23.3 (intergenic, OR(male) = 0.75 [95% CI = 0.68-0.83], P(male) = 1.59 × 10(−8) compared with OR(female) = 0.93 [95% CI = 0.82–1.06], P(female) = 0.21) that attained genome-wide significance in the joint meta-analysis. Herein, we provide evidence of sex-specific associations in RCC genetic susceptibility and advocate the necessity of larger genetic and genomic studies to unravel the endogenous causes of sex bias in sexually dimorphic traits and diseases like RCC.

The risk of pancreatic cancer (PC) is increased in melanoma-prone families but the causal relationship between germline CDKN2A mutations and PC risk is uncertain, suggesting the existence of non-CDKN2A factors. One genetic possibility involves patients having mutations in multiple high-risk PC-related genes; however, no systematic examination has yet been conducted. We used next-generation sequencing data to examine 24 putative PC-related genes in 43 PC patients with and 23 PC patients without germline CDKN2A mutations and 1001 controls. For each gene and the four pathways in which they occurred, we tested whether PC patients (overall or CDKN2A+ and CDKN2A- cases separately) had an increased number of rare nonsynonymous variants. Overall, we identified 35 missense variants in PC patients, 14 in CDKN2A+ and 21 in CDKN2A- PC cases. We found nominally significant associations for mismatch repair genes (MLH1, MSH2, MSH6, PMS2) in all PC patients and for ATM, CPA1, and PMS2 in CDKN2A- PC patients. Further, nine CDKN2A+ and four CDKN2A- PC patients had rare potentially deleterious variants in multiple PC-related genes. Loss-of-function variants were only observed in CDKN2A- PC patients, with ATM having the most pathogenic variants. Also, ATM variants (n = 5) were only observed in CDKN2A- PC patients with a family history that included digestive system tumors. Our results suggest that a subset of PC patients may have increased risk because of germline mutations in multiple PC-related genes.

Background Oesophageal cancer is the fourth leading cause of cancer death in China where essentially all cases are histologically oesophageal squamous cell carcinoma (ESCC). Agnostic pathway-based analyses of genome-wide association study (GWAS) data combined with tissue-specific expression quantitative trait loci (eQTL) analysis and publicly available functional data can identify biological pathways and/or genes enriched with functionally-relevant disease-associated variants. Method We used the adaptive multilocus joint test to analyse 1827 pathways containing 6060 genes using GWAS data from 1942 ESCC cases and 2111 controls with Chinese ancestry. We examined the function of risk alleles using in silico and eQTL analyses in oesophageal tissues. Results Associations with ESCC risk were observed for 36 pathways predominantly involved in apoptosis, cell cycle regulation and DNA repair and containing known GWAS-associated genes. After excluding genes with previous GWAS signals, candidate pathways (and genes) for ESCC risk included taste transduction (KEGG_hsa04742; TAS2R13, TAS2R42, TAS2R14, TAS2R46,TAS2R50), long-patch base excision repair (Reactome_pid; POLD2) and the metabolics pathway (KEGG_hsa01100; MTAP, GAPDH, DCTD, POLD2, AMDHD1). We identified and validated CASP8 rs13016963 and IDH2 rs11630814 as eQTLs, and CASP8 rs3769823 and IDH2 rs4561444 as the potential functional variants in high-linkage disequilibrium with these single nucleotide polymorphisms (SNPs), respectively. Further, IDH2 mRNA levels were down-regulated in ESCC (tumour:normal-fold change = 0.69, P = .75E-14). Conclusion Agnostic pathway-based analyses and integration of multiple types of functional data provide new evidence for the contribution of genes in taste transduction and metabolism to ESCC susceptibility, and for the functionality of both established and new ESCC risk-related SNPs.

Metastasis is the leading cause of death in patients with osteosarcoma, the most common pediatric bone malignancy. We conducted a multistage genome-wide association study of osteosarcoma metastasis at diagnosis in 935 osteosarcoma patients to determine whether germline genetic variation contributes to risk of metastasis. We identified an SNP, rs7034162, in NFIB significantly associated with metastasis in European osteosarcoma cases, as well as in cases of African and Brazilian ancestry (meta-analysis of all cases: P = 1.2 × 10−9; OR, 2.43; 95% confidence interval, 1.83–3.24). The risk allele was significantly associated with lowered NFIB expression, which led to increased osteosarcoma cell migration, proliferation, and colony formation. In addition, a transposon screen in mice identified a significant proportion of osteosarcomas harboring inactivating insertions in Nfib and with lowered NFIB expression. These data suggest that germline genetic variation at rs7034162 is important in osteosarcoma metastasis and that NFIB is an osteosarcoma metastasis susceptibility gene. Significance: Metastasis at diagnosis in osteosarcoma is the leading cause of death in these patients. Here we show data that are supportive for the NFIB locus as associated with metastatic potential in osteosarcoma. Cancer Discov; 5(9); 920–31. ©2015 AACR. This article is highlighted in the In This Issue feature, p. 893

Populations in north central China are at high risk for oesophageal squamous cell carcinoma (ESCC) and gastric cancer (GC), and genetic variation in epigenetic machinery genes and pathways may contribute to this risk.We used the adaptive multilocus joint test to analyse 192 epigenetic genes involved in chromatin remodelling, DNA methylation and microRNA biosynthesis in 1942 ESCC and 1758 GC cases [1126 cardia (GCA) and 632 non-cardia adenocarcinoma (GNCA)] and 2111 controls with Chinese ancestry. We examined potential function of risk alleles using in silico and expression quantitative trait loci (eQTLs) analyses.Suggestive pathway-based associations were observed for the overall epigenetic (P-value(PATH) = 0.034) and chromatin remodelling (P-value(PATH) = 0.039) pathways with risk of GCA, but not GC, GNCA or ESCC. Overall, 37 different epigenetic machinery genes were associated with risk of one or more upper gastrointestinal (UGI) cancer sites (P-value(GENE ) 0.05), including 14 chromatin remodelling genes whose products are involved in the regulation of HOX genes. We identified a gastric eQTL (rs12724079; rho = 0.37; P = 0.0006) which regulates mRNA expression of ASH1L. Several suggestive eQTLs were also found in oesophageal (rs10898459 in EED), gastric cardia (rs7157322 in DICER1; rs8179271 in ASH1L), and gastric non-cardia (rs1790733 in PPP1CA) tissues.Results of our analyses provide limited but suggestive evidence for a role of epigenetic gene variation in the aetiology of UGI cancer.

Recurrent large-scale somatic copy number alterations (SCNAs), and somatic point mutations can be analysed to stratify patients with chronic lymphocytic leukaemia (CLL) into distinct prognostic groups. To investigate the relationship between SCNAs and somatic mutations, we performed whole-exome sequencing and single nucleotide polymorphism microarray analyses on 98 CLL patients from 40 families with a high burden of CLL. Overall, 69 somatic mutations in 29 CLL driver genes were detected among 45 subjects (46%), with the most frequently mutated genes being TP53 (8·2%), NOTCH1 (8·2%) and ATM (5·1%). Additionally, 142 SCNAs from 54 subjects (57%) were detected, including losses of chromosome 13q14 (28·9%), 11q (5·6%), 17p (2·1%), and gain of chromosome 12 (4·2%). We found that patients having both an adverse point mutation in a CLL driver gene and an unfavourable SCNA tended to have poorer survival (Hazard ratio [HR] = 3·17, 95% confidence interval [CI] = 0·97-10·35; P = 0·056) than patients having either a point mutation (HR = 1·34, 95%CI = 0·66-2·71; P = 0·42) or SCNAs (HR = 2·65, 95%CI = 0·77-9·13; P = 0·12). TP53 mutation carriers were associated with the poorest overall survival (HR = 4·39, 95%CI = 1·28-15·04; P = 0·018). Our study suggests that combining SCNA and mutational data could contribute to predicting outcome in familial CLL.

Human papillomavirus type 16 (HPV16) is the most carcinogenic HPV, causing >50% of cervical cancers. Its unique oncogenicity is unsolved, e.g., the second most carcinogenic type (HPV18, causing ~16% of cancers) is relatively distantly related to HPV16. Mirabello et al. recently reported on viral genome data from 3,215 HPV positive specimens from women undergoing cervical cancer screening at Kaiser Permanente Northern California in the Persistence and Progression (PaP) cohort. Results demonstrated profound differences in disease risk by histologic subtype among the HPV16 sublineages (e.g., A1), with over 100-fold risk differences. Here, we expanded this analysis to examine the genetic and evolutionary underpinnings of HPV16 carcinogenicity through an exhaustive analysis of viral nucleotide distance. The HPV16 genome displays significant evidence of purifying selection (dN/dS = 0.267; PdN/dS = 0.049) than in controls (dN/dS = 0.27; PdN/dS values up to 59.44. Twelve (12) of these regions overlap phylogenetic sublineage-defining residues, and 9 match known HPV epitope sequences obtained from the Immune Epitope Database and/or the Human Papillomavirus T Cell Antigen Database. In particular, two regions of E6 were identified independently in cases and controls: codons 20-27 and codons 75-90, both of which overlap with experimentally verified HLA epitopes, and the second of which overlaps known antibody epitopes. Moreover, the 75-90 region exhibits a dN/dS ratio of 59.44 in cases and 28.05 in controls, shows substantial case vs. control divergence (between – within group dN/dS of 55.60), and includes a L83V variant that has been reported to contribute to persistence in European populations. We conclude that diversifying positive selection likely played a key role in the historical divergence of HPV16 variants, possibly as the result of inter-host environmental heterogeneity based on host immune genotype. Moreover, since our data imply that positive selection is targeted to many of the same loci that are diagnostic as sublineage-defining residues, it is likely that similar evolutionary pressures have operated throughout the evolutionary diversification of HPV16. In particular, we suggest that host immune genotype (e.g., HLA) may play a key role in disease outcome, and must be prioritized in future studies of HPV evolution and its link to cervical cancer.

Recent evidence from several relatively small nested case-control studies in prospective cohorts shows an association between longer telomere length measured phenotypically in peripheral white blood cell (WBC) DNA and increased lung cancer risk. We sought to further explore this relationship by examining a panel of seven telomere-length associated genetic variants in a large study of 5,457 never-smoking female Asian lung cancer cases and 4,493 never-smoking female Asian controls using data from a previously reported genome-wide association study. Using a group of 1,536 individuals with phenotypically measured telomere length in WBCs in the prospective Shanghai Women's Health study, we demonstrated the utility of a genetic risk score (GRS) of seven telomere-length associated variants to predict telomere length in an Asian population. We then found that GRSs used as instrumental variables to predict longer telomere length were associated with increased lung cancer risk (OR = 1.51 (95% CI = 1.34-1.69) for upper vs. lower quartile of the weighted GRS, p value = 4.54 × 10(-14) ) even after removing rs2736100 (p value = 4.81 × 10(-3) ), a SNP in the TERT locus robustly associated with lung cancer risk in prior association studies. Stratified analyses suggested the effect of the telomere-associated GRS is strongest among younger individuals. We found no difference in GRS effect between adenocarcinoma and squamous cell subtypes. Our results indicate that a genetic background that favors longer telomere length may increase lung cancer risk, which is consistent with earlier prospective studies relating longer telomere length with increased lung cancer risk.

CYBB encodes the gp91-phox protein of the phagocytic NADPH oxidase; the innate immunity-related enzymatic complex responsible for the respiratory burst. Mutations in CYBB can cause chronic granulomatous disease (CGD), a primary immunodeficiency characterized by ineffective microbicidal activity, for which over 150 family-specific mutations have been described. It is also plausible that common SNPs in CYBB alter the expression or function of gp91-phox, determining differences in susceptibility to complex disorders such as autoimmune or infectious diseases. We have resequenced the exons, UTRs, and intronic regions of CYBB in 102 ethnically diverse individuals and genotyped nine tag-SNPs in 942 individuals from 52 worldwide populations. The 28 observed SNPs (none of which nonsynonymous) reside on 28 haplotypes that can be collapsed into five clades. CYBB shows lower diversity than other X-chromosome genes and most of the between-population genetic variance was observed among Africans and non-Africans. The African population shows the highest diversity and the lowest linkage disequilibrium (LD). Because there is extensive shared LD among non,Africans, tag-SNPs can be effectively employed in gene,centric association studies and are portable across Eurasian and Native American populations. Comparison of CYBB coding sequences among mammals evidences the action of long-term purifying selection, which is stronger on the C-terminal cytosolic domain than on the N-terminal transmembrane domain of gp91-phox.

Introduction. Acquired severe aplastic anemia (SAA) is a life-threatening disorder characterized by severe progressive pancytopenia and hypocellular bone marrow. The etiology of acquired SAA is not understood but believed to be related to abnormal immune responses to environmental exposures. We conducted a genome-wide association study (GWAS) to identify common germline variants associated with SAA. Methods. We identified 895 patients with SAA who underwent related or unrelated hematopoietic cell transplant (HCT) with clinical data and pre-HCT blood samples available in the Center for International Blood and Marrow Transplant Research (CIBMTR) database and biorepository. Pre-HCT DNA was extracted from blood of patients with SAA and genome-wide genotyping was conducted using the Illumina OmniExpress array. We excluded 93 inherited bone marrow failure cases. The SAA cases were grouped cases into discovery and validation sets based on time of batch sample receipt. Analyses were limited to patients of European ancestry based on principal component analyses to minimize the potential effect of population stratification. Controls were genomically matched and selected from previously scanned cancer-free subjects at the Cancer Genomics Research Laboratory, Division of Cancer Epidemiology and Genetics, NCI. The final analysis included 534 acquired SAA cases (359 in the discovery set and 175 in the validation set), and 2,455 controls (1,396 in the discovery set, and 1,059 in the validation set). Results. Patients with SAA in this study received HCT between 1989-2015 at a median age of 21 years, 56% were male, and the median time between SAA diagnosis and HCT was 11 months. Strong genome-wide association signals were identified across the human leukocyte antigen (HLA) genes encoded at the major histocompatibility complex (MHC) on chromosome 6p21 (Figure 1). The top SNP was located in the P4 binding pocket of the HLA class II gene HLA-DPB1(rs1042151A>G, p.Met105Val, pooled-odds ratio [OR]=1.75, 95% confidence interval [CI]=1.50-2.03, p=1.94x10-13). The expression of HLA-DP in CD19+ cells from 175 healthy donors was significantly different by rs1042151 A>G genotype (p=2.04x10-6) (Figure 2). A second SNP near HLA-B, rs28367832G>A, also reached genome-wide significance (pooled-OR=1.49, 95% CI=1.22-1.78, p=7.27x10-9). Copy-number variant analysis and next generation sequencing also identified somatic, clonal copy-neutral loss-of-heterozygosity affecting class I HLA genes in 8.6% of the SAA cases and none of the controls. Conclusion. This SAA GWAS identified strong association signals between common germline genetic variants in HLA class I and II genes and SAA. The main SNP is associated with changes in HLA-DP expression suggesting a key role for this locus in SAA etiology. This study adds further evidence to the connection between SAA and immune dysregulation. Disclosures Cerhan: Janssen: Membership on an entity's Board of Directors or advisory committees; NanoString: Research Funding; Celgene: Research Funding. Lee:Incyte: Research Funding; Syndax: Research Funding; Amgen: Research Funding; Novartis: Research Funding; Takeda: Research Funding; Kadmon: Consultancy, Research Funding; Pfizer: Consultancy, Research Funding; AstraZeneca: Research Funding.

Populations in north central China are at high risk for gastric cancers (GC), and altered FAS-mediated cell signaling and/or apoptosis may contribute to this risk. We examined the association of 554 single nucleotide polymorphisms (SNPs) in 53 Fas signaling-related genes using a pathway-based approach in 1758 GC cases (1126 gastric cardia adenocarcinomas (GCA) and 632 gastric noncardia adenocarcinomas (GNCA)), and 2111 controls from a genome-wide association study (GWAS) of GC in ethnic Chinese. SNP associations with risk of overall GC, GCA and GNCA were evaluated using unconditional logistic regressions controlling for age, sex and study. Gene- and pathway-based associations were tested using the adaptive rank-truncated product (ARTP) method. Statistical significance was evaluated empirically by permutation. Significant pathway-based associations were observed for Fas signaling with risk of overall GC (p = 5.5E-04) and GCA (p = 6.3E-03), but not GNCA (p= 8.1E-02). Among examined genes in the Fas signaling pathway, MAP2K4, FAF1, MAPK8, CASP10, CASP8, CFLAR, MAP2K1, CAP8AP2, PAK2 and IKBKB were associated with risk of GC (nominal p < 0.05), and FAF1 and MAPK8 were significantly associated with risk of both GCA and GNCA (nominal p< 0.05). Our examination of genetic variation in the Fas signaling pathway is consistent with an association of altered Fas signaling and/or apoptosis with risk of GC. As one of the first attempts to investigate a pathway-level association, our results suggest that these genes and the Fas signaling pathway warrant further evaluation in relation to GC risk in other populations.

The DNA repair pathways help to maintain genomic integrity and therefore genetic variation in the pathways could affect the propensity to develop cancer. Selected germline single nucleotide polymorphisms (SNPs) in the pathways have been associated with esophageal cancer and gastric cancer (GC) but few studies have comprehensively examined the pathway genes. We aimed to investigate associations between DNA repair pathway genes and risk of esophageal squamous cell carcinoma (ESCC) and GC, using data from a genome-wide association study in a Han Chinese population where ESCC and GC are the predominant cancers. In sum, 1942 ESCC cases, 1758 GC cases and 2111 controls from the Shanxi Upper Gastrointestinal Cancer Genetics Project (discovery set) and the Linxian Nutrition Intervention Trials (replication set) were genotyped for 1675 SNPs in 170 DNA repair-related genes. Logistic regression models were applied to evaluate SNP-level associations. Gene- and pathway-level associations were determined using the resampling-based adaptive rank-truncated product approach. The DNA repair pathways overall were significantly associated with risk of ESCC (P = 6.37 × 10(-4)), but not with GC (P = 0.20). The most significant gene in ESCC was CHEK2 (P = 2.00 × 10(-6)) and in GC was CLK2 (P = 3.02 × 10(-4)). We observed several other genes significantly associated with either ESCC (SMUG1, TDG, TP53, GTF2H3, FEN1, POLQ, HEL308, RAD54B, MPG, FANCE and BRCA1) or GC risk (MRE11A, RAD54L and POLE) (P < 0.05). We provide evidence for an association between specific genes in the DNA repair pathways and the risk of ESCC and GC. Further studies are warranted to validate these associations and to investigate underlying mechanisms.

To the editor: Chronic lymphocytic leukemia (CLL) is one of the most common B-cell malignancies in individuals of European ancestry in the United States.[1][1] Although advanced age, white ancestry, and family history of hematologic malignancies are risk factors, the etiology of CLL is unknown.[2][

Genome-wide association studies (GWASs) of follicular lymphoma (FL) have previously identified human leukocyte antigen (HLA) gene variants. To identify additional FL susceptibility loci, we conducted a large-scale two-stage GWAS in 4,523 case subjects and 13,344 control subjects of European ancestry. Five non-HLA loci were associated with FL risk: 11q23.3 (rs4938573, p = 5.79 × 10 -20) near CXCR5; 11q24.3 (rs4937362, p = 6.76 × 10 -11) near ETS1; 3q28 (rs6444305, p = 1.10 × 10 -10) in LPP; 18q21.33 (rs17749561, p = 8.28 × 10 -10) near BCL2; and 8q24.21 (rs13254990, p = 1.06 × 10 -8) near PVT1. In an analysis of the HLA region, we identified four linked HLA-DRß1 multiallelic amino acids at positions 11, 13, 28, and 30 that were associated with FL risk (pomnibus = 4.20 × 10 -67 to 2.67 × 10 -70). Additional independent signals included rs17203612 in HLA class II (odds ratio [ORper-allele] = 1.44; p = 4.59 × 10 -16) and rs3130437 in HLA class I (ORper-allele = 1.23; p = 8.23 × 10 -9). Our findings further expand the number of loci associated with FL and provide evidence that multiple common variants outside the HLA region make a significant contribution to FL risk. © 2014 by The American Society of Human Genetics. All rights reserved.

Experimental and epidemiologic investigations suggest that certain pesticides may alter sex steroid hormone synthesis, metabolism or regulation, and the risk of hormone-related cancers. Here, we evaluated whether single-nucleotide polymorphisms (SNPs) involved in hormone homeostasis alter the effect of pesticide exposure on prostate cancer risk. We evaluated pesticide–SNP interactions between 39 pesticides and SNPs with respect to prostate cancer among 776 cases and 1,444 controls nested in the Agricultural Health Study cohort. In these interactions, we included candidate SNPs involved in hormone synthesis, metabolism or regulation (N = 1,100), as well as SNPs associated with circulating sex steroid concentrations, as identified by genome-wide association studies (N = 17). Unconditional logistic regression was used to estimate odds ratios (ORs) and 95% confidence intervals (CIs). Multiplicative SNP–pesticide interactions were calculated using a likelihood ratio test. We translated p-values for interaction into q-values, which reflected the false discovery rate, to account for multiple comparisons. We observed a significant interaction, which was robust to multiple comparison testing, between the herbicide dicamba and rs8192166 in the testosterone metabolizing gene SRD5A1 (p-interaction = 4.0 × 10−5; q-value = 0.03), such that men with two copies of the wild-type genotype CC had a reduced risk of prostate cancer associated with low use of dicamba (OR = 0.62 95% CI: 0.41, 0.93) and high use of dicamba (OR = 0.44, 95% CI: 0.29, 0.68), compared to those who reported no use of dicamba; in contrast, there was no significant association between dicamba and prostate cancer among those carrying one or two copies of the variant T allele at rs8192166. In addition, interactions between two organophosphate insecticides and SNPs related to estradiol metabolism were observed to result in an increased risk of prostate cancer. While replication is needed, these data suggest both agonistic and antagonistic effects on circulating hormones, due to the combination of exposure to pesticides and genetic susceptibility, may impact prostate cancer risk.

International audience; Evidence from a small number of studies suggests that longer telomere length measured in peripheral leukocytes is associated with an increased risk of non-Hodgkin lymphoma (NHL). However, these studies may be biased by reverse causation, confounded by unmeasured environmental exposures and might miss time points for which prospective telomere measurement would best reveal a relationship between telomere length and NHL risk. We performed an analysis of genetically inferred telomere length and NHL risk in a study of 10 102 NHL cases of the four most common B-cell histologic types and 9562 controls using a genetic risk score (GRS) comprising nine telomere length-associated single-nucleotide polymorphisms. This approach uses existing genotype data and estimates telomere length by weighing the number of telomere length-associated variant alleles an individual carries with the published change in kb of telomere length. The analysis of the telomere length GRS resulted in an association between longer telomere length and increased NHL risk [four B-cell histologic types combined; odds ratio (OR) = 1.49, 95% CI 1.22-1.82,P-value = 8.5 x 10(-5)]. Subtype-specific analyses indicated that chronic lymphocytic leukemia or small lymphocytic lymphoma (CLL/SLL) was the principal NHL subtype contributing to this association (OR = 2.60, 95% CI 1.93-3.51,P-value = 4.0 x 10(-10)). Significant interactions were observed across strata of sex for CLL/SLL and marginal zone lymphoma subtypes as well as age for the follicular lymphoma subtype. Our results indicate that a genetic background that favors longer telomere length may increase NHL risk, particularly risk of CLL/SLL, and are consistent with earlier studies relating longer telomere length with increased NHL risk

Genome-wide association studies have identified susceptibility loci for esophageal squamous cell carcinoma (ESCC). We conducted a meta-analysis of all single-nucleotide polymorphisms (SNPs) that showed nominally significant P-values in two previously published genome-wide scans that included a total of 2961 ESCC cases and 3400 controls. The meta-analysis revealed five SNPs at 2q33 with P5 × 10(-8), and the strongest signal was rs13016963, with a combined odds ratio (95% confidence interval) of 1.29 (1.19-1.40) and P= 7.63 × 10(-10). An imputation analysis of 4304 SNPs at 2q33 suggested a single association signal, and the strongest imputed SNP associations were similar to those from the genotyped SNPs. We conducted an ancestral recombination graph analysis with 53 SNPs to identify one or more haplotypes that harbor the variants directly responsible for the detected association signal. This showed that the five SNPs exist in a single haplotype along with 45 imputed SNPs in strong linkage disequilibrium, and the strongest candidate was rs10201587, one of the genotyped SNPs. Our meta-analysis found genome-wide significant SNPs at 2q33 that map to the CASP8/ALS2CR12/TRAK2 gene region. Variants in CASP8 have been extensively studied across a spectrum of cancers with mixed results. The locus we identified appears to be distinct from the widely studied rs3834129 and rs1045485 SNPs in CASP8. Future studies of esophageal and other cancers should focus on comprehensive sequencing of this 2q33 locus and functional analysis of rs13016963 and rs10201587 and other strongly correlated variants.

Genetic variation on the Y chromosome has not been convincingly implicated in prostate cancer risk. To comprehensively analyze the role of inherited Y chromosome variation in prostate cancer risk in individuals of European ancestry, we genotyped 34 binary Y chromosome markers in 3,995 prostate cancer cases and 3,815 control subjects drawn from four studies. In this set, we identified nominally significant association between a rare haplogroup, E1b1b1c, and prostate cancer in stage I (P = 0.012, OR = 0.51; 95% confidence interval 0.30–0.87). Population substructure of E1b1b1c carriers suggested Ashkenazi Jewish ancestry, prompting a replication phase in individuals of both European and Ashkenazi Jewish ancestry. The association was not significant for prostate cancer overall in studies of either Ashkenazi Jewish (1,686 cases and 1,597 control subjects) or European (686 cases and 734 control subjects) ancestry (Pmeta = 0.078), but a meta-analysis of stage I and II studies revealed a nominally significant association with prostate cancer risk (Pmeta = 0.010, OR = 0.77; 95% confidence interval 0.62–0.94). Comparing haplogroup frequencies between studies, we noted strong similarities between those conducted in the US and France, in which the majority of men carried R1 haplogroups, resembling Northwestern European populations. On the other hand, Finns had a remarkably different haplogroup distribution with a preponderance of N1c and I1 haplogroups. In summary, our results suggest that inherited Y chromosome variation plays a limited role in prostate cancer etiology in European populations but warrant follow-up in additional large and well characterized studies of multiple ethnic backgrounds. Electronic supplementary material The online version of this article (doi:10.1007/s00439-012-1139-5) contains supplementary material, which is available to authorized users.

In follow-up of a recent genome-wide association study (GWAS) that identified a locus in chromosome 2p21 associated with risk for renal cell carcinoma (RCC), we conducted a fine mapping analysis of a 120 kb region that includes EPAS1. We genotyped 59 tagged common single-nucleotide polymorphisms (SNPs) in 2278 RCC and 3719 controls of European background and observed a novel signal for rs9679290 [P=5 5.75× 3 10 -8, per-allele odds ratio (OR) = 1.27, 95% confidence interval (CI): 1.17-1.39]. Imputation of common SNPs surrounding rs9679290 using HapMap 3 and 1000 Genomes data yielded two additional signals, rs4953346 (P = 4.09 × 10 -14) and rs12617313 (P = 7.48 × 10 -12), both highly correlated withrs9679290 (r 2 > 0.95), but interestingly not correlated with the two SNPs reported in the GWAS: rs11894252 and rs7579899 (r 2 < 0.1 with rs9679290). Genotype analysis of rs12617313 confirmed an association with RCC risk (P = 1.72 × 10 -9, per-allele OR = 1.28, 95% CI: 1.18-1.39) In conclusion, we report that chromosome 2p21 harbors a complex genetic architecture for common RCC risk variants. Published by Oxford University Press 2011.