The Impairment of MAGMAS Function in Human Is Responsible for a Severe Skeletal Dysplasia

Summary

Skeletal dysplasias (SD) refer to a complex group of rare genetic disorders affecting the growth and development of the skeleton. The identification of the molecular basis of a considerable number of SD has greatly expanded our knowledge of the ossification process. Among SD, spondylodysplastic dysplasia is a generic term describing different conditions characterized by severe vertebral abnormalities and distinct by additional specific features. Several entities within this group are well defined. However, a few cases remain unclassified, of which a novel autosomal recessive spondylometaphyseal dysplasia recently reported by Mégarbané et al. in two Lebanese families. Here, we identified MAGMAS as a candidate gene responsible for this severe SD. MAGMAS, also referred to as PAM16, is a mitochondria-associated protein, involved in pre-proteins import into mitochondria and essential for cell growth and development. We demonstrated that MAGMAS is expressed in bone and cartilage in early developmental stages underlining its specific role in skeletogenesis. We also give strong evidence of the deleterious effect of the identified mutation on the in-vivo activity of Magmas and the viability of yeast strains. Reporting deleterious MAGMAS mutation in a SD supports a key and specific role for this mitochondrial protein in ossification.

Introduction

During embryonic development and postnatal growth period, most longitudinal bone growth occurs via endochondral ossification, a process in which chondrocytes produce a cartilage template that will ultimately be replaced by bone [1]. Disturbance of this tightly regulated process results in skeletal dysplasias (SD), which are an extremely diverse and complex group of rare genetic diseases [2]. Causative mutations have now been identified in over 230 different genes in more than 400 unique skeletal phenotypes [3], [4]. However, the genetic basis of over 100 different entities remains to be determined. Among SD, spondylodysplastic dysplasia is a generic descriptive term for a heterogeneous group of SD characterized by severe vertebral abnormalities and distinctive skeletal or extra-skeletal features [3], [4]. The outcome is variable from early lethality to long-term survival. The last international classification recognizes at least five distinct entities within this group of rare disorders, i) achondrogenesis type 1A [MIM 200600] due to mutations in TRIP11 [MIM 604505], which encodes the thyroid hormone receptor interactor 11, ii) Schneckenbecken dysplasia [MIM 296250] due to mutations in SLC35D1 [MIM 610804], encoding the solute carrier family 35 member D1, also known as UDP-Galactose transporter-related 7 iii) spondylometaphyseal dysplasia (SMD) Sedaghatian type [MIM 250220], a heterogeneous condition that is due in some cases, to mutations in SBDS [MIM 607444] involved in the processing of ribosomal RNA, iv) fibrochondrogenesis [MIM 228520], due in some cases to COL11A1 mutations and v) opsismodysplasia [MIM 258480], due to mutations in INPPL1 [MIM 600829] that encodes the inositol polyphosphate phosphatase like-1 [5]. The wide range of proteins involved in this group of SD points out the complexity of the ossification process.

Although different entities within the group of spondylodysplastic dysplasia are well defined, a few cases remain unclassified. Recently, Mégarbané et al. have reported on two unrelated consanguineous Lebanese families with four affected cases presenting a novel type of early lethal spondylometaphyseal dysplasia [6], [7]. To elucidate the molecular basis of this entity, we performed exome sequencing in these two distinct families and identified a homozygous mutation in MAGMAS (Mitochondria-associated granulocyte macrophage colony stimulating factor-signaling molecule) in patients from both families.

Results

Patients

Two unrelated consanguineous Lebanese families with four affected cases presenting a novel type of early lethal spondylodysplastic dysplasia, recently reported by Mégarbané et al. [6], [7], were included in this study. The main clinical features were pre- and postnatal growth retardation, developmental delay, dysmorphic features, prominent abdomen, narrow thorax with short ribs and short limbs. Skeletal features included severe platyspondyly at birth that improved somewhat over time, square iliac bones, horizontal acetabulae, short long bones with abnormal modeling, widening of the distal femoral metaphyses and delayed epiphyseal ossification (Figure 1). In the first reported family (Family F1, Figure 2) [6], both affected sibs died at an early age from respiratory insufficiency (F1-IV.1 at age 9 months and F1-IV.3 at age 2 years) while heart failure was the cause of the early death, at age 2 years, in patients (F2-IV.2 and F2-IV.3) of the second family (Family F2: individuals IV.2 and IV.3) [7] (Figure 2).

Molecular findings

Exome sequencing has been performed in Family 2, in three individuals: the proband (F2-IV.3) and her two unaffected parents (F2-III.1 and F2-III.2) (Figure 2). Considering the recessive inheritance of the disease in the pedigree and consanguinity (F = 1/16), we assumed identity by Descent [8]. Using an in-house bioinformatics tool, we realized identity by descent filtering of exome data and identified 1037 homozygous variants in the proband, which were also heterozygous in both parents. Of these 1037 variants, 916 were exonic and/or splice site mutations, of which 434 were either non-synonymous and nonsense single nucleotide variants, or small insertions/deletions leading to a shift in the reading frame.

In order to further refine the list of candidate nucleotide variations, additional filtering was performed on what we called the “dbSNP13700” database (i.e. dbSNP137, from which we have removed mutations that are reported as pathogenic alleles, as well as all referenced variants with frequencies lower than 1%, or without any frequency information). Additional filtering was also performed using our in-house exome database that includes all homozygous variants detected by exome sequencing of 14 Lebanese unaffected individuals. We ended up with a list of 16 candidate homozygous by descent variants in the proband F2-IV.3.

In family 1, exome sequencing has been performed only in the proband (F1-IV.3, Figure 2). Using the same analysis pipeline and after all filtering steps, we were able to identify 59 homozygous variants in the genes coding regions of this patient.

In total, patients F1-IV.3 and F2-IV.3 shared 13872 homozygous variants, of which only 3 exonic variations were both homozygous by descent in the proband from family 2 (F2-IV.3), and homozygous in the proband from family 1 (F1-IV.3), after applying the different filters (Table S1). Of those three variants, only one missense mutation (NM_016069: c.226A>G; p.Asn76Asp) in the PAM16 gene (presequence translocase-associated motor 16), that is more commonly referred to as MAGMAS (Mitochondria-associated granulocyte macrophage CSF-signaling molecule), was predicted to have pathogenic effects using different bioinformatics prediction softwares, such as: Mutation Taster (disease causing variation with a probability value of 0.99999957), SIFT (damaging effect, score 0.006) [9], PROVEAN (deleterious effect, score −4.021) [10] and PolyPhen-2 (possibly damaging effect, score 0.728 HumDiv). By Sanger sequencing, we confirmed the segregation of the c.226A>G transition in MAGMAS with the disease in both families (Figure 3A): the mutation was homozygous in the patients, heterozygous in the parents and in the unaffected sibs in both families. We also demonstrated that the variation was absent in a set of 550 Lebanese control chromosomes. Protein sequence alignments showed that the Asparagine at codon 76 in MAGMAS is well conserved among mammals and vertebrates (Figure 3B).

The two patients described here are not related. However, they display the same homozygous missense mutation, thus suggesting the existence of a founder mutation. In order to verify this hypothesis, we used both exome data and genotyping of STR markers. By analyzing exome data, we were able to identify 11 regions of shared homozygosity on chromosome 16. However, only three regions displayed a size over 100 kb, of whom two were containing or contiguous to the MAGMAS locus (Chr16: 4015729–7102036 (∼3 Mb) and Chr16: 8728764–8840025 (∼1.1 Mb)). Also, by using PLINK [11], we were able to reconstruct haplotypes from exome data of both patients, and here again, we identified a 3.6 Mb shared homozygous region located between SNP markers rs3730119 and rs4616299 on chromosome 16p13.3 (Chr16:4,057,603–7,657,432). Genotyping of microsatellite markers allowed to refine this region and define a minimal ancestral common homozygous haplotype spanning 1.9 Mb between markers D16S758 and D16S243 containing MAGMAS (Figure S1). Altogether, our results demonstrate that the c.226A>G transition in MAGMAS (p.Asn76Asp) is likely to be pathogenic in both families. MAGMAS is located on chromosome 16p13.3 and is composed of five exons. It encodes MAGMAS, an essential component of the pre-sequence translocase-associated protein import motor (PAM) that regulates pre-protein translocation into the mitochondrial matrix [12].

Gene expression analysis

A naturally occurring read-through transcription between the neighboring CORO7 (coronin 7) and MAGMAS genes on chromosome 16 is stated in the Genome UCSC browser (NM_001201479). To determine if the fusion protein is involved in the pathogenesis of the severe skeletal manifestations observed in the affected cases, we checked the expression of CORO7, MAGMAS and CORO7-MAGMAS fusion transcript in human cultured control cells (fetal chondrocytes, osteoblasts and skin fibroblasts) by Reverse Transcriptase PCR. CORO7, MAGMAS and CORO7-MAGMAS were all expressed in skin fibroblasts. MAGMAS expression in bone and cartilage was confirmed. On the other hand, neither CORO7 nor CORO7-MAGMAS was expressed in significant amounts in chondrocytes and osteoblasts (Figure 4A). This result was confirmed by quantitative real-time RT-PCR (Figure 4B): The ratio of MAGMAS in chondrocytes versus in fibroblasts is 0.82 whereas the ratios of CORO7 and CORO7-MAGMAS in chondrocytes versus in fibroblasts are 0.24 and 0.12, respectively.

Immunolocalization of MAGMAS protein in wild type growth plates

To further confirm the specific involvement of MAGMAS in skeletogenesis, we performed immunohistochemical staining of fixed femurs isolated from wild type (WT) mice at different developmental stages from E16.5 to 2 weeks of age. Type X collagen staining was also performed in order to mark the terminally differentiated hypertrophic chondrocytes. MAGMAS protein was expressed in bone and cartilage. Expression of the protein was mostly observed in the hypertrophic zone of the cartilage mainly at postnatal ages and in the primary spongiosa of trabecular bone. Its expression in proliferative chondrocytes and in the resting zone of the cartilage was not detected (Figure 5).

Functional complementation

Human and murine forms of MAGMAS present a very high degree of homology. Of the 125 amino acids comprising the mature protein, 120 are identical [13]. The evolutionary conservation of MAGMAS from yeast to mammals [14] and the lethal phenotype of its deletion (pam16Δ) in Saccharomyces cerevisae [15] and Drosophila [16] imply a fundamental role of this protein in cell growth and development. Previous studies have shown that human MAGMAS exhibits a complete growth complementation of deleted pam16 yeast cells at all temperatures [17].

In order to obtain in vivo evidence of the deleterious nature of the identified mutation, we performed a yeast complementation assay using genetic approaches as described previously by Sinha et al. [17]. We first found that at 28°C (301.15°K), both WT and mutant MAGMAS genes were able to complement the growth defect of the pam16Δ cells since spores expressing MAGMAS_WT and MAGMAS_Asn76Asp were viable when spotted onto G418 selective medium (Figure 6A, 6B). To figure out if the Asn76Asp mutant confers a deleterious growth phenotype at non-permissive temperatures, tenfold serial dilutions of pam16Δ spores expressing WT or mutant MAGMAS were subsequently spotted onto fermentable (YPD) and non-fermentable (YPG) growth media and incubated at different temperatures for 3 days. As shown in Figure 6C, the mutant strain pam16Δ-MAG_Asn76Asp that presented a normal growth phenotype on both YPD and YPG media at 28°C (301.15°K), grew more slowly than pam16Δ-MAG_WT on YPD at 34°C (307.15°K) and 36°C (309.15°K). A more severe growth defect of pam16Δ-MAG_Asn76Asp on YPG was evidenced by a very slow-growth at 34°C (307.15°K), and inviability at 36°C (309.15°K). To confirm this temperature-sensitive growth phenotype, we secondly tested the ability of MAGMAS_Asn76Asp to rescue the inviability of the pam16Δ strain by plasmid shuffling using the haploid YPH499pam16Δ strain transformed with a plasmid encoding the yeast WT PAM16 (yPAM16). Unlike WT MAGMAS expressing strains, the mutant pam16Δ-MAG_Asn76Asp exhibited an obvious slow-growth at 34°C and was unable to grow at 36°C (309.15°K) on YPD (Figure 6D). On the other hand, pam16Δ-PAM16-MAG_Asn76Asp cells carrying a WT copy of yPAM16 in addition to the mutant MAGMAS expressing plasmid, showed a normal growth phenotype at all temperatures.

Protein import and MAGMAS expression in yeast extracts

MAGMAS is an essential component of PAM, a peripheral mitochondrial inner membrane complex, that is crucial for preprotein translocation into the mitochondrial matrix [12], [18]. Preproteins targeted to the matrix cross the outer mitochondrial membrane through the TOM complex then pass through the inner membrane via the TIM23^MOTOR complex which is composed of the TIM23 translocase associated with the PAM complex [19] (Figure S2). To assess the transport function of the mutant MAGMAS_Asn76Asp, we evaluated the protein translocation activity by following the mitochondrial import of the Hsp60 protein. Non-processed mitochondrial Hsp60 was accumulated in significant amounts in the pam16Δ-MAG_Asn76Asp strain when the mutant phenotype was induced (37°C, 310.15°K), but was undetectable in WT pam16Δ-PAM16 and pam16Δ-MAGMAS expressing cells (Figure 7A). We then evaluated the expression levels of MAGMAS protein in the different strains, by performing a Western blot using Human anti-MAGMAS antibodies. Immunoblotting detected a lower level of mutant MAGMAS when the strain was shifted at the non permissive temperature, compared to the amount of MAGMAS_WT in the same condition (Figure 7B).

Mitochondrial morphology and peroxisomal biogenesis defects in the MAGMAS_Asn76Asp expressing strain

Previous studies in mammalian cells, yeast and drosophila have shown the crucial role of MAGMAS and its respective homologs in mitochondria biogenesis [12], 17,20. Altered mitochondria morphology observed in pam16 yeast mutants and in Blp-knowdown drosophila cells led us to study the aspect of these organelles in our mutant strains [20], [21]. Microscopic examination of mitochondria revealed an unhealthy feature in the mutant strain consisting in punctuate and highly fragmented mitochondria that differed from the normal reticulated mitochondrial network observed in WT MAGMAS expressing cells. This unusual aspect was more obviously seen in the mutant strains shifted at the non permissive temperature (Figure 8A).

Furthermore, the recent finding of impaired peroxisomes biogenesis in pam16 yeast mutants [21] prompted us to check the morphology of these organelles in our mutant strains. A remarkable aspect (vacuole labeling), indicative of pexophagic degradation, was only observed in pam16Δ-MAG_Asn76Asp cells grown at the non permissive temperature compared to pam16Δ-MAG_WT cells (Figure 8B).

Discussion

Here, we report the identification of a homozygous missense MAGMAS mutation, (p.Asn76Asp), in patients from two unrelated Lebanese families, affected with a rare lethal spondylometaphyseal dysplasia recently described by Mégarbané et al. [6], [7]. By reconstructing haplotypes from exome and STR genotyping data from patients F1-IV.3 and F2-IV.3, we were able to identify a minimal common ancestral homozygous haplotype around MAGMAS, spanning 1.9 Mb on chromosome 16p13.3, suggesting the existence of a founder mutation in gene. MAGMAS, also referred to as PAM16, was first identified in 2001 by Jubinsky et al. [13], as a gene encoding a novel mitochondria-associated protein involved in granulocyte macrophage colony stimulating factor (GM-CSF) signal transduction and expressed by nearly all mammalian cells. The observation of a MAGMAS mutation in a severe chondrodysplasia supports a specific role of MAGMAS protein in the endochondral ossification process. Its expression in the notochord and the somites has been reported at day 12.5 embryo and in the cartilage of the developing ribs at day 14.5 embryo [22]. However, its involvement in skeletogenesis has never been described to date. Our findings demonstrate that MAGMAS is expressed in trabecular bone and cartilage and more specifically by differentiated chondrocytes localized in the hypertophic zone and by osteoblasts at early developmental stages underlining its key role in skeletogenesis. On the other hand, the implication of the fusion protein CORO7-MAGMAS in the pathogenesis of the skeletal dysplasia has been ruled out, as the related transcript was not well expressed in human fetal chondrocytes and human osteoblasts.

Besides the skeletal features, respiratory insufficiency (F1-IV.1 and F1-IV.3) or heart failure (F2-IV.2 and F2-IV.3) were observed, leading to early deaths of the patients [6], [7]. Of note, an elevated MAGMAS level was detected in murine heart and lung in previous studies [22]. These features may also support a key role of MAGMAS in heart and lung.

MAGMAS consists of two main functional domains: a N-terminal hydrophobic region including a membrane association domain (M) and a predicted mitochondria targeting region (T), and a C-terminal J-like domain (J) characterized by three α-helical segments [12], [23]. Importantly, the p.Asn76Asp mutation lies in the helix II of MAGMAS J-like domain. The latter domain forms, with the J-domain of DNAJC19 (ortholog of yeast Pam18 which is another essential component of PAM), a stable heterodimeric subcomplex [19], [24], [25] (Figure S2). Mutations in DNAJC19 preventing MAGMAS/DNAJC19 dimerization have been shown to lead to dilated cardiomyopathy with ataxia syndrome (DCMA) [26], [27]. Considering the localization of the p.Asn76Asp mutation in the J-like domain of MAGMAS, a similar phenotype was likely to be observed in our reported patients. Cardiomyopathy was observed exclusively in patients of the second family and ataxia was not present but all 4 cases died early before 2 years of age. On the other hand, the absence of DCMA in affected cases may also suggest a redundancy in the function of the subcomplex MAGMAS/DNAJC19 in some tissues.

We further demonstrated that the p.Asn76Asp mutation confers to yeast strains a temperature-sensitive phenotype characterized by a very slow-growth at 34°C (307.15°K) and inviability at 36°C (309.15°K), on both fermentable and non-fermentable growth media. Remarkably, the deleterious effect of the mutation was only observed at a homozygous state which is in correlation with the recessive mode of inheritance of the spondylodysplastic dysplasia reported in this study.

Our findings of impaired import of mitochondrial-matrix proteins in the mutant strains shifted at non-permissive temperature confirmed the deleterious nature of the identified mutation. Likewise, it has been demonstrated that aa alterations in MAGMAS J-like domain resulted in temperature sensitivity and in vivo protein import defects in yeast cells [17], [28]. The low expression level of MAGMAS detected in the mutant strain at the non permissive temperature indicates that the p.Asn76Asp mutation confers protein instability. Additionally, we observed unusual punctuated and highly fragmented mitochondria together with enhanced pexophagy in the mutant MAGMAS strain at the non permissive temperature, supporting a cell death induction in the yeast mutant strain. These data are consistent with the finding of Roy et al., in Drosophila [20] and confirm the crucial role of MAGMAS for cell integrity and viability.

The association between a mitochondrial protein and a skeletal dysplasia is striking as skeletal manifestations are rarely seen in mitochondrial disorders. To date, the only severe skeletal dysplasia identified as having a possible mitochondrial basis is the anauxetic dysplasia [MIM 607095], which is allelic to two milder forms of SD, namely metaphyseal dysplasia Mc Kusick type [MIM 250250] and metaphyseal dysplasia without hypotrichosis [MIM 250460]). The gene responsible for these SD is RMRP [MIM 157660], which encodes the RNA component of the mitochondrial RNA processing endoribonuclease (MRP) complex [29]–[31]. Of note, RNase MRP is not exclusively found in the mitochondria but also in the nucleus. Previous studies have shown that mitochondrial and nucleolar RNase MRP have identical RNA components but distinct enzymatic activities and protein components [32]. The recent identification of mutations in POP1 (Processing of precursor 1) [MIM 602486], which encodes a protein component of the nuclear ribonucleoprotein complex in two siblings with anauxetic dysplasia [33] may support the involvement of the nuclear ribonucleoprotein complex rather than the mitochondrial complex in the etiology of these severe skeletal dysplasias. Similarly, besides its role in preprotein import into mitochondria, Magmas may also be involved in additional important mechanisms as suggested by Short et al. [21]. The detection of different human transcripts of MAGMAS, with some lacking the mitochondrial transit peptide (mTP) leader sequence [14] may support a novel extra mitochondrial activity of this protein.

In conclusion, this is the first report linking a mutation in MAGMAS with a human developmental disorder and supporting a key role for this mitochondrial protein in the ossification process.

Materials and Methods

Patients

Informed consents for participation, sample collection and photographs publication were obtained from both families in compliance with the ethics guidelines of “Conseil de recherche de l'Université Saint-Joseph” (number PTS5). Peripheral blood samples were collected from all available family members.

DNA extraction and exome sequencing

Genomic DNA was extracted from peripheral blood samples by standard salt-precipitation methods [34]. Exome sequencing was carried-out on 3 individuals from family 2 (F2-III.1, F2-III.2 and F2-IV.3) and the proband from family 1 (individual F1-IV.3) (Figure 2). Targeted exome sequencing, library preparation, capture and sequencing were performed by the French National Genotyping Center (CNG). Exomes were captured and enriched using the in solution Agilent SureSelect Human All Exon kit v3.0 [35] and then sequenced on an Illumina HiSeq2000, using a paired-end 100-bp read sequencing protocol. Image analysis and base calling were performed using the Illumina Data Analysis Pipeline Software 1.5 with default parameters. Raw data were mapped to the current built of the human genome (hg19) by using BWA 0.75 [36]. Variant calling was subsequently performed using GATK 2.5.2 [37] and annotation was done with ANNOVAR [38]. In order to predict the deleterious effect of the identified sequence variations, different bioinformatics tools were applied; such as MutationTaster (http://www.mutationtaster.org/), SIFT (http://sift.bii.a-star.edu.sg/) [9], PROVEAN (http://provean.jcvi.org/index.php/) [10] and PolyPhen-2 (http://genetics.bwh.harvard.edu/pph2/).

Linkage analysis

Searching for shared homozygous regions using exome data and reconstructing haplotypes were realized using the free software package PLINK (http://pngu.mgh.harvard.edu/purcell/plink/) [11].

Genotyping of STR markers

The following STR markers flanking the PAM16 gene and spanning a 8.7 Mb region on chromosome 16p13.3, were amplified using fluorescently labelled primers: D16S521 (AFMA139WG1), D16S3024 (AFMA134XC1), D16S3070 (AFMA33YH1), D16S3027 (AFMA154WC9), D16S758 (CHLC.ATA24C09), D16S3084 (AFMB080YH5), D16S3072 (AFMB015WA9), D16S3134 (AFMA059WF9), D16S510 (AFM312VD5), D16S423 (AFM249YC5), D16S423 (AFM249YC5), AFM339XG1 (RH4035) and D16S3135 (AFMA059XC1). Marker heterozygosity rate ranged from 0.71 (D16S521 and D16S510) to 0.87 (D16S3082 and D16S3027). PCR amplifications were performed in a 30-µl reaction mix, containing 50 ng of genomic DNA, 0.5 units of Taq polymerase (Invitrogen), 1.5 mM MgCl₂, 5 pmol of each primer, and 0.2 mM of each dNTP. After an initial denaturation step of 5 min at 96°C, 35 cycles of amplification (94°C for 30 sec; 55°C for 30 sec; 72°C for 30 sec) were performed, followed by a final elongation step of 5 min at 72°C. After amplification, PCR products were subjected to denaturing electrophoresis on an ABI PRISM 3100 DNA Sequencer (Applied Biosystems).

Capillary sequencing

Genomic and cDNA sequences of MAGMAS (also known as PAM16) were obtained from UCSC Genomic Browser on Human February 2009 (NM_016069). Primers used for PCR amplification were designed using Primer3 software (http://frodo.wi.mit.edu) to amplify the region surrounding the mutation detected by exome sequencing in exon 4. PCR products were purified by exonuclease I/Shrimp Alakaline Phosphatase treatment (ExoSAP-IT; Fisher Scientific SAS, Illkirch, France) according to the manufacturer's instructions and both strands were sequenced using the Big Dye Terminator v1.1 Cycle Sequencing Kit (Applied Biosystems). Sequence reactions were purified on Sephadex G50 (Amersham Pharmacia Biotech, Foster City, CA) and capillary electrophoresis was performed on Genetic Analyser 3100 (Applied Biosystems). Electropherograms were analyzed on the Sequence Analysis Software version 5.2 (Applied Biosystems) and aligned with the WT MAGMAS gene sequence using ChromasPro version 1.22 (Technelysium, Queensland, Australia).

Reverse Transcriptase PCR

Total RNAs were extracted from human primary cultured control cells (fetal chondrocytes, osteoblasts and skin fibroblasts) via the RNeasy Mini Kit (Qiagen). Complementary DNA (cDNA) was synthesized using the SuperScript II Reverse Transcriptase (Invitrogen Life Technologies, Carlsbad, CA, USA) according to the manufacturer's protocol. Sequences of MAGMAS (NM_016069), CORO7 (NM_001201472), CORO7-MAGMAS (NM_001201479) and GAPDH (NM_002046) were obtained from UCSC Genomic Browser on Human February 2009 and specific primers for cDNA amplification were designed through the Primer3 software as follows:

MAGMAS-Forward: 5′-TACCTGGCCCAGATCATTGT-3′, MAGMAS-Reverse: 5′-GCATCTGCCCTTTTTCTCTG-3′; CORO7-Forward: 5′-GTCCTGGTGTACTGGGCATT-3′, CORO7-Reverse: 5′-AGAGTCAGGGTCCAGCAGAG-3′; CORO7-MAGMAS-Forward: 5′-CTCAGCGCAGTACCTGGAAG-3′, CORO7-MAGMAS-Reverse: 5′-GCCACCCACGGATTTATCAT-3′ and GAPDH-F: 5′-ATGTTCGTCATGGGTGTGAA-3′, GAPDH-R: 5′- TTCAGCTCAGGGATGACCTT-3′.

Quantitative Real–Time RT-PCR

Reactions were conducted in a 96-well plate with the ABI 7500 Sequence Detection System (Applied Biosystems). Primers against MAGMAS, CORO7 and CORO7-MAGMAS were designed using Primer Express software (Applied Biosystems) (Figure S3). PCR was performed in a 20-µl reaction volume containing 10 µl Power SYBR-GreenPCR Master Mix of SYBR Green PCR Master Mix buffer (Applied Biosystems, Foster City, CA, USA), 10 pmol forward and reverse primers and 5 ng of RNAse H-treated (Invitrogen Life Technologies) cDNAs obtained from reverse transcription of DNAse I-treated total RNA extracted from control fibroblasts and chondrocytes cDNA. The reaction cycling conditions were 95°C for 10 min, followed by 40 cycles of 95°C for 15 s and at 60°C for 1 min. Sequence Detection Software was used for exporting the threshold cycle data and further analyzing differences in threshold cycle values (ΔCt) between the test locus and the control locus. The expression values are normalized to that of GAPDH. For each gene, the expression ratio for a sample was calculated as the ratio between the average mRNA signal in the sample and the signal from fibroblasts

Immunohistochemical staining

Femurs isolated from WT mice (at different developmental stages including E 16.5-day, birth and 2 weeks of age) were fixed with 4% paraformaldehyde and embedded in paraffin. Femur sections were stained with Hematoxylin and eosin (H&E) using standard protocol for histological analysis or were subjected to immunohistochemical staining. For immunohistochemistry, sections were stained with Ab specific to MAGMAS (Abcam plc., UK) at 1/200 dilution or with Ab specific to Col X (Quartett) at dilution 1/20 using DAKO EnVision kit. Images were captured with an Olympus PD70- IX2-UCB microscope. All experimental procedures were approved by the French Animal Care and Use Committee.

Plasmids, yeast strains and growth conditions

Plasmids and yeast strains used in this study are described in table 1. The yeast expression vector pRS415-WT-MAGMAS, encoding full-length human WT (WT) MAGMAS under the TEF promoter was provided by D'Silva et al. [17]. The pRS415-Asn76Asp-MAGMAS plasmid carrying the point mutation p.Asn76Asp was generated from pRS415-WT-MAGMAS through a site directed mutagenesis using the QuickChange protocol (Stratagene) according to the manufacturer's recommendations. Direct sequencing was performed to confirm the successful introduction of the mutation. For construction of yeast strains expressing the human WT/mutated MAGMAS, we used Saccharomyces cerevisiae Δpam16 strains (provided by Delahodde A). For the yeast complementation assay, the diploid strain BY4743 carrying the deleted pam16 allele (pam16Δ::kanMX6) in a heterozygous state was first transformed with pRS415-WT-MAGMAS or pRS415-Asn76Asp-MAGMAS plasmids, accordingly and then subjected to sporulation and tetrad dissection using standard techniques [39]. Spores from a single ascus were first grown at 28°C (301.15°K) in YPD media then replica-plated to G418 plates for the selection of spores carrying the deleted pam16 allele and to Leucine drop-out plate for the selection of pRS415MAGMAS positive spores. The haploids BY-pam16Δ-MAG_WT and BY-pam16Δ-MAG_Asn76Asp were thus generated and spotted onto YPD and YPG to analyze their phenotype. For the construction of YPH499-pam16Δ-pPAM16-MAG_WT and YPH499-pam16Δ-pPAM16-MAG_Asn76Asp, the haploid YPH499pam16Δ::ADE2 strain containing a plasmid encoding the yeast PAM16 gene (kindly provided by Guiard B), was transformed either with pRS415-WT-MAGMAS or pRS415-Asn76Asp-MAGMAS plasmids, respectively then patched onto the selection plates (lacking leucine and uracile). Transformants were subsequently replicated to 5-Fluoroorotic (5-FOA) containing plates, recovered and grown at 28°C (301.15°K), thus generating YPH499-pam16Δ-MAG_WT and YPH499-pam16Δ-MAG_Asn76Asp strains that have lost the plasmid containing the yeast PAM16 gene and the URA3 marker gene [40]. Ten fold serial dilutions of the different generated strains were subsequently spotted onto YPD plates then incubated at different temperatures for 3 days. Unless otherwise noted, yeast cultures were grown at 28°C (301.15°K). YPD (1% bactopeptone, 1% yeast extract, and 2% glucose) and YPG (1% bactopeptone, 1% yeast extract, and 2% glycerol) were used as rich fermentable and non-fermentable growth media respectively. W0 (0.67% yeast nitrogen base without aa and 2% glucose) supplemented with appropriate nutritional requirements according to the strains was used as a minimal medium. For solid media, 2% bacto agar (Difco) was added.

Protein import and western blotting

For the import activity analysis, YPH499pam16Δ-PAM16, YPH-pam16Δ-MAG_WT and YPH-pam16Δ-MAG_Asn76Asp cells were first grown to OD_{600 nm} of 1.0 at 28°C (301.15°K) in rich medium then shifted to 37°C (310.15°K), for 8 h (28800 s) to induce non-permissive conditions. Cell pellets were suspended in ice-cold TCA buffer (TrisHCl pH 8 20 mM, Ammonium Acetate 50 mM, EDTA 2 mM, protease inhibitor solution, and TCA10%) then lysed by adding glass beads and shaking vigorously for 300 s (5 fold 60 s in ice between). After centrifugation of 600 s (13,000 rpm) at 4°C (277.15°K), the proteins were solubilized in TCA-Laemmli loading buffer (SDS/glycerol stock solution 48%, Tris EDTA 40%, β-mercaptoethanol 0.5%, PMSF 0.2%, protease inhibitor solution 0.2%) and separated by SDS-PAGE. The proteins were subsequently transferred to a nitrocellulose membrane that was first blocked, then incubated 1 h at room temperature with a blocking solution containing the monoclonal HSP60 Ab (provided by Dujardin G.) at 1/15000 dilution. After several washes with TBS1X (Tris Buffered Saline 1X), the membrane was immersed for 1 h at room temperature with HRP – conjugated secondary anti-mouse Ab (at 2/15000 dilution in a blocking solution) then washed again with TBS1X. Secondary antibodies bound to the membrane were detected by chemiluminescence using the Thermo Scientific PierceECL 2 according to the manufacturer's instructions. The same protocol was used for the immunodetection of MAGMAS protein in YPH499, YPH-pam16Δ-MAG_WT and YPH-pam16Δ-MAG_Asn76Aspstrains. For this purpose, we incubated the membrane with a blocking solution containing the polyclonal anti-MAGMAS Ab (Abcam plc., UK) at a 1/1000 final concentration. Therefore, the secondary Ab used in this case was HRP – conjugated secondary anti-rabbit Ab (at 2/15000 dilution in a blocking solution). Phosphoglycerate kinase (PGK) and Porin were detected using anti-PGK (Invitrogen) and anti-Porin (Invitrogen), accordingly.

Fluorescence microscopy

In order to label mitochondria, YPH-pam16Δ-MAG_WT and YPH-pam16Δ-MAG_Asn76Asp strains were transformed with the pYX232-mtGFP plasmid that encodes the mitochondrial matrix-targeted GFP [41] and grown for 2 days at 28°C (301.15°K). Cells were secondly incubated overnight at the non-permissive temperature to induce the phenotype. For the visualization of peroxisomes, BY-pam16Δ-MAG_WT and BY-pam16Δ-MAG_Asn76Asp were transformed with pUG34DsRed.SKL encoding the red fluorescent protein DsRed containing the C-terminal peroxisomal targeting signal type 1 (PTS1, DsRed-SKL). Transformants were grown for 2 days at 28°C (301.15°K) in minimal medium. Harvested cells were washed in Phosphate-buffered saline 1X (PBS1X) then transferred to microscopic slides for imaging. For observations at high temperature, cells grown 2 days at 28°C (301.15°K) in minimal medium, were shifted to 37°C (310.15°K), for 6 hours then harvested, washed and transferred to slides for visualization. The slides were examined with a DMIRE2 microscope (Leica, Deerfield, IL). For GFP fluorescence, filters used were 450/490 nm for excitation and 500/550 nm for emission whereas for Ds.Red fluorescence, filters were 542/582 nm for excitation and 604/644 nm for emission. Images were captured using a CCD camera (Roper Scientific, Tucson, AZ) and Metamorph software (Universal Imaging, WestChester, PA) was used to deconvolute Z-series and treat the images.

Web resources

The URLs for data presented herein are as follows:

Online Mendelian Inheritance in Man (OMIM), http://www.omim.org

Polyphen, http://genetics.bwh.harvard.edu/pph2/

MutationTaster http://www.mutationtaster.org/

SIFT http://sift.bii.a-star.edu.sg/

PROVEAN http://provean.jcvi.org/index.php/

PLINK http://pngu.mgh.harvard.edu/purcell/plink/

Primer3 software http://frodo.wi.mit.edu/

Supporting Information

Attachment 1

Attachment 2

Attachment 3

Attachment 4