Inferring the effect of genomic variation in the new era of genomics

Utilization of Targeted RNA-Seq for the Resolution of Variant Pathogenicity and Enhancement of Diagnostic Yield in Dysferlinopathy

For inherited diseases, obtaining a definitive diagnosis is critical for proper disease management, family planning, and participation in clinical trials. This can be challenging for dysferlinopathy due to the significant clinical overlap between the 30+ subtypes of limb–girdle muscular dystrophy (LGMD) and the large number of variants of unknown significance (VUSs) that are identified in the dysferlin gene, DYSF. We performed targeted RNA-Seq using a custom gene-panel in 77 individuals with a clinical/genetic suspicion of dysferlinopathy and evaluated all 111 identified DYSF variants according to the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) guidelines. This evaluation identified 11 novel DYSF variants and allowed for the classification of 87 DYSF variants as pathogenic/likely pathogenic, 8 likely benign, while 16 variants remained VUSs. By the end of the study, 60 of the 77 cases had a definitive diagnosis of dysferlinopathy, which was a 47% increase in diagnostic yield over the rate at study onset. This data shows the ability of RNA-Seq to assist in variant pathogenicity classification and diagnosis of dysferlinopathy and is, therefore, a type of analysis that should be considered when DNA-based genetic analysis is not sufficient to provide a definitive diagnosis.

Identification of new variants and candidate genes in women with familial premature ovarian insufficiency using whole-exome sequencing

PURPOSE

To identify candidate variants in genes possibly associated with premature ovarian insufficiency (POI).

METHODS

Fourteen women, from 7 families, affected by idiopathic POI were included. Additionally, 98 oocyte donors of the same ethnicity were enrolled as a control group. Whole-exome sequencing (WES) was performed in 14 women with POI to identify possibly pathogenic variants in genes potentially associated with the ovarian function. The candidate genes selected in POI patients were analysed within the exome results of oocyte donors.

RESULTS

After the variant filtering in the WES analysis of 7 POI families, 23 possibly damaging genetic variants were identified in 22 genes related to POI or linked to ovarian physiology. All variants were heterozygous and five of the seven families carried two or more variants in different genes. We have described genes that have never been associated to POI pathology; however, they are involved in important biological processes for ovarian function. In the 98 oocyte donors of the control group, we found no potentially pathogenic variants among the 22 candidate genes.

CONCLUSION

WES has previously shown as an efficient tool to identify causative genes for ovarian failure. Although some studies have focused on it, and many genes are identified, this study proposes new candidate genes and variants, having potentially moderate/strong functional effects, associated with POI, and argues for a polygenic etiology of POI in some cases.

Genome interpretation using in silico predictors of variant impact

Estimating the effects of variants found in disease driver genes opens the door to personalized therapeutic opportunities. Clinical associations and laboratory experiments can only characterize a tiny fraction of all the available variants, leaving the majority as variants of unknown significance (VUS). In silico methods bridge this gap by providing instant estimates on a large scale, most often based on the numerous genetic differences between species. Despite concerns that these methods may lack reliability in individual subjects, their numerous practical applications over cohorts suggest they are already helpful and have a role to play in genome interpretation when used at the proper scale and context. In this review, we aim to gain insights into the training and validation of these variant effect predicting methods and illustrate representative types of experimental and clinical applications. Objective performance assessments using various datasets that are not yet published indicate the strengths and limitations of each method. These show that cautious use of in silico variant impact predictors is essential for addressing genome interpretation challenges.

Peptide Model of the Mutant Proinsulin Syndrome. II. Nascent Structure and Biological Implications

Toxic misfolding of proinsulin variants in β-cells defines a monogenic diabetes syndrome, designated mutant INS-gene induced diabetes of the young (MIDY). In our first study (previous article in this issue), we described a one-disulfide peptide model of a proinsulin folding intermediate and its use to study such variants. The mutations (LeuB15→Pro, LeuA16→Pro, and PheB24→Ser) probe residues conserved among vertebrate insulins. In this companion study, we describe 1H and 1H-13C NMR studies of the peptides; key NMR resonance assignments were verified by synthetic 13C-labeling. Parent spectra retain nativelike features in the neighborhood of the single disulfide bridge (cystine B19-A20), including secondary NMR chemical shifts and nonlocal nuclear Overhauser effects. This partial fold engages wild-type side chains LeuB15, LeuA16 and PheB24 at the nexus of nativelike α-helices α1 and α3 (as defined in native proinsulin) and flanking β-strand (residues B24-B26). The variant peptides exhibit successive structural perturbations in order: parent (most organized) > SerB24 >> ProA16 > ProB15 (least organized). The same order pertains to (a) overall α-helix content as probed by circular dichroism, (b) synthetic yields of corresponding three-disulfide insulin analogs, and (c) ER stress induced in cell culture by corresponding mutant proinsulins. These findings suggest that this and related peptide models will provide a general platform for classification of MIDY mutations based on molecular mechanisms by which nascent disulfide pairing is impaired. We propose that the syndrome's variable phenotypic spectrum-onsets ranging from the neonatal period to later in childhood or adolescence-reflects structural features of respective folding intermediates.

Spinal health in 40 patients in the initial stage of laboural life. Morfogeometric, biological and environmental study

Introduction

Assuming that spinal shape is a genetic expression, its analysis and acquired factors could assess their respective contribution to early spine deterioration.

Material and methods

A geometric morphometric analysis was retrospectively performed on sagittal lumbar MRI of young patients with back pain to identify lumbar spine shape changes. Using Geometric Morphometrics, findings were analyzed with anthropometric, radiological, and clinical variables.

Results

80 cases under 26 years of age were collected, 55 men (mean age 22.81) and 25 women (mean age 23.24). MRI abnormalities were reported in 57.5%: single altered disc (N=17), root compromises (N=8), and transition anomalies (35%).In the non-normal MRI subgroup, shape variation included: increased lordosis, enlarged vertebral body, canal stenosis, and lumbarization of S1. In non-Spanish origin patients, lumbar straightening and segmental deformities were prevalent. Morphometrics findings showed that lumbosacral transition anomalies are frequently underreported.

Conclusions

Genetic factors could be the main determinants of abnormality in MRIs under 26 years. The primary markers are transitional abnormalities, segmental deformities, and canal stenosis. In foreign populations, shape changes could suggest spine overload at an early age.

Challenges in variant interpretation in prenatal exome sequencing

The use of exome sequencing (ES) in the prenatal setting improves the diagnostic yield of genetic testing for fetuses with ultrasound anomalies. However, while the purpose of ES is to explain the fetal phenotype, secondary or incidental findings unrelated to the observed abnormalities might be detected. Recently, requests for ES in fetuses with no sonographic abnormalities have been increasing, raising serious ethical and medico-legal concerns. Variant interpretation is complex even in the postnatal setting and performing broad genomic data analyses in the prenatal setting presents additional dilemmas. This article discusses challenges and questions related to prenatal ES, including variant interpretation of incidental findings in cases of indicated prenatal ES, as well as in situations where ES is performed in asymptomatic fetuses.

Laniakea@ReCaS: exploring the potential of customisable Galaxy on-demand instances as a cloud-based service

BACKGROUND

Improving the availability and usability of data and analytical tools is a critical precondition for further advancing modern biological and biomedical research. For instance, one of the many ramifications of the COVID-19 global pandemic has been to make even more evident the importance of having bioinformatics tools and data readily actionable by researchers through convenient access points and supported by adequate IT infrastructures. One of the most successful efforts in improving the availability and usability of bioinformatics tools and data is represented by the Galaxy workflow manager and its thriving community. In 2020 we introduced Laniakea, a software platform conceived to streamline the configuration and deployment of "on-demand" Galaxy instances over the cloud. By facilitating the set-up and configuration of Galaxy web servers, Laniakea provides researchers with a powerful and highly customisable platform for executing complex bioinformatics analyses. The system can be accessed through a dedicated and user-friendly web interface that allows the Galaxy web server's initial configuration and deployment.

RESULTS

"Laniakea@ReCaS", the first instance of a Laniakea-based service, is managed by ELIXIR-IT and was officially launched in February 2020, after about one year of development and testing that involved several users. Researchers can request access to Laniakea@ReCaS through an open-ended call for use-cases. Ten project proposals have been accepted since then, totalling 18 Galaxy on-demand virtual servers that employ ~ 100 CPUs, ~ 250 GB of RAM and ~ 5 TB of storage and serve several different communities and purposes. Herein, we present eight use cases demonstrating the versatility of the platform.

CONCLUSIONS

During this first year of activity, the Laniakea-based service emerged as a flexible platform that facilitated the rapid development of bioinformatics tools, the efficient delivery of training activities, and the provision of public bioinformatics services in different settings, including food safety and clinical research. Laniakea@ReCaS provides a proof of concept of how enabling access to appropriate, reliable IT resources and ready-to-use bioinformatics tools can considerably streamline researchers' work.

Chemical map-based prediction of nucleosome positioning using the Bioconductor package nuCpos

Background

Assessing the nucleosome-forming potential of specific DNA sequences is important for understanding complex chromatin organization. Methods for predicting nucleosome positioning include bioinformatics and biophysical approaches. An advantage of bioinformatics methods, which are based on in vivo nucleosome maps, is the use of natural sequences that may contain previously unknown elements involved in nucleosome positioning in vivo. The accuracy of such prediction attempts reflects the genomic coordinate resolution of the nucleosome maps applied. Nucleosome maps are constructed using micrococcal nuclease digestion followed by high-throughput sequencing (MNase-seq). However, as MNase has a strong preference for A/T-rich sequences, MNase-seq may not be appropriate for this purpose. In addition to MNase-seq-based maps, base pair-resolution chemical maps of in vivo nucleosomes from three different species (budding and fission yeasts, and mice) are currently available. However, these chemical maps have yet to be integrated into publicly available computational methods.

Results

We developed a Bioconductor package (named nuCpos) to demonstrate the superiority of chemical maps in predicting nucleosome positioning. The accuracy of chemical map-based prediction in rotational settings was higher than that of the previously developed MNase-seq-based approach. With our method, predicted nucleosome occupancy reasonably matched in vivo observations and was not affected by A/T nucleotide frequency. Effects of genetic alterations on nucleosome positioning that had been observed in living yeast cells could also be predicted. nuCpos calculates individual histone binding affinity (HBA) scores for given 147-bp sequences to examine their suitability for nucleosome formation. We also established local HBA as a new parameter to predict nucleosome formation, which was calculated for 13 overlapping nucleosomal DNA subsequences. HBA and local HBA scores for various sequences agreed well with previous in vitro and in vivo studies. Furthermore, our results suggest that nucleosomal subsegments that are disfavored in different rotational settings contribute to the defined positioning of nucleosomes.

Conclusions

Our results demonstrate that chemical map-based statistical models are beneficial for studying nucleosomal DNA features. Studies employing nuCpos software can enhance understanding of chromatin regulation and the interpretation of genetic alterations and facilitate the design of artificial sequences.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04240-2.

Inborn errors of immunity-recent advances in research on the pathogenesis

Primary immunodeficiency (PID) is a genetic disorder with a defect of one of the important components of our immune system. Classical PID has been recognized as a disorder with loss of function of the immune system. Recent studies have unveiled disorders with immune dysfunction with autoimmunity, autoinflammation, allergy, or predisposition to malignancy. Some of them were caused by an augmented immune function or a defect in immune regulation. With this background, the term inborn errors of immunity (IEI) is now used to refer to PID in the International Union of Immunological Societies (IUIS) classification. More than 400 responsible genes have been identified in patients with IEI so far, and importantly, many of them identified lately were caused by a heterologous mutation. Moreover, the onset is not necessarily in childhood, and we started seeing more and more IEI patients diagnosed in adulthood in the clinical settings. Recent advances in genetic analysis, including whole-exome analysis, whole-genome analysis, and RNA-seq have contributed to the identification of the disease-causing gene mutation. We also started to find heterogeneity of phenotype even in the patients with the same mutation in the same family, leading us to wonder if modifier gene or epigenetic modification is involved in the pathogenesis. In contrast, we accumulated many cases suggesting genetic heterogeneity is associated with phenotypic homogeneity. It has thus become difficult to deduce a responsible gene only from the phenotype in a certain type of IEI. Current curative therapy for IEI includes hematopoietic cell transplantation and gene therapy. Other curative therapeutic modalities have been long waited and are to be introduced in the future. These include a small molecule that inhibits the gain-of-function of the molecule- and genome-editing technology. Research on IEI will surely lead to a better understanding of other immune-related disorders including rheumatic diseases and atopic disorders.

Clinical and Genomic Evaluation of 207 Genetic Myopathies in the Indian Subcontinent

Objective: Inherited myopathies comprise more than 200 different individually rare disease-subtypes, but when combined together they have a high prevalence of 1 in 6,000 individuals across the world. Our goal was to determine for the first time the clinical- and gene-variant spectrum of genetic myopathies in a substantial cohort study of the Indian subcontinent.

Methods: In this cohort study, we performed the first large clinical exome sequencing (ES) study with phenotype correlation on 207 clinically well-characterized inherited myopathy-suspected patients from the Indian subcontinent with diverse ethnicities.

Results: Clinical-correlation driven definitive molecular diagnosis was established in 49% (101 cases; 95% CI, 42–56%) of patients with the major contributing pathogenicity in either of three genes, GNE (28%; GNE-myopathy), DYSF (25%; Dysferlinopathy), and CAPN3 (19%; Calpainopathy). We identified 65 variant alleles comprising 37 unique variants in these three major genes. Seventy-eight percent of the DYSF patients were homozygous for the detected pathogenic variant, suggesting the need for carrier-testing for autosomal-recessive disorders like Dysferlinopathy that are common in India. We describe the observed clinical spectrum of myopathies including uncommon and rare subtypes in India: Sarcoglycanopathies (SGCA/B/D/G), Collagenopathy (COL6A1/2/3), Anoctaminopathy (ANO5), telethoninopathy (TCAP), Pompe-disease (GAA), Myoadenylate-deaminase-deficiency-myopathy (AMPD1), myotilinopathy (MYOT), laminopathy (LMNA), HSP40-proteinopathy (DNAJB6), Emery-Dreifuss-muscular-dystrophy (EMD), Filaminopathy (FLNC), TRIM32-proteinopathy (TRIM32), POMT1-proteinopathy (POMT1), and Merosin-deficiency-congenital-muscular-dystrophy-type-1 (LAMA2). Thirteen patients harbored pathogenic variants in >1 gene and had unusual clinical features suggesting a possible role of synergistic-heterozygosity/digenic-contribution to disease presentation and progression.

Conclusions: Application of clinically correlated ES to myopathy diagnosis has improved our understanding of the clinical and genetic spectrum of different subtypes and their overlaps in Indian patients. This, in turn, will enhance the global gene-variant-disease databases by including data from developing countries/continents for more efficient clinically driven molecular diagnostics.

Whole exome sequencing identifies multiple novel candidate genes in familial gastroschisis

Background

Genetic association studies for gastroschisis have highlighted several candidate variants. However, genetic basis in gastroschisis from noninvestigated heritable factors could provide new insights into the human biology for this birth defect. We aim to identify novel gastroschisis susceptibility variants by employing whole exome sequencing (WES) in a Mexican family with recurrence of gastroschisis.

Methods

We employed WES in two affected half‐sisters with gastroschisis, mother, and father of the proband. Additionally, functional bioinformatics analysis was based on SVS–PhoRank and Ensembl–Variant Effect Predictor. The latter assessed the potentially deleterious effects (high, moderate, low, or modifier impact) from exome variants based on SIFT, PolyPhen, dbNSFP, Condel, LoFtool, MaxEntScan, and BLOSUM62 algorithms. The analysis was based on the Human Genome annotation, GRCh37/hg19. Candidate genes were prioritized and manually curated based on significant phenotypic relevance (SVS–PhoRank) and functional properties (Ensembl–Variant Effect Predictor). Functional enrichment analysis was performed using ToppGene Suite, including a manual curation of significant Gene Ontology (GO) biological processes from functional similarity analysis of candidate genes.

Results

No single gene‐disrupting variant was identified. Instead, 428 heterozygous variations were identified for which SPATA17, PDE4DIP, CFAP65, ALPP, ZNF717, OR4C3, MAP2K3, TLR8, and UBE2NL were predicted as high impact in both cases, mother, and father of the proband. PLOD1, COL6A3, FGFRL1, HHIP, SGCD, RAPGEF1, PKD1, ZFHX3, BCAS3, EVPL, CEACAM5, and KLK14 were segregated among both cases and mother. Multiple interacting background modifiers may regulate gastroschisis susceptibility. These candidate genes highlight a role for development of blood vessel, circulatory system, muscle structure, epithelium, and epidermis, regulation of cell junction assembly, biological/cell adhesion, detection/response to endogenous stimulus, regulation of cytokine biosynthetic process, response to growth factor, postreplication repair/protein K63‐linked ubiquitination, protein‐containing complex assembly, and regulation of transcription DNA‐templated.

Conclusion

Considering the likely gene‐disrupting prediction results and similar biological pattern of mechanisms, we propose a joint “multifactorial model” in gastroschisis pathogenesis.

[Rational use of genetic tests in internal medicine : Possibilities and limitations of next generation sequencing diagnostics]

BACKGROUND

New methods of molecular genetic diagnostics enable a more comprehensive genetic analysis of patients.

OBJECTIVES

Rational use and benefits of molecular genetic testing in patients with various internal diseases.

METHOD

Evaluation of topic-related literature, discussion of own experiences, as well as consideration of current guidelines.

RESULTS

New genetic tests, such as next generation sequencing (NGS), improve the diagnosis of hereditary diseases; however, the use of this technology also leads to additional findings, which must be carefully considered.

CONCLUSION

The rational use of genetic tests is a benefit for patients and can significantly influence the prevention and treatment of a disease. The increasing complexity of genetic findings requires interdisciplinary approaches involving human genetics, internal medicine, and other disciplines.

NQO1: A target for the treatment of cancer and neurological diseases, and a model to understand loss of function disease mechanisms

NAD(P)H quinone oxidoreductase 1 (NQO1) is a multi-functional protein that catalyses the reduction of quinones (and other molecules), thus playing roles in xenobiotic detoxification and redox balance, and also has roles in stabilising apoptosis regulators such as p53. The structure and enzymology of NQO1 is well-characterised, showing a substituted enzyme mechanism in which NAD(P)H binds first and reduces an FAD cofactor in the active site, assisted by a charge relay system involving Tyr-155 and His-161. Protein dynamics play important role in physio-pathological aspects of this protein. NQO1 is a good target to treat cancer due to its overexpression in cancer cells. A polymorphic form of NQO1 (p.P187S) is associated with increased cancer risk and certain neurological disorders (such as multiple sclerosis and Alzheimer´s disease), possibly due to its roles in the antioxidant defence. p.P187S has greatly reduced FAD affinity and stability, due to destabilization of the flavin binding site and the C-terminal domain, which leading to reduced activity and enhanced degradation. Suppressor mutations partially restore the activity of p.P187S by local stabilization of these regions, and showing long-range allosteric communication within the protein. Consequently, the correction of NQO1 misfolding by pharmacological chaperones is a viable strategy, which may be useful to treat cancer and some neurological conditions, targeting structural spots linked to specific disease-mechanisms. Thus, NQO1 emerges as a good model to investigate loss of function mechanisms in genetic diseases as well as to improve strategies to discriminate between neutral and pathogenic variants in genome-wide sequencing studies.

Bioinformatic Analysis of Gene Variants from Gastroschisis Recurrence Identifies Multiple Novel Pathogenetic Pathways: Implication for the Closure of the Ventral Body Wall

We investigated whether likely pathogenic variants co-segregating with gastroschisis through a family-based approach using bioinformatic analyses were implicated in body wall closure. Gene Ontology (GO)/Panther functional enrichment and protein-protein interaction analysis by String identified several biological networks of highly connected genes in UGT1A3, UGT1A4, UGT1A5, UGT1A6, UGT1A7, UGT1A8, UGT1A9, UGT1A10, AOX1, NOTCH1, HIST1H2BB, RPS3, THBS1, ADCY9, and FGFR4. SVS–PhoRank identified a dominant model in OR10G4 (also as heterozygous de novo), ITIH3, PLEKHG4B, SLC9A3, ITGA2, AOX1, and ALPP, including a recessive model in UGT1A7, UGT1A6, PER2, PTPRD, and UGT1A3. A heterozygous compound model was observed in CDYL, KDM5A, RASGRP1, MYBPC2, PDE4DIP, F5, OBSCN, and UGT1A. These genes were implicated in pathogenetic pathways involving the following GO related categories: xenobiotic, regulation of metabolic process, regulation of cell adhesion, regulation of gene expression, inflammatory response, regulation of vascular development, keratinization, left-right symmetry, epigenetic, ubiquitination, and regulation of protein synthesis. Multiple background modifiers interacting with disease-relevant pathways may regulate gastroschisis susceptibility. Based in our findings and considering the plausibility of the biological pattern of mechanisms and gene network modeling, we suggest that the gastroschisis developmental process may be the consequence of several well-orchestrated biological and molecular mechanisms which could be interacting with gastroschisis predispositions within the first ten weeks of development.

Genetic landscape and novel disease mechanisms from a large LGMD cohort of 4656 patients

Objective

Limb‐girdle muscular dystrophies (LGMDs), one of the most heterogeneous neuromuscular disorders (NMDs), involves predominantly proximal‐muscle weakness with >30 genes associated with different subtypes. The clinical‐genetic overlap among subtypes and with other NMDs complicate disease‐subtype identification lengthening diagnostic process, increases overall costs hindering treatment/clinical‐trial recruitment. Currently seven LGMD clinical trials are active but still no gene‐therapy‐related treatment is available. Till‐date no nation‐wide large‐scale LGMD sequencing program was performed. Our objectives were to understand LGMD genetic basis, different subtypes’ relative prevalence across US and investigate underlying disease mechanisms.

Methods

A total of 4656 patients with clinically suspected‐LGMD across US were recruited to conduct next‐generation sequencing (NGS)‐based gene‐panel testing during June‐2015 to June‐2017 in CLIA‐CAP‐certified Emory‐Genetics‐Laboratory. Thirty‐five LGMD‐subtypes‐associated or LGMD‐like other NMD‐associated genes were investigated. Main outcomes were diagnostic yield, gene‐variant spectrum, and LGMD subtypes’ prevalence in a large US LGMD‐suspected population.

Results

Molecular diagnosis was established in 27% (1259 cases; 95% CI, 26–29%) of the patients with major contributing genes to LGMD phenotypes being: CAPN3(17%), DYSF(16%), FKRP(9%) and ANO5(7%). We observed an increased prevalence of genetically confirmed late‐onset Pompe disease, DNAJB6‐associated LGMD subtype1E and CAPN3‐associated autosomal‐dominant LGMDs. Interestingly, we identified a high prevalence of patients with pathogenic variants in more than one LGMD gene suggesting possible synergistic heterozygosity/digenic/multigenic contribution to disease presentation/progression that needs consideration as a part of diagnostic modality.

Interpretation

Overall, this study has improved our understanding of the relative prevalence of different LGMD subtypes, their respective genetic etiology, and the changing paradigm of their inheritance modes and novel mechanisms that will allow for improved timely treatment, management, and enrolment of molecularly diagnosed individuals in clinical trials.

Insight into the specificity and severity of pathogenic mechanisms associated with missense mutations through experimental and structural perturbation analyses

Most pathogenic missense mutations cause specific molecular phenotypes through protein destabilization. However, how protein destabilization is manifested as a given molecular phenotype is not well understood. We develop here a structural and energetic approach to describe mutational effects on specific traits such as function, regulation, stability, subcellular targeting or aggregation propensity. This approach is tested using large-scale experimental and structural perturbation analyses in over thirty mutations in three different proteins (cancer-associated NQO1, transthyretin related with amyloidosis and AGT linked to primary hyperoxaluria type I) and comprising five very common pathogenic mechanisms (loss-of-function and gain-of-toxic function aggregation, enzyme inactivation, protein mistargeting and accelerated degradation). Our results revealed that the magnitude of destabilizing effects and, particularly, their propagation through the structure to promote disease-associated conformational states largely determine the severity and molecular mechanisms of disease-associated missense mutations. Modulation of the structural perturbation at a mutated site is also shown to cause switches between different molecular phenotypes. When very common disease-associated missense mutations were investigated, we also found that they were not among the most deleterious possible missense mutations at those sites, and required additional contributions from codon bias and effects of CpG sites to explain their high frequency in patients. Our work sheds light on the molecular basis of pathogenic mechanisms and genotype–phenotype relationships, with implications for discriminating between pathogenic and neutral changes within human genome variability from whole genome sequencing studies.

Myofibrillar myopathy in the genomic context

Myofibrillar myopathy (MFM) is a group of inherited muscular disorders characterized by myofibril dissolution and abnormal accumulation of degradation products. The diagnosis of muscular disorders based on clinical presentation is difficult due to phenotypic heterogeneity and overlapping symptoms. In addition, precise diagnosis does not always explain the disease etiopathology or the highly variable clinical course even among patients diagnosed with the same type of myopathy. The advent of high-throughput next-generation sequencing (NGS) has provided a successful and cost-effective strategy for identification of novel causative genes in myopathies, including MFM. So far, pathogenic mutations associated with MFM phenotype, including atypical MFM-like cases, have been identified in 17 genes: DES, CRYAB, MYOT, ZASP, FLNC, BAG3, FHL1, TTN, DNAJB6, PLEC, LMNA, ACTA1, HSPB8, KY, PYROXD1, and SQSTM + TIA1 (digenic). Most of these genes are also associated with other forms of muscle diseases. In addition, in many MFM patients, numerous genomic variants in muscle-related genes have been identified. The various myopathies and muscular dystrophies seem to form a single disease continuum; therefore, gene identification in one disease impacts the genetic etiology of the others. In this review, we describe the heterogeneity of the MFM genetic background focusing on the role of rare variants, the importance of whole genome sequencing in the identification of novel disease-associated mutations, and the emerging concept of variant load as the basis of the phenotypic heterogeneity.