Research article | Open | Published:
Systematic protein-protein interaction and pathway analyses in the idiopathic inflammatory myopathies
Arthritis Research & Therapyvolume 18, Article number: 156 (2016)
The idiopathic inflammatory myopathies (IIM) are autoimmune diseases characterised by acquired proximal muscle weakness, inflammatory cell infiltrates in muscle and myositis-specific/associated autoantibodies. It is unclear which pathways are involved in IIM, and the functional relationship between autoantibody targets has not been systematically explored. Protein-protein interaction and pathway analyses were conducted to identify pathways relevant to disease, using autoantibody targets and gene products of IIM-associated single nucleotide polymorphism (SNP) loci.
Protein-protein interactions were analysed using Disease Association Protein-Protein Link Evaluator (DAPPLE). Gene ontology and pathway analyses were conducted using Database for Annotation Visualisation and Integrated Discovery (DAVID) and Gene Relationships Across Implicated Loci (GRAIL). Analyses were undertaken including the targets of published autoantibodies, significant and suggestive SNPs from an IIM association study and autoantibody targets plus SNPs combined.
The protein-protein interaction networks formed by autoantibody targets and associated SNPs showed significant direct and/or indirect connectivity (p < 0.05). Autoantibody targets plus associated SNPs combined resulted in more significant indirect and common interactor connectivity, suggesting autoantibody targets and proteins encoded by IIM-associated loci may be involved in common pathways. Tumour necrosis factor receptor-associated factor 6 (TRAF6) was identified as a hub protein, and UBE3B, HSPA1A, HSPA1B and PSMD3 also were identified as genes with significant connectivity. Pathway analysis identified that autoantibody targets and associated SNP regions are significantly interconnected (p < 0.01), and confirmed autoantibody target involvement in translational and post-translational processes. ‘Ubiquitin’ was the only keyword strongly linking significant genes across regions in all three GRAIL analyses of autoantibody targets and IIM-associated SNPs.
Autoantibody targets and IIM-associated loci show significant connectivity and inter-relatedness, and identify several key genes and pathways in IIM pathogenesis, possibly mediated via the ubiquitination pathway.
The idiopathic inflammatory myopathies (IIM) are a group of rare autoimmune diseases whereby patients typically present with weakness and inflammation of skeletal muscle, and extra-muscular manifestations including interstitial lung disease, malignancy and skin rashes . Approximately 80 % of patients with IIM possess serum autoantibodies; these may be myositis-specific autoantibodies (MSAs), which are found predominantly in IIM, or myositis-associated autoantibodies (MAAs), which are found also in other connective tissue diseases (CTD) or as part of a myositis/CTD-overlap condition. MSAs are associated with a particular clinical profile and are usually mutually exclusive .
The most common MSA is anti-Jo-1 which is found in about 20 % of patients with IIM. This MSA targets histidyl tRNA synthetase, one of eight aminoacyl tRNA synthetases (ARS) known to act as autoantigens for MSAs . Anti-ARS are associated with a distinct phenotype known as antisynthetase syndrome, which involves myopathy, non-erosive arthritis, interstitial lung disease, mechanic’s hands and fever .
Modification of the structure of an autoantigen, for example, by somatic mutation or post-translational modification, in a pro-immune context such as muscle injury may result in an autoimmune response . In IIM, regenerating muscle cells produce higher levels of autoantigens such as ARS molecules . There is some evidence that ARS molecules may act as chemoattractants inducing migration of immune cells to the affected muscle . Alongside elevated expression of major histocompatibility complex (MHC)-class I on cell surfaces, this may lead to a specific autoimmune response resulting in further muscle damage and regeneration; thereby more autoantigen is produced, resulting in a sustained immune response, disease propagation and amplification .
There are currently no disease-specific treatments for IIM; instead anti-inflammatory, immunosuppressive and immunodulatory treatments are borrowed from other rheumatic diseases. Furthermore, despite treatment some patients develop progressive muscle weakness and extramuscular manifestions with ensuing disability and poor quality of life. There is thus a clear need for new targeted treatments in IIM. It is currently unclear which pathways are involved in IIM and whether autoantibodies contribute to the disease process. The functional relationship between autoantibody targets has also not been systematically explored. Protein-protein interaction (PPI) analyses have been performed in a number of complex diseases, including breast cancer, schizophrenia and multiple sclerosis [5–8]. Genome-wide association studies (GWAS) identify single nucleotide polymorphisms (SNPs) associated with particular traits or diseases and PPI analyses have been used to identify the genes most likely to be functionally significant to pathological change [6, 7]; the hypothesis is that associated genes will be involved in a common set of biological pathways or processes, which may be perturbed in the disease. It is possible that the targeting of autoantigens in an autoimmune disease setting like IIM also could be related to these altered pathways.
In this study we analysed PPI and pathways using autoantibody targets and proteins encoded by loci associated with IIM in the largest genetic association study to date , to identify pathways involved in IIM. This is the first time a systematic PPI and pathway analysis has been performed in IIM and suggests particular proteins and pathways that may have important roles in disease.
Protein-protein interaction analysis
Analysis was performed using Disease Association Protein-Protein Link Evaluator (DAPPLE) to investigate physical connections between proteins . DAPPLE searches the InWeb database for PPI that have been reported in the literature and assigned a probabilistic score. The InWeb database compiles PPI data from numerous sources including Reactome, IntAct, the Molecular Interaction Database (MINT), the Biomolecular Interaction Network Database (BIND) and the Kyoto Encyclopaedia of Genes and Genomes (KEGG) . DAPPLE is designed to analyse disease-associated single nucleotide polymorphisms (SNPs) on the basis that disease-causing genetic variation is likely to affect common pathways that may be revealed by PPI .
Based on these interactions, DAPPLE forms networks of physical protein-protein connectivity where proteins are nodes connected by edges that represent interactions in the InWeb database. A direct connection may be present between two seed proteins (proteins derived from the input SNPs or genes, also referred to as ‘associated proteins’) or an indirect connection may be formed via a common interactor protein (a protein not derived from the input but which interacts with two or more seed proteins). Networks between proteins are evaluated based on four parameters: (1) number of direct interactions (number of edges in the direct network); (2) mean associated protein direct connectivity (the average number of proteins with which each seed protein directly interacts); (3) mean associated protein indirect connectivity (the average number of proteins with which each seed protein indirectly interacts); (4) mean common interactor connectivity (the average number of seed proteins bound by common interactor proteins) (Fig. 1).
The statistical significance of each parameter is assessed using 20,000 permutations, comparing the network formed against networks created by randomly reassigning proteins of the same binding degree as the protein in the original network to each node (binding degree was defined as the total number of interactions a protein has according to the InWeb database). In a random list of proteins any connectivity could simply be a function of the binding degree of each protein. This method controls for bias introduced by the differing extent to which certain PPI have been studied or behave in vitro.
DAPPLE accepts named genes or SNPs as input. From SNP input genes are identified as those in linkage disequilibrium (LD) (r 2 > 0.5), extended to the closest recombination hotspots plus 50 kb. Protein products of genes in associated loci are then scored based on their participation in direct and indirect networks. These scores are Bonferroni corrected for two tests if a protein participates in both networks and further Bonferroni corrected for the number of possible candidates in that locus to identify the ‘genes to prioritise’ (p corr < 0.05). . Analysis was performed using DAPPLE, version 2, Hg18, and HapMap SNPs.
Parallel pathway and gene ontology analyses were conducted using Gene Relationships Across Implicated Loci (GRAIL)  and the Database for Annotation, Visualisation and Integrated Discovery (DAVID) [12, 13]. GRAIL uses an automated text-based strategy to identify significantly inter-related genes based on a list of disease- associated regions, and identifies the most likely candidate gene in each region . Candidate genes within each region are selected based on a significance score for their relationships to genes in other associated regions, corrected for multiple hypothesis testing for the number of genes in the region. Candidate genes with uncorrected GRAIL text scores (p text < 0.2) are used to identify the top twenty highest scoring keywords based on the average frequency of that term for candidate genes versus all genes. These keywords must appear in at least 500 documents, be four or more letters long and contain no numbers . GRAIL makes no assumptions about the phenotype being studied or underlying pathways presumed to be relevant to the disease. Analysis was conducted using gene size correction and the October 2014 text database.
DAVID analysis was conducted using biological processes from the Gene ontology consortium (GOTERM_BP_FAT) and pathways from KEGG using DAVID version 6.7. P values were Benjamini-Hochberg corrected for multiple hypothesis testing.
Input for protein-protein interaction and pathway analyses
Analyses were performed using the following inputs: (1) all published MSA targets and anti-PMScl-75/100 targets; (2) the most significantly associated SNPs from each region from a recently reported IIM Immunochip association study ; (3) the combination of SNPs and autoantibody targets (Table 1; these data are also available in Additional file 1). Anti-PMScl-75/100 was included despite being classified as an MAA, as it is almost always mutually exclusive to other MSA/MAAs.
The Uniprot IDs of all published MSA and anti-PMScl-75/100 targets were identified (n = 28). SNPs reaching GWAS thresholds for significance (p < 5 × 10-8) or suggestive significance (p < 2.25 × 10-5, determined using the genetic type I error calculator ) in the IIM Immunochip association study  were identified from each clinical subgroup (combined IIM, polymyositis, adult/juvenile dermatomyositis) (n = 22). The genetic type 1 error calculator estimates the effective number of independent tests based on the LD between SNPs contained on the genotyping array. The most significantly associated SNP was used for each region, or a proxy was used if this SNP was not recognised by DAPPLE (proxy SNP in LD r 2 ≥ 0.9 with associated SNP, identified using SNAP version 2.2, hg18, HapMap 22 ). Input genes (encoding proteins that may participate in these networks) were identified based on LD (r 2 > 0.5 with the associated SNP) and extension to the closest recombination hotspot for GRAIL, or extension to the closest recombination hotspot +50 kb for DAPPLE , resulting in 185 and 151 genes used as input for GRAIL and DAPPLE, respectively. For DAPPLE the analysis was also run using a subset of autoantibody targets selecting one protein from each complex (n = 21), to remove bias from known interactions; this subset was used in the combined autoantibody targets-SNP analysis (n = 43). The genes from associated SNP regions identified by DAPPLE were used as input for DAVID; MHC region SNPs were excluded from the DAVID analysis to remove bias introduced by the strong LD and large number of immune-related genes across this region, resulting in an input of 62 genes from 20 SNPs plus the 28 autoantibody targets (n = 90).
Protein-protein interaction analysis identifies significant direct and indirect connectivity
The PPI networks formed by autoantibody targets and IIM-associated SNPs using DAPPLE had significant direct and/or indirect connectivity (p corr < 0.05, Fig. 2, Additional file 2), including when only one protein was selected from each known autoantibody target complex. These results indicate that these proteins interact more than a group of randomly selected proteins in 20,000 permutations of the networks and therefore could be involved in common molecular networks linked to IIM. Due to strong direct PPI, 17/28 autoantibody targets were identified as genes to prioritize as they participate in the network more than would be expected by chance (p corr values <0.05, Table 2). SNP loci alone had significant indirect and common interactor connectivity but direct connectivity was not significant, suggesting that the proteins encoded by associated loci do not directly interact more than a randomly selected group of proteins but may be involved in common pathways (indirect networks and common interactor p values may be found in Additional file 3). The SNP analysis identified UBE3B (p corr = 8.03 × 10-4) and TRAF6 (pcorr = 2.11 × 10-3) as genes to prioritise. Inclusion of both autoantibody targets and IIM-associated SNPs resulted in greater and more significant indirect and common interactor connectivity than either network alone, suggesting indirect interaction between autoantibody targets and proteins encoded by IIM-associated loci (Fig. 2b and c). Including both SNPs and MSA/MAA targets additionally identified PSMD3, HSPA1A and HSPA1B alongside UBE3B and TRAF6 and seven myositis-specific autoantibody (MSA) targets as genes to prioritise (Table 2). From visual assessment of the direct networks tumour necrosis factor receptor-associated factor 6 (TRAF6) appears to act as a hub protein, interacting with many proteins for both SNP and SNP-autoantibody target networks (Fig. 2a and b).
Text-based pathway analysis of autoantibody targets and associated SNP loci identifies significant inter-relatedness
GRAIL analysis identified 24/28 autoantibody targets as significant based on their connections to the other autoantibody targets in the literature (p corr < 0.01 for each autoantibody target). Only PRKDC, ABTB1, HMGCR and NT5C1A were not significant (Additional file 4). The twenty keywords identified by GRAIL confirmed previous understanding of MSA target functions (‘tRNA’, ‘synthetase’) and suggested autoantibody target involvement in post-translational modification (‘sumoylation’, ‘ubiquitin’) (Additional file 5). The VIZ-GRAIL (Visualising GRAIL connections) plot showed strong connections (p < 0.05), especially between the anti-synthetase targets and within known complexes such as CHD3-CHD4 and SAE1-UBA2 (Additional file 6).
In analysis of SNPs, only 4/19 regions were identified as significant based on their connections to the other associated regions (p corr < 0.01 for each candidate gene); the most likely candidate genes STAT4, CD27, CCL17 and CD28 (Additional file 4). The majority of keywords identified were immune-related, but also included ‘ubiquitin’ (Additional file 5). The VIZ-GRAIL plot showed many strong pairwise connections between candidate genes from associated loci such as C6orf21-GRB7, BLK-STAT4, CD27-RAG2 and CD28-RAG1 (Fig. 3).
GRAIL analysis of SNPs and autoantibody targets combined identified an increased number of autoantibody targets (25/27) and candidate genes (12/19) as significantly interconnected to the other associated SNP regions or autoantibody targets (p corr < 0.01 for each autoantibody target or candidate gene), suggesting inter-relatedness between autoantibody targets and gene products of IIM-associated loci (Additional files 4 and 7). Notably, STAT1 was identified as the most likely candidate gene in the combined analysis from the rs4853540 association, rather than STAT4. Keywords included ‘ubiquitin’ and ‘sumoylation’ and immune and tRNA synthetase-related terms (Additional file 6).
Gene ontology analyses confirm known pathways and suggest novel pathways involved in IIM pathogenesis
DAVID gene ontology analysis of 28 autoantibody targets confirmed previous knowledge of MSA target involvement in translation and RNA metabolism, including amino acid activation (p corr = 7.6 × 10-12), tRNA aminoacylation (p corr = 7.6 × 10-12), tRNA metabolic process (p corr = 6.5 × 10-9), non-coding RNA metabolic process (p corr = 2.1 × 10-9) and translation (p corr = 8.5 × 10-7) (Additional files 8 and 9). Autoantibody targets were shown also to be significantly involved in somatic recombination of immunoglobulin gene segments (p corr = 4.3 × 10-2). The only significant KEGG pathway was aminoacyl-tRNA biosynthesis (p corr = 9.8 × 10-12), which is unsurprising as nine of the autoantibody targets are aminoacyl tRNA synthetases. Gene ontology analysis of the associated SNP loci identified only immune response as significant (p corr = 0.015). Combined analyses of autoantibody targets and associated SNP loci identified additional significant immune-related gene ontologies such as somatic diversification of immune receptors (p corr = 2.9 × 10-3).
We report the first systematic protein-protein interaction and pathway analysis of myositis autoantibody targets and proteins encoded by IIM-associated loci . Our analyses identify significant direct and indirect connectivity and inter-related genes/proteins. Several genes are highlighted by PPI analysis including TRAF6, HSPA1A/B, UBE3B and PSMD3. Pathway analysis confirmed the over-representation of autoantibody target involvement in translation and suggests that ubiquitination may play an important role in the IIM disease processes.
The significant direct and indirect connectivity identified between autoantibody targets indicates physical interactions between these antigens and with common interactors and suggests related biological pathways. The significant indirect connectivity between associated SNP loci suggests that although proteins encoded by these loci do not directly interact, they may influence common pathways, in keeping with potential long-range regulatory effects. The more significant indirect and common interactor connectivity observed in the combined analysis suggests autoantibody targets and proteins encoded by IIM-associated loci act in the same pathways.
The direct networks produced by DAPPLE suggest an important role for TRAF6 interactions with both autoantibody targets and SNP-associated gene products in IIM. TRAF6 is an E3 ligase which acts as an intermediate signalling adaptor between signals such as pathogen-associated molecular patterns (PAMPs) and downstream pathways including nuclear factor-kB (NF-kb), c-Jun N-terminal kinase (JNK), mitogen-activated protein kinase (MAPK) and AMP-activated protein kinase (AMPK). TRAF6 is involved in immune responses via the regulation of inflammation , differentiation and proliferation of B cells  and development of medullary thymic epithelial cells . TRAF6 also is implicated in muscle atrophy . These roles support the hypothesis that TRAF6 may have an important role within muscle as part of the autoimmune, inflammatory process in IIM, and TRAF6 with associated pathways may represent potential therapeutic targets.
Of the seven autoantibody target-derived ‘genes to prioritise’ in the DAPPLE analysis of SNPs and autoantibody targets, five were aminoacyl tRNA synthetases (ARS). In the direct network four of these (TARS, HARS, GARS and YARS) directly interacted with TRAF6. Some of the statistical significance of these proteins may be derived from interactions within the ARS group. However, a proportion will be due to their interaction with TRAF6. Further study of these interactions may be valuable in understanding the link between presence of anti-ARS and the specific phenotype observed in anti-synthetase syndrome.
A further ‘gene to prioritise’, UBE3B, is not present in the direct networks, and derives significant connectivity indirectly via common interactor proteins. UBE3B encodes ubiquitin protein ligase E3B, which is involved in target recognition in the ubiquitination pathway. Another ‘gene to prioritise’, PSMD3, encodes a 26 s proteasome regulatory subunit involved in degradation of ubiquitinated proteins via the ubiquitin proteasome pathway (UPP). In the Immunochip study two associations were reported as intergenic of PRR5L/TRAF6 and UBE3B/MMAB . Our analyses could suggest TRAF6 and UBE3B as the causal genes in their respective regions. UBE2L3, an E2 ubiquitin-conjugating enzyme, was also highlighted in the Immunochip study as it has previously been implicated in other autoimmune diseases. Four of five prioritised SNP-derived genes identified in the DAPPLE analyses encode proteins related to the UPP (Fig. 4). Furthermore, ‘ubiquitin’ was the only keyword identified in all three GRAIL analyses (Additional file 6). In IIM, ubiquitination is involved in activation of NF-kB which in turn leads to upregulation of MHC class I, promotion of inflammatory cytokines leading to muscle fibre damage, and MyoD inhibition resulting in reduced myoblast differentiation (Fig. 4) .
The UPP is thought to be the most significant degradation system in myogenesis . Inhibition or knockdown of the proteasome prevents myoblast fusion and inhibits differentiation, whereas proteasomal activity increases during differentiation of mouse myoblast C212 cells  and has been found to be increased in IIM muscle tissue .
In inclusion body myositis (IBM), abnormalities of muscle protein homeostasis may lead to accumulation of proteins in muscle fibres, possibly as a result of impaired proteasome function. Gene mutations resulting in disrupted ubiquitination have been identified in several inherited myopathies such as nemaline and polyglucosan body myopathy [23, 24].
As illustrated by Fig. 4, enhanced proteasome activity generates MHC class I restricted antigens, leading to activation of the immune response. Stress-inducible heat shock 70 kDa (Hsp70) proteins, such as HSPA1A and HSPA1B, upregulated when cells are under stress, act as chaperones for transferring these antigens from the proteasome to be presented by MHC class I or II . Hsp70 proteins play a critical role in protecting against muscle damage in the innate immune response, promoting muscle regeneration and recovery, and maintaining skeletal muscle mass and integrity . This response decreases with increased age. As IIM is generally a disease of middle and old age it is thought that restoration of the heat-shock response may be a good therapeutic approach, especially in IBM where protein mishandling is most clearly seen in the form of protein aggregates in muscle fibres. Arimoclomol is a pharmacological agent that induces an increased heat-shock response in stressed cells but does not appear to affect unstressed cells, thus avoiding off-target responses . Arimoclomol has shown promise in ameliorating disease in both in vitro rat myoblast and in vivo mouse models of IBM and has passed safety testing in IBM patients . A further pharmacological agent, 17-N-allylamino-17-demethoxygeldanamycin (17AAG), upregulates Hsp70. Treatment with 17AAG has been shown to preserve mitochondrial function and prevent atrophy of C2C12 mouse myotubes in a tunicamycin-induced endoplasmic reticulum (ER)-stress model of IIM . This PPI analysis highlighting heat-shock proteins in association with IIM lends further support to targeting the heat shock response in development of future IIM treatments.
A limitation of this study is that MSA/MAAs target intracellular antigens and it is unclear whether the autoantibodies or the pathways in which the antigens are involved influence the pathogenesis of IIM, and how these antigens come to be exposed to the immune system in IIM break of tolerance. The IIM-associated SNPs are from an Immunochip study, therefore, the identification of immune-related keywords is likely influenced by the biased content of the Immunochip. However, these results are still informative and highlight non-immune pathways that could be investigated as potential therapeutic targets in IIM.
In conclusion, through PPI and pathway analyses, this study highlights several genes including TRAF6, HSPA1A/B, UBE3B and PSMD3, and pathways including ubiquitination, which may play important roles in IIM disease processes.
ARS, aminoacyl tRNA synthetase; BIND, Biomolecular Interaction Network Database; CTD, connective tissue disease; DAPPLE, Disease Association Protein-Protein Link Evaluator; DAVID, Database for Annotation Visualisation and Integrated Discovery; ER, endoplasmic reticulum; GRAIL, Gene Relationships Across Implicated Loci; GWAS, genome-wide association studies; IBM, inclusion body myositis; IIM, idiopathic inflammatory myopathies; KEGG, Kyoto Encyclopaedia of Genes and Genomes; LD, linkage disequilibrium; MAA, myositis-associated autoantibody; MHC, major histocompatibility complex; MINT, Molecular Interaction Database; MSA, myositis-specific autoantibody; PAMPs, Pathogen associated molecular patterns; PPI, protein-protein interactions; SNP, single nucleotide polymorphism; TRAF, tumour necrosis factor receptor-associated factor; UPP, ubiquitin proteasome pathway.
Betteridge ZE, Gunawardena H, Mchugh NJ. Novel autoantibodies and clinical phenotypes in adult and juvenile myositis. Arthritis Res Ther. 2011;13:1–7.
Rosen A, Casciola-rosen L. Autoantigens as partners in initiation and propagation of autoimmune rheumatic diseases. Annu Rev Immunol. 2016;15:1–15.26. doi:10.1146/annurev-immunol-032414-112205.
Casciola-rosen L, Nagaraju K, Plotz P, et al. Enhanced autoantigen expression in regenerating muscle cells in idiopathic inflammatory myopathy. J Exp Med. 2005;201:591–601. doi:10.1084/jem.20041367.
Levine SM, Rosen A, Casciola-rosen LA. Anti-aminoacyl tRNA synthetase immune responses: insights into the pathogenesis of the idiopathic inflammatory myopathies. Curr Opin Rheumatol. 2003;15:708–13.
Chuang H, Lee E, Liu Y, et al. Network-based classification of breast cancer metastasis. Mol Syst Biol. 2007;3:1–10. doi:10.1038/msb4100180.
Luo X, Huang L, Jia P, et al. Protein-protein interaction and pathway analyses of top schizophrenia genes reveal schizophrenia susceptibility genes converge on common molecular networks and enrichment of nucleosome (chromatin) assembly genes in schizophrenia susceptibility loci. Schizophr Bull. 2014;40:39–49. doi:10.1093/schbul/sbt066.
Ragnedda G, Disanto G, Giovannoni G, et al. Protein-protein interaction analysis highlights additional loci of interest for multiple sclerosis. PLoS One. 2012;7:6–12. doi:10.1371/journal.pone.0046730.
Tuller T, Atar S, Ruppin E, et al. Global map of physical interactions among differentially expressed genes in multiple sclerosis relapses and remissions. Hum Mol Genet. 2011;20:3606–19. doi:10.1093/hmg/ddr281.
Rothwell S, Cooper RG, Lundberg IE, et al. Dense genotyping of immune-related loci in idiopathic inflammatory myopathies confirms HLA alleles as the strongest genetic risk factor and suggests different genetic background for major clinical subgroups. Ann Rheum Dis. 2015:1–9. doi:10.1136/annrheumdis-2015-208119.
Rossin EJ, Lage K, Raychaudhuri S, et al. Proteins encoded in genomic regions associated with immune-mediated disease physically interact and suggest underlying biology. PLoS Genet. 2011;7:e1001273. doi:10.1371/journal.pgen.1001273.
Raychaudhuri S, Plenge RM, Rossin EJ, et al. Identifying relationships among genomic disease regions: predicting genes at pathogenic SNP associations and rare deletions. PLoS Genet. 2009;5:e1000534. doi:10.1371/journal.pgen.1000534.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57. doi:10.1038/nprot.2008.211.
Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13. doi:10.1093/nar/gkn923.
Miao-xin L, Yeung JMY, Cherny SS, et al. Evaluating the effective numbers of independent tests and significant p -value thresholds in commercial genotyping arrays and public imputation reference datasets. Hum Genet. 2012;131:747–56. doi:10.1007/s00439-011-1118-2.
Johnson AD, Handsaker RE, Pulit SL, et al. SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics. 2008;24:2938–9. doi:10.1093/bioinformatics/btn564.
Jiao S, Zhang Z, Li C, et al. The kinase MST4 limits inflammatory responses through direct phosphorylation of the adaptor TRAF6. Nat Immunol. 2015;16. doi:10.1038/ni.3097.
Iwata S, Yamaoka K, Niiro H, et al. Increased Syk phosphorylation leads to overexpression of TRAF6 in peripheral B cells of patients with systemic lupus erythematosus. Lupus. 2015;40:1–10.
Akiyama T, Maeda S, Yamane S, et al. Dependence of self-tolerance on TRAF6-directed development of thymic stroma. Science. 2005;308:248–52.
Kumar A, Bhatnagar S, Paul PK. TWEAK and TRAF6 regulate skeletal muscle atrophy. Curr Opin Clin Nutr Metab Care. 2012;15:233–9. doi:10.1097/MCO.0b013e328351c3fc.TWEAK.
Nagaraju K, Casciola-Rosen L, Lundberg I, et al. Activation of the endoplasmic reticulum stress response in autoimmune myositis: potential role in muscle fiber damage and dysfunction. Arthritis Rheum. 2005;52:1824–35. doi:10.1002/art.21103.
Cui Z, Hwang SM, Gomes AV. Identification of the immunoproteasome as a novel regulator of skeletal muscle differentiation. Mol Cell Biol. 2014;34:96–109. doi:10.1128/MCB.00622-13.
Ghannam K, Martinez-Gamboa L, Spengler L, et al. Upregulation of immunoproteasome subunits in myositis indicates active inflammation with involvement of antigen presenting cells, CD8 T-cells and IFNγ. PLoS One. 2014;9:1–12. doi:10.1371/journal.pone.0104048.
Gupta V, Ravenscroft G, Shaheen R, et al. Identification of KLHL41 mutations implicates BTB-Kelch-mediated ubiquitination as an alternate pathway to myofibrillar disruption in nemaline myopathy. Am J Hum Genet. 2013;93:1108–17. doi:10.1016/j.ajhg.2013.10.020.
Nilsson J, Schoser B, Laforet P, et al. Polyglucosan body myopathy caused by defective ubiquitin ligase RBCK1. Ann Neurol. 2013;74:914–9. doi:10.1002/ana.23963.
Stocki P, Morris NJ, Preisinger C, et al. Identification of potential HLA class I and class II epitope precursors associated with heat shock protein 70 (HSPA). Cell Stress Chaperones. 2010;15:729–41.
Senf SM, Howard TM, Ahn B, et al. Loss of the inducible Hsp70 delays the inflammatory response to skeletal muscle injury and severely impairs muscle regeneration. PLoS One. 2013;8:e62687. doi:10.1371/journal.pone.0062687.
Ahmed M, Machado PM, Miller A, et al. Targeting protein homeostasis in sporadic inclusion body myositis. Sci Transl Med. 2016;8:331-41.
Lightfoot AP, Jackson MJ, Mcardle A, et al. Endoplasmic reticulum (ER) stress-induced mitochondrial dysfunction and atrophy can be prevented by pharmacological upregulation of heat shock protein 70 (HSP) in cultured murine myotubes. Arthritis Rheumatol. 2014;66:S966.
Paul Martin (Centre for Musculoskeletal Research, University of Manchester, UK) assisted with the construction of VIZ-GRAIL plots. All genotyping was generated through the Myositis Genetics Consortium (MYOGEN), the members of which contributed the samples analysed in the Immunochip study. Study investigators of MYOGEN, in addition to the authors of this article, are as follows: Ingrid E. Lundberg (Rheumatology Unit, Department of Medicine, Karolinska Institutet, Stockholm, Sweden), Frederick W. Miller (Environmental Autoimmunity Group, Clinical Research Branch, National Institute of Environmental Health Science, National Institutes of Health, Bethesda, Maryland, USA), Peter K. Gregersen (The Robert S Boas Center for Genomics and Human Genetics, Feinstein Institute for Medical Research, Manhasset, NY, USA), Jiri Vencovsky (Institute of Rheumatology and Department of Rheumatology, 1st Faculty of Medicine, Charles University, Prague, Czech Republic), Katalin Danko (Division of Clinical Immunology, Department of Internal Medicine, University of Debrecen, Debrecen, Hungary), Vidya Limaye (Rheumatology Unit, Royal Adelaide Hospital and University of Adelaide, Adelaide, Australia), Albert Selva-O'Callaghan (Department of Internal Medicine, Vall d’Hebron General Hospital, Barcelona, Spain), Pedro M. Machado, Michael G. Hanna (MRC Centre for Neuromuscular Diseases, UCL Institute of Neurology, London, UK), Lauren M. Pachman (Northwestern University Feinberg School of Medicine and Ann & Robert H. Lurie Children’s Hospital of Chicago, Chicago, Illinois, USA), Ann M. Reed (Department of Pediatrics, Duke University, Durham, North Carolina, USA), Lisa G. Rider (Environmental Autoimmunity Group, Clinical Research Branch, National Institute of Environmental Health Science, National Institutes of Health, Bethesda, Maryland, USA), Hazel Platt (Centre for Integrated Genomic Medical Research, University of Manchester, UK), Øyvind Molberg (Department of Rheumatology, Oslo University Hospital, Oslo, Norway), Olivier Benveniste (Pitié-Salpêtrière Hospital, UPMC, APHP, Paris, France), Pernille Mathiesen (Paediatric Department, Naestved Hospital, Næstved, Denmark), Timothy Radstake (Department of Rheumatology and Clinical Immunology, University Medical Center Utrecht, Utrecht, The Netherlands), Andrea Doria (Department of Medicine, University of Padova, Padova, Italy), Jan De Bleecker, Boel De Paepe (Department of Neurology, Neuromuscular Reference Centre, Ghent University Hospital, Ghent, Belgium), Britta Maurer (Department of Rheumatology and Center of Experimental Rheumatology, University Hospital Zurich, Zurich, Switzerland), Leonid Padyukov (Rheumatology Unit, Department of Medicine, Karolinska University Hospital, Karolinska Institutet, Stockholm, Sweden), Terrance P. O'Hanlon (Environmental Autoimmunity Group, Clinical Research Branch, National Institute of Environmental Health Science, National Institutes of Health, Bethesda, Maryland, USA), Annette Lee (The Robert S Boas Center for Genomics and Human Genetics, Feinstein Institute for Medical Research, Manhasset, New York, USA), Christopher I. Amos (Geisel School of Medicine, Dartmouth College, Hanover, New Hampshire, USA), Christian Gieger (Helmholtz Zentrum München, Deutsches Forschungszentrum für Gesundheit und Umwelt (GmbH), Neuherberg, Germany), Thomas Meitinger (Institute of Human Genetics, Technische Universität München, Munich, Germany, Institute of Human Genetics, Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Germany), Juliane Winkelmann (Neurologische Klinik und Poliklinik, Klinikum rechts der Isar, Technische Universität München, Munich, Germany, Institute of Neurogenomics, Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Germany) Lucy R Wedderburn (Arthritis Research UK Centre for Adolescent Rheumatology, and Institute of Child Health, University College London, London, UK), Christopher Denton (Royal Free Hospital, London, UK), Herman Mann (Institute of Rheumatology, Prague), David Hilton-Jones (John Radcliffe Hospital, Oxford, UK), Patrick Kiely (St. George's Hospital, London, UK), Paul H. Plotz (National Institute of Arthritis and Musculoskeletal and Skin Diseases, National Institutes of Health, Bethesda, MD, USA), Mark Gourley (National Institute of Arthritis and Musculoskeletal and Skin Diseases, National Institutes of Health, Bethesda, MD, USA), Kelly Rouster-Stevens (Emory University School of Medicine, Atlanta, GA, USA), Adam M Huber (Dalhousie University, Halifax, Nova Scotia, Canada), Galina Marder (North Shore Univeristy Hospital, Great Neck, NY, USA) and Mazen Dimachkie (University of Kansas Medical Center, Kansas City, KS, USA).
This study was supported in part by: Myositis UK; Arthritis Research UK (18474); Association Francaise Contre Les Myopathies (AFM);The European Union Sixth Framework Programme (project AutoCure; LSH-018661); European Science Foundation (ESF) in the framework of the Research Networking Programme European Myositis Network (EUMYONET); The Swedish Research Council and the regional agreement on medical training and clinical research (ALF) between Stockholm County Council and Karolinska Institutet; the intramural research programme of the National Institute of Environmental Health Sciences (NIEHS); the National Institutes of Health (NIH); European Community’s FP6, AutoCure LSHB CT-2006-018661; The Cure JM Foundation; the European Science Foundation; the Wellcome Trust; the Henry Smith Charity UK; Action Medical UK; and the Swedish Research Council. The Czech cohort was supported by Project for Conceptual Development of Research Organization 00023728 from the Ministry of Health in the Czech Republic. JP is supported by a University of Manchester alumni “Research Impact” PhD studentship and President’s Doctoral Scholar Award. This report includes independent research supported by the National Institute for Health Research. The views expressed in this publication are those of the author(s) and not necessarily those of the National Health Service (NHS), the National Institute for Health Research or the Department of Health. The funding bodies had no role in the design of the study or collection, analysis, or interpretation of data or in writing the manuscript.
Availability of data and materials
All of the data used in this study are available in published papers.
All authors meet the criteria for authorship. JP contributed to the design of the study, performed the analyses, interpreted the data and wrote the manuscript. SR provided the Immunochip data, and assisted with data interpretation and critical revision of the manuscript. NM and ZB contributed to study design by advising on autoantibodies to include and critically revised the manuscript. PD, RC and WO contributed to the conception and design of the study and revision of the manuscript. HC and JL designed the study, assisted with acquisition of data, supervised the study and critically revised the manuscript. All authors read and approved the manuscript.
None to add.
The authors declare that they have no competing interests.
Consent for publication
Ethical approval and consent to participate
Input for analyses. Excel file of the information in Table 1 about myositis-specific and associated autoantibody targets and significant or suggestive SNPs used in analyses (XLSX 12 kb)
Protein-protein interaction analysis showing direct networks among myositis autoantibody targets. a Myositis autoantibody targets. b Subset of myositis autoantibody targets (one selected from each complex). Colours indicate significance of node p value (TIF 97 kb)
DAPPLE output. Excel file of the output provided by http://www.broadinstitute.org/mpg/dapple/dappleTMP.php from each DAPPLE analysis on SNPs, autoantibody targets, a subset of autoantibody targets and SNPs and autoantibody target subset combined (XLSX 6162 kb)
Most likely candidate genes in each significantly interconnected region identified by text-based pathway analysis, using GRAIL. Inputs were 19 myositis-associated SNPs, 28 autoantibody targets and 19 SNPs and 27 autoantibody targets combined. Due to the same candidate genes being selected from multiple inputs, CHD4, rs2286896, rs11724804 and rs3094013 were removed from the analysis. NS not significant (p > 0.01), ND not determined (DOCX 23 kb)
Keywords strongly linking the significant genes in each region identified by text-based pathway analysis, using GRAIL. Inputs were 19 myositis-associated SNPs, 28 autoantibody targets and 19 SNPs and 27 auotantibody targets combined. Due to the same candidate genes being selected from multiple inputs, CHD4, rs2286896, rs11724804 and rs3094013 were removed from the analysis (DOCX 15.1 KB)
VIZ-GRAIL plot of 28 myositis autoantibody targets. Outer circle labelled with regions, inner circle with genes. Thickness and brightness of lines represents the strength of the connections between genes (TIF 672 kb)
VIZ-GRAIL plot of 27 myositis autoantibody targets and 19 SNPs. Outer circle labelled with regions/SNPs, inner circle with genes. Thickness and brightness of lines represents the strength of the connections between genes (TIF 958 kb)
Significant DAVID results for gene ontology and pathway analysis of autoantibody targets and genes from associated loci. Genes were nominated by DAPPLE for 20 associated SNPs (the two MHC SNPs were excluded). Only significant results are shown (Benjamini-Hochberg false discovery rate <0.05) (DOCX 16 kb)