Identification of HAVCR1 gene haplotypes associated with mRNA expression levels and susceptibility to autoimmune diseases

Human HAVCR1 gene maps on 5q33.2, a region linked with susceptibility to allergic and autoimmune diseases. The aims of the present study were to define the haplotypes of HAVCR1 gene taking into account both HapMap Project SNP haplotypes and exon 4 variants, to investigate a possible relationship between these haplotypes and mRNA expression levels, and to assess whether HAVCR1 gene is involved in susceptibility to rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE). Genotyping of three ins/del variants in the exon 4 was performed by fragment length analysis. Five tag SNPs genotypes and mRNA levels were determined using TaqMan assays. We defined four major haplotypes in our population: the two major haplotypes (named haplotypes A and B) bear both the 5383_5397del variant and the two most common SNP sets found in the CEU population. Quantification analysis revealed that genotype B/B had the highest median of mRNA expression levels (vs. BX + XX, p < 0.0001). Additionally, frequency of the genotype BB was significantly higher in RA patients than in controls (12.3 vs. 5.9% in controls, p = 0.0046, p _c = 0.014, OR = 2.23, 95% CI 1.23–4.10). Our results support a relationship between HAVCR1 haplotypes and mRNA expression levels, and suggest an association of this gene with autoimmune diseases.     

Targeted exon sequencing fails to identify rare coding variants with large effect in rheumatoid arthritis

Introduction

Although it has been suggested that rare coding variants could explain the substantial missing heritability, very few sequencing studies have been performed in rheumatoid arthritis (RA). We aimed to identify novel functional variants with rare to low frequency using targeted exon sequencing of RA in Korea.

Methods

We analyzed targeted exon sequencing data of 398 genes selected from a multifaceted approach in Korean RA patients (n = 1,217) and controls (n = 717). We conducted a single-marker association test and a gene-based analysis of rare variants. For meta-analysis or enrichment tests, we also used ethnically matched independent samples of Korean genome-wide association studies (GWAS) (n = 4,799) or immunochip data (n = 4,722).

Results

After stringent quality control, we analyzed 10,588 variants of 398 genes from 1,934 Korean RA case controls. We identified 13 nonsynonymous variants with nominal association in single-variant association tests. In a meta-analysis, we did not find any novel variant with genome-wide significance for RA risk. Using a gene-based approach, we identified 17 genes with nominal burden signals. Among them, VSTM1 showed the greatest association with RA (P = 7.80 × 10⁻⁴). In the enrichment test using Korean GWAS, although the significant signal appeared to be driven by total genic variants, we found no evidence for enriched association of coding variants only with RA.

Conclusions

We were unable to identify rare coding variants with large effect to explain the missing heritability for RA in the current targeted resequencing study. Our study raises skepticism about exon sequencing of targeted genes for complex diseases like RA.

Electronic supplementary material

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

A genome-wide association study of severe teenage acne in European Americans

Despite the family aggregation of severe teenage acne, the genetic basis of this common skin condition remains unclear. We conducted a genome-wide association study (GWAS) on severe teenage acne in 928 European Americans. The SNP rs4133274 on chromosome 8q24 (72 kb upstream of MYC) revealed the most significant association with severe teenage acne (p value = 1.7 × 10^?6). The variant allele of this SNP (G allele) was associated with an increased risk of severe teenage acne with odds ratio of 4.01 (95 % confidence interval = 2.37–6.82). Upon further replication, our findings suggest new genetic basis of acne and may explain the association between acne and cancer risk observed in the epidemiological studies.     

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 studies identified novel loci for non-high-density lipoprotein cholesterol and its postprandial lipemic response

Non-high-density lipoprotein cholesterol(NHDL) is an independent and superior predictor of CVD risk as compared to low-density lipoprotein alone. It represents a spectrum of atherogenic lipid fractions with possibly a distinct genomic signature. We performed genome-wide association studies (GWAS) to identify loci influencing baseline NHDL and its postprandial lipemic (PPL) response. We carried out GWAS in 4,241 participants of European descent. Our discovery cohort included 928 subjects from the Genetics of Lipid-Lowering Drugs and Diet Network Study. Our replication cohorts included 3,313 subjects from the Heredity and Phenotype Intervention Heart Study and Family Heart Study. A linear mixed model using the kinship matrix was used for association tests. The best association signal was found in a tri-genic region at RHOQ-PIGF-CRIPT for baseline NHDL (lead SNP rs6544903, discovery p = 7e?7, MAF = 2 %; validation p = 6e?4 at 0.1 kb upstream neighboring SNP rs3768725, and 5e?4 at 0.7 kb downstream neighboring SNP rs6733143, MAF = 10 %). The lead and neighboring SNPs were not perfect surrogate proxies to each other (D′ = 1, r ² = 0.003) but they seemed to be partially dependent (likelihood ration test p = 0.04). Other suggestive loci (discovery p < 1e?6) included LOC100419812 and LOC100288337 for baseline NHDL, and LOC100420502 and CDH13 for NHDL PPL response that were not replicated (p > 0.01). The current and first GWAS of NHDL yielded an interesting common variant in RHOQ-PIGF-CRIPT influencing baseline NHDL levels. Another common variant in CDH13 for NHDL response to dietary high-fat intake challenge was also suggested. Further validations for both loci from large independent studies, especially interventional studies, are warranted.     

Hierarchical modeling identifies novel lung cancer susceptibility variants in inflammation pathways among 10,140 cases and 11,012 controls

Recent evidence suggests that inflammation plays a pivotal role in the development of lung cancer. In this study, we used a two-stage approach to investigate associations between genetic variants in inflammation pathways and lung cancer risk based on genome-wide association study (GWAS) data. A total of 7,650 sequence variants from 720 genes relevant to inflammation pathways were identified using keyword and pathway searches from Gene Cards and Gene Ontology databases. In Stage 1, six GWAS datasets from the International Lung Cancer Consortium were pooled (4,441 cases and 5,094 controls of European ancestry), and a hierarchical modeling (HM) approach was used to incorporate prior information for each of the variants into the analysis. The prior matrix was constructed using (1) role of genes in the inflammation and immune pathways; (2) physical properties of the variants including the location of the variants, their conservation scores and amino acid coding; (3) LD with other functional variants and (4) measures of heterogeneity across the studies. HM affected the priority ranking of variants particularly among those having low prior weights, imprecise estimates and/or heterogeneity across studies. In Stage 2, we used an independent NCI lung cancer GWAS study (5,699 cases and 5,818 controls) for in silico replication. We identified one novel variant at the level corrected for multiple comparisons (rs2741354 in EPHX2 at 8q21.1 with p value = 7.4 × 10^?6), and confirmed the associations between TERT (rs2736100) and the HLA region and lung cancer risk. HM allows for prior knowledge such as from bioinformatic sources to be incorporated into the analysis systematically, and it represents a complementary analytical approach to the conventional GWAS analysis.     

Resequencing Three Candidate Genes for Major Depressive Disorder in a Dutch Cohort

Major depressive disorder (MDD) is a psychiatric disorder, characterized by periods of low mood of more than two weeks, loss of interest in normally enjoyable activities and behavioral changes. MDD is a complex disorder and does not have a single genetic cause. In 2009 a genome wide association study (GWAS) was performed on the Dutch GAIN-MDD cohort. Many of the top signals of this GWAS mapped to a region spanning the gene PCLO, and the non-synonymous coding single nucleotide polymorphism (SNP) rs2522833 in the PCLO gene became genome wide significant after post-hoc analysis. We performed resequencing of PCLO, GRM7, and SLC6A4 in 50 control samples from the GAIN-MDD cohort, to detect new genomic variants. Subsequently, we genotyped these variants in the entire GAIN-MDD cohort and performed association analysis to investigate if rs2522833 is the causal variant or simply in linkage disequilibrium with a more associated variant. GRM7 and SLC6A4 are both candidate genes for MDD from literature. We aimed to gather more evidence that rs2522833 is indeed the causal variant in the GAIN-MDD cohort or to find a previously undetected common variant in either PCLO, GRM7, or SLC6A4 with a higher association in this cohort. After next generation sequencing and association analysis we excluded the possibility of an undetected common variant to be more associated. For neither PCLO nor GRM7 we found a more associated variant. For SLC6A4, we found a new SNP that showed a lower P-value (P = 0.07) than in the GAIN-MDD GWAS (P = 0.09). However, no evidence for genome-wide significance was found. Although we did not take into account rare variants, we conclude that our results provide further support for the hypothesis that the non-synonymous coding SNP rs2522833 in the PCLO gene is indeed likely to be the causal variant in the GAIN-MDD cohort.     

The majority of the genetic risk for Paget’s disease of bone is explained by genetic variants close to the CSF1, OPTN, TM7SF4, and TNFRSF11A genes

Paget’s disease of bone (PDB) is one of the most frequent metabolic bone disorders (1–5%), next to osteoporosis, affecting individuals above age 55. Sequestosome1 mutations explain a part of the PDB patients, but still the disease pathogenesis in the remaining PDB patients is largely unknown. Therefore, association studies investigating the relationship between genetic polymorphisms and sporadic PDB have been performed to find the genetic risk variants. Previously such studies indicated a role of the OPG and RANK gene. The latter was recently confirmed in a genome-wide association study (GWAS) which also indicated the involvement of chromosomal regions harbouring the CSF1 and OPTN gene. In this study, we sought to replicate these findings in a Belgian and a Dutch population. Similar significant results were obtained for the single nucleotide polymorphisms and the haplotypes. The most significant results are found in the CSF1 gene region, followed by the OPTN and TNFRSF11A gene region (p values ranging from 1.3 × 10^?4 to 3.8 × 10^?8, OR = 1.523–1.858). We next obtained significant association with a polymorphism from the chromosomal region around the TM7SF4 gene (p = 2.7 × 10^?3, OR = 1.427), encoding DC-STAMP which did not reach genome-wide significance in the GWAS, but based on its function in osteoclasts it can be considered a strong candidate gene. After meta-analysis with the GWAS data, p values ranged between 2.6 × 10^?4 and 8.8 × 10^?32. The calculated cumulative population attributable risk of these four loci turned out to be about 67% in our two populations, indicating that most of the genetic risk for PDB is coming from genetic variants close to these four genes.     

Integrated analysis of direct and proxy genome wide association studies highlights polygenicity of Alzheimer’s disease outside of the APOE region

Recent meta-analyses combining direct genome-wide association studies (GWAS) with those of family history (GWAX) have indicated very low SNP heritability of Alzheimer’s disease (AD). These low estimates may call into question the prospects of continued progress in genetic discovery for AD within the spectrum of common variants. We highlight dramatic downward biases in previous methods, and we validate a novel method for the estimation of SNP heritability via integration of GWAS and GWAX summary data. We apply our method to investigate the genetic architecture of AD using GWAX from UK Biobank and direct case-control GWAS from the International Genomics of Alzheimer’s Project (IGAP). We estimate the liability scale common variant SNP heritability of Clinical AD outside of APOE region at ~7–11%, and we project the corresponding estimate for AD pathology to be up to approximately 23%. We estimate that nearly 90% of common variant SNP heritability of Clinical AD exists outside the APOE region. Rare variants not tagged in standard GWAS may account for additional variance. Our results indicate that, while GWAX for AD in UK Biobank may result in greater attenuation of genetic effects beyond that conventionally assumed, it does not introduce appreciable contamination of signal by genetically distinct traits relative to direct case-control GWAS in IGAP. Genetic risk for AD represents a strong effect of APOE superimposed upon a highly polygenic background.     

A genome-wide association study on copy-number variation identifies a 11q11 loss as a candidate susceptibility variant for colorectal cancer

Colorectal cancer (CRC) is a complex disease, and therefore its development is determined by the combination of both environmental factors and genetic variants. Although genome-wide association studies (GWAS) of SNP variation have conveniently identified 20 genetic variants so far, a significant proportion of the observed heritability is yet to be explained. Common copy-number variants (CNVs) are one of the most important genomic sources of variability, and hence a potential source to explain part of this missing genetic fraction. Therefore, we have performed a GWAS on CNVs to explore the relationship between common structural variation and CRC development. Phase 1 of the GWAS consisted of 881 cases and 667 controls from a Spanish cohort. Copy-number status was validated by quantitative PCR for each of those common CNVs potentially associated with CRC in phase I. Subsequently, SNPs were chosen as proxies for the validated CNVs for phase II replication (1,342 Spanish cases and 1,874 Spanish controls). Four common CNVs were found to be associated with CRC and were further replicated in Phase II. Finally, we found that SNP rs1944682, tagging a 11q11 CNV, was nominally associated with CRC susceptibility (p value = 0.039; OR = 1.122). This locus has been previously related to extreme obesity phenotypes, which could suggest a relationship between body weight and CRC susceptibility.     

Low Frequency Variants,Collapsed Based on Biological Knowledge,Uncover Complexity of Population Stratification in 1000 Genomes Project Data

Analyses investigating low frequency variants have the potential for explaining additional genetic heritability of many complex human traits. However, the natural frequencies of rare variation between human populations strongly confound genetic analyses. We have applied a novel collapsing method to identify biological features with low frequency variant burden differences in thirteen populations sequenced by the 1000 Genomes Project. Our flexible collapsing tool utilizes expert biological knowledge from multiple publicly available database sources to direct feature selection. Variants were collapsed according to genetically driven features, such as evolutionary conserved regions, regulatory regions genes, and pathways. We have conducted an extensive comparison of low frequency variant burden differences (MAF<0.03) between populations from 1000 Genomes Project Phase I data. We found that on average 26.87% of gene bins, 35.47% of intergenic bins, 42.85% of pathway bins, 14.86% of ORegAnno regulatory bins, and 5.97% of evolutionary conserved regions show statistically significant differences in low frequency variant burden across populations from the 1000 Genomes Project. The proportion of bins with significant differences in low frequency burden depends on the ancestral similarity of the two populations compared and types of features tested. Even closely related populations had notable differences in low frequency burden, but fewer differences than populations from different continents. Furthermore, conserved or functionally relevant regions had fewer significant differences in low frequency burden than regions under less evolutionary constraint. This degree of low frequency variant differentiation across diverse populations and feature elements highlights the critical importance of considering population stratification in the new era of DNA sequencing and low frequency variant genomic analyses.     

Association of CD247 Polymorphisms with Rheumatoid Arthritis: A Replication Study and a Meta-Analysis

Given the role of CD247 in the response of the T cells, its entailment in autoimmune diseases and in order to better clarify the role of this gene in RA susceptibility, we aimed to analyze CD247 gene variants previously associated with other autoimmune diseases (rs1052237, rs2056626 and rs864537) in a large independent European Caucasian population. However, no evidence of association was found for the analyzed CD247 single-nucleotide polymorphisms (SNPs) with RA and with the presence/absence of anti-cyclic citrullinated polypeptide. We performed a meta-analysis including previously published GWAS data from the rs864537 variant, revealing an overall genome-wide significant association between this CD247 SNP and RA with anti-CCP (OR = 0.90, CI 95% = 0.87–0.93, P_overall = 2.1×10⁻¹⁰). Our results show for first time a GWAS-level association between this CD247 polymorphism and RA risk.     

Identification of rare genetic variants in novel loci associated with Paget’s disease of bone

In genome-wide association studies, single nucleotide polymorphisms located in five novel loci were associated with PDB. We aimed at identifying rare genetic variants of candidate genes located in these loci and search for genetic association with PDB in the French-Canadian population. Exons, promoter and exon–intron junctions from patients with familial PDB and healthy individuals were sequenced in candidate genes, located within novel loci associated with PDB in our population. Rare variant was defined by a minor allele frequency <0.05 or absent from dbSNP (NCBI). We sequenced seven genes in 1p13 locus, three genes in 7q33, three genes in 8q22, and five genes in 15q24 locus. We identified 126 rare variants in at least one patient with PDB of whom 55 were located in 1p13 locus, 32 in 7q33, 10 in 8q22 and 29 in 15q24 locus. We located 71 of these 126 rare variants in an intron, 30 in an exon and 9 in an untranslated region. 60 % of these variants were located in functionally relevant gene regions. Among the 12 missense rare variants in PDB, two (rs62620995 in TM7SF4; rs62641691 in CD276) were predicted to be damaging by in silico analysis tools. Rs62620995, which altered a conserved amino acid (p.Leu397Phe) in the TM7SF4 gene, encoding the DC-STAMP protein involved in osteoclastogenesis through RANK signaling pathway, was found to have a marginal association with PDB (p = 0.09). Rs35500845, located in the CTHRC1 gene, which encodes a regulator of collagen matrix deposition, was also associated with PDB in the French-Canadian population (p = 0.046).     

A Functional TNFAIP2 3'-UTR rs8126 Genetic Polymorphism Contributes to Risk of Esophageal Squamous Cell Carcinoma

Background

Accumulated evidences demonstrated that single nucleotide polymorphisms (SNPs) in mRNA 3''-untranslated region (3''-UTR) may impact microRNAs (miRNAs)-mediated expression regulation of oncogenes and tumor suppressors. There is a TNFAIP2 3''-UTR rs8126 T>C genetic variant which has been proved to be associated with head and neck cancer susceptibility. This SNP could disturb binding of miR-184 with TNFAIP2 mRNA and influence TNFAIP2 regulation. However, it is still unclear how this polymorphism is involved in development of esophageal squamous cell carcinoma (ESCC). Therefore, we hypothesized that the functional TNFAIP2 rs8126 SNP may affect TNFAIP2 expression and, thus, ESCC risk.

Methods

We investigated the association between the TNFAIP2 rs8126 variant and ESCC risk as well as the functional relevance on TNFAIP2 expression in vivo. Genotypes were determined in a case-control set consisted of 588 ESCC patients and 600 controls. The allele-specific regulation on TNFAIP2 expression by the rs8126 SNP was examined in normal and cancerous tissue specimens of esophagus.

Results

We found that individuals carrying the rs8126 CC or CT genotype had an OR of 1.89 (95%CI  = 1.23–2.85, P = 0.003) or 1.38 (95%CI  = 1.05–1.73, P = 0.017) for developing ESCC in Chinese compared with individual carrying the TT genotype. Carriers of the rs8126 CC and CT genotypes had significantly lower TNFAIP2 mRNA levels than those with the TT genotypes in normal esophagus tissues (P<0.05).

Conclusions

Our data demonstrate that functional TNFAIP2 rs8126 genetic variant is a ESCC susceptibility SNP. These results support the hypothesis that genetic variants interrupting miRNA-mediated gene regulation might be important genetic modifiers of cancer risk.     

BET1L and TNRC6B associate with uterine fibroid risk among European Americans

Uterine fibroid (UFs) affect 77 % of women by menopause and account for $9.4 billion in healthcare costs each year. Although UFs are heritable, genetic risk is poorly understood. The first genome-wide association study (GWAS) of UFs was recently performed in a Japanese population, with reported genome-wide significance for single nucleotide polymorphisms (SNPs) across three chromosomal regions. We tested these SNPs for association with UFs in US cohorts. Women were enrolled in the Right from the Start (RFTS) cohort and the BioVU DNA repository. UF status in both cohorts was determined by pelvic imaging. We tested 65 candidate and haplotype-tagging SNPs for association with UFs presence using logistic regression in RFTS and the top three GWAS-associated SNPs in BioVU. We also combined association results from both cohorts using meta-analysis. 1,086 European American (EA) cases and 1,549 controls were examined. Two SNP associations replicated [blocked early in transport 1 homolog (BET1L) rs2280543, RFTS–BioVU meta-odds ratio (OR) = 0.67 95 % confidence interval (CI) 0.38–0.96, Q = 0.70, I = 0, p = 6.9 × 10^?3; trinucleotide repeat containing 6B (TNRC6B) rs12484776, RFTS–BioVU meta-OR = 1.21, 95 % CI 1.07–1.35, Q = 0.24, I = 28.37, p = 8.7 × 10^?3). Meta-analyses combining evidence from RFTS, BioVU, and prior GWAS showed little heterogeneity in effect sizes across studies, with meta-p values between 7.45 × 10^?8 and 3.89 × 10^?9, which were stronger than prior GWAS and supported associations observed for all previously identified loci. These data suggest common variants increase risk for UF in both EA and Japanese populations. However, further research is needed to assess the role of these genes across other racial groups.     

Genome-wide association study of aspirin-exacerbated respiratory disease in a Korean population

Aspirin-exacerbated respiratory disease (AERD) is a nonallergic clinical syndrome characterized by a severe decline in forced expiratory volume in one second (FEV1) following the ingestion of non-steroidal anti-inflammatory drugs (NSAIDs) such as aspirin. The effects of genetic variants have not fully explained all of the observed individual differences to an aspirin challenge despite previous attempts to identify AERD-related genes. In the present study, we performed genome-wide association study (GWAS) and targeted association study in Korean asthmatics to identify new genetic factors associated with AERD. A total of 685 asthmatic patients without AERD and 117 subjects with AERD were used for the GWAS of the first stage, and 996 asthmatics without AERD and 142 subjects with AERD were used for a follow-up study. A total of 702 SNPs were genotyped using the GoldenGate assay with the VeraCode microbead. GWAS revealed the top-ranked variants in 3′ regions of the HLA-DPB1 gene. To investigate the detailed genetic effects of an associated region with the risk of AERD, a follow-up targeted association study with the 702 single nucleotide polymorphisms (SNPs) of 14 genes was performed on 802 Korean subjects. In a case–control analysis, HLA-DPB1 rs1042151 (Met105Val) shows the most significant association with the susceptibility of AERD (p = 5.11 × 10^?7; OR = 2.40). Moreover, rs1042151 also shows a gene dose for the percent decline of FEV1 after an aspirin challenge (p = 2.82 × 10^?7). Our findings show that the HLA-DPB1 gene polymorphism may be the most susceptible genetic factor for the risk of AERD in Korean asthmatics and confirm the importance of HLA-DPB1 in the genetic etiology of AERD.     

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.     

The −319C/+49G/CT60G Haplotype of CTLA-4 Gene Confers Susceptibility to Rheumatoid Arthritis in Mexican Population

Several single nucleotide polymorphisms (SNPs) within the CTLA-4 gene and elevated serum levels of soluble CTLA-4 (sCTLA-4) have been associated with autoimmunity including rheumatoid arthritis (RA). In this case–control study, we evaluated the relationship between the ?319C/T (rs5742909) and CT60 G/A (rs3087243) SNPs and sCTLA-4 levels in 200 RA patients and 200 control subjects (CS) from Western Mexico. Both SNPs were genotyped with the polymerase chain reaction–restriction fragment length polymorphism technique and the sCTLA-4 levels were quantified using an enzyme-linked immunosorbent assay kit. In addition, we performed a haplotype analysis, including our previous data of the +49A/G (rs231775) SNP. The G/A genotype of the rs3087243 SNP was associated with a decreased risk of RA [odd ratio (OR) 0.61, 95 % confidence interval (CI) 0.38–0.96, p = 0.024]. This protection was also observed in the negative anti-cyclic citrullinated peptide group of RA carriers of the A allele (OR 0.48, 95 % CI 0.22–1.05, p = 0.042). On the contrary, we identified the ?319C/+49G/CT60G haplotype of CTLA-4 gene as a risk factor for RA (OR 1.69, 95 % CI 1.13–2.52, p = 0.01). The sCTLA-4 levels were not associated with RA (p = 0.377), but were correlated with the functional disability of these patients (r = 0.282, p = 0.012). However, in CS the C/T genotype of the rs5742909 SNP, as well as the G/G and G/A genotypes of the rs3087243 SNP were associated with higher sCTLA-4 levels (p < 0.001). In conclusion, our results suggest that the ?319C/+49G/CT60G haplotype of CTLA-4 gene is a genetic marker of susceptibility to RA in Western Mexico, whereas the rs3087243 SNP confers protection against this disease. Moreover, both SNPs showed an effect on the sCTLA-4 production in our control population. However, further studies are required to evaluate the role of sCTLA-4 in RA, as well as the molecular and functional basis of the association between both CTLA-4 gene SNPs and soluble levels of CTLA-4 in CS.     

Use of Allele-Specific FAIRE to Determine Functional Regulatory Polymorphism Using Large-Scale Genotyping Arrays

Following the widespread use of genome-wide association studies (GWAS), focus is turning towards identification of causal variants rather than simply genetic markers of diseases and traits. As a step towards a high-throughput method to identify genome-wide, non-coding, functional regulatory variants, we describe the technique of allele-specific FAIRE, utilising large-scale genotyping technology (FAIRE-gen) to determine allelic effects on chromatin accessibility and regulatory potential. FAIRE-gen was explored using lymphoblastoid cells and the 50,000 SNP Illumina CVD BeadChip. The technique identified an allele-specific regulatory polymorphism within NR1H3 (coding for LXR-α), rs7120118, coinciding with a previously GWAS-identified SNP for HDL-C levels. This finding was confirmed using FAIRE-gen with the 200,000 SNP Illumina Metabochip and verified with the established method of TaqMan allelic discrimination. Examination of this SNP in two prospective Caucasian cohorts comprising 15,000 individuals confirmed the association with HDL-C levels (combined beta = 0.016; p = 0.0006), and analysis of gene expression identified an allelic association with LXR-α expression in heart tissue. Using increasingly comprehensive genotyping chips and distinct tissues for examination, FAIRE-gen has the potential to aid the identification of many causal SNPs associated with disease from GWAS.     

Multivariate eQTL mapping uncovers functional variation on the X-chromosome associated with complex disease traits

Very few studies have investigated the associations between genetic polymorphisms and gene expression on the X-chromosome. This is a major bottleneck when conducting functional follow-up studies of trait-associated variants, as those identified in genome-wide association studies (GWAS). We used a multivariate approach to test the association between individual single nucleotide polymorphisms (SNPs) and exon expression levels measured in 356 Epstein–Barr virus-transformed lymphoblastoid cell lines (LCLs) from the Geuvadis RNA sequencing project to identify SNPs associated with variation in gene expression on the X-chromosome, which we refer to as eSNPs. At an FDR of 5 %, we discovered 548 independent [linkage disequilibrium (LD) r ² < 0.1] eSNPs on the X-chromosome. Of these, 35 were in LD (r ² > 0.8) with previously published disease- or trait-associated variants identified through GWAS. One of the strongest eSNPs identified was rs35975601, which was associated with F8A1 expression (p value = 3 × 10^?20) and was in LD with a type 1 diabetes risk variant. Additionally, we identified a number of genes for which eSNPs were in LD with multiple diseases or traits, including DNASE1L1 which was mapped to bilirubin levels, type 1 diabetes and schizophrenia. Our results also indicate that multivariate exon-level analysis provides a more powerful approach than univariate gene-level analysis, particularly when SNPs influence the expression of different exons with different magnitude and/or direction of effect. The associations identified in our study may provide new insights into the molecular process by which gene expression may contribute to trait variation or disease risk in humans.