Inferred expression regulator activities suggest genes mediating cardiometabolic genetic signals.

Expression QTL (eQTL) analyses have suggested many genes mediating genome-wide association study (GWAS) signals but most GWAS signals still lack compelling explanatory genes. We have leveraged an adipose-specific gene regulatory network to infer expression regulator activities and phenotypic master regulators (MRs), which were used to detect activity QTLs (aQTLs) at cardiometabolic trait GWAS loci. Regulator activities were inferred with the VIPER algorithm that integrates enrichment of expected expression changes among a regulator's target genes with confidence in their regulator-target network interactions and target overlap between different regulators (i.e., pleiotropy). Phenotypic MRs were identified as those regulators whose activities were most important in predicting their respective phenotypes using random forest modeling. While eQTLs were typically more significant than aQTLs in cis, the opposite was true among candidate MRs in trans. Several GWAS loci colocalized with MR trans-eQTLs/aQTLs in the absence of colocalized cis-QTLs. Intriguingly, at the 1p36.1 BMI GWAS locus the EPHB2 cis-aQTL was stronger than its cis-eQTL and colocalized with the GWAS signal and 35 BMI MR trans-aQTLs, suggesting the GWAS signal may be mediated by effects on EPHB2 activity and its downstream effects on a network of BMI MRs. These MR and aQTL analyses represent systems genetic methods that may be broadly applied to supplement standard eQTL analyses for suggesting molecular effects mediating GWAS signals.

Germline variants in hereditary breast cancer genes are associated with early age at diagnosis and family history in Guatemalan breast cancer.

PURPOSE: Mutations in hereditary breast cancer genes play an important role in the risk for cancer. METHODS: Cancer susceptibility genes were sequenced in 664 unselected breast cancer cases from Guatemala. Variants were annotated with ClinVar and VarSome. RESULTS: A total of 73 out of 664 subjects (11%) had a pathogenic variant in a high or moderate penetrance gene. The most frequently mutated genes were BRCA1 (37/664, 5.6%) followed by BRCA2 (15/664, 2.3%), PALB2 (5/664, 0.8%), and TP53 (5/664, 0.8%). Pathogenic variants were also detected in the moderate penetrance genes ATM, BARD1, CHEK2, and MSH6. The high ratio of BRCA1/BRCA2 mutations is due to two potential founder mutations: BRCA1 c.212 + 1G > A splice mutation (15 cases) and BRCA1 c.799delT (9 cases). Cases with pathogenic mutations had a significantly earlier age at diagnosis (45 vs 51 years, P < 0.001), are more likely to have had diagnosis before menopause, and a higher percentage had a relative with any cancer (51% vs 37%, P = 0.038) or breast cancer (33% vs 15%, P < 0.001). CONCLUSIONS: Hereditary breast cancer mutations were observed among Guatemalan women, and these women are more likely to have early age at diagnosis and family history of cancer. These data suggest the use of genetic testing in breast cancer patients and those at high risk as part of a strategy to reduce breast cancer mortality in Guatemala.

Characterising cis-regulatory variation in the transcriptome of histologically normal and tumour-derived pancreatic tissues.

OBJECTIVE: To elucidate the genetic architecture of gene expression in pancreatic tissues. DESIGN: We performed expression quantitative trait locus (eQTL) analysis in histologically normal pancreatic tissue samples (n=95) using RNA sequencing and the corresponding 1000 genomes imputed germline genotypes. Data from pancreatic tumour-derived tissue samples (n=115) from The Cancer Genome Atlas were included for comparison. RESULTS: We identified 38 615 cis-eQTLs (in 484 genes) in histologically normal tissues and 39 713 cis-eQTL (in 237 genes) in tumour-derived tissues (false discovery rate <0.1), with the strongest effects seen near transcriptional start sites. Approximately 23% and 42% of genes with significant cis-eQTLs appeared to be specific for tumour-derived and normal-derived tissues, respectively. Significant enrichment of cis-eQTL variants was noted in non-coding regulatory regions, in particular for pancreatic tissues (1.53-fold to 3.12-fold, p≤0.0001), indicating tissue-specific functional relevance. A common pancreatic cancer risk locus on 9q34.2 (rs687289) was associated with ABO expression in histologically normal (p=5.8×10-8) and tumour-derived (p=8.3×10-5) tissues. The high linkage disequilibrium between this variant and the O blood group generating deletion variant in ABO (exon 6) suggested that nonsense-mediated decay (NMD) of the 'O' mRNA might explain this finding. However, knockdown of crucial NMD regulators did not influence decay of the ABO 'O' mRNA, indicating that a gene regulatory element influenced by pancreatic cancer risk alleles may underlie the eQTL. CONCLUSIONS: We have identified cis-eQTLs representing potential functional regulatory variants in the pancreas and generated a rich data set for further studies on gene expression and its regulation in pancreatic tissues.

Identification of a novel susceptibility locus at 13q34 and refinement of the 20p12.2 region as a multi-signal locus associated with bladder cancer risk in individuals of European ancestry.

Candidate gene and genome-wide association studies (GWAS) have identified 15 independent genomic regions associated with bladder cancer risk. In search for additional susceptibility variants, we followed up on four promising single-nucleotide polymorphisms (SNPs) that had not achieved genome-wide significance in 6911 cases and 11 814 controls (rs6104690, rs4510656, rs5003154 and rs4907479, P < 1 × 10(-6)), using additional data from existing GWAS datasets and targeted genotyping for studies that did not have GWAS data. In a combined analysis, which included data on up to 15 058 cases and 286 270 controls, two SNPs achieved genome-wide statistical significance: rs6104690 in a gene desert at 20p12.2 (P = 2.19 × 10(-11)) and rs4907479 within the MCF2L gene at 13q34 (P = 3.3 × 10(-10)). Imputation and fine-mapping analyses were performed in these two regions for a subset of 5551 bladder cancer cases and 10 242 controls. Analyses at the 13q34 region suggest a single signal marked by rs4907479. In contrast, we detected two signals in the 20p12.2 region-the first signal is marked by rs6104690, and the second signal is marked by two moderately correlated SNPs (r(2) = 0.53), rs6108803 and the previously reported rs62185668. The second 20p12.2 signal is more strongly associated with the risk of muscle-invasive (T2-T4 stage) compared with non-muscle-invasive (Ta, T1 stage) bladder cancer (case-case P ≤ 0.02 for both rs62185668 and rs6108803). Functional analyses are needed to explore the biological mechanisms underlying these novel genetic associations with risk for bladder cancer.

Meta-analysis of genome-wide association studies identifies multiple lung cancer susceptibility loci in never-smoking Asian women.

Genome-wide association studies (GWAS) of lung cancer in Asian never-smoking women have previously identified six susceptibility loci associated with lung cancer risk. To further discover new susceptibility loci, we imputed data from four GWAS of Asian non-smoking female lung cancer (6877 cases and 6277 controls) using the 1000 Genomes Project (Phase 1 Release 3) data as the reference and genotyped additional samples (5878 cases and 7046 controls) for possible replication. In our meta-analysis, three new loci achieved genome-wide significance, marked by single nucleotide polymorphism (SNP) rs7741164 at 6p21.1 (per-allele odds ratio (OR) = 1.17; P = 5.8 × 10(-13)), rs72658409 at 9p21.3 (per-allele OR = 0.77; P = 1.41 × 10(-10)) and rs11610143 at 12q13.13 (per-allele OR = 0.89; P = 4.96 × 10(-9)). These findings identified new genetic susceptibility alleles for lung cancer in never-smoking women in Asia and merit follow-up to understand their biological underpinnings.

Identification of new susceptibility loci for gastric non-cardia adenocarcinoma: pooled results from two Chinese genome-wide association studies.

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.

Common genetic variation and risk of gallbladder cancer in India: a case-control genome-wide association study.

BACKGROUND: Gallbladder cancer is highly lethal, with notable differences in incidence by geography and ethnic background. The aim of this study was to identify common genetic susceptibility alleles for gallbladder cancer. METHODS: In this case-control genome-wide association study (GWAS), we did a genome-wide scan of gallbladder cancer cases and hospital visitor controls, both of Indian descent, followed by imputation across the genome. Cases were patients aged 20-80 years with microscopically confirmed primary gallbladder cancer diagnosed or treated at Tata Memorial Hospital, Mumbai, India, and enrolled in the study between Sept 12, 2010, and June 8, 2015. We only included patients who had been diagnosed less than 1 year before the date of enrolment and excluded patients with any other malignancies. We recruited visitor controls aged 20-80 years with no history of cancer visiting all departments or units of Tata Memorial Hospital during the same time period and frequency matched them to cases on the basis of age, sex, and current region of residence. We estimated association using logistic regression, adjusting for age, sex, and five eigenvectors. We recruited samples for a replication cohort from patients visiting Tata Memorial Hospital between Aug 4, 2015, and May 17, 2016, and patients visiting the Sanjay Gandhi Postgraduate Institute of Medical Sciences, Lucknow, India, between July, 2010, and May, 2015. We used the same inclusion and exclusion criteria for the replication set. We examined three of the most significant single-nucleotide polymorphisms (SNPs) in the replication cohort and did a meta-analysis of the GWAS discovery and replication sets to get combined estimates of association. FINDINGS: The discovery cohort comprised 1042 gallbladder cancer cases and 1709 controls and the replication cohort contained 428 gallbladder cancer cases and 420 controls. We observed genome-wide significant associations for several markers in the chromosomal region 7q21.12 harbouring both the ABCB1 and ABCB4 genes, with the most notable SNPs after replication and meta-analysis being rs1558375 (GWAS p=3·8 × 10-9; replication p=0·01; combined p=2·3 × 10-10); rs17209837 (GWAS p=2·0 × 10-8; replication p=0·02; combined p=2·3 × 10-9), and rs4148808 (GWAS p=2·4 × 10-8; replication p=0·008; combined p=2·7 × 10-9). Combined estimates of per-allele trend odds ratios were 1·47 (95% CI 1·30-1·66; p=2·31 × 10-10) for rs1558375, 1·61 (1·38-1·89; p=2·26 × 10-9) for rs17209837, and 1·57 (1·35-1·82; p=2·71 × 10-9) for rs4148808. GWAS heritability analysis suggested that common variants are associated with substantial variation in risk of gallbladder cancer (sibling relative risk 3·15 [95% CI 1·80-5·49]). INTERPRETATION: To our knowledge, this study is the first report of common genetic variation conferring gallbladder cancer risk at genome-wide significance. This finding, along with in-silico and biological evidence indicating the potential functional significance of ABCB1 and ABCB4, underlines the likely importance of these hepatobiliary phospholipid transporter genes in the pathology of gallbladder cancer. FUNDING: The Tata Memorial Centre and Department of Biotechnology.

Integration of multiethnic fine-mapping and genomic annotation to prioritize candidate functional SNPs at prostate cancer susceptibility regions.

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 genome-encoded 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 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.

Genome-wide association study identifies five susceptibility loci for follicular lymphoma outside the HLA region.

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 [OR(per-allele)] = 1.44; p = 4.59 × 10(-16)) and rs3130437 in HLA class I (OR(per-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.

Pathway-based analysis of GWAs data identifies association of sex determination genes with susceptibility to testicular germ cell tumors.

Genome-wide association (GWA) studies of testicular germ cell tumor (TGCT) have identified 18 susceptibility loci, some containing genes encoding proteins important in male germ cell development. Deletions of one of these genes, DMRT1, lead to male-to-female sex reversal and are associated with development of gonadoblastoma. To further explore genetic association with TGCT, we undertook a pathway-based analysis of SNP marker associations in the Penn GWAs (349 TGCT cases and 919 controls). We analyzed a custom-built sex determination gene set consisting of 32 genes using three different methods of pathway-based analysis. The sex determination gene set ranked highly compared with canonical gene sets, and it was associated with TGCT (FDRG = 2.28 × 10(-5), FDRM = 0.014 and FDRI = 0.008 for Gene Set Analysis-SNP (GSA-SNP), Meta-Analysis Gene Set Enrichment of Variant Associations (MAGENTA) and Improved Gene Set Enrichment Analysis for Genome-wide Association Study (i-GSEA4GWAS) analysis, respectively). The association remained after removal of DMRT1 from the gene set (FDRG = 0.0002, FDRM = 0.055 and FDRI = 0.009). Using data from the NCI GWA scan (582 TGCT cases and 1056 controls) and UK scan (986 TGCT cases and 4946 controls), we replicated these findings (NCI: FDRG = 0.006, FDRM = 0.014, FDRI = 0.033, and UK: FDRG = 1.04 × 10(-6), FDRM = 0.016, FDRI = 0.025). After removal of DMRT1 from the gene set, the sex determination gene set remains associated with TGCT in the NCI (FDRG = 0.039, FDRM = 0.050 and FDRI = 0.055) and UK scans (FDRG = 3.00 × 10(-5), FDRM = 0.056 and FDRI = 0.044). With the exception of DMRT1, genes in the sex determination gene set have not previously been identified as TGCT susceptibility loci in these GWA scans, demonstrating the complementary nature of a pathway-based approach for genome-wide analysis of TGCT.

Genome-wide association study identifies multiple loci associated with bladder cancer risk.

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.

Imputation and subset-based association analysis across different cancer types identifies multiple independent risk loci in the TERT-CLPTM1L region on chromosome 5p15.33.

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 × 10(-39); Region 3: rs2853677, P = 3.30 × 10(-36) and PConditional = 2.36 × 10(-8); Region 4: rs2736098, P = 3.87 × 10(-12) and PConditional = 5.19 × 10(-6), Region 5: rs13172201, P = 0.041 and PConditional = 2.04 × 10(-6); and Region 6: rs10069690, P = 7.49 × 10(-15) and PConditional = 5.35 × 10(-7)) and one in the neighboring CLPTM1L gene (Region 2: rs451360; P = 1.90 × 10(-18) and PConditional = 7.06 × 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.

Further Confirmation of Germline Glioma Risk Variant rs78378222 in TP53 and Its Implication in Tumor Tissues via Integrative Analysis of TCGA Data.

We confirmed strong association of rs78378222:A>C (per allele odds ratio [OR] = 3.14; P = 6.48 × 10(-11) ), a germline rare single-nucleotide polymorphism (SNP) in TP53, via imputation of a genome-wide association study of glioma (1,856 cases and 4,955 controls). We subsequently performed integrative analyses on the Cancer Genome Atlas (TCGA) data for GBM (glioblastoma multiforme) and LUAD (lung adenocarcinoma). Based on SNP data, we imputed genotypes for rs78378222 and selected individuals carrying rare risk allele (C). Using RNA sequencing data, we observed aberrant transcripts with ∼3 kb longer than normal for those individuals. Using exome sequencing data, we further showed that loss of haplotype carrying common protective allele (A) occurred somatically in GBM but not in LUAD. Our bioinformatic analysis suggests rare risk allele (C) disrupts mRNA termination, and an allelic loss of a genomic region harboring common protective allele (A) occurs during tumor initiation or progression for glioma.

Testicular germ cell tumor susceptibility associated with the UCK2 locus on chromosome 1q23.

Genome-wide association studies (GWASs) have identified multiple common genetic variants associated with an increased risk of testicular germ cell tumors (TGCTs). A previous GWAS reported a possible TGCT susceptibility locus on chromosome 1q23 in the UCK2 gene, but failed to reach genome-wide significance following replication. We interrogated this region by conducting a meta-analysis of two independent GWASs including a total of 940 TGCT cases and 1559 controls for 122 single-nucleotide polymorphisms (SNPs) on chromosome 1q23 and followed up the most significant SNPs in an additional 2202 TGCT cases and 2386 controls from four case-control studies. We observed genome-wide significant associations for several UCK2 markers, the most significant of which was for rs3790665 (PCombined = 6.0 × 10(-9)). Additional support is provided from an independent familial study of TGCT where a significant over-transmission for rs3790665 with TGCT risk was observed (PFBAT = 2.3 × 10(-3)). Here, we provide substantial evidence for the association between UCK2 genetic variation and TGCT risk.

A subset-based approach improves power and interpretation for the combined analysis of genetic association studies of heterogeneous traits.

Pooling genome-wide association studies (GWASs) increases power but also poses methodological challenges because studies are often heterogeneous. For example, combining GWASs of related but distinct traits can provide promising directions for the discovery of loci with small but common pleiotropic effects. Classical approaches for meta-analysis or pooled analysis, however, might not be suitable for such analysis because individual variants are likely to be associated with only a subset of the traits or might demonstrate effects in different directions. We propose a method that exhaustively explores subsets of studies for the presence of true association signals that are in either the same direction or possibly opposite directions. An efficient approximation is used for rapid evaluation of p values. We present two illustrative applications, one for a meta-analysis of separate case-control studies of six distinct cancers and another for pooled analysis of a case-control study of glioma, a class of brain tumors that contains heterogeneous subtypes. Both the applications and additional simulation studies demonstrate that the proposed methods offer improved power and more interpretable results when compared to traditional methods for the analysis of heterogeneous traits. The proposed framework has applications beyond genetic association studies.

The chromosome 2p21 region harbors a complex genetic architecture for association with risk for renal cell carcinoma.

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.75 × 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 with rs9679290 (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.

Genotypic variants at 2q33 and risk of esophageal squamous cell carcinoma in China: a meta-analysis of genome-wide association studies.

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 P< 5 × 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.

A comprehensive resequence-analysis of 250 kb region of 8q24.21 in men of African ancestry.

BACKGROUND: Genome-wide association studies (GWAS) have identified that a ∼1 M region centromeric to the MYC oncogene on chromosome 8q24.21 harbors at least five independent loci associated with prostate cancer risk and additional loci associated with cancers of breast, colon, bladder, and chronic lymphocytic leukemia (CLL). Because GWAS identify genetic markers that may be indirectly associated with disease, fine-mapping based on sequence analysis provides important insights into patterns of linkage disequilibrium (LD) and is critical in defining the optimal variants to nominate for biological follow-up. METHODS: To catalog variation in individuals of African ancestry, we resequenced a region (250 kb; chr8:128,050, 768­128, 300,801, hg19) containing several prostate cancer susceptibility loci as well as a locus associated with CLL. Our samples included 78 individuals from Ghana and 47 of African-Americans from Johns Hopkins University. RESULTS: After quality control metrics were applied to next-generation sequence data, 1,838 SNPs were identified. Of these, 285 were novel and not yet reported in any public database. Using genotypes derived from sequencing, we refined the LD and recombination hotspots within the region and determined a set of tag SNPs to be used in future fine-mapping studies. Based on LD, we annotated putative risk loci and their surrogates using ENCODE data, which should help guide laboratory studies. CONCLUSIONS: In comparison to the 1000 Genome Project data, we have identified additional variants that could be important in establishing priorities for future functional work designed to explain the biological basis of associations between SNPs and both prostate cancer and CLL.

A genome-wide association study of prostate cancer in West African men.

Age-adjusted mortality rates for prostate cancer are higher for African-American men compared with those of European ancestry. Recent data suggest that West African men also have elevated risk for prostate cancer relative to European men. Genetic susceptibility to prostate cancer could account for part of this difference. We conducted a genome-wide association study (GWAS) of prostate cancer in West African men in the Ghana Prostate Study. Association testing was performed using multivariable logistic regression adjusted for age and genetic ancestry for 474 prostate cancer cases and 458 population-based controls on the Illumina HumanOmni-5 Quad BeadChip. The most promising association was at 10p14 within an intron of a long non-coding RNA (lncRNA RP11-543F8.2) 360 kb centromeric of GATA3 (p = 1.29E-7). In sub-analyses, SNPs at 5q31.3 were associated with high Gleason score (≥7) cancers, the strongest of which was a missense SNP in PCDHA1 (rs34575154, p = 3.66E-8), and SNPs at Xq28 (rs985081, p = 8.66E-9) and 6q21 (rs2185710, p = 5.95E-8) were associated with low Gleason score (<7) cancers. We sought to validate our findings in silico in the African Ancestry Prostate Cancer GWAS Consortium, but only one SNP, at 10p14, replicated at p < 0.05. Of the 90 prostate cancer loci reported from studies of men of European, Asian or African-American ancestry, we were able to test 81 in the Ghana Prostate Study, and 10 of these replicated at p < 0.05. Further genetic studies of prostate cancer in West African men are needed to confirm our promising susceptibility loci.