==== Front Orphanet J Rare Dis Orphanet J Rare Dis Orphanet Journal of Rare Diseases 1750-1172 BioMed Central London 1629 10.1186/s13023-020-01629-9 Research Gene expression profiling of orbital muscles in treatment-resistant ophthalmoplegic myasthenia gravis Europa Tarin A. taradr008@gmail.com 12 Nel Melissa nlxmel004@myuct.ac.za 12 http://orcid.org/0000-0002-8999-718XHeckmann Jeannine M. jeanine.heckmann@uct.ac.zajmheckmann01@gmail.com 12 1 grid.7836.a0000 0004 1937 1151Neurology Research Group, Division of Neurology, Department of Medicine, Faculty of Health Sciences, E8-74, New Groote Schuur Hospital, University of Cape Town, Cape Town, 7925 South Africa 2 UCT Neuroscience Institute, Cape Town, South Africa 11 12 2020 11 12 2020 2020 15 34625 8 2020 27 11 2020 © The Author(s) 2020Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.Background Unbiased in silico approaches applied to genome-wide data prioritized putative functional gene variants associating with treatment-resistant ophthalmoplegic myasthenia gravis (OP-MG). Although altered expression of genes harbouring these variants, or associated pathways, were shown in patient-derived transdifferentiated-myocyte models, gene expression in orbital-derived muscle was required to test the validity of the predictions. Methods We sampled orbicularis oculi muscle (OOM) and one paralysed extraocular muscle (EOM) from six individuals with OP-MG during blepharoptosis and re-alignment surgeries, respectively. For controls, the OOMs were sampled from four individuals without myasthenia undergoing surgery for non-muscle causes of ptosis, and one non-paralysed EOM. Using a qPCR array, expression of 120 genes was compared between OP-MG and control OOMs, profiling putative “OP-MG” genes, genes in related biological pathways and genes reported to be dysregulated in MG cases or experimental MG models, and in EOMs of cases with strabismus. Normalization was performed with two stable reference genes. Differential gene expression was compared between OP-MG and control samples using the ΔΔCT method. Co-expression was analysed by pairwise correlation of gene transcripts to infer expression networks. Results Overall, transcript levels were similar in OOMs and EOMs (p = 0.72). In OOMs, significant downregulated expression of eight genes was observed in OP-MG cases compared with controls (> twofold; p ≤ 0.016), including TFAM, a mitochondrial transcription factor, and genes related to the following pathways: atrophy signalling; muscle regeneration and contraction; glycogen synthesis; and extracellular matrix remodelling. Several microRNAs, known to be highly expressed in EOMs, are predicted to regulate some of these genes. Co-expression analyses of gene-pairs suggested high interconnectedness of gene expression networks in OP-MG muscle, but not controls (r > 0.96, p < 0.01). Significant inverse directions of gene-pair correlations were noted in OP-MG versus controls OOM networks (r ≥ 0.92, p < 0.001) involving most OP-MG genes overlapping prominently with muscle atrophy/contractility and oxidative metabolism genes. Conclusions The gene expression in orbital muscles derived from OP-MG individuals compared with normal controls, support the pathogenic hypothesis previously generated from whole genome sequence analyses. Repression of gene transcripts in OP-MG orbital muscle implicate tissue-specific regulatory mechanisms, which may inform future biomarker discovery approaches. Keywords Orbicularis oculiMuscle atrophyOxidative metabolismOphthalmoplegiaExtraocular musclesContractilityMyasthenia gravisTreatment-resistantGene expressionhttp://dx.doi.org/10.13039/501100001321National Research Foundation113416109084NRF-DAADEuropa Tarin A. Heckmann Jeannine M. http://dx.doi.org/10.13039/501100007112University of Cape TownHealth Sciences and Neuroscience Institute fellowshipsNel Melissa issue-copyright-statement© The Author(s) 2020 ==== Body Introduction Previously we identified individuals with otherwise characteristic myasthenia gravis (MG) and who responded to immune therapies as expected in their non-ocular muscles, but who remained with treatment-resistant ophthalmoplegia with/without ptosis [1, 2]. These individuals are most frequently MG subjects with juvenile onset symptoms, African genetic ancestry (see methods) and who have generalised disease with detectable circulating acetylcholine receptor (AChR) antibodies [1, 3]. Clinically, the treatment-resistant weakness of the extraocular muscles (EOMs) may range from moderate to severe. When severe, in which case most of the EOMs are paralysed, we refer to the subphenotype of treatment-resistant ophthalmoplegic MG, or OP-MG. In an African setting we anticipate up to 20% of cases may develop treatment-resistant EOM weakness [4]. To understand the underlying genetic pathogenesis of the OP-MG subphenotype we previously used next generation sequencing to dissect the molecular genetic landscape using an extreme phenotype approach i.e. OP-MG cases vs control MG. Genes associated with OP-MG were identified using single variant and gene-based cumulative variant statistical association analyses [5, 6]. The pathogenetic pathways postulated to be altered in OP-MG, and supported by expression studies in patient-specific dermal-derived transdifferentiated myocyte models, related to muscle regeneration and atrophy signalling [7]. The aim of this work was to investigate the expression of these putative OP-MG susceptibility genes/pathways in tissues relevant to the subphenotype, i.e. patient-derived orbital muscles, and to verify the hypotheses generated by the unbiased genomic studies in the tissues of OP-MG cases vs controls without MG. Briefly, a genome-wide single nucleotide variant analysis followed by prioritization based on skeletal muscle gene expression (from GTExPortal) [8], identified putative modifier variants in FAM92A1 and PEF1, which were more frequent in OP-MG vs control MG genomes (p < 1 × 10–5) [6]. As an adjunct, we sought genes with a cumulative variant burden at gene-level associating with OP-MG compared to control MG genomes, and after filtering by relevance in muscle (GTExPortal), identified several genes which may be relevant in OP-MG pathogenesis [6]. For simplicity we will refer to genes from both single-variant and gene-based analyses as putative “OP-MG genes” (Fig. 1, OP-MGWGS genes). Although the aim of the limb muscle-based prioritization was to select genes with biological relevance in this tissue, a limitation of this approach was that gene expression in EOM, as an unique muscle group, was not specifically considered due to the limited availability of gene expression data.Fig. 1 Genes profiled in the array organised by category. Ophthalmoplegic myasthenia gravis (OP-MG). “OP-MG genes” refers to genes identified by unbiased next generation sequencing analyses in OP-MG and control-MG subjects, either by whole genome sequencing (WGS) or whole exome sequencing (WES). “OP-MG pathways” refers to the main biological pathways associated with the known physiological functions of the putative OP-MG genes. The “atrophy signalling” genes were selected from biological pathways in reference 7 using associated KEGG pathways as resources. Ubiquitin proteasome system (UPS); experimental autoimmune myasthenia gravis (EAMG); acetylcholine receptor (AChR); muscle-specific kinase (MuSK); extracellular matrix (ECM). “Strabismus pathways” refers to genes with altered expression in human strabismic extraocular muscles. References in figure: (1) [5]; (2) [6]; (3) [50]; (4) [51]; (5) [12]; (6) [52]; (7) [32]; (8) [53]; (9) [18]; (10) [21]; (11) [22]; (12) [7]; (13) [54] The EOMs constitute a unique muscle allotype [9–11]. However, we have used biopsies of orbicularis oculi muscle (OOM) as a surrogate orbital muscle to profile transcripts which may closely resemble the EOM transcriptome. We have done so for several reasons; the relative inaccessibility of EOMs; the ultra-rare indications to offer oculoplastic surgery in MG and therefore the possibility of EOM biopsy; and the paralysed EOMs in OP-MG would reflect the gene expression profile of this ‘burnt-out’ stage which may differ from the active evolving pathophysiological process. As craniofacial muscles, both EOMs and OOMs share similar embryological origins and many of the anatomical, physiological and molecular characteristics which distinguish EOMs from skeletal limb muscle [12]. Both EOMs and OOMs have fast firing rates [13], high mitochondrial content and oxidative capacities [14, 15] and type II fibre predominance [9, 12]. Therefore, in contrast to previous work using transdifferentiated cell models which originated from dermal fibroblasts, the use of OOMs as a surrogate for EOM will provide important validation of gene expression tissue-specificity of the craniofacial myotranscriptome in OP-MG vs normal controls. Methods Patient selection and tissue collection OP-MG cases and controls selected for ocular re-alignment or lid elevation procedures were invited to participate and provided biopsy specimens with informed consent. The diagnostic criteria of OP-MG cases have been previously described [1, 6, 7], and those included in this report all had circulating anti-acetylcholine receptor (AChR) antibodies (> 0.4 nmol/L) (Table 1). Importantly, for OP-MG patients to be eligible for ptosis surgeries (lid elevation), OOM contractility had to be sufficient to minimize the risk of post-operative corneal exposure. In contrast, the EOMs of the OP-MG case had been paralysed for several years by the time the medial rectus was sampled. Control samples were derived from cases with strabismus or ptosis from non-MG causes. All subjects had self-reported African-genetic ancestry i.e. either black African genetic ancestry, or mixed-African ancestry comprising predominantly KhoiSan ancestry, with lesser genetic contributions from black Africans, Europeans and Asians (see [6]). The study was approved by the institutional Human Research Ethics Committee (HREC 257/2012).Table 1 Clinical information of patient-derived orbicularis oculi muscle and medial recti samples Age at surgery (years) Orbital muscle Diagnosis MG duration (years) 42 OOM OP-MG 18 15 OOMa OP-MG 7 27 OOM OP-MG 10 34 OOM OP-MG 14 55 OOMa Levator dehiscenceb NA 51 OOM Levator traumab NA 56 OOM Canalicular obstructionb NA 27 EOM OP-MG 15 77 EOM Secondary esotropia-blind from glaucoma NA Myasthenia gravis (MG) duration was taken from onset of symptoms. Extraocular muscle (EOM) refers to medial rectus; ophthalmoplegic MG (OP-MG); orbicularis oculi muscle (OOM) aSamples obtained from both eyes at independent surgeries bNormal OOM biopsied Biopsies of the orbital part of the orbicularis oculi (OOM) were received from those undergoing lid elevation procedures and medial recti muscle samples (EOM) from ocular realignment surgery. Both OOM and EOM samples were dissected by the ophthalmic surgeon without cautery or clamping, immediately placed in RNAlater (ThermoFisher) and stored at -80 °C. RNA isolation The orbital muscle samples were homogenized in TRI Reagent (ThermoFisher) using a handheld Pellet Pestle (Fisherbrand). After chloroform precipitation, RNA was isolated from the aqueous phase using a spin-column method (Direct-zol Miniprep kit, Zymo) according to the protocol which included an in-column DNAse step. RNA was quantified using Nanodrop 1000 software (v3.5.2, Inqaba Biotechnical industries) and RNA integrity (RIN) was evaluated using the PicoChip Bioanalyzer (Agilent 2100 Bioanalyzer). The muscle samples weighed 3–26 mg yielding 5–49 ng/µL RNA. Spectrophometric ratios suggested some impurities despite optimised RNA extraction protocols using phenol/chloroform precipitation, but were in the expected range for RNA isolated from muscle [16] with A260/280 range 1.6–1.9 and A260/230 range 0.4–1.6. RIN values ranged between 5.2 and 7.6. Custom array plate A total of 120 target and 5 reference genes were profiled in a custom qPCR array (Qiagen) using proprietary assays. The target genes were categorized into OP-MG genes, OP-MG pathways, MG pathways and strabismus pathways based on the source of evidence for their relationship to OP-MG, MG or strabismus, although there is overlap between categories (Fig. 1). OP-MG genes broadly refer to those identified by WGS [6] and the preceding whole exome sequencing (WES) analysis [5], although the prime objective of this study was to validate the genes identified by the WGS analysis. OP-MG genes also served as “seed queries” which informed the broader enquiry of additional gene candidates in related pathways (OP-MG pathways). Genes previously shown to be dysregulated by MG in human limb muscles [17], in vitro models or experimental autoimmune myasthenia gravis (EAMG) [18] were also included. As reduced contractility (as occurs in OP-MG EOMs) [1, 19] potentially alters the regulation of biological pathways in ocular muscle, we also profiled genes previously shown to be dysregulated in the EOMs of cases with ocular misalignment/strabismus due to other non-myasthenic causes. We postulated that functional variants in muscle contractility pathways may render EOMs vulnerable in the setting of MG weakness. The reference genes were selected from genes that have previously been validated in human [9, 20–22] and rodent EOMs [18, 23] and showed stability in a patient derived muscle model [7]. Reverse transcription, qPCR and normalization For each muscle RNA sample 40 ng was reverse transcribed using the RT2 First Strand Kit (Qiagen) according to the manufacturer’s protocol and qPCR performed using the RT2 SYBR Green Mastermix (Qiagen) on the QuantStudio 12K Flex thermocycler. No inhibition of reverse transcription or PCR was evident and no genomic DNA contamination was detected. Three control EOM samples with the lowest RNA concentrations and RIN values were excluded from data analysis for quantitation cycle (Cq) values > 35 for at least 15% of genes in the array (see Additional file 1: Table S1). Reference genes for normalization were selected based on the geNorm [24] and Bestkeeper methods for OOMs [25] and chosen based on stability of expression levels for the two EOMs. Gene expression analysis and statistical methods Raw Cq values from the array were analysed using the ΔΔCt method [26]. Differences between OP-MG and controls were expressed as fold change or log2 fold change (Prism Graphpad 8). For OP-MG vs control comparisons in OOMs, unpaired Student’s t tests for normally distributed 2−∆Cq values and Mann–Whitney U tests for non-parametric data, were performed. Shapiro Wilk tests were performed to assess normality. The Benjamini–Hochberg correction was performed on the results of the t tests and Mann–Whitney U tests in the OOM comparisons. Due to the small sample, the false discovery rate (FDR) correction was set at 15%. Statistical significance of the difference in gene expression by phenotype could not be calculated for EOMs. Using OOM data, pairwise gene co-expression by phenotype was performed by Pearson’s correlations. Genes not expressed in > 1 sample (Cq > 35) were excluded from these analyses (n = 5). Correlation plots using gene pairs with significant correlations (p < 0.01) were generated using the Corrplot R package. Hierarchical clustering of the OP-MG correlation plot was performed using Hmisc (R package) and the controls correlation plot was reordered to match for visual comparison. Co-expression networks were generated using Cytoscape (v.3.71). Potentially functional modules were identified using the MCODE Cytoscape plug-in. To determine which gene pairs were significantly different between OP-MG and controls, differential correlation analysis was performed using the DiffCorr R package [27]. Gene pairs with a local false discovery rate > 0.6 were excluded. Results Six OP-MG muscle (1 EOM; 5 OOM) and five control muscle (1 EOM; 4 OOM) samples passed data quality analysis (Table 1). Differential gene expression in orbicularis oculi muscles To assess differences in relative gene expression between 5 OP-MG and 4 control OOMs, data normalization was performed using the average values for RPLP0 and ACTN2 as reference genes (Additional file 1: Table S2). Figure 2 shows the genes that showed > twofold difference in gene expression levels between OP-MG vs control OOMs (FDR = 15%; p < 0.01). All the significantly differentially expressed genes were downregulated in the OP-MG muscles. Interestingly, the three most differentially repressed gene transcripts in the OP-MG samples were those genes included in the profiling array because of their altered expression levels reported in post-mortem strabismic EOMs compared with controls [21]. MYH2, which encodes the myosin heavy chain 2a isoform, found in fast-twitch muscle fibres and abundant in both OOMs [13] and EOMs [28], showed fivefold downregulation compared to controls (p = 0.008). VCAN and ANK1, genes related to extracellular matrix remodelling pathways, as well as genes related to muscle atrophy signalling, CAPN3 and MAPK6 (p < 0.009) were downregulated in the OP-MG OOMs. VCAN and ANK1 encode components of the extracellular matrix and muscle, respectively, that regulate myoblast fusion/differentiation and attachment to the sarcoplasmic reticulum in regenerating muscle [29]. CAPN3 (calpain 3) is involved in preserving sarcomeric integrity [30]. MAPK6/ERK3 encodes mitogen-activated protein kinase-6, and acts at integration points for multiple processes including myogenesis [31]. TFAM, encoding a critical mitochondrial DNA transcription factor [32], was differentially downregulated (p = 0.005) in OP-MG.Fig. 2 Genes differentially expressed in orbicularis oculi muscles. The differentially expressed genes between the OP-MG and control orbicularis oculi muscles represented as a heatmap (log2 fold change) PPP1R2 which was prioritized as an OP-MG gene and showed significant repression in OP-MG orbital muscles (2.9-fold; p = 0.016), encodes a key enzyme in glycogen synthesis, phosphatase-1 regulatory inhibitor subunit R-2. The repression of PPP1R2 transcripts together with DYRK2 (2.7-fold; p = 0.003), which phosphorylates glycogen synthase, suggest dysregulation of muscle energy metabolism in OP-MG muscles, but not control muscle. Overall, we found differential repression of genes involved in the maintenance of the extracellular matrix and sarcomeric stability, as well as mitochondrial biogenesis and muscle metabolism in orbital muscles from OP-MG cases vs normal controls. Interestingly, genes in the latter categories were included in the array based on postulated pathway involvement informed by genomic studies, whereas the former were considered “strabismus-pathway” genes at the time and included as a controls for the clinical phenotype of altered EOMs contractility in OP-MG. Differential gene-co-expression in orbicularis oculi muscle by phenotype Pairwise correlations identified 232 strongly correlated gene pairs in OP-MG OOMs and 92 gene pairs in controls (r > 0.96, p < 0.01). OP-MG and control correlation plots were generated using this data and hierarchical clustering (k-median) identified potentially functional modules [33] visualised in the OP-MG plot that were not seen in the control plot (Fig. 3a vs b). Figure 3c shows the visual network of co-expressed gene pairs in the OOMs of OP-MG samples (from Fig. 3a) and the most interconnected genes (10–15 gene–gene interactions each) involved putative OP-MG genes from the WGS analysis (FAM92A1) and from the WES analysis (CANX, DDX17, FAM69A) (Fig. 3c; unlabelled OP-MG genes are from the WES analysis). DDX17 (Fig. 3c solid arrow) is a master transcriptional regulator involved in muscle gene expression [34] and appears central to a subnode. CANX encodes calnexin, a chaperone protein that facilitates the assembly of AChR subunits [35], FAM92A1 has recently been recognized for its importance in mitochondrial functioning [36] and both FAM69A and PEF1 have poorly characterized biological function. PEF1 encoding Peflin, which may be involved with endoplasmic reticulum-golgi transport and/or calcium binding [37], connected with a ‘module’ of genes (Table 2) which were profiled because they showed altered muscle expression levels in EAMG [18], and they are known to regulate mitochondrial biogenesis and oxidative metabolism. PARK2, encoding parkin, a major regulator of mitophagy [32], appeared to form a subnode of negatively correlated genes (Fig. 3c, dashed arrow), several of which are involved in muscle regeneration (e.g. DDX17) and mitochondrial biogenesis. PPP1R2, which was differentially repressed in OP-MG muscle (Fig. 2), did not correlate with other transcripts in OP-MG muscle yet showed significant correlations with four genes in control OOMs (data not shown).Fig. 3 Visual comparison of gene co-expression between OP-MG and control orbicularis oculi muscles. a Ophthalmoplegic myasthenia gravis (OP-MG) orbicularis oculi muscles (OOMs) correlation plot with gene pairs (r > 0.96; p < 0.01). b Control correlation matrix (r > 0.98; p < 0.01) ordered to match hierarchical clustering configuration of A for comparison. c Gene co-expression network of OP-MG orbicularis oculi muscles (r > 0.96; p < 0.01); the OP-MG genes, identified by whole genome sequencing analysis, are labelled (intermodular gene in bold). Degree of interconnectivity is shown by increased node size. Interconnecting lines are weighted by strength of correlation. Dashed lines indicate negative correlations. The solid arrow indicates DDX17, which is central to a subnode and the most highly prioritized gene from the whole exome sequencing analysis [5]. The dashed arrow indicates PARK2 which is central to an inhibitory subnode Table 2 Genes in putative co-expression network with PEF1 transcripts in OP-MG orbicularis oculi muscle Gene Function (biological pathway) ANGPTL4 Oxidative metabolism PPARG Mitochondrial biogenesis/metabolism FMO2 Oxidative metabolism BCL2 Mitochondrial biogenesis/metabolism TIMP4 Inhibitor of metalloproteinases (extracellular matrix remodelling) MMP16 Metalloproteinase (extracellular matrix remodelling) Next, a differential co-expression algorithm was used to statistically identify gene pairs differing significantly and inversely between the two phenotypes (Fisher’s Z tests, p < 0.002). An example is showed in Fig. 4a with significant positive correlations in a gene pair in one phenotype (OP-MG or controls) but inverse correlations in the other (p < 0.01) (see Additional file 1: Figure S1). Figure 4b displays a correlation network of the significant inversely correlated gene pairs for the OP-MG network (implying that the control network showed opposite correlations). The importance of this analysis is that it shows inverse gene–gene interconnectivity (correlations) in OP-MG orbital muscle vs controls, which included seven (of nine) genes prioritized as OP-MG genes in the previous unbiased genomic comparisons viz. PEF1, FAM92A1, AKT1, MYL12B, PPP1R12C, ZFP36L2 and SH3BGR [6]. Interestingly, the gene product of SH3BGR is involved with myosin heavy chain isoforms and sarcomere function and is highly expressed in EOMs.Fig. 4 Differential and inverse co-expression of gene–gene pairs in orbicularis oculi muscles. a Examples of gene–gene interactions which are significantly correlated in both ophthalmoplegic myasthenia gravis (OP-MG) muscle and control muscle, either positively or negatively, but inversely (i.e. opposite directions; raw data in supplementary Fig. 1a). b Visualization of differentially co-expressed gene–gene pairs in orbicularis oculi muscles in the OP-MG network (p < 0.01). Solid lines indicate positive associations in OP-MG and dashed lines, negative association in OP-MG. The genes in bold are OP-MG genes identified in the whole genome sequencing analysis [6]. c A Venn diagram showing the broad biological categories represented by the most differentially correlated gene–gene transcript levels which are inverse by muscle phenotype (p ≤ 0.0001; supplementary Figure A). d Potential messenger RNA (mRNA) and microRNAs (miRNAs) regulatory networks in extraocular muscle. Differentially downregulated genes in OP-MG orbicularis oculi muscle (Fig. 2) and differentially co-expressed genes (Figs. 3, 4) are listed in the columns to show potential mRNA-miRNA interaction with miRs known to be highly expressed in extraocular muscles (EOMs). OP-MGwes refers to genes identified by whole exome sequencing analysis [5] Figure 4c depicts the representative biological categories of strongly correlated gene pairs in the OP-MG network and inverse in controls. Of the most differentially correlated gene pairs which were inverse by phenotype (n = 13; p ≤ 0.001), > 50% involved a gene whose transcripts were substantially altered in EAMG [17] suggesting altered MG pathways in the muscle transcriptome of OP-MG vs controls. The broad category of ‘muscle atrophy’ signalling was the most representative functional category accounting for two-thirds of the most differential and inversely correlated gene pairs by phenotype with oxidative metabolism the next most frequent category (p ≤ 0.001) (genes may be implicated in 2 categories). Downregulated genes in OP-MG orbital muscle may be regulated by microRNAs In OOMs the significant differential downregulation of transcripts in OP-MG suggest post-transcriptional gene repression by microRNAs (miRs) binding to 3′ regulatory variants in OP-MG as a potential regulatory mechanism. Therefore, the miRTarBase database [38] was interrogated for interactions between miRs previously shown to be highly expressed in EOMs [39], and genes differentially downregulated ≥ 1.5-fold in OP-MG, compared with controls. MiR-499 and miR-206, the most highly expressed miRs in EOM (> sevenfold higher than in skeletal muscle), potentially interact with three of the genes which showed differential repression in OOMs (Fig. 4d). In addition, CANX and IL6R, which feature prominently in the differential gene–gene cross-correlations (Fig. 4b), and predicted to have EOM-miRs interactions, were prioritized previously to have 3′ regulatory region variants present in the OP-MG, but not controls [5, 7]. Extraocular muscle Although only two EOMs (medical recti) passed quality control we briefly mention the results due to scarcity of data on these tissues from live donors. For the EOMs, ACTB was used for normalization (SD = 0.06 between samples; Additional file 1: Table S2). Of relevance to the interpretation of the OOM results, the normalized expression levels of the 125 genes were similar between the two EOMs and nine OOMs (p = 0.72), although the EOM-specific isoform, MYH3, differed significantly (p < 0.0001; Additional file 1: Table S3, Figure S2). As the comparison of gene expression between only two samples required cautious interpretation a quantile–quantile plot of the transcripts of the two samples showed a similar distribution of the majority of the genes with only isolated values > tenfold difference (Fig. 5a, b). We briefly discuss WIF1, IL6 and MT1A transcript levels which were ≥ 25-fold higher in OP-MG compared to control.Fig. 5 Genes differentially expressed in the medial rectus of an ophthalmoplegic myasthenia gravis (OP-MG) patient. a Quantile–quantile plot showing the distribution of log2 fold change values between the OP-MG and control extraocular muscle (medial recti) and b outlying values with > tenfold differences in gene expression values are represented in the heatmap WIF1 encodes Wnt inhibitory factor-1, and as an inhibitor of Wnt signalling may affect muscle regeneration and endplate AChR clustering [40]. MT1A encodes metallothionein-1A and protects cells against oxidative stress [in [18]]. Increased expression of metallothionein 1 has been reported in atrophic muscle, possibly by preventing the upregulation of the insulin growth factor (IGF)-1 pathway [41], and in the EOMs of the EAMG rodent models [18]. IL6 encodes interleukin-6, a pivotal myokine in satellite cell/myoblast proliferation amongst others, but with persistently elevated levels it facilitates muscle atrophy and inhibits IGF-1 signalling [42]. Therefore, the higher transcript levels of IL6 and MTIA in the chronically paralysed EOM of an OP-MG vs a non-paralytic control, could be compatible with severely atrophic muscle and/or MG-induced effects [17]. Discussion Gene expression profiling in orbicularis oculi muscles showed the genes which were previously prioritized as OP-MG-associated variants, were differently regulated in the muscles of OP-MG cases compared with normal controls, and supported their involvement in OP-MG pathogenesis. The original gene discovery strategy used an ‘extreme phenotype’ in which all the cases had MG, but who differed only by their EOM resistance to treatment observed in the clinic. In the current study we aimed to validate expression differences in these putative OP-MG genes/pathways using orbital muscles from affected cases vs controls without MG. Differential gene repression and gene–gene correlations differed significantly between the OOMs from OP-MG cases and control muscle. Importantly, genes such as PEF1 and FAM92A1, which featured prominently in the dysregulated gene correlation networks by phenotype, were genes with unknown function when prioritized and included in the custom panel. FAM92A1, has recently been found to be critical for mitochondrial ultrastructural integrity [36]. Although PEF1 remains with poorly characterised function [43], the strong correlations with the transcripts of several genes related to oxidative metabolism and extracellular matrix remodelling, suggest “guilt-by-association”. These findings are highly suggestive of perturbed crosstalk in the OP-MG orbital muscles between genes representing the muscle atrophy/regeneration and mitochondrial/oxidative metabolism functional categories. In addition to assessing the putative OP-MG gene transcripts in orbital muscles, we were also interested in whether the transcriptional profile of genes known to represent MG/EAMG and strabismus-related pathways, may cross correlate with OP-MG gene transcripts. Here we included strabismus-related genes based on reports showing dysregulation in maligned EOMs (strabismus) in a non-MG context. However, recently RNA-sequencing of cultured human myoblasts revealed that AChR-antibodies significantly impacted genes related to extracellular matrix and actin/myosin cytoskeleton pathways [44], which overlap with the genes categorised in Fig. 1 as ‘strabismus pathways’, such as myosin genes. Although only selected myosin genes were included in the OOM profiling, they featured as dysregulated genes (Figs. 2, 4), and underscore the previous hypothesis based on clinical and histopathological reports in MG muscles, that loss of muscle contractility is critical in driving OP-MG pathogenesis [19, 45]. These results support the postulate that the EOMs of susceptible individuals, carrying gene variants in the myosin genes/pathways amongst others, are particularly vulnerable during the active MG attack to extended periods of poor contractility. Genes related to ‘muscle atrophy signalling’ pathways at cellular level including those regulating muscle fibre size by protein synthesis (IGF1-Akt/mTOR signalling pathway), degradation (ubiquitin-proteosome and autophagy-lysosome pathways) and the balance maintained by the TGFβ-BMP superfamily [46], were common to differential gene-pair cross-correlations summarised in Fig. 4c. Interestingly, MYH2 which showed significant repression in the OP-MG OOMs, also featured prominently in the muscle genes regulated by AChR-antibodies [44]. MYH2 is a substrate of MuRF1 and therefore directly responsive to atrophy signalling [47]. MuRF1, in turn, is moderated by the IGF1/PI3K/Akt signalling pathway which inhibits induction of muscle atrophy pathways [47]. The observed gene-pair correlations, in which transcript levels in genes from these pathways in OP-MG muscle were in direct opposite to those observed in controls, suggest perturbations between the muscle atrophy, cytoskeletal components and regeneration pathways in OP-MG OOMs. A concern for the translation of these results to the clinic is that under experimental conditions using dexamethasone, myosin heavy chain transcripts are downregulated resulting in steroid-linked muscle atrophy [47]. Steroid therapy, such as prednisone, is an important adjunct to immune therapies in MG, and in the early management of myasthenic ocular muscle paralysis, we have observed improved clinical outcomes linked with higher doses of prednisone [4]. However, it is not known whether higher doses of prednisone will benefit all individuals with ocular muscle involvement such as those carrying functional variants in high-risk pathways [6, 48, 49]. The differential gene downregulation of gene transcripts in OP-MG muscle compared to controls, implicate regulatory mechanisms mediated by miRs. Bioinformatic tools predict several miR-mRNA interactions by miRs which are highly expressed in EOMs and genes which showed differential repression in OP-MG muscle. Furthermore, OP-MG associated genes with putative functional 3′ regulatory variants, such as DDX17 [5, 7] (Figs. 3c, 4b), featured prominently in the differential and inverse cross-correlations by phenotype in the orbital muscles, suggesting its involvement in dysregulated crosstalk. DDX17, which was highly expressed in both EOMs and OOMs, is a master regulator in muscle splicing events and miRNA biogenesis [34]. The study of circulating miRs as possible prognostic biomarkers for treatment-responsiveness in EOMs, may prove useful for future clinical treatment trials. We acknowledge a number of limitations to this work. The sample sizes are limited due to opportunistic sampling of rare events. As the opportunity to evaluate the expression of putative OP-MG genes in the EOM of a well-characterized OP-MG case is a rarity, even the comparison of a single OP-MG medial rectus sample to a control was considered a valuable addition despite years of paralysis in the OP-MG EOM which would have impacted the gene-transcript snapshot. Although the study conclusions are based on the expression levels within the larger sample of OOMs and despite the similar origins of EOMs and craniofacial muscles (orbital orbicularis oculi), the EOM allotype differs from OOMs in some respects notably fatigue resistance, multi-innervated fibres (compared to singly innervated fibres in OOMs) and EOMs are constantly regenerating [12, 15]. Nevertheless, the orbital-derived orbicularis oculi muscles were not clinically weak and therefore a good comparator to the normal control muscle transcripts, and the gene transcript levels did not differ significantly between EOMs and OOMs. Conclusion The profiling of patient-derived orbital muscle supports the previous unbiased, genome-driven hypotheses of candidate genes and pathways involved in the pathogenesis of OP-MG. The dysregulated gene expression in the orbital muscles of treatment-resistant ophthalmoplegic MG cases implicate pathways related to muscle contractility and mitochondrial homeostasis and strongly suggests altered extraocular muscle-specific regulatory events. Supplementary Information Additional file 1. Supplementary Table S1 present the results of the quality of RNA extracted from the orbital muscles. Table S2 depicts the raw data informing the reference gene selection and Table S3 the genes which were highly expressed in extraocular muscle and orbicular oculi muscles. Supplementary Figure S1. Differentially co-expressed gene pairs between OP-MG and control orbicularis oculi muscle ranked by significance derived from Fisher’s Z test. Figure S2. Scatter plot of gene expression levels in two extraocular muscles and nine orbicularis oculi muscles. Abbreviations MGMyasthenia gravis OP-MGOphthalmoplegic myasthenia gravis EOMsExtraocular muscles OOMsOrbicularis oculi muscles EAMGExperimental autoimmune MG AChRAcetylcholine receptor RINRNA integrity number qPCRQuantitative PCR CqQuantitation cycle value SDStandard deviation CVCoefficient of variation ∆Delta FDRFalse discovery rate Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Supplementary Information The online version contains supplementary material available at 10.1186/s13023-020-01629-9. Acknowledgements The authors wish to thank Dr M Lenake and Prof A Murray from the University of Cape Town (UCT) Ophthalmology division for their assistance in obtaining the samples with such care, and to the patients who agreed to the muscle biopsies. Authors’ contributions TAE generated, analyzed and interpreted the data and wrote the first draft. MN and JMH supervised the project and were both major contributors in editing the manuscript. All authors read and approved the final manuscript. Funding The National Research Foundation (NRF) of South Africa provided funding (113416, 109084, JMH) and the NRF-DAAD scholarship (TAE). MN received postdoctoral fellowships from the University of Cape Town (UCT) Health Sciences Faculty and UCT-Neuroscience Institute awards. Availability of data and materials The datasets during and/or analysed during the current study available from the corresponding author on reasonable request. Ethics approval and consent to participate The study was approved by the University of Cape Town Human Research Ethics Committee (HREC 257/2012) and all participants provided informed written consent. Consent for publication All persons participating in this study agreed to their anonymous data being published and shared. Competing interests The authors declare that they have no competing interests. ==== Refs References 1. Heckmann JM Nel M A unique subphenotype of myasthenia gravis Ann N Y Acad Sci 2018 1412 14 20 10.1111/nyas.13471 28984362 2. Heckmann JM Owen EP Little F Myasthenia gravis in South Africans: racial differences in clinical manifestations Neuromuscul Disord 2007 17 929 934 10.1016/j.nmd.2007.07.002 17720497 3. Heckmann JM Hansen P Van Toorn R Lubbe E Janse van Rensburg E Wilmshurst JM The characteristics of juvenile myasthenia gravis among South Africans S Afr Med J 2012 102 532 536 10.7196/SAMJ.5390 22668957 4. Europa TA Nel M Heckmann JM Myasthenic ophthalmoparesis: time to resolution after initiating immune therapies Muscle Nerve 2018 58 542 549 10.1002/mus.26172 29790193 5. Nel M Jalali Sefid Dashti M Gamieldien J Heckmann JM Exome sequencing identifies targets in the treatment-resistant ophthalmoplegic subphenotype of myasthenia gravis Neuromuscul Disord 2017 27 816 825 10.1016/j.nmd.2017.06.009 28673556 6. Nel M Mulder N Europa TA Heckmann JM Using whole genome sequencing in an African subphenotype of myasthenia gravis to generate a pathogenetic hypothesis Front Genet 2019 10 136 10.3389/fgene.2019.00136 30881381 7. Nel M Prince S Heckmann JM Profiling of patient-specific myocytes identifies altered gene expression in the ophthalmoplegic subphenotype of myasthenia gravis Orphanet J Rare Dis 2019 14 24 10.1186/s13023-019-1003-y 30696470 8. Abecasis GR Auton A Brooks LD DePristo MA Durbin RM Handsaker RE Kang HM Marth GT McVean GA An integrated map of genetic variation from 1,092 human genomes Nature 2012 491 56 65 10.1038/nature11632 23128226 9. Fischer MD Budak MT Bakay M Gorospe JR Kjellgren D Pedrosa-Domellof F Hoffman EP Khurana TS Definition of the unique human extraocular muscle allotype by expression profiling Physiol Genom 2005 22 283 291 10.1152/physiolgenomics.00158.2004 10. Porter JD Khanna S Kaminski HJ Rao JS Merriam AP Richmonds CR Leahy P Li J Andrade FH Extraocular muscle is defined by a fundamentally distinct gene expression profile Proc Natl Acad Sci U S A 2001 98 12062 12067 10.1073/pnas.211257298 11572940 11. Porter JD Israel S Gong B Merriam AP Feuerman J Khanna S Kaminski HJ Distinctive morphological and gene/protein expression signatures during myogenesis in novel cell lines from extraocular and hindlimb muscle Physiol Genom 2006 24 264 275 10.1152/physiolgenomics.00234.2004 12. Verma M Fitzpatrick K McLoon LK Extraocular muscle repair and regeneration Curr Ophthalmol Rep 2017 5 207 215 10.1007/s40135-017-0141-4 29109908 13. Cheng NC Liao SL Wang IJ Lin IC Tang YB Fiber type and myosin heavy chain compositions of adult pretarsal orbicularis oculi muscle J Mol Histol 2007 38 177 182 10.1007/s10735-007-9085-6 17370042 14. Porter JD Baker RS Muscles of a different 'color': the unusual properties of the extraocular muscles may predispose or protect them in neurogenic and myogenic disease Neurology 1996 46 30 37 10.1212/WNL.46.1.30 8559415 15. Goodmurphy CW Ovalle WK Morphological study of two human facial muscles: orbicularis oculi and corrugator supercilii Clin Anat 1999 12 1 11 10.1002/(SICI)1098-2353(1999)12:1<1::AID-CA1>3.0.CO;2-J 9890723 16. Kuang J Yan X Genders AJ Granata C Bishop DJ An overview of technical considerations when using quantitative real-time PCR analysis of gene expression in human exercise research PLoS ONE 2018 13 30196438 10.1371/journal.pone.0196438 17. Maurer M Bougoin S Feferman T Frenkian M Bismuth J Mouly V IL-6 and Akt are involved in muscular pathogenesis in myasthenia gravis Acta Neuropathol Commun 2015 3 1 10.1186/s40478-014-0179-6 25627031 18. Kaminski HJ Himuro K Alshaikh J Gong B Cheng G Kusner LL Differential RNA expression profile of skeletal muscle induced by experimental autoimmune myasthenia gravis in rats Front Physiol 2016 7 524 10.3389/fphys.2016.00524 27891095 19. Rautenbach RM Pillay K Murray ADN Heckmann JM Extraocular muscle findings in myasthenia gravis associated treatment-resistant ophthalmoplegia J Neuroophthalmol 2017 37 414 417 10.1097/WNO.0000000000000534 28742638 20. Andrade FH McMullen CA Lactate is a metabolic substrate that sustains extraocular muscle function Pflugers Arch 2006 452 102 108 10.1007/s00424-005-0010-0 16328456 21. Altick AL Feng CY Schlauch K Johnson LA von Bartheld CS Differences in gene expression between strabismic and normal human extraocular muscles Investig Ophthalmol Vis Sci 2012 53 5168 5177 10.1167/iovs.12-9785 22786898 22. Agarwal AB Feng CY Altick AL Quilici DR Wen D Johnson LA von Bartheld CS Altered protein composition and gene expression in strabismic human extraocular muscles and tendons Investig Ophthalmol Vis Sci 2016 57 5576 5585 10.1167/iovs.16-20294 27768799 23. Fischer MD Gorospe JR Felder E Bogdanovich S Pedrosa-Domellof F Ahima RS Rubinstein NA Hoffman EP Khurana TS Expression profiling reveals metabolic and structural components of extraocular muscles Physiol Genom 2002 9 71 84 10.1152/physiolgenomics.00115.2001 24. Vandesompele J De Preter K Pattyn F Poppe B Van Roy N De Paepe A Speleman F Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes Genome Biol 2002 3 RESEARCH0034 10.1186/gb-2002-3-7-research0034 12184808 25. Pfaffl MW Tichopad A Prgomet C Neuvians TP Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper–Excel-based tool using pair-wise correlations Biotechnol Lett 2004 26 509 515 10.1023/B:BILE.0000019559.84305.47 15127793 26. Schmittgen TD Livak KJ Analyzing real-time PCR data by the comparative C(T) method Nat Protoc 2008 3 1101 1108 10.1038/nprot.2008.73 18546601 27. Fukushima A DiffCorr: an R package to analyze and visualize differential correlations in biological networks Gene 2013 518 209 214 10.1016/j.gene.2012.11.028 23246976 28. Park KA Lim J Sohn S Oh SY Myosin heavy chain isoform expression in human extraocular muscles: longitudinal variation and patterns of expression in global and orbital layers Muscle Nerve 2012 45 713 720 10.1002/mus.23240 22499099 29. Stupka N Kintakas C White JD Fraser FW Hanciu M Aramaki-Hattori N Versican processing by a disintegrin-like and metalloproteinase domain with thrombospondin-1 repeats proteinases-5 and -15 facilitates myoblast fusion J Biol Chem 2013 288 1907 1917 10.1074/jbc.M112.429647 23233679 30. Beckmann JS Spencer M Calpain 3, the "gatekeeper" of proper sarcomere assembly, turnover and maintenance Neuromuscul Disord 2008 18 913 921 10.1016/j.nmd.2008.08.005 18974005 31. Coulombe P Rodier G Pelletier S Pellerin J Meloche S Rapid turnover of extracellular signal-regulated kinase 3 by the ubiquitin-proteasome pathway defines a novel paradigm of mitogen-activated protein kinase regulation during cellular differentiation Mol Cell Biol 2003 23 4542 4558 10.1128/MCB.23.13.4542-4558.2003 12808096 32. Tryon LD Vainshtein A Memme JM Crilly MJ Hood DA Recent advances in mitochondrial turnover during chronic muscle disuse Integr Med Res 2014 3 161 171 10.1016/j.imr.2014.09.001 28664093 33. van Dam S Vosa U van der Graaf A Franke L de Magalhaes JP Gene co-expression analysis for functional classification and gene-disease predictions Brief Bioinform 2018 19 575 592 28077403 34. Dardenne E Polay Espinoza M Fattet L Germann S Lambert MP Neil H RNA helicases DDX5 and DDX17 dynamically orchestrate transcription, miRNA, and splicing programs in cell differentiation Cell Rep 2014 7 1900 1913 10.1016/j.celrep.2014.05.010 24910439 35. Chang W Gelman MS Prives JM Calnexin-dependent enhancement of nicotinic acetylcholine receptor assembly and surface expression J Biol Chem 1997 272 28925 28932 10.1074/jbc.272.46.28925 9360963 36. Wang L Yan Z Vihinen H Eriksson O Wang W Soliymani R FAM92A1 is a BAR domain protein required for mitochondrial ultrastructure and function J Cell Biol 2019 218 97 111 10.1083/jcb.201806191 30404948 37. McGourty CA Akopian D Walsh C Gorur A Werner A Schekman R Regulation of the CUL3 ubiquitin ligase by a calcium-dependent co-adaptor Cell 2016 167 525 538 10.1016/j.cell.2016.09.026 27716508 38. Chou CH Shrestha S Yang CD Chang NW Lin YL Liao KW miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions Nucleic Acids Res 2018 46 D296 302 10.1093/nar/gkx1067 29126174 39. Zeiger U Khurana TS Distinctive patterns of microRNA expression in extraocular muscles Physiol Genom 2010 41 289 296 10.1152/physiolgenomics.00169.2009 40. Cisternas P Henriquez JP Brandan E Inestrosa NC Wnt signaling in skeletal muscle dynamics: myogenesis, neuromuscular synapse and fibrosis Mol Neurobiol 2014 49 574 589 10.1007/s12035-013-8540-5 24014138 41. Summermatter S Bouzan A Pierrel E Melly S Stauffer D Gutzwiller S Blockade of metallothioneins 1 and 2 increases skeletal muscle mass and strength Mol Cell Biol 2017 37 e00305 e316 10.1128/MCB.00305-16 27956698 42. Munoz-Canoves P Scheele C Pedersen BK Serrano AL Interleukin-6 myokine signaling in skeletal muscle: a double-edged sword? Febs J 2013 280 4131 4148 10.1111/febs.12338 23663276 43. Consortium GT The genotype-tissue expression (GTEx) project Nat Genet 2013 45 580 585 10.1038/ng.2653 23715323 44. Hong Y Liang X Gilhus NE AChR antibodies show a complex interaction with human skeletal muscle cells in a transcriptomic study Sci Rep 2020 10 11230 10.1038/s41598-020-68185-x 32641696 45. Europa TA Nel M Heckmann JM A review of the histopathological findings in myasthenia gravis: clues to the pathogenesis of treatment-resistance in extraocular muscles Neuromuscul Disord 2019 29 381 387 10.1016/j.nmd.2019.03.009 31029532 46. Sartori R Gregorevic P Sandri M TGFbeta and BMP signaling in skeletal muscle: potential significance for muscle-related disease Trends Endocrinol Metab 2014 25 464 471 10.1016/j.tem.2014.06.002 25042839 47. Glass DJ PI3 kinase regulation of skeletal muscle hypertrophy and atrophy Curr Top Microbiol Immunol 2010 346 267 278 20593312 48. Auret J Abrahams A Prince S Heckmann JM The effects of prednisone and steroid-sparing agents on decay accelerating factor (CD55) expression: implications in myasthenia gravis Neuromuscul Disord 2014 24 499 508 10.1016/j.nmd.2014.02.010 24703255 49. Nel M Buys JM Rautenbach R Mowla S Prince S Heckmann JM The African-387 C>T TGFB1 variant is functional and associates with the ophthalmoplegic complication in juvenile myasthenia gravis J Hum Genet 2016 61 307 316 10.1038/jhg.2015.146 26632886 50. McLoon LK Wirtschafter J Activated satellite cells in extraocular muscles of normal adult monkeys and humans Investig Ophthalmol Vis Sci 2003 44 1927 1932 10.1167/iovs.02-0673 12714625 51. Anderson BC Christiansen SP McLoon LK Myogenic growth factors can decrease extraocular muscle force generation: a potential biological approach to the treatment of strabismus Investig Ophthalmol Vis Sci 2008 49 221 229 10.1167/iovs.07-0600 18172096 52. McLoon LK Harandi VM Brannstrom T Andersen PM Liu JX Wnt and extraocular muscle sparing in amyotrophic lateral sclerosis Investig Ophthalmol Vis Sci 2014 55 5482 5496 10.1167/iovs.14-14886 25125606 53. Zhou Y Kaminski HJ Gong B Cheng G Feuerman JM Kusner L RNA expression analysis of passive transfer myasthenia supports extraocular muscle as a unique immunological environment Invest Ophthalmol Vis Sci 2014 55 4348 4359 10.1167/iovs.14-14422 24917137 54. Sekulic-Jablanovic M Ullrich ND Goldblum D Palmowski-Wolfe A Zorzato F Treves S Functional characterization of orbicularis oculi and extraocular muscles J Gen Physiol 2016 147 395 406 10.1085/jgp.201511542 27069119