[Skip to Content]
Access to paid content on this site is currently suspended due to excessive activity being detected from your IP address 54.159.129.152. Please contact the publisher to request reinstatement.
Sign In
Individual Sign In
Create an Account
Institutional Sign In
OpenAthens Shibboleth
[Skip to Content Landing]
Download PDF
Figure 1.
Myasthenia Gravis Cases and P Values From the Genome-Wide Association Study of Myasthenia Gravis
Myasthenia Gravis Cases and P Values From the Genome-Wide Association Study of Myasthenia Gravis

A, Quartile-quartile plot showing the distribution of expected vs observed P values for the US discovery cohort (972 myasthenia gravis cases and 1977 control individuals; λ = 1.036). B, Manhattan plot of genome-wide association results for all myasthenia gravis cases (972 cases and 1977 control individuals). C, Density plot of myasthenia gravis cases (n = 972) showing bimodal frequency distribution with disease being more common in younger women and older men. D, Plot of genome-wide association results for late-onset myasthenia gravis cases (737 cases and 1977 control individuals). E, Plot of genome-wide association results for early-onset myasthenia gravis cases (235 cases and 1977 control individuals). CTLA4 indicates cytotoxic T-lymphocyte–associated protein 4; HLA, human leukocyte antigen; and TNFRSF11A, tumor necrosis factor receptor 4 superfamily, member 11a, NFKB activator gene. The orange horizontal line depicts the Bonferroni-adjusted significance threshold (5.0 × 10−8). Single-nucleotide polymorphisms exceeding this threshold are indicated with red dots, whereas single-nucleotide polymorphisms with a P < 10−4 and >5.0 × 10−8 are indicated with yellow dots.

Figure 2.
Regional Association Plots Across Loci Identified by the Genome-Wide Association Study of Myasthenia Gravis
Regional Association Plots Across Loci Identified by the Genome-Wide Association Study of Myasthenia Gravis

Regional association plots of the signals at the major histocompatibility complex (MHC) region (A) and CTLA4 (B) in all myasthenia gravis cases (972 cases and 1977 control individuals), as well as the MHC region in early-onset myasthenia gravis cases (C; 235 cases and 1977 control individuals), the MHC region in late-onset cases (D), and TNFRSF11A in late-onset myasthenia gravis cases (E; 737 cases and 1977 control individuals). Single-nucleotide polymorphisms are colored on the basis of their linkage disequilibrium with the labeled hit single-nucleotide polymorphism. Recombination rates estimated from the CEU HapMap population are represented by the blue line and genes are marked as arrows. Red diamonds represent the most associated single-nucleotide polymorphism in the discovery cohort and blue diamonds, P values for the combined discovery and replication cohorts. F, CTLA4 regulates T-cell activation by antigen-presenting cells (APCs). TCR indicates T-cell receptor.

Table 1.  
Loci Showing Genome-Wide Association With Myasthenia Gravis
Loci Showing Genome-Wide Association With Myasthenia Gravis
Table 2.  
Heritability Estimates Based on the Discovery Cohort of Myasthenia Gravis
Heritability Estimates Based on the Discovery Cohort of Myasthenia Gravis
1.
Drachman  DB.  Myasthenia gravis. N Engl J Med. 1994;330(25):1797-1810.
PubMedArticle
2.
Spillane  J, Higham  E, Kullmann  DM.  Myasthenia gravis. BMJ. 2012;345:e8497.
PubMedArticle
3.
Higuchi  O, Hamuro  J, Motomura  M, Yamanashi  Y.  Autoantibodies to low-density lipoprotein receptor-related protein 4 in myasthenia gravis. Ann Neurol. 2011;69(2):418-422.
PubMedArticle
4.
Zhang  B, Shen  C, Bealmear  B,  et al.  Autoantibodies to agrin in myasthenia gravis patients. PLoS One. 2014;9(3):e91816.
PubMedArticle
5.
Carr  AS, Cardwell  CR, McCarron  PO, McConville  J.  A systematic review of population based epidemiological studies in Myasthenia Gravis. BMC Neurol. 2010;10:46.
PubMedArticle
6.
Alshekhlee  A, Miles  JD, Katirji  B, Preston  DC, Kaminski  HJ.  Incidence and mortality rates of myasthenia gravis and myasthenic crisis in US hospitals. Neurology. 2009;72(18):1548-1554.
PubMedArticle
7.
Jani-Acsadi  A, Lisak  RP.  Myasthenic crisis: guidelines for prevention and treatment. J Neurol Sci. 2007;261(1-2):127-133.
PubMedArticle
8.
Cavalcante  P, Bernasconi  P, Mantegazza  R.  Autoimmune mechanisms in myasthenia gravis. Curr Opin Neurol. 2012;25(5):621-629.
PubMedArticle
9.
Vincent  A, Palace  J, Hilton-Jones  D.  Myasthenia gravis. Lancet. 2001;357(9274):2122-2128.
PubMedArticle
10.
Sanders  DB, Andrews  I, Howard  JF, Massey  JM.  Seronegative myasthenia gravis. Neurology. 1997;48(suppl 5):S40-S45. doi:10.1212/WNL.48.Suppl_5.40S.
PubMedArticle
11.
Somnier  FE.  Myasthenia gravis. Dan Med Bull. 1996;43(1):1-10.
PubMed
12.
LeMaoult  J, Szabo  P, Weksler  ME.  Effect of age on humoral immunity, selection of the B-cell repertoire and B-cell development. Immunol Rev. 1997;160:115-126.
PubMedArticle
13.
Carlsson  B, Wallin  J, Pirskanen  R, Matell  G, Smith  CI.  Different HLA DR-DQ associations in subgroups of idiopathic myasthenia gravis. Immunogenetics. 1990;31(5-6):285-290.
PubMedArticle
14.
Pirskanen  R.  Genetic aspects in myasthenia gravis: a family study of 264 Finnish patients. Acta Neurol Scand. 1977;56(5):365-388.
PubMedArticle
15.
Gregersen  PK, Kosoy  R, Lee  AT,  et al.  Risk for myasthenia gravis maps to a (151) Pro→Ala change in TNIP1 and to human leukocyte antigen-B*08. Ann Neurol. 2012;72(6):927-935.
PubMedArticle
16.
Wang  XB, Pirskanen  R, Giscombe  R, Lefvert  AK.  Two SNPs in the promoter region of the CTLA-4 gene affect binding of transcription factors and are associated with human myasthenia gravis. J Intern Med. 2008;263(1):61-69.
PubMedArticle
17.
Christensen  PB, Jensen  TS, Tsiropoulos  I,  et al.  Associated autoimmune diseases in myasthenia gravis: a population-based study. Acta Neurol Scand. 1995;91(3):192-195.
PubMedArticle
18.
Tsao  CY, Mendell  JR, Lo  WD, Luquette  M, Rennebohm  R.  Myasthenia gravis and associated autoimmune diseases in children. J Child Neurol. 2000;15(11):767-769.
PubMedArticle
19.
Szobor  A.  Myasthenia gravis: familial occurrence: a study of 1100 myasthenia gravis patients. Acta Med Hung. 1989;46(1):13-21.
PubMed
20.
Padua  L, Stalberg  E, LoMonaco  M, Evoli  A, Batocchi  A, Tonali  P.  SFEMG in ocular myasthenia gravis diagnosis. Clin Neurophysiol. 2000;111(7):1203-1207.
PubMedArticle
21.
Hoch  W, McConville  J, Helms  S, Newsom-Davis  J, Melms  A, Vincent  A.  Auto-antibodies to the receptor tyrosine kinase MuSK in patients with myasthenia gravis without acetylcholine receptor antibodies. Nat Med. 2001;7(3):365-368.
PubMedArticle
22.
Hamza  TH, Zabetian  CP, Tenesa  A,  et al.  Common genetic variation in the HLA region is associated with late-onset sporadic Parkinson’s disease. Nat Genet. 2010;42(9):781-785.
PubMedArticle
23.
R Core Team. A language and environment for statistical computing. 2012. http://www.R-project.org/. Accessed September 30, 2013.
24.
Purcell  S, Neale  B, Todd-Brown  K,  et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559-575.
PubMedArticle
25.
Li  Y, Willer  C, Sanna  S, Abecasis  G.  Genotype imputation. Annu Rev Genomics Hum Genet. 2009;10:387-406.
PubMedArticle
26.
Yang  J, Benyamin  B, McEvoy  BP,  et al.  Common SNPs explain a large proportion of the heritability for human height. Nat Genet. 2010;42(7):565-569.
PubMedArticle
27.
Pe’er  I, Yelensky  R, Altshuler  D, Daly  MJ.  Estimation of the multiple testing burden for genomewide association studies of nearly all common variants. Genet Epidemiol. 2008;32(4):381-385.
PubMedArticle
28.
Lee  SH, Wray  NR, Goddard  ME, Visscher  PM.  Estimating missing heritability for disease from genome-wide association studies. Am J Hum Genet. 2011;88(3):294-305.
PubMedArticle
29.
Yang  J, Manolio  TA, Pasquale  LR,  et al.  Genome partitioning of genetic variation for complex traits using common SNPs. Nat Genet. 2011;43(6):519-525.
PubMedArticle
30.
Djabiri  F, Caillat-Zucman  S, Gajdos  P,  et al.  Association of the AChRalpha-subunit gene (CHRNA), DQA1*0101, and the DR3 haplotype in myasthenia gravis. Evidence for a three-gene disease model in a subgroup of patients. J Autoimmun. 1997;10(4):407-413.
PubMedArticle
31.
Schneider  H, Downey  J, Smith  A,  et al.  Reversal of the TCR stop signal by CTLA-4. Science. 2006;313(5795):1972-1975.
PubMedArticle
32.
Waterhouse  P, Penninger  JM, Timms  E,  et al.  Lymphoproliferative disorders with early lethality in mice deficient in Ctla-4. Science. 1995;270(5238):985-988.
PubMedArticle
33.
Qureshi  OS, Zheng  Y, Nakamura  K,  et al.  Trans-endocytosis of CD80 and CD86: a molecular basis for the cell-extrinsic function of CTLA-4. Science. 2011;332(6029):600-603.
PubMedArticle
34.
Dubois  PC, Trynka  G, Franke  L,  et al.  Multiple common variants for celiac disease influencing immune gene expression. Nat Genet. 2010;42(4):295-302.
PubMedArticle
35.
Barrett  JC, Clayton  DG, Concannon  P,  et al; Type 1 Diabetes Genetics Consortium.  Genome-wide association study and meta-analysis find that over 40 loci affect risk of type 1 diabetes. Nat Genet. 2009;41(6):703-707.
PubMedArticle
36.
Chu  X, Pan  CM, Zhao  SX,  et al; China Consortium for Genetics of Autoimmune Thyroid Disease.  A genome-wide association study identifies two new risk loci for Graves’ disease. Nat Genet. 2011;43(9):897-901.
PubMedArticle
37.
Gregersen  PK, Amos  CI, Lee  AT,  et al.  REL, encoding a member of the NF-kappaB family of transcription factors, is a newly defined risk locus for rheumatoid arthritis. Nat Genet. 2009;41(7):820-823.
PubMedArticle
38.
Schiff  M, Keiserman  M, Codding  C,  et al.  Efficacy and safety of abatacept or infliximab vs placebo in ATTEST: a phase III, multi-centre, randomised, double-blind, placebo-controlled study in patients with rheumatoid arthritis and an inadequate response to methotrexate. Ann Rheum Dis. 2008;67(8):1096-1103.
PubMedArticle
39.
Vincenti  F, Charpentier  B, Vanrenterghem  Y,  et al.  A phase III study of belatacept-based immunosuppression regimens versus cyclosporine in renal transplant recipients (BENEFIT Study). Am J Transplant. 2010;10(3):535-546.
PubMedArticle
40.
McIntosh  KR, Linsley  PS, Bacha  PA, Drachman  DB.  Immunotherapy of experimental autoimmune myasthenia gravis: selective effects of CTLA4Ig and synergistic combination with an IL2-diphtheria toxin fusion protein. J Neuroimmunol. 1998;87(1-2):136-146.
PubMedArticle
41.
McIntosh  KR, Linsley  PS, Drachman  DB.  Immunosuppression and induction of anergy by CTLA4Ig in vitro: effects on cellular and antibody responses of lymphocytes from rats with experimental autoimmune myasthenia gravis. Cell Immunol. 1995;166(1):103-112.
PubMedArticle
42.
Anderson  DM, Maraskovsky  E, Billingsley  WL,  et al.  A homologue of the TNF receptor and its ligand enhance T-cell growth and dendritic-cell function. Nature. 1997;390(6656):175-179.
PubMedArticle
43.
Dougall  WC, Glaccum  M, Charrier  K,  et al.  RANK is essential for osteoclast and lymph node development. Genes Dev. 1999;13(18):2412-2424.
PubMedArticle
44.
Hughes  AE, Shearman  AM, Weber  JL,  et al.  Genetic linkage of familial expansile osteolysis to chromosome 18q. Hum Mol Genet. 1994;3(2):359-361.
PubMedArticle
45.
Guerrini  MM, Sobacchi  C, Cassani  B,  et al.  Human osteoclast-poor osteopetrosis with hypogammaglobulinemia due to TNFRSF11A (RANK) mutations. Am J Hum Genet. 2008;83(1):64-76.
PubMedArticle
46.
Freedman  ML, Reich  D, Penney  KL,  et al.  Assessing the impact of population stratification on genetic association studies. Nat Genet. 2004;36(4):388-393.
PubMedArticle
47.
Buhler  S, Sanchez-Mazas  A.  HLA DNA sequence variation among human populations: molecular signatures of demographic and selective events. PLoS One. 2011;6(2):e14643.
PubMedArticle
48.
Namba  T, Brunner  NG, Brown  SB, Muguruma  M, Grob  D.  Familial myasthenia gravis: report of 27 patients in 12 families and review of 164 patients in 73 families. Arch Neurol. 1971;25(1):49-60.
PubMedArticle
49.
Keller  MF, Ferrucci  L, Singleton  AB,  et al.  Genome-wide analysis of the heritability of amyotrophic lateral sclerosis. JAMA Neurol. 2014;71(9):1123-1134.
PubMedArticle
50.
Keller  MF, Saad  M, Bras  J,  et al; International Parkinson’s Disease Genomics Consortium (IPDGC); Wellcome Trust Case Control Consortium 2 (WTCCC2).  Using genome-wide complex trait analysis to quantify ‘missing heritability’ in Parkinson’s disease. Hum Mol Genet. 2012;21(22):4996-5009.
PubMedArticle
51.
Chan  KH, Lachance  DH, Harper  CM, Lennon  VA.  Frequency of seronegativity in adult-acquired generalized myasthenia gravis. Muscle Nerve. 2007;36(5):651-658.
PubMedArticle
52.
Aarli  JA.  Myasthenia gravis in the elderly: is it different? Ann N Y Acad Sci. 2008;1132:238-243.
PubMedArticle
53.
Evoli  A, Tonali  PA, Padua  L,  et al.  Clinical correlates with anti-MuSK antibodies in generalized seronegative myasthenia gravis. Brain. 2003;126(pt 10):2304-2311.
PubMedArticle
  • Download PDF
  • Cite This

    Citation

    Renton AE, Pliner HA, Provenzano C, Evoli A, Ricciardi R, Nalls MA, Marangi G, Abramzon Y, Arepalli S, Chong S, Hernandez DG, Johnson JO, Bartoccioni E, Scuderi F, Maestri M, Gibbs JR, Errichiello E, Chiò A, Restagno G, Sabatelli M, Macek M, Scholz SW, Corse A, Chaudhry V, Benatar M, Barohn RJ, McVey A, Pasnoor M, Dimachkie MM, Rowin J, Kissel J, Freimer M, Kaminski HJ, Sanders DB, Lipscomb B, Massey JM, Chopra M, Howard JF, Koopman WJ, Nicolle MW, Pascuzzi RM, Pestronk A, Wulf C, Florence J, Blackmore D, Soloway A, Siddiqi Z, Muppidi S, Wolfe G, Richman D, Mezei MM, Jiwa T, Oger J, Drachman DB, Traynor BJ. A Genome-Wide Association Study of Myasthenia Gravis. JAMA Neurol. 2015;72(4):396-404. doi:10.1001/jamaneurol.2014.4103

    Download citation file:

Original Investigation
April 2015

A Genome-Wide Association Study of Myasthenia Gravis

Author Affiliations
  • 1Neuromuscular Diseases Research Unit, Laboratory of Neurogenetics, National Institute on Aging, National Institutes of Health, Porter Neuroscience Research Center, Bethesda, Maryland
  • 2Institute of General Pathology, Catholic University, Rome, Italy
  • 3Institute of Neurology, Catholic University, Rome, Italy
  • 4Department of Neuroscience, Cisanello Hospital, University of Pisa, Pisa, Italy
  • 5Molecular Genetics Section, Laboratory of Neurogenetics, National Institute on Aging, National Institutes of Health, Porter Neuroscience Research Center, Bethesda, Maryland
  • 6Institute of Medical Genetics, Catholic University, Rome, Italy
  • 7Genomics Technology Group, Laboratory of Neurogenetics, National Institute on Aging, National Institutes of Health, Porter Neuroscience Research Center, Bethesda, Maryland
  • 8Computational Biology Core, Laboratory of Neurogenetics, National Institute on Aging, National Institutes of Health, Porter Neuroscience Research Center, Bethesda, Maryland
  • 9Rita Levi Montalcini Department of Neuroscience, University of Turin, Turin, Italy
  • 10Molecular Genetics Unit, Department of Clinical Pathology, ASO OIRM-S Anna, Turin, Italy
  • 11Department of Neurology, Johns Hopkins School of Medicine, Baltimore, Maryland
  • 12Department of Neurology, University of Miami, Miami, Florida
  • 13Department of Neurology, University of Kansas Medical Center, Kansas City
  • 14Department of Neurology, University of Illinois College of Medicine, Chicago
  • 15Department of Neurology, Ohio State University Medical Center, Columbus
  • 16Department of Neurology, George Washington University, Washington, DC
  • 17Department of Neurology, Duke University Medical Center, Durham, North Carolina
  • 18Department of Neurology, University of North Carolina, Chapel Hill
  • 19Department of Clinical Neurosciences, London Health Sciences Centre, London, Ontario, Canada
  • 20Department of Neurology, Indiana University–Purdue University, Indianapolis
  • 21Department of Neurology, Washington University School of Medicine, St Louis, Missouri
  • 22Department of Medicine, University of Alberta Hospital, Edmonton, Alberta, Canada
  • 23Department of Neurology, University at Buffalo SMBS, State University of New York, Buffalo
  • 24Department of Neurology, University of California, Davis Medical Center
  • 25Division of Neurology, University of British Columbia, Vancouver, British Columbia, Canada
JAMA Neurol. 2015;72(4):396-404. doi:10.1001/jamaneurol.2014.4103
Abstract

Importance  Myasthenia gravis is a chronic, autoimmune, neuromuscular disease characterized by fluctuating weakness of voluntary muscle groups. Although genetic factors are known to play a role in this neuroimmunological condition, the genetic etiology underlying myasthenia gravis is not well understood.

Objective  To identify genetic variants that alter susceptibility to myasthenia gravis, we performed a genome-wide association study.

Design, Setting, and Participants  DNA was obtained from 1032 white individuals from North America diagnosed as having acetylcholine receptor antibody–positive myasthenia gravis and 1998 race/ethnicity-matched control individuals from January 2010 to January 2011. These samples were genotyped on Illumina OmniExpress single-nucleotide polymorphism arrays. An independent cohort of 423 Italian cases and 467 Italian control individuals were used for replication.

Main Outcomes and Measures  We calculated P values for association between 8 114 394 genotyped and imputed variants across the genome and risk for developing myasthenia gravis using logistic regression modeling. A threshold P value of 5.0 × 10−8 was set for genome-wide significance after Bonferroni correction for multiple testing.

Results  In the overall case-control cohort, we identified association signals at CTLA4 (rs231770; P = 3.98 × 10−8; odds ratio, 1.37; 95% CI, 1.25-1.49), HLA-DQA1 (rs9271871; P = 1.08 × 10−8; odds ratio, 2.31; 95% CI, 2.02 - 2.60), and TNFRSF11A (rs4263037; P = 1.60 × 10−9; odds ratio, 1.41; 95% CI, 1.29-1.53). These findings replicated for CTLA4 and HLA-DQA1 in an independent cohort of Italian cases and control individuals. Further analysis revealed distinct, but overlapping, disease-associated loci for early- and late-onset forms of myasthenia gravis. In the late-onset cases, we identified 2 association peaks: one was located in TNFRSF11A (rs4263037; P = 1.32 × 10−12; odds ratio, 1.56; 95% CI, 1.44-1.68) and the other was detected in the major histocompatibility complex on chromosome 6p21 (HLA-DQA1; rs9271871; P = 7.02 × 10−18; odds ratio, 4.27; 95% CI, 3.92-4.62). Association within the major histocompatibility complex region was also observed in early-onset cases (HLA-DQA1; rs601006; P = 2.52 × 10−11; odds ratio, 4.0; 95% CI, 3.57-4.43), although the set of single-nucleotide polymorphisms was different from that implicated among late-onset cases.

Conclusions and Relevance  Our genetic data provide insights into aberrant cellular mechanisms responsible for this prototypical autoimmune disorder. They also suggest that clinical trials of immunomodulatory drugs related to CTLA4 and that are already Food and Drug Administration approved as therapies for other autoimmune diseases could be considered for patients with refractory disease.

Introduction

Autoimmune myasthenia gravis is a disorder of neuromuscular transmission clinically characterized by muscle fatigability manifested by diplopia, ptosis, and bulbar and limb weakness.1,2 The disorder is typically mediated by antibodies against nicotinic acetylcholine receptors (AChRs) or against related proteins located at the neuromuscular junction such as muscle-specific tyrosine kinase (MuSK), lipoprotein receptor–related protein 4, and agrin.14 Although myasthenia gravis is relatively uncommon, the apparent incidence has increased in the white population over time owing, at least in part, to improved recognition of the disorder among elderly individuals.5 Acute respiratory failure requiring mechanical ventilation (myasthenic crisis) occurs in up to 20% of patients and is associated with significant morbidity and mortality.6,7

There is increasing recognition that myasthenia gravis is not a monolithic disease.8,9 Epidemiological studies have shown a bimodal pattern of incidence, with early-onset cases (defined as initial symptoms occurring before age 40 years) being predominantly women and late-onset patients being mostly men.911 Advanced age is associated with an increased response to autoantigens, although the implications of the age- and sex-specific frequency distribution of myasthenia gravis regarding pathogenesis remain unclear.10,12

Genetic factors contribute to the susceptibility to develop myasthenia gravis. Identified more than 30 years ago, the human leukocyte antigen (HLA) locus remains the most strongly associated risk factor for the disease.9,13,14 A genome-wide association study (GWAS) involving 649 early-onset cases drawn from the Scandinavian, British, French, Dutch, German, and American populations identified variants in the major histocompatibility complex (MHC) class II locus, protein tyrosine phosphatase nonreceptor type 22 (PTPN22), and TNFAIP3 interacting protein 1 (TNIP1).15 The cytotoxic T-lymphocyte–associated protein 4 gene (CTLA4) has also been previously suggested as a susceptibility factor for myasthenia gravis based on candidate gene studies.16 Furthermore, patients with myasthenia gravis frequently have a personal or family history of other autoimmune disorders, most notably autoimmune thyroid disease, rheumatoid arthritis, and type 1 diabetes mellitus,17,18 although the genetic basis underlying this predisposition to autoimmunity is unknown. Finally, approximately 5% of patients report a family history of myasthenia gravis that usually follows an autosomal dominant pattern of inheritance.14,19

To identify susceptibility loci operating in AChR antibody–positive myasthenia gravis, we undertook a GWAS of 1032 patients diagnosed as having myasthenia gravis (case cohort) and 1998 neurologically normal individuals (control cohort) using Illumina HumanOmniExpress Beadchips.

Methods
Samples

From January 2010 to January 2011, DNA was collected from patients attending myasthenia gravis clinics at 14 centers throughout North America using Oragene DNA Saliva Collection kits (DNA Genotek Inc). All patients included in the case cohort had been diagnosed by a neurologist specializing in myasthenia gravis, had onset of symptoms after age 18 years, and were of non-Hispanic white race/ethnicity (eTable 1 in the Supplement). The diagnosis of myasthenia gravis was based on standard clinical criteria of characteristic weakness, fatigue, and electrophysiological and/or pharmacological abnormalities,1 and it was confirmed by the presence of anti-AChR antibodies.20 Patients with positive test results for antibodies to muscle-specific kinase (anti-MuSK) were excluded from enrollment.21 For the control cohort, we downloaded genotype data from dbGAP (http://www.ncbi.nlm.nih.gov/gap; phs000196.v2.p1) for 1998 US neurologically normal individuals. The control cohort was matched to the case cohort for race/ethnic group but not for age or sex.

The replication cohort consisted of DNA obtained from 423 Italian patients diagnosed as having AChR-positive myasthenia gravis and 467 Italian neurologically normal control individuals (eTable 1 in the Supplement). Blood samples were collected from January 2010 to January 2011 at the Catholic University Rome and at Cisanello Hospital, Pisa, Italy. The control cohort was matched to the case cohort for race/ethnic group but not for age or sex.

Written consent was obtained from all patients enrolled in this study. Institutional review board approval was obtained at all participating institutions including Johns Hopkins University and the National Institute on Aging.

Genome-Wide Genotyping

The case cohort was genotyped in the Laboratory of Neurogenetics, National Institute on Aging, using HumanOmniExpress BeadChips (Illumina), which assay 730 525 single-nucleotide polymorphisms (SNPs) across the genome. Individual-level genotypes for the myasthenia gravis case cohort are available on the dbGAP web portal (phs000726.v1.p1). The control cohort had been previously genotyped at the Center for Inherited Disease Research at Johns Hopkins University on HumanOmni1-Quad BeadChips (Illumina) as part of a GWAS of Parkinson disease.22 Analyses were confined to the 677 673 autosomal SNPs that were common to both chips.

Genotyping Assays in the Replication Cohort

Rs231770, rs4263037, and rs9270986 were analyzed using Taqman genotyping assays and scanned on an ABI 7900HT Real-Time polymerase chain reaction system (Applied Biosystems) according to the manufacturer’s protocol. Rs601006 and rs9271850 in the MHC region were analyzed by sequencing using the Big-Dye Terminator version 3.1 sequencing kit (Applied Biosystems), run on an ABI 3730xl DNA analyzer, and analyzed with Sequencher software (version 4.2; Gene Codes) and Mutation Surveyor (version 4.0.9; Softgenetics). Primers and polymerase chain reaction conditions are listed in eTable 2 in the Supplement.

Genomic DNA from the Italian replication cohort was used for the Taqman assays. Owing to low quantities, the remaining genomic DNA was amplified using the REPLI-g kit according to the manufacturer’s protocol (Qiagen), and the amplified DNA was used in the Sanger sequencing experiments.

Statistical Analysis
Genome-Wide Association

Statistical analyses were performed using R statistical software (version 2.15.2).23 We applied standard quality-control procedures to the data: exclusion of samples with SNP call rates of less than 95%, non-European ancestry, or cryptic relatedness defined as identity-by-descent proportion of inheritance (pi_hat from the PLINK software toolset version 1.0624) greater than 0.1, as well as exclusion of SNPs with call rates of less than 95%, minor allele frequency less than 0.01 in the control cohort, and Hardy-Weinberg equilibrium P < .001 in the control cohort. The cryptic-relatedness threshold led to the exclusion of individuals who shared more than 10% of their genome, which meant that related individuals down to third- or fourth-degree relatives were not included in the final analysis. The index individual whose sample had the better call rate from each related pair was excluded from the analysis (17 related pairs).

After quality control, we used a Markov chain–based Haplotyper (version 1.0.16) to impute genotypes for all participants.25 A total of 8 114 394 SNPs (consisting of 513 081 genotyped SNPs and 7 601 313 imputed SNPs) were available for analysis.

We calculated P values using logistic regression modeling incorporating the first 2 principal components as covariates to compensate for any residual population stratification. Principal components were generated using Genome-wide Complex Trait Analysis software package implementation of eigenstrat.26 A threshold of 5.0 × 10−8 was set for genome-wide significance after Bonferroni correction for multiple testing.27

Probability Analysis and Heritability Estimates

We used density estimation to generate posterior probabilities of developing myasthenia gravis based on age and sex.

The Genome-wide Complex Trait Analysis package uses covariance matrices and mixed modeling to estimate the heritability of a trait in an ostensibly outbred population.26,28,29 We applied this method to compare each case series (all cases, early-onset cases, and late-onset cases) to all control individuals. For each model, we compared 2 separate subsets of SNPs, one inclusive of all genotyped SNPs and another including only SNPs within ±1 MB from the loci identified as genome-wide significant in the discovery phase of analyses (Table 1). Only SNPs and samples passing quality control were used to evaluate the heritability and all analyses were adjusted for population substructure by using C1 and C2 from principle component analysis as covariates. Further details are provided in the eAppendix in the Supplement.

Results
GWAS Identifies Loci Influencing Risk for Myasthenia Gravis

After imputation and quality-control measures, 8 114 394 genotyped and imputed SNPs from 972 North American patients diagnosed as having AChR antibody–positive myasthenia gravis and 1977 North American control individuals were available for analysis. Quartile-quartile plots did not show evidence of significant population stratification (Figure 1A; genomic inflation factor λ = 1.036).

In addition to the previously reported MHC class II region signal near the MHC class II DQ α 1 (HLA-DQA1) gene,30 we observed a strong association signal on chromosome 2q33 in the CTLA4 locus with a P value that reached genome-wide significance (rs231770; P = 3.98 × 10−8; odds ratio [OR], 1.37; 95% CI, 1.25-1.49; Figure 1B). We also identified a strong association peak on chromosome 18q21.33 within the tumor necrosis factor receptor 4 superfamily, member 11a, NFKB activator (TNFRSF11A) gene (rs4263037; P = 1.60 × 10−9;OR, 1.41; 95% CI, 1.29-1.53). Markers that achieved genome-wide significance in the analysis of the whole cohort are shown in Table 1, and SNPs associated at P < 1.0 × 10−6 are listed in eTable 3 in the Supplement.

Early- and Late-Onset Myasthenia Gravis Possess Distinct Genetic Architecture

Density estimate plots confirmed a bimodal sex distribution within our case cohort, with myasthenia gravis being more common in younger women and older men (Figure 1C). To elucidate the genetic architecture underlying this age difference, we performed GWAS of early-onset cases (235 cases with onset at age <40 years and 1977 control individuals) and late-onset cases (737 cases with onset at age >40 years and 1977 control individuals). The TNFRSF11A and HLA-DQA1 association signals were significantly enhanced among late-onset cases (rs4263037, P = 1.32 × 10−12; OR, 1.56; 95% CI, 1.44-1.68 and rs9271871, P = 7.02 × 10−18; OR, 4.27; 95% CI, 3.92-4.62, respectively; Figure 1D; Table 1). In contrast, there was no evidence of association within TNFRSF11A in the young-onset cases (rs4263037; P = .83; Figure 1E). Furthermore, although an association signal was present near HLA-DQA1 in the younger patients, the most significantly associated SNPs were different compared with those observed in the late-onset cohort (rs601006; P = 2.52 × 10−11; OR, 4.0; 95% CI, 3.57-4.43; eTables 4 and 5 in the Supplement).

Myasthenia Gravis Risk Loci Replicated in an Independent Cohort

We selected 3 SNPs from the risk loci identified on chromosomes 2, 6, and 18 in the overall cohort for genotyping in an independent replication collection consisting of 423 Italian AChR antibody–positive myasthenia gravis cases and 467 Italian control individuals. We found the strongest signals at rs9270986 in the intergenic region between HLA-DRB1 and HLA-DQA1 and at rs231770 located 3.3 kb upstream of CTLA4 (Figure 2A and B). In a combined analysis, we integrated the US GWAS and the Italian replication data using fixed effects modeling. Two loci (CTLA4 and HLA-DQA1) were associated with myasthenia gravis at genome-wide significance (P = 9.25 × 10−11 and 6.17 × 10−8, respectively; Table 1).

We also attempted to replicate the most highly associated SNPs observed in the early- and late-onset cases using the same Italian replication collection. Rs601006 near HLA-DQA1 was significantly associated among early-onset cases (169 Italian myasthenia gravis cases and 467 control individuals; P = 5.53 × 10−5; Table 1; Figure 2C). This SNP was not associated among late-onset cases (250 cases and 467 control individuals; P = .22), recapitulating the pattern observed in the discovery cohort.

Similarly, we replicated the SNPs observed near HLA-DQA1 and within TNFRSF11A in the late-onset cases. Rs9271850 near HLA-DQA1 was significantly associated in the late-onset Italian cohort (P = 8.15 × 10−4; Table 1; Figure 2D) but was not associated among the early-onset cases (P = .53). Rs4263037 located on chromosome 18 showed a trend toward significance (P for patients aged ≥40 years at symptom onset = .09; Figure 2E). This SNP became progressively more associated with increasing age of the replication cohort (P for patients aged ≥50 years at symptom onset = .04; P for patients aged ≥ 70 years = 8.75 × 10−3).

Myasthenia Gravis Is Highly Heritable

The relative importance of genetic factors in a complex disease like myasthenia gravis can be estimated by quantifying heritability. We generated heritability estimates for AChR antibody–positive myasthenia gravis using a sophisticated algorithm (Genome-wide Complex Trait Analysis) that analyzes all of the genome-wide SNPs simultaneously (Table 2). Our genotype data accounted for 25.6% (95% CI, 18.6-32.6) of the phenotype variance observed in all myasthenia gravis cases. Analysis of more precise phenotypes (early- and late-onset cases) generated even higher estimates of heritability (37.9%; 95% CI, 16.8-59.0 and 35.3%; 95% CI, 27.1-43.5, respectively). These estimates are substantially higher than heritability identified using genome-wide–associated SNPs alone (1%-3%), indicating that additional risk loci remain to be identified.

Discussion

Our GWAS of AChR antibody–positive myasthenia gravis identified several new loci that drive susceptibility to develop disease. We also provide strong evidence that early- and late-onset myasthenia gravis have distinct, but overlapping, genetic architecture, thereby confirming previous epidemiological and clinical data suggesting a division between these 2 forms of disease. Genetic variation within the TNFSRF11A locus drives susceptibility to disease only among older cases, whereas different haplotypes across the same HLA region on chromosome 6 were identified in early- and late-onset cases. CTLA4 exerts a significant effect regardless of age at symptom onset, suggesting that it plays a central role in generating the aberrant autoimmune response that leads to neuromuscular junction dysfunction.

CTLA4 is a 45-kD immunoglobulin protein expressed by activated T cells and sharing significant sequence identity with CD28 (Figure 2F). CTLA4 increases T-cell motility and reduces contact periods between T cells and antigen-presenting cells leading to decreased cytokine production and proliferation. In this way, CTLA4 is thought to downregulate T-cell activation, terminate T-cell responses, and protect against autoimmunity.31 CTLA4-deficient mice have high serum immunoglobulin concentrations and develop massive lymphocyte infiltration of lymph nodes and spleen, suggesting that CTLA4 plays an inhibitory role in regulating lymphocyte expansion.32 The mechanism of CTLA4’s regulatory action is controversial; however, evidence suggests that cells expressing CTLA4 avidly endocytose the costimulatory factors CD80 and CD86 on antigen-presenting cells and degrade them.33 Genetic variant in the CTLA4 locus has been implicated in other autoimmune disorders including celiac disease,34 type 1 diabetes mellitus,35 thyroiditis,36 and rheumatoid arthritis.37 Indeed, abatacept and belatacept are commercially available proteins consisting of CTLA4 rendered soluble by fusion to antibodies (CTLA4-IgG) approved for use in the treatment of refractory rheumatoid arthritis38 and as first-line immunosuppressants for renal transplant patients immune to Epstein-Barr virus.39 Clinical trials of these agents are also under way in inflammatory bowel disease, systemic lupus erythematosus, and other autoimmune diseases. They have not previously been used therapeutically in myasthenia gravis, although our experimental studies demonstrated the effectiveness of CTLA4-Ig in the treatment of rats with experimental autoimmune myasthenia gravis.40,41

Genetic analysis of additional case-control cohorts will be required to confirm the association observed for the TNFRSF11A locus. Nevertheless, the pattern of association uncovered for this locus, with a large effect in late-onset disease that is absent in younger patients, led us to speculate that age-related changes in local gene expression predispose toward an aberrant immune response to autoantigens. TNFRSF11A encodes the 4.5-kDa receptor activator of nuclear factor-κ B expressed on the surface of antigen-presenting dendritic cells.42 This receptor is an important regulator of the interaction between T cells and dendritic cells that is essential for immune surveillance and the regulation of specific immunity.42 Furthermore, TNFRSF11A is critical for lymph node organogenesis and osteoclast differentiation.43 Mutations in this gene are responsible for a form of familial Paget disease of the bone44 and for autosomal recessive osteopetrosis associated with defective immunoglobulin production.45 Receptor activator of nuclear factor-κ B knockout mice also have profound osteopetrosis resulting from a complete lack of osteoclasts, as well as near total absence of peripheral lymph nodes.43,45

The association signal that we identified in the MHC region was distinct from that reported by Gregersen et al15 to increase susceptibility for myasthenia gravis. Furthermore, our study did not show any evidence of association in the PTPN22 or TNIP1 loci found by the same study (eTable 6 in the Supplement).15 These disparities likely arise from differences in the populations studied. In their study, patients and control individuals from 7 countries were included.15 Such population stratification impedes the identification of true association and increases the rate of false-positive association46 particularly in the HLA region that is known to be highly divergent across races/ethnicities.47 In contrast, we attempted to reduce clinical heterogeneity by selecting only cases from North America for the discovery cohort and from Italy for the replication cohort.

We considered the possibility that our study may have had population stratification, especially because we used an outbred American population as our original discovery cohort. The lack of age and sex matching of control samples to cases also may have introduced artifacts into our analysis. To avoid this, we applied strict filters to exclude outlying samples and incorporated principal components generated from the genome-wide data into the analysis model to compensate for residual population stratification between case and control samples. Furthermore, we replicated our results in a case-control cohort drawn from a distinct single population and found identical patterns of association among early- and late-onset cases.

Heritability can be considered the proportion of a disease that arises from genetics. Prior studies have estimated that only between 3.8% and 7.1% of individuals who receive a diagnosis of myasthenia gravis have a family history of disease.14,48 The genome-wide nature of our data set yielded a higher estimate of the heritability associated with myasthenia gravis risk by capturing more of the polygenic variance. This pattern is consistent with other neurological diseases, such as amyotrophic lateral sclerosis and Parkinson disease, where heritability estimates are significantly higher than the rate of familial disease.49,50 Furthermore, our heritability estimates suggest that there are genetic factors contributing to the disease that have yet to be identified and that GWASs involving larger cohorts will yield additional loci. To that end, we have made our raw genotype data publicly available on dbGaP to facilitate meta-analysis with existing and future studies.

A strength of our study was that genetic analysis was confined to patients diagnosed as having AChR antibody–positive myasthenia gravis. This approach decreased disease heterogeneity and improved the power of the study to identify association signals. Approximately 85% of patients with myasthenia gravis have detectable AChR antibodies and approximately 40% of the remainder are MuSK antibody positive.2,51,52 Furthermore, AChR antibody–negative myasthenia gravis is less common among elderly individuals,52 and MuSK-positive patients are more like to be women.53 In light of this, care should be taken in generalizing our results to the AChR antibody–negative myasthenia gravis population. Nevertheless, our data support the existence of distinct, but overlapping, genetic profiles for patients with early- and late-onset myasthenia gravis with detectable AChR antibodies.

Conclusions

Our GWAS identified susceptibility loci for AChR antibody–positive myasthenia gravis at CTLA4, HLA-DRB1/HLA-DQA1, and TNFRSF11A, and it showed distinct, but overlapping, genetic risk factors underlying early- and late-onset disease. Although future studies involving deep resequencing of these loci will be required to identify the functional variants underlying these association signals, our data also suggest that clinical trials of CTLA4 immunomodulatory therapies could be considered in patients with myasthenia gravis with refractory disease failing to respond to standard therapies.

Back to top
Article Information

Corresponding Author: Bryan J. Traynor, MD, PhD, Neuromuscular Diseases Research Unit, Laboratory of Neurogenetics, Porter Neuroscience Research Center, 35 Convent Dr, Bethesda, MD 20892 (traynorb@mail.nih.gov).

Accepted for Publication: November 6, 2014.

Published Online: February 2, 2015. doi:10.1001/jamaneurol.2014.4103.

Author Contributions: Drs Drachman and Traynor had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. Drs Renton, Provenzano, Evoli, Ricciardi, Drachman, and Traynor, and Ms Pliner contributed equally to this work.

Study concept and design: Drachman, Traynor.

Acquisition, analysis, or interpretation of data: All authors.

Drafting of the manuscript: Drachman, Traynor.

Critical revision of the manuscript for important intellectual content: All authors.

Statistical analysis: Renton, Nalls, Gibbs, Drachman, Traynor.

Obtained funding: Drachman, Traynor.

Administrative, technical, or material support: Provenzano, Evoli, Ricciardi, Hernandez, Bartoccioni, Scuderi, Maestri, Errichiello, Macek, Scholz, Corse, Benatar, Barohn, Pasnoor, Rowin, Kissel, Massey, Chopra, Koopman, Nicolle, Pestronk, Wulf, Florence, Soloway, Wolfe, Jiwa, Oger, Drachman, Traynor.

Study supervision: Drachman, Traynor.

Conflict of Interest Disclosures: None reported.

Funding/Support: The work was supported by the Myasthenia Gravis Foundation (Drs Drachman and Traynor), a generous bequest by Geraldine Weinrib, and a gift from Philip Swift. This work was also supported in part by the Intramural Research Programs of the National Institutes of Health National Institute on Aging (Z01-AG000949-02).

Role of the Funder/Sponsor: The funders had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Additional Contributions: We thank the patients and research participants who contributed samples for this study.

Additional Information: This study used the high-performance computational capabilities of the Biowulf Linux cluster at the National Institutes of Health, Bethesda, Maryland (http://biowulf.nih.gov).

References
1.
Drachman  DB.  Myasthenia gravis. N Engl J Med. 1994;330(25):1797-1810.
PubMedArticle
2.
Spillane  J, Higham  E, Kullmann  DM.  Myasthenia gravis. BMJ. 2012;345:e8497.
PubMedArticle
3.
Higuchi  O, Hamuro  J, Motomura  M, Yamanashi  Y.  Autoantibodies to low-density lipoprotein receptor-related protein 4 in myasthenia gravis. Ann Neurol. 2011;69(2):418-422.
PubMedArticle
4.
Zhang  B, Shen  C, Bealmear  B,  et al.  Autoantibodies to agrin in myasthenia gravis patients. PLoS One. 2014;9(3):e91816.
PubMedArticle
5.
Carr  AS, Cardwell  CR, McCarron  PO, McConville  J.  A systematic review of population based epidemiological studies in Myasthenia Gravis. BMC Neurol. 2010;10:46.
PubMedArticle
6.
Alshekhlee  A, Miles  JD, Katirji  B, Preston  DC, Kaminski  HJ.  Incidence and mortality rates of myasthenia gravis and myasthenic crisis in US hospitals. Neurology. 2009;72(18):1548-1554.
PubMedArticle
7.
Jani-Acsadi  A, Lisak  RP.  Myasthenic crisis: guidelines for prevention and treatment. J Neurol Sci. 2007;261(1-2):127-133.
PubMedArticle
8.
Cavalcante  P, Bernasconi  P, Mantegazza  R.  Autoimmune mechanisms in myasthenia gravis. Curr Opin Neurol. 2012;25(5):621-629.
PubMedArticle
9.
Vincent  A, Palace  J, Hilton-Jones  D.  Myasthenia gravis. Lancet. 2001;357(9274):2122-2128.
PubMedArticle
10.
Sanders  DB, Andrews  I, Howard  JF, Massey  JM.  Seronegative myasthenia gravis. Neurology. 1997;48(suppl 5):S40-S45. doi:10.1212/WNL.48.Suppl_5.40S.
PubMedArticle
11.
Somnier  FE.  Myasthenia gravis. Dan Med Bull. 1996;43(1):1-10.
PubMed
12.
LeMaoult  J, Szabo  P, Weksler  ME.  Effect of age on humoral immunity, selection of the B-cell repertoire and B-cell development. Immunol Rev. 1997;160:115-126.
PubMedArticle
13.
Carlsson  B, Wallin  J, Pirskanen  R, Matell  G, Smith  CI.  Different HLA DR-DQ associations in subgroups of idiopathic myasthenia gravis. Immunogenetics. 1990;31(5-6):285-290.
PubMedArticle
14.
Pirskanen  R.  Genetic aspects in myasthenia gravis: a family study of 264 Finnish patients. Acta Neurol Scand. 1977;56(5):365-388.
PubMedArticle
15.
Gregersen  PK, Kosoy  R, Lee  AT,  et al.  Risk for myasthenia gravis maps to a (151) Pro→Ala change in TNIP1 and to human leukocyte antigen-B*08. Ann Neurol. 2012;72(6):927-935.
PubMedArticle
16.
Wang  XB, Pirskanen  R, Giscombe  R, Lefvert  AK.  Two SNPs in the promoter region of the CTLA-4 gene affect binding of transcription factors and are associated with human myasthenia gravis. J Intern Med. 2008;263(1):61-69.
PubMedArticle
17.
Christensen  PB, Jensen  TS, Tsiropoulos  I,  et al.  Associated autoimmune diseases in myasthenia gravis: a population-based study. Acta Neurol Scand. 1995;91(3):192-195.
PubMedArticle
18.
Tsao  CY, Mendell  JR, Lo  WD, Luquette  M, Rennebohm  R.  Myasthenia gravis and associated autoimmune diseases in children. J Child Neurol. 2000;15(11):767-769.
PubMedArticle
19.
Szobor  A.  Myasthenia gravis: familial occurrence: a study of 1100 myasthenia gravis patients. Acta Med Hung. 1989;46(1):13-21.
PubMed
20.
Padua  L, Stalberg  E, LoMonaco  M, Evoli  A, Batocchi  A, Tonali  P.  SFEMG in ocular myasthenia gravis diagnosis. Clin Neurophysiol. 2000;111(7):1203-1207.
PubMedArticle
21.
Hoch  W, McConville  J, Helms  S, Newsom-Davis  J, Melms  A, Vincent  A.  Auto-antibodies to the receptor tyrosine kinase MuSK in patients with myasthenia gravis without acetylcholine receptor antibodies. Nat Med. 2001;7(3):365-368.
PubMedArticle
22.
Hamza  TH, Zabetian  CP, Tenesa  A,  et al.  Common genetic variation in the HLA region is associated with late-onset sporadic Parkinson’s disease. Nat Genet. 2010;42(9):781-785.
PubMedArticle
23.
R Core Team. A language and environment for statistical computing. 2012. http://www.R-project.org/. Accessed September 30, 2013.
24.
Purcell  S, Neale  B, Todd-Brown  K,  et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559-575.
PubMedArticle
25.
Li  Y, Willer  C, Sanna  S, Abecasis  G.  Genotype imputation. Annu Rev Genomics Hum Genet. 2009;10:387-406.
PubMedArticle
26.
Yang  J, Benyamin  B, McEvoy  BP,  et al.  Common SNPs explain a large proportion of the heritability for human height. Nat Genet. 2010;42(7):565-569.
PubMedArticle
27.
Pe’er  I, Yelensky  R, Altshuler  D, Daly  MJ.  Estimation of the multiple testing burden for genomewide association studies of nearly all common variants. Genet Epidemiol. 2008;32(4):381-385.
PubMedArticle
28.
Lee  SH, Wray  NR, Goddard  ME, Visscher  PM.  Estimating missing heritability for disease from genome-wide association studies. Am J Hum Genet. 2011;88(3):294-305.
PubMedArticle
29.
Yang  J, Manolio  TA, Pasquale  LR,  et al.  Genome partitioning of genetic variation for complex traits using common SNPs. Nat Genet. 2011;43(6):519-525.
PubMedArticle
30.
Djabiri  F, Caillat-Zucman  S, Gajdos  P,  et al.  Association of the AChRalpha-subunit gene (CHRNA), DQA1*0101, and the DR3 haplotype in myasthenia gravis. Evidence for a three-gene disease model in a subgroup of patients. J Autoimmun. 1997;10(4):407-413.
PubMedArticle
31.
Schneider  H, Downey  J, Smith  A,  et al.  Reversal of the TCR stop signal by CTLA-4. Science. 2006;313(5795):1972-1975.
PubMedArticle
32.
Waterhouse  P, Penninger  JM, Timms  E,  et al.  Lymphoproliferative disorders with early lethality in mice deficient in Ctla-4. Science. 1995;270(5238):985-988.
PubMedArticle
33.
Qureshi  OS, Zheng  Y, Nakamura  K,  et al.  Trans-endocytosis of CD80 and CD86: a molecular basis for the cell-extrinsic function of CTLA-4. Science. 2011;332(6029):600-603.
PubMedArticle
34.
Dubois  PC, Trynka  G, Franke  L,  et al.  Multiple common variants for celiac disease influencing immune gene expression. Nat Genet. 2010;42(4):295-302.
PubMedArticle
35.
Barrett  JC, Clayton  DG, Concannon  P,  et al; Type 1 Diabetes Genetics Consortium.  Genome-wide association study and meta-analysis find that over 40 loci affect risk of type 1 diabetes. Nat Genet. 2009;41(6):703-707.
PubMedArticle
36.
Chu  X, Pan  CM, Zhao  SX,  et al; China Consortium for Genetics of Autoimmune Thyroid Disease.  A genome-wide association study identifies two new risk loci for Graves’ disease. Nat Genet. 2011;43(9):897-901.
PubMedArticle
37.
Gregersen  PK, Amos  CI, Lee  AT,  et al.  REL, encoding a member of the NF-kappaB family of transcription factors, is a newly defined risk locus for rheumatoid arthritis. Nat Genet. 2009;41(7):820-823.
PubMedArticle
38.
Schiff  M, Keiserman  M, Codding  C,  et al.  Efficacy and safety of abatacept or infliximab vs placebo in ATTEST: a phase III, multi-centre, randomised, double-blind, placebo-controlled study in patients with rheumatoid arthritis and an inadequate response to methotrexate. Ann Rheum Dis. 2008;67(8):1096-1103.
PubMedArticle
39.
Vincenti  F, Charpentier  B, Vanrenterghem  Y,  et al.  A phase III study of belatacept-based immunosuppression regimens versus cyclosporine in renal transplant recipients (BENEFIT Study). Am J Transplant. 2010;10(3):535-546.
PubMedArticle
40.
McIntosh  KR, Linsley  PS, Bacha  PA, Drachman  DB.  Immunotherapy of experimental autoimmune myasthenia gravis: selective effects of CTLA4Ig and synergistic combination with an IL2-diphtheria toxin fusion protein. J Neuroimmunol. 1998;87(1-2):136-146.
PubMedArticle
41.
McIntosh  KR, Linsley  PS, Drachman  DB.  Immunosuppression and induction of anergy by CTLA4Ig in vitro: effects on cellular and antibody responses of lymphocytes from rats with experimental autoimmune myasthenia gravis. Cell Immunol. 1995;166(1):103-112.
PubMedArticle
42.
Anderson  DM, Maraskovsky  E, Billingsley  WL,  et al.  A homologue of the TNF receptor and its ligand enhance T-cell growth and dendritic-cell function. Nature. 1997;390(6656):175-179.
PubMedArticle
43.
Dougall  WC, Glaccum  M, Charrier  K,  et al.  RANK is essential for osteoclast and lymph node development. Genes Dev. 1999;13(18):2412-2424.
PubMedArticle
44.
Hughes  AE, Shearman  AM, Weber  JL,  et al.  Genetic linkage of familial expansile osteolysis to chromosome 18q. Hum Mol Genet. 1994;3(2):359-361.
PubMedArticle
45.
Guerrini  MM, Sobacchi  C, Cassani  B,  et al.  Human osteoclast-poor osteopetrosis with hypogammaglobulinemia due to TNFRSF11A (RANK) mutations. Am J Hum Genet. 2008;83(1):64-76.
PubMedArticle
46.
Freedman  ML, Reich  D, Penney  KL,  et al.  Assessing the impact of population stratification on genetic association studies. Nat Genet. 2004;36(4):388-393.
PubMedArticle
47.
Buhler  S, Sanchez-Mazas  A.  HLA DNA sequence variation among human populations: molecular signatures of demographic and selective events. PLoS One. 2011;6(2):e14643.
PubMedArticle
48.
Namba  T, Brunner  NG, Brown  SB, Muguruma  M, Grob  D.  Familial myasthenia gravis: report of 27 patients in 12 families and review of 164 patients in 73 families. Arch Neurol. 1971;25(1):49-60.
PubMedArticle
49.
Keller  MF, Ferrucci  L, Singleton  AB,  et al.  Genome-wide analysis of the heritability of amyotrophic lateral sclerosis. JAMA Neurol. 2014;71(9):1123-1134.
PubMedArticle
50.
Keller  MF, Saad  M, Bras  J,  et al; International Parkinson’s Disease Genomics Consortium (IPDGC); Wellcome Trust Case Control Consortium 2 (WTCCC2).  Using genome-wide complex trait analysis to quantify ‘missing heritability’ in Parkinson’s disease. Hum Mol Genet. 2012;21(22):4996-5009.
PubMedArticle
51.
Chan  KH, Lachance  DH, Harper  CM, Lennon  VA.  Frequency of seronegativity in adult-acquired generalized myasthenia gravis. Muscle Nerve. 2007;36(5):651-658.
PubMedArticle
52.
Aarli  JA.  Myasthenia gravis in the elderly: is it different? Ann N Y Acad Sci. 2008;1132:238-243.
PubMedArticle
53.
Evoli  A, Tonali  PA, Padua  L,  et al.  Clinical correlates with anti-MuSK antibodies in generalized seronegative myasthenia gravis. Brain. 2003;126(pt 10):2304-2311.
PubMedArticle
×