Article Text

Download PDFPDF

The human GIMAP5 gene has a common polyadenylation polymorphism increasing risk to systemic lupus erythematosus
  1. Anna Hellquist1,2,*,
  2. Marco Zucchelli1,*,
  3. Katja Kivinen1,
  4. Ulpu Saarialho-Kere3,4,
  5. Sari Koskenmies3,5,
  6. Elisabeth Widen6,
  7. Heikki Julkunen7,
  8. Andrew Wong8,
  9. Marja-Liisa Karjalainen-Lindsberg9,
  10. Tiina Skoog1,2,
  11. Johanna Vendelin5,
  12. Deborah S Cunninghame-Graham8,
  13. Timothy J Vyse8,
  14. Juha Kere1,2,5,
  15. Cecilia M Lindgren1
  1. 1Department of Biosciences at Novum, Karolinska Institute, Stockholm, Sweden
  2. 2Clinical Research Centre, Karolinska University Hospital, Huddinge, Sweden
  3. 3Department of Dermatology, University of Helsinki, Helsinki, Finland
  4. 4Department of Dermatology, Karolinska Institutet at Stockholm Söder Hospital, Stockholm, Sweden
  5. 5Department of Medical Genetics, University of Helsinki, Finland
  6. 6Finnish Genome Center, University of Helsinki, Helsinki, Finland
  7. 7Department of Medicine, Peijas Hospital and Helsinki University Hospital, Helsinki, Finland
  8. 8Rheumatology Section, Imperial College, Hammersmith Hospital, London, UK
  9. 9Department of Pathology, Helsinki University Central Hospital, Helsinki, Finland
  1. Correspondence to:
 Professor Juha Kere
 Karolinska Institutet, Department of Biosciences and Nutrition, S-141 57 Huddinge, Sweden; juha.kere{at}biosci.ki.se

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.

  • BSA, bovine serum albumin
  • PBS, phosphate-buffered solution
  • RT–PCR, polymerase chain reaction with reverse transcription
  • SLE, systemic lupus erythematosus
  • genetic association
  • autoimmune
  • apoptosis
  • susceptibility gene
View Full Text

Statistics from Altmetric.com

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 1A) was found to cause lymphopenia and autoimmune insulin-dependent diabetes in a rat model.11–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–20

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) contains five motifs conserved in guanosine triphosphate (GTP)-binding proteins and a putative C-terminal transmembrane region.19 B, representation of LD structure in the GIMAP5 gene using D′ with 95% CI.28 Each square depicts the magnitude of linkage disequilibrium (LD) for a pair of markers, red colour indicating high D′ and white colour low D′. C, The difference in ratio of amounts between the longer (526 bp) and the shorter (457 bp) bands compared between individuals either homozygous or heterozygous for the common AATAAA polyA signal and individuals homozygous for the AATAGA signal in the 3′ untranslated region (UTR) of GIMAP5 (p<0.005). D, Positive complementary DNA (cDNA) controls for the polyA polymerase chain reaction with reverse transcription (RT–PCR) showing a difference in bands independent of genotype (all non-pooled samples are homozygous GG for rs6598 and are thereby homozygous for the rare AATAGA polyA signal) in different tissues, suggesting that the GIMAP5 polyA signal functions in a cell-type-specific manner; (1) placenta, (2) U937 monocytes, (3) A549 lung epithelial cells, (4) neuroblasts BE (2)-C cells, (5) Jurkat E6-1 T lymphocytes, (6) normal pooled skin, (7) pooled brain and (8) pooled leucocytes.

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 1 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 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 1.

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 1). 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 1A). 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 1B). 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 1; rs1046355 p⩽0.06) and to the most common haplotype (table 2; rs759011–rs2286898, CCGGTC, p⩽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 2).

Table 2

 Analysis of the common six marker haplotypes spanning GIMAP5 in five independent sample sets

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 95% CI for the six-marker haplotype of GIMAP5 is given. For each study, • represents the estimated OR for the risk haplotype. Bars indicate the 95% CIs.

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 2). 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 2).

GIMAP5 is a small protein, suggested to reside in the mitochondrial membrane19,34,35 (fig 1A). 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 1A). 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 4C, 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.

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 1A). 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 1D), 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 3).

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), pancreas (6), liver (7), lung (8), spleen (9), colon (10), spleen (11), lymph node (12), thymus (13), peripheral blood leucocytes (14), bone marrow (15) and fetal liver (16) were used. After exposing the blot with a 32P-random labelled 80 bp probe for 72 hours, we detected a 1.6 kilobase (kb) transcript and after additional 72 h a weak band of approximately 5.5 kb appears, as reported previously.35 For control, the northern blots were probed with βactin. C, RT-PCR on GIMAP5 and GAPDH on cDNA from neuroblasts BE (2)-C cells (1), Jurkat E6-1 T lymphocytes (2), 15 HL-60 eosinophils (3), A549 lung epithelial cells (4), U937 monocytes (5), normal pooled skin (6), placenta (7).

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 2).

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,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.

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 in a similar fashion in both tested cell lines (IFN-γ 5.7/6.7-fold, p values 10−8/10−5; IL-1β 4.9/3.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). C, Differences in expression between patients with SLE and controls shown in log scale (p = 0.074, power for α<0.05 significance is 76%), showing a trend for lower expression in patients. Expression levels are given as relative fold changes of GIMAP5 (y axis).

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/).

Acknowledgments

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).

REFERENCES

View Abstract

Supplementary materials

Footnotes

  • * These authors contributed equally to the manuscript

  • Published Online First 12 January 2007

  • Competing interests: None.

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.