|
|
||||||||


* Department of Clinical Neuroscience, Neuroimmunology Unit, and
Department of Medicine, Rheumatology Unit, Center for Molecular Medicine, Karolinska Institutet, Stockholm, Sweden
| Abstract |
|---|
|
|
|---|
30-Mb large region. Fine-mapping was then performed in an advanced intercross line consisting of a (DA x PVG.1AV1)F7 intercross, resulting in two adjacent EAE-regulating QTLs designated Eae18a and Eae18b. The two QTLs span 5.5 and 3 Mb, respectively, and the 3-Mb Eae18b contains as few as 10 genes, including a cluster of chemokine genes (CCL1, CCL2, CCL7, and CCL11). Eae18a and Eae18b are syntenic to human chromosome 17p13 and 17q11, respectively, which both display linkage to multiple sclerosis. Thus, Eae18 consists of at least two EAE-regulating genes, providing additional evidence that clustering of disease-regulating genes in QTLs is an important phenomenon. The overlap between Eae18a and Eae18b with previously identified QTLs in humans and mice further supports the notion that susceptibility alleles in inflammatory disease are evolutionary conserved between species. | Introduction |
|---|
|
|
|---|
We here report the first application of an advanced intercross line (AIL) in the genetic dissection of a complex experimental inflammatory disease. To date, gene mapping in rodent inflammatory diseases has relied on whole genome scans in backcross or F2 populations followed by fine-mapping in congenic strains. Although
35 quantitative trait loci (QTLs) regulating rodent EAE have been identified with this approach (17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31), there are certain limitations. Confidence intervals in F2 crosses are typically
20 cM and contain, on average, as many as 500 genes. Further complications using linkage analysis in F2 populations are: 1) a single QTL may comprise several adjacent QTLs (32, 33); 2) nearby QTLs with opposing effects may not be detected (34); and 3) a detected QTL may be incorrectly positioned (33, 35). Congenic mapping may provide formal evidence for the role of a particular gene, although it requires laborious testing because initially even a larger fragment than the one representing the F2 confidence interval is introduced into the recipient background genome (36). High resolution QTL mapping can instead be achieved using an AIL. An AIL is created by random intercross breeding of two inbred strains for several generations, resulting in genetically unique individuals with a mixture of founder chromosomal fragments. The increase in recombinations increases the genomic resolutions, which allows separation of linked QTLs, given that t (t = number of generations) is large enough (35). Furthermore, the relative distance between recombinations is 100/t cM; thus, the AIL theoretically gives at least a t/2-fold reduction in the confidence interval compared with a F2 cross (35, 37). The estimation of the QTL localization will therefore be considerably improved.
Our interest in the region on rat chromosome 10, hereby designated Eae18 (http://ratmap.gen.gu.se/), stems from previous observations of linkage to MOG-EAE (17). We have demonstrated that a rat congenic strain with a large ACI chromosome 10 fragment (
50 cM, corresponding to
60 Mb) introduced into the susceptible DA background conferred protection from MOG-EAE (38). In the present study, we reduced the size of Eae18 in recombinant DA.ACI congenic strains. Eae18 was then resolved into two narrow QTLs using an F7 (DA x PVG.1AV1) AIL.
| Materials and Methods |
|---|
|
|
|---|
Full length (DA.ACI-D10Rat2-D10Rat29) congenic animals were established using a speed congenization approach (38). They were then intercrossed with parental DA rats to create heterozygous animals that were subsequently intercrossed to obtain recombinations within the region. Rats with recombinations were made homozygous before experiments. The following recombinant strains were tested in MOG-EAE (Fig. 1): R1 (DA.ACI-D10Rat15-D10Rat29) and R2 (DA.ACI-D10Rat2-D10Rat6) that did not overlap, R9 (DA.ACI-D10Rat10-D10Rat142) that covered the region between R1 and R2 overlapping with both, R10 (DA.ACI-D10Rat12-D10Rat144) that overlapped with R1 and R9 and finally R5 (DA.ACI-D10Rat219-D10Rat29) that covered the central
15 cM (
30 Mb) large chromosome 10 fragment of the R1 recombinant. Between 7 and 16 male rats were used in each experiment.
|
60% according to the Whitehead Institute (http://www.genome.wi.mit.edu/rat/public/). To create the F1 generation, two breeding pairs with DA female founders and two breeding pairs with PVG.1AV1 female founders were bred. The F2 generation was produced from seven couples each of F1 rats with DA and PVG.1AV1 as female founders, respectively. The F3 generation originated from 50 breeding couples with both types of female founders. Random breeding of 50 males and females, consistently avoiding brother-sister mating, produced all subsequent generations. Thereafter, three F7 litters were produced for MOG-EAE experiments. The litters were similar in size and contained almost equal numbers of females and males. Inbred DA (Dark Agouti) and PVG.1AV1 rats were originally obtained from Zentralinstitut für Versuchstierzucht (Hannover, Germany). Rats were bred and kept at the Karolinska Hospital (Stockholm, Sweden) in a 12-h light/dark cycle, housed in polystyrene cages containing aspen wood shavings and with free access to water and standard rodent chow. They were routinely tested for specific pathogens according to a health-monitoring program for rats at the National Veterinary Institute in Uppsala, Sweden. The local ethical committee in Northern Stockholm approved the experiments.
Phenotyping
Recombinant rat MOG (aa 1125 from the N terminus) was expressed in Escherichia coli and purified to homogeneity by chelate chromatography (39). Rats were anesthetized with isoflurane (Forene; Abbott Laboratories, Abbott Park, IL) and immunized by a single intradermal injection in the dorsal base of the tail with 200 µl of inoculum containing rMOG in saline emulsified (1:1) with IFA (Sigma-Aldrich, St. Louis, MO). Rats were weighed and monitored daily for clinical signs of EAE, from day 10 until day of sacrifice at day 3140 postimmunization (p.i.). Serum for anti-MOG IgG assessment was collected on day 12 p.i. The clinical scoring was: 0, no clinical signs of EAE; 1, tail weakness or tail paralysis; 2, hind leg paraparesis or hemiparesis; 3, hind leg paralysis or hemiparalysis; 4, tetraplegy or moribund; and 5, death. The rMOG dose used for disease induction in the congenic mapping experiments was 15, 20, or 7080 µg/rat, depending on the batch of rMOG protein, and 40 µg for AIL animals. The dose was selected upon titration in the susceptible parental DA rats. Age-matched rats between 9 and 15 wk of age were used in all experiments.
Genotyping and linkage analysis
A total of 152 clinically affected rats and 162 randomly selected unaffected rats were genotyped. Affected animals were selected on the basis of displaying unambiguous signs of the disease (minimum score 1 for >2 days accompanied with weight loss). Rats in the unaffected group did not display any signs of disease, including a steady increase in weight. DNA was extracted from the tail tip according to a standard protocol (40). The region analyzed in the AIL covered the fragment in the R5 congenic strain and the centromeric fragment of the R10 congenic strain (Fig. 1). Both of these fragments displayed phenotypic differences. The fragment in the R9 congenic strain did not display any phenotypic differences and therefore the telomeric part of the R10 fragment, which overlaps with the R9 fragment, was not included in the AIL analysis (Fig. 1). The 35-cM (50 Mb) large region, extending from D10Rat78 to D10Rat142, was genotyped with 56 microsatellite markers (<1 cM spacing in average) (GENSET, Paris, France). Primers for additional microsatellite markers within the region were designed using Oligo version 5.0 (National Biosciences, Plymouth, MN). The following primers were used: MJ082 (forward 5'-AGAGACCAGACAGACTTTAAT; reverse 5'-TGTACGTGTGGTACGTAAATA); MJ411 (forward 5'-GCATCTGTATTCATATTT; reverse 5'-TATACATGTTTCACTGCAAAT), MJ770 (forward 5'-AGTAATTGCTGTCCTTTTCT; reverse 5'-GGATCAACAATTTTGAAAGT); MJ546 (forward 5'-TTTGCTCTTCAGGTCTAGTA; reverse 5'-TCCTCCATGTCTTATTTGTA) and MJ310 (forward 5'-GATCTCTGTAGTCCCATAAAT; reverse 5'-GGAAAATAGAGACTTAGAGAA). PCR amplification was performed as previously described (41) with [
-33P]ATP end-labeling of the forward primer. The PCR products were size fractionated on 6% polyacrylamide gels and visualized by autoradiography. All genotypes were evaluated manually and double-checked.
The map, defined using the publicly available rat genome sequence (http://www.ensembl.org/), corresponded to the linkage map of the present data generated using MAPMAKER/EXP. Accordingly, the marker order from the physical map was applied, and linkage analysis was performed using R/qtl software (42). Data were analyzed implementing the nonparametric model for maximum EAE score, cumulative EAE score, duration and onset of EAE, and the binary model for incidence of EAE. Permutation tests using R/qtl were performed to determine the threshold levels for significant linkage. The permutation procedure based on the investigated material is empirical and reflects the characteristics of the particular experiment to which it is applied (43). To further evaluate the identified QTLs, we have implemented a multiple QTL model test using R/qtl software. This analysis is based on the creation of an initial model of phenotypic variance, comprising all identified QTLs displaying significant linkage. Each QTL is then excluded from the initial model in subsequent steps, and the influence on the phenotypic variance by the excluded QTL is estimated. The outcome of the multiple QTL model analysis therefore provides information on whether certain QTL has a significant effect on the phenotype (p
0.05) or not (p > 0.05). The region was also tested for epistatic interactions by implementing a two-dimensional scan with a two-QTL model using R/qtl (42).
The results were confirmed by nonparametric pointwise analysis (JMP, 4.0.2; SAS Institute, Cary, NC). Fishers exact test was used to analyze whether there was a difference in the observed genotype distribution between affected and unaffected rats. The Kruskal-Wallis ranking test was used to analyze whether different genotypes were associated with maximum EAE score, cumulative EAE score, duration of EAE, onset of EAE, and anti-MOG Ab levels, respectively. The analysis using R/qtl and MAPMAKER/QTL resulted in similar levels of significance and peak intervals (44). We therefore present the outcome of R/qtl alone, because this software was used to determine the significance levels using permutation tests, for the fit-multiple QTL model and for the analysis of epistatic interactions.
Anti-MOG Ab determination
Blood samples for Ab measurement were collected at day 12 p.i. ELISA plates with 96 wells (Nunc, Roskilde, Denmark) were coated with 10 µg/ml rat recombinant MOG diluted in 0.1M NaHCO3 (100 µl/well) and stored overnight at 4°C. Plates were washed with PBS, 0.05% Tween 20, and free binding sites were blocked with 5% fat-free milk in PBS, 0.05% Tween 20 for 1 h at room temperature. After washing, diluted serum samples were added in duplicates and incubated for 1 h at room temperature. Plates were washed and incubated for 1 h at room temperature with rabbit anti-rat IgG, IgG1, IgG2b, or IgG2c (Nordic, Tilburg, The Netherlands). Unbound Abs were removed by washing before addition of a peroxidase-conjugated goat anti-rabbit antiserum (Nordic) diluted in PBS, 0.05% Tween 20 (1/10,000). After 30 min of incubation, plates were washed extensively, and bound Abs were visualized with 3,3',5,5'-tetramethylbenzidine (Sigma-Aldrich). The reaction was stopped by addition of 1 M HCl after 15 min in darkness and OD was read at 450 nm. Nine 2-fold dilutions in duplicates of serum containing high levels of anti-MOG IgG were included on each plate to construct a standard curve and to permit interplate comparisons. Levels were presented in arbitrary units deduced from the standard curve.
| Results |
|---|
|
|
|---|
A DA.ACI -D10Rat2-D10Rat29 congenic rat strain in which the inserted ACI congenic fragment is
50 cM (60 Mb) reproducibly demonstrated protection against MOG-EAE (38). Congenic dissection of recombinant strains demonstrated that clinical EAE features differed between the full length congenic and the recombinant R1 and R5 strains compared with the parental DA strain (Table I and Fig. 1.) The difference in clinical EAE signs varied between experiments, but the full length congenic and the recombinant R1 and R5 strains consistently displayed lower mean cumulative scores compared with the susceptible DA parental strain (Table I). In contrast, R2 and R9 did not differ from the parental DA strain (Table I). The absence of any up- and/or down-regulatory effects in the R2 and R9 recombinant congenic strains strongly argues against effects by contaminating ACI genes outside the congenic fragment. The R10 strain displayed an intermediate phenotype in a single experiment and would need further testing before concluding additional gene regulatory influences from this region (Table I). We have thus isolated and narrowed down a 50-cM (60-Mb) large MOG-EAE-regulating region to a 15-cM (30-Mb) region present in R1 and R5, but not in R2 and R9.
|
A total of 1068 MOG immunized F7 rats were examined for clinical EAE signs for 31 days p.i. (Table II). Clinical signs defined as a minimum of 2 days of EAE score 1 (tail paralysis) were recorded in 14.8% of animals (158 of 1068). Affected animals displayed different degrees of severity and a variety of disease courses (i.e., monophasic, relapsing-remitting, and chronic progressive EAE). The disease course could not be determined unambiguously due to a late onset in a subgroup of 33 rats. An additional group of 92 animals displayed minor disease signs with tail tip paralysis, ragged fur, weight loss, or arrested weight gain. These signs may potentially reflect a subclinical disease. We therefore excluded this group from genetic analyses, because it was not possible to separate them from the healthy rats using the current EAE scoring scale.
|
|
5.5-Mb-large region. The DA allele at the peak marker (OT18.06) is disease-enhancing in an additive manner. Sequence alignments and map comparisons (http://www.ensembl.org/) revealed that Eae18a is syntenic to human 17p13.1-p13.2 (Fig. 2 and Table III).
|
3 Mb. The DA/PVG allele distribution tested for the peak marker (D10Rat123) suggests a recessive DA allele disease-promoting effect. Eae18b is syntenic to the human chromosome 17q11.2-q12 region (Fig. 2 and Table III). A summary of all genes residing within Eae18b is presented in Table IV.
|
0.05). This test confirms a regulatory effect of Eae18a and Eae18b on EAE phenotypes. The third peak displaying LOD scores above the threshold for significance did not fit in the multiple QTL model described (p > 0.05). Moreover, it was supported by only one microsatellite marker (D10Rat81) and may therefore be a false positive. In addition, we wanted to evaluate any possible gene interactions within this region. However, no significant epistatic interactions were detected between Eae18a and Eae18b, or any other two positions within the tested region, when using a two-dimensional scan with a two-QTL model (using R/qtl; data not included).
In addition to the clinical EAE phenotypes, we assessed the anti-MOG autoimmune response by measurement of anti-MOG Ab levels in serum. Clinical disease was associated with elevated levels of total anti-MOG IgG, IgG2b, and IgG2c isotype levels, whereas there was no correlation with IgG1 isotype levels (Fig. 3). Nonetheless, Eae18a and Eae18b did not display any linkage to total anti-MOG IgG or to anti-MOG IgG1, IgG2b, or IgG2c isotypes when performing linkage analysis (data not included).
|
50-Mb-large region on rat chromosome 10 in rats that contained at least two recombinations within this region. The average fragment size was
13 Mb, corresponding to
9 cM. However, this is an underestimation because 6 of the 20 selected rats carried only one or no recombination, although the fragment size roughly fits the theoretically expected value in the seventh generation, i.e., 14 cM (100/7). | Discussion |
|---|
|
|
|---|
In our study, we combine the AIL and the recombinant congenic mapping approaches. An
60-Mb region previously isolated in the DA.ACI congenic strain was narrowed down to
30 Mb in the recombinant congenic strain DA.ACI-R5. Further congenic fine-mapping was, however, hampered by the lack of polymorphic microsatellite markers between the DA and the ACI strains, because the polymorphism rate is only
10% (17). Furthermore, the lack of polymorphic markers is not evenly distributed within this region but instead covers certain blocks. The DA.ACI-R5 region is one of those, with 3 polymorphic of a total of 70 markers in a 40-Mb region. We therefore currently define polymorphic SNPs to select intraregional recombinant strains. The low degree of polymorphisms between the DA and the ACI strains may, however, also be advantageous, because polymorphisms responsible for phenotypic differences and disease regulation are easier to identify during sequencing of candidate genes. The DA x PVG.1AV1 strain combination provides a substantially higher degree of polymorphism (
60% according to http://www-genome.wi.mit.edu/rat/public/), which is required for a successful high-resolution mapping in an AIL. Obviously, the use of the DA x PVG.1AV1 strain combination may unravel different genes than those in DA x ACI, although it is highly unlikely that none of the DA disease-promoting genes would be the same.
Eae18b was refined to a 3-Mb region, which contains 10 genes, including a cluster of chemokine genes, CCL1, CCL2, CCL7, and CCL11 (Table IV). These genes are particularly attractive candidates in EAE and MS, because chemokines are the most important mediators of leukocyte migration to inflammatory lesions. Blocking of chemokine receptors has been a successful therapeutic approach in EAE and is now also tested in MS (48, 49). Any polymorphism(s) in a certain chemokine gene that can be linked to EAE susceptibility will provide a more specific therapeutic target. Multiple sequence polymorphisms are already identified in mouse CCL1, CCL2, and CCL12 (50). These polymorphisms were suggested to alter the chemokine activity and the capability of binding to tissue matrices. Interestingly, two polymorphisms in the promoter/enhancer region of CCL7 were associated with MS in two case-control studies (51, 52). Sequencing of all genes identified in Eae18b is currently under way in our laboratory.
Eae18a was refined to a 5.5-Mb interval, situated within a very gene-dense genome region, consisting of
130 genes. A current list of genes located in Eae18a can easily be retrieved using the flanking marker information and the publicly available rat genome map (http://www.ensembl.org). Eae18a contains several relevant candidate genes, such as arachidonate 12/15-lipoxygenase, two translation initiation factors (EIF-4A and 5A), MAPK, and C1q.
Both Eae18a and Eae18b display an interesting overlap with loci identified in several human inflammatory diseases and in their experimental counterparts. Eae18b overlaps with two QTLs, which displayed linkage to murine EAE (27, 28, 29, 31). Moreover, Eae18b is syntenic to human 17q11.2-q12, which has displayed linkage to MS in several studies (4, 6, 53). Eae18a overlaps with Eae3 that was identified in an LEW x BN intercross (19). Furthermore, it is syntenic to the human chromosome 17p13.117p13.2, which displayed linkage in the British and Australian MS materials (4, 6, 12). The overlap between Eae18a and Eae18b and the regions implicated in human and mouse linkage studies strengthens the approach of comparative genomics. Eae18a and Eae18b also overlap with regions that displayed linkage to other experimental organ-specific inflammatory diseases such as neuritis and arthritis (54, 55). In addition, the syntenic mouse region on chromosome 11B3-B5 has been linked to type 1 diabetes and systemic lupus erythematosus (56, 57), whereas the human 17p13 and 17q11 regions have been linked to rheumatoid arthritis (58, 59). Accordingly, Eae18a and Eae18b are an example of clustering of genes of importance for different inflammatory diseases (60).
We conclude that for the first time in the setting of an experimental inflammatory disease, an AIL is successfully used to define two highly significant, closely situated QTLs. The data support that clustering of immunoregulatory genes in QTLs is a rule rather than an exception. Furthermore, we demonstrate synteny between the identified rat MOG-EAE-regulatory QTLs with MS QTLs on human chromosome 17. Eae18b comprises only 10 genes, which makes possible both positional cloning attempts using recombinant congenic strains, expression and sequence polymorphism analysis of the genes within the region, and association studies in large cohorts of persons with inflammatory disease and controls. Interesting candidate genes also emerge from the
130 genes located in Eae18a. Positional identification of evolutionary conserved key immunoregulatory genes is crucial for the understanding of the organ-specific chronic inflammatory diseases.
| Acknowledgments |
|---|
| Footnotes |
|---|
2 Address correspondence and reprint requests to Dr. Maja Jagodic, Neuroimmunology Unit, Center for Molecular Medicine, L8:04, S-17176 Stockholm, Sweden. E-mail address: Maja.Jagodic{at}cmm.ki.se ![]()
3 Abbreviations used in this paper: MS, multiple sclerosis; EAE, experimental autoimmune encephalomyelitis; LOD, base 10 logarithm of the likelihood ratio; MOG, myelin oligodendrocyte glycoprotein; p.i., postimmunization; QTL, quantitative trait locus. ![]()
Received for publication February 9, 2004. Accepted for publication May 11, 2004.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. Callander, S. Haghighi, A.-M. Landtblom, C.E. Ahlgren, S.I. Nilsson, L. Rydberg, H. Al Khoury, L. Rosengren, and O. Andersen Multiple sclerosis immunopathic trait and HLA-DR(2)15 as independent risk factors in multiple sclerosis Multiple Sclerosis, May 1, 2007; 13(4): 441 - 445. [Abstract] [PDF] |
||||
![]() |
X. Yu, K. Bauer, P. Wernhoff, D. Koczan, S. Moller, H.-J. Thiesen, and S. M. Ibrahim Fine Mapping of Collagen-Induced Arthritis Quantitative Trait Loci in an Advanced Intercross Line J. Immunol., November 15, 2006; 177(10): 7042 - 7049. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Ockinger, P. Serrano-Fernandez, S. Moller, S. M. Ibrahim, T. Olsson, and M. Jagodic Definition of a 1.06-Mb Region Linked to Neuroinflammation in Humans, Rats and Mice Genetics, July 1, 2006; 173(3): 1539 - 1545. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Becanovic, M. Jagodic, J. R. Sheng, I. Dahlman, F. Aboul-Enein, E. Wallstrom, P. Olofsson, R. Holmdahl, H. Lassmann, and T. Olsson Advanced Intercross Line Mapping of Eae5 Reveals Ncf-1 and CLDN4 as Candidate Genes for Experimental Autoimmune Encephalomyelitis J. Immunol., May 15, 2006; 176(10): 6055 - 6064. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. R. Sheng, M. Jagodic, I. Dahlman, K. Becanovic, R. Nohra, M. Marta, E. Iacobaeus, T. Olsson, and E. Wallstrom Eae19, a New Locus on Rat Chromosome 15 Regulating Experimental Autoimmune Encephalomyelitis Genetics, May 1, 2005; 170(1): 283 - 289. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Jagodic, M. Marta, K. Becanovic, J. R. Sheng, R. Nohra, T. Olsson, and J. C. Lorentzen Resolution of a 16.8-Mb Autoimmunity-Regulating Rat Chromosome 4 Region into Multiple Encephalomyelitis Quantitative Trait Loci and Evidence for Epistasis J. Immunol., January 15, 2005; 174(2): 918 - 924. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |