Genome-wide meta-analysis identifies novel gender specific loci associated with thyroid antibodies level in Croatians

Autoimmune thyroid diseases (AITD) are multifactorial endocrine diseases most frequently accompanied by Tg and TPO autoantibodies. Both antibodies have a higher prevalence in females and act under a strong genetic influence.To identify novel variants underlying thyroid antibody levels, we performed GWAS meta-analysis on the plasma levels of TgAb and TPOAb in three Croatian cohorts, as well as gender specific GWAS and a bivariate analysis.No significant association was detected with the level of TgAb and TPOAb in the meta-analysis of GWAS or bivariate results for all individuals. The bivariate analysis in females only revealed a genome-wide significant association for the locus near GRIN3A (rs4457391, P = 7.76 × 10^?9). The same locus had borderline association with TPOAb levels in females (rs1935377, P = 8.58 × 10^?8).In conclusion, we identified a novel gender specific locus associated with TgAb and TPOAb levels. Our findings provide a novel insight into genetic and gender differences associated with thyroid antibodies.     

Genome-wide association study identifies a single major locus contributing to survival into old age; the APOE locus revisited

By studying the loci that contribute to human longevity, we aim to identify mechanisms that contribute to healthy aging. To identify such loci, we performed a genome-wide association study (GWAS) comparing 403 unrelated nonagenarians from long-living families included in the Leiden Longevity Study (LLS) and 1670 younger population controls. The strongest candidate SNPs from this GWAS have been analyzed in a meta-analysis of nonagenarian cases from the Rotterdam Study, Leiden 85-plus study, and Danish 1905 cohort. Only one of the 62 prioritized SNPs from the GWAS analysis (P<1×10(-4) ) showed genome-wide significance with survival into old age in the meta-analysis of 4149 nonagenarian cases and 7582 younger controls [OR=0.71 (95% CI 0.65-0.77), P=3.39 × 10(-17) ]. This SNP, rs2075650, is located in TOMM40 at chromosome 19q13.32 close to the apolipoprotein E (APOE) gene. Although there was only moderate linkage disequilibrium between rs2075650 and the ApoE ε4 defining SNP rs429358, we could not find an APOE-independent effect of rs2075650 on longevity, either in cross-sectional or in longitudinal analyses. As expected, rs429358 associated with metabolic phenotypes in the offspring of the nonagenarian cases from the LLS and their partners. In addition, we observed a novel association between this locus and serum levels of IGF-1 in women (P=0.005). In conclusion, the major locus determining familial longevity up to high age as detected by GWAS was marked by rs2075650, which tags the deleterious effects of the ApoE ε4 allele. No other major longevity locus was found.     

Genome-wide association meta-analysis identifies pleiotropic risk loci for aerodigestive squamous cell cancers

Squamous cell carcinomas (SqCC) of the aerodigestive tract have similar etiological risk factors. Although genetic risk variants for individual cancers have been identified, an agnostic, genome-wide search for shared genetic susceptibility has not been performed. To identify novel and pleotropic SqCC risk variants, we performed a meta-analysis of GWAS data on lung SqCC (LuSqCC), oro/pharyngeal SqCC (OSqCC), laryngeal SqCC (LaSqCC) and esophageal SqCC (ESqCC) cancers, totaling 13,887 cases and 61,961 controls of European ancestry. We identified one novel genome-wide significant (P_meta<5x10^-8) aerodigestive SqCC susceptibility loci in the 2q33.1 region (rs56321285, TMEM273). Additionally, three previously unknown loci reached suggestive significance (P_meta<5x10^-7): 1q32.1 (rs12133735, near MDM4), 5q31.2 (rs13181561, TMEM173) and 19p13.11 (rs61494113, ABHD8). Multiple previously identified loci for aerodigestive SqCC also showed evidence of pleiotropy in at least another SqCC site, these include: 4q23 (ADH1B), 6p21.33 (STK19), 6p21.32 (HLA-DQB1), 9p21.33 (CDKN2B-AS1) and 13q13.1(BRCA2). Gene-based association and gene set enrichment identified a set of 48 SqCC-related genes rel to DNA damage and epigenetic regulation pathways. Our study highlights the importance of cross-cancer analyses to identify pleiotropic risk loci of histology-related cancers arising at distinct anatomical sites.     

Genome-Wide Association and Trans-ethnic Meta-Analysis for Advanced Diabetic Kidney Disease: Family Investigation of Nephropathy and Diabetes (FIND)

Diabetic kidney disease (DKD) is the most common etiology of chronic kidney disease (CKD) in the industrialized world and accounts for much of the excess mortality in patients with diabetes mellitus. Approximately 45% of U.S. patients with incident end-stage kidney disease (ESKD) have DKD. Independent of glycemic control, DKD aggregates in families and has higher incidence rates in African, Mexican, and American Indian ancestral groups relative to European populations. The Family Investigation of Nephropathy and Diabetes (FIND) performed a genome-wide association study (GWAS) contrasting 6,197 unrelated individuals with advanced DKD with healthy and diabetic individuals lacking nephropathy of European American, African American, Mexican American, or American Indian ancestry. A large-scale replication and trans-ethnic meta-analysis included 7,539 additional European American, African American and American Indian DKD cases and non-nephropathy controls. Within ethnic group meta-analysis of discovery GWAS and replication set results identified genome-wide significant evidence for association between DKD and rs12523822 on chromosome 6q25.2 in American Indians (P = 5.74x10^-9). The strongest signal of association in the trans-ethnic meta-analysis was with a SNP in strong linkage disequilibrium with rs12523822 (rs955333; P = 1.31x10^-8), with directionally consistent results across ethnic groups. These 6q25.2 SNPs are located between the SCAF8 and CNKSR3 genes, a region with DKD relevant changes in gene expression and an eQTL with IPCEF1, a gene co-translated with CNKSR3. Several other SNPs demonstrated suggestive evidence of association with DKD, within and across populations. These data identify a novel DKD susceptibility locus with consistent directions of effect across diverse ancestral groups and provide insight into the genetic architecture of DKD.     

Genome-wide pathway analysis of genome-wide association studies on systemic lupus erythematosus and rheumatoid arthritis

The aim of this study was to explore candidate single nucleotide polymorphisms (SNPs) and candidate mechanisms of systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA). Two SLE genome-wide association studies (GWASs) datasets were included in this study. Meta-analysis was conducted using 737,984 SNPs in 1,527 SLE cases and 3,421 controls of European ancestry, and 4,429 SNPs that met a threshold of p?<?0.01 in a Korean RA GWAS dataset was used. ICSNPathway (identify candidate causal SNPs and pathways) analysis was applied to the meta-analysis results of the SLE GWAS datasets, and a RA GWAS dataset. The most significant result of SLE GWAS meta-analysis concerned rs2051549 in the human leukocyte antigen (HLA) region (p?=?3.36E?22). In the non-HLA region, meta-analysis identified 6 SNPs associated with SLE with genome-wide significance (STAT4, TNPO3, BLK, FAM167A, and IRF5). ICSNPathway identified five candidate causal SNPs and 13 candidate causal pathways. This pathway-based analysis provides three hypotheses of the biological mechanism involved. First, rs8084 and rs7192?→?HLA-DRA?→?bystander B cell activation. Second, rs1800629?→?TNF?→?cytokine network. Third, rs1150752 and rs185819?→?TNXB?→?collagen metabolic process. ICSNPathway analysis identified three candidate causal non-HLA SNPs and four candidate causal pathways involving the PADI4, MTR, PADI2, and TPH2 genes of RA. We identified five candidate SNPs and thirteen pathways, involving bystander B cell activation, cytokine network, and collagen metabolic processing, which may contribute to SLE susceptibility, and we revealed candidate causal non-HLA SNPs, genes, and pathways of RA.     

Meta-analysis of two Chinese populations identifies an autoimmune disease risk allele in 22q11.21 as associated with systemic lupus erythematosus

IntroductionSystemic lupus erythematosus (SLE) is a heterogeneous disease with a diverse spectrum of clinical symptoms, ranging from skin rash to end-organ damage. 22q11.21 has been identified as a susceptibility region for several autoimmune diseases, including SLE. However, detailed information for SLE association and the underlying functional mechanism(s) is still lacking.MethodsThrough meta-analysis of two genome-wide association studies (GWAS) on Han Chinese populations, comprising a total of 1,659 cases and 3,398 controls matched geographically, we closely examined the 22q11.21 region, especially on the reported single-nucleotide polymorphisms (SNPs) associated with different autoimmune diseases and their relationships. We further replicated the most significant associations of SNPs with SLE using 2,612 cases and 2,323 controls of Asian ancestry.ResultsAll reported SNPs in the 22q11.21 region with different autoimmune diseases were examined using the two GWAS data and meta-analysis results, and supportive evidence of association with SLE was found (meta-analysis: P_meta ≤ 7.27E-05), which might require further investigation. SNP rs2298428 was identified as the most significant SNP associated with SLE in this region (P_meta =2.70E-09). It showed independent effects through both stepwise and conditional logistic regression, and there is no evidence of other independent association signals for SLE in this region. The association of rs2298428 was further replicated in three cohorts from Hong Kong, Anhui and Thailand comprising a total of 2,612 cases and 2,323 controls (joint analysis of GWAS and replication result: P_all =1.31E-11, odds ratio =1.23). SNP rs2298428 was shown to be an expression quantitative locus for UBE2L3 gene in different cell types, with the risk allele (T) being correlated with higher expression of UBE2L3. This is consistent with earlier reports on higher expression of UBE2L3 in patients with SLE.ConclusionsAssociation with distinct autoimmune diseases highlights the significance of this region in autoreactive responses and potentially shared functional mechanisms in these diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0577-6) contains supplementary material, which is available to authorized users.     

Common Variation at 1q24.1 (ALDH9A1) Is a Potential Risk Factor for Renal Cancer

So far six susceptibility loci for renal cell carcinoma (RCC) have been discovered by genome-wide association studies (GWAS). To identify additional RCC common risk loci, we performed a meta-analysis of published GWAS (totalling 2,215 cases and 8,566 controls of Western-European background) with imputation using 1000 Genomes Project and UK10K Project data as reference panels and followed up the most significant association signals [22 single nucleotide polymorphisms (SNPs) and 3 indels in eight genomic regions] in 383 cases and 2,189 controls from The Cancer Genome Atlas (TCGA). A combined analysis identified a promising susceptibility locus mapping to 1q24.1 marked by the imputed SNP rs3845536 (P _combined =2.30x10^-8). Specifically, the signal maps to intron 4 of the ALDH9A1 gene (aldehyde dehydrogenase 9 family, member A1). We further evaluated this potential signal in 2,461 cases and 5,081 controls from the International Agency for Research on Cancer (IARC) GWAS of RCC cases and controls from multiple European regions. In contrast to earlier findings no association was shown in the IARC series (P=0.94; P _combined =2.73x10^-5). While variation at 1q24.1 represents a potential risk locus for RCC, future replication analyses are required to substantiate our observation.     

A Genome-Wide Association Study of a Biomarker of Nicotine Metabolism

Individuals with fast nicotine metabolism typically smoke more and thus have a greater risk for smoking-induced diseases. Further, the efficacy of smoking cessation pharmacotherapy is dependent on the rate of nicotine metabolism. Our objective was to use nicotine metabolite ratio (NMR), an established biomarker of nicotine metabolism rate, in a genome-wide association study (GWAS) to identify novel genetic variants influencing nicotine metabolism. A heritability estimate of 0.81 (95% CI 0.70–0.88) was obtained for NMR using monozygotic and dizygotic twins of the FinnTwin cohort. We performed a GWAS in cotinine-verified current smokers of three Finnish cohorts (FinnTwin, Young Finns Study, FINRISK2007), followed by a meta-analysis of 1518 subjects, and annotated the genome-wide significant SNPs with methylation quantitative loci (meQTL) analyses. We detected association on 19q13 with 719 SNPs exceeding genome-wide significance within a 4.2 Mb region. The strongest evidence for association emerged for CYP2A6 (min p = 5.77E-86, in intron 4), the main metabolic enzyme for nicotine. Other interesting genes with genome-wide significant signals included CYP2B6, CYP2A7, EGLN2, and NUMBL. Conditional analyses revealed three independent signals on 19q13, all located within or in the immediate vicinity of CYP2A6. A genetic risk score constructed using the independent signals showed association with smoking quantity (p = 0.0019) in two independent Finnish samples. Our meQTL results showed that methylation values of 16 CpG sites within the region are affected by genotypes of the genome-wide significant SNPs, and according to causal inference test, for some of the SNPs the effect on NMR is mediated through methylation. To our knowledge, this is the first GWAS on NMR. Our results enclose three independent novel signals on 19q13.2. The detected CYP2A6 variants explain a strikingly large fraction of variance (up to 31%) in NMR in these study samples. Further, we provide evidence for plausible epigenetic mechanisms influencing NMR.     

Genome-wide association study and meta-analysis of intraocular pressure

Elevated intraocular pressure (IOP) is a major risk factor for glaucoma and is influenced by genetic and environmental factors. Recent genome-wide association studies (GWAS) reported associations with IOP at TMCO1 and GAS7, and with primary open-angle glaucoma (POAG) at CDKN2B-AS1, CAV1/CAV2, and SIX1/SIX6. To identify novel genetic variants and replicate the published findings, we performed GWAS and meta-analysis of IOP in >6,000 subjects of European ancestry collected in three datasets: the NEI Glaucoma Human genetics collaBORation, GLAUcoma Genes and ENvironment study, and a subset of the Age-related Macular Degeneration-Michigan, Mayo, AREDS and Pennsylvania study. While no signal achieved genome-wide significance in individual datasets, a meta-analysis identified significant associations with IOP at TMCO1 (rs7518099-G, p = 8.0 × 10^?8). Focused analyses of five loci previously reported for IOP and/or POAG, i.e., TMCO1, CDKN2B-AS1, GAS7, CAV1/CAV2, and SIX1/SIX6, revealed associations with IOP that were largely consistent across our three datasets, and replicated the previously reported associations in both effect size and direction. These results confirm the involvement of common variants in multiple genomic regions in regulating IOP and/or glaucoma risk.     

Use of a multiethnic approach to identify rheumatoid- arthritis-susceptibility loci, 1p36 and 17q12

We have previously shown that rheumatoid arthritis (RA) risk alleles overlap between different ethnic groups. Here, we utilize a multiethnic approach to show that we can effectively discover RA risk alleles. Thirteen putatively associated SNPs that had not yet exceeded genome-wide significance (p < 5 × 10(-8)) in our previous RA genome-wide association study (GWAS) were analyzed in independent sample sets consisting of 4,366 cases and 17,765 controls of European, African American, and East Asian ancestry. Additionally, we conducted an overall association test across all 65,833 samples (a GWAS meta-analysis plus the replication samples). Of the 13 SNPs investigated, four were significantly below the study-wide Bonferroni corrected p value threshold (p < 0.0038) in the replication samples. Two SNPs (rs3890745 at the 1p36 locus [p = 2.3 × 10(-12)] and rs2872507 at the 17q12 locus [p = 1.7 × 10(-9)]) surpassed genome-wide significance in all 16,659 RA cases and 49,174 controls combined. We used available GWAS data to fine map these two loci in Europeans and East Asians, and we found that the same allele conferred risk in both ethnic groups. A series of bioinformatic analyses identified TNFRSF14-MMEL1 at the 1p36 locus and IKZF3-ORMDL3-GSDMB at the 17q12 locus as the genes most likely associated with RA. These findings demonstrate empirically that a multiethnic approach is an effective strategy for discovering RA risk loci, and they suggest that combining GWASs across ethnic groups represents an efficient strategy for gaining statistical power.     

Genome-wide association study of african and European americans implicates multiple shared and ethnic specific Loci in sarcoidosis susceptibility

Sarcoidosis is a systemic inflammatory disease characterized by the formation of granulomas in affected organs. Genome-wide association studies (GWASs) of this disease have been conducted only in European population. We present the first sarcoidosis GWAS in African Americans (AAs, 818 cases and 1,088 related controls) followed by replication in independent sets of AAs (455 cases and 557 controls) and European Americans (EAs, 442 cases and 2,284 controls). We evaluated >6 million SNPs either genotyped using the Illumina Omni1-Quad array or imputed from the 1000 Genomes Project data. We identified a novel sarcoidosis-associated locus, NOTCH4, that reached genome-wide significance in the combined AA samples (rs715299, P _AA-meta = 6.51×10⁻¹⁰) and demonstrated the independence of this locus from others in the MHC region in the same sample. We replicated previous European GWAS associations within HLA-DRA, HLA-DRB5, HLA-DRB1, BTNL2, and ANXA11 in both our AA and EA datasets. We also confirmed significant associations to the previously reported HLA-C and HLA-B regions in the EA but not AA samples. We further identified suggestive associations with several other genes previously reported in lung or inflammatory diseases.     

Genome-wide association study identified novel genetic variant on SLC45A3 gene associated with serum levels prostate-specific antigen (PSA) in a Chinese population

Prostate-specific antigen (PSA) is a commonly used cancer biomarker for prostate cancer, and is often included as part of routine physical examinations in China. Serum levels of PSA may be influenced by genetic factors as well as other factors. A genome-wide association study (GWAS) conducted in a European population successfully identified six genetic loci that were significantly associated with PSA level. In this study, we aimed to identify common genetic variants that are associated with serum level of PSA in a Chinese population. We also evaluated the effects of those variants by creating personalized PSA cutoff values. A two-stage GWAS of PSA level was performed among men age 20–69 years and self-reported cancer-free participants that underwent routine physical examinations at several hospitals in Guangxi Province, China. Single nucleotide polymorphisms (SNPs) significantly associated with PSA levels in the first stage of sample (N = 1,999) were confirmed in the second stage of sample (N = 1,496). Multivariate linear regression was used to assess the independent contribution of confirmed SNPs and known covariates, such as age, to the level of PSA. SNPs in three regions were significantly associated with levels of PSA in this two-stage GWAS, and had combined P values between 4.62 × 10^?17 and 6.45 × 10^?37. The three regions are located on 1q32.1 at SLC45A3, 10q11.23 at MSMB, and 19q13.33 at KLK3. The region 1q32.1 at SLC45A3 was identified as a novel locus. Genetic variants contributed significantly more to the variance of PSA level than known covariates such as age. Personalized cutoff values of serum PSA, calculated based on the inheritance of these associated SNPs, differ considerably among individuals. Identification of these genetic markers provides new insight into the molecular mechanisms of PSA. Taking individual variation into account, these genetic variants may improve the performance of PSA to predict prostate cancer.     

Genome-wide association analysis of canine atopic dermatitis and identification of disease related SNPs

In humans, genome-wide association studies (GWAS) have been shown to be an effective and thorough approach for identifying polymorphisms associated with disease phenotypes. Here, we describe the first study to perform a genome-wide association study in canine atopic dermatitis (cAD) using the Illumina Canine SNP20 array, containing 22,362 single-nucleotide polymorphisms (SNPs). The aim of the study was to identify SNPs associated with cAD using affected and unaffected Golden Retrievers. Further validation studies were performed for potentially associated SNPs using Sequenom genotyping of larger numbers of cases and controls across eight breeds (Boxer, German Shepherd Dog, Labrador, Golden Retriever, Shiba Inu, Shih Tzu, Pit Bull, and West Highland White Terriers). Using meta-analysis, two SNPs were associated with cAD in all breeds tested. RS22114085 was identified as a susceptibility locus (p?=?0.00014, odds ratio?=?2) and RS23472497 as a protective locus (p?=?0.0015, odds ratio?=?0.6). Both of these SNPs were located in intergenic regions, and their effects have been demonstrated to be independent of each other, highlighting that further fine mapping and resequencing is required of these areas. Further, 12 SNPs were validated by Sequenom genotyping as associated with cAD, but these were not associated with all breeds. This study suggests that GWAS will be a useful approach for identifying genetic risk factors for cAD. Given the clinical heterogeneity within this condition and the likelihood that the relative genetic effect sizes are small, greater sample sizes and further studies will be required.     

How Genome-Wide SNP-SNP Interactions Relate to Nasopharyngeal Carcinoma Susceptibility

This study is the first to use genome-wide association study (GWAS) data to evaluate the multidimensional genetic architecture underlying nasopharyngeal cancer. Since analysis of data from GWAS confirms a close and consistent association between elevated risk for nasopharyngeal carcinoma (NPC) and major histocompatibility complex class 1 genes, our goal here was to explore lesser effects of gene-gene interactions. We conducted an exhaustive genome-wide analysis of GWAS data of NPC, revealing two-locus interactions occurring between single nucleotide polymorphisms (SNPs), and identified a number of suggestive interaction loci which were missed by traditional GWAS analyses. Although none of the interaction pairs we identified passed the genome-wide Bonferroni-adjusted threshold for significance, using independent GWAS data from the same population (Stage 2), we selected 66 SNP pairs in 39 clusters with P<0.01. We identified that in several chromosome regions, multiple suggestive interactions group to form a block-like signal, effectively reducing the rate of false discovery. The strongest cluster of interactions involved the CREB5 gene and a SNP rs1607979 on chromosome 17q22 (P = 9.86×10⁻¹¹) which also show trans-expression quantitative loci (eQTL) association in Chinese population. We then detected a complicated cis-interaction pattern around the NPC-associated HLA-B locus, which is immediately adjacent to copy-number variations implicated in male susceptibility for NPC. While it remains to be seen exactly how and to what degree SNP-SNP interactions such as these affect susceptibility for nasopharyngeal cancer, future research on these questions holds great promise for increasing our understanding of this disease’s genetic etiology, and possibly also that of other gene-related cancers.     

A meta-analysis of two genome-wide association studies to identify novel loci for maximum number of alcoholic drinks

Maximum number of alcoholic drinks consumed in a 24-h period (maxdrinks) is a heritable (>50 %) trait and is strongly correlated with vulnerability to excessive alcohol consumption and subsequent alcohol dependence (AD). Several genome-wide association studies (GWAS) have studied alcohol dependence, but few have concentrated on excessive alcohol consumption. We performed two GWAS using maxdrinks as an excessive alcohol consumption phenotype: one in 118 extended families (N = 2,322) selected from the Collaborative Study on the Genetics of Alcoholism (COGA), and the other in a case–control sample (N = 2,593) derived from the Study of Addiction: Genes and Environment (SAGE). The strongest association in the COGA families was detected with rs9523562 (p = 2.1 × 10^?6) located in an intergenic region on chromosome 13q31.1; the strongest association in the SAGE dataset was with rs67666182 (p = 7.1 × 10^?7), located in an intergenic region on chromosome 8. We also performed a meta-analysis with these two GWAS and demonstrated evidence of association in both datasets for the LMO1 (p = 7.2 × 10^?7) and PLCL1 genes (p = 4.1 × 10^?6) with maxdrinks. A variant in AUTS2 and variants in INADL, C15orf32 and HIP1 that were associated with measures of alcohol consumption in a meta-analysis of GWAS studies and a GWAS of alcohol consumption factor score also showed nominal association in the current meta-analysis. The present study has identified several loci that warrant further examination in independent samples. Among the top SNPs in each of the dataset (p ≤ 10^?4) far more showed the same direction of effect in the other dataset than would be expected by chance (p = 2 × 10^?3, 3 × 10^?6), suggesting that there are true signals among these top SNPs, even though no SNP reached genome-wide levels of significance.     

Genome-wide association study of body weight in Wenshang Barred chicken based on the SLAF-seq technology

Chicken body weight (BW) is an economically important trait, and many studies have been conducted on genetic selection for BW. However, previous studies have detected functional chromosome mutations or regions using gene chips. The present study used the specific-locus amplified fragment sequencing (SLAF-seq) technology to perform a genome-wide association study (GWAS) on purebred Wengshang Barred chicken. A total of 1,286,715 single-nucleotide polymorphisms (SNPs) were detected, and 175,211 SNPs were selected as candidate SNPs for genome-wide association analysis using TASSEL general linear models. Six SNP markers reached genome-wide significance. Of these, rs732048524, rs735522839, rs738991545, and rs15837818 were significantly associated with body weight at 28 days (BW28), while rs314086457 and rs315694878 were significantly associated with BW120. These SNPs are close to seven genes (PRSS23, ME3, FAM181B, NABP1, SDPR, TSSK6L2, and RBBP8). Moreover, 24 BW-associated SNPs reached “suggestive” genome-wide significance. Of these, 6, 13, 1, and 4 SNPs were associated with BW28, BW56, BW80, and BW120, respectively. These results would enrich the studies on BW and promote the use of Chinese chicken, especially the Wenshang Barred chicken.     

Identification of a Potential Regulatory Variant for Colorectal Cancer Risk Mapping to Chromosome 5q31.1: A Post-GWAS Study

Large-scale genome-wide association studies (GWAS) have established chromosome 5q31.1 as a susceptibility locus for colorectal cancer (CRC), which was still lack of causal genetic variants. We searched potentially regulatory single nucleotide polymorphisms (SNPs) in the overlap region between linkage disequilibrium (LD) block of 5q31.1 and regulatory elements predicted by histone modifications, then tested their association with CRC via a case-control study. Among three candidate common variants, we found rs17716310 conferred significantly (heterozygous model: OR = 1.273, 95% confidence interval (95%CI) = 1.016–1.595, P = 0.036) and marginally (dominant model: OR = 1.238, 95%CI = 1.000–1.532, P = 0.050) increase risk for CRC in a Chinese population including 695 cases and 709 controls. This variation was suggested to be regulatory altering the activity of enhancer that control PITX1 expression. Using epigenetic information such as chromatin immunoprecipitation-sequencing (ChIP-seq) data might help researchers to interpret the results of GWAS and locate causal variants for diseases in post-GWAS era.