- Open Access
Next generation sequencing in a large cohort of patients presenting with neuromuscular disease before or at birth
Orphanet Journal of Rare Diseasesvolume 10, Article number: 148 (2015)
Fetal akinesia/hypokinesia, arthrogryposis and severe congenital myopathies are heterogeneous conditions usually presenting before or at birth. Although numerous causative genes have been identified for each of these disease groups, in many cases a specific genetic diagnosis remains elusive. Due to the emergence of next generation sequencing, virtually the entire coding region of an individual’s DNA can now be analysed through “whole” exome sequencing, enabling almost all known and novel disease genes to be investigated for disorders such as these.
Genomic DNA samples from 45 patients with fetal akinesia/hypokinesia, arthrogryposis or severe congenital myopathies from 38 unrelated families were subjected to next generation sequencing. Clinical features and diagnoses for each patient were supplied by referring clinicians. Genomic DNA was used for either whole exome sequencing or a custom-designed neuromuscular sub-exomic supercapture array containing 277 genes responsible for various neuromuscular diseases. Candidate disease-causing variants were investigated and confirmed using Sanger sequencing. Some of the cases within this cohort study have been published previously as separate studies.
A conclusive genetic diagnosis was achieved for 18 of the 38 families. Within this cohort, mutations were found in eight previously known neuromuscular disease genes (CHRND, CHNRG, ECEL1, GBE1, MTM1, MYH3, NEB and RYR1) and four novel neuromuscular disease genes were identified and have been published as separate reports (GPR126, KLHL40, KLHL41 and SPEG). In addition, novel mutations were identified in CHRND, KLHL40, NEB and RYR1. Autosomal dominant, autosomal recessive, X-linked, and de novo modes of inheritance were observed.
By using next generation sequencing on a cohort of 38 unrelated families with fetal akinesia/hypokinesia, arthrogryposis, or severe congenital myopathy we therefore obtained a genetic diagnosis for 47 % of families. This study highlights the power and capacity of next generation sequencing (i) to determine the aetiology of genetically heterogeneous neuromuscular diseases, (ii) to identify novel disease genes in small pedigrees or isolated cases and (iii) to refine the interplay between genetic diagnosis and clinical evaluation and management.
Fetal akinesia deformation sequence (FADS) or Pena Shokeir syndrome, characterized by intrauterine growth retardation, contractures, craniofacial anomalies, limb anomalies, pulmonary hypoplasia and polyhydramnios, results from reduced movement in utero [1, 2]. A number of other fetal akinesia syndromes overlap phenotypically with FADS. These include the lethal congenital contracture syndromes, multiple pterygium syndromes, and arthrogryposis multiplex congenita , in which the clinical findings are dependent upon the time of onset of the dyskinesia, earlier onset being associated with a more severe phenotype . It’s thought that more than 50 % of all causes of fetal akinesia are of neuromuscular origin ; at least 30 causative genes have been identified, involving all points along the neuromuscular axis (motor neurons, peripheral nerves, neuromuscular junction and the skeletal muscle regulatory and contractile apparatus) [5–7].
Arthrogryposis refers to non-progressive congenital joint contractures in >1 area of the body, and has been described in more than 300 specific disorders [6, 8]. Arthrogryposis is thought to result from reduced fetal movement, and affects approximately 1 in 3,000 live births [8, 9]. There is a range of disease severity: severe cases present with arthrogryposis multiplex congenita, which is lethal prior to or at birth, while milder cases with a longer life expectancy may have predominantly distal involvement [8, 9]. The distal arthrogryposes are a group of disorders with contractures primarily involving the extremities of the body, often associated with camptodactyly, hypoplastic or absent flexion creases, and talipes equinovarus [10, 11]. There are ten distinct subtypes of distal arthrogryposis, for which seven causative genes have been identified: ECEL1 (OMIM 605896), MYH3 (OMIM 160720), MYH8 (OMIM 160741), PIEZO2 (OMIM 613692), TNNI2 (OMIM 191043), TNNT3 (OMIM 600692) and TPM2 (OMIM 190990) [11–13].
The congenital myopathies are a diverse group of disorders, characterised by skeletal muscle dysfunction (most often weakness and hypotonia), with specific morphological features on skeletal muscle biopsies [14, 15]. Three distinct major groups are recognized based upon the presence of one or more major histopathological features: centronuclear myopathy, core myopathy and nemaline myopathy (NEM), although there is extensive overlap in both genotype and phenotype within and between these groups [16, 17]. While muscle biopsy remains critical for diagnosis, there can be overlap in the morphological abnormalities seen in these conditions, and marked variability in their clinical progression and severity [14, 15]. The clinical spectrum of the congenital myopathies ranges from severe fetal akinesia to adult-onset progressive weakness. Typical features of these conditions include proximal weakness, respiratory insufficiency, facial weakness, skeletal deformities such as hip dislocation and deformities of the feet, feeding difficulties, hypotonia and delayed motor milestones , however hypertonic cases are also encountered [15, 18].
More than 15 disease genes are known to cause congenital myopathies. However, many cases remain genetically unresolved, suggesting further heterogeneity [5, 7, 12, 19–21]. This study aimed to assess the potential of next generation sequencing technologies to identify causative genes in small families or isolated probands presenting with fetal hypokinesia, arthrogryposis or a severe congenital myopathy.
Subject information and study ethics approval
Informed consent was given for participation in this study, which was approved by the Human Research Ethics Committee of the University of Western Australia, Perth, Western Australia, Australia.
Exome sequencing for this study was performed at the Lotterywest State Biomedical Facility Genomics Node (LSBFG) in Perth, Australia. Exome sequencing was performed on the 5500XL SOLiD™ system (Applied Biosystems), as described elsewhere [20, 22–24], and the Ion Proton™ (Ampliseq chemistry, Life Technologies) (Family 16 and 38). For AmpliSeq exome sequencing, 100 ng of DNA from the probands was amplified in 12 PCR pools and sequencing adaptors ligated. The library was then purified using AMPure beads (Beckman Coulter), and amplified using Platinum® High-Fidelity Taq Polymerase. The amplified library was again purified with AMPure beads and analysed on a 2100 Bioanalyser (Agilent Technologies Genomics). Libraries were diluted to 18-26pM and attached to Ion Sphere™ Particles using an Ion Proton™ Template 200 v3 kit and sequenced on a P1 sequencing chip on an Ion Proton sequencer™ (Ion Sequencing 200 kit v3) in pools of two.
Targeted capture and sequencing of neuromuscular disease genes by next generation sequencing
Neuromuscular sub-exomic sequencing (NSES) was also performed at the LSBFG. The NSES panel comprised those genes listed within the December 2012 freeze of Neuromuscular Disorders gene table  in which the disease-causing mutations could be identified by next generation sequencing, some unpublished candidate disease genes identified by our group and others and 59 cardiomyopathy genes. NSES analysis was performed on DNA from the probands using the Ion Proton™ sequencer (Life Technologies), as previously described . For NSES, 2 μg of DNA was captured in pools of 16 DNA samples using a custom TargetSeq™ (Life Technologies) capture system, enriching for the 336 known and candidate neuromuscular and cardiomyopathy disease genes. These captured pools were then sequenced in batches of 16 using an Ion P1 200 V2 sequencing kit (Life Technologies) for 520 flows.
Variant calling was performed against the GRCh37 human reference genome, using LifeScope™ 2.5 (exome sequencing) and Torrent Suite V 3.6.2 (NSES) (Life Technologies). Data was filtered using an ANNOVAR annotation software suite. Variants were annotated using the EncodeGencode gene annotation set. Variants were filtered against the 1000 Genomes database (2012 release, ) and the dbSNP137 common database, and variants with a frequency of >0.5 % were excluded. Variants were then filtered against an in-house common variant list and were checked against the HGMD professional database to identify any known disease-causing mutations. The frequencies of candidate disease variants in the 1000 Genomes Project, Exome Variant Server (http://evs.gs.washington.edu/EVS/) and ExAC Browser (http://exac.broadinstitute.org) were also determined. Pathogenicity predictions were made using online prediction software programs: SIFT, PolyPhen , and MutationTaster .
The LSBFG has a cut-off of 90 % of on-target regions covered to 20-fold or greater for the neuromuscular panel (NSES) and 80 % covered to 20-fold or greater for exome sequencing, however some samples, especially early samples, did not achieve these cut-offs (Additional file 1: Table S1). There was no significant difference in the average coverage (mean ± SEM) of exome sequencing data for genetically resolved (80 ± 14-fold; n = 15) versus unresolved cases (70 ± 8-fold; n = 23). For the NSES panel, average coverages were 220 ± 23-fold (n = 6) for resolved cases versus unresolved cases (195 ± 13-fold, n = 9). Thus coverage is unlikely to contribute to the lack of a genetic diagnosis in most cases.
Sanger confirmation and co-segregation studies
PCR amplification and Sanger sequencing was performed to verify potential mutations identified by next generation sequencing. Co-segregation was also verified for all existing family members where available. Primers were based on genomic and cDNA sequences obtained from the UCSC Human Genome Browser (http://genome.ucsc.edu/) and Ensembl (http://www.ensembl.org/). Primer sequences and conditions are available upon request. Sanger sequencing data was processed by LSBFG and results viewed using CodonCode Aligner software.
Functional studies of the CHRND missense substitution
The mutation CHRND p.Cys257Arg was directly introduced into the wild-type human delta subunit cDNA in the vector pcDNA3.1/hygro (−) by site-directed mutagenesis (QuikChange® Site-Directed Mutagenesis Kit, Stratagene, Amsterdam, The Netherlands). Primer sequences can be obtained on request. To confirm the presence of the introduced mutation, and to rule out any errors, the construct was subjected to Sanger sequencing.
Wild-type and mutant human AChR δ-subunits cDNAs in the vector pcDNA3.1/hygro (−) (Life Technologies, V875–20) were used for transfection studies.
Wild-type and mutant AChR δ-subunit cDNAs, in combination with wild-type α-, β-and ɛ-subunit cDNAs, were transfected into HEK 293 cells grown on six-well tissue culture plates using polyethyleneimine. Surface AChR expression was determined 2 days post-transfection by incubating cells in 10 nM 125I-α-bungarotoxin (125I-α-BuTx) with 1 mg/ml BSA for 30 minutes. Cells were washed three times with PBS and extracted in 1 % Triton X-100, in 10 mM Tris–HCl (pH 7.4), 100 mM NaCl, 1 mM EDTA and 125I-α-BuTx binding determined using a gamma counter.
Results and discussion
A total of 45 subjects from 38 families (including ten consanguineous pedigrees) diagnosed with FADS, arthrogryposis, or a severe congenital myopathy were included in this study. Of these seven probands were submitted for NSES only, eight families had probands sequenced using both NSES and exome sequencing, and 23 families underwent only exome sequencing (Additional file 1: Table S1). Families were grouped into three disease entities: FADS (n = 9), arthrogryposis (n = 13), and severe congenital myopathies (n = 16). Clinical details for the genetically resolved families are summarized in Table 1.
A conclusive genetic diagnosis was achieved for 18/38 families (47 %, Table 2). This included two kindreds with FADS, six with arthrogryposis and 10 presenting with a congenital myopathy. From these results, autosomal dominant (n = 1), autosomal recessive (n = 15), de novo (n = 1) and X-linked (n = 1) modes of inheritance were identified. Mutations were identified in eight previously known neuromuscular disease genes. As part of this cohort study, four then novel disease genes were initially identified from five families (Families 3, 4, 5, 12 and 38) in the cohort and these families have been previously published: GPR126 (Family 3) , KLHL40 (Family 10 and 17; OMIM 615340) , KLHL41 (Patient ID: D12-203; OMIM 607701)  and SPEG (Patient ID: P3; OMIM 615950) .
Mutations in fetal hypokinesia and congenital myopathy genes
Since our initial publication of KLHL40 as a novel NEM gene, two further families within our cohort were shown to have mutations in KLHL40 (Families 16 and 20). A previously-unpublished homozygous nonsense mutation in KLHL40 (exon 1, c.46C>T, p.Gln16*) was identified in a proband from consanguineous parents (Family 16, Fig. 1a). This proband was born by emergency Caesarean section at 35/40 weeks gestation and presented with severe arthrogryposis, congenital fractures, respiratory insufficiency and complete akinesia. An initial clinical diagnosis of spinal muscular atrophy type 0 was made, but both light and electron microscopy of the child’s muscle biopsy demonstrated miliary nemaline bodies (Fig. 2), adding to the body of evidence suggesting that miliary nemaline bodies are a good indicator suggesting KLHL40 as the causative gene.
The proband in Family 20 was born to consanguineous parents (Fig. 1b) by Caesarean section at 37 weeks gestation. He had profound hypotonia, an absent gag reflex, myopathic facies, and was ventilated from birth, but survived only a few days. His muscle biopsy showed numerous nemaline bodies. The family history included two previous miscarriages, two neonatal deaths and a sibling who died at seven months of age with suspected NEM (light microscopy indicated rods, but electron microscopy was not performed). No mutations were found on Sanger sequencing of ACTA1, but NSES showed a novel homozygous missense mutation in KLHL40 (exon 1, c.931C>A, p.Arg311Ser) affecting the same highly conserved amino acid residue as that in Family 3 (Fig. 3a).
The second male child of a non-consanguineous family (Family 14, Fig. 1c) was born after an uncomplicated pregnancy, by emergency Caesarean section for failure to progress. The baby was weak and hypotonic at birth, was very long (reported >90th percentile for length with weight 10-25th percentile), had advanced bone age, and initially required intubation. By age 5 months the infant’s strength and spontaneous movement improved markedly, but he had significant residual weakness and bulbar dysfunction. NSES identified a known missense mutation (exon 8, c.595C>T, p.Pro199Ser) in the myotubularin gene (MTM1; OMIM 300415)  associated with myotubular myopathy. A muscle biopsy taken at 10 weeks of age revealed hypoplastic myofibres, some with internal nuclei, typical features of myotubular myopathy (OMIM 310400) . However, enzyme staining showed reduced central staining in some myofibres, while electron microscopy showed foci of sarcomeric dissolution, suggestive of cores. IHC for myosin confirmed the preservation of type II/fast myofibres and numerous small type I myofibres. Thus a diagnosis of congenital myopathy with fibre-type disproportion and occasional minicores had been suggested (Fig. 2d-e). This highlights that MTM1 cases can present with congenital weakness and muscle biopsies displaying features of fibre type disproportion and minicores.
The proband and affected sibling of Family 6 (Fig. 1d), were born to consanguineous parents. The proband was born at 30 weeks gestation with profound hypotonia, facial weakness, dysmorphic features and ambiguous genitalia, after a pregnancy complicated by fetal hypokinesia. He died at 3 weeks of age. A subsequent pregnancy with a female sibling was complicated by polyhydramnios. At birth there was minimal limb movement, respiratory distress necessitating mechanical ventilation, subcutaneous oedema, contractures of the hips and knees and camptodactyly of the fingers. She died at 4 weeks of age. Maternal testing for myotonic dystrophy (DM1) was negative. Vastus lateralis biopsies from both babies showed non-specific abnormalities of myofibre typing, with type II myofibre predominance and numerous small myofibres. Occasional minicores and cores were seen in the proband but not his sibling. Neither had nemaline bodies or histologic features of myotubular myopathy. Exome sequencing performed on the proband revealed two mutations in the ryanodine receptor gene (RYR1, OMIM 180901): a novel heterozygous missense mutation affecting a highly conserved amino acid (Fig. 3b) (exon 36, c.5989G>A, p.Glu1997Lys) and a heterozygous previously-reported nonsense mutation (exon 20, c.2455C>T, p.Arg819* ). Sanger sequencing confirmed these mutations and showed co-segregation with disease. Thus in this instance, the consanguinity does not appear to be a contributing factor in the siblings’ disease. The nonsense mutation was previously identified in a 49-year old ambulant patient with a moderate form of slowly-progressive myopathy with cores . That patient also harboured a previously identified heterozygous missense mutation (p.Arg4558Gln) [33, 34]. Thus the same nonsense mutation, in combination with different missense mutations, can result in variable phenotypes, from fetal hypokinesia and death in the perinatal period, to a mild delay in motor milestones and normal life expectancy.
The proband in non-consanguineous Family 8 (Fig. 1e) presented with non-immune hydrops fetalis and arthrogryposis, and was stillborn at 26 weeks gestation. Autopsy showed multiple contractures and reduced muscle bulk. Microscopically, there were marked dystrophic changes in all muscles examined (Fig. 2g). The contactin-1 gene (CNTN1, OMIM 600016) was Sanger sequenced but no mutations were identified. A subsequent pregnancy with a male fetus was complicated by polyhyhydramnios, contractures, and in utero fetal demise at 29 weeks gestation. Both affected individuals were diagnosed with FADS and congenital muscular dystrophy. Exome sequencing of the proband identified two previously reported heterozygous mutations in the RYR1 gene; a nonsense mutation (exon 41, c.6721C>T, p.Arg2241* ) and a missense mutation (exon 101, c.14645C>T, p.Thr4882Met ) associated with multiminicore disease and core rod disease, respectively. Sanger sequencing confirmed compound heterozygosity in both affected individuals, and showed that both parents were carriers and that the unaffected sibling did not harbour either mutation.
The affected individual in Family 13 was born to non-consanguineous Turkish parents (Fig. 1f) after reports of reduced intrauterine movement. He was delivered at term by Caesarean section, due to poor positioning. At birth bilateral humeral fractures were noted. He had a diagnosis of osteogenesis imperfecta, and followed up with alendronate treatment. He was referred to pediatric neurology outpatient clinic at the age of 19 months when the parents had concerns in terms of hypotonia and delay in motor developmental milestones. At the time, he remained hypotonic with a myopathic face and high-arched palate. He had axial and vertical hypotonia, head lag, facial weakness, and absence of deep tendon reflexes. He could sit but not stand. The muscle biopsy showed muscle tissue embedded in fibro-adipose tissue with severe non-specific myopathic changes (Fig. 2f). There were hypertrophic and atrophic myofibres, central nuclei, type II myofibre predominance and some core-like regions on oxidative enzyme stains. Exome sequencing of the proband revealed two pathogenic missense mutations in the RYR1 gene, (exon 15, c.1589G>A, p.Arg530His , exon 101, c.14645C>T, p.Thr4882Met ) which had previously been associated with central core disease/malignant hyperthermia (MH) and core rod myopathy, respectively. The p.Arg530His substitution was inherited paternally, thus the presence of this MH (OMIM 145600) susceptibility mutation in both the proband and asymptomatic father changes their clinical management.
Thus, affected individuals in three families (Family 6, 8, and 13), harbored compound heterozygous mutations in RYR1. Disease severity was much greater in the two families possessing a nonsense (null) mutation as well as a missense mutation (Family 6 and 8), resulting in death at or soon after birth. The affected individual in the third RYR1 family, (Family 13), possessed two missense mutations, and survived infancy, albeit with severe muscle weakness and motor delay. He had a rather static improving course with physiotherapy. These findings mirror those of recent publications expanding the phenotypes associated with recessive RYR1 disease to include arthrogryposis multiplex congenita and fetal akinesia [9, 38, 39]. Despite RYR1 originally being described as a disease gene for central core disease and minicore disease, cores are seen in only a minority of recessive RYR1 cases, and are less likely to be seen in cases with hypomorphic (null) muations [38, 39]. In this study, cores were not a prominent feature in two of the recessive RYR1 families, both of which harboured a hypomorphic mutation.
A consanguineous family (Family 9) presented early in pregnancy with monoamniotic male twins (Fig. 1g) and a history of a previous fetus therapeutically aborted due to hydrops fetalis at 19 weeks gestation. Ultrasound scanning revealed severe hydrops in both fetuses, and the pregnancy was terminated at 16 weeks gestation. Post-mortem analysis of both twins showed bilateral joint contractures, bilateral talipes, multiple pterygia, hypertelorism and cystic hygromas. Muscle biopsies were not taken. A diagnosis of fetal akinesia with lethal multiple pterygia syndrome was made. Karyotyping showed a normal 46XY karyotype, with no apparent genomic imbalance. Exome sequencing was performed on one twin, and a novel homozygous nonsense mutation (exon 29, c.2920C>T, p.Arg974*) in the nebulin gene (NEB; OMIM 161650) was identified. Sanger sequencing confirmed that both twins were homozygous for this mutation and that each parent was a carrier (Fig. 1). This mutation was included in the recent NEB mutation update . Although this case was diagnosed as FADS/lethal multiple pterygia syndrome, recessive mutations in the NEB gene are a known cause of NEM, which in severe cases can have a FADS phenotype . Without a muscle biopsy however, it cannot be determined whether these cases had nemaline myopathy.
In three additional families, diagnosed with NEM presenting with fetal akinesia, single heterozygous pathogenic mutations were identified in NEB by either exome sequencing or NSES (Table 3). In Family 17 a known splice-site mutation (intron 5, c.78+1G > A, ) was identified, and in Family 19 a known frameshift mutation (exon 55, c.7523_7526del4, p.Ile2508Thrfs*14, ), was identified, both of which are associated with NEM. In Family 18, a previously unpublished nonsense mutation (exon 29, c.2864G>A, p.Trp955*) was identified. A common deletion of exon 55 of NEB, originating within the Ashkenazi Jewish population, is known to cause a severe NEM phenotype . A heterozygous deletion of this exon would not be identifiable through next generation sequencing techniques. Deletion analysis was performed on the affected individuals of Family 17 and Family 18, which confirmed they did not have a deletion of this exon. The proband in Family 19 could not have harboured a deletion of exon 55, since the exon 55 variant identified in this proband was heterozygous. Although only single heterozygous mutations were identified in these three severe NEM cases, given their severity and the absence of likely pathogenic variants in the other known NEM genes, it is likely that they are harbouring a second pathogenic NEB variant that was not identified by next generation sequencing. In support of AR NEM, Family 18 and 19 both had a previously affected fetus. In further support that these cases (three of nine NEM families, 33 %) are harbouring an additional pathogenic NEB variant, only one truncating NEB variant was identified by next generation sequencing in non-NEM cases, of which we have sequenced and analysed in excess of >500 probands (~0.2 %). Due to the highly repetitive nature of exons 83–105 of NEB, next generation sequencing is unable to accurately sequence and map this region; in addition, next generation sequencing data is not reliable for the detection of small CNVs. However a targeted NEB array CGH has been developed as an adjunct to overcome these limitations  and has recently identified a recurrent CNV within this triplicated repeat .
Therefore, of the nine NEM cases in our cohort, five cases had mutations in the newly described genes KLHL40 and KLHL41, and an additional three cases are thought likely to harbour a second pathogenic mutation in NEB. It is likely that many undiagnosed NEM cases are due to mutations in NEB, however due to its size it has not been routinely screened. With the introduction of next generation sequencing techniques, more NEB-related NEM cases are beginning to be identified. This may mean that there are not as many new NEM genes to find as might have been thought.
A non-consanguineous family (Family 2) presented with recurrent fetal akinesia and multiple pterygium syndrome . We identified compound heterozygous mutations in the gene GBE1, a known splice site mutation (intron 5, c.691+2T>C) associated with a non-lethal neonatal glycogenosis type IV, and a missense mutation (exon 7, c.956A>G, p.His319Arg). This report extended the phenotypic spectrum of GBE1 disease to include lethal multiple pterygium syndrome .
Mutations in known disease genes for arthrogryposes
The proband in Family 10 was the first child to non-consanguineous parents, born following an IVF pregnancy, (Fig. 1h). A routine 20-week ultrasound identified bilateral fetal talipes. Chromosome microarray was normal. The fetal phenotype evolved with polyhydramnios, fetal micrognathia and an absence of hand movements noted at 32 weeks. The polyhydramnios required three amnioreduction procedures. The male infant was delivered by elective Caesarian section for placenta praevia at 37 weeks gestation and weighed 2.82 kg. He was intubated and ventilated at 10 minutes for apnoea and poor respiratory effort after APGARS of 51, 65 and 710. He had micrognathia, cryptorchidism, a left single palmar crease, bilateral talipes, moderate large joint contractions, hypotonia, an absent gag/suck, and paucity of movement. He developed a weak suck and infrequent antigravity movement of the fingers after a week. Prader-Willi syndrome, SMA and myotonic dystrophy were excluded. Endocrine and metabolic investigations were normal as was the ophthalmologic examination. Brain MRI showed a right MCA infarct in the context of positive maternal serology for SLE. Multiple attempts to extubate the patient to CPAP failed. Ongoing ventilatory support was considered futile and was withdrawn at 4 weeks of age. Exome sequencing was performed and two novel heterozygous mutations were identified in CHRND (OMIM 100720) that encodes the delta-subunit of the acetylcholine receptor (AChR) . A missense mutation (c.769T>C) in exon 7 that resulted in substitution of a highly conserved amino acid (p.Cys257Arg, Fig. 3c) and a frameshift mutation in exon 11 (c.1315delG, p.Val439Trpfs*11). To our knowledge neither of these mutations have been previously reported and are not listed in the CHRND locus-specific database (http://www.dmd.nl/nmdb/home.php?select_db=CHRND). Sanger sequencing confirmed the presence of the mutations in the affected individual and showed each parent was a carrier of one of the variants. The mother conceived a second time, naturally. Prenatal diagnosis was performed and the fetus had both variants. The pregnancy was terminated.
Studies in HEK cells, found that cell surface expression levels of AChRs harbouring the δC257R subunit to be approximately 20 % of wild-type (Fig. 4). This result is consistent with the c.769T>C mutation (in combination with c.1315delG, p.Val439Trpfs*11 on the second allele) underlying a congenital myasthenic syndrome due to AChR deficiency. The mother is currently pregnant and is approaching term with a healthy fetus following PGD. Mutations of CHNRD typically result in congenital myasthenic syndromes (OMIM 608930 (fast-channel) and 601462 (slow-channel) [48, 49]). but have also more recently been associated with lethal multiple pterygium syndrome . In two families presenting with recurrent lethal multiple pterygium syndrome, resulting in terminations during the second trimester of pregnancy, null mutations of CHRND were identified (one consanguineous family with a homozygous p.Trp57* mutation and one with compound heterozygous p. Phe74Leu and p.Arg464* mutations). Substitutions of amino acids in close proximity to Cys257 have been shown to cause congenital myasthenia and impaired channel function (p.Pro250Gln  and p.Ser268Phe ).
The affected female individual in Family 15 was born to unrelated parents (Fig. 1i). At birth there was arthrogryposis with distinctive shin dimples. The clinical picture of this patient is presented in Hall et al., (Patient 10) . NSES was performed on the proband and revealed a known frequent homozygous frameshift mutation (exon 5, c.459dupA, p. Val154Serfs*24) in the gene encoding the gamma-subunit of the AChR (CHRNG; OMIM 100730) [53, 54]. Sanger sequencing confirmed the presence of the mutation in the affected individual, as well as showing each parent had the mutation in the heterozygous state. Given the unique presentation of arthrogryposis with shin dimples in this case and others harbouring CHRNG mutations , CHRNG should be considered in individuals presenting with this particular phenotype.
The proband in Family 1 was born from unaffected parents and has an unaffected sibling (Fig. 1j). He presented with Freeman-Sheldon Syndrome (DA2A, OMIM 193700)  and on examination at 2 years of age he showed some facial features and proximal weakness. Exome sequencing of the proband demonstrated heterozygosity for a mutation in MYH3 (exon 18, c.2015G>A, p.Arg672His; OMIM 160720 ;) previously associated with Freeman-Sheldon syndrome. Sanger sequencing confirmed the presence of this mutation in the proband as well as its absence in the unaffected sibling and both unaffected parents, confirming the mutation was de novo.
The male proband from Family 11, was born from an unaffected mother, but affected father (Fig. 1k). Both proband and the father were born with a very typical Sheldon-Hall distal arthrogryposis (DA2B) phenotype. Sheldon-Hall syndrome can be caused by autosomal dominant or de novo mutations in a number of genes. In this kindred, screening of TPM2 and TNNI2 identified no mutations. On NSES, however, a known heterozygous MYH3 mutation (exon 12, c.1123G>A, p.Glu375Lys ) previously associated with Freeman-Sheldon syndrome, was identified. Sanger sequencing confirmed the mutation in both the proband and his affected father, confirming autosomal dominant inheritance.
Non-consanguineous Family 7, previously described in , was also part of this cohort study. The proband was born from an uncomplicated pregnancy with extended hips, multiple arthrogrypotic features, multiple pterygium, adducted thumbs and bilateral ptosis. The couple presented when pregnant again, and on ultrasound at 20 weeks the fetus appeared to have similar features to those of the proband. The pterygia and ptosis led to consideration of multiple pterygium syndrome (Table 1). Exome sequencing revealed compound heterozygous mutations in ECEL1 (OMIM 605896), a missense substitution (c.1531G>A, p.Gly511Ser) and a essential splice-site mutation (c.1797-1G>A). Mutations in ECEL1 are associated with distal arthrogryposis type 5D (OMIM 615065), and the clinical presentation was in keeping with those recently described for DA5D [13, 57], although pterygia was a more prominent feature in this family.
In another cohort study, Laquerriere et al. identified two novel genes (CNTNAP1 and ADCY6) for severe arthrogryposis multiplex congenita (AMC) by exome sequencing, and achieved a genetic diagnosis for 24 of 31 multiplex and/or consanguineous AMC families studied (>75 %). This highlights the importance of working with well-phenotyped cohorts . Mutations in CNTNAP1 were identified in four of their 31 families, suggesting that mutations in this gene underlie a significant proportion of recessive AMC cases.
Results from our study, and that of Laquerriere et al., suggest that there are further arthrogryposis disease genes to be identified . ADCY6 and CNTNAP1 are both involved in axonal function , as is ECEL1 . GPR126 is critical for myelination of peripheral nerves  and we identified AMC patients with loss-of-function mutations in GPR126 . Genes involved in axonal function should therefore be considered as candidates for arthrogryposis, in addition to skeletal muscle contractile proteins.
Our study highlights the widening spectrum of phenotypes associated with mutations in known fetal akinesia, arthrogryposis and myopathy genes, as is increasingly demonstrated for other neuromuscular disorders [26, 59, 60]. As sequencing of targeted gene panels or exome sequencing becomes the mainstay of genetic diagnostics [61, 62], it is likely that there will be greater broadening of genotype-phenotype correlations for neuromuscular diseases. With the overwhelming amount of genetic information obtained via next generation sequencing, the reliability of meticulously curated locus-specific databases, the availability of large exome datasets from ethnically matched reference populations and appropriate functional and/or protein studies will be critical to obtaining accurate genetic diagnosis. Given that numerous novel disease genes and mutations are being described in non-Caucasian inbred populations  and genetic isolates , there is a real need for exome sequencing of healthy individuals within these populations.
Within our cohort, three novel disease genes were initially identified by exome sequencing of single probands (GPR126, KLHL41 and SPEG). The success of disease gene discovery in NEM (KLHL40, KLHL41, LMOD3, MYO18B) and centronuclear myopathies (SPEG) is likely due to the ability to identify patients with a very similar presentation (clinically and based upon very specific muscle biopsy findings) such as to enable screening of candidate genes in patients with the same disease [21, 65].
For fetal hypokinesia and arthrogryposis cases, it is more difficult to deeply phenotype the patients, due in many cases to the poor preservation of fetal tissue and the lack of specific pathological hallmarks from biopsy or autopsy material. A recent study describes exome sequencing of 143 multiplex consanguineous families, in which 33 novel candidate neurogenic disease genes were identified , highlighting the value of studying consanguineous families. As a comparison, only three of the 20 (15 %) genetically-unresolved cases were consanguineous whereas seven of 18 of the genetically-diagnosed cases were consanguineous (39 %, Additional file 1: Table S1), thus one is 2.5-times more likely to identify the causative disease gene in consanguineous families. A genetic diagnosis was achieved in ten of 16 congenital myopathy cases (63 %) and six of 13 arthrogryposis cases (46 %) but only 22 % of fetal akinesia cases (two of nine). It is also possible that the cause of disease, in some of the isolated cases (particularly those diagnosed with fetal akinesia), is not due to a monogenic disorder but may be environmental and/or polygenic. In families with multiple affected siblings and normal CGH arrays, we will pursue whole genome sequencing and/or RNA-seq of target tissue cDNA to try to identify novel disease genes and/or mechanisms.
In summary, this study highlights the use of next generation sequencing to genetically diagnose 47 % of cases within a heterogeneous severe neuromuscular disease cohort. The study has also resulted in the identification of four novel neuromuscular disease genes, and has led to the identification of a novel mechanism of sarcomere assembly and muscle dysfunction involving KLHL40, KLHL41 and LMOD3 [21, 66, 67]. Finally, this study has contributed to extending the phenotypic spectrum of CHRNG, ECEL1, GBE1 and RYR1.
Arthrogryposis multiplex congenita
Fetal akinesia deformation sequence
Lotterywest state biomedical facility genomics node
Neuromuscular sub-exomic sequencing
Hall JG. Analysis of Pena Shokeir phenotype. Am J Med Genet. 1986;25:99–117.
Hall JG. Pena-Shokeir phenotype (fetal akinesia deformation sequence) revisited. Birth Defects Res A Clin Mol Teratol. 2009;85:677–94.
Hall JG. Arthrogryposis multiplex congenita: etiology, genetics, classification, diagnostic approach, and general aspects. J Pediatr Orthop B. 1997;6:159–66.
Quinn CM, Wigglesworth JS, Heckmatt J. Lethal arthrogryposis multiplex congenita: a pathological study of 21 cases. Histopathology. 1991;19:155–62.
Ravenscroft G, Sollis E, Charles AK, North KN, Baynam G, Laing NG. Fetal akinesia: review of the genetics of the neuromuscular causes. J Med Genet. 2011;48:793–801.
Filges I, Hall JG. Failure to identify antenatal multiple congenital contractures and fetal akinesia--proposal of guidelines to improve diagnosis. Prenat Diagn. 2013;33:61–74.
Wilbe M, Ekvall S, Eurenius K, Ericson K, Casar-Borota O, Klar J, et al. MuSK: a new target for lethal fetal akinesia deformation sequence (FADS). J Med Genet. 2015;52(3):195–202.
Haliloglu G, Topaloglu H. Arthrogryposis and fetal hypomobility syndrome. Handb Clin Neurol. 2013;113:1311–9.
Laquerriere A, Maluenda J, Camus A, Fontenas L, Dieterich K, Nolent F, et al. Mutations in CNTNAP1 and ADCY6 are responsible for severe arthrogryposis multiplex congenita with axoglial defects. Hum Mol Genet. 2014;23:2279–89.
Bamshad M, Jorde LB, Carey JC. A revised and extended classification of the distal arthrogryposes. Am J Med Genet. 1996;65:277–81.
Bamshad M, Van Heest AE, Pleasure D. Arthrogryposis: a review and update. J Bone Joint Surg Am. 2009;91 Suppl 4:40–6.
McMillin MJ, Beck AE, Chong JX, Shively KM, Buckingham KJ, Gildersleeve HI, et al. Mutations in PIEZO2 cause Gordon syndrome, Marden-Walker syndrome, and distal arthrogryposis type 5. Am J Hum Genet. 2014;94:734–44.
McMillin MJ, Below JE, Shively KM, Beck AE, Gildersleeve HI, Pinner J, et al. Mutations in ECEL1 cause distal arthrogryposis type 5D. Am J Hum Genet. 2013;92:150–6.
Romero NB, Clarke NF. Congenital myopathies. Handb Clin Neurol. 2013;113:1321–36.
Ravenscroft G, Laing NG, Clarke NF. Chapter 28-Congenital/ultrastructural myopathies. In: Hilton-Jones D, editor. Oxford Textbook of Neuromuscular Disorders. 2014.
Maggi L, Scoto M, Cirak S, Robb SA, Klein A, Lillis S, et al. Congenital myopathies--clinical features and frequency of individual subtypes diagnosed over a 5-year period in the United Kingdom. Neuromuscul Disord. 2013;23:195–205.
Nance JR, Dowling JJ, Gibbs EM, Bonnemann CG. Congenital myopathies: an update. Curr Neurol Neurosci Rep. 2012;12:165–74.
Jain RK, Jayawant S, Squier W, Muntoni F, Sewry CA, Manzur A, et al. Nemaline myopathy with stiffness and hypertonia associated with an ACTA1 mutation. Neurology. 2012;78:1100–3.
Majczenko K, Davidson AE, Camelo-Piragua S, Agrawal PB, Manfready RA, Li X, et al. Dominant mutation of CCDC78 in a unique congenital myopathy with prominent internal nuclei and atypical cores. Am J Hum Genet. 2012;91:365–71.
Agrawal PB, Pierson CR, Joshi M, Liu X, Ravenscroft G, Moghadaszadeh B, et al. SPEG Interacts with Myotubularin, and Its Deficiency Causes Centronuclear Myopathy with Dilated Cardiomyopathy. Am J Hum Genet. 2014;95:218–26.
Yuen M, Sandaradura SA, Dowling JJ, Kostyukova AS, Moroz N, Quinlan KG, et al. Leiomodin-3 dysfunction results in thin filament disorganization and nemaline myopathy. J Clin Invest. 2014;124:4693–708.
Gupta VA, Ravenscroft G, Shaheen R, Todd EJ, Swanson LC, Shiina M, et al. Identification of KLHL41 Mutations Implicates BTB-Kelch-Mediated Ubiquitination as an Alternate Pathway to Myofibrillar Disruption in Nemaline Myopathy. Am J Hum Genet. 2013;93:1108–17.
Ravenscroft G, Thompson EM, Todd EJ, Yau KS, Kresoje N, Sivadorai P, et al. Whole exome sequencing in foetal akinesia expands the genotype-phenotype spectrum of GBE1 glycogen storage disease mutations. Neuromuscul Disord. 2013;23:165–9.
Ravenscroft G, Miyatake S, Lehtokari VL, Todd EJ, Vornanen P, Yau KS, et al. Mutations in KLHL40 Are a Frequent Cause of Severe Autosomal-Recessive Nemaline Myopathy. Am J Hum Genet. 2013;93:6–18.
Kaplan JC. The 2012 version of the gene table of monogenic neuromuscular disorders. Neuromuscul Disord. 2011;21:833–61.
Cabrera-Serrano M, Ghaoui R, Ravenscroft G, Johnsen RD, Davis MR, Corbett A, et al. Expanding the phenotype of GMPPB mutations. Brain. 2015;138:836–44.
Clarke L, Zheng-Bradley X, Smith R, Kulesha E, Xiao C, Toneva I, et al. The 1000 Genomes Project: data management and community access. Nat Methods. 2012;9:459–62.
Flanagan SE, Patch AM, Ellard S. Using SIFT and PolyPhen to predict loss-of-function and gain-of-function mutations. Genet Test Mol Biomarkers. 2010;14:533–7.
Schwarz JM, Cooper DN, Schuelke M, Seelow D. MutationTaster2: mutation prediction for the deep-sequencing age. Nat Methods. 2014;11:361–2.
Ravenscroft G, Nolent F, Rajagopalan S, Meireles AM, Paavola KJ, Gaillard D, et al. Mutations of GPR126 Are Responsible for Severe Arthrogryposis Multiplex Congenita. Am J Hum Genet. 2015;96:955–61.
Flex E, De Luca A, D’Apice MR, Buccino A, Dallapiccola B, Novelli G. Rapid scanning of myotubularin (MTM1) gene by denaturing high-performance liquid chromatography (DHPLC). Neuromuscul Disord. 2002;12:501–5.
Laporte J, Hu LJ, Kretz C, Mandel JL, Kioschis P, Coy JF, et al. A gene mutated in X-linked myotubular myopathy defines a new putative tyrosine phosphatase family conserved in yeast. Nat Genet. 1996;13:175–82.
Monnier N, Marty I, Faure J, Castiglioni C, Desnuelle C, Sacconi S, et al. Null mutations causing depletion of the type 1 ryanodine receptor (RYR1) are commonly associated with recessive structural congenital myopathies with cores. Hum Mutat. 2008;29:670–8.
Kossugue PM, Paim JF, Navarro MM, Silva HC, Pavanello RC, Gurgel-Giannetti J, et al. Central core disease due to recessive mutations in RYR1 gene: is it more common than described? Muscle Nerve. 2007;35:670–4.
Zhou H, Lillis S, Loy RE, Ghassemi F, Rose MR, Norwood F, et al. Multi-minicore disease and atypical periodic paralysis associated with novel mutations in the skeletal muscle ryanodine receptor (RYR1) gene. Neuromuscul Disord. 2010;20:166–73.
von der Hagen M, Kress W, Hahn G, Brocke KS, Mitzscherling P, Huebner A, et al. Novel RYR1 missense mutation causes core rod myopathy. Eur J Neurol. 2008;15:e31–2.
Robinson R, Carpenter D, Shaw MA, Halsall J, Hopkins P. Mutations in RYR1 in malignant hyperthermia and central core disease. Hum Mutat. 2006;27:977–89.
Amburgey K, Bailey A, Hwang JH, Tarnopolsky MA, Bonnemann CG, Medne L, et al. Genotype-phenotype correlations in recessive RYR1-related myopathies. Orphanet J Rare Dis. 2013;8:117.
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:1584–9.
Lehtokari VL, Kiiski K, Sandaradura SA, Laporte J, Repo P, Frey JA, et al. Mutation update: the spectra of nebulin variants and associated myopathies. Hum Mutat. 2014;35:1418–26.
Wallgren-Pettersson C, Donner K, Sewry C, Bijlsma E, Lammens M, Bushby K, et al. Mutations in the nebulin gene can cause severe congenital nemaline myopathy. Neuromuscul Disord. 2002;12:674–9.
Lehtokari VL, Pelin K, Sandbacka M, Ranta S, Donner K, Muntoni F, et al. Identification of 45 novel mutations in the nebulin gene associated with autosomal recessive nemaline myopathy. Hum Mutat. 2006;27:946–56.
Lehtokari VL, Greenleaf RS, DeChene ET, Kellinsalmi M, Pelin K, Laing NG, et al. The exon 55 deletion in the nebulin gene--one single founder mutation with world-wide occurrence. Neuromuscul Disord. 2009;19:179–81.
Anderson SL, Ekstein J, Donnelly MC, Keefe EM, Toto NR, LeVoci LA, et al. Nemaline myopathy in the Ashkenazi Jewish population is caused by a deletion in the nebulin gene. Hum Genet. 2004;115:185–90.
Kiiski K, Laari L, Lehtokari VL, Lunkka-Hytonen M, Angelini C, Petty R, et al. Targeted array comparative genomic hybridization--a new diagnostic tool for the detection of large copy number variations in nemaline myopathy-causing genes. Neuromuscul Disord. 2013;23:56–65.
Kiiski, K., Lehtokari, V.L., Loytynoja, A., Ahlsten, L., Laitila, J., Wallgren-Pettersson, C., and Pelin, K. (2015). A recurrent copy number variation of the NEB triplicate region: only revealed by the targeted nemaline myopathy CGH array. Eur J Hum Genet. Jul 22. doi:10.1038/ejhg.2015.166. [Epub ahead of print].
Brownlow S, Webster R, Croxen R, Brydson M, Neville B, Lin JP, et al. Acetylcholine receptor delta subunit mutations underlie a fast-channel myasthenic syndrome and arthrogryposis multiplex congenita. J Clin Invest. 2001;108:125–30.
Gomez CM, Maselli RA, Vohra BP, Navedo M, Stiles JR, Charnet P, et al. Novel delta subunit mutation in slow-channel syndrome causes severe weakness by novel mechanisms. Ann Neurol. 2002;51:102–12.
Muller JS, Baumeister SK, Schara U, Cossins J, Krause S, von der Hagen M, et al. CHRND mutation causes a congenital myasthenic syndrome by impairing co-clustering of the acetylcholine receptor with rapsyn. Brain. 2006;129:2784–93.
Michalk A, Stricker S, Becker J, Rupps R, Pantzar T, Miertus J, et al. Acetylcholine receptor pathway mutations explain various fetal akinesia deformation sequence disorders. Am J Hum Genet. 2008;82:464–76.
Shen XM, Ohno K, Fukudome T, Tsujino A, Brengman JM, De Vivo DC, et al. Congenital myasthenic syndrome caused by low-expressor fast-channel AChR delta subunit mutation. Neurology. 2002;59:1881–8.
Hall JG. Pretibial linear vertical creases or indentations (shin dimples) associated with arthrogryposis. Am J Med Genet A. 2013;161A:737–44.
Morgan NV, Brueton LA, Cox P, Greally MT, Tolmie J, Pasha S, et al. Mutations in the embryonal subunit of the acetylcholine receptor (CHRNG) cause lethal and Escobar variants of multiple pterygium syndrome. Am J Hum Genet. 2006;79:390–5.
Robinson KG, Viereck MJ, Margiotta MV, Gripp KW, Abdul-Rahman OA, Akins RE. Neuromotor synapses in Escobar syndrome. Am J Med Genet A. 2013;161A:3042–8.
Toydemir RM, Rutherford A, Whitby FG, Jorde LB, Carey JC, Bamshad MJ. Mutations in embryonic myosin heavy chain (MYH3) cause Freeman-Sheldon syndrome and Sheldon-Hall syndrome. Nat Genet. 2006;38:561–5.
Barnett CP, Todd EJ, Ong R, Davis MR, Atkinson V, Allcock R, et al. Distal arthrogryposis type 5D with novel clinical features and compound heterozygous mutations in ECEL1. Am J Med Genet A. 2014;164A:1846–9.
Dieterich K, Quijano-Roy S, Monnier N, Zhou J, Faure J, Smirnow DA, et al. The neuronal endopeptidase ECEL1 is associated with a distinct form of recessive distal arthrogryposis. Hum Mol Genet. 2013;22:1483–92.
Monk KR, Oshima K, Jors S, Heller S, Talbot WS. Gpr126 is essential for peripheral nerve development and myelination in mammals. Development. 2011;138:2673–80.
Chong JX, McMillin MJ, Shively KM, Beck AE, Marvin CT, Armenteros JR, et al. De Novo Mutations in NALCN Cause a Syndrome Characterized by Congenital Contractures of the Limbs and Face, Hypotonia, and Developmental Delay. Am J Hum Genet. 2015;96:462–73.
Rossor AM, Oates EC, Salter HK, Liu Y, Murphy SM, Schule R, et al. Phenotypic and molecular insights into spinal muscular atrophy due to mutations in BICD2. Brain. 2015;138:293–310.
Xue, Y., Ankala, A., Wilcox, W.R., and Hegde, M.R. (2014). Solving the molecular diagnostic testing conundrum for Mendelian disorders in the era of next-generation sequencing: single-gene, gene panel, or exome/genome sequencing. Genet Med. Jun;17 (6):444-51. doi:10.1038/gim.2014.122. Epub 2014 Sep 18.
Ankala A, Da Silva C, Gualandi F, Ferlini A, Bean LJ, Collins C, et al. A comprehensive genomic approach for neuromuscular diseases gives a high diagnostic yield. Ann Neurol. 2015;77:206–14.
Alazami AM, Patel N, Shamseldin HE, Anazi S, Al-Dosari MS, Alzahrani F, et al. Accelerating novel candidate gene discovery in neurogenetic disorders via whole-exome sequencing of prescreened multiplex consanguineous families. Cell Rep. 2015;10:148–61.
Tan-Sindhunata, M.B., Mathijssen, I.B., Smit, M., Baas, F., de Vries, J.I., van der Voorn, J.P., Kluijt, I., Hagen, M.A., Blom, E.W., Sistermans, E., et al. (2014). Identification of a Dutch founder mutation in MUSK causing fetal akinesia deformation sequence. Eur J Hum Genet. Sep;23 (9):1151-7. doi:10.1038/ejhg.2014.273. Epub 2014 Dec 24.
Alazami AM, Kentab AY, Faqeih E, Mohamed JY, Alkhalidi H, Hijazi H, et al. A novel syndrome of Klippel-Feil anomaly, myopathy, and characteristic facies is linked to a null mutation in MYO18B. J Med Genet. 2015;52:400–4.
Garg A, O’Rourke J, Long C, Doering J, Ravenscroft G, Bezprozvannaya S, et al. KLHL40 deficiency destabilizes thin filament proteins and promotes nemaline myopathy. J Clin Invest. 2014;124:3529–39.
Cenik BK, Garg A, McAnally JR, Shelton JM, Richardson JA, Bassel-Duby R, et al. Severe myopathy in mice lacking the MEF2/SRF-dependent gene leiomodin-3. J Clin Invest. 2015;125:1569–78.
We are grateful to the patients and their families for their participation in this research. This research was supported by the National Health and Medical Research Council of Australia (Early Career Researcher Fellowship #1035955 to GR, Research Fellowship APP1002147 to NGL and Project Grant APP1022707; EU Collaborative grant APP1055295); the Association Francaise contre les Myopathies (#15734) and a UWA Collaborative Research Award. EJT and KSY were supported by University Postgraduate Awards.
The authors declare that they have no competing interests.
GR and NGL conceived and coordinated the project. EJT, KSY, RO, PS, RJNA, SM, DB, MRD and GR performed the experiments and analysed data. JS, GM, CPB, GH, ZA, AK, AC, NFC, M-LF, DW, CAS, AC, MMR and CK-B cared for the patients and co-ordinated patient enrolment in the study. BT, NBR and CAS performed histology and immunohistochemistry. EJT, CAS, MMR, NGL and GR drafted the manuscript. All authors read and approved the manuscript.
For each family the affected individuals sequenced are shown as well as the type of NGS that was performed for each patient. The disease group and consanguinity for each family is also indicated as are the coverage statistics for each sample and NGS. (DOCX 21 kb)