|Home | About | Journals | Submit | Contact Us | Français|
Juvenile rheumatoid arthritis (JRA) is mediated by Th1-immune responses. In children with JRA, synovial T-cells express high levels of the Th1-chemokine receptor CCR5, which has been implicated in susceptibility to rheumatoid arthritis. To test the hypothesis that genetic variation in CCR5 is associated with susceptibility to JRA, we analyzed patterns of variation in the 5'cis-regulatory region of CCR5 in 124 multiplex families from a JRA affected-sibpair registry. After sequencing the upstream region of CCR5, variants were tested for association with JRA by transmission disequilibrium testing. A single nucleotide polymorphism, C-1835T, was significantly under-transmitted to children with early-onset JRA (p<0.01). C-1835T was genotyped in 424 additional simplex and multiplex families. CCR5-1835T allele was under-transmitted in the cohort of all probands with JRA (p<0.02), as well as in those with early-onset (p<0.01) or pauciarticular JRA (p<0.05). Another variant, a 32-bp deletion in the open reading frame of CCR5, (CCR5-Δ32), was also tested in ~700 simplex and multiplex families. CCR5-Δ32 was also significantly under-transmitted to probands with early-onset JRA (p<0.05). Both variants are in regions under natural selection, and result in functional consequences. Our results suggest these CCR5 variants are protective against early-onset JRA.
Juvenile rheumatoid arthritis (JRA) is a heterogeneous collection of chronic arthritides in children.1 The subtypes of JRA are similar to the major subtypes of Juvenile Idiopathic Arthritis (JIA) described by the International League of Associations for Rheumatology.2 Both genetic and environmental factors are believed to play a role in the pathogenesis of JRA. Numerous associations have been reported between JRA and polymorphisms in genes that encode the human leukocyte antigens (HLA). However, HLA associations account only in part for the overall genetic susceptibility to JRA3, suggesting that other, non-HLA variants also influence susceptibility to JRA.4 The search for non-HLA variants that influence susceptibility to JRA has been challenging for several reasons, including the large number of candidate genes and access to relatively small JRA cohorts that has limited statistical power to detect associations. These challenges can be overcome, in part, by performing adequately powered studies of candidate genes that play a compelling role in the pathophysiology of JRA.
Direct and indirect evidence suggests that the pathogenesis of JRA is mediated, in part, by T-cells. The inflamed synovial tissue in JRA is characterized by a prominent infiltrate of T-cells, plasma cells, macrophages and proliferating synoviocytes.5,6 The recruitment of pro-inflammatory cells into the synovium of children with JRA is mediated by chemokines that selectively attract Th1 T-cells.7–9 These T-cells are characterized by the production of interleukin-2, interferon-γ, and tumor necrosis factor-β. Several studies have demonstrated that Th1 cytokines also predominate in the synovial fluid samples from patients with JRA.10–12
A complex network of chemokines and chemokine receptors orchestrates the homing of T-cells to sites of inflammation. CC chemokine receptor 5 (CCR5) is a receptor preferentially expressed on the surface of Th1-cells. The predominance of CCR5-positive synovial mononuclear cells among patients with different types of inflammatory arthritis suggests that CCR5 plays an important role in synovial inflammation.13 Expression of CCR5 is also increased in synovial T-cells from children with JRA.7 Furthermore, two of the ligands for CCR5, CCL3 and CCL5, have been shown to be elevated in synovial fluid samples from adults with rheumatoid arthritis (RA), as well as in children with JIA.14 CCR5 blockade in collagen induced arthritis, an animal model of RA results in the inhibition of arthritis.15, 16 Together, these data suggest that genetic variation in CCR5 could affect the phenotype of inflammatory arthritis including JRA.
A 32-bp deletion in the open reading frame (exon 3) of the gene encoding CCR5 (CCR5-Δ32), has been tested for association with RA in adults with conflicting results.17–20 Homozygosity for CCR5-Δ32 results in absent expression of CCR5 on the cell surface. This allele is associated with protection against RA in adults.21 However, a systematic investigation of this and other CCR5 polymorphisms in children with JRA has not been reported to our knowledge. The objective of the present study was to test the hypothesis that genetic variants at the CCR5 locus are associated with susceptibility to JRA or subtypes of JRA.
The complete sequence of the 5’ cis-regulatory region of CCR5 was determined in 507 individuals from 124 multiplex JRA families. A total of 10 single nucleotide polymorphisms (SNPs) were identified (Table 1), most of which have been described previously.22 The frequencies of these alleles were comparable to published frequencies in individuals of European ancestry. Two SNPS (G-2248A, G-1939A) were seen in one family each, and one SNP (C-2132T) was observed in two families. All SNPs were in Hardy Weinberg Equilibrium (HWE).
The allele frequency of CCR5-Δ32 among unrelated, unaffected parents was ~10 %, similar to figures published for individuals of European ancestry23. The genotype frequencies of the wild type CCR5 and CCR5-Δ32 in unaffected, unrelated individuals did not deviate from HWE. The allele and genotype frequencies among unrelated, unaffected parents, siblings and probands are shown in Table 2. None of more than 2000 individuals from ~700 families that were typed for both CCR5 −1835 and CCR5-Δ32 were homozygous for the minor alleles of both loci.
Examination of pairwise linkage disequilibrium (LD) between the CCR5 variants revealed that although this region comprised a single haplotype block using the definition of Gabriel et al,24 several pairs of polymorphisms showed low levels of LD using r2, another measure of LD (Figure 1). A total of 9 haplotypes were observed, of which four were rare, being detected in < 1% of the samples. CCR5 -1835T and CCR5-Δ32 were not in LD (r2 = 0.012), and these variants were on two different haplotypes, each with a frequency of ~10%. (Figure 1).
Analysis of the 124 multiplex families revealed no statistically significant associations with any of the SNPs or haplotypes and the JRA cohort as a whole. However after stratification of the JRA cohort into the sub-phenotypes, one allele, CCR5 -1835T was significantly under-transmitted to probands with early-onset JRA (p < 0.01) by transmission disequilibrium test (TDT). A CCR5-haplotype that contained this SNP was also was under-transmitted to children with disease onset before 6 years of age, but this result was not statistically significant. CCR5 C-1835T was genotyped in 23 additional multiplex families and 401 simplex trios with JRA. Examination of the combined cohort confirmed that the −1835T allele was significantly under-transmitted in the cohort of all patients with JRA (p < 0.02), as well as patients with pauciarticular JRA (p < 0.05) and early-onset JRA (p < 0.01) (Table 3). The family-based estimate of the odds ratio for the level of protection conferred by the transmission of −1835T was 0.5 (95 % confidence interval 0.3–0.8, p <0.001) in 254 early-onset JRA trios where there were 37 transmissions (34 %) of −1835T, while the −1835 C allele was transmitted 52 % of the time. Because 21 of the simplex families were ascertained in Utah, whereas the other families were part of a Cincinnati Cohort, we repeated the analyses excluding the Utah cases. The result of the TDT was virtually the same and the frequency of the minor allele among unrelated, unaffected parents was not significantly different between the Cincinnati and Utah simplex cohorts.
Examination of the combined simplex and multiplex cohort of ~700 families by TDT revealed no significant deviation of transmission of CCR5-Δ32 in the combined cohort of probands with JRA (Table 4). However, after stratification, CCR5-Δ32 was significantly under-transmitted to probands with early onset JRA (p< 0.05). The family-based estimate of the odds ratio for the level of protection conferred by the transmission of CCR5-Δ32 was 0.8 (95 % confidence interval 0.5–1.2, p < 0.2) in 330 early-onset JRA trios where there were 57 transmissions (44.5%) of CCR5Δ32, while the CCR5 WT was transmitted 51 % of the time. The lack of statistical significance of the odds ratio might reflect lack of power since only ~70 % of the cohort had both parents’ genotypes available.
In addition, we applied TDT only to the unaffected siblings of the probands with JRA in order to detect if there was a possible transmission distortion of CCR5 −1835T and CCR5-Δ32 in these families. An allele-specific segregation bias, which is defined as an altered transmission of an allele independent of its role in disease, could falsely inflate the statistical significance. Hence, we tested if these alleles could also be preferentially transmitted to unaffected siblings of the JRA probands. There was no significant under or over-transmission of either the CCR5 -1835T or CCR5-Δ32 to the unaffected siblings in these families (Tables 3 and and4).4). Thus, these data provide evidence that our findings are not a result of segregation bias.
We also performed TDT analysis for the transmission of the CCR5 C-1835T and CCR5-Δ32 haplotypes (Table 5). The haplotype containing the rare alleles of both loci was not observed. However, the haplotype containing the major alleles at both loci (−1835C and CCR5 wild type in the open reading frame), was significantly over-transmitted to all probands with JRA (p < 0.03), and to those with early onset (<0.003), or pauciarticular subtypes (< 0.02). TDT analysis of the unaffected siblings showed no significant deviation of transmission of this haplotype. The results of the two other haplotypes containing one minor allele at each locus were similar to the results of the TDT analysis for the individual polymorphisms, and generally trended in the direction of under-transmission of the minor alleles to the probands with JRA, early onset or pauciarticular subtypes.
These results demonstrate that two variants (CCR5-1835T and CCR5-Δ32) in the gene encoding CCR5 are associated with JRA, especially in children with disease onset before the age of 6 years. Both variants are believed to have functional consequences, and it should be noted that both effects on JRA were in the same direction, i.e., protective against JRA.
Case-control association studies are commonly used to investigate associations between phenotypes of interest and genetic variants. However, case-control studies can produce spurious associations due to population stratification. Several approaches have been suggested to account for population stratification.25 Many of these approaches use family-based controls to reduce the effect of population stratification. The TDT provides a joint test of linkage and association and eliminates the effects of stratification when applied to probands and parents. Although TDT based methods are more robust to population stratification, they are thought to be less powerful than conventional case-control studies at detecting associations.25 While failure to detect statistically significant associations by TDT could be due to type 2 statistical errors, positive results obtained by TDT are meaningful. We think that the demonstration of significant deviation of transmission of two CCR5 variants to children with JRA likely represents a true association. This conclusion is further supported by the observation that there is no deviation of transmission of these alleles to unaffected children in these families, although ideally a larger cohort of unaffected siblings would have likely given more robust results.
To minimize statistical artifacts, our strategy was to first test the 5’ cis-regulatory CCR5 variants in only the multiplex cohort, and then to confirm the significant results by genotyping a large cohort of simplex families. Thus, the only SNP we tested in the entire cohort was C-1835T, and we demonstrated a significant association not only with the entire JRA cohort, but also some JRA subtypes, which had been designated a priori. The other polymorphism tested, CCR5-Δ32 was chosen based on its known functional consequence and confirmed associations in adults with RA. Again we had designated the subtypes to be tested a priori. While the JRA subtypes have unique clinical and genetic features, they cannot be thought of as being entirely independent. Similarly variants in LD with one another are also not independent. For these reasons, we did not think it appropriate to apply the conventional Bonferroni correction to our results. Moreover, the Bonferroni correction tends to be very conservative, increasing the possibility of a type-2 statistical error. Our combined JRA cohort is one of the largest reported among studies looking at genetic associations in JRA, likely further minimizing type-2 statistical errors. Our study had sufficient power (> 0.8) to detect an association under most models for a marker allele frequency of 0.1, both for JRA as a whole, and for most sub-phenotypes examined.26 Ultimately, like other genetic association studies the present observations need to be validated by replication in independent cohorts.
The SNP CCR5 C-1835T, located in the 5’ cis-regulatory region of CCR5 has been shown to result in loss of novel nuclear factor binding.27 Investigations of the transcriptional activity of different haplotypes in the cis-regulatory region of CCR5 have demonstrated that CCR5-haplotypes containing -1835T (HHF) have significantly different promoter efficiency compared to other haplotypes.27 Haplotypes containing −1835T are associated with different disease modifying effects among patients with acquired immune deficiency syndrome.28 These data suggest the possibility that −1835T is protective against developing early onset JRA, conceivably by altering the expression level of CCR5. Although one would expect these haplotypes to decrease the expression of CCR5, haplotypes containing −1835T appear to increase the efficiency of the CCR5 promoter in vitro. It is possible that −1835T may have other effects in vivo. Also, this effect of −1835T is on gene expression, and the effect if any, of −1835T on cell surface expression of CCR5 is unknown. Furthermore the effects observed could be due to complex interactions with other polymorphisms on these haplotypes. Ideally, experiments addressing the genotype-phenotype correlations of CCR5 variants and CCR5 expression in subjects with JRA would better address the role played by this polymorphism and others in the development or phenotype of inflammatory arthritis. It is also possible that another variant in LD with this SNP influences the protective effect observed here. For instance −1835T has been reported to be in LD with a SNP in the coding region of another chemokine receptor gene, CCR2 (CCR2-V64I), located ~9kb 5’ of the CCR5 gene on chromosome 3p21.31.22 The CCR2-64I allele has also been shown to have protective influences against human immunodeficiency virus (HIV) disease progression.29
It has been proposed that genetic variants subjected to natural selection are likely to harbor functional polymorphisms and hence are good candidates for association studies.30 The 5’ cis-regulatory region of CCR5 has been demonstrated to be under balancing selection, where diverse haplotypes have been maintained.22 The 5’cis-regulatory region of CCR5 has several polymorphisms that regulate the expression of CCR5, and different CCR5 SNPs and haplotypes have been associated with phenotypic differences in HIV infection.28 It has also been suggested that the CCR5-Δ32 allele has been subjected to recent positive selection.31, 32 (But see 33). The fact that our associations have been discovered in regions affected by natural selection strengthens the evidence that they might have functional significance.
Studies investigating the association between CCR5-Δ32 and RA in adults have reported conflicting results.17, 20 However, a meta-analysis of the published association studies between CCR5-Δ32 and RA suggests that some of these studies were underpowered, and confirms that this polymorphism is strongly and negatively associated with RA (p < 0.0001).21 This suggests that CCR5-Δ32 is associated with protection against the development of RA. Until now, this polymorphism has not been formally investigated in association studies of individuals with JRA, although two children with JRA who were homozygous for this mutation have been described.7, 9 The demonstration of under-transmission of CCR5-Δ32 to probands with early onset JRA is consistent with these observations and suggests that JRA and RA share some genetic risk factors.
To answer the question whether these association results could be due to LD between CCR5 - 1835T and CCR5-Δ32, we examined the haplotypes and the LD pattern in CCR5. The minor alleles, −1835T and CCR5-Δ32 were not in LD, and they occur on two different haplotypes. Examination of the published data on Great Apes, and Old and New world monkeys, shows that neither allele is present in these species and that both are newly derived in humans, suggesting that both of these variants are derived from the ancestral haplotype of CCR5-1835C-CCR5 WT.27 Thus, the protective effects observed at these loci appear to be independent of each other and not due to LD. As anticipated, we did not find haplotypes that had the minor alleles at both loci. However, the haplotype containing CCR5 −1835C and CCR5 WT (i.e., lacking the protective variants) did show evidence of excessive transmission to all children with JRA and those with early onset JRA by TDT, further strengthening our results. Together these observations support the notion that CCR5 might play important roles in the pathogenesis of inflammatory arthritis, suggesting that CCR5 blockade could modify the phenotype of inflammatory arthritis.
DNA samples were obtained from children with JRA and one or both parents. DNA from unaffected siblings was also available in several families. The diagnosis of JRA was made according to the American College of Rheumatology (ACR) criteria.1 As most subjects had been enrolled before adoption of the more recently proposed criteria for the diagnosis of JIA, the ACR criteria for a diagnosis of JRA was used. Samples from 147 multiplex families with two or more affected siblings with JRA from around the USA, and 540 simplex families with one affected child with JRA from Cincinnati OH were available from a JRA Registry sponsored by the National Institutes of Arthritis and Musculoskeletal and Skin diseases, and the Cincinnati Children’s Hospital Medical Center. The multiplex families were ascertained and enrolled from around the USA through a JRA Registry as described earlier.34 Additionally, 21 simplex JRA families were ascertained from Salt Lake City UT. In all, DNA from both parents was available in 105 of the multiplex families, 370 of the simplex families from Cincinnati, and all 21 families from Utah. DNA from at least one parent was available in most of the remaining families. DNA from unaffected siblings was available in 38 multiplex and 313 Cincinnati simplex families. Institutional review board approvals were obtained from participating institutions, and parents and/or children provided informed consent. In all, there were 419 children with pauciarticular onset (< 5 joints in the first 6 months of disease), 331 with polyarticular onset (5 or more joints in the first 6 months of disease), and 107 with systemic onset JRA. The median age of onset of JRA in the simplex cohort was 5.8 years, and the median age on onset in the multiplex cohort was 4.5 years. The children with JRA were predominantly female, similar to other studies, with 72% and 73 % of the cases from simplex and multiplex cohorts, respectively, being female. Most (85% of the simplex and 88% of the multiplex) families were of European Ancestry.
A 1122 bp region corresponding to human CCR5 −2867 to −1745 was amplified by polymerase chain reaction (PCR) in 507 individuals from 124 multiplex families, and sequenced on both strands as previously described.22 Briefly, cycle sequencing was carried out in 5μl reaction volumes using ABI BigDye (v3.1) Terminator chemistry. The sequence ladders were precipitated with 62.5% EtOH/1M KOAc, resuspended in 15 μl of dH2O and electrophoresed on an ABI3700 DNA analyzer prepared with POP-5 capillary gel matrix. The sequenced region includes part of exon 1, intron 1, exons 2A and 2B, and the 5’ end of intron 2. All of these exons are non-coding. This region, defined in earlier publications as the 5’ cis regulatory region, contains several unique CCR5-SNPs and haplotypes that have been shown to have functional consequences.22, 27 The primers used for PCR amplification were CCR5L2 (5’ - CCA AAC TGT GAC CCT TTC C – 3’) and CCR5R3: (5’ - AGT AGC TCT CTG CTG TCT TCT CA – 3’). The programs PHRED, PHRAP and CONSED were used to evaluate sequence trace files35. Potential heterozygotes were identified by using the program POLYPHRED version 3.5. Polymorphisms were verified by examining the individual sequence trace files manually, and for most polymorphisms both the forward and reverse sequences were evaluated. CCR5 re-sequencing allowed us to successfully identify ~ 99% of the polymorphisms in this region.
DNA samples from 23 multiplex families and 401 simplex families were genotyped for the SNP C-1835T, using a PCR restriction fragment length polymorphism based assay using EcoRV as described previously.28 This SNP results in a C to T change at position −1835 in the 5’ cis-regulatory region of CCR5, and has been shown to result in loss of novel nuclear factor binding.27 The restriction site for EcoRV was created by changing an A to G in the antisense primer at position −1832. The enzyme digests the amplicons that contain the T allele at position − 1835. The following primers were used for PCR amplification: CCR5927_S: (5’ GTT GGT TTA AGT TGG CTT 3’) and CCR5927_AS: (5’ ATC TTA AAG ATT ATA TTT TAA GAT AAT TGT ATG AGC ACT TGG TGT TTG CCA GAT -3’). The genotyping success rate was 95%. Several representative samples that had been sequenced were also genotyped by the RFLP method for quality control, and the genotyping calls were consistent.
Genotyping for the CCR5-Δ32 polymorphism was carried out by performing a PCR reaction in 697 multiplex and simplex families with JRA. PCR was performed with the following primers spanning the region of the CCR5-Δ32 deletion: CCR5_d32_R: 5’ TTC AGG AGA AGG ACA ATG TTG TAG G 3’ and CCR5_d32_L: 5’ ATT ACA CCT GCA GCT CTC ATT TTC 3’. The observed wild type and deletion fragments were 250 bp and 218 bp respectively. The genotyping success rate was 99%.
Pairwise LD between the different SNPs in the promoter region and the CCR5-Δ32 polymorphism was assessed using the program Haploview, which also was used to infer haplotypes and haplotype blocks.36 Only variants with frequency > 5 % were used for inferring haplotypes and haplotype blocks. The haplotype block definition proposed by Gabriel et al was used, wherein pairs of SNPs in strong LD are defined as those in which the one-sided upper 95% confidence bound on D′ (a normalized measure of allelic variation) is >98% and the lower bound is >70%.24 LD was also assessed using an alternate measure, r2, which unlike D' is not strongly biased by low frequency variants or missing haplotypes.
All SNPs were tested for HWE using unrelated, unaffected individuals (i.e., parents). SNPs and haplotypes were tested for association with JRA by TDT, which tests for preferential transmission of specified alleles from heterozygous parents to the affected offspring. TDT analysis was performed using TRANSMIT 126.96.36.199 TRANSMIT infers missing parental haplotypes using HWE, thus allowing the inclusion of all families in the analysis.
The multiplex JRA cohort was analyzed for an association with the CCR5 SNPs as well as the haplotypes. The combined JRA cohort was analyzed for an association with C-1835T by TDT. Additionally, 697 simplex and multiplex families were analyzed for an association with CCR5- Δ32 polymorphism by TDT. All analyses were repeated after stratifying the JRA cohort by onset type (pauciarticular, polyarticular, or systemic), as well as by age of onset (early onset < 6 years, or late onset > 6 years). The division into early/late onset JRA included all subtypes. These JRA sub-phenotypes were designated a priori based on reports suggesting that JRA comprises a collection of clinically and genetically heterogeneous phenotypes38, 39, and by the demonstration of age-specific effects of JRA associated HLA alleles.40 Furthermore the various subtypes of juvenile arthritis differ in their outcomes.41 The approach of stratifying the JRA cohorts in the present study was similar to the stratification scheme employed in the genome-wide scan for JRA, which identified several regions of suggestive linkage on stratification of the JRA cohort into sub-phenotypes.42 These observations suggest that the various subtypes of JRA likely represent distinct phenotypes that share the feature of inflammatory arthritis in children. In addition TDT was also performed using only the unaffected siblings to determine if there was preferential transmission of CCR5 alleles or haplotypes to unaffected siblings as well in these families.
We also calculated the family-based estimate of the odds ratio for the level of protection conferred by CCR5 −1835T and by CCR5-D32 in families of probands with early onset JRA using the methods described by Ackerman et al, where transmitted versus non-transmitted alleles from the parental matings are compared.43 They proposed that in the absence of segregation distortion, and in families in HWE, the non-transmitted chromosomes can be used as population controls to estimate relative risk. Both of these variants were in HWE in our cohort, and there was no significant deviation of transmission of these variants to unaffected siblings, suggesting a lack of segregation distortion. Only trios where both parental genotypes where available were used for these calculations. Thus, 254 and 330 trios with early onset JRA were used for the analyses of C-1835T and CCR5-Δ32 respectively.
Supported in part by grants from the National Institute of Arthritis and Musculoskeletal and Skin Diseases (K23-AR-50177, N01-AR-42272, P01-AR-048929, P30-AR-47363), National Center for Research Resources (M01-RR-00064); The Cincinnati Pediatric Rheumatology Tissue Repository, and The Cincinnati Children’s Hospital Medical Center, Cincinnati, OH; The Val A Browning Charitable Foundation, The Primary Children’s Medical Center Foundation, The Clinical Genetics Research Program, and The Children’s Health Research Center Salt Lake City, UT. We would also like to thank the families that participated in this project.