Prenatal diagnosis identifies compound heterozygous variants in RYR1 that causes ultrasound abnormalities in a fetus
BMC Medical Genomics volume 15, Article number: 202 (2022)
We presented a non-consanguineous healthy Chinese couple with five pregnancies, three early miscarriages, the fetus II-2 and II-5 with similar abnormal phenotypes of fetal hydrops, scoliosis, fetal akinesia and polyhydramnios. This study aimed to uncover the molecular etiology of this family with a history of multiple adverse pregnancies.
Materials and methods
DNA extracted from the fifth fetal umbilical cord and parents’ peripheral blood were subjected to SNP-array and whole exome sequencing. The result was verified by Sanger sequencing. Functional characterization of the c.2682G > C (p.Ile860_Pro894del) variant was completed by minigene splicing assay.
Trio whole-exome sequencing has identified compound heterozygous variants in RYR1 (c.2682G > C; p.Ile860_Pro894del and c.12572G > A; p.Arg4191His) in fetus II-5. The variant c.2682G > C (p.Ile860_Pro894del) comes from the father and the c.12572G > A (p.Arg4191His) comes from the mother. The c.2682G > C (p.Ile860_Pro894del) affects the splice site resulting in exon 21 skipping, therefore is classified as likely pathogenic. The c.12572G > A (p.Arg4191His) locates in the C-terminal hot spots region of the RYR1, classified as of uncertain significance.
We report the first prenatal case of RYR1-related disorders in Chinese population, expanding the variant spectrum of RYR1 in fetuses.
The skeletal muscle ryanodine receptor (RyR1) locates in the sarcoplasmic reticulum (SR) of skeletal muscle and releases calcium (Ca2+) into the cytoplasm upon depolarization of the sarcolemma, playing a significant role in excitation–contraction coupling . As one of the largest genes in the human genome, the RYR1 gene includes 106 exons and covers 15.4 kb on chromosome 19q13.1 .
Variant in the RYR1 gene is first described in a Canadian family with malignant hyperthermia (MH) . RYR1 variants have subsequently been identified in various myopathies including central core disease (CCD), multi-minicore disease (MmD), centronuclear myopathy (CNM), congenital fiber-type disproportion (CFTD), core-rod myopathy (CRM), King-Denborough syndrome, dusty core disease (DuCD), samaritan myopathy, a late-onset axial myopathy, fetal akinesia and lethal multiple pterygium syndrome (LMPS) [4,5,6,7,8,9]. The age of onset, clinical presentation, disease progression and severity of these myopathies vary widely. With the constant advancement of next generation sequencing, the RYR1 disease spectrum has further expanded. Clinical and pathological overlaps make the differential diagnosis of RYR1-related myopathies difficult. Lawal et al. propose RYR1-related disorders (RYR1-RD) as a unified nomenclature to describe this complex and evolving disease spectrum .
In this study, we report a Chinese couple with five pregnancies who had three early miscarriages and two fetal losses because of a primary myopathy characterized by fetal akinesia, polyhydramnios, umbilical cord cyst, abnormality of the vertebral column, abnormality of the scrotum, rocker bottom foot, absence of stomach bubble and hydrops fetalis. However, no further genetic tests were carried out on the first four fetuses as their samples were not available. Abnormal karyotype and pathogenic CNV of the fifth fetus and parents had been excluded. Trio whole-exome sequencing was performed to identify the genetic causes of these symptoms. Lethality of the condition and potential phenotype recurrence based on the reproductive history suggested an autosomal recessive condition.
Material and methods
We report a non-consanguineous Chinese family that had five pregnancies affected by RYR1-related disorders (Fig. 1A). The mother was a 27 year-old woman. The couple was healthy and had no personal or family history of the neuromuscular disorder or malignant hyperthermia. The first pregnancy (II-1) ended in abortion.
The second pregnancy (II-2) was complicated by fetal ultrasound findings that polyhydramnios, single umbilical artery, intrauterine growth retardation, and absence of stomach bubble at 28 weeks of gestation. At post mortem examination there was spina bifida and scoliosis associated with vertebral segmentation defects. An additional genetic analysis was not selected.
Their third and fourth pregnancy (II-3 and II-4) were terminated, because of no yolk sac in early pregnancy.
Subsequently, fetal ultrasound of the fifth pregnancy (II-5) revealed increased nuchal translucency (NT: 8.8 mm), hand clenching, adducted thumb, absence of stomach bubble on fetal sonography, polyhydramnios, umbilical cord cysts, abnormal conus medullaris morphology, and thickened heel subcutaneous tissue (Fig. 1B–D) at 25 weeks’ gestation. Magnetic resonance imaging (MRI) showed abnormality of the scrotum, thoracic scoliosis, abnormal conus medullaris morphology, hand clenching, rocker bottom foot, fetal immobility, and subcutaneous oedema (Fig. 1E–K). The fetal and the couple's karyotypes were normal. Based on the multiple malformation outcomes, the pregnancy was terminated 3 weeks later. Whole exome sequencing analysis was performed using the fetal umbilical cord and parents’ blood. Unfortunately, the II-5 muscle samples were not retained for histological studies.
All participants provided written informed consent. The study protocols were approved by the ethics committee from the Second Affiliated Hospital of Chongqing Medical University.
SNP-array and exome sequencing
Genomic DNA was extracted from the umbilical cord of II-5 and peripheral blood of parents, using DNeasy Blood & Tissue kit (Qiagen, Hilden, Germany). Single nucleotide polymorphism array analysis was performed on the Affymetrix CytoScan platform (Affymetrix, Santa Clara, CA, USA) following the protocol. Exome capture was performed on 3 µg of genomic DNA sample using the Agilent SureSelectXT V5 capture kit (Agilent Technologies, Santa Clara, CA), and then sequenced on Illumina HiSeq 2000 (Illumina, San Diego, CA) with 100-bp paired-end reads.
Variant annotation and selection
SNP-array data were analyzed by Chromosome Analysis Suite 3.2 (Affymetrix, Santa Clara, CA, USA). Raw data were mapped to the human reference genome (GRCh37/hg19) using Burrows‐Wheeler Aligner  and variants calling using Genomic Analysis Tool Kit best practices . Variants were annotated by ANNOVAR using ExAC, 1000 Genomes, Genome Aggregation Database (gnomAD) and other public databases . Stepwise filtering included the removal of variants with quality scores < 20 and allele counts < 10X, minor allele frequency (MAF) > 1%, common SNPs, intergenic and 3’/5’ UTR variants, non-splice-related intronic variants, benign variants and likely benign variants, and synonymous variants. Only variants matched the inheritance model of the family and fetuses’ phenotypes were kept. The pathogenicity of variants was interpreted and classified according to the American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) .
Parental blood samples as well as amniotic fluid of the fetus II-5 were used for sanger sequencing after WES analysis. Primers were designed using Primer3 online (http://primer3.ut.ee/). Primers RYR1_1 (Forward, 5’-AGAACATCCACGAGCTCTGG-3’ and Reverse, 5’-CCAGACTATGACCCCTGACC-3’) for variant c.2682G > C (p.Ile860_Pro894del) and RYR1_2 (Forward, 5’-ATCCTTGAGTACTTCCGCCC-3’ and Reverse, 5’-TGAGGGTGCAGGAAGTGAG-3’) for variant c.12572G > A (p.Arg4191His) were used to span the sequences. Sequencing was performed in both directions by BGI (Chongqing, China). The reference sequence of RYR1 used was GenBank NM_000540.3.
In silico analysis
In silico prediction tools (SIFT, PolyPhen-2, Mutation Taster, and LRT) were used to assess variant pathogenicity and potential functional effects. Online analysis tools (Human Splicing Finder, and NetGene2) were used to evaluated potential effects on splicing. NCBI Protein blast was used to analyze the conservation of amino acid residue alteration by comparing across different species.
Minigene splicing assay
We performed functional analysis of the c.2682G > C (p.Ile860_Pro894del) variant with an in vitro minigene assay. Genomic DNA of the parents was amplified using specific primers (RYR1_Ex20_F: 5′-CCCAAGCTTCCGCCATGGTGAATTCAAGT-3′, RYR1_Ex22_R: 5′-ATAGTTTAGCGGCCGCCTCTGGAAGGCTGTGGAAGT-3′) for exons 20, 21, 22 and intron sequences (3432 bp). The amplified minigene products were cloned into the pcDNA3.1(+) plasmid (Invitrogen) at the Hind III and Not I sites. After verifying by sequencing, we obtained a homozygous mutant-type hybrid minigene from the father (pcDNA3.1-MT) and a wild-type hybrid minigene from the mother (pcDNA3.1-WT). Then, hybrid minigenes were transfected into HEK293T cells transiently with the Lipofectamine 2000 reagent (Invitrogen). After 48 h, RNA was extracted from each plate of cells and RT‐PCR was performed with the same primers. RT‐PCR products were analyzed by electrophoresis analysis and Sanger sequencing.
G-banding analysis showed that the karyotype of the fetus II-5 and the couple were normal. No CNVs were identified by SNP-array in fetus II-5 and parents. Trio WES analysis found compound heterozygous variants in RYR1, a c.2682G > C (p.Ile860_Pro894del) variant that was inherited from the father and a c.12572G > A (p.Arg4191His) variant that was inherited from the mother. Subsequently, sanger sequencing revealed the same results (Fig. 2).
The c.2682G > C (p.Ile860_Pro894del) variant is novel (i.e., not reported in ExAC/gnomAD, ClinVar, 1000 Genomes, or HGMD databases and not published to date). This variant is located at the last nucleotide of exon 21 and is predicted to cause RYR1 alternative splicing by Mutation Taster, NetGene2 and HSF 3.0. To characterize the abnormal splicing, minigene Splicing Assay was performed.
Splicing products of c.2682G > C (p.Ile860_Pro894del) was analyzed by RT-PCR. Agarose gel electrophoresis showed that pcDNA3.1-WT plasmid produced a 361 bp band, while pcDNA3.1-MT generated a 256 bp band (Fig. 3A and Additional file 1). Sequencing analysis confirmed that the product of pcDNA3.1-WT was consistent with the reference sequence. In contrast, sequencing pcDNA3.1-MT indicated a complete skipping of exon 21 (Fig. 3A). That is, the c.2682 G > C (p.Ile860_Pro894del) variant caused a deletion of 105 bp in exon 21. According to the AMP/ACMG guidelines for the interpretation of sequence variants, the novel c.2682 G > C (p.Ile860_Pro894del) variant was assessed to be likely pathogenic.
The population frequency of this c.12572G > A (p.Arg4191His) variant is extremely low (MAF of 0.000008). The variant was previously detected (rs765730833), but had not been reported in individuals with RYR1-related disease. It was highly conserved and was predicted deleterious by SIFT, Polyphen‐2, and LRT. The c.12572G > A (p.Arg4191His) variant located in the C-terminal domain on exon 90, was highly conserved across species (Fig. 3B). According to the ACMG/AMP standards and guidelines, the variant is uncertain significance.
With advances in genetic analysis techniques, the spectrum of RYR1-RD has gradually expanded. RYR1 has been repeatedly reported to be associated with severe congenital onset in the neonatal period [9, 15,16,17,18,19,20]. A severe form of central core disease (CCD) had been reported in seven cases with fetal akinesia, multiple arthrogryposis and hydramnion during pregnancy. Ultimately, two fetuses died before birth and three infants died before their first birthday . In 2016, Lethal multiple pterygium syndrome (LMPS) was considered at the extreme end of the RYR1-RD spectrum. It usually presents as prenatal growth failure with pterygium present in multiple areas, severe arthrogryposis, akinesia and fetal hydrops . In this report, we describe two fetuses with fetal akinesia, polyhydramnios, multi-malformation and hydrops fetalis presented with some features of CCD and LMPS. Since there is no typical multiple pterygium phenotype and no fetal tissue remains for histopathological studies, we can’t identify the disease subtype. Based on fetuses’ phenotype and WES test results, we use RYR1-RD to describe the disease we found.
In these cases, there are both dominant inheritance and recessive inheritance. Dominant variants are found to concentrate at three “hot spots” : N-terminal (residues p.Met1-p.Arg614), central (p.Arg2163–p.Arg2458) and C-terminal (p.Arg4136-p.Pro4973). While recessive variants are distributed through the entire RYR1 gene. Dominant variants are clinically characterized by mild hypotonia and non-progressive weakness, while recessive variants may cause early-onset and severe myopathy. And, patients with at least one hypomorphic allele (nonsense, frameshift, indel, or splice) had increased disease severity [22, 23]. To date, the molecular disease patho-mechanisms of RYR1-RD caused by dominant versus recessive variants are not completely clear.
Here, we reported compound heterozygous variants (splicing and missense) in RYR1 in a fetus. It was consistent with recessive inheritance pattern in this family and leads to more severe phenotype. The variant (c.2682G > C) was the skipping of the entire exon and consequently the loss of the c.2578_2682 in the mRNA, resulting in an in-frame deletion (p.Ile860_Pro894del). The largest RyR1 domain was the myoplasmic domain, also known as the RyR1 foot region, which consisted of the first 3613 amino acid residues . The myoplasmic domain formed key inter-subunit interactions and housed binding sites for channel activity modulator proteins calmodulin, S100A1, FK506-binding protein (FKBP12), and dihydropyrine receptor [25,26,27]. Another variant c.2682G > A (p.Pro894 =) at the same site had been reported in an individual affected with multiminicore disease. Furthermore, previously missense substitutions within this deletion (c.2654G > T; p.Arg885Leu and c.2603G > A; p.Arg868His) had been reported in dominantly inherited central core myopathy . Thirty five amino acid deletions may result in an alteration of the protein’s secondary structure, which then may affect protein interactions. So, the novel splicing variant was classified as “likely pathogenic”.
The second missense variant (c.12572G > A) located in C-terminal ‘hot spots’ had been previously reported. The variant c.12572G > A resulted in an aromatic amino acid (His) was substituted for an aliphatic amino acid (Arg) at amino acid position 4191. The maternal p.Arg4191His variant affected a highly conserved amino acid located in the M2–M3 cytoplasmic loop of the channel release domain according to the RYR1 10-transmembrane domain-based model . In addition, the variation was located in a previously predicted EF-hand pair that constituted the conserved Ca2+-binding domain in RyRs . However, the p.Arg4191His had also been reported to be located in the thumb-and-forefinger domain of RyR1, which connects to the C-terminal extension of the S6 helix. The TaF domain is an ATP ligand binding site and close to putative Ca2+ binding site. It is part of the activation module for conformational changes upon activation and opening of RyR1 [24, 31]. Furthermore, the central domain contained amino acids 3668–4251 and constituted a docking station for a variety of allosteric regulatory factors including proteins and small molecules . Therefore, the p.Arg4191His directly or indirectly altered Ca2+ sensing or the interaction with RyR1 agonists or antagonists, thereby interfering with channel activity. According to the ACMG and AMP criteria, this variant was classified as “uncertain significance”.
Taken together, we describe a ‘severe neonatal form’ of RYR1-RD with the antenatal onset of symptoms, two fetuses in a non-consanguineous healthy Chinese family presented with an abnormality of the vertebral column, absence of stomach bubble, polyhydramnios, and multi-malformation during pregnancy. They were terminated at 28 and 25 weeks of gestational age, respectively. Trio-WES to identified missense variant and splicing variant in RYR1 and we confirmed the pathogenicity of the splicing variant using a simple in vitro reporter minigene assay. In our case, clinical diagnosis was difficult as we lacked fetal muscle biopsy specimens to confirm the disease and the parental biopsies were not available for morphological investigation. Further genetic tests could not be carried out on the other four fetuses in this family, as no samples were available. We can only infer the same molecular etiology by phenotypic recurrence. This is the limitation of this study. We advocate that future fetuses with similar clinical presentation are examined with muscle biopsy of affected muscle. Moreover, genome-wide sequencing, like WES, would be a helpful and effective testing in prenatal diagnosis, especially for those fetuses without definite clinical differential diagnosis. Although some RYR1-RD cases of the severe neonatal phenotype had been reported abroad, this is the first report of variants in the RYR1 gene involved in a severe form of RYR1-RD in China. We reveal new findings expanding the disease spectrum of RYR1 gene variants.
Availability of data and materials
The clinical case records of our patients are available from the corresponding author on reasonable request. The whole exome sequencing data generated during the current study have been deposited in the Genome Sequence Archive (GSA). The direct web links to this dataset is https://ngdc.cncb.ac.cn/gsa-human/s/1YQ2331a.
Santulli G, Lewis D, des Georges A, Marks AR, Frank J. Ryanodine receptor structure and function in health and disease. Sub-cell Biochem. 2018;87:329–52.
Phillips MS, Fujii J, Khanna VK, DeLeon S, Yokobata K, de Jong PJ, et al. The structural organization of the human skeletal muscle ryanodine receptor (RYR1) gene. Genomics. 1996;34(1):24–41.
Gillard EF, Otsu K, Fujii J, Khanna VK, de Leon S, Derdemezi J, et al. A substitution of cysteine for arginine 614 in the ryanodine receptor is potentially causative of human malignant hyperthermia. Genomics. 1991;11(3):751–5.
Lawal TA, Todd JJ, Meilleur KG. Ryanodine receptor 1-related myopathies: diagnostic and therapeutic approaches. Neurother J Am Soc Exp Neurother. 2018;15(4):885–99.
D’Arcy CE, Bjorksten A, Yiu EM, Bankier A, Gillies R, McLean CA, et al. King-denborough syndrome caused by a novel mutation in the ryanodine receptor gene. Neurology. 2008;71(10):776–7.
Garibaldi M, Rendu J, Brocard J, Lacene E, Fauré J, Brochier G, et al. “Dusty core disease” (DuCD): expanding morphological spectrum of RYR1 recessive myopathies. Acta Neuropathol Commun. 2019;7(1):3.
Böhm J, Leshinsky-Silver E, Vassilopoulos S, Le Gras S, Lerman-Sagie T, Ginzberg M, et al. Samaritan myopathy, an ultimately benign congenital myopathy, is caused by a RYR1 mutation. Acta Neuropathol. 2012;124(4):575–81.
Løseth S, Voermans NC, Torbergsen T, Lillis S, Jonsrud C, Lindal S, et al. A novel late-onset axial myopathy associated with mutations in the skeletal muscle ryanodine receptor (RYR1) gene. J Neurol. 2013;260(6):1504–10.
McKie AB, Alsaedi A, Vogt J, Stuurman KE, Weiss MM, Shakeel H, et al. Germline mutations in RYR1 are associated with foetal akinesia deformation sequence/lethal multiple pterygium syndrome. Acta Neuropathol Commun. 2014;2:148.
Lawal TA, Todd JJ, Witherspoon JW, Bönnemann CG, Dowling JJ, Hamilton SL, et al. Ryanodine receptor 1-related disorders: an historical perspective and proposal for a unified nomenclature. Skeletal muscle. 2020;10(1):32.
Li H, Durbin R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinform (Oxford, England). 2010;26(5):589–95.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.
Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16): e164.
Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American college of medical genetics and genomics and the association for molecular pathology. Genet Med Off J Am Coll Med Genet. 2015;17(5):405–24.
Romero NB, Monnier N, Viollet L, Cortey A, Chevallay M, Leroy JP, et al. Dominant and recessive central core disease associated with RYR1 mutations and fetal akinesia. Brain J Neurol. 2003;126(Pt 11):2341–9.
Monnier N, Laquerrière A, Marret S, Goldenberg A, Marty I, Nivoche Y, et al. First genomic rearrangement of the RYR1 gene associated with an atypical presentation of lethal neonatal hypotonia. Neuromuscul Disord NMD. 2009;19(10):680–4.
Bevilacqua JA, Monnier N, Bitoun M, Eymard B, Ferreiro A, Monges S, et al. Recessive RYR1 mutations cause unusual congenital myopathy with prominent nuclear internalization and large areas of myofibrillar disorganization. Neuropathol Appl Neurobiol. 2011;37(3):271–84.
Bharucha-Goebel DX, Santi M, Medne L, Zukosky K, Dastgir J, Shieh PB, et al. Severe congenital RYR1-associated myopathy: the expanding clinicopathologic and genetic spectrum. Neurology. 2013;80(17):1584–9.
Kariminejad A, Ghaderi-Sohi S, Hossein-Nejad Nedai H, Varasteh V, Moslemi AR, Tajsharghi H. Lethal multiple pterygium syndrome, the extreme end of the RYR1 spectrum. BMC Musculoskelet Disord. 2016;17:109.
Mauri E, Piga D, Govoni A, Brusa R, Pagliarani S, Ripolone M, et al. Early findings in neonatal cases of RYR1-related congenital myopathies. Front Neurol. 2021;12: 664618.
Treves S, Anderson AA, Ducreux S, Divet A, Bleunven C, Grasso C, et al. Ryanodine receptor 1 mutations, dysregulation of calcium homeostasis and neuromuscular disorders. Neuromuscul Disord NMD. 2005;15(9–10):577–87.
Chang X, Wei R, Wei C, Liu J, Qin L, Yan H, et al. Correlation of phenotype-genotype and protein structure in RYR1-related myopathy. Front Neurol. 2022;13: 870285.
Todd JJ, Sagar V, Lawal TA, Allen C, Razaqyar MS, Shelton MS, et al. Correlation of phenotype with genotype and protein structure in RYR1-related disorders. J Neurol. 2018;265(11):2506–24.
des Georges A, Clarke OB, Zalk R, Yuan Q, Condon KJ, Grassucci RA, et al. Structural basis for gating and activation of RyR1. Cell. 2016;167(1):145–57.
Hernández-Ochoa EO, Pratt SJP, Lovering RM, Schneider MF. Critical role of intracellular RyR1 calcium release channels in skeletal muscle function and disease. Front Physiol. 2015;6:420.
Wright NT, Prosser BL, Varney KM, Zimmer DB, Schneider MF, Weber DJ. S100A1 and calmodulin compete for the same binding site on ryanodine receptor. J Biol Chem. 2008;283(39):26676–83.
Zalk R, Lehnart SE, Marks AR. Modulation of the ryanodine receptor and intracellular calcium. Annu Rev Biochem. 2007;76:367–85.
Snoeck M, van Engelen BG, Küsters B, Lammens M, Meijer R, Molenaar JP, et al. RYR1-related myopathies: a wide spectrum of phenotypes throughout life. Eur J Neurol. 2015;22(7):1094–112.
Du GG, Sandhu B, Khanna VK, Guo XH, MacLennan DH. Topology of the Ca2+ release channel of skeletal muscle sarcoplasmic reticulum (RyR1). Proc Natl Acad Sci USA. 2002;99(26):16725–30.
Xiong L, Zhang JZ, He R, Hamilton SL. A Ca2+-binding domain in RyR1 that interacts with the calmodulin binding site and modulates channel activity. Biophys J. 2006;90(1):173–82.
Woll KA, Van Petegem F. Calcium-release channels: structure and function of IP(3) receptors and ryanodine receptors. Physiol Rev. 2022;102(1):209–68.
Bai XC, Yan Z, Wu J, Li Z, Yan N. The central domain of RyR1 is the transducer for long-range allosteric gating of channel opening. Cell Res. 2016;26(9):995–1006.
We would like to thank the family members who provided samples and clinical information for this study.
The study was supported by grants from Chongqing Natural Science Foundation (cstc2019jcyj-msxmx0318); “Kuanren talents” project of the Second Affiliated Hospital of Chongqing Medical University (13-003-003); Scientific and Technological Research Program of Chongqing Municipal Education Commission (Grant No. KJQN202200420).
Ethics approval and consent to participate
The study was conducted in accordance with the principles of the Declaration of Helsinki. All participants provided written informed consent. The study protocols were approved by the ethics committee from the Second Affiliated Hospital of Chongqing Medical University.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhao, Q., Li, X., Liu, L. et al. Prenatal diagnosis identifies compound heterozygous variants in RYR1 that causes ultrasound abnormalities in a fetus. BMC Med Genomics 15, 202 (2022). https://doi.org/10.1186/s12920-022-01358-x
- Ryanodine receptor type 1 gene
- RYR1-related disorders
- Ultrasound abnormalities
- Splice-site variant
- Whole-exome sequencing