We’re sorry, something doesn't seem to be working properly.
Please try refreshing the page. If that doesn't work, please contact support so we can address the problem.
An investigation of polymorphisms in the 17q11.2-12 CC chemokine gene cluster for association with multiple sclerosis in Australians
BMC Medical Geneticsvolume 7, Article number: 64 (2006)
Multiple sclerosis (MS) is a disorder of the central nervous system (CNS) characterised by inflammation and neuronal degeneration. It is believed to result from the complex interaction of a number of genes, each with modest effect. Chemokines are vital to the migration of cells to sites of inflammation, including the CNS, and many are implicated in MS pathogenesis. Most of the CC chemokine genes are encoded in a cluster on chromosome 17q11.2-12, which has been identified in a number of genome wide screens as being potentially associated with MS.
We conducted a two-stage analysis to investigate the chemokine gene cluster for association with MS. After sequencing the chemokine genes in several DNA pools to identify common polymorphisms, 12 candidate single-nucleotide polymorphisms (SNPs) were genotyped in a cohort of Australian MS trio families.
Marginally significant (uncorrected) transmission distortion was identified for four of the SNPs after stratification for several factors. We also identified marginally significant (uncorrected) transmission distortion for haplotypes encompassing the CCL2 and CCL11 genes, using two independent cohorts, which was consistent with recent reports from another group.
Our results implicate several chemokines as possibly being associated with MS susceptibility, and given that chemokines and their receptors are suitable targets for therapeutic agents, further investigation is warranted in this region.
Multiple sclerosis (MS) is the most common chronic neurological disease in young adults. It is characterised by inflammation of the central nervous system (CNS), believed to be the result of an autoimmune reaction resulting in demyelination and destruction of neural supporting cells . Epidemiological studies suggest a multifactorial aetiology for MS, implicating a complex interplay between environmental and genetic factors . In the past ten years, a large number of genome wide screens have been conducted, including the recent GAMES collaboration (Genetic Analysis of Multiple sclerosis in EuropeanS) . Multiple regions of potential linkage and association with MS have been identified, suggesting that genetic predisposition to MS might result from the modest contribution of many genetic factors, which, if identified, may present important new therapeutic targets .
The inflammatory response that is characteristic of MS requires the targetted migration of leukocytes into the CNS, which is under the control of chemokines. Over 40 members of the human chemokine family have been identified, which act upon a variety of leukocytes via interactions with almost 20 seven-transmembrane domain chemokine receptors . Chemokines are small molecules of approximately 8–10 kDa in size, and are primarily classified on the basis of the relative positioning of two conserved cysteines. In the CC chemokines, the two cysteines are adjacent, and in the CXC chemokines, a single amino acid residue separates the cysteines. Two smaller subgroups have also been identified; the CX3C chemokines, and the C chemokines.
Functionally, CC chemokines chemoattract a wide range of cells, including lymphocytes, dendritic cells, monocytes and some granulocytes [6, 7], whilst CXC chemokines are chemoattractant for neutrophils and lymphocytes [6, 7]. There is increasing evidence for chemokines possessing abilities beyond that of migration, including T-helper cell subset differentiation , T cell costimulation [9, 10], and macrophage and natural killer cell maturation [11, 12].
Substantial evidence supports the involvement of CC chemokines in the pathogenesis of MS. In the mouse model of MS, experimental autoimmune encephalomyelitis (EAE), knockout of CCL2 leads to resistance to disease induction , whilst CCL3, CCL4 and CCL5 have all been implicated in EAE development [14, 15]. However, CCL3-knockout mice were found to be fully susceptible to myelin oligodendrocyte glycoprotein (MOG)-induced EAE . In MS lesions, expression of CC chemokines, including CCL2, CCL3, CCL4, CCL5, CCL7 and CCL8, and their receptors, have been identified on a wide variety of cells, such as astrocytes, microglia and perivascular T cells [17–22]. In addition, altered levels of CC chemokines and receptors have been identified in the serum and cerebrospinal fluid (CSF) of MS patients; some are elevated (including CCL5), whilst CCL2 is decreased in the CSF [22–26], possibly due to removal by CCR2-positive migrating cells as they cross the blood-brain barrier .
Fourteen of the 28 CC chemokine genes are clustered on chromosome 17q11.2-12 (Figure 1) . This cluster spans slightly less than 2 Mb, and is split into two sub-clusters separated by a gap of 1.5 Mb. The 17q11.2-12 region has been implicated in genome-wide screens for linkage and association with MS [29–34], and in a meta-analysis of three genome screens, the most significant nonparametric linkage score was obtained for this region . The 17q region is also syntenic to an EAE quantitative trait locus on chromosome 10, which includes a chemokine gene cluster , and non-synonymous polymorphisms in murine CCL1, CCL2 and CCL12 were identified as candidates for the eae7 quantitative trait locus .
Candidate gene studies of CC chemokines in MS have not been extensive. Rather, the majority of studies have focussed on other diseases. Polymorphisms from several CC chemokine genes have been variably associated with diseases such as tuberculosis , HIV , arthritis , and asthma and atopy [41, 42]. The CCR5Δ32 mutation has been widely studied in MS. Whilst it was found not to be a general susceptibility factor for MS , it may have effects on age of onset or disease progression [44, 45]. More recently, several moderate single-loci and haplotypic associations with MS were identified for single nucleotide polymorphisms (SNPs) from the CC chemokine gene cluster .
It was the aim of our study to conduct candidate gene analyses for the CC chemokines in order to identify polymorphisms and/or haplotypes associated with MS. We undertook a two-stage analysis. The first stage involved the scanning of the genes for published (online) and novel polymorphisms in several DNA pools by DNA sequencing. From these, 12 candidate SNPs were then individually genotyped in a cohort of MS trio families in the second stage. We also sought confirmation of our findings in an independent cohort of Australian MS families. Marginally significant (uncorrected) transmission distortion was identified for four of the SNPs, as well as for haplotypes encompassing the CCL2 and CCL11 genes.
For the main study, all sporadic MS cases and parents (including 373 MS trio families) were recruited by our facility at the Institute for Immunology and Allergy Research, Westmead Hospital (Australia). The familial DNA pool was derived from MS probands from multicase families obtained from the National Register of Multiple Sclerosis Families (Rex Simmons, Canberra Hospital, Australia). Control individuals were composed of local staff members and spouses of sporadic MS cases. MS patients and controls had the same ethnic composition, and were of similar average age (patients = 50yo; controls = 52.8yo). 90% of MS patients were of northern European origin, while the remainder were of southern European origin. The ratio of female to male patients was 4:1; the control ratio was 1.5:1. Approximately 60% of MS patients were HLA-DRB1*1501 positive; 73% had relapsing-remitting (RR)-MS, 20% secondary-progressive (SP)-MS, 5% primary-progressive (PP)-MS, and 2% progressive relapsing (PR)-MS. An additional cohort of 208 Australian MS trio families was obtained from the Southern MS Genetics Consortium for independent validation of our initial results. All cases were verified as having MS as defined by the Poser criteria , and provided written informed consent.
Pooled DNA sequencing
DNA was extracted from whole blood using a rapid salting out method . Methods for the construction of the DNA pools has been described elsewhere . Four DNA pools were constructed: 217 HLA-DRB1*1501 positive sporadic MS patients; 155 HLA-DRB1*1501 negative sporadic MS patients; 169 MS probands from multicase families; 185 unrelated, unaffected controls. HLA-DRB1*1501 status was genotyped using a SYBR Green assay (Applied Biosystems (ABI)), described elsewhere .
DNA sequencing primers were designed to encompass the exons of all genes analysed, as well as segments of the putative promoter region, and often part of the 3' untranslated region (UTR) (see additional file 1: 17qSuppTable1). Where possible, primers were also designed to maximise coverage of published SNPs, which were identified from online SNP databases.
DNA sequencing of the chemokine genes was performed on each of the DNA pools and a single individual, included for SNP allele frequency correction purposes and for potential identification of rare SNPs. PCR products for each gene segment were amplified using the ABI 2× PCR Master Mix, with 50 ng of genomic DNA and appropriate amplification primer pairs (Sigma Genosys) to a final concentration of 5 ng/μL. PCR reactions followed a standard PCR program: 1 cycle of 95°C for 10 minutes; 35 cycles of 94°C for 30 seconds, 62°C for 30 seconds, and 72°C for 1 minute and 10 seconds; a final elongation step of 72°C for 10 minutes. PCR products were purified from the PCR cocktail by use of either the QIAQuick® PCR purification kit (QIAGEN), or ExoSAP-IT (USB Corporation). Sequencing reactions were performed using the ABI Big Dye® Terminator v3.1 Cycle Sequencing Kit. Final samples were sequenced using an ABI 3100 sequence analyser. Sequences were viewed using the ABI PRISM™ EditView software.
Genotyping in Trio Families
An initial set of seven CC chemokine gene SNPs were selected for individual genotyping in a cohort of 204 MS trio families (Table 1): CCL2 -2581, CCL2 -2138, CCL11 67, CCL5 -471, CCL15 -1284, CCL15 136+88 and CCL23 -289. These SNPs were genotyped using the SNaPshot assay, which allows the genotyping of multiple SNPs simultaneously in a single-tube multiplexed reaction, using primer-extension methodology  (see additional file 3: 17qSuppTable3 for SNaPshot primer details for each of the seven SNPs).
SNaPshot genotyping was performed on the pooled PCR products of each individual, amplified using either the ABI 2× PCR master mix, or the Fermentas 2× PCR master mix (Progen). SNaPshot reactions were optimised in individuals of known genotype. Between 100–200 ng of DNA was used from each individual with the appropriate amplification primer pair (Sigma Genosys) to a final concentration of 5 ng/μL. PCR products were amplified using the standard conditions, as described above. The successful amplification of PCR products was verified on 1% agarose gels. The PCR products for each individual were then pooled, and an aliquot purified using ExoSAP-IT (USB Corporation). To this aliquot was then added 2.5 μL of SNaPshot Multiplex reagent (ABI), and pooled SNaPshot primers as defined in Supplementary Table 3, and the SNaPshot procedure was followed as per the manufacturer's instructions (ABI). Samples were then scanned using an ABI 3100. SNaPshot results were viewed using the GeneScan™ software (ABI) and genotypes determined using the GeneMapper™ software package (ABI).
An additional five CC chemokine gene SNPs were selected for genotyping in the full cohort of 373 MS trio families (Table 1): CCL11 -488, CCL8 205, CCL16 -595, CCL14 -649 and CCL23 316. These SNPs were genotyped using the SNPlex™ genotyping system (ABI), for which SUPAMAC (University of Sydney, Australia) were contracted. This system allows high-throughput genotyping of up to 48 SNPs simultaneously in a single tube using an oligonucleotide ligation assay. Initially, candidate SNP details were submitted, after which suitability for the assay was determined by running the SNPs through an algorithm. Oligonucleotides were then designed for each SNP, and applied to the DNA samples (500 ng genomic DNA). Samples were scanned using an ABI 3730. Genotyping results were determined using the GeneMapper™ software package (ABI), utilising an allelic discrimination cluster analysis to allocate genotype calls with maximum stringency.
Additional genotyping of validation cohort
Genotyping of the additional cohort of MS trio families was contracted to the Australian Genome Research Facility (AGRF; Brisbane, Australia). Only the CCL2 -2138A>T and CCL11 -488C>A SNPs were genotyped using the Sequenom Autoflex Mass Spectrometer. After submission of SNP details, PCR oligonucleotides were designed and applied to the DNA samples.
Analyses & statistics
Sequencing results from DNA pools and individuals were compared and scanned for published SNPs and novel polymorphisms. All SNPs were classified according to the nomenclature recommended by den Dunnen and Antonarakis , and are relative to the translation start site. For genes transcribed in the anti-sense direction (indicated in Table 1 with (-)), the complementary nucleotide for each allele was used, and is used throughout this manuscript. Minor allele frequencies (MAF) of SNPs were estimated by comparing the relative peak heights of the alleles. Where the individual was heterozygous for a particular SNP, the MAF in the DNA pools could be corrected for SNP-specific variation in peak-height intensity. The significance of differences between the MS patient pools and the control pool was measured by estimated relative risk (ERR) calculations.
The HelixTree genetics analysis software (Golden Helix Inc, Bozeman, USA) was used to calculate linkage disequilibrium (LD) for the 12 SNPs analysed and to determine haplotypes in parents only. By using the Expectation/Maximisation (EM) algorithm, the software was able to calculate probabilities of each haplotype occurring, based on multi-locus genotypes. It calculated D' and r2 values, which are accepted measures of LD, and performed a χ2 comparison for each pair of SNPs. HaploBlockFinder  was used to establish the haplotype block structure across the CC chemokine gene cluster.
TRANSMIT  was then used to analyse for haplotype transmission disequilibrium, estimating χ2 values for individual haplotypes, as well as global χ2 values for all haplotypes analysed. Mendelian transmission was checked using MERLIN .
Identification and analysis of SNPs in CC chemokine genes by DNA pool sequencing
The CC chemokine genes were scanned using DNA pool sequencing in order to identify common polymorphisms in the region, including the potential identification of novel SNPs, and to establish whether associations might exist for these SNPs. In all, 50 SNPs were identified across the CC chemokine gene cluster (a list of all variants identified is available in the additional file 2: 17qSuppTable2), 48 of which were known; novel SNPs were detected in CCL8 and CCL15. Five common SNPs within exons were confirmed, four of which cause codon changes. Twenty-nine SNPs had MAF >0.15 in the control pool, designated as 'common' here, and for which we had statistical power to detect associations for ERR >1.5.
The greater than 90% DNA sequence homology between CCL3 and CCL3L1, and between CCL4 and CCL4L1, made it impractical to design specific primers for each of these genes. Thus, these four genes were not sequenced.
Population genetics of individually genotyped markers
Based on their likely functional significance (codon changing, putative promoter region, reported functional significance) and ERR from the pooled DNA sequencing, 12 SNPs were chosen for genotyping in the MS cohort, either by SNaPshot or SNPlex™ (Table 1). The CCL8 205 SNP was originally chosen for analysis by SNPlex™, but failed to pass the SNPlex™ algorithm. Thus, it was replaced with the CCL8 -572 SNP, which is in LD with the original SNP (International HapMap Project ).
Genotype information retrieval for the seven SNPs genotyped using SNaPshot was 100%. Genotype information available for single-locus analysis for the five SNPs genotyped by SNPlex™ ranged between 84–87%; losses due to failed genotyping, and the exclusion of genotypes called with low stringency. The SNPlex™ genotyping, which is based on a highly multiplex PCR, was quite sensitive to DNA quality. Unambiguous genotypes could be determined using the SNPlex™ calling algorithm, which we have used at the highest stringency. The genotype frequencies for all SNPs conformed to Hardy-Weinberg equilibrium estimates, except for the CCL11 -488 SNP in mothers. The Mendelian error rate for the seven SNPs was ≤ 2%. Table 2 lists the MAF of the 12 markers as determined by individual genotyping.
Linkage disequilibrium and haplotypes
The HelixTree genetics analysis software (Golden Helix Inc, Bozeman, USA) was used to conduct a pair-wise analysis of LD across the 12 CC chemokine gene markers, using the parents only. This software calculated both D' and r2 measures of LD, and the data are represented in Figure 2. These results were confirmed using the Haploview software package . Results suggest a clear separation of the SNPs into two haplotype blocks of moderate to strong LD, confirmed by HaploBlockFinder  (represented in Figure 2). The first block extended from the CCL2 -2581 SNP to the CCL11 67 SNP; a distance of 33.1 kb. The second block extended from the CCL14 -649 SNP to the CCL23 -289 SNP, but did not include the CCL15 136+88 SNP, which interestingly was not in LD with any of its neighbouring SNPs. This block extended across 30.9 kb.
The HelixTree software was then used to construct haplotypes only across the two clear haplotype blocks, in parents only (350–650 individuals; actual numbers varied for each SNP); and to estimate haplotype frequencies. Table 2 lists the common haplotypes (frequency >0.05) across the two haplotype blocks and their frequencies within the unaffected parents. Four common haplotypes were identified in the 5' haplotype block, with frequencies between 0.18 and 0.33. In the 3' haplotype block, three common haplotypes were identified, but the TATA haplotype was the most common by a substantial margin (frequency = 0.76).
Analysis for transmission distortion across single markers and haplotypes
Single-locus transmission disequilibrium
Analyses for transmission distortion were undertaken for the 12 chemokine gene markers in all of the trios, as well as subgroups stratified based on HLA-DRB1*1501 status, disease course and gender. Results are listed in Table 3.
No significant transmission distortion was found for any of the SNPs when analysed in all individuals. Upon stratification, a trend towards excess transmission of the CCL5 -471T allele was found in the HLA-DRB1*1501 negative group (P(uncorrected) = 0.03). Interestingly, DNA sequencing had also suggested an increased ERR for the T allele in the HLA-DRB1*1501 negative pool. For the CCL14 -649T>A SNP, transmission distortion of the T allele was found in patients with RR-MS (P(uncorrected) = 0.03). In contrast, transmission distortion of the A allele was found in SP-MS patients (P(uncorrected) = 0.02). Transmission distortion in SP-MS patients was also found for the CCL23 316C allele (P(uncorrected) = 0.02). For the CCL23 -289A>C SNP, transmission distortion of the C allele was found in male MS patients (P(uncorrected) = 0.03). This distortion was accentuated in male patients with SP-MS (P(uncorrected) = 0.007), however this result was inconclusive given the small number of individuals in this group (n = 10 informative transmissions). None of these results would survive a conservative Bonferroni correction for multiple comparisons.
Whilst single-locus analyses did not suggest that any of the SNPs were significantly associated with MS susceptibility, it is conceivable that they may define MS-associated haplotypes. Using the TRANSMIT  software, haplotype transmission was evaluated for common haplotypes (>5% frequency) of decreasing size from either block, as well as pairwise analyses across all twelve SNPs. Analyses were conducted only where genotyping information was available for all individuals across all 12 SNPs (n = 162 families), with additional trio families removed after identification of Mendelian inconsistencies.
Table 4 lists common haplotypes for which transmission distortion was identified. All of these findings were identified within the 5' haplotype block (CCL2 -2581.CCL2 -2138.CCL11 -488.CCL11 67). Marginally significant transmission distortion was found for two four-marker haplotypes, and several three-marker haplotypes (P ≤ 0.05). Marginally significant transmission distortion was also identified for pairwise analyses for CCL2 -2581.CCL2 -2138 (P(uncorrected) = 0.05) and CCL2 -2138.CCL11 -488 haplotypes (P(uncorrected) = 0.04). No significant transmission distortion was found for haplotypes across the 3' haplotype block.
Verification of haplotype results in an independent cohort
Independent validation of the results discussed above was sought in an independent cohort of 208 Australian MS trio families obtained from the Southern MS Genetics Consortium. It was determined that genotyping the CCL2 -2138A>T and CCL11 -488C>A SNPs was sufficient for information extraction across the four markers within the 5' haplotype block.
Genotype and minor allele frequencies for both SNPs were equivalent to those obtained for our cohort, and were in Hardy-Weinberg equilibrium. LD properties between the two markers were also in agreement with our results. No evidence for single-locus transmission distortion was identified for either SNP, supporting our original findings. Full haplotype transmission data could be obtained for 169 trios; losses were due to genotyping failure for either SNP. Two-marker haplotype analysis revealed slight transmission distortion for the A.C and T.C haplotypes in the same direction as our original findings, but the distortion was not significant (Table 4(b)). However, combining the total datasets of the two cohorts restored the original trends for transmission distortion for the A.C and T.C haplotypes (P = 0.05) (Table 4(c)).
In this study, we have analysed the members of the CC chemokine gene cluster for association with MS. Variant chemokine expression could diminish or enhance the inflammatory response characteristic of MS. The majority of CC chemokine genes are located in a chromosomal region (17q11.2-12) that has shown suggestive linkage and association with MS in a number of genome wide screens [29–34]. Most recently, in a large genome wide linkage screen, no genome-wide significant results could be identified beyond the MHC . However, suggestive linkage was identified on chromosome 17q23. An important deduction from this work was that linkage studies are under-powered to detect modest associations even in large cohorts.
We utilised a two-stage approach for this study, the first of which was sequencing of four DNA pools, three composed of MS patients and one control DNA pool across the 17q11.2-12 CC chemokine gene cluster. This approach allowed estimation of MAF for all common (MAF >0.15) SNPs across the CC chemokine genes in our Australian Caucasian cohort, and identification of two novel SNPs. This technique has been validated by our group . Twelve SNPs were chosen for further analysis by individual genotyping, based on their likely functional significance (codon changing, location in the putative promoter region), prevalence (MAF >0.1), and/or uneven representation between the MS and control pools.
Overall, a large amount of genetic data has been gathered from pooled DNA sequencing and individual genotyping. The pooled sequencing provided a close estimate of MAF for SNPs across the 17q chemokine gene region, and verified the existence of common SNPs published in online databases. The individual genotyping data from 350–650 unaffected parents allowed the determination of MAF for the 12 SNPs chosen for further study, establishment of LD structure across the region, and calculation of haplotype frequencies in our unaffected Australian cohort. The MAF data facilitate power calculations for additional disease-susceptibility studies, and it and the haplotype data should prove useful in studies of population variation for these immunologically important genes. From individual genotyping, MAF were established for the 12 markers in a large unaffected cohort (up to 600 unaffected parents), whilst LD analysis revealed the haplotype block structure across the cluster, and haplotype frequencies were also established.
Measurement of transmission distortion for each of the 12 SNPs in MS trios provided some evidence for over-transmission of several of the SNPs after stratification, but these trends would not survive a conservative correction for multiple comparisons.
The CCL5 -471C>T SNP, which we found might be associated with MS in HLA-DRB1*1501 negative patients, is of potential functional relevance, as it creates a new transcription factor binding site , and has shown associations with atopic dermatitis, atopy and asthma [41, 42]. CCL5 is an important inflammatory chemokine, with a range of activities upon eosinophils , monocytes and activated T cells [63, 64], and has been identified repeatedly in the CNS and CSF of MS patients [19, 21, 22]. Our data also identify CCL14 and CCL23 as warranting further investigation, with marginally significant (uncorrected) trends towards transmission distortion found in RR-MS and SP-MS patients for CCL14 -649T>A, and in males for CCL23 -289A>C. Little is known for either of these chemokines, beyond basic functionality, and neither has been studied with regard to MS pathogenesis. The CCL23 -289A>C SNP was found to affect a potential SP-1 binding site using the TRANSFAC database .
Upon analysing for haplotype transmission distortion across the CC chemokine gene cluster, we found suggestive evidence for association of haplotypes encompassing CCL2 and CCL11 (Table 4), which lie within a haplotype block spanning 33.1 kb. CCL2 has been implicated in MS pathogenesis. It is chemotactic for T cells and monocytes [66, 67], is important in the induction of inflammation in EAE , has been identified in MS lesions [17–19], and in contrast to EAE studies has been associated with remission of MS .
Recently (and subsequent to the completion of our study), a similar study of polymorphisms across the CC chemokine gene cluster was conducted by Vyshkina et al . In this study, 31 SNPs derived from online databases were selected for genotyping in a variety of individuals. The basis of selection of SNPs was not discussed. This is in contrast to our approach, in which an informed decision for SNP selection was based on several criteria, discussed above, with a particular focus on SNPs with potential functional consequences. The SNPs assayed by Vyshkina et al  were predominantly non-coding, with no SNPs from putative promoter regions. Fifteen SNPs overlapped with those identified in our pooled DNA sequencing. Two SNPs, both exonic, were individually genotyped in our study; the remainder did not pass the first stage. One exonic SNP was excluded as it was not common, whilst an exonic SNP in CCL4 was not analysed as this gene was excluded from our study. Similar to our study, Vyshkina et al  found no strong evidence for association with any single locus. In addition, we identified similar LD structure surrounding CCL2.CCL11 and between CCL14 and CCL23, and identified haplotypic associations for the CCL2.CCL11 haplotype block. Whilst the SNPs analysed in this block were different between the two studies, it might be assumed that they are subject to the LD within the block. Thus, whilst we can not specifically say that we have replicated the haplotype associations of Vyshkina et al , we would suggest that the determination of haplotype tagging SNPs in this block is justified. It was interesting that our TDT analysis was conducted entirely in sporadic MS cases, whilst that of Vyshkina et al  had an emphasis on familial MS, yet both studies found weak association for the CCL2.CCL11 haplotype block with MS, suggestive of a general MS susceptibility factor within this block. This group has since verified their results in a second-phase study . Thus, in a combined population totalling almost 1000 MS families (331 Australian; 644 North American), the results implicate a haplotype association encompassing the CCL2.CCL11 markers.
We have conducted a two-stage analysis of polymorphisms across the CC chemokine gene cluster. We identified novel SNPs in this region, and added further information to the data available on LD structure and haplotypes across the cluster. Interestingly, our best single-locus findings were for promoter SNPs (CCL5 -471C>T, CCL14 -649T>A, CCL23 -289A>C), which might affect the relative expression of these chemokines and lead to downstream effects on leukocyte migration to the CNS, and therefore influence MS pathogenesis. We also confirmed the potential presence of a haplotypic association across CCL2 and CCL11. Further validation of the association of these SNPs in independent cohorts, and confirmation of their functional significance would support therapeutic targetting of these chemokines and their receptors.
Ewing C, Bernard CC: Insights into the aetiology and pathogenesis of multiple sclerosis. Immunol Cell Biol. 1998, 76 (1): 47-54. 10.1046/j.1440-1711.1998.00718.x.
Dyment DA, Sadovnick AD, Ebers GC, Sadnovich AD: Genetics of multiple sclerosis. Hum Mol Genet. 1997, 6 (10): 1693-1698. 10.1093/hmg/6.10.1693.
Sawcer S, Compston A: The Genetic Analysis of Multiple sclerosis in EuropeanS: concepts and design. J Neuroimmunol. 2003, 143 (1–2): 13-16. 10.1016/j.jneuroim.2003.08.035.
Stewart GJ: Multiple sclerosis, multiple genes. Medical Journal of Australia. 1997, 166: 573-574.
Brown KD, Zurawski SM, Mosmann TR, Zurawski G: A family of small inducible proteins secreted by leukocytes are members of a new superfamily that includes leukocyte and fibroblast-derived inflammatory agents, growth factors, and indicators of various activation processes. J Immunol. 1989, 142 (2): 679-687.
Baggiolini M: Chemokines and leukocyte traffic. Nature. 1998, 392 (6676): 565-568. 10.1038/33340.
Gale LM, McColl SR: Chemokines: extracellular messengers for all occasions?. Bioessays. 1999, 21 (1): 17-28. 10.1002/(SICI)1521-1878(199901)21:1<17::AID-BIES3>3.0.CO;2-4.
Karpus WJ, Lukacs NW, Kennedy KJ, Smith WS, Hurst SD, Barrett TA: Differential CC chemokine-induced enhancement of T helper cell cytokine production. J Immunol. 1997, 158 (9): 4129-4136.
Taub DD, Turcovski-Corrales SM, Key ML, Longo DL, Murphy WJ: Chemokines and T lymphocyte activation: I. Beta chemokines costimulate human T lymphocyte activation in vitro. J Immunol. 1996, 156 (6): 2095-2103.
Molon B, Gri G, Bettella M, Gomez-Mouton C, Lanzavecchia A, Martinez AC, Manes S, Viola A: T cell costimulation by chemokine receptors. Nat Immunol. 2005, 6 (5): 465-471. 10.1038/ni1191.
Cappello P, Caorsi C, Bosticardo M, De Angelis S, Novelli F, Forni G, Giovarelli M: CCL16/LEC powerfully triggers effector and antigen-presenting functions of macrophages and enhances T cell cytotoxicity. J Leukoc Biol. 2004, 75 (1): 135-142. 10.1189/jlb.0403146.
Taub DD, Sayers TJ, Carter CR, Ortaldo JR: Alpha and beta chemokines induce NK cell migration and enhance NK-mediated cytolysis. J Immunol. 1995, 155 (8): 3877-3888.
Huang DR, Wang J, Kivisakk P, Rollins BJ, Ransohoff RM: Absence of monocyte chemoattractant protein 1 in mice leads to decreased local macrophage recruitment and antigen-specific T helper cell type 1 immune response in experimental autoimmune encephalomyelitis. J Exp Med. 2001, 193 (6): 713-726. 10.1084/jem.193.6.713.
Karpus WJ, Lukacs NW, McRae BL, Strieter RM, Kunkel SL, Miller SD: An important role for the chemokine macrophage inflammatory protein-1 alpha in the pathogenesis of the T cell-mediated autoimmune disease, experimental autoimmune encephalomyelitis. J Immunol. 1995, 155 (10): 5003-5010.
Miyagishi R, Kikuchi S, Takayama C, Inoue Y, Tashiro K: Identification of cell types producing RANTES, MIP-1 alpha and MIP-1 beta in rat experimental autoimmune encephalomyelitis by in situ hybridization. J Neuroimmunol. 1997, 77 (1): 17-26. 10.1016/S0165-5728(97)00040-4.
Tran EH, Kuziel WA, Owens T: Induction of experimental autoimmune encephalomyelitis in C57BL/6 mice deficient in either the chemokine macrophage inflammatory protein-1alpha or its CCR5 receptor. Eur J Immunol. 2000, 30 (5): 1410-1415. 10.1002/(SICI)1521-4141(200005)30:5<1410::AID-IMMU1410>3.0.CO;2-L.
McManus C, Berman JW, Brett FM, Staunton H, Farrell M, Brosnan CF: MCP-1, MCP-2 and MCP-3 expression in multiple sclerosis lesions: an immunohistochemical and in situ hybridization study. J Neuroimmunol. 1998, 86 (1): 20-29. 10.1016/S0165-5728(98)00002-2.
Van Der Voorn P, Tekstra J, Beelen RH, Tensen CP, Van Der Valk P, De Groot CJ: Expression of MCP-1 by reactive astrocytes in demyelinating multiple sclerosis lesions. Am J Pathol. 1999, 154 (1): 45-51.
Simpson JE, Newcombe J, Cuzner ML, Woodroofe MN: Expression of monocyte chemoattractant protein-1 and other beta-chemokines by resident glia and inflammatory cells in multiple sclerosis lesions. J Neuroimmunol. 1998, 84 (2): 238-249. 10.1016/S0165-5728(97)00208-7.
Simpson J, Rezaie P, Newcombe J, Cuzner ML, Male D, Woodroofe MN: Expression of the beta-chemokine receptors CCR2, CCR3 and CCR5 in multiple sclerosis central nervous system tissue. J Neuroimmunol. 2000, 108 (1–2): 192-200. 10.1016/S0165-5728(00)00274-5.
Hvas J, McLean C, Justesen J, Kannourakis G, Steinman L, Oksenberg JR, Bernard CC: Perivascular T cells express the pro-inflammatory chemokine RANTES mRNA in multiple sclerosis lesions. Scand J Immunol. 1997, 46 (2): 195-203. 10.1046/j.1365-3083.1997.d01-100.x.
Sorensen TL, Tani M, Jensen J, Pierce V, Lucchinetti C, Folcik VA, Qin S, Rottman J, Sellebjerg F, Strieter RM, et al: Expression of specific chemokines and chemokine receptors in the central nervous system of multiple sclerosis patients. J Clin Invest. 1999, 103 (6): 807-815.
Sorensen TL, Sellebjerg F, Jensen CV, Strieter RM, Ransohoff RM: Chemokines CXCL10 and CCL2: differential involvement in intrathecal inflammation in multiple sclerosis. Eur J Neurol. 2001, 8 (6): 665-672. 10.1046/j.1468-1331.2001.00327.x.
Strunk T, Bubel S, Mascher B, Schlenke P, Kirchner H, Wandinger KP: Increased numbers of CCR5+ interferon-gamma- and tumor necrosis factor-alpha-secreting T lymphocytes in multiple sclerosis patients. Ann Neurol. 2000, 47 (2): 269-273. 10.1002/1531-8249(200002)47:2<269::AID-ANA23>3.0.CO;2-G.
Miyagishi R, Kikuchi S, Fukazawa T, Tashiro K: Macrophage inflammatory protein-1 alpha in the cerebrospinal fluid of patients with multiple sclerosis and other inflammatory neurological diseases. J Neurol Sci. 1995, 129 (2): 223-227. 10.1016/0022-510X(95)00004-L.
Teleshova N, Pashenkov M, Huang YM, Soderstrom M, Kivisakk P, Kostulas V, Haglund M, Link H: Multiple sclerosis and optic neuritis: CCR5 and CXCR3 expressing T cells are augmented in blood and cerebrospinal fluid. J Neurol. 2002, 249 (6): 723-729. 10.1007/s00415-002-0699-z.
Mahad D, Callahan MK, Williams KA, Ubogu EE, Kivisakk P, Tucky B, Kidd G, Kingsbury GA, Chang A, Fox RJ, et al: Modulating CCR2 and CCL2 at the blood-brain barrier: relevance for multiple sclerosis pathogenesis. Brain. 2006, 129 (Pt 1): 212-223.
Naruse K, Ueno M, Satoh T, Nomiyama H, Tei H, Takeda M, Ledbetter DH, Coillie EV, Opdenakker G, Gunge N, et al: A YAC contig of the human CC chemokine genes clustered on chromosome 17q11.2. Genomics. 1996, 34 (2): 236-240. 10.1006/geno.1996.0274.
Giedraitis V, Modin H, Callander M, Landtblom AM, Fossdal R, Stefansson K, Hillert J, Gulcher J: Genome-wide TDT analysis in a localized population with a high prevalence of multiple sclerosis indicates the importance of a region on chromosome 14q. Genes Immun. 2003, 4 (8): 559-563. 10.1038/sj.gene.6364024.
Hensiek AE, Roxburgh R, Smilie B, Coraddu F, Akesson E, Holmans P, Sawcer SJ, Compston DA: Updated results of the United Kingdom linkage-based genome screen in multiple sclerosis. J Neuroimmunol. 2003, 143 (1–2): 25-30. 10.1016/j.jneuroim.2003.08.007.
Kuokkanen S, Gschwend M, Rioux JD, Daly MJ, Terwilliger JD, Tienari PJ, Wikstrom J, Palo J, Stein LD, Hudson TJ, et al: Genomewide scan of multiple sclerosis in Finnish multiplex families. Am J Hum Genet. 1997, 61 (6): 1379-1387. 10.1086/301637.
Sawcer S, Jones HB, Feakes R, Gray J, Smaldon N, Chataway J, Robertson N, Clayton D, Goodfellow PN, Compston A: A genome screen in multiple sclerosis reveals susceptibility loci on chromosome 6p21 and 17q22. Nat Genet. 1996, 13 (4): 464-468. 10.1038/ng0896-464.
Sawcer S, Maranian M, Setakis E, Curwen V, Akesson E, Hensiek A, Coraddu F, Roxburgh R, Sawcer D, Gray J, et al: A whole genome screen for linkage disequilibrium in multiple sclerosis confirms disease associations with regions previously linked to susceptibility. Brain. 2002, 125 (Pt 6): 1337-1347. 10.1093/brain/awf143.
Dyment DA, Sadovnick AD, Willer CJ, Armstrong H, Cader ZM, Wiltshire S, Kalman B, Risch N, Ebers GC: An extended genome scan in 442 Canadian multiple sclerosis-affected sibships: a report from the Canadian Collaborative Study Group. Hum Mol Genet. 2004, 13 (10): 1005-1015. 10.1093/hmg/ddh123.
The Transatlantic Multiple Sclerosis Genetics Cooperative: A meta-analysis of genomic screens in multiple sclerosis. Mult Scler. 2001, 7 (1): 3-11. 10.1191/135245801669625359.
Jagodic M, Becanovic K, Sheng JR, Wu X, Backdahl L, Lorentzen JC, Wallstrom E, Olsson T: An advanced intercross line resolves eae18 into two narrow quantitative trait Loci syntenic to multiple sclerosis candidate Loci. J Immunol. 2004, 173 (2): 1366-1373.
Teuscher C, Butterfield RJ, Ma RZ, Zachary JF, Doerge RW, Blankenhorn EP: Sequence polymorphisms in the chemokines Scya1 (TCA-3), Scya2 (monocyte chemoattractant protein (MCP)-1), and Scya12 (MCP-5) are candidates for eae7, a locus controlling susceptibility to monophasic remitting/nonrelapsing experimental allergic encephalomyelitis. J Immunol. 1999, 163 (4): 2262-2266.
Jamieson SE, Miller EN, Black GF, Peacock CS, Cordell HJ, Howson JM, Shaw MA, Burgner D, Xu W, Lins-Lainson Z, et al: Evidence for a cluster of genes on chromosome 17q11-q21 controlling susceptibility to tuberculosis and leprosy in Brazilians. Genes Immun. 2004, 5 (1): 46-57. 10.1038/sj.gene.6364029.
Modi WS, Goedert JJ, Strathdee S, Buchbinder S, Detels R, Donfield S, O'Brien SJ, Winkler C: MCP-1-MCP-3-Eotaxin gene cluster influences HIV-1 transmission. Aids. 2003, 17 (16): 2357-2365. 10.1097/00002030-200311070-00011.
Gonzalez-Escribano MF, Torres B, Aguilar F, Rodriguez R, Garcia A, Valenzuela A, Nunez-Roldan A: MCP-1 promoter polymorphism in spanish patients with rheumatoid arthritis. Hum Immunol. 2003, 64 (7): 741-744. 10.1016/S0198-8859(03)00090-9.
Nickel RG, Casolaro V, Wahn U, Beyer K, Barnes KC, Plunkett BS, Freidhoff LR, Sengler C, Plitt JR, Schleimer RP, et al: Atopic dermatitis is associated with a functional mutation in the promoter of the C-C chemokine RANTES. J Immunol. 2000, 164 (3): 1612-1616.
Fryer AA, Spiteri MA, Bianco A, Hepple M, Jones PW, Strange RC, Makki R, Tavernier G, Smilie FI, Custovic A, et al: The -403 G-->A promoter polymorphism in the RANTES gene is associated with atopy and asthma. Genes Immun. 2000, 1 (8): 509-514. 10.1038/sj.gene.6363717.
Bennetts BH, Teutsch SM, Buhler MM, Heard RN, Stewart GJ: The CCR5 deletion mutation fails to protect against multiple sclerosis. Hum Immunol. 1997, 58 (1): 52-59. 10.1016/S0198-8859(97)00207-3.
Barcellos LF, Schito AM, Rimmler JB, Vittinghoff E, Shih A, Lincoln R, Callier S, Elkins MK, Goodkin DE, Haines JL, et al: CC-chemokine receptor 5 polymorphism and age of onset in familial multiple sclerosis. Multiple Sclerosis Genetics Group. Immunogenetics. 2000, 51 (4–5): 281-288. 10.1007/s002510050621.
Kantor R, Bakhanashvili M, Achiron A: A mutated CCR5 gene may have favorable prognostic implications in MS. Neurology. 2003, 61 (2): 238-240.
Vyshkina T, Shugart YY, Birnbaum G, Leist TP, Kalman B: Association of haplotypes in the beta-chemokine locus with multiple sclerosis. Eur J Hum Genet. 2005, 13 (2): 240-247. 10.1038/sj.ejhg.5201295.
Poser CM, Paty DW, Scheinberg L, McDonald WI, Davis FA, Ebers GC, Johnson KP, Sibley WA, Silberberg DH, Tourtellotte WW: New diagnostic criteria for multiple sclerosis: guidelines for research protocols. Ann Neurol. 1983, 13 (3): 227-231. 10.1002/ana.410130302.
Lahiri DK, Nurnberger JI: A rapid non-enzymatic method for the preparation of HMW DNA from blood for RFLP studies. Nucleic Acids Res. 1991, 19 (19): 5444-
Teutsch SM, Booth DR, Bennetts BH, Heard RN, Stewart GJ: Identification of 11 novel and common single nucleotide polymorphisms in the interleukin-7 receptor-alpha gene and their associations with multiple sclerosis. Eur J Hum Genet. 2003, 11 (7): 509-515. 10.1038/sj.ejhg.5200994.
Ban M, Sawcer SJ, Heard RN, Bennetts BH, Adams S, Booth D, Perich V, Setakis E, Compston A, Stewart GJ: A genome-wide screen for linkage disequilibrium in Australian HLA-DRB1*1501 positive multiple sclerosis patients. J Neuroimmunol. 2003, 143 (1–2): 60-64. 10.1016/j.jneuroim.2003.08.012.
Norton N, Williams NM, Williams HJ, Spurlock G, Kirov G, Morris DW, Hoogendoorn B, Owen MJ, O'Donovan MC: Universal, robust, highly quantitative SNP allele frequency measurement in DNA pools. Hum Genet. 2002, 110 (5): 471-478. 10.1007/s00439-002-0706-6.
den Dunnen JT, Antonarakis E: Nomenclature for the description of human sequence variations. Hum Genet. 2001, 109: 121-124. 10.1007/s004390100505.
Spielman RS, McGinnis RE, Ewens WJ: Transmission test for linkage disequilibrium: the insulin gene region and insulin-dependent diabetes mellitus (IDDM). Am J Hum Genet. 1993, 52 (3): 506-516.
Kruglyak L, Daly MJ, Reeve-Daly MP, Lander ES: Parametric and nonparametric linkage analysis: a unified multipoint approach. Am J Hum Genet. 1996, 58 (6): 1347-1363.
Zhang K, Jin L: HaploBlockFinder: haplotype block analyses. Bioinformatics. 2003, 19 (10): 1300-1301. 10.1093/bioinformatics/btg142.
Clayton D: A generalization of the transmission/disequilibrium test for uncertain-haplotype transmission. Am J Hum Genet. 1999, 65 (4): 1170-1177. 10.1086/302577.
Abecasis GR, Cherny SS, Cookson WO, Cardon LR: Merlin – rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet. 2002, 30 (1): 97-101. 10.1038/ng786.
Altshuler D, Brooks LD, Chakravarti A, Collins FS, Daly MJ, Donnelly P: A haplotype map of the human genome. Nature. 2005, 437 (7063): 1299-1320. 10.1038/nature04226.
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.
Sawcer S, Ban M, Maranian M, Yeo TW, Compston A, Kirby A, Daly MJ, De Jager PL, Walsh E, Lander ES, et al: A high-density screen for linkage in multiple sclerosis. Am J Hum Genet. 2005, 77 (3): 454-467. 10.1086/444547.
Bugeja MJ, Booth DR, Bennetts BH, Heard RNS, Stewart GJ: An investigation of polymorphisms in the 4q13.3-21.1 CXCchemokine gene cluster for association with multiple sclerosis inAustralians. Mult Scler. 2006.
Kameyoshi Y, Dorschner A, Mallet AI, Christophers E, Schroder JM: Cytokine RANTES released by thrombin-stimulated platelets is a potent attractant for human eosinophils. J Exp Med. 1992, 176 (2): 587-592. 10.1084/jem.176.2.587.
Schall TJ: Biology of the RANTES/SIS cytokine family. Cytokine. 1991, 3 (3): 165-183. 10.1016/1043-4666(91)90013-4.
Schall TJ, Bacon K, Toy KJ, Goeddel DV: Selective attraction of monocytes and T lymphocytes of the memory phenotype by cytokine RANTES. Nature. 1990, 347 (6294): 669-671. 10.1038/347669a0.
Wingender E, Chen X, Hehl R, Karas H, Liebich I, Matys V, Meinhardt T, Pruss M, Reuter I, Schacherer F: TRANSFAC: an integrated system for gene expression regulation. Nucleic Acids Res. 2000, 28 (1): 316-319. 10.1093/nar/28.1.316.
Taub DD, Proost P, Murphy WJ, Anver M, Longo DL, van Damme J, Oppenheim JJ: Monocyte chemotactic protein-1 (MCP-1), -2, and -3 are chemotactic for human T lymphocytes. J Clin Invest. 1995, 95 (3): 1370-1376.
Loetscher P, Seitz M, Clark-Lewis I, Baggiolini M, Moser B: Monocyte chemotactic proteins MCP-1, MCP-2, and MCP-3 are major attractants for human CD4+ and CD8+ T lymphocytes. Faseb J. 1994, 8 (13): 1055-1060.
Mahad DJ, Howell SJ, Woodroofe MN: Expression of chemokines in the CSF and correlation with clinical disease activity in patients with multiple sclerosis. J Neurol Neurosurg Psychiatry. 2002, 72 (4): 498-502.
Vyshkina T, Kalman B: Haplotypes within genes of beta-chemokines in 17q11 are associated with multiple sclerosis: a second phase study. Hum Genet. 2005, 1-9.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/7/64/prepub
The authors would like to thank the members of the Southern MS Genetics Consortium: Justin Rubio, Trevor Kilpatrick, Helmut Butzkueven, Niall Tubridy, Mark Mariott, Caron Chapman, John Carey, Jo Baker, Laura Johnson, Rachel Tan, Simon Foote, Stewart Huxtable, Melanie Bahlo, Jim Stankovich, and Terry Speed, Najwa Marmash for technical assistance, Susan Adams for recruitment of MS trio patients, Rex Simmons for providing familial MS cases, Vanessa Perich, Maria Ban, Natalie Hartley and Suzy Teutsch for preparation of the DNA samples, Mark Wheeler and Ilya Henner for running samples for sequencing and SNaPshot analysis, and Mike Payne for overseeing the genotyping of our samples by SNPlex™. This work was supported by a grant from the National Health and Medical Research Council of Australia (grant no. 153990). M.B. is the recipient of an Australian Postgraduate Award Scholarship and additional support was provided by Biogen Australia.
The author(s) declare that they have no competing interests.
MJB was responsible for the conception of the project, acquisition and analysis of the data, interpretation of results and writing of the manuscript. DB, BB and GS were responsible for the conception and design of the project, intellectual input, interpretation of results, and helped to draft the manuscript. RH was responsible for the clinical aspects of this project, including the assessment of MS patients. All authors read and approved the final manuscript.