Identification of a Dutch founder mutation in MUSK causing fetal akinesia deformation sequence.

Fetal akinesia deformation sequence (FADS) refers to a clinically and genetically heterogeneous group of disorders with congenital malformations related to impaired fetal movement. FADS can result from mutations in CHRNG, CHRNA1, CHRND, DOK7 and RAPSN; however, these genes only account for a minority of cases. Here we identify MUSK as a novel cause of lethal FADS. Fourteen affected fetuses from a Dutch genetic isolate were traced back to common ancestors 11 generations ago. Homozygosity mapping in two fetuses revealed MUSK as a candidate gene. All tested cases carried an identical homozygous variant c.1724T>C; p.(Ile575Thr) in the intracellular domain of MUSK. The carrier frequency in the genetic isolate was 8%, exclusively found in heterozygous carriers. Consistent with the established role of MUSK as a tyrosine kinase that orchestrates neuromuscular synaptogenesis, the fetal myopathy was accompanied by impaired acetylcholine receptor clustering and reduced tyrosine kinase activity at motor nerve endings. A functional assay in myocytes derived from human fetuses confirmed that the variant blocks MUSK-dependent motor endplate formation. Taken together, the results strongly support a causal role of this founder mutation in MUSK, further expanding the gene set associated with FADS and offering new opportunities for prenatal genetic testing.

Introduction

Fetal akinesia deformation sequence (FADS; MIM 208150) represents a group of disorders characterized by decreased or absent fetal movements. Affected infants may die in utero or shortly after birth due to severe pulmonary hypoplasia. FADS was first recognized in 1974 by Pena and Shokeir, and therefore was originally called the Pena–Shokeir syndrome.[1, 2] The etiology of FADS is heterogeneous, including neuromuscular disorders, metabolic diseases, syndromic and environmental causes.[3] Mutations in genes expressed at the neuromuscular junction (NMJ) are increasingly recognized as important causes of FADS, including CHRNA1 (MIM 100690), CHRND (MIM 100720), CHRNG (MIM 100730), DOK7 (MIM 610285) and RAPSN (MIM 601592).[3] These genes however explain a minority of cases: in a cohort of 59 kindreds with Lethal Multiple Pterygium Syndrome (LMPS [MIM 253290]) and FADS, it was estimated that 24% were positive for mutations in either of these genes.[4]

Here we describe a novel founder mutation in MUSK, encoding a muscle-specific tyrosine kinase essential for NMJ development[5] as a frequent cause of lethal FADS in a Dutch genetic isolate. We demonstrate that this variant causes a defect in neuromuscular synaptogenesis.

Materials and methods

Informed consent was given by all subjects included in this study. For homozygosity mapping in two fetuses affected with FADS (case 10 and 11 in Table 1), an Affymetrix 6.0 array was used. Overlapping regions of homozygosity were analyzed for candidate genes (‘Genomic Oligoarray and SNP array evaluation tool v1.0').[6] Sanger sequencing of the coding region and flanking intronic regions of MUSK (primer details available on request) was performed on a ABI 3730 type DNA analyzer (Applied Biosystems, Foster City, CA, USA). The variant was deposited in the Leiden Open Variation Database at www.lovd.nl (patient ID no. 0033104). After routine autopsy, various fetal muscle tissues (m. quadriceps, m. psoas, m. biceps femoris, diaphragm) were collected for vital and nonvital assays, including hematoxylin and eosin staining, elastic Van Gieson staining and ATPase assays to distinguish type I and type II muscle fibers. To visualize NMJs, cryosections were fixed with 4% paraformaldehyde and stained for synaptotagmin-1 (rabbit polyclonal W855, 1:500), phosphotyrosine (mouse monoclonal 4G10, 1:200, Millipore, Amsterdam-Zuidoost, Netherlands)[7] and acetylcholine receptors (α-bungarotoxin-Rhodamine, 1:1000, Sigma, St Louis, MO, USA). For a side-by-side comparison of wild-type and mutant motor endplates of similar developmental age, stainings of case 8 (23 weeks of gestation) were performed in parallel with that of a control fetus aborted due to cervical insufficiency (22 weeks of gestation).

Table 1

Clinical features of 14 fetuses with FADS

Fetus number	1 (F)	2 (F) Sib of 3	3 (F) Sib of 2	4 (M)	5 (M)	6 (F)	7 (M)	8 (M)	9 (M)	10 (M) Sib of 11	11 (F) Sib of 10	12 (F) Sib of 13	13 (F) Sib of 12	14 (M)
Molecular analysis	I575T/I575T	I575T/I575T	NA	I575T/I575T	I575T/I575T	I575T/I575T	I575T/I575T	I575T/I575T	NA	I575T/I575T	I575T/I575T	I575T/I575T	NA	I575T/I575T
Gestational age at birth	41w	33w 4d	34w	38w 1d	38w 4d	33w 1d	35w	23w	37w	32w 2d	22w	31w	35w 6d	23w
Pregnancy outcome	Died <2 h pp	Died <1 h pp	Died <1 h pp	Died 5d pp	Died <1 h pp	Stillborn	Died <1 h pp	TOP	Died <1 h pp	Died <1 h pp	TOP	Died 2 h pp	Died 4d pp	TOP
Prenatal findings
Polyhydramnios	NA	+	+	From 30w	+	+	+	−	From 31w	From 22w	−	+	From 31w	+
Fetal growth restriction	NA	+	NA	−	−	−	−	−	−	+	−	−	−	−
Nuchal translucency/cystic hygroma	NA	NA	NA	NA	NA	−	NA	+	−	−	−	NA	NA	+
Generalized edema/hydrops	NA	NA	NA	NA	NA	−	NA	−	−	−	−	NA	NA	+
Reduced/absent fetal movements	NA	NA	+	NA	−	+	NA	+	+	+	+	NA	NA	+
Joint contractures	NA	NA	NA	NA	NA	+	NA	+	NA	+	+	NA	NA	+
Club feet	NA	NA	NA	NA	NA	c)	NA	−	NA	+/+	+/+	NA	NA	+/+
Postnatal findings
Birth weight (g)	2700 (<p3)	1870 (<p3)	NA	3470 (p75)	3250 (p50)	1363 (p3-p10)	1910 (p10)	551 (p50)	2500 (p10)	1610 (p25-50)	408 (p10)	1250 (p10-p25)	2210 (p10)	490 (p25)
Head circumference (cm)	36 (p50)	NA	NA	36.5 (p97)	NA	28 (p10)	32.5 (p25-50)	22.5 (p75)	35.5 (p75-90)	30 (p50)	19.6 (p10)	NA	36 (p97)	20.8 (p50)
Hypertelorism	NA	+	+	+	−	+	NA	+	+	+	+	−	−	+
Low-set ears	NA	+	+	−	−	+	NA	+	+	+	+	+	−	−
Micrognathia	+	−	NA	+	−	+	−	+	+	+	+	−	−	+
Small mouth	NA	NA	NA	+	−	−	−	+	−	−	−	−	−	+
Cleft palate	−	−	NA	−	−	−	−	−	−	−	−	−	−	−
Pterygia/webbing	NA	−	NA	−	−	−	−	−	−	−	−	−	−	−
Reduced muscle bulk	NA	NA	NA	−	−	+	NA	+	+	+	+	−	−	+
Oedema/hydrops	NA	−	NA	−	+	−	−	+	−	−	−	−	−	+
Pulmonary hypoplasia	+	+	+	+	+	+	+	+	+	+	+	NA	+	+
Undescended testis	NR	NR	NR	+/−	+/+	NR	+/+	+/+	+/+	+/+	NR	NR	NR	+/+
Joint contractures
Shoulders	NA	NA	NA	−	−	+	NA	+	−	NA	NA	NA	−	+
Elbows	NA	+	+	−	−	+	−	+	−	+	−	NA	−	−
Wrists	NA	+	NA	−	−	+	−	+	−	+	+	NA	−	−
Fingers	NA	+	+	+	+	+	+	+	+	+	+	NA	+	+
Hips	+	+	+	−	+	+	−	+	+	+	−	NA	+	+
Knees	NA	+	+	−	+	+	−	+	+	+	−	NA	−	−
Club feet	+/+	+/−	+	−	+/+	−	+/+	−	NA	+/+	+/+	−	−	+/+
Rockerbottom feet	NA	−	−	−	−	−	−	+/+	NA	−	−	−	−	−
Other		a,e		b		c		d	a	f	g

Abbreviations: F, female; M, male; pp, postpartum; NA, not assessed; NR, not relevant; TOP, termination of pregnancy; +, present; −, absent; +/+, bilateral; +/−, unilateral.

All cases were from a Dutch ethnic origin.

(a)Gracile ribs; (b)hemivertebrae Th4, convulsions; (c)plantar flexion of the feet; (d)ventriculomegaly, choroid plexus cysts; (e)large liver and adrenals; (f)unilateral syndactyly digits 1–2 and agenesis digit 3 of the foot; (g)bilateral cervical ribs.

To assay MUSK-induced motor endplate maturation, fetal myocytes from a fetus homozygous for the c.1724T>C; p.(Ile575Thr) variant were maintained in DMEM supplemented with 20% fetal bovine serum, glutaMAX and penicillin/streptomycin at 37 °C in 5% CO2. Myotubes were induced by growing the myocytes in 2% serum on coverslips coated with 72 μg/ml Matrigel (BD Biosciences). After 1 day in vitro (DIV1), myocytes were transduced with lentiviral vectors. Wild-type or mutant MUSK was engineered to contain a Flag tag in the ectodomain, a strategy previously shown to preserve MUSK activity.[8] For efficient and safe cellular delivery, a third-generation lentiviral expression system was used.[9] At DIV6, acetylcholine receptor (AChR) cluster formation was induced by extracellular administration of 100 ng/ml neuronal agrin (splice variant 3,4,8; R&D Systems, Minneapolis, MN, USA).[10] At DIV8, the myocytes were fixed with 4% paraformaldehyde, washed in Tris-buffered saline containing 0.1% Tween-20 (TBST) and stained with DAPI (1:2000) and α-bungarotoxin-Rhodamine (1:2000) in TBST for 48 h at 4 °C. The cells were washed again in TBST and mounted in Mowiol, followed by image acquisition on a custom-build FAINT microscope (Chromaphor, Oberhausen, Germany). AChR clusters were counted in n=35 images from three coverslips per condition, and a Student's t-test was used to compare the mean cluster density.

Results

Homozygosity mapping in two fetuses with FADS from a single consanguineous family (case 10 and 11 in Table 1 and Figure 1a) identified MUSK as the prime candidate gene. Sanger sequencing revealed a novel homozygous missense variant (NM_005592.3: c.1724T>C; p.(Ile575Thr); Figure 1b) absent in relevant databases (dbSNP build 132, 1000 Genomes project, and Exome Variant Server data ESP5400). The amino-acid change was predicted to be potentially pathogenic by various bioinformatics tools: SIFT: Deleterious (score: 0, median: 4.32); Align GVGD: C65 (GV: 0.00–GD: 89.28); PhyloP: 4.97 (highly conserved); Mutation Taster: disease causing, prob: 0.9999999987; Polyphen: probably damaging (HumVar score 0.999). The variant was also absent in a cohort of 1326 Dutch controls, analyzed by Taqman SNP Genotyping (Life Technologies).

Figure 1

Co-segregation of a variant in MUSK with lethal FADS in a Dutch genetic isolate. (a) Schematic illustrating the common ancestry of 14 affected fetuses from 11 families. (b) Sanger sequencing chromatograms of fetal (top) and parental (bottom) DNA samples identifying the variant NM_005592.3: c.1724T>C; p.(Ile575Thr) and multispecies alignment of the affected residue (right). (c) Domain structure of MUSK depicting the extracellular domain containing IgG-like domains, a Frizzled domain (FZ, also named C6 box), a transmembrane domain (TM) and a cytoplasmic domain that contains a juxtamembrane (JM)/PTB-domain and the tyrosine kinase domain. Ile575 is located at the boundary between the JM and kinase domains, very close to the Dok7-interacting region and the autophosphorylation site Tyr554. Other mutations previously reported to cause CMS are also indicated. Mutations found in compound heterozygotes are color-matched[23, 24, 25, 26, 27, 28]. All amino-acid numbers refer to MUSK isoform 1 (GenPept accession NP_005583.1).

As the affected fetuses originated from a genetic isolate, we hypothesized the occurrence of a founder effect. Twelve other cases from 10 families from the same genetic isolate came to our clinic between 1980 and 2012. DNA was available from nine cases. All carried the same homozygous variant and all parents were heterozygous carriers. All 14 cases had common ancestors traced back 11 generations (Figure 1a). Consistent with a founder effect, the MUSK variant co-segregated with a shared haplotype around the MUSK gene (Supplementary Table 1). Genetic testing of 530 healthy controls from the genetic isolate revealed no homozygous individuals and 43 heterozygous carriers, resulting in a carrier frequency of 8% in this region. Subsequently, a cohort of 26 unresolved LMPS/FADS cases who did not come from the genetic isolate were tested by Sanger sequencing, revealing no likely pathogenic variants in MUSK.

The clinical symptoms available from 14 affected fetuses are summarized in Table 1. Pictures of four fetuses are shown in Figure 2. The pregnancy was terminated in three fetuses, whereas the natural course of the pregnancy was followed in 11 fetuses. From the naturally born fetuses, seven were born prematurely between 31–37 weeks of gestation and four were born at term. All 11 were stillborn or died shortly after birth, usually after failed resuscitation. Typical prenatal ultrasound findings were polyhydramnios, progressively reduced fetal movements and joint contractures. The polyhydramnios generally developed in the second or third trimesters, and became severely pronounced around 30 weeks, sometimes leading to premature rupture of membranes or premature contractions. Post-mortem examination revealed severe pulmonary hypoplasia in all cases. Hypertelorism, low-set ears, micrognathia and reduced muscle bulk were seen in most fetuses. The contractures could affect all joints with variable severity; however, all fetuses had contractures of the fingers (camptodactyly) and 10 out of 13 analyzed fetuses had severely affected feet (club feet, plantar flexion or rocker-bottom feet). Two fetuses who died 4 and 5 days after birth were severely hypotonic. All male fetuses had cryptorchidism, although this may not be representative because most children were born prematurely. Three of twelve analyzed fetuses were mildly edematous. Gracile ribs were observed in two cases. Some additional findings were observed only in a single fetus and are considered to be incidental: cervical ribs, mild ventriculomegaly, choroid plexus cysts, a solitary hemivertebra, syndactyly of the first and second toes and agenesis of the third toe. None of the fetuses had pterygia or a cleft palate.

Figure 2

Macroscopic appearance of four fetuses homozygous for the c.1724T>C; p.(Ile575Thr) variant in MUSK. (a) Case 6, stillborn at a gestational age of 33 weeks; (b) case 8, termination of pregnancy at 23 weeks; (c) case 10, died immediately after birth at 32 weeks; (d) case 11, termination of pregnancy at 22 weeks. Note the multiple contractures, facial appearance of hypertelorism and micro-retrognathia, club feet (b–d), oligodactyly of the right foot (c), plantar flexion of the feet (a). For details see Table 1. Bars indicate 5 cm.

Microscopic examination of tissues from case 8 (a terminated pregnancy at 23 weeks of gestation) showed no abnormalities of the cerebrum, cerebellum, cranial nerve motor nuclei and spinal motor neurons. In contrast, muscle tissues (m. quadriceps, m. psoas, m. biceps femoris, diaphragm) showed several myopathic features consistent with a denervating or neurogenic cause (Figure 3). There was clear variation in muscle fiber diameter with small and rounded shaped atrophic fibers and there was a mild increase in endomysial fibrosis. The number of internal nuclei was increased consistent with maturational delay. The ratio between type 1 and 2 fibers was 1:10 suggesting a loss of type 1 (slow) fibers (Figure 3b and c).

Figure 3

Histopathological findings indicative for defects in neuromuscular synaptogenesis. (a) Skeletal muscle biopsies showed small and rounded atrophic fibers and an increased number of intracellular nuclei (arrowheads). (b and c) Immunostaining of type I (b) and type II fibers (c) indicated a relative loss of type I fibers. (d–g) Combined staining for nerve terminals (synaptotagmin-1) and AChR clusters (α-bungarotoxin-rhodamine) in quadriceps muscle indicated fewer and smaller neuromuscular junctions in case 8 (arrowheads in d, e; 23 weeks of gestation) compared with a control fetus (f, g; 22 weeks). (h–k) Combined staining for nerve terminals, AChR clusters and anti-phosphotyrosine (pTyr) monoclonal 4G10 in diaphragm indicated fewer clusters (arrowheads) and a loss of tyrosine kinase activity in MUSKI575T/I575T endplates (h, i) compared with control tissue (j, k). Note that residual tyrosine kinase activity was detected in some mutant endplates (arrow), which may also result from unrelated tyrosine kinases.

MUSK is a receptor tyrosine kinase expressed on the muscle cell membrane.[5, 11] It consists of an extracellular domain containing IgG-like domains, a Frizzled-like cysteine-rich domain (C6 box), a transmembrane domain and a cytoplasmic domain composed of a juxtamembrane/PTB domain and a tyrosine kinase domain (Figure 1c).[12, 13] The variant c.1724T>C; p.(Ile575Thr) is located in a highly conserved region at the boundary of the juxtamembrane and kinase domains. The important role of both domains in tyrosine kinase activation and Dok-7 signaling suggests that the variant may affect MUSK signaling and AChR clustering.[14, 15] To investigate this, local tyrosine kinase activity and AChR clusters were visualized at nascent motor endplates in muscle biopsies from an affected fetus (case 8, 23 weeks of gestation) and a control fetus aborted at 22 weeks due to cervical insufficiency (Figure 3d–k). Ingrowing axons from motor neurons were visualized with the synaptic vesicle marker Synaptotagmin-1. In wild-type tissue, developing motor endplates were closely associated with AChR clusters, stained by α-bungarotoxin-Rhodamine. In MUSK mutant muscle, these clusters were strongly reduced in number and appeared very small (Figure 3d, e and h). Whereas motor endplates in the control fetus showed increased tyrosine kinase activity (Figure 3j and k), the small AChR clusters in the FADS patient did not accumulate phosphotyrosines (Figure 3h and i). Thus, histological examination was consistent with a defect in the MUSK signaling pathway.

To corroborate the above findings in fetal muscle tissue, we next investigated the AChR-clustering activity of mutant and wild-type MUSK in a direct comparison. Wild-type and mutant MUSK were expressed in HEK293 cells with similar molecular mass (Figure 4c), migrating in SDS-PAGE as a 110-kDa protein consistent with previous studies.[16, 17] In myocytes derived from a fetus homozygous for the c.1724T>C variant, overexpression of wild-type MUSK significantly increased the density of agrin-induced AChR clusters (P<0.001, n=35), whereas mutant MUSKI575T did not (P=0.65, n=35; see Figure 4a and b). Taken together, our observations demonstrate that the c.1724T>C; p.(Ile575Thr) variant strongly affects MUSK autoactivation, tyrosine phosphorylation and AChR clustering in both fetal tissue and cultured myocytes.

Figure 4

Wild-type but not mutant MUSK augments AChR clustering in MUSKI575T myocytes. (a) Primary myocytes derived from an affected fetus homozygous for MUSKI575T were plated on Matrigel-coated coverslips and transduced with lentiviral vectors encoding either MUSKWT or MUSKI575T. As a control, non-transduced MUSKI575T cells were used. Myocytes were treated with neuronal Agrin(3,4,8) and stained for nuclei (DAPI, displayed in green) and AChR clusters (Rhodamine-bungarotoxin, red). Representative examples from each experimental group are shown. (b) Mean±SEM values for the AChR cluster density, quantitated from 35 images from three coverslips of each group. ***P<0.001 (Student's t-test). (c) Immunoblot analysis of Flag-tagged MUSKWT and MUSKI575T constructs expressed in HEK293 cell lysates, suggesting similar yield and molecular mass.

Discussion

Mutations in MUSK, RAPSN, CHRNA1, CHRNB1, CHRND, CHRNE, CHRNG and DOK7 can also cause Congenital Myasthenic Syndromes (CMSs; MIM 608931, 608930, 601462).[18] The severity of the mutation likely determines the severity of the phenotype ranging from CMS to lethal FADS.[19, 20, 21, 22] For CHRNG, a recent study has shown intra- and interfamilial variability for the same mutation and the severity of the phenotype (LMPS/FADS versus the nonlethal Escobar Variant of Multiple Pterygium Syndrome; MIM 265000), suggesting that there may be genetic and/or environmental modifiers.[4]

Previously, recessive mutations in MUSK were found in 11 individuals from five families with CMS. Affected subjects in the first family were compound heterozygous for a frameshift (p.Arg74fs) and missense mutation (p.Val790Met).[23] In the second consanguineous family a homozygous missense mutation was found (p.Pro344Arg in the cysteine-rich extracellular domain).[24] The third family included an individual with heteroallelic missense mutations (p.Met605Ile and p.Ala727Val) both located in the kinase domain.[25] An Iranian case from consanguineous parents was homozygous for a missense mutation (p.Met835Val).[26] Finally, two brothers from Turkish origin with a salbutamol-responsive limb-girdle CMS were compound heterozygous for a missense mutation (p.Asp38Glu) and a genomic deletion of exons 2–3.[27] The clinical symptoms of these 11 CMS cases consisted of a variable, fatigable, predominantly proximal muscle weakness, partial ophthalmoparesis, mild ptosis and facial weakness. None of the patients had congenital joint contractures. Although there could be severe respiratory failure and hypotonia in the neonatal period, the condition seemed to improve during childhood and there was no delay of initial motor milestones.

Comparing the clinical and histological findings, the c.1724T>C; p.(Ile575Thr) variant appears to halt neuromuscular synaptogenesis at an early stage, and cause a more severe loss of function than the previously observed mutations, which cause the milder CMS phenotype. The relatively severe effect of the mutation may be explained by its location very close to the Dok7 interaction site, critically involved in MUSK dimerization and its tyrosine kinase activation.[15] Consistent with this hypothesis, the clinical phenotype of the MUSKI575T/I575T fetuses is reminiscent of that of Musk knockout mice that fail to form stable NMJs, resulting in paralyzed embryos that die perinatally before taking a single breath.[5]

In conclusion, the novel founder mutation in MUSK is a major cause of lethal FADS in the studied genetic isolate and caused early defects in motor endplate development including decreased tyrosine kinase activity and a lack of AChR cluster formation. Although the absence of pterygia is noteworthy, this feature is shared with RAPSN- and DOK7-associated FADS phenotypes.[4] As clinical heterogeneity further obscures genotype–phenotype correlations, a molecular genetic analysis of the gene set involved in FADS/LMPS is probably the most effective classification method for this disease spectrum. The addition of MUSK to this gene set strengthens the existing idea that genes encoding neuromuscular synaptogenesis can cause phenotypes with variable severity, ranging from late-onset and neonatal CMS to the peri- or prenatal lethality observed in fetal akinesia.