ORIGINAL RESEARCH article

Front. Neurol., 10 July 2018

Sec. Neuromuscular Disorders and Peripheral Neuropathies

Volume 9 - 2018 | https://doi.org/10.3389/fneur.2018.00529

IL-4Rα Polymorphism Is Associated With Myasthenia Gravis in Chinese Han Population

  • 1. Department of Neurology, Affiliated Hospital of Qingdao University, Qingdao, China

  • 2. Department of Neurology, Qilu Hospital of Shandong University, Jinan, China

  • 3. Department of Clinical Medicine, University of Bergen, Bergen, Norway

  • 4. Department of Neurology, Beijing Friendship Hospital, Capital Medical University, Beijing, China

  • 5. ICU, Affiliated Hospital of Qingdao University, Qingdao, China

  • 6. Department of Neurology, Peking University First Hospital, Beijing, China

Abstract

Interleukin-4 (IL-4) is a potent growth and differentiation factor for B cells which play a vital role in the pathogenesis of myasthenia gravis (MG). IL-4 exerts its function by binding to three types of IL-4 receptor (IL-4R) complexes. IL-4Rα is the key component of the IL-4R complex. We hypothesize that polymorphism of IL-4Rα gene may be associated with the susceptibility and severity of MG. A Chinese cohort of 480 MG patients and 487 healthy controls were recruited. Polymorphisms of IL-4Rα gene were determined with SNPscan™ methods and compared between MG and control groups, as well as among MG subgroups. Rs2107356 and rs1805010 were found to be associated with adult thymoma associated MG, and rs1801275 was found to be associated with adult non-thymoma AChR-Ab positive MG. We did not found association between IL-4Rα polymorphism and the severity of MG. Genetic variations of IL-4Rα were found associated with the susceptibility of MG in Chinese Han population.

Introduction

Myasthenia gravis (MG) is a T cell-dependent antibody-mediated autoimmune neuromuscular disorder and characterized by weakness and fatigability of skeletal muscles (1). Interleukin-4 (IL-4) is one of the Th2 cytokines which plays a vital role in the pathogenesis of MG. IL-4 is a potent growth and differentiation factor for B cells and stimulates class-switching and autoantibody production. IL-4 also takes part in the positive modulation of regulatory T cell (Treg) functions (2). IL-4 levels were found lower in MG patients compared with healthy controls (3). In mice genetically deficient in IL-4, long-lasting muscle weakness developed after a single immunization with acetylcholine receptor (AChR) (4).

IL-4 binds to three types of IL-4 receptor (IL-4R) complexes to exert its function. The type I IL-4R complex is present on lymphoid T cells, NK cells, basophils and mast cells; type II complex is found on non-lymphoid tumor cells, and type III complex is present on B cells and monocytes. Type I IL-4R is composed of two subunits, the IL-4Rα and the IL-2Rγc-common (γC) chain. Type II IL-4R is composed of two subunits, the IL-4Rα and the IL-13Rα1. In type III IL-4R, all three chains are present. Hence, IL-4Rα is the key component of all three complexes. When IL-4 binds to IL-4Rα, it recruits one of the other two chains in the complexes depending on the cell type and forms a heterodimer complex that initiates signal transduction (5). In experimental MG, IL-4R is upregulated in the myocytes following exposure to anti-AChR antibody, and demonstrated an increased responsiveness to IL-4, which implies the role of IL-4R in disease progression of experimental MG (6).

Till now, there was no evidence on the association of IL-4 gene polymorphism with MG (7). Polymorphisms of IL-4Rα were found to be associated with the susceptibility of autoimmune or inflammatory diseases such as rheumatoid arthritis (810), systemic lupus erythematosus (11), Graves' disease (12), eczema (13) and asthma (1416), and with the severity and treatment effects of some of the diseases (9, 14). Pal et al. reported the frequency of rs1805010 GG genotype was significantly higher in MG patients than in healthy controls (17). The association between the polymorphism of IL-4Rα and MG remains to be confirmed independently. In the current study, we explore the potential association between polymorphism of IL-4Rα and the susceptibility and severity of MG in a Chinese cohort.

Subjects and methods

Study population

Four hundred and eighty MG patients and 487 healthy controls (HC) were recruited from Affiliated Hospital of Qingdao University and Beijing Friendship Hospital, Capital Medical University. All MG patients and healthy controls were Han Chinese population in origin and non-consanguineous. All patients met the diagnostic criteria of MG: the typical symptoms of fluctuating muscle weakness and positive result of neostigmine test as the essential conditions, and at least one of the supporting conditions: anti-AChR antibody (AChR-Ab) positive and/or the amplitude decrement of low frequency repetitive nerve stimulation (RNS) >10%. AChR-Ab was detected with ELISA kits (RSR Limited, Cardiff, UK). Patients were followed at least twice a year, with additional follow-ups on the worsening and within 2–3 months thereafter. Informed consent was obtained from the MG patients or their guardian and the controls. The study protocol was approved by ethical committees in the two hospitals.

MG patients were classified into subgroups according to essential clinical variables in MG: gender, onset age (<15/15–50/ > 50 years) (18), thymoma (CT and/or pathology identified), AChR-Ab (+/−), muscle involvement at onset (ocular/generalized), and severity at the maximal worsening (modified Oosterhuis score 0–2/3–5) (19) during the first 2 years after onset. Previous study suggested that 82% MG patients reached the maximum worsening within 2 years after onset (20), thus, we analyzed the maximum Oosterhuis scores only in patients with clinical duration of 2 years or more. Due to potential interaction among the essential clinical variables (21, 22), a new classification for subgroup (23) was also applied in analysis (Figure 1).

Figure 1

SNPs selection and genotyping

Based on previous reports till July, 2013 and Pubmed SNP database, six SNPs (rs2107356, rs1805010, rs1805011, rs1805015, rs1801275, and rs8832) were selected in IL-4Rα gene. All SNPs met the requirement that the minor allele frequencies (MAF) > 5% (Table 1). Peripheral blood was collected in ethylenediamine tetraacetic acid (EDTA) vials, genomic DNA was isolated from peripheral blood using the Wizard genomic DNA purification kit (Biochain, Newark, CA, USA). SNP genotyping was performed with a custom-designed 2 × 48-Plex SNPscan™ Kit (Cat#: G0104, Genesky Biotechnologies Inc., Shanghai, China), which was based on double ligation and multiplex fluorescence PCR. Forty randomly selected samples from the patients and controls were typed for double-blinded quality control and were consistent with the original genotyping results.

Table 1

GeneSNPsAllelesPosition(bp)FunctionPHWMAFa
IL-4Rrs2107356C > T27312083promoter0.3020.380
rs1805010G > A27344882missense0.3530.484
rs1805011A > C27362551missense0.5360.083
rs1805015T > C27362859missense0.180.084
rs1801275A > G27363079missense0.2740.166
rs8832A > G273644663′ UTR0.1940.486

General characteristic of SNPs in IL-4R gene.

a

MAF value of East Asian [EAS] in the 1000 Genomes Project (http://www.internationalgenome.org/).

Statistical analysis

The online SNPstats software (http://bioinfo.iconcologia.net/snpstats/start.htm) was used to test the Hardy-Weinberg equilibrium (HWE) in HC. The χ2-test and Fisher exact test (SPSS 19.0) were used to compare allele frequencies and the SNPstats was used to compare genotype frequencies between MG and HC, between each MG subgroup and HC, between each pairs of MG subgroups (e.g., with thymoma vs. without thymoma), and among subgroups of the new classification and HC group. Bonferonni correction of the p-value was used for adjusting multiple tests. When there were significant differences in frequencies of alleles between MG subgroups and HC or among MG subgroups, logistic regression (SPSS 19.0) was used to adjust the potential confounders under codominant or additive models of inheritance. The criterion of significant difference was P < 0.05.

Result

General characteristics

Four hundred and eighty MG patients (291 females and 189 males) and 487 HC (238 females and 249 males) were enrolled in our study. The successful genotyping rates of six SNPs were 97.93 ~ 99.48%. The genotype distributions of six SNPs in HC were consistent with Hardy–Weinberg genetic equilibrium. The onset age of MG ranged 1 ~ 86 (median 40, interquartile range [IQR] 32). The age of HC group ranged 14 ~ 78 (median 45, IQR 24). 107 patients were with thymoma, 367 patients were without thymoma and 6 patients were not defined; 338 patients were AChR-Ab positive, 124 patients were AChR-Ab negative and 18 patients were not defined; 342 patients were ocular presenting at onset, 135 were generalized presenting and 3 patients were not defined. The disease duration of MG patients ranged 8–220 months (median 43, IQR 61). The mean number of follow-up each year was 6.5. The data of Oosterhuis score on the maximum worsening was available in 370 patients (76.9%). 216 patients were in mild subgroup (Oosterhuis score 0–2) and 154 patients were in severe subgroup (Oosterhuis score 3–5).

Comparison between the whole MG group and healthy control group

Frequency of rs1801275 G allele was higher in MG than that in HC (P = 0.036, OR = 1.293, 95%CI 1.017–1.643, PBon = 0.216). There were significant differences in genotype frequencies of rs1801275 between MG and HC under codominant (P = 0.048) and additive (P = 0.015) models of inheritance (Table 2).

Table 2

SNPsMAFGenotypesGenotype frequenciesP-value
MGHCMGHCCodominantAdditive
rs21073560.380.37CC174 (37.6)200 (41.3)0.460.43
CT222 (48.0)214 (44.2)
TT67 (14.5)70 (14.5)
rs18050100.50.48GG127 (26.7)135 (27.7)0.710.42
GA222 (46.7)233 (47.8)
AA126 (26.5)119 (24.4)
rs18050110.070.08AA409 (86.3)410 (84.2)0.510.52
CA62 (13.1)75 (15.4)
CC3 (0.6)2 (0.4)
rs18050150.070.08TT409 (86.3)409 (84.0)0.260.49
CT62 (13.1)77 (15.8)
CC3 (0.6)1 (0.2)
rs18012750.190.15aAA305 (64.8)348 (71.5)0.0480.015
GA156 (33.1)131 (26.9)
GG10 (2.1)8 (1.6)
rs88320.450.42AA146 (30.7)157 (32.4)0.240.18
GA233 (49.0)250 (51.5)
GG96 (20.2)78 (16.1)

Distribution of SNPs among all 480 MG patients without any subgroup distinction and 487 controls.

a

MG vs. HC, P = 0.036, OR = 1.293 95%CI (1.017–1.643), PBon = 0.216.

Comparison between MG subgroups and healthy control group and among MG subgroups

Frequencies of rs2107356 T allele and rs1805010 A allele were significantly higher in thymoma (+) subgroup than those in HC and in thymoma (−) subgroup both after Bonferonni correction. There were significant differences in genotype frequencies of rs2107356 and rs1805010 between thymoma (+) subgroup and HC and between thymoma (+) and thymoma (−) subgroups under codominant and additives models. There were no significant differences in frequencies of rs2107356 T allele and rs1805010 A allele between thymoma (−) subgroup and HC (Table 3).

Table 3

ControlMaleFemaleAnti-AchR AbThymoma
MGcontrolMGcontrolPositiveNegativepositivenegative
Number487189249291238338124107367
Rs2107356484*180247283237325120102355
CC200(41.3)66(36.7)105(42.5)108(38.2)95(40.1)122(37.5)44(36.7)31(30.4)141(39.7)
CT214(44.2)88(48.9)109(44.1)134(47.4)105(44.3)154(47.4)59(49.2)45(44.1)174(49.0)
TT70(14.5)26(14.4)33(13.4)41(14.5)37(15.6)49(15.1)17(14.2)26(25.5)40(11.3)
T allele0.370.390.350.380.380.390.390.48a0.36b
Rs1805010487186249289238334124106363
GG135(27.7)43(23.1)64(25.7)84(29.1)71(29.8)88(26.4)30(24.2)20(18.9)105(28.9)
GA233(47.8)91(48.9)125(50.2)131(45.3)108(45.4)161(48.2)53(42.7)45(42.5)175(48.2)
AA119(24.4)52(28.0)60(24.1)74(25.6)59(24.8)85(25.4)41(33.1)41(38.7)83(22.9)
A allele0.480.520.490.480.470.50.540.6c0.47d
Rs1801275487183249288238331122104361
AA348(71.5)108(59.0)175(70.3)197(68.4)173(72.7)202(61.0)87(71.3)66(63.5)236(65.4)
AG131(26.9)68(37.2)69(27.7)88(30.6)62(26.1)120(36.3)34(27.9)36(34.6)117(32.4)
GG8(1.6)7(3.8)5(2.0)3(1.0)3(1.3)9(2.7)1(0.8)2(1.9)8(2.2)
G allele0.150.220.160.160.140.21e0.150.190.18
ControlOnset age (year)Onset involvementOosterhuis score
<1515–50>50OcularGeneralized0–23–5
Number48771253156342135216154
Rs210735648468247148330130209149
CC200(41.3)29(42.6)93(37.7)52(35.1)126(38.2)47(36.1)81(38.8)62(40.3)
CT214(44.2)32(47.1)113(45.7)77(52.0)160(48.5)60(46.1)103(49.3)65(43.6)
TT70(14.5)7(10.3)41(16.6)19(12.8)44(13.3)23(17.7)25(12.0)24(16.1)
T allele0.370.340.390.390.380.410.370.38
Rs180501048770252153337135216152
GG135(27.7)21(30.0)69(27.4)37(24.2)88(26.1)39(28.9)56(25.9)45(29.6)
GA233(47.8)27(38.6)109(43.2)86(56.2)161(48.4)57(42.2)100(46.3)74(48.7)
AA119(24.4)22(31.4)74(29.4)30(19.6)86(25.5)39(28.9)60(27.8)33(21.7)
A allele0.480.510.510.480.50.50.510.46
Rs180127548769252150337131214152
AA348(71.5)51(73.9)163(64.7)91(60.7)216(64.1)88(67.2)141(65.9)98(65.1)
AG131(26.9)17(24.6)83(32.9)56(37.3)114(33.8)40(30.5)70(32.7)49(32.2)
GG8(1.6)1(1.4)6(2.4)3(2.0)7(2.1)3(2.3)3(1.4)5(3.3)
G allele0.150.140.190.210.190.180.180.19

Frequencies of alleles and genotypes of rs1805010, rs1801275 and rs2107356 in MG subgroups and controls [case (%)].

*

Italic bold refers to the number of successfully genotyped patients with adequate clinical data in each group.

a

. Thymoma (+) MG vs. HC, P = 0.003, OR = 1.572 (1.160–2.132), PBon = 0.018

codominant model, thymoma (+) MG vs. HC, P = 0.017

additive model, thymoma (+) MG vs. HC, P = 0.0055

b

. Thymoma (+) MG vs. Thymoma (−) MG, P = 0.002, OR = 1.627 (1.188–2.229), PBon = 0.012

codominant model, thymoma (+) MG vs. thymoma (−) MG, P = 0.0025

additive model, thymoma (+) MG vs. thymoma (−) MG, P = 0.0023

c

. Thymoma (+) MG vs. HC, P = 0.002, OR = 1.596 (1.180-2.158), PBon = 0.012

codominant model, thymoma (+) MG vs. HC, P = 0.008

additive model, thymoma (+) MG vs. HC, P = 0.003

d

. Thymoma (+) MG vs. Thymoma (−) MG, P = 0.001, OR = 1.687 (1.236–2.302), PBon = 0.006

codominant model, thymoma (+) MG vs. thymoma (−) MG, P = 0.004

additive model, thymoma (+) MG vs. thymoma (−) MG, P = 0.0012

e

. AChR-Ab (+) MG vs. HC, P = 0.003, OR = 1.482 (1.146–1.915), PBon = 0.018

codominant model, AChR Ab (+) MG vs. HC, P = 0.0032

additive model, AChR Ab (+) MG vs. HC, P = 0.0007.

Frequency of rs1801275 G allele was significantly higher in AChR-Ab (+) subgroup than that in HC after Bonferonni correction. There were significant differences in genotype frequencies of rs1801275 between AChR-Ab (+) subgroup and HC under codominant and additive models of inheritance. Frequency of rs1801275 G allele was significantly higher in male MG patients than that in male controls (P = 0.017, OR = 1.531, 95%CI 1.086–2.160), and higher in AChR-Ab (+) subgroup than that in AChR-Ab (−) subgroup (P = 0.039, OR = 1.522, 95%CI 1.020–2.271). After Bonferonni correction, the differences were not significant. There were no significant differences in allele and genotype frequencies between AChR-Ab (−) subgroup and HC (Table 3).

Adjustment of potential confounders in clinical variable based subgroup analysis

Logistic regression analysis was performed with thymoma (±) as dependent variables, and with onset age (<15/15–50/ > 50 years, <15 as reference), gender (male/female), AChR-Ab (±), muscle involvement at onset (ocular/generalized) and genotypes of rs2107356 or rs1805010 separately (under codominant or additives model) as independent variables. Onset age, AChR-Ab (+) and genotypes of rs2107356 and rs1805010 were found as independent risk factors for the presence of thymoma (Tables 4).

Table 4

VariablesRegression coefficientStandard errorWald x2-valueP-valueOR95%CI
Onset age (< 15 as reference)10.4480.005
15–50 year1.9640.62110.0170.0027.1282.112–24.056
>50 year1.6820.6366.9910.0085.3751.545–18.699
AChR-Ab positive1.1460.33111.9940.0013.1451.644–6.015
rs18050100.5350.16810.1980.0011.7081.230–2.372
Onset age (< 15 as reference)10.5880.005
15–50 year1.8460.6208.8630.0036.3351.879–21.358
>50 year1.3920.6384.7640.0294.0281.153–14.077
AChR-Ab positive1.0240.3319.5960.0022.7851.457–5.323
rs21073560.4400.1726.4980.0111.5521.107–2.176

Logistic regression analysis in subgroups classified by thymoma.a

a

The results were same under both codominant model and additive model.

Logistic regression analysis was performed with AChR-Ab (±) as dependent variables, and with onset age, gender, thymoma, muscle involvement at onset and genotypes (under codominant or additive model) as independent variables. Thymoma, ocular presenting and rs1801275 genotypes were found as independent risk factors for AChR-Ab positivity (Table 5).

Table 5

VariablesRegression coefficientStandard errorWald x2-valueP-valueOR95%CI
Onset age (< 15 as reference)1.3290.515
15–500.3580.3141.3040.2541.4310.774–2.645
>500.3030.3340.8270.3631.3540.704–2.604
Thymoma1.0540.32910.2830.0012.8711.507–5.469
Ocular presenting0.6930.2865.8810.0150.5000.286–0.876
rs18012750.5150.2255.2460.0221.6741.077–2.600

Logistic regression analysis in subgroups classified by AChR-Ab.a

a

The results were same under both codominant model and additive model.

Comparison between MG subgroups and healthy control group and among MG subgroups in the new classification

The new classification (Figure 1) used the combination of essential clinical features of MG to sort biological and clinical meaningful subgroups. It ensured any MG patients with sufficient data to be allocated into one subgroup and only in one subgroup. According to the new classification, 71 patients were juvenile onset MG (JMG). In 409 adult onset MG patients, 104 were with thymoma. In 300 adult MG patients without thymoma, 84 were AChR-Ab negative. In 204 non-thymoma AChR-Ab positive patients, 126 patients were early onset MG (EOMG, onset age 15–50) and 78 were late onset MG (LOMG, onset age > 50); 135 patients were ocular presenting and 67 were generalized presenting. The other patients had not relevant information to be defined into subgroups in this new classification.

Frequencies of rs2107356 T allele and rs1805010 A allele in adult thymoma subgroup was significantly higher than those in HC and in adult non-thymoma subgroup after Bonferonni correction. There were significant differences in genotype frequencies of rs2107356 and rs1805010 between adult thymoma subgroup and HC and between adult thymoma subgroup and adult non-thymoma subgroup under codominant and additive models (p < 0.05, Table 6). There was no difference in allele and genotype frequencies of the two SNPs between adult non-thymoma subgroup and HC.

Table 6

HCJuvenile (< 15)Adult (≥15)
ThymomaNon-thymoma
AchRAb (−)AchRAb (+)
15–50>50OcularGeneralized
Number48771409104300842041267813567
rs2107356484*6839599291811981237513363
CC200 (41.3)29 (42.6)145 (36.7)30 (30.3)114 (39.2)29 (35.8)79 (39.9)52 (42.3)27 (36.0)57 (42.9)19 (30.2)
CT214 (44.2)32 (47.1)190 (48.1)43 (43.4)144 (49.5)40 (49.5)98 (49.5)58 (47.2)40 (53.3)58 (43.6)37 (58.7)
TT70 (14.5)7 (10.3)60 (15.2)26 (26.3)33 (11.3)12 (14.8)21(10.6)13 (10.6)8 (10.7)18 (13.5)7 (11.1)
T allele0.370.340.390.48a0.36b0.400.350.340.370.350.40
rs180501048770405103297842011267513466
GG135 (27.7)21 (30.0)106 (26.2)20 (19.4)85 (28.6)18 (21.4)60 (29.9)42 (33.3)18 (24.0)38 (28.4)13 (19.7)
GA233 (47.8)27 (38.6)195 (48.1)43 (41.7)150 (50.5)42 (50.0)103 (51.2)58 (46.0)45 (60.0)65 (48.5)35 (53.0)
AA119 (24.4)22 (31.4)104 (25.7)40 (38.8)62 (20.9)24 (28.6)38 (18.9)26 (20.6)12 (16.0)1 (23.1)18 (27.3)
A allele0.480.510.500.60c0.46d0.540.450.440.460.470.43
rs180127548769402101296832011267513366
AA348 (71.5)51 (73.9)254 (63.2)63 (62.4)189 (63.9)57 (68.7)122 (60.7)77 (61.1)45 (60.0)86 (64.7)39 (59.1)
GA131 (26.9)17 (24.6)139 (34.6)36(35.6)100(33.8)25(30.1)73 (36.3)45 (35.7)28 (37.3)47 (35.3)24 (36.4)
GG8 (1.6)1 (1.4)9 (2.2)2(2.0)7(2.4)1 (1.2)6 (3.0)4 (3.2)2 (2.7)0 (0)3 (4.5)
G allele0.150.140.200.200.190.170.21e0.210.210.180.23

Frequencies of alleles and genotypes of rs1805010, rs1801275, and rs2107356 in MG subgroups of a new classification and the control group [case (%)].

*

Italic bold refers to the number of successfully genotyped patients with adequate clinical data in each group.

a

. adult thymoma MG vs. HC, P = 0.003, OR = 1.600 (1.176-2.177), PBon = 0.018

codominant model, adult thymoma MG vs. HC, P = 0.0056

additive model, adult thymoma MG vs. HC, P = 0.0017

b

. adult thymoma MG vs. adult non-thymoma MG, P = 0.003, OR = 1.634 (1.179–2.264), PBon = 0.018

codominant model, adult thymoma MG vs. non-thymoma MG, P = 0.0025

additive model, adult thymoma MG vs. non-thymoma MG, P = 0.0032

c

. adult thymoma MG vs. HC, P = 0.003, OR = 1.583 (1.166–2.148), PBon = 0.018

codominant model, adult thymoma MG vs. HC, P = 0.0034

additive model, adult thymoma MG vs. HC, P = 0.0013

d

. adult thymoma MG vs. adult non-thymoma MG, P = 0.001, OR = 1.731 (1.255–2.388), PBon = 0.006

codominant model, adult thymoma MG vs. non-thymoma MG, P = 0.0025

additive model, adult thymoma MG vs. adult non-thymoma MG, P = 0.0012

e

. adult non-thymoma AChR-Ab positive MG vs. HC, P = 0.006, OR = 1.509 (1.121-2.030), PBon = 0.036

codominant model, adult non-thymoma AChR-Ab positive MG vs. HC, P = 0.011

additive model, adult non-thymoma AChR-Ab positive MG vs. HC, P = 0.0026.

Frequency of rs1801275 G allele was significantly higher in adult non-thymoma AChR-Ab positive subgroup than that in HC after Bonferonni correction. There were significant differences in genotype frequencies of rs1801275 between adult non-thymoma AChR-Ab positive subgroup and HC group under codominant and additive models (Table 6). Frequency of rs1801275 G allele was significantly higher in adult non-thymoma AChR-Ab positive MG than in adult non-thymoma AChR-Ab negative MG (P = 0.023). After Bonferonni correction, the difference was not significant. There were no significantly difference in rs1801275 G allele and genotype frequencies between adult non-thymoma AChR-Ab negative subgroup and HC group. Frequency of rs1801275 G allele was higher in EOMG, LOMG, ocular presenting and generalized presenting subgroups than those in HC, but there were no allelic differences between EOMG subgroup and LOMG subgroup, and between ocular presenting subgroup and generalized presenting subgroup (Table 6).

Discussion

The current study tried to explore the association between candidate SNPs and the susceptibility and severity of MG in a representing MG cohort. Significant association was found between rs2107356 and rs1805010 and thymoma subgroup, and between rs1801275 and AChR-Ab positivie subgroup. No association was found among these SNPs and severity of MG. In previous subgroup analysis (21, 22), we found any single clinical variable of MG may be confound with other clinical variables, therefore, we adjusted these variables to analysis whether genotypes are independent risk factors for the association with a specific clinical feature. We found the association of rs2107356 and rs1805010 with thymoma and association of rs1801275 with AChR-Ab positivity. We then used the new classification which was based on the current insight of pathogenic mechanism and subgroup reasoning of MG to minimize the interaction among different clinical variables. We found association of rs2107356 and rs1805010 with adult thymoma associated MG; and found association of rs1801275 with adult non-thymoma AChR-Ab positivie MG.

Although allele and genotype frequencies of rs2107356 and rs1805010 were significantly higher in the adult thymoma MG subgroup than those in HC and in adult non-thymoma MG subgroup, there was no difference between adult non-thymoma MG subgroup and HC. Thymoma is a tumor of thymic epithelium. The association of thymoma and MG is well-known (24). As far as we known, this is the first time to find the association of these two SNPs with thymoma.

Rs2107356 is a promoter SNP which might increase the expression of IL-4R, which should be studied in the future. Rs2107356 was also found associated with gastric cancer (25) and multiple myeloma (26). Rs1805010 is a gain-of-function mutation in the extracellular domain of IL-4Rα subunit. This variant enhances IL-4 signaling by keeping STAT6 transcription factor abnormally active, which was thought not due to modulation of intracellular mechanisms (27). Rs1805010 was also found associated with renal cell cancer (28) and cervical cancer (29). Most of the association was reported in epithelic tumors. IL-4 supports cellular proliferation of human tumor cell lines in a concentration-dependent manner. IL-4R blockade dramatically reduces lymphatic and hematogenous metastases of some kinds of tumors in vivo. IL-4 also acts as a survival factor on tumor cells by increasing the expression of anti-apoptotic molecules. These effects of IL-4 on tumor cells predominantly signals through the type II IL-4R complex (5). The role of IL-4R on the pathogenesis in thymoma might resemble with those in other tumors.

In our study, frequency of G allele in rs1801275 was higher in MG patients than that in HC group, especially in adult non-thymoma AChR-Ab positive MG subgroup. Genotypes of rs1801275 were found to be independent factors for the presence of AChR-Ab in MG. The rs1801275 is another gain-of-function mutation in the intracellular domain of IL-4Rα, which enhancs IL-4Rα signal transduction via STAT6 pathway, which has important effects on immune modulating (30, 31). This variant was also found associated with IL-4-directed and IL-6-dependent subversion of Treg cells into Th17-like cells (16). Th17 cell is critical in the pathogenesis of MG (32, 33), The effect of rs1801275 on the balance between Treg and Th17 may explain our results on the association of this mutant with AChR-Ab positive MG.

Our study had limitations. We did not enroll thymoma patients without MG, therefore, the association between IL-4Rα gene polymorphism and thymoma should be further confirmed. Moreover, previous studies have shown that the T cell receptors (TCRs) of MG with thymoma are different from those of MG without thymoma (34). Whether IL-4Rα plays a role in the growth of thymoma and induces the T cells expressing specific TCRs to metastasize to periphery to participate in the production of AChR-Ab wait for further researches.

In conclusion, in Chinese Han population, rs2107356 and rs1805010 were found to be associated with adult thymoma associated MG, and rs1801275 was found to be associated with adult non-thymoma AChR-Ab positive MG.

Statements

Ethics statement

This study was carried out in accordance with the recommendations of ethical committees of Affiliated Hospital of Qingdao University, and Beijing Friendship Hospital, Capital Medical University. The protocol was approved by the ethical committees of Affiliated Hospital of Qingdao University, and Beijing Friendship Hospital, Capital Medical University. All subjects gave written informed consent in accordance with the Declaration of Helsinki.

Author contributions

PJ, Y-XY, and H-FL conceptualized and designed the study, interpreted the data, and wrote the manuscript. Y-XY and YH designed the genotyping experiment. H-JH performed the AChR antibody testing. PJ, Y-XY, and YH performed statistical analysis. YQ, X-JD, and MS contributed to discussion. H-FL, YX, XG, XZ, and C-KG diagnosed and treated patients in this study. YX, XG, and Y-XY maintained the database.

Funding

The present study was supported by the National Natural Science Foundation of China (No. 81070963 and 81771362 to H-FL, No. 81400068 to YX), Natural Science Foundation of Shandong Province (No. ZR2010HM019 to H-FL), Key Research and Development Plan of Shandong Province (No. 2016GSF201210 to H-FL), Beijing Natural Science Foundation (No. 7142042 to YX), Qingdao Technology Program for Health and Welfare (No. 17-3-3-26-nsh to H-FL) and Research Grant from Qilu Hospital of Shandong University (Qingdao) (No. QDKY20152D01 to H-FL).

Acknowledgments

Participants in this study were enrolled and followed from Affiliated Hospital of Qingdao University and Beijing Friendship Hospital, Capital Medical University. The authors thank Genesky Biotechnologies Inc., Shanghai, China, for technical support.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1.

    MeriggioliMNSandersDBReményiVRózsaCFalusA. Autoimmune myasthenia gravis: emerging clinical and biological heterogeneity. Lancet Neurol. (2009) 8:47590. 10.1016/S1474-4422(09)70063-8

  • 2.

    YangWCHwangYSChenYYLiuCLShenCNHongWHet al. Interleukin-4 supports the suppressive immune responses elicited by regulatory T cells. Front Immunol. (2017) 14:1508. 10.3389/fimmu.2017.01508

  • 3.

    UzawaAKawaguchiNHimuroKKanaiTKuwabaraS. Serum cytokine and chemokine profiles in patients with myasthenia gravis. Clin Exp Immunol. (2014) 176:23237. 10.1111/cei.12272

  • 4.

    OstlieNMilaniMWangWOkitaDConti-FineBM. Absence of IL-4 facilitates the development of chronic autoimmune myasthenia gravis in C57BL/6 mice. J Immunol. (2003) 170:60412. 10.4049/jimmunol.170.1.604

  • 5.

    SuzukiALelandPJoshiBHPuriRK. Targeting of IL-4 and IL-13 receptors for cancer therapy. Cytokine (2015) 75:7988. 10.1016/j.cyto.2015.05.026

  • 6.

    ShandleySMartinezSKrolickK. IL-4 receptor as a bridge between the immune system and muscle in experimental myasthenia gravis I: up-regulation of muscle IL-15 by IL-4. Clin Immunol. (2009) 132:24656. 10.1016/j.clim.2009.03.523

  • 7.

    HuangDXiaSZhouYPirskanenRLiuLLefvertAK. No evidence for interleukin-4 gene conferring susceptibility to myasthenia gravis. J Neuroimmunol. (1998) 92:20811. 10.1016/S0165-5728(98)00186-6

  • 8.

    WallisSKCooneyLAEndresJLLeeMJRyuJSomersECet al. A polymorphism in the interleukin-4 receptor affects the ability of interleukin-4 to regulate Th17 cells: a possible immunoregulatory mechanism for genetic control of the severity of rheumatoid arthritis. Arthritis Res Ther. (2011) 13:R15. 10.1186/ar3239

  • 9.

    KrabbenAWilsonAGde RooyDPZhernakovaABrouwerELindqvistEet al. Association of genetic variants in the IL4 and IL4R genes with the severity of joint damage in rheumatoid arthritis: a study in seven cohorts. Arthritis Rheum. (2013) 65:305157. 10.1002/art.38141

  • 10.

    LeipeJSchrammMAProtsISchulze-KoopsHSkapenkoA. Increased Th17 cell frequency and poor clinical outcome in rheumatoid arthritis are associated with a genetic variant in the IL4R gene, Rs1805010. Arthritis Rheumatol. (2014) 66:116575. 10.1002/art.38343

  • 11.

    XuYChenZQLiYMGongJQLiASChenMet al. Correlation between some Th1 and Th2 cytokine receptor gene polymorphisms and systemic lupus erythematosus in Chinese patients. Int J Dermatol. (2007) 46:112935. 10.1111/j.1365-4632.2007.03258.x

  • 12.

    YabikuKHayashiMKomiyaIYamadaTKinjoYOhshiroYet al. Polymorphisms of interleukin (IL)-4 receptor alpha and signal transducer and activator of transcription-6 (Stat6) are associated with increased IL-4Rα-Stat6 signalling in lymphocytes and elevated serum IgE in patients with Graves' disease. Clin Exp Immunol. (2007) 148:42531. 10.1111/j.1365-2249.2007.03366.x

  • 13.

    MiyakeYTanakaKArakawaM. Case-control study of eczema in relation to IL4Rα genetic polymorphisms in Japanese women: The Kyushu Okinawa Maternal and Child Health Study. Scand J Immunol. (2013) 77:41318. 10.1111/sji.12043

  • 14.

    SlagerREOtulanaBAHawkinsGAYenYPPetersSPWenzelSEet al. IL-4 receptor polymorphisms predict reduction in asthma exacerbations during response to an anti-IL-4 receptor α antagonist. J Allergy Clin Immunol. (2012) 130:516-22.e4. 10.1016/j.jaci.2012.03.030

  • 15.

    HuangZYChengBJCaiGJZhangBF. Association of the IL-4R Q576R polymorphism and asthma in the Chinese Han population: a meta-analysis. Genet Mol Res. (2015) 14:290011. 10.4238/2015.March.31.21

  • 16.

    MassoudAHCharbonnierLMLopezDPellegriniMPhipatanakulWChatilaTA. An asthma-associated IL4R variant exacerbates airway inflammation by promoting conversion of regulatory T cells to TH17-like cells. Nat Med. (2016) 22:101322. 10.1038/nm.4147

  • 17.

    PálZVargaZSemseiÁReményiVRózsaCFalusAet al. Interleukin-4 receptor alpha polymorphisms in autoimmune myasthenia gravis in a Caucasian population. Hum Immunol. (2012) 73:19395. 10.1016/j.humimm.2011.11.001

  • 18.

    AarliJA. Myasthenia gravis in the elderly: is it different. Ann N Y Acad Sci. (2008) 1132:23843. 10.1196/annals.1405.040

  • 19.

    SkeieGOPandeyJPAarliJAGilhusNE. TNFA and TNFB polymorphisms in myasthenia gravis. Arch Neurol. (1999) 56:45761. 10.1001/archneur.56.4.457

  • 20.

    GrobDBrunnerNNambaTPagalaM. Lifetime course of myasthenia gravis. Muscle Nerve (2008) 37:14149. 10.1002/mus.20950

  • 21.

    YueYXHongYXieYHaoHJSuiYGuCKet al. Association study between IL-17A and IL-17F gene polymorphism and myasthenia gravis in Chinese patients. Neurol Sci. (2016) 37:12330. 10.1007/s10072-015-2375-y

  • 22.

    ZhangXDingXJWangQYueYXXieYHaoHJet al. Rs3761389 polymorphism in autoimmune regulator (AIRE) gene is associated with susceptibility of myasthenia gravis in Chinese patients. J Clin Neurosci. (2017) 40:18084. 10.1016/j.jocn.2017.02.049

  • 23.

    LiHFXieYYueYX. Myasthenia gravis: subgroup classifications. Lancet Neurol. (2016) 15:3556. 10.1016/S1474-4422(16)00032-6

  • 24.

    RomiF. Thymoma in myasthenia gravis: from diagnosis to treatment. Autoimmune Dis. (2011) 2011:474512. 10.4061/2011/474512

  • 25.

    XiaHZDuWDWuQChenGZhouYTangXFet al. E-selectin rs5361 and FCGR2A rs1801274 variants were associated with increased risk of gastric cancer in a Chinese population. Mol Carcinog. (2012) 51:597607. 10.1002/mc.20828

  • 26.

    BrownEELanQZhengTZhangYWangSSHoar-ZahmSet al. Common variants in genes that mediate immunity and risk of multiple myeloma. Int J Cancer (2007) 120:271522. 10.1002/ijc.22618

  • 27.

    FordAQHellerNMStephensonLBoothbyMRKeeganAD. An atopy-associated polymorphism in the ectodomain of the IL-4Rα chain (V50) regulates the persistence of STAT6 phosphorylation. J Immunol. (2009) 183:160716. 10.4049/jimmunol.0803266

  • 28.

    LuoYYeZLiKChenRLiSPangJ. Associations between polymorphisms in the IL-4 and IL-4 receptor genes and urinary carcinomas: a meta-analysis. Int J Clin Exp Med. (2015) 8:122733.

  • 29.

    IvanssonELGustavssonIMMagnussonJJSteinerLLMagnussonPKErlichHAet al. Variants of chemokine receptor 2 and interleukin 4 receptor, but not interleukin 10 or Fas ligand, increase risk of cervical cancer. Int J Cancer (2007) 121:24517. 10.1002/ijc.22989

  • 30.

    HersheyGKFriedrichMFEssweinLAThomasMLChatilaTA. The association of atopy with a gain-of-function mutation in the alpha subunit of the interleukin-4 receptor. N Engl J Med. (1997) 337:17205. 10.1056/NEJM199712113372403

  • 31.

    RismaKAWangNAndrewsRPCunninghamCMEricksenMBBernsteinJAet al. V75R576 IL-4 receptor alpha is associated with allergic asthma and enhanced IL-4 receptor function. J Immunol. (2002) 169:160410. 10.4049/jimmunol.169.3.1604

  • 32.

    SchaffertHPelzASaxenaALosenMMeiselAThielAet al. IL-17-producing CD4(+) T cells contribute to the loss of B-cell tolerance in experimental autoimmune myasthenia gravis. Eur J Immunol. (2015) 45:133947. 10.1002/eji.201445064

  • 33.

    VillegasJAVan WassenhoveJLe PanseRBerrih-AkninSDraginN. An imbalance between regulatory T cells and T helper 17 cells in acetylcholine receptor-positive myasthenia gravis patients. Ann N Y Acad Sci. (2018) 1413:15462. 10.1111/nyas.13591

  • 34.

    InadaKOkumuraMShionoHInoueMKadotaYOhtaMet al. Role of positive selection of thymoma-associated T cells in the pathogenesis of myasthenia gravis. J Surg Res. (2005) 126:344010.1016/j.jss.2005.01.023

Summary

Keywords

myasthenia gravis, interleukin-4 receptor, polymorphism, susceptibility, severity

Citation

Jiang P, Yue Y-X, Hong Y, Xie Y, Gao X, Gu C-K, Hao H-J, Qin Y, Ding X-J, Song M, Li H-F and Zhang X (2018) IL-4Rα Polymorphism Is Associated With Myasthenia Gravis in Chinese Han Population. Front. Neurol. 9:529. doi: 10.3389/fneur.2018.00529

Received

16 February 2018

Accepted

13 June 2018

Published

10 July 2018

Volume

9 - 2018

Edited by

Genevieve Gourdon, Institut National de la Santé et de la Recherche Médicale (INSERM), France

Reviewed by

Rozen Le Panse, INSERM U974 Thérapie des Maladies du Muscle Strié (Institut de Myologie), France; Raghav Govindarajan, University of Missouri, United States

Updates

Copyright

*Correspondence: Hai-Feng Li Xu Zhang

This article was submitted to Neuromuscular Diseases, a section of the journal Frontiers in Neurology

†These authors have contributed equally to this work.

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics