PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of jmedgeneJournal of Medical GeneticsVisit this articleSubmit a manuscriptReceive email alertsContact usBMJ
 
J Med Genet. 2007 May; 44(5): 314–321.
Published online 2007 January 12. doi:  10.1136/jmg.2006.046185
PMCID: PMC2597989

The human GIMAP5 gene has a common polyadenylation polymorphism increasing risk to systemic lupus erythematosus

Abstract

Background

Several members of the GIMAP gene family have been suggested as being involved in different aspects of the immune system in different species. Recently, a mutation in the GIMAP5 gene was shown to cause lymphopenia in a rat model of autoimmune insulin‐dependent diabetes. Thus it was hypothesised that genetic variation in GIMAP5 may be involved in susceptibility to other autoimmune disorders where lymphopenia is a key feature, such as systemic lupus erythematosus (SLE).

Material and methods

To investigate this, seven single nucleotide polymorphisms in GIMAP5 were analysed in five independent sets of family‐based SLE collections, containing more than 2000 samples.

Result

A significant increase in SLE risk associated with the most common GIMAP5 haplotype was found (OR 1.26, 95% CI 1.02 to 1.54, p = 0.0033). In families with probands diagnosed with trombocytopenia, the risk was increased (OR 2.11, 95% CI 1.09 to 4.09, p = 0.0153). The risk haplotype bears a polymorphic polyadenylation signal which alters the 3′ part of GIMAP5 mRNA by producing an inefficient polyadenylation signal. This results in higher proportion of non‐terminated mRNA for homozygous individuals (p<0.005), a mechanism shown to be causal in thalassaemias. To further assess the functional effect of the polymorphic polyadenylation signal in the risk haplotype, monocytes were treated with several cytokines affecting apoptosis. All the apoptotic cytokines induced GIMAP5 expression in two monocyte cell lines (1.5–6 times, p<0.0001 for all tests).

Conclusion

Taken together, the data suggest the role of GIMAP5 in the pathogenesis of SLE.

Keywords: genetic association, autoimmune, apoptosis, susceptibility gene

Systemic lupus erythematosus (SLE, OMIM 152700) is a multisystem autoimmune disease with largely unknown aetiology affecting predominantly women. Its prevalence in the Western population is about 0.05%,1,2 a number that differs widely in other ethnic groups.3 Although many candidate genes have been suggested as being involved in susceptibility to SLE,4,5 few have been replicable.6

SLE commonly shows familial co‐segregation with other autoimmune disorders: 10–20% of SLE probands have at least one first‐ or second‐degree relative with another autoimmune disorder.3 Disturbances of apoptosis of regulatory T‐cells, which normally facilitates maintenance of immune tolerance, has been suggested as causing autoimmunity7 and is thought to contribute to lymphopenia in several autoimmune disorders such as autoimmune insulin‐dependent diabetes,8 Sjogren's syndrome9 and SLE.1,10

Recently, GIMAP5 (OMIM 608086; previously called Immunity‐associated nucleotide 4‐like1, IAN4L1, fig 1A1A)) was found to cause lymphopenia and autoimmune insulin‐dependent diabetes in a rat model.11,12,13 Several members of the GIMAP gene family have been suggested as being involved in different aspects of the immune system in different species.11,12,14,15 These events appear to be associated with cell proliferation and apoptosis, suggesting a role for GIMAP proteins in control of cell death/survival.16 In line with this, several studies have shown that GIMAP5 normally protects cells from apoptosis and that loss of GIMAP5 function results in decreased mitochondrial integrity with lymphocyte apoptosis of primarily T cells.11,12,17,18,19,20

figure mg46185.f1
Figure 1 Schematic of the GIMAP5 gene. A, GIMAP5 spans 6286 base pairs (bp) and contains three exons; exon 1 is untranslated and exon 3 contains more than 94% of the coding sequence. The deduced 307 amino‐acid protein (34.8 kilodaltons) ...

Thus, based on current research, GIMAP5 is a potential candidate gene for autoimmune disorders where apoptosis and lymphopenia are key features, such as SLE. Thus, we set out to test association between SLE and GIMAP5 and designed our study to control as many common confounders of genetic association studies as possible.6

Material and methods

Subjects and samples

All subjects gave written informed consent for participation in genetic studies on SLE. The study protocols for all samples were reviewed and approved by local ethical committees, for the Finnish samples at University of Helsinki and at Karolinska Institutet, and for the UK samples at the Multi‐Centre Research Ethics Committee.

Finnish SLE samples

Collection of this sample has been described in detail previously.1 All patients with SLE were interviewed by the same physician, and the case records from the hospitals where the patients were treated were reviewed. All patients met the American College of Rheumatology criteria for the diagnosis of SLE.21 According to the prevalence of SLE in Finland,22 approximately 1200 of 1500 available patients were contacted during two phases of recruitment, accounting for more than 80% of all Finnish patients with SLE requiring hospital‐based treatment. Blood samples from patients and their relatives were obtained, comprising a total of 192 families; of which 86 were multiply affected by SLE (multiplex families); the rest were families with sporadic cases. Clinical characteristics of the patients included in this cohort are described in supplementary table 11 and have been described previously.1 Of these families, 76 parent‐affected offspring trios and 84 discordant sib pairs were included in association studies.

Table thumbnail
Table 1 Minor allele frequencies and nominal pointwise association in the different samples

British samples

In the UK collection, diagnosis of SLE was established by telephone interview, health questionnaire and details from clinical notes. All collected probands conformed to the American College of Rheumatology criteria for SLE21 and written consent was obtained from all participants. The study cohort comprised 378 parent‐affected offspring trios and 171 discordant sib pairs of European Caucasian descent together with 23 parent‐affected offspring trios of Indian–Asian descent. An additional 14 Indian–Asian discordant sib pairs, 16 parent‐affected offspring trios and 15 discordant sib pairs of Afro–Caribbean descent and 12 families with mixed ethnicity were excluded from the analysis because of their low information content. The families have all been described in detail previously23 and the clinical characteristics of the patients in the cohort are described in supplementary table 11.

Genotyping

First, we sequenced (see supplementary material for details) the entire GIMAP5 gene in one patient homozygous for all of the associated single nucleotide polymorphisms (SNPs) within GIMAP5 to identify and validate SNPs in the GIMAP5 region, but only identified previously reported SNPs.24 We chose 10 markers in the GIMAP5 gene for genotyping: rs3807383, rs759011, rs1046355, rs9657881, rs13229808, rs13229806, rs10361, rs6598, rs2286899 and rs2286898. Primers and probes for genotyping were designed using Primer3 (see electronic information).

The Finnish samples were genotyped using the SnuPe kit according to the manufacturer's instructions (Amersham Pharmacia Biotech, Uppsala, Sweden). Three of the tested assays, rs9657881, rs13229806 and rs13229808, failed to meet our quality criteria (as described below) and were subsequently excluded from further genotyping. We found the remaining seven SNPs to identify all common variation within GIMAP5 redundantly.

Subsequently, DNA from our UK replication samples were genotyped for the same seven markers using matrix‐assisted laser desorption/ionisation time‐of‐flight mass spectrometry according to the manufacturer's instructions (Sequenom Inc., San Diego, California, USA).25 PCR assays and associated extension reactions were designed using the SpectroDESIGNER software (Sequenom Inc.).

Using either genotyping method, genotypes for each SNP were confirmed by two independent scorers. For all data, PedCheck26 was used to detect Mendelian inconsistencies and markers were excluded if they yielded more than 10% errors. Hardy–Weinberg calculations were performed to ensure that each marker was in equilibrium (p>0.05). Success rate of all assays was required to be at least 85%. To ensure consistency between the two genotyping methods, three markers were regenotyped in the Finnish material with the matrix‐assisted laser desorption/ionisation time‐of‐flight mass spectrometry method (rs6598, rs2286899 and rs2286898). The concordance levels between the methods were greater than 98%. The concordance within the UK sample sets (repeated in two different labs on two different Sequenom platforms) was greater than 99%.

Expression analysis of GIMAP5

Peripheral blood mononuclear cells were isolated from fresh blood samples from nine patients with SLE and nine healthy controls from Finland using Ficoll‐Paque PLUS (Amersham Biosciences, Uppsala, Sweden). Total RNA was extracted using RNeasy Mini‐kit (Qiagen, Chatsworth, California, USA) as instructed by the manufacturer and 800 ng were reverse transcribed to cDNA using TaqMan Reverse transcription reagents (Applied Biosystems, Warrington, UK).

Primers were designed to cover a 80 bp fragment, spanning from exon 2 to exon 3, specific for the coding domain (359–438 bp in GIMAP5 mRNA; NM_018384) using Primer3 (electronic information). Quantitative real‐time PCR (TaqMan) was performed on an ABI PRISM 7700 Sequence Detection System according to the manufacturer's instructions (Applied Biosystems, see supplementary material for detailed information).

Allelic expression analysis of rs759011, rs1046355, rs10361, rs6598 and rs2286899

Genomic DNA from blood and cDNA from the same individuals as described above were sequenced (see supplementary material) for SNP detection and investigation of allelic expression levels.27

Northern Blot analysis of GIMAP5 and βactin

An 80 bp RT–PCR product from brain (Clontech, Palo Alto, California, USA), specific for the coding domain of GIMAP5 (359–438 bp in GIMAP5 mRNA; accession number NM_018384) was labelled with P32‐dCTP (Amersham Biosciences) by random priming (High Prime, Roche, Penzberg, Germany) and used as probe on Human Multiple Tissue Northern Blots (Ambion, Austin, Texas, USA) and Human Immune System MTN® Blot II (Clontech) with βactin as control probe.

RT–PCR on GIMAP5 and GAPDH

The same primers used for the GIMAP5 expression (TaqMan) analysis were used to study the expression of GIMAP5 at mRNA level. Total RNA (0.5 μg) from neuroblast BE (2)‐C cells, Jurkat E6‐1 T‐cells, 15 HL‐60 eosinophils, A549 lung carcinoma epithelial cells, U937 histocytic lymphoma monocytes, normal pooled human skin (BD Biosciences, Palo Alto, California, USA) and human placenta (Ambion) were converted to cDNA using SuperScript™ III Reverse Transcriptase reagents (Invitrogen). Amplification reactions were run in a total volume of 25 μl with 2 μl of cDNA, 0.8 μM of each primer, 0.2 mM dNTP, 1.5 mM MgCl2 and 0.03 U/μl of HotStarTaq (Qiagen) using the following programme; 95°C for 15 min followed by 40 cycles (45 s at 95°C, 30 s at 55°C, and 1 min 15 s at 72°C) and 10 min extension at 72°C and run on an agarose gel.

Regulation of GIMAP5 expression by SLE‐associated cytokines

Two monocytic cell lines, THP‐1 and U937, were plated on six‐well‐plates (1.4×106 cells per well) and grown overnight. The cells were depleted of serum overnight before stimulation with tumor necrosis factor alpha (TNF‐α, 10 ng/ml, Sigma, St Louis, Missouri, USA), transforming growth factor beta 1 (TGF‐β1, 10 ng/ml, Sigma), basic fibroblast growth factor (FGF2, 10 ng/ml, Sigma), interleukin‐1 beta (IL‐1β, 0.5 and 5 U/ml, Roche Molecular Biochemicals, Indiananpolis, Indiana, USA) and gamma‐interferon (IFN‐γ, 1 and 10 ng/ml, Sigma). Stimulation was allowed to proceed for 24 h after which the cells were lysed and total RNA was extracted using RNeasy Mini‐kit (Qiagen). Cells grown in serum‐free medium were used as controls. All experiments were performed at least twice, in triplicate wells.

Total RNA (0.5 μg) from cultured cells was reverse transcribed to cDNA (Invitrogen) and used as a template for TaqMan PCR. Reactions were performed with the 7500 Fast Real‐Time PCR system using standard protocols (Applied Biosystems).

Differential polyA RT–PCR on GIMAP5

One forward primer located 412 bp upstream (F) and two reverse primers located 132 bp upstream (R1) or 42 bp downstream (R2) of rs6598 were designed. Fresh peripheral blood mononuclear cells were collected from 17 individuals, and lymphocytes were immortalised by Epstein‐Barr virus transformation. RNA extraction and cDNA transcription were performed as above.

As positive control, cDNA from placenta (Ambion), U937 monocytes, A549 lung epithelial cells, neuroblasts BE (2)‐C cells, Jurkat E6‐1 T lymphocytes, normal pooled skin (BD Biosciences), pooled human fetal brain (Clontech) and pooled human leucocytes (Clontech) were used. Genotypes for rs6598 in all non‐pooled samples were obtained by genotyping, sequencing or both.

Amplification reactions were run in duplicates in 25 μl with 20 ng of cDNA, 1.6 μM of forward primer (Proligo, France) and 0.8 μM each reverse primer (Proligo), 0.2 mM dNTP, 1.5 mM MgCl2 and 0.075 U/μl of HotStarTaq (Qiagen). Reactions were first heated to 95°C for 15 min, and then subjected to 39 cycles (45 s at 95°C, 30 s at 59°C and 1 min 15 s at 72°C) followed by a 10 min extension at 72°C and run on an agarose gel.

The signal of the long and short band for the different genotypes was quantified by measuring the area of the band and the number of pixels in the area after maximising the intensity on the gel. A mean of the ratio of the long band over the short band from two measurements on each individual was used to adjust for laboratory variation in the signal amounts. The individuals homozygous for the rs6598 canonical polyA signal (n = 2) were combined with the individuals heterozygous for the rs6598 canonical polyA signal (n = 6) for statistical analysis and compared with the individuals homozygous for the rare AATAGA polyA signal (n = 9).

GIMAP5 protein analysis

We examined the expression of GIMAP5 protein (accession number Q96F15) in formalin‐fixed paraffin embedded tissue sections from healthy individuals, using exon‐3‐specific polyclonal antibodies (fig 11).). The antibodies were generated against the peptide LGREREGSFHSNDLF in two rabbits (Sigma Genosys, Suffolk, UK) and affinity purified to final concentrations of 0.4 mg/ml and 1.0 mg/ml, respectively.

The tissue sections were stained using the avidin–biotin–peroxidase complex technique (DakoCytomation Norden, Solna, Sweden). Diaminobenzidine (Vector Laboratories, Immunkemi F&D AB, Jarfalla, Sweden) was used as the chromogenic substrate and Mayer's haematoxylin as the counterstain (Merck, Darmstadt, Germany). Incubation with primary antibody, diluted to a final concentration of 0.02 mg/ml in phosphate‐buffered solution (PBS) containing 3% bovine serum albumin (BSA, Sigma‐Aldrich), was performed overnight at 4°C. Corresponding preimmune serum was used as a negative control.

To test the specificity of the polyclonal antibodies, a Western blot analysis was performed. U937 monocytes were grown to a density of 1×106 cells per milliltre, washed twice in ice‐cold PBS and harvested in 1 ml RIPA buffer. Aliquots from cleared lysates containing 10, 20 or 30 μg protein were subjected to electrophoresis on a 10% SoftGel precast polyacrylamide gel (Mirador DNA Design INC, Montreal, Canada) and transferred to a Hybond‐P PVDF transfer membrane (Amersham Biosciences). Protein was detected using an ECL Advance Western Blotting Detection Kit (Amersham Biosciences), and then visualised on an ECL HyperfilmTM (Amersham Biosciences).

Statistics

To evaluate haplotype block patterns over the GIMAP5 locus, linkage disequilibrium between the SNPs was analysed using Haploview 2.128 Data for 119 SNPs (between rs952825 and rs717658) from the International HapMap project spanning the GIMAP5 locus was used as a reference. Two of our seven SNPs (rs759011 and rs2286898) were also genotyped in the Centre d'Etude du Polymorphisme Humain (CEPH) reference material from the HapMap data (greater than 99% concordance with genotyping in our laboratory).

Haplotypes were reconstructed using maximum‐likelihood estimation in nuclear families by FAMHAP.29 Subsequent family‐based, nominal pointwise and haplotype associations were first assessed by the transmission disequilibrium test30 as implemented in Genehunter 2.1 R2.31 Global significance of the distribution of haplotypes was assessed by a randomisation test on the transmitted alleles/haplotypes using 50 000 iterations. For association analysis using discordant sib pairs, we used the discordant allele test.32 Gender and age of onset as covariates were non‐significant and therefore disregarded from further analyses.

Meta‐analysis33 of the independent sample sets was performed in two steps, first for SLE and then for thrombocytopenia; first the Breslow–Day test for non‐compatibility of ORs was performed (p>0.2) and then the pooled estimate of ORs was obtained using the Mantel–Haenzel method as implemented in the “R” software (see electronic information).

The differences in expression levels between individuals with SLE compared with healthy controls as well as by haplotype, independent of affection status, were calculated as ΔCt values between (GIMAP5−GAPDH) and subsequently compared using a t test as implemented in StatsDirect software version 2.3.4.

The differences in expression levels between untreated and cytokine treated monocytes were calculated as ΔCt values between (GIMAP5−GAPDH) and subsequently compared using a t test as implemented in StatsDirect version 2.3.4. Fold changes in GIMAP5 expression levels were calculated according to the manufacturer's instructions (Applied Biosystems).

The differences in the signal amount of the long and short band between the different genotypes in the polyA differential RT–PCR was assessed by t test as implemented in StatsDirect version 2.3.4.

All p values are given as nominal, two sided p values if not stated otherwise.

Results

To capture all common variation in GIMAP5, we genotyped all individuals using seven SNPs distributed over the gene (fig 1A1A).). All markers except one (rs3807383) belong to one haplotype block (as defined by Gabriel et al28), covering the protein coding parts as well as the 3′ untranslated region (UTR) (fig 1B1B).). GIMAP5 spans a 360 kilobase (kb) region on chromosome 7q36.1, also containing a cluster of 10 additional GIMAP gene family members. To ensure that the genetic effect was specific to GIMAP5, we evaluated the linkage disequilibrium structure in the region, and GIMAP5 is clearly separated from the nearby genes.

In Finnish trios we found suggestive association of SLE to a single SNP (table 11;; rs1046355 p[less-than-or-eq, slant]0.06) and to the most common haplotype (table 22;; rs759011–rs2286898, CCGGTC, p[less-than-or-eq, slant]0.05). To evaluate these findings more extensively,6 we genotyped the same seven SNPs in four additional family‐based sample sets and tested for GIMAP5 association with SLE. In 84 discordant Finnish sib pairs, we found nominal significance for two markers (rs759011 p<0.03 and rs2286898 p<0.04). The suggested risk haplotype showed a trend in the same direction as the initial finding (OR 1.09, 95% CI 0.68 to 1.73, p = 0.13). In 378 UK Caucasian trios, the most common haplotype showed a trend towards an association with SLE (OR 1.25, 95% CI 0.96 to 1.64, p = 0.1, fig 22).

Table thumbnail
Table 2 Analysis of the common six marker haplotypes spanning GIMAP5 in five independent sample sets
figure mg46185.f2
Figure 2 Meta‐analysis of the association between the GIMAP5 risk haplotype and SLE (A) and the association of the GIMAP5 risk haplotype thrombocytopenia in patients with SLE (B). In both cases the OR of each analysis surrounded by the ...

In 21 UK Asian trios, two SNPs suggested pointwise association (rs759011; p = 0.039 and rs6598 p = 0.033) and the common risk‐haplotype was trending towards an association with SLE (OR 3.67, 95% CI 1.02 to 13.2, p = 0.03). However, the 95% CI for the genetic effect in this population was broad, reflecting the small sample size (fig 22).). In an additional UK sample, 171 discordant sib pairs were tested for association but no point wise or haplotype associations were detected (OR 0.92, 95% CI 0.69 to 1.24, p = 0.27).

The most exact estimate of the suggested haplotype association was obtained by combining the ORs in a Mantel–Haenzel test. This analysis reveals a significant association between the most common GIMAP5 haplotype and SLE (OR 1.26, 95% CI 1.02 to 1.54, p = 0.0033, one sided, fig 22).

GIMAP5 is a small protein, suggested to reside in the mitochondrial membrane19,34,35 (fig 1A1A).). None of the genotyped SNPs of the risk haplotype changed the amino‐acid sequence of GIMAP5; two were silent polymorphisms, one altered the first canonical polyadenylation (polyA) signal (AATAAA) and one was a putative splice variant 22 bp from the end of the mRNA occurring after the canonical polyA signal (fig 1A1A).). We therefore compared total GIMAP5 mRNA expression level between patients with SLE (n = 8) and controls (n = 9), or by their individual haplotype, independent of affection status. We saw a nearly significant difference in expression between patients with SLE and controls (fig 4C4C,, p = 0.074, power for α<0.05 significance is 76%), with lower expression in patients. No significant difference was seen between non‐risk homozygous + risk‐heterozygous (n = 9) and risk‐homozygous (n = 10) (p = 0.7, data not shown). We also assessed the allele‐specific expression of all genotyped SNPs in exon 3 (rs759011, rs1046355, rs10361, rs6598, rs2286899), but did not find any significant differences between patients and controls.

figure mg46185.f4
Figure 4 Cytokine effect on GIMAP5 mRNA expression in (A) THP‐1 monocytes and (B) U937 monocytes. Data were collected for untreated cells (controls) and cells treated with five different cytokines. All cytokines significantly upregulated ...

GIMAP5 contains two consecutive polyA signals, and one SNP (rs6598) that marks the risk haplotype resides at the fifth base of the first canonical polyA signal (fig 1A1A).). This SNP causes a switch from the canonical AATAAA polyA signal to a rare polyA signal (AATAGA), present in only 0.7% of mature mRNA species.36 We hypothesised that the alternative polyA signal might cause a differential termination of the transcription of GIMAP5, an mechanism that has previously been described for the globulin gene causing α‐ and βthalassaemia.37,38 Thus, we investigated differences in mRNA lengths between individuals with different genotypes at rs6598. Individuals homozygous for the rarer AATAGA signal have a higher amount of longer mRNA transcripts than individuals heterozygous or homozygous for the conservative polyA signal AATAAA (comparing the ratio of long bands/short bands in individuals, n = 9, with either AA/AG genotypes = 0.22±0.10 or individuals, n = 8, with GG genotypes = 0.37±0.16, p<0.005, power for α<0.05 significance is 76%). The AATAGA polyA signal efficiency has been suggested as cellular‐context dependent,36 and there is emerging evidence that different polyA signals in the same gene could react differently on environmental factors in a tissue, differentiation or development specific way.36,39 Supporting this, we observed either one (suggesting predominant use of the first polyA signal) or two bands (suggesting leaky recognition and use of both polyA signals) in our polyA RT–PCR when using cDNA from placenta, U937 monocytes, A549 lung epithelial cells, neuroblasts BE (2)‐C cells, Jurkat E6‐1 T lymphocytes, normal pooled skin, pooled human fetal brain and pooled human leucocytes (fig 1D1D),), with all non‐pooled cDNA samples being homozygous GG for rs6598, and thereby also homozygous for the rarer AATAGA polyA signal. The genotype for the pooled samples could not be determined because of mixed genotypes in these samples.

Next, we investigated whether the identified GIMAP5 risk haplotype might affect the haematological profile of patients with SLE, as could be expected from the initial results from the animal model of autoimmune type 1 diabetes.11,12 As a systemic, autoimmune disorder, SLE typically displays several different haematological disturbances, such as thrombocytopenia or lymphopenia, in a large proportion of patients (20–40%).40,41 Thrombocytopenia was available as the only haematological phenotype in all our patient cohorts. We analysed our data stratified by the presence of thrombocytopenia and found evidence that in patients with SLE, carriage of the risk haplotype increased the likelihood of developing thrombocytopenia (OR 2.11, 95%, CI 1.09 to 4.09, p = 0.0153). GIMAP5 mRNA was strongly expressed in the spleen (the major site of thrombocyte clearance), lymph nodes and lung, but less abundantly in other tissues (fig 33).

figure mg46185.f3
Figure 3 Multiple tissues (A) and immune‐related tissue northern blot analysis (B). Commercial northern blots containing 2 μg polyA+ RNA from brain (1), placenta (2), skeletal muscle (3), heart (4), kidney (5), ...

To characterise the human GIMAP5 protein expression in humans more precisely, we studied relevant tissues using immunohistochemistry. Interestingly, GIMAP5 was detected in endothelial cells in blood vessels and kidney, lung septal cells, parenchymal cells in the thymus, and polymorphonuclear leucocytes, in agreement with mRNA expression patterns (supplementary fig 22).

To investigate what effect several cytokines (IFN‐γ, IL‐1β, TNF‐α, TGF‐β and FGF2), well‐known to affect apoptosis and inflammation, would have on GIMAP5 expression, we treated two different monocyte cell lines (U937 and THP‐1) with these cytokines and measured changes in mRNA levels. We found consistent and significant upregulation of GIMAP5 expression in both monocyte cell lines (fig 4A,B4A,B;; IFN‐γ 5.7/6.7‐fold, p values 10−8/10−5; IL‐1β 3.9/4.9‐fold, p values 10−7/10−6, TNF‐α 2.1/2.2‐fold, p values 10−4/10−4, TGF‐β 1.4/2.4‐fold, p values 0.05/10−7 and FGF2 1.5/1.8‐fold, p values 0.02/10−6 ; power for α<0.05 significance is 80%) after cytokine treatment.

Discussion

An association between SLE and GIMAP5 is consistent in all independent family‐based samples (combined OR 1.26, 95% CI 1.02 to .54, p = 0.0033). The suggested risk haplotype is common (about 50% in control chromosomes) as seen in several associations to common complex disorders.42,43 This may explain why no genome‐wide scans have identified this region as a significant candidate for harbouring an SLE‐susceptibility gene,4 because most genome‐wide studies are underpowered to detect common risk variants.44 However, GIMAP5 is located about 2.5 centimorgans away from the suggestive scores of linkage in the combined mapping results for Minnesota SLE cohorts 1 and 2,45 and the region on chromosome 7q34–36.3 is ranked as one of the most interesting regions in a recent meta‐analysis of all SLE genome‐wide linkage studies.5

GIMAP5 was associated to thrombocytopenia in patients with SLE with the risk haplotype. Because thrombocytopenia predicts severe disease progression and premature death in patients with SLE, this finding may potentially have clinical relevance. It would be important to investigate the GIMAP5 association with the other haematological features of SLE, but larger collaborations will be needed for sufficient power for the less common subphenotypes.

The expression of GIMAP5 was not significantly altered, but the detected trend towards lower levels of mRNA in patients is in line with findings of increased leucocyte apoptosis in patients with SLE and findings of reduced T‐lymphocyte survival in a mouse knockdown for GIMAP5.46 The two consecutive polyA signals in GIMAP5 are a rather unusual feature in human‐expressed sequence tags.36 A similar A→G replacement of the fifth base of the polyA signal as rs6598 is observed in the human globin mRNA and is known to cause α‐ and βthalassaemia.37,38 The generated globin polyA signal is inefficient and allows transcription to continue past it, producing a high proportion of non‐terminated, long globin mRNA causing thalassaemia.37,38 Our data suggested similarly that the shorter GIMAP5 transcript is the predominant mRNA species, but that the A→G substitution causes an inefficient polyA site and results in a longer transcript, as hypothesised and shown previously.37,38 Our data also support the idea that different polyA sites in the same gene might react differently on environmental factors in a tissue‐, differentiation‐ or development‐specific way.36,39 Thus, taken together, we suggest that the alternative polyA sites created by the disease‐associated, risk haplotype‐tagging SNP may be functionally relevant for disease susceptibility in SLE.

An in vivo role of GIMAP5 is supported by recent observations in the BB rat. In this animal model, a frameshift mutation creates a stop codon that contributes to autoimmune diabetes type 1 with lymphopenia.11,12 The actual effect of the mutation on the immune cells of the BB rat remains unsolved, but the strongly decreased expression of mRNA for rGimap511 suggest that the cells either lack the rGimap5 protein completely or produce a truncated protein. Also, using transgenic rescue in lymphopenic BB‐derived congenic F344.lyp/lyp rats, Michalkiewicz et al have demonstrated that the Gimap5 protein is essential for maintaining normal T‐cell levels.13

Sandal et al suggested that the antiapoptotic action of rGimap5 is executed through its central coiled‐coil‐containing domain, responsible for centrosomal anchoring and resistance towards death, and that if expressed in the BB rat, the mutated rGimap5 will be nonfunctional for this apoptosis antagonism.20 They show that Jurkat T‐cells transfected with antisense‐GIMAP5 cDNA became hypersensitive to general apoptosis inducers (γ‐radiation and okadaic acid), suggesting that the endogenous level of GIMAP5 in Jurkat cells normally protect against cell death.20 The members of the GIMAP family of antiapoptotic proteins have been reported to protect against a broad spectrum of apoptogens (including anti‐CD95 (Fas/Apo‐1), TNF‐α, staurosporine and etoposide),47 but the findings of Sandal et al contrast this.20

We examined what effect several pleiotropic cytokines commonly implicated in SLE pathogenesis had on GIMAP5 expression and saw a consistent, significant upregulation of GIMAP5 in two monocyte cell lines. IFN‐γ, IL‐1β and TNF‐α are all proinflammatory cytokines that are suggested to play a part in apoptosis.48,49 TGF‐β and FGF2 are growth factors and disturbances in their activation and signalling may result in apoptosis. This supports the notion that GIMAP5 is affected by cytokines important in the immune system and for the development of autoimmune disease, and is upregulated in response to cytokines known to be elevated in SLE, perhaps to protect against apoptosis. Our identified, potentially functional polyA variation is more common, not fully penetrant and would subsequently be expected to have a milder effect on the GIMAP5 antiapoptotic effects than the mutation described in the BB rat. It is also possible that the genetic variation in GIMAP5 predisposing to SLE susceptibility is regulating its action in a tissue‐ or cell‐specific manner.

In summary, by analysing five independent family base samples with over 2000 individuals, we have found a significant (OR 1.26, p = 0.0033) increase in the risk of SLE associated with a common GIMAP5 risk haplotype. We also confirm a functional variant on the SLE risk haplotype, which introduces a non‐optimal polyA site as a possible susceptibility mechanism. Our data imply that inherited variation in GIMAP5 is a susceptibility factor for SLE.

Electronic database information

HapMap, public data release #11 September 2004 (http://www.hapmap.org/); NCBI (http://www.ncbi.nlm.nih.gov/); Online Mendelian Inheritance in Man (OMIM) (http://www.ncbi.nih.gov/entrez/query.fcgi?db = OMIM); Primer3 (http://frodo.wi.mit.edu/primer3/primer3_code.html); Saccharomyces Genome Database (http://seq.yeastgenome.org/cgi‐bin/SGD/web‐primer); “R” (http://www.r‐project.org/).

Acknowledgements

We thank Linda Berglind and Ville‐Veikko Mäkelä for their contributions in genotyping the samples. We thank all patients and family members participating in the Finnish and UK SLE studies. This study was supported by the Sigrid Jusélius Foundation, the Academy of Finland and the Swedish Research Council (to JK and US‐K); the Ulla & Gustaf af Ugglas Foundation and the Alex & Eva Wallströms Foundation and travel grants from the Royal Physiographic Society and the Novartis Foundation (to CML); the Edvard Welanders Foundation and the Finsen Foundation (to TS); and a Wellcome Trust Senior Fellowship (to TJV, supporting AW and DCG).

Abbreviations

BSA - bovine serum albumin

PBS - phosphate‐buffered solution

RT–PCR - polymerase chain reaction with reverse transcription

SLE - systemic lupus erythematosus

Footnotes

Competing interests: None.

References

1. Koskenmies S, Widen E, Kere J. et al Familial systemic lupus erythematosus in Finland. J Rheumatol 2001. 28758–760.760 [PubMed]
2. Prokunina L, Castillejo‐Lopez C, Oberg F. et al A regulatory polymorphism in PDCD1 is associated with susceptibility to systemic lupus erythematosus in humans. Nat Genet 2002. 32666–669.669 [PubMed]
3. Nath S K, Kelly J A, Harley J B. et al Mapping the systematic lupus erythematosus susceptibility genes. Methods Mol Med 2004. 10211–29.29 [PubMed]
4. Tsao B P. Update on human systemic lupus erythematosus genetics. Curr Opin Rheumatol 2004. 16513–521.521 [PubMed]
5. Lee Y H, Nath S K. Systemic lupus erythematosus susceptibility loci defined by genome scan meta‐analysis. Hum Genet 2005. 118434–443.443 [PubMed]
6. Hirschhorn J N, Altshuler D. Once and again‐issues surrounding replication in genetic association studies. J Clin Endocrinol Metab 2002. 874438–4441.4441 [PubMed]
7. Marleau A M, Sarvetnick N. T cell homeostasis in tolerance and immunity. J Leukoc Biol 2005. 78575–584.584 [PubMed]
8. Kaaba S A, Al‐Harbi S A. Abnormal lymphocyte subsets in Kuwaiti patients with type‐1 insulin‐dependent diabetes mellitus and their first‐degree relatives. Immunol Lett 1995. 47209–213.213 [PubMed]
9. Mandl T, Bredberg A, Jacobsson L T. et al CD4+ T‐lymphocytopenia–a frequent finding in anti‐SSA antibody seropositive patients with primary Sjogren's syndrome. J Rheumatol 2004. 31726–728.728 [PubMed]
10. Mirzayan M J, Schmidt R E, Witte T. Prognostic parameters for flare in systemic lupus erythematosus. Rheumatology (Oxford) 2000. 391316–1319.1319 [PubMed]
11. Hornum L, Romer J, Markholst H. The diabetes‐prone BB rat carries a frameshift mutation in Ian4, a positional candidate of Iddm1. Diabetes 2002. 511972–1979.1979 [PubMed]
12. MacMurray A J, Moralejo D H, Kwitek A E. et al Lymphopenia in the BB rat model of type 1 diabetes is due to a mutation in a novel immune‐associated nucleotide (Ian)‐related gene. Genome Res 2002. 121029–1039.1039 [PubMed]
13. Michalkiewicz M, Michalkiewicz T, Ettinger R A. et al Transgenic rescue demonstrates involvement of the Ian5 gene in T cell development in the rat. Physiol Genomics 2004. 19228–232.232 [PubMed]
14. Reuber T L, Ausubel F M. Isolation of Arabidopsis genes that differentiate between resistance responses mediated by the RPS2 and RPM1 disease resistance genes. Plant Cell 1996. 8241–249.249 [PubMed]
15. Krucken J, Schmitt‐Wrede H P, Markmann‐Mulisch U. et al Novel gene expressed in spleen cells mediating acquired testosterone‐resistant immunity to Plasmodium chabaudi malaria. Biochem Biophys Res Commun 1997. 230167–170.170 [PubMed]
16. Stamm O, Krucken J, Schmitt‐Wrede H P. et al Human ortholog to mouse gene imap38 encoding an ER‐localizable G‐protein belongs to a gene family clustered on chromosome 7q32–36. Gene 2002. 282159–167.167 [PubMed]
17. Cambot M, Aresta S, Kahn‐Perles B. et al Human immune associated nucleotide 1: a member of a new guanosine triphosphatase family expressed in resting T and B cells. Blood 2002. 993293–3301.3301 [PubMed]
18. Zenz T, Roessner A, Thomas A. et al hIan5: the human ortholog to the rat Ian4/Iddm1/lyp is a new member of the Ian family that is overexpressed in B‐cell lymphoid malignancies. Genes Immun 2004. 5109–116.116 [PubMed]
19. Pandarpurkar M, Wilson‐Fritch L, Corvera S. et al Ian4 is required for mitochondrial integrity and T cell survival. Proc Natl Acad Sci U S A 2003. 10010382–10387.10387 [PubMed]
20. Sandal T, Aumo L, Hedin L. et al Irod/Ian5: an inhibitor of gamma‐radiation‐ and okadaic acid‐induced apoptosis. Mol Biol Cell 2003. 143292–3304.3304 [PMC free article] [PubMed]
21. Tan E M, Cohen A S, Fries J F. et al The 1982 revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheum 1982. 251271–1277.1277 [PubMed]
22. Helve T. Prevalence and mortality rates of systemic lupus erythematosus and causes of death in SLE patients in Finland. Scand J Rheumatol 1985. 1443–46.46 [PubMed]
23. Russell A I, Cunninghame Graham D S, Shepherd C. et al Polymorphism at the C‐reactive protein locus influences gene expression and predisposes to systemic lupus erythematosus. Hum Mol Genet 2004. 13137–147.147 [PubMed]
24. Payne F, Smyth D J, Pask R. et al Haplotype tag single nucleotide polymorphism analysis of the human orthologues of the rat type 1 diabetes genes Ian4 (Lyp/Iddm1) and Cblb. Diabetes 2004. 53505–509.509 [PubMed]
25. Jurinke C, van den Boom D, Cantor C R. et al Automated genotyping using the DNA MassArray technology. Methods Mol Biol 2002. 187179–192.192 [PubMed]
26. O'Connell J R, Weeks D E. PedCheck: a program for identification of genotype incompatibilities in linkage analysis. Am J Hum Genet 1998. 63259–266.266 [PubMed]
27. Pastinen T, Sladek R, Gurd S. et al A survey of genetic and epigenetic variation affecting human gene expression. Physiol Genomics 2004. 16184–193.193 [PubMed]
28. Gabriel S B, Schaffner S F, Nguyen H. et al The structure of haplotype blocks in the human genome. Science 2002. 2962225–2229.2229 [PubMed]
29. Becker T, Knapp M. Maximum‐likelihood estimation of haplotype frequencies in nuclear families. Genet Epidemiol 2004. 2721–32.32 [PubMed]
30. Spielman R S, McGinnis R E, Ewens W J. Transmission test for linkage disequilibrium: the insulin gene region and insulin‐dependent diabetes mellitus (IDDM). Am J Hum Genet 1993. 52506–516.516 [PubMed]
31. Kruglyak L, Daly M J, Reeve‐Daly M P. et al Parametric and nonparametric linkage analysis: a unified multipoint approach. Am J Hum Genet 1996. 581347–1363.1363 [PubMed]
32. Boehnke M, Langefeld C D. Genetic association mapping based on discordant sib pairs: the discordant‐alleles test. Am J Hum Genet 1998. 62950–961.961 [PubMed]
33. Lohmueller K E, Pearce C L, Pike M. et al Meta‐analysis of genetic association studies supports a contribution of common variants to susceptibility to common disease. Nat Genet 2003. 33177–182.182 [PubMed]
34. Andersen U N, Markholst H, Hornum L. The antiapoptotic gene Ian4l1 in the rat: genomic organization and promoter characterization. Gene 2004. 341141–148.148 [PubMed]
35. Krucken J, Schroetel R M, Muller I U. et al Comparative analysis of the human gimap gene cluster encoding a novel GTPase family. Gene 2004. 341291–304.304 [PubMed]
36. Beaudoing E, Gautheret D. Identification of alternate polyadenylation sites and analysis of their tissue distribution using EST data. Genome Res 2001. 111520–1526.1526 [PubMed]
37. van Solinge W W, Lind B, van Wijk R. et al Clinical expression of a rare beta‐globin gene mutation co‐inherited with haemoglobin E‐disease. Eur J Clin Chem Clin Biochem 1996. 34949–954.954 [PubMed]
38. Jankovic L, Efremov G D, Petkov G. et al Two novel polyadenylation mutations leading to beta(+)‐thalassemia. Br J Haematol 1990. 75122–126.126 [PubMed]
39. Edwalds‐Gilbert G, Veraldi K L, Milcarek C. Alternative poly(A) site selection in complex transcription units: means to an end? Nucleic Acids Res 1997. 252547–2561.2561 [PMC free article] [PubMed]
40. Hochberg M C. Updating the American College of Rheumatology revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheum 1997. 401725 [PubMed]
41. Kao A H, Manzi S, Ramsey‐Goldman R. Review of ACR hematologic criteria in systemic lupus erythematosus. Lupus 2004. 13865–868.868 [PubMed]
42. Corder E H, Saunders A M, Strittmatter W J. et al Gene dose of apolipoprotein E type 4 allele and the risk of Alzheimer's disease in late onset families. Science 1993. 261921–923.923 [PubMed]
43. Altshuler D, Hirschhorn J N, Klannemark M. et al The common PPARgamma Pro12Ala polymorphism is associated with decreased risk of type 2 diabetes. Nat Genet 2000. 2676–80.80 [PubMed]
44. Hirschhorn J N, Daly M J. Genome‐wide association studies for common diseases and complex traits. Nat Rev Genet 2005. 695–108.108 [PubMed]
45. Gaffney P M, Ortmann W A, Selby S A. et al Genome screening in human systemic lupus erythematosus: results from a second Minnesota cohort and combined analyses of 187 sib‐pair families. Am J Hum Genet 2000. 66547–556.556 [PubMed]
46. Baechler E C, Gregersen P K, Behrens T W. The emerging role of interferon in human systemic lupus erythematosus. Curr Opin Immunol 2004. 16801–807.807 [PubMed]
47. Deveraux Q L, Reed J C. IAP family proteins–suppressors of apoptosis. Genes Dev 1999. 13239–252.252 [PubMed]
48. Firme L, Bush A B. FGF signaling inhibits the proliferation of human myeloma cells and reduces c‐myc expression. BMC Cell Biol 2003. 417 [PMC free article] [PubMed]
49. Davidson A, Diamond B. Autoimmune diseases. N Engl J Med 2001. 345340–350.350 [PubMed]

Articles from Journal of Medical Genetics are provided here courtesy of BMJ Group