Next-generation sequencing and bioinformatics to identify genetic causes of malignant hyperthermia

Background/purpose: Malignant hyperthermia (MH) is a life-threatening pharmacogenetic disease with only two known causative genes, RYR1 and CACNA1S. Both are huge genes containing numerous exons, and they reportedly only account for 50-70% of known MH patients. Next-generation sequencing (NGS) technology and bioinformatics could help delineate the genetic diagnosis of MH and several MH-like clinical presentations. Methods: We established a capture-based targeted NGS sequencing framework to examine the whole genomic regions of RYR1, CACNA1S and the 16.6 Kb mitochondrial genome, as well as 12 other genes related to excitation-contraction coupling and/or skeletal muscle calcium homeostasis. We applied bioinformatics analyses to the variants identified in this study and also to the 48 documented RYR1 pathogenic variants. Results: The causative variants were identified in seven of the eight (87.5%) MH families, but in none of the 10 individuals classified as either normal controls (N = 2) or patients displaying MH-like clinical features later found to be caused by other etiologies (N = 8). We showed that RYR1 c.1565A>G (p.Tyr522Cys)(rs118192162) could be a genetic hot spot in the Taiwanese population. Bioinformatics analyses demonstrated low population frequencies and predicted damaging effects from all known pathogenic RYR1 variants. We estimated that more than one in 1149 individuals worldwide carry MH pathogenic variants at RYR1. Background/purpose: Malignant hyperthermia (MH) is a life-threatening pharmacogenetic disease with only two known causative genes, RYR1 and CACNA1S. Both are huge genes containing numerous exons, and they reportedly only account for 50-70% of known MH patients. Next generation sequencing (NGS) technology and bioinformatics could help delineate the genetic diagnosis of MH and several MH-like clinical presentations. Methods: We established a capture-based targeted NGS sequencing framework to examine the whole genomic regions of RYR1, CACNA1S and the 16.6 Kb mitochondrial genome, as well as 12 other genes related to excitation-contraction coupling and/or skeletal muscle calcium homeostasis. We applied bioinformatics analyses to the variants identified in this study and also to the 48 documented RYR1 pathogenic variants. Results: The causative variants were identified in seven of the eight (87.5%) MH families, but in none of the 10 individuals classified as either normal controls (N = 2) or patients displaying MH-like clinical features later found to be caused by other etiologies (N = 8). We showed that RYR1 c.1565A>G (p.Tyr522Cys)(rs118192162) could be a genetic hot spot in the Taiwanese population. Bioinformatics analyses demonstrated low population frequencies and predicted damaging effects from all known pathogenic RYR1 variants. We estimated that more than one in 1149 individuals worldwide carry MH pathogenic variants at RYR1. Conclusion: NGS and bioinformatics are sensitive and specific tools to examine RYR1 and CACNA1S for the genetic diagnosis of MH. Pathogenic variants in RYR1 can be found in the majority of MH patients in Taiwan. Copyright (C) 2020, Formosan Medical Association. Published by Elsevier Taiwan LLC. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).