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.     

Exonic Re-Sequencing of the Chromosome 2q24.3 Parkinson’s Disease Locus

Genome-wide association studies (GWAS) in Parkinson’s disease (PD) have identified over 20 genomic regions associated with disease risk. Many of these loci include several candidate genes making it difficult to pinpoint the causal gene. The locus on chromosome 2q24.3 encompasses three genes: B3GALT1, STK39, and CERS6. In order to identify if the causal variants are simple missense changes, we sequenced all 31 exons of these three genes in 187 patients with PD. We identified 13 exonic variants including four non-synonymous and three insertion/deletion variants (indels). These non-synonymous variants and rs2102808, the GWAS tag SNP, were genotyped in three independent series consisting of a total of 1976 patients and 1596 controls. Our results show that the seven identified 2q24.3 coding variants are not independently responsible for the GWAS association signal at the locus; however, there is a haplotype, which contains both rs2102808 and a STK39 exon 1 6bp indel variant, that is significantly associated with PD risk (Odds Ratio [OR] = 1.35, 95% CI: 1.11–1.64, P = 0.003). This haplotype is more associated than each of the two variants independently (OR = 1.23, P = 0.005 and 1.10, P = 0.10, respectively). Our findings suggest that the risk variant is likely located in a non-coding region. Additional sequencing of the locus including promoter and regulatory regions will be needed to pinpoint the association at this locus that leads to an increased risk to PD.     

Influence of the CCR2-V64I Polymorphism on Human Immunodeficiency Virus Type 1 Coreceptor Activity and on Chemokine Receptor Function of CCR2b, CCR3, CCR5, and CXCR4

The chemokine receptors CCR5 and CXCR4 are used by human immunodeficiency virus type 1 (HIV-1) in conjunction with CD4 to infect cells. In addition, some virus strains can use alternative chemokine receptors, including CCR2b and CCR3, for infection. A polymorphism in CCR2 (CCR2-V64I) is associated with a 2- to 4-year delay in the progression to AIDS. To investigate the mechanism of this protective effect, we studied the expression of CCR2b and CCR2b-V64I, their chemokine and HIV-1 coreceptor activities, and their effects on the expression and receptor activities of the major HIV-1 coreceptors. CCR2b and CCR2b-V64I were expressed at similar levels, and neither molecule affected the expression or coreceptor activity of CCR3, CCR5, or CXCR4 in cotransfected cell lines. Peripheral blood mononuclear cells (PBMCs) from CCR2-V64I heterozygotes had normal levels of CCR2b and CCR5 but slightly reduced levels of CXCR4. CCR2b and CCR2b-V64I functioned equally well as HIV-1 coreceptors, and CCR2-V64I PBMCs were permissive for HIV-1 infection regardless of viral tropism. The MCP-1-induced calcium mobilization mediated by CCR2b signaling was unaffected by the polymorphism, but MCP-1 signaling mediated by either CCR2b- or CCR2-V64I-encoded receptors resulted in heterologous desensitization (i.e., limiting the signal response of other receptors) of both CCR5 and CXCR4. The heterologous desensitization of CCR5 and CXCR4 signaling by both CCR2 allele receptor types provides a mechanistic link that might help explain the in vivo effects of CCR2 gene variants on progression to AIDS as well as the reported antiviral activity of natural CCR2 ligands.     

Association study of CCR6 rs3093024 with Rheumatoid Arthritis in a Pakistani cohort

C-C Chemokine receptor 6 (CCR6), an important protein in inflammatory and immunological responses, has been previously reported to be associated with rheumatoid arthritis (RA). Therefore, in order to replicate these findings, a case-control study was conducted on 500 subjects (including 250 RA patients and 250 healthy controls) of Pakistani origin. The aim of this study was to determine the association of CCR6 rs3093024 variant with RA and identify its role in splicing events using bioinformatics tools. The clinical and demographic characteristics of the patients were collected using a well-designed questionnaire. The genotype frequencies of CCR6 rs3093024 variant were determined using tetra-primer ARMS-PCR (amplification of refractory mutation system-polymerase chain reaction) method. A significant difference was found between CCR6 rs3093024 genotype frequencies [P = 0.0016, χ² = 12.915]. Similarly, a significant difference in the allele frequencies between RA patients and healthy controls was also observed [P = 0.0003 and OR (95% CI) = 0.63 (0.49–0.80)]. The stratification of patients showed that there was a significant increase in AA genotype against AG + GG in patients [P = 0.0014, OR (95% CI) = 2.0 (1.32–3.02)]. Furthermore, using bioinformatics analysis, it was found that CCR6 rs3093024 variant might create a potential splicing enhancer motif (SF2/ASF (IgM-BRCA1) with score of 77.92; Threshold 70.53), which might have important impact on the product of this gene. This study suggests that the A variant of CCR6 rs3093024 variant is significantly associated with RA-risk and its G variant is protective in Pakistani population but a multi-cohort large sized population study is needed to elucidate these results. Moreover, functional studies are needed to highlight the effects of this polymorphism on the function of CCR6 gene.     

Increased frequency of CCR4+ and CCR6+ memory T-cells including CCR7+CD45RAmed very early memory cells in granulomatosis with polyangiitis (Wegener's)

Introduction

Chemokine receptors play an important role in mediating the recruitment of T cells to inflammatory sites. Previously, small proportions of circulating Th1-type CCR5+ and Th2-type CCR3+ cells have been shown in granulomatosis with polyangiitis (GPA). Wondering to what extent CCR4 and CCR6 expression could also be implicated in T cell recruitment to inflamed sites in GPA, we investigated the expression of CCR4 and CCR6 on T cells and its association with T cell diversity and polarization.

Methods

Multicolor flow cytometry was used to analyze CCR4, CCR6, and intracellular cytokine expression of T cells from whole blood of GPA-patients (n = 26) and healthy controls (n = 20). CCR7 and CD45RA were included for phenotypic characterization.

Results

We found a significant increase in the percentages of circulating CCR4+ and CCR6+ cells within the total CD4+ T cell population in GPA. In contrast, there was no difference in the percentages of CD8+CCR4+ and CD8+CCR6+ T cells between GPA and healthy controls. CCR4 and CCR6 expression was largely confined to central (T_CM) and effector memory T cells (T_EM, T_EMRA). A significant increase in the frequency of CCR4+ and CCR6+ T_EMRA and CCR6+ T_CM was shown in GPA. Of note, we could dissect CCR4 and CCR6 expressing CCR7+CD45RA^med very early memory T cells (T_VEM) from genuine CCR7+CD45RA^high naïve T cells lacking CCR4 and CCR6 expression for peripheral tissue-migration within the CCR7+CD45RA+ compartment. The frequencies of CCR4+ and CCR6+ T_VEM were also significantly increased in GPA. An increased percentage of IL-17+ and IL-22+ cells was detected in the CCR6+ cell subsets and IL-4+ cells in the CRR4+ cell subset when compared with CD4+ cells lacking CCR4 and CCR6 expression.

Conclusions

Increased frequencies of circulating CCR4+ and CCR6+ memory T cell subsets including hitherto unreported T_VEM suggest persistent T cell activation with the accumulation of CCR4+ and CCR6+ cells in GPA. CCR4 and CCR6 could be involved in the recruitment of T cells including cytokine-producing subsets to inflamed sites in GPA.     

CCR6, the Sole Receptor for the Chemokine CCL20, Promotes Spontaneous Intestinal Tumorigenesis

Interactions between the inflammatory chemokine CCL20 and its receptor CCR6 have been associated with colorectal cancer growth and metastasis, however, a causal role for CCL20 signaling through CCR6 in promoting intestinal carcinogenesis has not been demonstrated in vivo. In this study, we aimed to determine the role of CCL20-CCR6 interactions in spontaneous intestinal tumorigenesis. CCR6-deficient mice were crossed with mice heterozygous for a mutation in the adenomatous polyposis coli (APC) gene (APC^MIN/+ mice) to generate APC^MIN/+ mice with CCR6 knocked out (CCR6KO-APC^MIN/+ mice). CCR6KO-APC^MIN/+ mice had diminished spontaneous intestinal tumorigenesis. CCR6KO-APC^MIN/+ also had normal sized spleens as compared to the enlarged spleens found in APC^MIN/+ mice. Decreased macrophage infiltration into intestinal adenomas and non-tumor epithelium was observed in CCR6KO-APC^MIN/+ as compared to APC^MIN/+ mice. CCL20 signaling through CCR6 caused increased production of CCL20 by colorectal cancer cell lines. Furthermore, CCL20 had a direct mitogenic effect on colorectal cancer cells. Thus, interactions between CCL20 and CCR6 promote intestinal carcinogenesis. Our results suggest that the intestinal tumorigenesis driven by CCL20-CCR6 interactions may be driven by macrophage recruitment into the intestine as well as proliferation of neoplastic epithelial cells. This interaction could be targeted for the treatment or prevention of malignancy.     

Phenotype Refinement Strengthens the Association of AHR and CYP1A1 Genotype with Caffeine Consumption

Two genetic loci, one in the cytochrome P450 1A1 (CYP1A1) and 1A2 (CYP1A2) gene region (rs2472297) and one near the aryl-hydrocarbon receptor (AHR) gene (rs6968865), have been associated with habitual caffeine consumption. We sought to establish whether a more refined and comprehensive assessment of caffeine consumption would provide stronger evidence of association, and whether a combined allelic score comprising these two variants would further strengthen the association. We used data from between 4,460 and 7,520 women in the Avon Longitudinal Study of Parents and Children, a longitudinal birth cohort based in the United Kingdom. Self-report data on coffee, tea and cola consumption (including consumption of decaffeinated drinks) were available at multiple time points. Both genotypes were individually associated with total caffeine consumption, and with coffee and tea consumption. There was no association with cola consumption, possibly due to low levels of consumption in this sample. There was also no association with measures of decaffeinated drink consumption, indicating that the observed association is most likely mediated via caffeine. The association was strengthened when a combined allelic score was used, accounting for up to 1.28% of phenotypic variance. This was not associated with potential confounders of observational association. A combined allelic score accounts for sufficient phenotypic variance in caffeine consumption that this may be useful in Mendelian randomization studies. Future studies may therefore be able to use this combined allelic score to explore causal effects of habitual caffeine consumption on health outcomes.     

CYP2B6 Non-Coding Variation Associated with Smoking Cessation Is Also Associated with Differences in Allelic Expression,Splicing, and Nicotine Metabolism Independent of Common Amino-Acid Changes

The Cytochrome P450 2B6 (CYP2B6) enzyme makes a small contribution to hepatic nicotine metabolism relative to CYP2A6, but CYP2B6 is the primary enzyme responsible for metabolism of the smoking cessation drug bupropion. Using CYP2A6 genotype as a covariate, we find that a non-coding polymorphism in CYP2B6 previously associated with smoking cessation (rs8109525) is also significantly associated with nicotine metabolism. The association is independent of the well-studied non-synonymous variants rs3211371, rs3745274, and rs2279343 (CYP2B6*5 and *6). Expression studies demonstrate that rs8109525 is also associated with differences in CYP2B6 mRNA expression in liver biopsy samples. Splicing assays demonstrate that specific splice forms of CYP2B6 are associated with haplotypes defined by variants including rs3745274 and rs8109525. These results indicate differences in mRNA expression and splicing as potential molecular mechanisms by which non-coding variation in CYP2B6 may affect enzymatic activity leading to differences in metabolism and smoking cessation.     

Analysis of a genome-wide association study-linked locus (CCR6) in Asian rheumatoid arthritis

A genome-wide association study in Japan identified the C-C chemokine receptor type 6 gene (CCR6) as associated with rheumatoid arthritis (RA). This finding has not been validated in other Asian populations. A case-control study involving 996 subjects, comprising 440 controls and 556 RA patients, was done to determine their anticyclic citrullinated peptide (anti-CCP) antibody status and CCR6 polymorphism (rs3093024) genotype. Three hundred eighty-seven patients were anti-CCP positive and 153 anti-CCP negative. Logistic regression showed that allele A was likely to increase the risk of developing RA among females via a recessive model (odds ratio [OR]=1.55, 95% confidence interval [CI]=1.01, 2.39), whereas the risk effect appeared to be reduced among males via an additive model (OR=0.60, 95% CI=0.42, 0.85). Considering only subjects who are anti-CCP positive, allele A increased RA risk among females via a recessive model (OR=1.68, 95% CI=1.07, 2.64) but decreased the risk among males via an additive model (OR=0.59, 95% CI=0.39, 0.89). We showed that CCR6 polymorphism was a risk factor among females but a protective factor among males. Functional studies are warranted to unravel the pathophysiological relevance of the gene variant and other linked variants with RA.