# Integrative Data Mining Highlights Candidate Genes for Monogenic Myopathies

Osorio Abath Neto¹,², Olivier Tassy³, Valérie Biancalana¹,⁴, Edmar Zanoteli², Olivier Pourquié³, Jocelyn Laporte¹*

1. Dept. of Translational Medicine and Neurogenetics, IGBMC, INSERM U964, CNRS UMR7104, University of Strasbourg, Collège de France, Illkirch, Strasbourg, France,
2. Departamento de Neurologia, Faculdade de Medicina de São Paulo (FMUSP), São Paulo, Brazil, 3. Dept. of Development & Stem Cells, IGBMC, INSERM U964, CNRS UMR7104, University of Strasbourg, Collège de France, Illkirch, Strasbourg, France,
4. Faculté de Médecine, Laboratoire de Diagnostic Génétique, Nouvel Hopital Civil, Strasbourg, France

<figure><img src="image_1.png" /></figure>

## Abstract

Inherited myopathies are a heterogeneous group of disabling disorders with still barely understood pathological mechanisms. Around 40% of afflicted patients remain without a molecular diagnosis after exclusion of known genes. The advent of high-throughput sequencing has opened avenues to the discovery of new implicated genes, but a working list of prioritized candidate genes is necessary to deal with the complexity of analyzing large-scale sequencing data. Here we used an integrative data mining strategy to analyze the genetic network linked to myopathies, derive specific signatures for inherited myopathy and related disorders, and identify and rank candidate genes for these groups. Training sets of genes were selected after literature review and used in Manteia, a public web-based data mining system, to extract disease group signatures in the form of enriched descriptor terms, which include functional annotation, human and mouse phenotypes, as well as biological pathways and protein interactions. These specific signatures were then used as an input to mine and rank candidate genes, followed by filtration against skeletal muscle expression and association with known diseases. Signatures and identified candidate genes highlight both potential common pathological mechanisms and allelic disease groups. Recent discoveries of gene associations to diseases, like *B3GALNT2*, *GMPPB* and *B3GNT1* to congenital muscular dystrophies, were prioritized in the ranked lists, suggesting a *posteriori* validation of our approach and predictions. We show an example of how the ranked lists can be used to help analyze high-throughput sequencing data to identify candidate genes, and highlight the best candidate genes matching genomic regions linked to myopathies without known causative genes. This strategy can be automatized to generate fresh candidate gene lists, which help cope with database annotation updates as new knowledge is incorporated.

**Citation:** Neto OA, Tassy O, Biancalana V, Zanoteli E, Pourquié O, et al. (2014) Integrative Data Mining Highlights Candidate Genes for Monogenic Myopathies. PLoS ONE 9(10): e110888. doi:10.1371/journal.pone.0110888

**Editor:** Francisco J. Esteban, University of Jaén, Spain

**Received:** June 15, 2014; **Accepted:** September 18, 2014; **Published:** October 29, 2014

**Copyright:** © 2014 Neto et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability:** The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was supported by Institut National de la Santé et de la Recherche Médicale (INSERM), Centre National de la Recherche Scientifique (CNRS), Strasbourg University, Agence Nationale de la Recherche (ANR-11-BSV1-026), Association Francaise contre les Myopathies (15352), Muscular Dystrophy Association (186985), and Myotubular Trust (http://www.myotubulartrust.com/). OAN was supported by a fellowship from CAPES Foundation, Ministry of Education of Brazil, process number 1286/51-2 (http://www.capes.gov.br/). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have read the journal's policy and have the following competing interests: OT and OP are the intellectual proprietors of the Manteia data mining tool. This does not alter the authors' adherence to PLOS ONE policies on sharing data and materials.

* Email: jocelyn@igbmc.fr

## Background

A large number of disorders affecting skeletal muscle have a genetic basis, with multiple modes of inheritance. They are classified based on phenotype and histopathological features into several groups, which include muscular dystrophies, congenital myopathies and myotonic syndromes, among others (Table 1) [1]. Muscular dystrophies and congenital muscular dystrophies, for example, are characterized by dystrophic changes on muscle biopsy, as opposed to congenital myopathies, which have non-dystrophic peculiar histopathologic findings [2–5]. Despite being rare, most inherited myopathies impose a heavy burden on the life of affected persons, and have a strong impact on the health care system. The identification of the causative gene and mutations is often a pre-requisite for genetic counseling and potentially

prenatal diagnosis, improved disease care, and access to more specific therapies or inclusion in clinical trials. A lot of advances have been made in the last few decades on the molecular bases of inherited myopathies, which included the discovery of about 130 genes associated with different disorders [1]. Still, it is estimated that around 40% of patients afflicted with myopathies remain without a molecular diagnosis, supporting the implication of additional genes [6,7]. Further identification of these genes is the focus of a tremendous research effort at present, and will help understand pathological mechanisms and defining novel drug targets.

Next-generation sequencing (NGS) is a relatively new technology that enables massive parallel sequencing of a huge number of bases. It has revolutionized molecular diagnosis and genetic
<table><caption>**Table 1.** Breakdown of disease groups and known associated genes.</caption><thead><tr><td>Disease group</td><td>Main diseases</td><td>Associated genes</td></tr></thead><tbody><tr><td>Muscular dystrophies</td><td>Duchenne and Becker muscular dystrophies, Emery-Dreifuss muscular dystrophy, Limb-girdle muscular dystrophies</td><td>ANO5, CAPN3, CAV3, DAG1, DES, DMD, DNAJB6, DPM3, DUX4, DYSF, EMD, FHL1, FKRP, FKTN, LMNA, MYOT, PABPN1, PLEC, POMGNT1, POMT1, POMT2, PTRF, SGCA, SGCB, SGCD, SGCG, SYNE1, SYNE2, TCAP, TMEM43, TNPO3, TRAPPC11, TRIM32, TTN</td></tr><tr><td>Congenital muscular dystrophies</td><td>Merosin-deficient CMD, Dystroglycanopathies, Ulrich and Bethlem myopathies</td><td>CHKB, COL6A1, COL6A2, COL6A3, DNM2, DPM2, FHL1, FKRP, FKTN, GTDC2, ISPD, ITGA7, LAMA2, LARGE, LMNA, POMGNT1, POMT1, POMT2, SEPN1, TCAP</td></tr><tr><td>Congenital myopathies</td><td>Centronuclear myopathy, Nemaline myopathy, Central core disease</td><td>ACTA1, BIN1, CCDC78, CFL2, CNTN1, DNM2, KBTBD13, KLHL40, MEGF10, MTM1, MTMR14, MYH2, MYH7, NEB, RYR1, SEPN1, STIM1, TNNT1, TPM2, TPM3, TRIM32, TTN</td></tr><tr><td>Metabolic myopathies</td><td>Glycogen storage diseases (Pompe, McArdle), Lipid storage diseases (CPTII deficiency)</td><td>ACADVL, AGL, CPT2, ENO3, GAA, GBE1, GYG1, GYS1, LDHA, LPIN1, PFKM, PGAM2, PGK1, PGM1, PHKA1, PNPLA2, PYGM, SLC22A5, SLC25A20</td></tr><tr><td>Congenital myastenic syndromes</td><td>Acetylcholine receptor deficiency, Choline acetyl transferase deficiency, Escobar syndrome</td><td>AGRN, CHAT, CHRNA1, CHRNB1, CHRND, CHRNE, CHRG, COLQ, DOK7, DPAGT1, GFPT1, LAMB2, MUSK, RAPSN, SCN4A</td></tr><tr><td>Myotonic syndromes</td><td>Myotonic dystrophy type 1 (Steinert disease), Schwartz-Jampel disease</td><td>ATP2A1, CAV3, CNBP, DMPK, HSPG2</td></tr><tr><td>Ion channel muscle diseases</td><td>Myotonia congenita, Hyperkalemic periodic paralysis, Paramyotonia congenita</td><td>CACNA1S, CLCN1, SCN4A</td></tr><tr><td>Vacuolar myopathies</td><td>Myopathy with excessive autophagia, Danon disease, Inclusion body myopathy with Paget disease of bone and frontotemporal dementia</td><td>EPG5, GNE, LAMP2, VCP, VMA21</td></tr><tr><td>Myofibrillar myopathies</td><td>Alpha-B-crystallin related myofibrillar myopathy, Desmin related myofibrillar myopathy</td><td>BAG3, CRYAB, DES, FLNC, LDB3, MYOT, SEPN1</td></tr></tbody></table>

doi:10.1371/journal.pone.0110888.t001

research, as it represents a cost-effective way of testing several genes at once in disorders with genetic heterogeneity, such as myopathies [8–10]. Moreover, exome sequencing (ES) or genome sequencing (GS) aid in the discovery of new genes associated with various diseases [11,12]. There has recently been a surge in publications that use NGS to discover new genes associated with diseases, including myopathies [13–20].

The biggest challenge of NGS is to cope with the complexity of analyzing the massive amount of variants generated by the approach. Indeed, comparing two unrelated individuals may lead to about 3 million matching variants in their genomes or about 20,000 in their exomes, but only one of these variants can cause a monogenic disease. The resolution of this issue demands good filtering pipelines to exclude common or meaningless variants, based on the biochemical function of genome location as studied through the ENCODE project [21], and on relationships between human variations and phenotype as in ClinVar and in locus specific databases [22,23]. In addition, ranking systems can help prioritize validation of the most promising variants. It makes sense to focus on genes presumably implicated in the disease process via functional, structural or phenotypical links with known genes. One of the approaches to collect and compare these data is via *in silico* analysis using a multitude of open-access knowledge information sources. This approach has been recently done successfully for some disorders but not yet for myopathies [24,25]. Lists of candidate genes thus generated can be ranked and used to prioritize variants resulting from NGS analysis.

Here, we propose ranked lists of candidate genes for individualized groups of inherited myopathies and related diseases that were obtained via data mining of online information databases.

These lists can be coupled to NGS analyses pipelines to help filter and prioritize variants aiming at the discovery of novel genes. We also put forward a number of genetic and functional insights taken from the generation of signatures for such disease groups to suggest common pathological pathways between them that can be subject of further scrutiny.

## Methods

### Classification of myopathy genes into 9 overlapping disease groups

The disease groups and associated known genes were based on a modified version of the Gene Table of Neuromuscular Disorders (GTNMD) [26]. We selected the following disease groups, which are primarily related to skeletal muscle pathology: Muscular Dystrophies, Congenital Muscular Dystrophies, Congenital Myopathies, Myotonic Syndromes, Ion Channel Muscle Diseases, Metabolic Myopathies, and Congenital Myasthenic Syndromes. To cope with an ill-defined classification of “Other Myopathies” in the GTNMD, we decided to cluster genes from this group into two new disease groups, Myofibrillar Myopathies and Vacuolar Myopathies. A literature search was performed to find recently published genes not yet listed in the Gene Table version that was used in our present study, which resulted in the addition of the following genes: *VMA21* [16] and *EPG5* [14] to the Vacuolar Myopathies group; *TRAPPC11* [13] and *TNPO3* [19,27] to the Muscular Dystrophies group; and *STIM1* [20], *CCDC78* [15] and *KLHL40* [17] to the Congenital Myopathies group.

The disease groups have some degree of overlap due to phenotypic heterogeneity of certain genes. For example, *SEPN1* is
implicated in multi-minicore disease (a congenital myopathy), and in rigid-spine muscular dystrophy (a congenital muscular dystrophy); CAV3 both causes limb-girdle muscular dystrophy 1C (a muscular dystrophy), and rippling muscle disease (a myotonic syndrome). The largest overlap is found between Muscular Dystrophies and Congenital Muscular Dystrophies, with 8 genes out of 34 muscular dystrophy-associated genes also found among 20 congenital muscular dystrophy-associated genes.

GTNMD's disease group "Distal Myopathies" was not included as a separate class in this work due to the lack of a gene uniquely associated with it - all genes were also found in other disease groups. Non-myopathy disease groups, such as ataxias, neuropathies, and motor neuron diseases, were also not included, as well as genes that, although listed in the GTNMD under included disease groups, do not lead to a skeletal muscle phenotype. This was the case of *MYBPC3*, implicated in cardiomyopathies, removed from the Congenital Myopathies group; *PRKAG2*, which causes a glycogen storage disease of the heart, not included in the Metabolic Myopathies group; and genes excluded from the Ion Channel Muscle Diseases group because they lead to various ataxia and cardiac arrhythmia syndromes, while not resulting in periodic muscle paralysis. The full list of genes and disease groups used in this work can be found in Table 1.

## Data-mining from online databases to address complex biological questions

We used the data mining system Manteia [28], a public resource available online (manteia.igbmc.fr) that retrieves and combines data from freely available online data sources such as Ensembl, Reactome, OMIM, NCBI, Human Phenotype Ontology (HPO), Gene Ontology (GO), Mouse Genome Informatics, and InterPro. Manteia makes it possible to address complex biological questions by running several queries at the same time to mine and statistically analyze gene sets to highlight their annotation specificities compared to the rest of the genome. This study was conducted with Manteia version 2 with data downloaded in June 2013 from the different databases used in the system.

Using Manteia's orthology module, we analyzed human gene sets and their mouse orthologs to find an enrichment on statistically significant terms within several annotation categories, including Gene Ontology (GO), Human Phenotype Ontology (HPO), Mammalian Phenotype Ontology (MPO), pathways (Reactome), protein motifs (Interpro) and interacting complexes (Reactome). Gene length was not taken into account as there is no clear enrichment of large genes mutated in myopathies; while some large genes are indeed implicated (*TTN*, *NEB*), smaller genes were found to accumulate mutations along their sequence (e.g. *ACTA1*). Manteia calculates the enrichment of each term in the gene set compared to all genes in the genome, and sorts the terms according to individual statistical significances.

The list of specific terms for each data set can then be used to screen the genome looking for genes that have similar properties. This is achieved using a query builder, which outputs a list of candidate genes ranked according to their similarity with the data set annotation signature and the weight given for each term.

## Extraction of specific signatures for each disease group based on known genes

Statistical analysis for human genes in each disease group was individually performed for Gene Ontology (GO) terms, Human Phenotype Ontology (HPO) terms, pathways, complexes and protein motifs. Mouse orthologs were additionally used to get

statistical breakdowns of Mammalian Phenotype Ontology (MPO) terms.

Signatures were represented by a weighted combination of GO terms, HPO terms, MPO terms, and what were collectively called "Interactions Annotation" (IA) terms - pathways and protein complexes descriptors from Reactome and descriptors of protein motifs from InterPro. For each disease group, terms were chosen from the various domains in order to obtain a signature of the disease group. We used the following criteria to select GO terms, HPO terms and MPO terms for each group: 1) significance *p*-value less than 0.05 (corrected using the Benjamini-Hochberg false discovery rate (FDR) procedure); 2) occurrence in the disease group gene set greater than 1; 3) occurrence in the genome <800; 4) GO level (or HPO level or MPO level) >2.

The FDR-BH correction of the *p*-value was chosen because it reduced the large size of the list of resulting terms while not being as stringent as the Bonferroni correction. Terms with only one occurrence in the gene set were deemed not representative of the set. Criteria 3 and 4 enrich for specificity and are closely associated owing to less specific terms (higher ontology level) being associated with a large number of genes; such general terms would not only be unproductive in compounding a signature for a disease group, but also could degrade the performance of a complex query.

For IA terms, the restriction on the occurrence of only one gene in the set was dropped with the aim of improving the scores of new genes related to protein function linked to single known myopathy gene. Indeed, a large proportion of significant terms have only one occurrence in any given disease group gene set. Finally, criteria 4 does not apply to Interpro and Reactome data, which are not structured in defined hierarchies as gene and phenotype ontologies.

## Ranking formula based on weighted scores of signature terms

After experimenting with different signature definitions, we decided to define a signature as having an equal contribution of GO terms, phenotype terms (HPO and MPO terms) and IA terms (Figure 1), so as not to *a priori* give more importance to any term set. A signature with a stronger component of phenotype terms, for example, yields a list of purported candidates strongly biased to genes implicated in known diseases or for which mouse models have been extensively phenotyped. Likewise, if GO terms are the main component of the signature, genes with functional links are preferentially ranked. Finally, IA terms boost the interactome of the known genes to the top of the ranked lists.

Thus, in our approach, for each disease group, the weight of individual terms was calculated so that the added weights of all GO terms was the same as the total PO score (added HPO terms weights combined with added MPO terms weights) and as the added weights of all IA terms, which was arbitrarily set as 1000.

In the GO domain, we defined strata of term weights with percentile cutoffs so that terms with a higher significance would respond for a larger share of the total GO score. The top 20% of terms (p80) contributed to 40% of the total GO score, a middle tier comprising 40% of terms (p40-p80) contributed to an additional 40% of the total GO score, and the lower 40% of terms (p40) provided 20% of the total GO score.

A similar approach was used to calculate individual weights for HPO and MPO terms, with the exception that, as the total PO score reflects the combination of equal shares of HPO and MPO terms, the maximum score of either HPO terms or MPO terms was set at 500.

The weight of each IA term, on the other hand, was the same no matter its position on the corresponding list, and was simply
<figure><img src="image_2.png" /><figcaption>**Figure 1. Integrated data mining workflow.** A signature of a disease group, composed of weighted terms, is generated from statistical analyses of genes already implicated in diseases of the group. Terms come from the three main annotation groups, GO (Gene Ontology), PO (Phenotype Ontology, an aggregate of Human Phenotype Ontology and Mammalian Phenotype Ontology) and IA (Interactions Annotation), are mined using Manteia and receive weights proportional to the their enrichment in the set of genes implicated in the disease group, as compared to the set of all genes in the human genome. Weights are attributed to terms so that annotation groups contribute equally to the composition of the signature. The signature of the disease group is then used to mine the genome for additional genes. Every gene in the genome receives a score equal to the sum of weights of terms that describe the gene if they match terms that define the disease group signature, for a maximum possible score of 3000. Further filtering steps mark genes that have low relative skeletal muscle expression or are annotated with known diseases.
doi:10.1371/journal.pone.0110888.g001</figcaption></figure>

calculated as 1000 divided by the number of significant IA terms for each disease group. This approach helped mine all genes that interacted with any single gene in the training set, provided pathways and interactions were statistically significant.

While the choice of percentile cutoffs that define the strata of weights was arbitrary, we observed that the modification of the cutoffs did not result in substantially different ranked candidate lists for each disease group, as long as the signature definition is the same. All terms for every domain in each disease group, with their corresponding calculated weights, can be found in Table S1.

## Generation of ranked lists of candidate genes for each disease group

Manteia's query builder feature was used to filter genes in the human genome that matched the signature defined for each disease group. Queries combining terms that constitute the signature were run to obtain a list of genes ranked by a score represented by the sum of all matched term weights (Figure 1). More specifically, a gene score results as the sum of a GO score (sum of weights of the disease group's signature GO terms that match the gene's GO terms), a PO score (sum of weights of matching HPO and MPO terms), and an IA score (sum of weights
of matching IA terms). The maximum total score a gene can receive is thus 3000 (1000 for GO score+1000 for PO score+1000 for IA score). To deal with MPO terms applied to murine orthologs, the predicted best human orthologs were selected using Manteia's ORTHO function after the ranking process. The ranked lists for each of the 9 disease groups, including annotation detailed in the next subsection, can be found in Table S2.

## Additional filtering of ranked candidate genes using expression data and association with human diseases

The ranked lists for each disease group include the known genes of the group, which were used to create the signatures, genes known to be associated with myopathies but implicated in other disease groups, and genes that are not listed in any disease group and thus represent potentially good candidate genes for myopathies. Among those, additional filtering was performed using tissue expression databases. Data from Illumina Body Map E-MAT513, established from mRNA-Seq of 16 human tissues, was downloaded for every gene, and genes with no expression in skeletal muscle or with an expression in skeletal muscle that was less than a cutoff of 10% of the maximum expression found in any other tissue were excluded. The rationale behind this filtering is that if a gene is expressed in a tissue other than skeletal muscle at a much higher level, one expects such a gene to be implicated in disorders primarily involving that tissue. The 10% cutoff was empirically determined due to the fact that all genes already implicated in myopathies have skeletal muscle expression levels above this cutoff. To deal with missing expression data and eventual heterogeneity in Illumina Body Map's expression database, genes ruled out by the 10% threshold and candidate genes within the 100 first positions in the rankings for each disease group were double checked using expression data from the Genotype-Tissue Expression Project (GTEx) [29].

For the Congenital Myasthenic Syndromes disease group, which includes diseases primarily related to neuromuscular junction protein defects, but also some peripheral nerve terminal protein defects, we decided to disregard the muscle expression filtering due to the fact that the implicated genes *AGRN*, *CHAT*, and *CHRNE* do not have significant skeletal muscle expression (they are instead expressed in the nerve terminal).

The lists of candidate genes after skeletal muscle expression filtration was further annotated with Online Mendelian Inheritance in Men (OMIM) data on existing human phenotype in the form of well-characterized diseases or syndromes, in order to easily identify genes biased by phenotype, such as *SMN1*, which results in a phenotype very similar to many myopathies, characterized by flaccid proximal limb weakness, but which gives rise instead to a motor neuron disease; or biased by interactions, as occurs to a number of carbohydrate metabolism genes that share common pathways to metabolic myopathies but cause instead inborn errors of metabolism without a muscle phenotype.

# Results

## Myopathy groups are clustered by gene ontology and protein function

To identify novel candidate genes for myopathies, we established an integrated data mining approach aiming first to extract specific signatures for disease groups encompassing previously implicated genes, and then to use these signatures to search for additional matching genes in the human genome. As detailed in Figure 1 and the methods section, this approach consists of a weighted ranking of three main sets of data: gene ontology, human

and mouse phenotypes ontologies, and “interactions annotation” incorporating pathways and protein motifs and complexes.

To test this approach and better visualize signature composition analysis, we first analyzed a training set that consisted in all myopathy-associated genes using the data mining system Manteia [28]. Figure 2 shows graphs with relationships between all known genes of the nine chosen disease groups. In particular, the combination of GO, PO and IA terms aggregate most genes that are part of the same myopathy group for metabolic myopathies, the congenital myasthenic syndromes, and the glycosylation components of congenital muscular dystrophies (Figure 2A). Of note, the gene *GFPT1*, which causes a congenital myasthenic syndrome with tubular aggregates, has mainly relationships with genes in the metabolic myopathy cluster, presumably because it codes for an enzyme in the metabolism of glycoproteins. Another large cluster in the graph encompasses the main genes implicated in muscular dystrophies and congenital or myofibrillar myopathies, without subdivision, suggesting a strong overlap in the function of the related genes and potentially in the pathogenesis. This approach can thus retrieve several phenotypic and pathologic clusters. However, applying only the human phenotype ontology analysis generates a single large, highly connected graph (Figure 2B), even when the threshold for representing an edge in the graph - number of matching HPO term between two genes - is increased or decreased, or when HPO term hierarchy is taken into account. This means that genes implicated in myopathies share a common hierarchy of phenotype ontology terms, e.g. with most genes annotated with muscle weakness or abnormal muscle physiology related terms. While they do not help separate genes into disease groups, HPO terms are important to help emerge genes with phenotype annotation associated to skeletal muscle. GO terms and IA terms, on the other hand, are responsible for the final clustering. Different myopathy groups appear using only GO terms (Figure 2C), while IA terms, even considering a lower threshold of 5 terms shared between genes, create smaller clusters of genes that interact closely by sharing the same pathways, interactions complexes or motifs, such as constituents of collagen VI, genes responsible for the assembly of nicotinic cholinergic receptors, or conglomerated proteins involved with the sarcomere (Figure 2D). Only the combination of the different GO, PO and IA terms reaches the most precise clustering.

## Characterization of disease groups via biological processes annotation

We next aimed to extract specific signatures for each disease group, classified based on the Gene Table of Neuromuscular Disorders [1]. Statistical analysis of known genes was conducted for each disease group.

GO terms include three types of ontologies: cellular components indicate the localization of gene products; molecular function refers to the normal roles of genes at the molecular level; and biological processes describe the higher-order roles of genes from a biological perspective. Four main general skeletal muscle-related biological processes were extrapolated from the hierarchy of GO terms: muscle contraction, calcium homeostasis, muscle development, and muscle intracellular organization (Table 2). Analysis of the breakdown of biological process-related GO terms that make up the signatures of different disease groups reveals differences in the implicated skeletal muscle processes and hints on other important biological processes that do not primarily involve the skeletal muscle.

Vacuolar myopathies and myofibrillar myopathies did not receive in their signature GO terms associated with biological processes, because the training set genes for these groups were
<figure><img src="image_3.png" /><figcaption>**Figure 2. Graph representation of relationships of known genes.** All known genes for the different disease groups were concurrently analyzed for matching terms in different ontologies. Nodes represent genes, and edges between two given nodes are depicted when the number of terms shared by the two connected genes is greater than a certain threshold. Edge width is proportional to the number of terms shared between two genes, and node size and color in a scale from green (lowest) to red (highest) is proportional to the number of associations of a gene in the graph. Closely related genes appear clustered together, and hubs in the graph appear centrally located. A: graph for combined terms from Gene Ontology (GO), Human Phenotype Ontology (HPO) and Interactions Annotation (IA), with a threshold of 30 matching terms. The cluster with a yellow background includes genes implicated in metabolic myopathies, the one with a red background groups congenital muscular dystrophy genes, and the cluster with a gray background represents genes associated with congenital myasthenic syndromes. B: graph for HPO terms with a threshold of 20 matching terms. C: graph for GO terms, with a threshold of 10 matching terms. Background colors correspond to clusters represented in A. D: IA terms with a threshold of 5 matching terms. The gray background highlights a cluster with gene that code subunits of cholinergic receptors, implicated in congenital myasthenic syndromes, the green one groups components of collagen VI, and the cluster with a blue background links elements of the contractile apparatus.</figcaption></figure>

doi:10.1371/journal.pone.0110888.g002

annotated with heterogeneous terms that did not attain statistical significance. Metabolic myopathies and congenital myasthenic syndromes inferred biological processes were, as expected, not primarily muscle related, but had mostly to do with glycogen

metabolism and neuromuscular junction, respectively. Congenital muscular dystrophies, while having two GO terms associated with muscle development, were also primarily annotated with non-muscle-specific biological processes, especially protein glycosyla-
<table><caption>**Table 2.** Composition of biological processes GO terms that make up the signature of each disease group.</caption><thead><tr><th>Disease Group</th><th>Muscle contraction</th><th>Calcium homeostasis</th><th>Muscle development</th><th>Muscle intracellular organization</th><th>Other terms</th><th>Total</th><th>Other term categories</th></tr></thead><tbody><tr><td>Congenital myopathy</td><td>7</td><td>10</td><td>12</td><td>3</td><td>13</td><td>45</td><td>Cardiac development, catabolism of nucleotides</td></tr><tr><td>Muscular dystrophy</td><td>9</td><td>1</td><td>14</td><td>23</td><td>30</td><td>77</td><td>Glycosylation, cardiac development, cardiac contraction</td></tr><tr><td>Congenital muscular dystrophy</td><td>0</td><td>0</td><td>2</td><td>0</td><td>20</td><td>22</td><td>Cardiac development, glycosylation</td></tr><tr><td>Metabolic myopathy</td><td>1</td><td>1</td><td>0</td><td>0</td><td>73</td><td>75</td><td>Glycogen metabolism</td></tr><tr><td>Congenital myasthenic syndrome</td><td>2</td><td>7</td><td>18</td><td>2</td><td>64</td><td>93</td><td>Neuromuscular junction, synapses</td></tr><tr><td>Myotonic myopathy</td><td>7</td><td>21</td><td>11</td><td>2</td><td>9</td><td>50</td><td>Heart contraction, circulation</td></tr><tr><td>Ion channel muscle disease</td><td>1</td><td>4</td><td>0</td><td>0</td><td>2</td><td>7</td><td>Ion transport</td></tr></tbody></table>

doi:10.1371/journal.pone.0110888.t002

tion. Myotonic myopathies and ion channel muscle diseases take the larger contribution from calcium homeostasis-related terms. Muscular dystrophies mainly involve muscle intracellular organization terms, but also receive some contribution from muscle development and muscle contraction terms. Other important biological processes for muscular dystrophies are associated with heart muscle contraction and development. Finally, for congenital myopathies, muscle development and calcium homeostasis seem to be the most significant processes, but muscle contraction-related terms also play a role, as well as processes not specific for skeletal muscle, such as catabolism of nucleotides - these appear enriched due to the association of *DNM2* to the catabolism of GTP, as well as *MYH7* and *TPM2* to the catabolism of ATP.

## Training set genes appear at the top of the ranked lists of the disease groups

We used the signature specific to each disease group to screen the whole set of human genes and identify candidate genes for myopathies. The breakdown of the gene score for these training set genes shows that a similar contribution of the different term domains can be consistently found throughout the various disease groups. GO score, PO score and IA score respond for approximately 30 to 45%, 40 to 55%, and 10 to 15% of the gene score, respectively. Table 3 shows the ranked lists of known genes for three disease groups (Congenital Myopathies, Muscular Dystrophies and Metabolic Myopathies), along with each gene score and breakdown of partial scores. Table S3 shows similar additional data for all disease groups.

As expected, genes already known to be mutated in the various disease groups, which were used as the training set to create the mining signatures, appear at the top of the ranked lists of data mining. Considering congenital myopathies, out of the 22 genes chosen as the training set, 19 genes appeared in the data mining, while genes *CCDC78*, *KBTBD13*, and *KLHL40* did not have annotation in the databases used at the time of this work. Thirteen of these genes were ranked within the first 100 genes, a coverage of 13/19 (68.4%). The muscular dystrophy group had 31 out of 34 genes of its training set appearing in the data mining list, and of these 23 were found within the first 100 ranked genes (79.3%). In the metabolic myopathy disease group, all 19 genes of the training set were ranked, and 18/19 were found within the top 100 genes in the rank (94.7%). Outliers among the known genes are mostly poorly annotated genes, and genes with no score are actually not annotated at all (see discussion). Thus, the high ranking of most previously implicated genes supports the signature choice having adequately defined the disease group.

## Proposed candidate genes after filtration

A number of candidate genes sharing disease group signatures with known myopathy genes are barely expressed in skeletal muscle or sometimes mutated in other diseases not affecting skeletal muscle. We thus added filtering steps based on tissue expression and known implication in diseases (see methods for details). Table 4 shows the top 8 ranked genes for each disease group after filtration on skeletal muscle expression and absence of link with diseases in OMIM (Online Mendelian Inheritance in Men, omim.org) database. Table S2 lists the full ranked list of genes for each disease group without filtration, but annotated with skeletal muscle expression and OMIM diseases, and can be linked to NGS filtering pipelines to help prioritization of novel gene discovery, as shown in the discussion. In the following paragraphs, we discuss a few genes found as candidates in some of the disease groups, to illustrate the connections between the integrated data mining results and evidence from the literature.
<table><caption>Table 3. Score breakdown for training sets of congenital myopathies, muscular dystrophies and metabolic myopathies disease groups.</caption><thead><tr><th colspan="11">Muscular dystrophies</th></tr><tr><th>Rank</th><th>Gene</th><th>Gene Score</th><th>GO score</th><th>%GO</th><th>HPO score</th><th>%HPO</th><th>MPO score</th><th>%MPO</th><th>%PO</th><th>IA score</th><th>%IA</th></tr></thead><tbody><tr><td>1</td><td>DMD</td><td>1115</td><td>352</td><td>31.57</td><td>180</td><td>16.14</td><td>375</td><td>33.63</td><td>49.78</td><td>208</td><td>18.65</td></tr><tr><td>2</td><td>TTN</td><td>1064</td><td>520</td><td>48.87</td><td>224</td><td>21.05</td><td>208</td><td>19.55</td><td>40.60</td><td>112</td><td>10.53</td></tr><tr><td>3</td><td>LMNA</td><td>997</td><td>208</td><td>20.86</td><td>361</td><td>36.21</td><td>284</td><td>28.49</td><td>64.69</td><td>144</td><td>14.44</td></tr><tr><td>4</td><td>TCAP</td><td>915</td><td>508</td><td>55.52</td><td>209</td><td>22.84</td><td>86</td><td>9.40</td><td>32.24</td><td>112</td><td>12.24</td></tr><tr><td>5</td><td>DES</td><td>814</td><td>252</td><td>30.96</td><td>195</td><td>23.96</td><td>223</td><td>27.40</td><td>51.35</td><td>144</td><td>17.69</td></tr><tr><td>6</td><td>CAV3</td><td>777</td><td>384</td><td>49.42</td><td>181</td><td>23.29</td><td>196</td><td>25.23</td><td>48.52</td><td>16</td><td>2.06</td></tr><tr><td>11</td><td>SYNE2</td><td>599</td><td>340</td><td>56.76</td><td>83</td><td>13.86</td><td>16</td><td>2.67</td><td>16.53</td><td>160</td><td>26.71</td></tr><tr><td>13</td><td>SYNE1</td><td>577</td><td>236</td><td>40.90</td><td>87</td><td>15.08</td><td>110</td><td>19.06</td><td>34.14</td><td>144</td><td>24.96</td></tr><tr><td>17</td><td>FKRP</td><td>560</td><td>48</td><td>8.57</td><td>333</td><td>59.46</td><td>163</td><td>29.11</td><td>88.57</td><td>16</td><td>2.86</td></tr><tr><td>18</td><td>DAG1</td><td>553</td><td>288</td><td>52.08</td><td>0</td><td>0.00</td><td>233</td><td>42.13</td><td>42.13</td><td>32</td><td>5.79</td></tr><tr><td>19</td><td>PLEC</td><td>552</td><td>140</td><td>25.36</td><td>65</td><td>11.78</td><td>171</td><td>30.98</td><td>42.75</td><td>176</td><td>31.88</td></tr><tr><td>20</td><td>CAPN3</td><td>539</td><td>176</td><td>32.65</td><td>123</td><td>22.82</td><td>176</td><td>32.65</td><td>55.47</td><td>64</td><td>11.87</td></tr><tr><td>22</td><td>TRIM32</td><td>518</td><td>92</td><td>17.76</td><td>299</td><td>57.72</td><td>95</td><td>18.34</td><td>76.06</td><td>32</td><td>6.18</td></tr><tr><td>23</td><td>SGCB</td><td>518</td><td>112</td><td>21.62</td><td>245</td><td>47.30</td><td>145</td><td>27.99</td><td>75.29</td><td>16</td><td>3.09</td></tr><tr><td>24</td><td>SGCA</td><td>503</td><td>136</td><td>27.04</td><td>163</td><td>32.41</td><td>156</td><td>31.01</td><td>63.42</td><td>48</td><td>9.54</td></tr><tr><td>25</td><td>EMD</td><td>487</td><td>196</td><td>40.25</td><td>180</td><td>36.96</td><td>79</td><td>16.22</td><td>53.18</td><td>32</td><td>6.57</td></tr><tr><td>28</td><td>SGCG</td><td>471</td><td>112</td><td>23.78</td><td>136</td><td>28.87</td><td>207</td><td>43.95</td><td>72.82</td><td>16</td><td>3.40</td></tr><tr><td>32</td><td>MYOT</td><td>446</td><td>176</td><td>39.46</td><td>270</td><td>60.54</td><td>0</td><td>0.00</td><td>60.54</td><td>0</td><td>0.00</td></tr><tr><td>36</td><td>DYSF</td><td>416</td><td>36</td><td>8.65</td><td>155</td><td>37.26</td><td>161</td><td>38.70</td><td>75.96</td><td>64</td><td>15.38</td></tr><tr><td>80</td><td>POMTT1</td><td>278</td><td>96</td><td>34.53</td><td>116</td><td>41.73</td><td>18</td><td>6.47</td><td>48.20</td><td>48</td><td>17.27</td></tr><tr><td>89</td><td>FHL1</td><td>270</td><td>32</td><td>11.85</td><td>238</td><td>88.15</td><td>0</td><td>0.00</td><td>88.15</td><td>0</td><td>0.00</td></tr><tr><td>92</td><td>SGCD</td><td>266</td><td>112</td><td>42.11</td><td>0</td><td>0.00</td><td>138</td><td>51.88</td><td>51.88</td><td>16</td><td>6.02</td></tr><tr><td>94</td><td>POMGNT1</td><td>265</td><td>76</td><td>28.68</td><td>92</td><td>34.72</td><td>81</td><td>30.57</td><td>65.28</td><td>16</td><td>6.04</td></tr><tr><td>114</td><td>ANOS</td><td>243</td><td>0</td><td>0.00</td><td>227</td><td>93.42</td><td>0</td><td>0.00</td><td>93.42</td><td>16</td><td>6.58</td></tr><tr><td>130</td><td>FKTN</td><td>226</td><td>84</td><td>37.17</td><td>112</td><td>49.56</td><td>14</td><td>6.19</td><td>55.75</td><td>16</td><td>7.08</td></tr><tr><td>193</td><td>PABPN1</td><td>181</td><td>40</td><td>22.10</td><td>77</td><td>42.54</td><td>0</td><td>0.00</td><td>42.54</td><td>64</td><td>35.36</td></tr><tr><td>262</td><td>POMT2</td><td>154</td><td>88</td><td>57.14</td><td>0</td><td>0.00</td><td>18</td><td>11.69</td><td>11.69</td><td>48</td><td>31.17</td></tr><tr><td>282</td><td>PTRF</td><td>149</td><td>16</td><td>10.74</td><td>124</td><td>83.22</td><td>8</td><td>5.37</td><td>88.59</td><td>0</td><td>0.00</td></tr><tr><td>443</td><td>DPM3</td><td>112</td><td>80</td><td>71.43</td><td>0</td><td>0.00</td><td>0</td><td>0.00</td><td>0.00</td><td>32</td><td>28.57</td></tr><tr><td colspan="2">Average</td><td>502.24</td><td>170.21</td><td>32.68</td><td>154.31</td><td>34.51</td><td>115.90</td><td>20.30</td><td>54.81</td><td>61.79</td><td>12.48</td></tr></tbody></table><table><thead><tr><th colspan="11">Congenital myopathies</th></tr><tr><th>Rank</th><th>Gene</th><th>Gene Score</th><th>GO score</th><th>%GO</th><th>HPO score</th><th>%HPO</th><th>MPO score</th><th>%MPO</th><th>%PO</th><th>IA score</th><th>%IA</th></tr></thead><tbody><tr><td>1</td><td>TTN</td><td>1309</td><td>641</td><td>48.97</td><td>206</td><td>15.74</td><td>315</td><td>24.06</td><td>39.80</td><td>147</td><td>11.23</td></tr></tbody></table>
<table><caption>Table 3. Cont.</caption><thead><tr><th colspan="11">Muscular dystrophies</th></tr><tr><th>Rank</th><th>Gene</th><th>Gene Score</th><th>GO score</th><th>%GO</th><th>HPO score</th><th>%HPO</th><th>MPO score</th><th>%MPO</th><th>%PO</th><th>IA score</th><th>%IA</th></tr></thead><tbody><tr><td>2</td><td>ACTA1</td><td>1291</td><td>555</td><td>42.99</td><td>441</td><td>34.16</td><td>169</td><td>13.09</td><td>47.25</td><td>126</td><td>9.76</td></tr><tr><td>3</td><td>RYR1</td><td>1162</td><td>349</td><td>30.03</td><td>292</td><td>25.13</td><td>332</td><td>28.57</td><td>53.70</td><td>189</td><td>16.27</td></tr><tr><td>4</td><td>NEB</td><td>1146</td><td>364</td><td>31.76</td><td>370</td><td>32.29</td><td>160</td><td>13.96</td><td>46.25</td><td>252</td><td>21.99</td></tr><tr><td>5</td><td>TPM3</td><td>1090</td><td>425</td><td>38.99</td><td>371</td><td>34.04</td><td>0</td><td>0.00</td><td>34.04</td><td>294</td><td>26.97</td></tr><tr><td>7</td><td>TPM2</td><td>998</td><td>393</td><td>39.38</td><td>311</td><td>31.16</td><td>0</td><td>0.00</td><td>31.16</td><td>294</td><td>29.46</td></tr><tr><td>11</td><td>MYH7</td><td>917</td><td>635</td><td>69.25</td><td>198</td><td>21.59</td><td>0</td><td>0.00</td><td>21.59</td><td>84</td><td>9.16</td></tr><tr><td>16</td><td>TRIM32</td><td>758</td><td>234</td><td>30.87</td><td>279</td><td>36.81</td><td>203</td><td>26.78</td><td>63.59</td><td>42</td><td>5.54</td></tr><tr><td>22</td><td>MTM1</td><td>702</td><td>53</td><td>7.55</td><td>211</td><td>30.06</td><td>354</td><td>50.43</td><td>80.48</td><td>84</td><td>11.97</td></tr><tr><td>27</td><td>TNNI1</td><td>663</td><td>333</td><td>50.23</td><td>141</td><td>21.27</td><td>0</td><td>0.00</td><td>21.27</td><td>189</td><td>28.51</td></tr><tr><td>33</td><td>DNM2</td><td>618</td><td>89</td><td>14.40</td><td>200</td><td>32.36</td><td>224</td><td>36.25</td><td>68.61</td><td>105</td><td>16.99</td></tr><tr><td>44</td><td>MYH2</td><td>527</td><td>443</td><td>84.06</td><td>0</td><td>0.00</td><td>0</td><td>0.00</td><td>0.00</td><td>84</td><td>15.94</td></tr><tr><td>55</td><td>BIN1</td><td>454</td><td>101</td><td>22.25</td><td>209</td><td>46.04</td><td>60</td><td>13.22</td><td>59.25</td><td>84</td><td>18.50</td></tr><tr><td>103</td><td>SEPN1</td><td>333</td><td>0</td><td>0.00</td><td>0</td><td>0.00</td><td>333</td><td>100.00</td><td>100.00</td><td>0</td><td>0.00</td></tr><tr><td>209</td><td>CNTN1</td><td>216</td><td>34</td><td>15.74</td><td>141</td><td>65.28</td><td>20</td><td>9.26</td><td>74.54</td><td>21</td><td>9.72</td></tr><tr><td>313</td><td>STM1</td><td>170</td><td>122</td><td>71.76</td><td>33</td><td>19.41</td><td>15</td><td>8.82</td><td>28.24</td><td>0</td><td>0.00</td></tr><tr><td>354</td><td>MTMR14</td><td>160</td><td>53</td><td>33.13</td><td>0</td><td>0.00</td><td>65</td><td>40.63</td><td>40.63</td><td>42</td><td>26.25</td></tr><tr><td colspan="2">Average</td><td>736.12</td><td>283.76</td><td>37.14</td><td>200.18</td><td>26.20</td><td>132.35</td><td>21.47</td><td>47.67</td><td>119.82</td><td>15.19</td></tr><tr><th colspan="11">Metabolic myopathies</th></tr><tr><th>Rank</th><th>Gene</th><th>Gene Score</th><th>GO score</th><th>%GO</th><th>HPO score</th><th>%HPO</th><th>MPO score</th><th>%MPO</th><th>%PO</th><th>IA score</th><th>%IA</th></tr><tr><td>1</td><td>GAA</td><td>875</td><td>410</td><td>46.86</td><td>136</td><td>15.54</td><td>309</td><td>35.31</td><td>50.86</td><td>20</td><td>2.29</td></tr><tr><td>2</td><td>PFKM</td><td>840</td><td>225</td><td>26.79</td><td>162</td><td>19.29</td><td>353</td><td>42.02</td><td>61.31</td><td>100</td><td>11.90</td></tr><tr><td>3</td><td>PHKA1</td><td>831</td><td>395</td><td>47.53</td><td>236</td><td>28.40</td><td>150</td><td>18.05</td><td>46.45</td><td>50</td><td>6.02</td></tr><tr><td>4</td><td>GBE1</td><td>822</td><td>330</td><td>40.15</td><td>59</td><td>7.18</td><td>333</td><td>40.51</td><td>47.69</td><td>100</td><td>12.17</td></tr><tr><td>6</td><td>ACADVL</td><td>747</td><td>140</td><td>18.74</td><td>275</td><td>36.81</td><td>222</td><td>29.72</td><td>66.53</td><td>110</td><td>14.73</td></tr><tr><td>7</td><td>GYS1</td><td>746</td><td>370</td><td>49.60</td><td>0</td><td>0.00</td><td>296</td><td>39.68</td><td>39.68</td><td>80</td><td>10.72</td></tr><tr><td>8</td><td>PGM1</td><td>709</td><td>450</td><td>63.47</td><td>159</td><td>22.43</td><td>0</td><td>0.00</td><td>22.43</td><td>100</td><td>14.10</td></tr><tr><td>12</td><td>AGL</td><td>666</td><td>480</td><td>72.07</td><td>136</td><td>20.42</td><td>0</td><td>0.00</td><td>20.42</td><td>50</td><td>7.51</td></tr><tr><td>14</td><td>PYGM</td><td>624</td><td>385</td><td>61.70</td><td>159</td><td>25.48</td><td>0</td><td>0.00</td><td>25.48</td><td>80</td><td>12.82</td></tr><tr><td>19</td><td>GYG1</td><td>590</td><td>480</td><td>81.36</td><td>0</td><td>0.00</td><td>0</td><td>0.00</td><td>0.00</td><td>110</td><td>18.64</td></tr><tr><td>24</td><td>LDHA</td><td>575</td><td>200</td><td>34.78</td><td>190</td><td>33.04</td><td>95</td><td>16.52</td><td>49.57</td><td>90</td><td>15.65</td></tr><tr><td>25</td><td>CPT2</td><td>571</td><td>205</td><td>35.90</td><td>316</td><td>55.34</td><td>0</td><td>0.00</td><td>55.34</td><td>50</td><td>8.76</td></tr><tr><td>28</td><td>PGAM2</td><td>551</td><td>265</td><td>48.09</td><td>216</td><td>39.20</td><td>0</td><td>0.00</td><td>39.20</td><td>70</td><td>12.70</td></tr><tr><td>37</td><td>PGK1</td><td>477</td><td>260</td><td>54.51</td><td>157</td><td>32.91</td><td>0</td><td>0.00</td><td>32.91</td><td>60</td><td>12.58</td></tr></tbody></table>
<table><caption>Table 3. Cont.</caption><thead><tr><th colspan="11">Muscular dystrophies</th></tr><tr><th>Rank</th><th>Gene</th><th>Gene Score</th><th>GO score</th><th>%GO</th><th>HPO score</th><th>%HPO</th><th>MPO score</th><th>%MPO</th><th>%PO</th><th>IA score</th><th>%IA</th></tr></thead><tbody><tr><td>47</td><td>ENO3</td><td>441</td><td>245</td><td>55.56</td><td>106</td><td>24.04</td><td>0</td><td>0.00</td><td>24.04</td><td>90</td><td>20.41</td></tr><tr><td>77</td><td>PNPLA2</td><td>354</td><td>10</td><td>2.82</td><td>0</td><td>0.00</td><td>334</td><td>94.35</td><td>94.35</td><td>10</td><td>2.82</td></tr><tr><td>84</td><td>LPIN1</td><td>344</td><td>25</td><td>7.27</td><td>141</td><td>40.99</td><td>158</td><td>45.93</td><td>86.92</td><td>20</td><td>5.81</td></tr><tr><td>87</td><td>SLC25A20</td><td>340</td><td>190</td><td>55.88</td><td>120</td><td>35.29</td><td>0</td><td>0.00</td><td>35.29</td><td>30</td><td>8.82</td></tr><tr><td>116</td><td>SLC22A5</td><td>311</td><td>60</td><td>19.29</td><td>123</td><td>39.55</td><td>108</td><td>34.73</td><td>74.28</td><td>20</td><td>6.43</td></tr><tr><td colspan="2">Average</td><td>600.74</td><td>269.74</td><td>43.28</td><td>141.63</td><td>25.05</td><td>124.11</td><td>20.89</td><td>45.93</td><td>65.26</td><td>10.78</td></tr></tbody></table>

Gene score is the sum of GO, HPO, MPO and IA scores. Relative contributions of GO, HPO, MPO and IA scores to the gene score are shown in the columns %GO, %HPO, %MPO and %IA, respectively. Training set genes without database annotation received a gene score of 0 and are not shown.
doi:10.1371/journal.pone.0110888.t003

Candidate genes for muscular dystrophies display strong links with muscle development, contraction and intracellular organization, expected subcomponents of skeletal muscle-related biological processes terms from the breakdown of GO terms. *ITGB1* codes for a subunit of ubiquitous fibronectin receptors and has a number of suggested functions in different tissues. In skeletal muscle, it has been proposed as a possible target for myostatin in mice myoblast differentiation [30] and is also critical for the development of neuromuscular junctions [31]. *TMOD1* encodes for tropomodulin, a protein that regulates tropomyosin and F-actin organization. Knockout mice present with age-dependent sarcomere misalignment and sarcoplasmic reticulum morphological defects [32]. *MYL1* is involved with early differentiation of fast muscle cells [33] and *TNNI1* codes for the slow-twitch skeletal muscle isoform of troponin I, which has yet to be associated with human diseases even though the fast-twitch isoform is responsible for a subtype of arthrogryposis and the cardiac isoform causes cardiomyopathy syndromes.

Candidate genes for congenital myopathies have a significant overlap with genes proposed for muscular dystrophies, for example for *TMOD1* and *TNNI1*. Among those with a high rank in congenital myopathies are *RYR3* and *MYH1*. *RYR3* codes for a ryanodine calcium release channel with a low $Ca^{2+}$ sensitivity that has a physiologic role in the excitation-contraction coupling of neonatal skeletal muscles and is up regulated in steroid-associated muscle damage [34], while *MYH1* is one of the adult skeletal muscle isoforms of myosin heavy chain that predominates in 2B myofibers. *RYR3* high ranking is boosted by a strong contribution of calcium homeostasis terms, explaining why *RYR3* received a similarly high score in the Myotonic Syndromes and Ion Channel Muscle Diseases groups, which also have a strong component of calcium homeostasis terms. Also in the group of Ion Channel Muscle Diseases, the gene *CACNB1* encodes both the brain and skeletal muscle isoforms of the calcium channel beta subunit, and its loss in mouse is associated to a phenotype similar to that seen in mice with mutations in the known genes *CACNA1S* or *RYR1* [35].

Within congenital muscular dystrophies, the gene *B3GALNT2*, ranked in the 97th place out of 4841 genes with annotation for this group's signature, was recently found to be associated with hypoglycosylation of alpha-dystroglycan and a congenital muscular dystrophy phenotype in humans [18]. Two other genes, *GMPPB*, ranked in 225th, and *B3GNT1*, ranked in 479th, were also implicated in a form of congenital muscular dystrophy with hypoglycosylation of alpha-dystroglycan and Walker-Warburg disease, respectively [36,37]. These genes had not been used in the training set of genes for congenital muscular dystrophies, and have since been annotated in OMIM, but their high placement in the ranking list validate the proposed data mining strategy and subsequent filtering steps.

## Candidate genes within genomic regions linked to myopathies

A number of neuromuscular diseases have mapped loci awaiting gene identification [1]. Matching the genomic positions of the top 100 candidate genes of each disease group with such mapped loci reveals some interesting candidates (Table 5).

The gene *XIRP1*, matching the locus for hyalin body myopathy and congenital muscular dystrophy with joint hyperlaxity, was originally studied in relation to murine cardiac morphogenesis and later shown to bind skeletal muscle actin in *in vitro* assays [38]. Its product, the Xin protein, is skeletal muscle-specific and has recently been put forward as a potentially useful biomarker of muscle damage, which can be used to monitor disease progression
<table><caption>**Table 4.** Top 8 ranked candidate genes for each disease group.</caption><thead><tr><th>Rank</th><th>Gene</th><th>Name</th><th>Score</th></tr></thead><tbody><tr><td>33</td><td>ITGB1</td><td>integrin, beta 1 (fibronectin receptor, beta polypeptide, antigen CD29 includes MDF2, MSK12)</td><td>437</td></tr><tr><td>42</td><td>TMOD1</td><td>troponulin 1</td><td>391</td></tr><tr><td>48</td><td>MYL1</td><td>myosin, light chain 1, alkali; skeletal, fast</td><td>368</td></tr><tr><td>53</td><td>TNNI1</td><td>troponin I type 1 (skeletal, slow)</td><td>356</td></tr><tr><td>62</td><td>MYH4</td><td>myosin, heavy chain 4, skeletal muscle</td><td>332</td></tr><tr><td>67</td><td>UTRN</td><td>utrophin</td><td>325</td></tr><tr><td>72</td><td>TNNC2</td><td>troponin C type 2 (fast)</td><td>304</td></tr><tr><td>81</td><td>SRF</td><td>serum response factor (c-fos serum response element-binding transcription factor)</td><td>278</td></tr></tbody></table><table><thead><tr><th>Rank</th><th>Gene</th><th>Name</th><th>Score</th></tr></thead><tbody><tr><td>36</td><td>GCNT4</td><td>glucosaminyl (N-acetyl) transferase 4, core 2</td><td>458</td></tr><tr><td>44</td><td>GALNT1</td><td>UDP-N-acetyl-alpha-D-galactosamine:polypeptide N-acetylgalactosaminyltransferase 1 (GalNAc-T1)</td><td>444</td></tr><tr><td>47</td><td>ST8SIA2</td><td>ST8 alpha-N-acetyl-neuraminide alpha-2,8-sialyltransferase 2</td><td>444</td></tr><tr><td>51</td><td>OGT</td><td>O-linked N-acetylglucosamine (GlcNAc) transferase</td><td>444</td></tr><tr><td>53</td><td>GALNT2</td><td>UDP-N-acetyl-alpha-D-galactosamine:polypeptide N-acetylgalactosaminyltransferase 2 (GalNAc-T2)</td><td>444</td></tr><tr><td>55</td><td>SDF2</td><td>stromal cell-derived factor 2</td><td>443</td></tr><tr><td>57</td><td>ST8SIA6</td><td>ST8 alpha-N-acetyl-neuraminide alpha-2,8-sialyltransferase 6</td><td>442</td></tr><tr><td>62</td><td>MGAT1</td><td>mannosyl (alpha-1,3-)-glycoprotein beta-1,2-N-acetylglycosaminyltransferase</td><td>426</td></tr></tbody></table><table><thead><tr><th>Rank</th><th>Gene</th><th>Name</th><th>Score</th></tr></thead><tbody><tr><td>17</td><td>MYH4</td><td>myosin, heavy chain 4, skeletal muscle</td><td>755</td></tr><tr><td>29</td><td>MYL1</td><td>myosin, light chain 1, alkali; skeletal, fast</td><td>640</td></tr><tr><td>31</td><td>TNNI1</td><td>troponin I type 1 (skeletal, slow)</td><td>634</td></tr><tr><td>37</td><td>RYR3</td><td>ryanodine receptor 3</td><td>587</td></tr><tr><td>38</td><td>TMOD1</td><td>troponulin 1</td><td>570</td></tr><tr><td>43</td><td>TNNC2</td><td>troponin C type 2 (fast)</td><td>527</td></tr><tr><td>51</td><td>MYH1</td><td>myosin, heavy chain 1, skeletal muscle, adult</td><td>458</td></tr><tr><td>58</td><td>MYL6B</td><td>myosin, light chain 6B, alkali, smooth muscle and non-muscle</td><td>438</td></tr></tbody></table><table><thead><tr><th>Rank</th><th>Gene</th><th>Name</th><th>Score</th></tr></thead><tbody><tr><td>10</td><td>PRKAA2</td><td>protein kinase, AMP-activated, alpha 2 catalytic subunit</td><td>672</td></tr><tr><td>18</td><td>PPP1R3C</td><td>protein phosphatase 1, regulatory subunit 3C</td><td>593</td></tr><tr><td>21</td><td>MTOR</td><td>mechanistic target of rapamycin (serine/threonine kinase)</td><td>588</td></tr><tr><td>26</td><td>PRKAB2</td><td>protein kinase, AMP-activated, beta 2 non-catalytic subunit</td><td>560</td></tr><tr><td>39</td><td>ACACB</td><td>acetyl-CoA carboxylase beta</td><td>470</td></tr><tr><td>40</td><td>PHKG1</td><td>phosphorylase kinase, gamma 1 (muscle)</td><td>470</td></tr><tr><td>44</td><td>PPARGC1A</td><td>peroxisome proliferator-activated receptor gamma, coactivator 1 alpha</td><td>459</td></tr><tr><td>48</td><td>GSK3A</td><td>glycogen synthase kinase 3 alpha</td><td>441</td></tr></tbody></table><table><thead><tr><th>Rank</th><th>Gene</th><th>Name</th><th>Score</th></tr></thead><tbody></tbody></table>
<table><caption>**Table 4.** Cont.</caption><thead><tr><th>Rank</th><th>Gene</th><th>Name</th><th>Score</th></tr></thead><tbody><tr><td colspan="4"><strong>Muscular dystrophies</strong></td></tr><tr><td>11</td><td>CHRNB4</td><td>cholinergic receptor, nicotinic, beta 4 (neuronal)</td><td>817</td></tr><tr><td>15</td><td>CHRNA6</td><td>cholinergic receptor, nicotinic, alpha 6 (neuronal)</td><td>680</td></tr><tr><td>19</td><td>CHRNB3</td><td>cholinergic receptor, nicotinic, beta 3 (neuronal)</td><td>636</td></tr><tr><td>21</td><td>CACNA2D2</td><td>calcium channel, voltage-dependent, alpha 2/delta subunit 2</td><td>628</td></tr><tr><td>22</td><td>CHRNA9</td><td>cholinergic receptor, nicotinic, alpha 9 (neuronal)</td><td>616</td></tr><tr><td>40</td><td>ITGB1</td><td>integrin, beta 1 (fibronectin receptor, beta polypeptide, antigen CD29 includes MDF2, MSK12)</td><td>441</td></tr><tr><td>52</td><td>CHRNA10</td><td>cholinergic receptor, nicotinic, alpha 10 (neuronal)</td><td>413</td></tr><tr><td>61</td><td>HTR3B</td><td>5-hydroxytryptamine (serotonin) receptor 3B, ionotropic</td><td>399</td></tr><tr><td colspan="4"><strong>Ion channel muscle diseases</strong></td></tr><tr><td>12</td><td>SCN3A</td><td>sodium channel, voltage-gated, type III, alpha subunit</td><td>1416</td></tr><tr><td>13</td><td>CACNB1</td><td>calcium channel, voltage-dependent, beta 1 subunit</td><td>1378</td></tr><tr><td>40</td><td>RYR3</td><td>ryanodine receptor 3</td><td>1019</td></tr><tr><td>52</td><td>CACNG1</td><td>calcium channel, voltage-dependent, gamma subunit 1</td><td>936</td></tr><tr><td>55</td><td>CACNA2D1</td><td>calcium channel, voltage-dependent, alpha 2/delta subunit 1</td><td>930</td></tr><tr><td>63</td><td>CACNA2D3</td><td>calcium channel, voltage-dependent, alpha 2/delta subunit 3</td><td>915</td></tr><tr><td>70</td><td>KCNQ5</td><td>potassium voltage-gated channel, KQT-like subfamily, member 5</td><td>888</td></tr><tr><td>71</td><td>KCNA7</td><td>potassium voltage-gated channel, shaker-related subfamily, member 7</td><td>888</td></tr><tr><td colspan="4"><strong>Myotonic syndromes</strong></td></tr><tr><td>7</td><td>CASQ1</td><td>calsequestrin 1 (fast-twitch, skeletal muscle)</td><td>1079</td></tr><tr><td>8</td><td>RYR3</td><td>ryanodine receptor 3</td><td>954</td></tr><tr><td>15</td><td>JPH1</td><td>junctophilin 1</td><td>806</td></tr><tr><td>21</td><td>MYL1</td><td>myosin, light chain 1, alkali; skeletal, fast</td><td>703</td></tr><tr><td>26</td><td>CAMK2D</td><td>calcium/calmodulin-dependent protein kinase II delta</td><td>658</td></tr><tr><td>32</td><td>SYPL2</td><td>synaptophysin-like 2</td><td>611</td></tr><tr><td>34</td><td>ITGB1</td><td>integrin, beta 1 (fibronectin receptor, beta polypeptide, antigen CD29 includes MDF2, MSK12)</td><td>610</td></tr><tr><td>47</td><td>MYH4</td><td>myosin, heavy chain 4, skeletal muscle</td><td>559</td></tr><tr><td colspan="4"><strong>Myofibrillar myopathies</strong></td></tr><tr><td>27</td><td>MYL1</td><td>myosin, light chain 1, alkali; skeletal, fast</td><td>984</td></tr><tr><td>33</td><td>MYH4</td><td>myosin, heavy chain 4, skeletal muscle</td><td>924</td></tr><tr><td>35</td><td>MYL12B</td><td>myosin, light chain 12B, regulatory</td><td>921</td></tr><tr><td>41</td><td>TNNI1</td><td>troponin I type 1 (skeletal, slow)</td><td>909</td></tr><tr><td>42</td><td>TNNC2</td><td>troponin C type 2 (fast)</td><td>909</td></tr><tr><td>50</td><td>PDLIM3</td><td>PDZ and LIM domain 3</td><td>862</td></tr><tr><td>51</td><td>MYO18B</td><td>myosin XVIIIB</td><td>845</td></tr><tr><td>52</td><td>PDLIM5</td><td>PDZ and LIM domain 5</td><td>844</td></tr><tr><td colspan="4"><strong>Vacuolar myopathies</strong></td></tr><tr><td>4</td><td>CD63</td><td>CD63 molecule</td><td>1160</td></tr><tr><td>8</td><td>AP1G1</td><td>adaptor-related protein complex 1, gamma 1 subunit</td><td>1006</td></tr></tbody></table>
<table><caption>**Table 4.** Cont.</caption><thead><tr><th colspan="5">Muscular dystrophies</th></tr><tr><th>Rank</th><th>Gene</th><th>Name</th><th>Score</th><th></th></tr></thead><tbody><tr><td>18</td><td>VAMP7</td><td>vesicle-associated membrane protein 7</td><td>1006</td><td></td></tr><tr><td>20</td><td>MARCH8</td><td>membrane-associated ring finger (C3HC4) 8, E3 ubiquitin protein ligase</td><td>1006</td><td></td></tr><tr><td>22</td><td>ZNRF1</td><td>zinc and ring finger 1, E3 ubiquitin protein ligase</td><td>1006</td><td></td></tr><tr><td>26</td><td>AP1M1</td><td>adaptor-related protein complex 1, mu 1 subunit</td><td>1006</td><td></td></tr><tr><td>27</td><td>AP1B1</td><td>adaptor-related protein complex 1, beta 1 subunit</td><td>1006</td><td></td></tr><tr><td>40</td><td>ABCC4</td><td>ATP-binding cassette, sub-family C (CFTR/MRP), member 4</td><td>817</td><td></td></tr></tbody></table>

Candidate genes are not associated with disease (as per annotation in OMIM) and are expressed in skeletal muscle with at least 10% of the maximum expression in any tissue, except for congenital myasthenic syndromes, where there was no expression filtering.

doi:10.1371/journal.pone.0110888.t004

and treatment effects in myopathies [39]. *OBSCN* encodes obscurin, a giant sarcomeric signaling protein similar to titin, which has a suspected role in myofibrilogenesis. It is also involved in dystrophin localization and maintenance of sarcolemma integrity [40], and is proposed here as a candidate for congenital muscular dystrophy with merosin deficiency (MDC1B). An additional candidate gene mapped in the linked region is *GALNT2*, a glycosylating enzyme similar to *B3GALNT2* recently found mutated in another form of congenital muscular dystrophy [18], and also involved in the O-glycosylation of peptides in the Golgi apparatus.

Although not directly analyzed in this work, malignant hyperthermia susceptibility regions encompass *CACNG1* and *CACNA2D1*, which are associated with calcium homeostasis and calcium channels, are highly ranked for Ion Channel Muscle Diseases, and are thus interesting candidates. *CACNA2D1*, specifically, has been suggested at least as a modifier of hyperthermia susceptibility in association to other genes [41]. These genes have been excluded in a limited number of families not linked to *RYR1* mutations [42], results which may be revisited with the advent of NGS data. Additionally, another candidate gene, *CACNB1*, has no associated human disease. However, *CACNG1*, *CACNB1* and *CACNA2D1* encode for subunits of the DHPR calcium channel, which is in direct contact and regulating *RYR1* in skeletal muscle, and one mutation in the channel subunit *CACNA1S* of DHPR was linked to malignant hyperthermia [43].

The *CALR* and *AP1M1* genes both map to 19q13, the locus associated to autosomal muscular dystrophy with rimmed vacuoles. In a recent work, the product of the *CALR* gene, calreticulin, has been shown to localize in cardiomyocyte mitochondria, and its content increases in mouse models with dilated cardiomyopathy [44]. Strikingly, calreticulin was found to be highly expressed in GNE myopathy, a distal myopathy associated with rimmed vacuoles [45]. Also in distal myopathies with rimmed vacuoles, though not necessarily GNE myopathy, adaptin related-proteins subunits, which are normally not marked in the immunohistochemistry of normal muscle, appear inside or on the rims of vacuoles. The *AP1M1* gene codes for the mu subunit of adaptin related-proteins [46].

## Discussion

In this study, an integrated data mining strategy was used to cluster and rank genes with known or potential importance for skeletal muscle, and to provide candidate genes for myopathies and some related diseases. Results from the clustering and ranking highlight pathological pathways specific for disease groups. The list

of candidate genes was further filtered based on expression data and association with other diseases, and the ensuing identification of mutations in high-ranked genes for congenital muscular dystrophies (*B3GALNT2*, *GMPPB* and *B3GNT1*) illustrated the validity of this approach.

## Gene clustering and ranking are dependent on database annotations

Proposed genes in the final ranked list have gene scores with a major contribution of GO and IA terms, and eventual contribution of MPO terms. Thus, they represent genes that have mostly functional links with known myopathy genes (IA terms and GO term ontologies for biological processes and molecular functions), but also some degree of product colocalization in the muscle cell, as expected from matching cellular component-related GO terms. When available, data on altered skeletal muscle function in mouse models also tend to contribute to higher scores for proposed candidate genes.

Database annotation can vary from one gene to another, as it is dependent on the history of research for each gene, including both the date when the gene was discovered and the amount of effort spent for its functional characterization. In addition, animal models are generally phenotyped with a specific organ system in mind.

To illustrate the effects of incomplete annotation, the genes *TRAPPC11* and *TNPO3*, recently implicated in muscular dystrophies, were used as components of this disease group training set, but did not impact the results of the gene ranking due to their poor database annotation. Likewise, they were themselves not captured by the signature used for the Muscular Dystrophies disease group. *TRAPPC11* does not appear in the ranking, as it was annotated with only two GO terms that are not significant for muscular dystrophies ("vesicle-mediated transport" and "Golgi apparatus"), it has no annotation for pathways or phenotypes, and its two protein motifs are unique in the genome. Annotation biases also account for higher placements of better-annotated genes that have some kind of overlap with myopathy genes. Such is the case for motor neuron disease-associated genes, which give rise to human and mouse phenotypes that present some degree of phenotypic overlap with myopathies and tend to share many HPO or MPO terms with myopathy phenotypes. In the ranked list of muscular dystrophies, high scores with a predominant contribution of HPO terms were given to the genes *SMN1*, *SMN2*, *ALS2*, *IGHMBP2* and *AR*. These genes are linked to different types of motor neuron disease, which ultimately manifest with muscle weakness and atrophy.
<table><caption>**Table 5.** Candidate genes within genomic regions linked to myopathies and related diseases.</caption><thead><tr><th>Linked region</th><th>Phenotypes and associated disease symbols</th><th>Candidate genes</th></tr></thead><tbody><tr><td>1q42</td><td>Congenital muscular dystrophy with merosin deficiency - MDC1B</td><td>OBSCN, GALNT2</td></tr><tr><td>3p22.2-p21.32</td><td>Hyalin body myopathy - HBM</td><td>XIRP1</td></tr><tr><td>3p23-21</td><td>Congenital muscle dystrophy with joint hyperlaxity</td><td>XIRP1</td></tr><tr><td>7q21-q22</td><td>Malignant hyperthermia susceptibility 3 - MHS3</td><td>CACNA2D1</td></tr><tr><td>17q11.2-q24</td><td>Malignant hyperthermia susceptibility 2 - MHS2</td><td>SDF2, SYNRG, CACNB1, CACNG1</td></tr><tr><td>19p13</td><td>Muscular dystrophy, autosomal dominant, with rimmed vacuoles - MDRV</td><td>CALR, PRKACA, AP1M1</td></tr></tbody></table>

doi:10.1371/journal.pone.0110888.t005

*In silico* approaches need thus periodic revisits to adjust candidate lists based on association of new genes that impact training sets and discovery of new pathways or interactions that change corresponding database annotation, such as the recently published interactome of skeletal muscle proteins centered on proteins that cause limb-girdle muscular dystrophies [47].

## Possible insights into pathological mechanisms

The integrated data mining identified gene signatures revealing common function within specific myopathy groups or between groups, and highlighting pathological mechanisms.

Known and candidate genes for metabolic myopathies, congenital myasthenic syndromes, myotonic syndromes, and ion channel muscle diseases define distinctive functions for each disease group (Table 2). Highly ranked genes for congenital myasthenic syndromes are associated to function not primarily linked to skeletal muscle but point as expected to the neuromuscular junction. The cellular basis of myotonic myopathies and ion channel muscle diseases consists in the alteration of ion homeostasis. Additional genes contributing to glycogen metabolism were identified as good candidate genes for the metabolic myopathies.

Muscle development, muscle contraction and calcium homeostasis are key pathways linked to congenital myopathies; indeed this myopathy group presents generally at or before birth, and is characterized by histological hallmarks reflecting alteration and aggregation of proteins implicated in muscle contraction (nemaline myopathies) or due to primary defects in the excitation-contraction coupling (core myopathies and potentially the centronuclear myopathies).

Muscular dystrophies mainly involve muscle intracellular organization terms that reflect the structural importance of most proteins already reported mutated. However, other pathways may have been overlooked because the way the first genes were discovered; once the *DMD* gene was found, investigators started seeking mutations on genes from the dystrophin-glycoprotein complex. Of note, based on the terms breakdown and candidate genes identified, muscular dystrophies may have a larger contribution from the contractile apparatus than previously assumed, which would bring this disease group closer to congenital myopathies.

## Allelic diseases

The integrated data mining reveals or confirms allelic diseases. Indeed, while proposed genes for metabolic myopathies or myasthenic syndromes are rather group-specific, a larger overlap occurred between congenital myopathies and muscular dystro-

phies than what was expected from the analysis of overlap between these groups' training sets. While only 2 genes out of 34 muscular dystrophy training set genes also appeared among the 22 congenital myopathy training set genes (*TTN* and *TRIM32*), the first positions on the ranked lists of candidate genes after filtering for known diseases encompass a large overlap of genes: 5 out of the top 8 candidate genes for muscular dystrophies are within the top 8 for congenital myopathies (Table 4), and 33 out of the top 50 candidate muscular dystrophy genes are also within the top 50 candidate genes for congenital myopathies (Figure 3 and Table S2). Overlaps are also substantial between both these lists and the one for myofibrillar myopathies, but in this case the overlap was expected as the training set for myofibrillar myopathies, despite being small (7 genes), included 2 genes also associated to muscular dystrophies and 1 gene associated to congenital myopathies. On the other hand, while the training sets of muscular dystrophies and congenital muscular dystrophies overlap with a significant share of 8 genes, only 3 genes within the top 50 candidate genes is the same for both groups. The reasons for these results lie in the signature of the disease groups: the breakdown of biological processes terms (depicted in Table 2), which represent the larger share of GO terms, is more comparable between congenital myopathies and muscular dystrophies, with similar contributions of terms involving muscle contraction and development, as opposed to the absence of resemblance between these disease groups and the congenital muscular dystrophies breakdown, which is enriched with mainly non-skeletal muscle-related terms, especially glycosylation. Taken together, gene clustering and candidate genes retrieval suggest that mutations in the same genes will eventually be linked to both muscular dystrophies and congenital myopathies.

Table 4 presents the top 8 ranked genes after excluding genes with known disease annotation. Genes with known disease annotation (listed in Table S2) might still be good candidate genes for myopathies, considering that phenotypic variability is more a rule than an exception for known myopathy genes. This is the case for genes linked to both myopathy and cardiomyopathy. It is thus expected that known cardiomyopathy-associated genes may be found associated with a skeletal muscle phenotype. Such phenotypic variability may even transcend the realm of muscle alteration. *DNM2*, for example, is associated both with centronuclear myopathy, a congenital myopathy, and Charcot-Marie Tooth disease, a hereditary neuropathy [48,49]. *LMNA*, in addition to multiple myopathic phenotypes, also causes Charcot-Marie Tooth disease or progeria [50,51], and *SYNE1* can cause one type of Emery-Dreifuss muscular dystrophy, a dilated cardiomyopathy syndrome, a form of autosomal recessive arthrogryposis, and autosomal recessive spinocerebellar ataxia
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Venn diagrams of gene set overlaps. A: Venn diagram showing the overlap of training set genes between muscular dystrophies (MD), congenital myopathies (CM) and congenital muscular dystrophies (CMD). B: Venn diagram showing the overlap of genes found within the top 50 ranked candidate genes in the three disease groups.
doi:10.1371/journal.pone.0110888.g003</figcaption></figure>

[52–55]. The origin of the variability may stem from varying impacts of mutations in different protein domains. For example, in *DNM2*, mutations giving rise to a centronuclear myopathy phenotype are enriched in the interface between the middle domain and the pleckstrin homology domain, while mutations implicated in Charcot-Marie-Tooth disease tend to cluster in other parts of the pleckstrin homology domain [56].

After we carried out this work, a form of vacuolar myopathy was
associated to *CLN3*, implicated in neuronal ceroid lipofuscinosis,
ranked 343rd for vacuolar myopathies [57]. In addition, *LRP4*, a
gene associated with Cenani-Lenz syndactyly syndrome and
sclerosteosis, and ranked high (54th) for congenital muscular
dystrophies, was implicated in a patient with a CMS disease [58].
We therefore suggest that genes with disease annotation in the
ranked lists should be considered with caution in analyzing NGS
results, but not *a priori* excluded when using filtering pipelines.

Example of the usage of the ranked lists

We found the ranked lists to be helpful in our own analysis of
exome data to prioritize the scrutiny of potential novel genes
implicated in myopathies. The gene ranks in the Excel file sheets
in Table S2 can be easily used as additional genomic annotation.

Consider this exome of a sporadic patient affected with nemaline myopathy, a congenital myopathy, from unaffected parents (Table S4). Out of an initial 86,333 variants called in the exome data, 250 remained after variants filtering to exclude purported sequencing errors and polymorphisms. The first variants we analyze closely are the ones found in genes with known implication in myopathies. The heterozygous variants in *DMD* and *CACNA1S* were subsequently found in the unaffected father, while the missense variant in *ANO5*, associated with autosomal recessive limb girdle dystrophy 2L, would require an association to a second mutation to cause disease. We can thus exclude the implication of these known genes.

We next focus on the candidate genes for congenital myopathies
(Table 4 and Table S2). If a gene has more than 10% expression
in skeletal muscle and is not associated to a disease, it receives a

flag as a “candidate”. Only 29 variants, a significant reduction
from the original 250, survive this additional filtering and are
shown in Table 6.

The gene with the highest rank was *PXM*, which codes for paxillin, a protein believed to have a function related to integrins and cytoskeletal localization in multiple tissues, skeletal muscle included [59]. The second best gene was *MYBPH*, which codes for myosin-binding protein H, the second most abundant protein of the family of myosin-associated proteins [60]. Except for its cloning, not much is known about its specific function. The third gene, *PPARGC1A*, has regulatory functions on glucose and fat oxidation in muscle cells and protects skeletal muscle fibers against atrophy in mouse models [61]. However, in all these genes, single mutations were found in a heterozygous state, thus a putative dominant negative effect or haploinsufficiency would be required for a pathogenicity call. The next gene in the list, *SYNPO*, produces synaptopodin, a protein whose name stems from its involvement in synapses involving dendritic spines, in addition to renal podocytes [62]. In spite of its name, skeletal muscle is actually the tissue where it is most strongly expressed. Furthermore, synaptopodin directly binds actin, one of the proteins known to be involved in nemaline myopathy. The missense variant found in this gene was homozygous, in a highly evolutionary conserved position, with prediction of pathogenicity in multiple tools, and was Sanger confirmed to be homozygous in the patient and heterozygous in her parents. We believe *SYNPO* is the best candidate gene for this family, based on a recessive scenario. While we cannot discard the other genes, the ranking of candidate genes based on our integrative data mining quickly highlights the best genes to proceed further in functional analysis.

Conclusions

The above integrated data mining approach was successfully
used to retrieve both specific signatures for different myopathy
groups and to uncover and rank interesting candidate genes for
myopathies. Recent discoveries of gene implications that were
<table><caption>**Table 6.** Resulting 29 variants after filtration of exome data of a patient affected with nemaline myopathy.</caption><thead><tr><th>CM Rank</th><th>Flag</th><th>VariantID</th><th>State</th><th>Gene</th><th>Spec%</th></tr></thead><tbody><tr><td>208</td><td>candidate</td><td>12_120660719_C_T</td><td>Heterozygous</td><td>PXN</td><td>22</td></tr><tr><td>355</td><td>candidate</td><td>1_203139425_T_A</td><td>Heterozygous</td><td>MYBPH</td><td>16</td></tr><tr><td>446</td><td>candidate</td><td>4_23803919_C_T</td><td>Heterozygous</td><td>PPARGC1A</td><td>30</td></tr><tr><td>586</td><td>candidate</td><td>5_150028613_T_C</td><td>Homozygous</td><td>SYNPO</td><td>100</td></tr><tr><td>610</td><td>candidate</td><td>5_138160333_G_A</td><td>Heterozygous</td><td>CTNNA1</td><td>34</td></tr><tr><td>786</td><td>candidate</td><td>10_115374035_A_T</td><td>Heterozygous</td><td>NRAP</td><td>100</td></tr><tr><td>856</td><td>candidate</td><td>1_87208087_A_G</td><td>Heterozygous</td><td>SH3GLB1</td><td>100</td></tr><tr><td>951</td><td>candidate</td><td>6_36076169_A_G</td><td>Heterozygous</td><td>MAPK14</td><td>18</td></tr><tr><td>1044</td><td>candidate</td><td>11_1901435_C_T</td><td>Heterozygous</td><td>LSP1</td><td>36</td></tr><tr><td>1044</td><td>candidate</td><td>11_1901435_C_T</td><td>Heterozygous</td><td>LSP1</td><td>36</td></tr><tr><td>1199</td><td>candidate</td><td>9_125863896_C_T</td><td>Heterozygous</td><td>RABGAP1</td><td>38</td></tr><tr><td>1758</td><td>candidate</td><td>12_95604081_G_A</td><td>Heterozygous</td><td>FGD6</td><td>17</td></tr><tr><td>1902</td><td>candidate</td><td>14_103420979_G_A</td><td>Heterozygous</td><td>CDC42BPB</td><td>32</td></tr><tr><td>1976</td><td>candidate</td><td>9_124522285_C_T</td><td>Heterozygous</td><td>DAB2IP</td><td>43</td></tr><tr><td>2066</td><td>candidate</td><td>22_19865895_A_C</td><td>Heterozygous</td><td>TXNRD2</td><td>16</td></tr><tr><td>2231</td><td>candidate</td><td>5_95116054_A_T</td><td>Heterozygous</td><td>RHOTB3</td><td>32</td></tr><tr><td>2245</td><td>candidate</td><td>22_41652800_A_C</td><td>Heterozygous</td><td>RANGAP1</td><td>25</td></tr><tr><td>2263</td><td>candidate</td><td>2_159477732_C_A</td><td>Heterozygous</td><td>PKP4</td><td>10</td></tr><tr><td>2360</td><td>candidate</td><td>2_152980460_G_T</td><td>Heterozygous</td><td>STAM2</td><td>26</td></tr><tr><td>2679</td><td>candidate</td><td>1_46472006_A_G</td><td>Heterozygous</td><td>MAST2</td><td>100</td></tr><tr><td>3075</td><td>candidate</td><td>10_68138967_C_T</td><td>Heterozygous</td><td>CTNNA3</td><td>19</td></tr><tr><td>3434</td><td>candidate</td><td>7_156976610_G_A</td><td>Heterozygous</td><td>UBE3C</td><td>100</td></tr><tr><td>3530</td><td>candidate</td><td>9_32407367_C_T</td><td>Heterozygous</td><td>ACO1</td><td>16</td></tr><tr><td>3627</td><td>candidate</td><td>1_19453077_C_T</td><td>Heterozygous</td><td>UBR4</td><td>100</td></tr><tr><td>4029</td><td>candidate</td><td>20_35632140_C_G</td><td>Heterozygous</td><td>RBL1</td><td>40</td></tr><tr><td>4235</td><td>candidate</td><td>22_50356432_A_T</td><td>Heterozygous</td><td>PIM3</td><td>47</td></tr><tr><td>4375</td><td>candidate</td><td>11_75115893_C_A</td><td>Heterozygous</td><td>RPS3</td><td>14</td></tr><tr><td>5062</td><td>candidate</td><td>1_204494668_G_A</td><td>Heterozygous</td><td>MDM4</td><td>47</td></tr><tr><td>5084</td><td>candidate</td><td>7_21469915_C_T</td><td>Heterozygous</td><td>SP4</td><td>36</td></tr></tbody></table>

An initial 86,333 variants were reduced to 250 using criteria on the variant level, which resulted in the 29 variants after exclusion of genes already ascribed to diseases and based on specificity of skeletal muscle expression. Variants are then sorted according to the gene ranking calculated for the congenital myopathy group.
doi:10.1371/journal.pone.0110888.t006

correctly identified by the disease group's signature validated this approach. *In silico* approaches allow for systematic, but modifiable criteria to be used in generating ranked candidate lists and have the added benefit of automation, whereby such lists can be updated on the fly as new knowledge is incorporated in genomic databases.

Signatures and candidate genes highlighted both potential common pathological mechanisms and overlap between several disease groups. In addition, the ranked candidate gene lists are helpful to prioritize functional validation of filtered variants from overwhelming NGS data.

## Supporting Information

**Table S1 Breakdown of descriptor terms for every domain of each disease group, with corresponding calculated weights.**

(XLSX)

**Table S2 Full ranked gene lists for each disease group.**

(XLSX)

**Table S3 Ranked lists of known genes for each disease group.**

(XLSX)

**Table S4 Filtered variants from an exome of a patient with nemaline myopathy, ordered according to the ranked gene lists for congenital myopathies.**

(XLSX)

## Author Contributions

Conceived and designed the experiments: OAN OT OP JL. Performed the experiments: OAN. Analyzed the data: OAN OT VB EZ OP JL. Contributed reagents/materials/analysis tools: OAN OT OP. Contributed to the writing of the manuscript: OAN JL.
References

1. Kaplan JC, Hamroun D (2013) The 2014 version of the gene table of monogenic neuromuscular disorders (nuclear genome). Neuromuscul Disord 23: 1081-1111.

2. Mercuri E, Muntoni F (2013) Muscular dystrophies. Lancet 381: 845-860.

3. Nance JR, Dowling JJ, Gibbs EM, Bonnemann CG (2012) Congenital myopathies: an update. Curr Neurol Neurosci Rep 12: 165-174.

4. Mercuri E, Muntoni F (2012) The ever-expanding spectrum of congenital muscular dystrophies. Ann Neurol 72: 9-17.

5. Heatwole CR, Statland JM, Logigian EL (2013) The diagnosis and treatment of myotonic disorders. Muscle Nerve 47: 632-648.

6. Amburgey K, McNamara N, Bennett LR, McCormick ME, Acsadi G, et al. (2011) Prevalence of congenital myopathies in a representative pediatric united states population. Ann Neurol 70: 662-665.

7. Cuisset JM, Maurage CA, Carpentier A, Briand G, Thevenon A, et al. (2013) [Muscle biopsy in children: Usefulness in 2012]. Rev Neurol (Paris) 169: 632-639.

8. Bohm J, Vasli N, Malfatti E, Le Gras S, Feger C, et al. (2013) An integrated diagnosis strategy for congenital myopathies. PLoS One 8: e67527.

9. Vasli N, Laporte J (2013) Impacts of massively parallel sequencing for genetic diagnosis of neuromuscular disorders. Acta Neuropathol 125: 173-185.

10. Nigro V, Piluso G (2012) Next generation sequencing (NGS) strategies for the genetic testing of myopathies. Acta Myol 31: 196-200.

11. Boycott KM, Vanstone MR, Bulman DE, MacKenzie AE (2013) Rare-disease genetics in the era of next-generation sequencing: discovery to translation. Nat Rev Genet 14: 681-691.

12. Ng SB, Buckingham KJ, Lee C, Bigham AW, Tabor HK, et al. (2010) Exome sequencing identifies the cause of a mendelian disorder. Nat Genet 42: 30-35.

13. Bogershausen N, Shahrzad N, Chong JX, von Kleist-Retzow JC, Stanga D, et al. (2013) Recessive TRAPP1C1 Mutations Cause a Disease Spectrum of Limb Girdle Muscular Dystrophy and Myopathy with Movement Disorder and Intellectual Disability. Am J Hum Genet 93: 181-190.

14. Cullup T, Kho AL, Dionisi-Vici C, Brandmeier B, Smith F, et al. (2013) Recessive mutations in EPG5 cause Vici syndrome, a multisystem disorder with defective autophagy. Nat Genet 45: 83-87.

15. Majczenko K, Davidson AE, Camelo-Piragua S, Agrawal PB, Manfready RA, et al. (2012) Dominant mutation of CCDC78 in a unique congenital myopathy with prominent internal nuclei and atypical cores. Am J Hum Genet 91: 365-371.

16. Ramachandran N, Munteanu I, Wang P, Ruggieri A, Rilstone JJ, et al. (2013) VMA21 deficiency prevents vacuolar ATPase assembly and causes autophagic vacuolar myopathy. Acta Neuropathol 125: 439-457.

17. Ravenscroft G, Miyatake S, Lehtokari VL, Todd EJ, Vornanen P, et al. (2013) Mutations in KLHL40 Are a Frequent Cause of Severe Autosomal-Recessive Nemaline Myopathy. Am J Hum Genet 93: 6-18.

18. Stevens E, Carss KJ, Cirak S, Foley AR, Torelli S, et al. (2013) Mutations in B3GALNT2 cause congenital muscular dystrophy and hypoglycosylation of alpha-dystroglycan. Am J Hum Genet 92: 354-365.

19. Torella A, Fanin M, Mutarelli M, Peterle E, Del Vecchio Blanco F, et al. (2013) Next-generation sequencing identifies transportin 3 as the causative gene for LGMD1F. PLoS One 8: e63536.

20. Bohm J, Chevessier F, Maues De Paula A, Koch C, Attarian S, et al. (2013) Constitutive activation of the calcium sensor STIM1 causes tubular-aggregate myopathy. Am J Hum Genet 92: 271-278.

21. Consortium EP (2012) An integrated encyclopedia of DNA elements in the human genome. Nature 499: 57-74.

22. Landrum MJ, Lee JM, Riley GR, Jang W, Rubinstein WS, et al. (2014) ClinVar: public archive of relationships among sequence variation and human phenotype. Nucleic Acids Res 42: D980-985.

23. Fokkema IF, Taschner PE, Schaafsma GC, Celli J, Laros JF, et al. (2011) LOVD v.2.0: the next generation in gene variant databases. Hum Mutat 32: 557-563.

24. Accetturo M, Creanza TM, Santoro C, Tria G, Giordano A, et al. (2010) Finding new genes for non-syndromic hearing loss through an in silico prioritization study. PLoS One 5:

25. Chang S, Zhang W, Gao L, Wang J (2012) Prioritization of candidate genes for attention deficit hyperactivity disorder by computational analysis of multiple data sources. Protein Cell 3: 526-534.

26. Kaplan JC. GeneTable of Neuromuscular Disorders. Available: http://musclegenetable.fr/. Accessed 2014 January 17.

27. Melia MJ, Kubota A, Ortolano S, Vilchez JJ, Gamez J, et al. (2013) Limb-girdle muscular dystrophy 1F is caused by a microdeletion in the transportin 3 gene. Brain 136: 1508-1517.

28. Tassy O, Pourquie O (2014) Manteia, a predictive data mining system for vertebrate genes and its applications to human genetic diseases. Nucleic Acids Res 42: D882-891.

29. Consortium GT (2013) The Genotype-Tissue Expression (GTEx) project. Nat Genet 45: 580-585.

30. Wicik Z, Sadkowski T, Jank M, Motyl T (2011) The transcriptomic signature of myostatin inhibitory influence on the differentiation of mouse C2C12 myoblasts. Pol J Vet Sci 14: 643-652.

31. Schwander M, Shirasaki R, Pfaff SL, Muller U (2004) Beta1 integrins in muscle, but not in motor neurons, are required for skeletal muscle innervation. J Neurosci 24: 8181-8191.

32. Gokhin DS, Fowler VM (2011) Cytoplasmic gamma-actin and tropomodulin isoforms link to the sarcoplasmic reticulum in skeletal muscle fibers. J Cell Biol 194: 105-120.

33. Burguiere AC, Nord H, von Hofsten J (2011) Alkali-like myosin light chain-1 (myl1) is an early marker for differentiating fast muscle cells in zebrafish. Dev Dyn 240: 1856-1863.

34. Mohaupt MG, Karas RH, Babiychuk EB, Sanchez-Freire V, Monastyrskaya K, et al. (2009) Association between statin-associated myopathy and skeletal muscle damage. CMAJ 181: E11-18.

35. Gregg RG, Messing A, Strube C, Beurg M, Moss R, et al. (1996) Absence of the beta subunit (cchb1) of the skeletal muscle dihydropyridine receptor alters expression of the alpha 1 subunit and eliminates excitation-contraction coupling. Proc Natl Acad Sci U S A 93: 13961-13966.

36. Carss KJ, Stevens E, Foley AR, Cirak S, Riemersma M, et al. (2013) Mutations in GDP-mannose pyrophosphorylase B cause congenital and limb-girdle muscular dystrophies associated with hypoglycosylation of alpha-dystroglycan. Am J Hum Genet 93: 29-41.

37. Buyse K, Riemersma M, Powell G, van Reeuwijk J, Chitayat D, et al. (2013) Missense mutations in beta-1,3-N-acetylglucosaminyltransferase 1 (B3GNT1) cause Walker-Warburg syndrome. Hum Mol Genet 22: 1746-1754.

38. Pacholsky D, Vakeel P, Himmel M, Lowe T, Stradal T, et al. (2004) Xin repeats define a novel actin-binding motif. J Cell Sci 117: 5257-5268.

39. Nilsson MI, Nissar AA, Al-Sajee D, Tarnopolsky MA, Parise G, et al. (2013) Xin is a marker of skeletal muscle damage severity in myopathies. Am J Pathol 183: 1703-1709.

40. Randazzo D, Giacomello E, Lorenzini S, Rossi D, Pierantozzi E, et al. (2013) Obscurin is required for ankyrinB-dependent dystrophin localization and sarcolemma integrity. J Cell Biol 200: 523-536.

41. Robinson RL, Curran JL, Ellis FR, Halsall PJ, Hall WJ, et al. (2000) Multiple interacting gene products may influence susceptibility to malignant hyperthermia. Ann Hum Genet 64: 307-320.

42. Iles DE, Segers B, Sengers RC, Monsieurs K, Heytens L, et al. (1993) Genetic mapping of the beta 1- and gamma-subunits of the human skeletal muscle L-type voltage-dependent calcium channel on chromosome 17q and exclusion as candidate genes for malignant hyperthermia susceptibility. Hum Mol Genet 2: 863-868.

43. Monnier N, Procaccio V, Stieglitz P, Lunardi J (1997) Malignant-hyperthermia susceptibility is associated with a mutation of the alpha 1-subunit of the human dihydropyridine-sensitive L-type voltage-dependent calcium-channel receptor in skeletal muscle. Am J Hum Genet 60: 1316-1325.

44. Zhang M, Wei J, Li Y, Shan H, Yan R, et al. (2014) Novel distribution of calreticulin to cardiomyocyte mitochondria and its increase in a rat model of dilated cardiomyopathy. Biochem Biophys Res Commun 449: 62-68.

45. Li H, Chen Q, Liu F, Zhang X, Li W, et al. (2013) Unfolded protein response and activated degradative pathways regulation in GNE myopathy. PLoS One 8: e58116.

46. Kumamoto T, Ito T, Horinouchi H, Ueyama H, Toyoshima I, et al. (2000) Increased lysosome-related proteins in the skeletal muscles of distal myopathy with rimmed vacuoles. Muscle Nerve 23: 1686-1693.

47. Blandin G, Marchand S, Charton K, Daniele N, Gicquel E, et al. (2013) A human skeletal muscle interactome centered on proteins involved in muscular dystrophies: LGMD interactome. Skelet Muscle 3: 3.

48. Bitoun M, Maugenre S, Jeannet PY, Lacene E, Ferrer X, et al. (2005) Mutations in dynamin 2 cause dominant centronuclear myopathy. Nat Genet 37: 1207-1209.

49. Zuchner S, Noureddine M, Kennerson M, Verhoeven K, Claeys K, et al. (2005) Mutations in the pleckstrin homology domain of dynamin 2 cause dominant intermediate Charcot-Marie-Tooth disease. Nat Genet 37: 289-294.

50. De Sandre-Giovannoli A, Chaouch M, Kozlov S, Vallat JM, Tazir M, et al. (2002) Homozygous defects in LMNA, encoding lamin A/C nuclear-envelope proteins, cause autosomal recessive axonal neuropathy in human (Charcot-Marie-Tooth disorder type 2) and mouse. Am J Hum Genet 70: 726-736.

51. De Sandre-Giovannoli A, Bernard R, Cau P, Navarro C, Amiel J, et al. (2003) Lamin A truncation in Hutchinson-Gilford progeria. Science 300: 2055.

52. Gros-Louis F, Dupre N, Dion P, Fox MA, Laurent S, et al. (2007) Mutations in SYNE1 lead to a newly discovered form of autosomal recessive cerebellar ataxia. Nat Genet 39: 80-85.

53. Zhang Q, Bethmann C, Worth NF, Davies JD, Wasner C, et al. (2007) Nesprin-1 and -2 are involved in the pathogenesis of Emery Dreifuss muscular dystrophy and are critical for nuclear envelope integrity. Hum Mol Genet 16: 2816-2833.

54. Attali R, Warwar N, Israel A, Gurt I, McNally E, et al. (2009) Mutation of SYNE-1, encoding an essential component of the nuclear lamina, is responsible for autosomal recessive arthrogryposis. Hum Mol Genet 18: 3462-3469.

55. Puckelwartz MJ, Kessler EJ, Kim G, Dewitt MM, Zhang Y, et al. (2010) Nesprin-1 mutations in human and murine cardiomyopathy. J Mol Cell Cardiol 48: 600-608.
56. Bohm J, Biancalana V, Dechene ET, Bitoun M, Pierson CR, et al. (2012) Mutation spectrum in the large GTPase dynamin 2, and genotype-phenotype correlation in autosomal dominant centronuclear myopathy. Hum Mutat 33: 949-959.

57. Cortese A, Tucci A, Piccolo G, Galimberti CA, Fratta P, et al. (2014) Novel CLN3 mutation causing autophagic vacuolar myopathy. Neurology.

58. Ohkawara B, Cabrera-Serrano M, Nakata T, Milone M, Asai N, et al. (2014) LRP4 third beta-propeller domain mutations cause novel congenital myasthenia by compromising agrin-mediated MuSK signaling in a position-specific manner. Hum Mol Genet 23: 1856-1868.

59. Fernandez-Valle C, Tang Y, Ricard J, Rodenas-Ruano A, Taylor A, et al. (2002) Paxillin binds schwannomin and regulates its density-dependent localization and effect on cell morphology. Nat Genet 31: 354-362.

60. Vaughan KT, Weber FE, Ried T, Ward DC, Reinach FC, et al. (1993) Human myosin-binding protein H (MyBP-H): complete primary sequence, genomic organization, and chromosomal localization. Genomics 16: 34-40.

61. Sandri M, Lin J, Handschin C, Yang W, Arany ZP, et al. (2006) PGC-1alpha protects skeletal muscle from atrophy by suppressing FoxO3 action and atrophy-specific gene transcription. Proc Natl Acad Sci U S A 103: 16260-16265.

62. Mundel P, Heid HW, Mundel TM, Kruger M, Reiser J, et al. (1997) Synaptopodin: an actin-associated protein in telencephalic dendrites and renal podocytes. J Cell Biol 139: 193-204.