Mutations in the TGF-β repressor SKI cause Shprintzen-Goldberg syndrome with aortic aneurysm.

Elevated transforming growth factor (TGF)-β signaling has been implicated in the pathogenesis of syndromic presentations of aortic aneurysm, including Marfan syndrome (MFS) and Loeys-Dietz syndrome (LDS). However, the location and character of many of the causal mutations in LDS intuitively imply diminished TGF-β signaling. Taken together, these data have engendered controversy regarding the specific role of TGF-β in disease pathogenesis. Shprintzen-Goldberg syndrome (SGS) has considerable phenotypic overlap with MFS and LDS, including aortic aneurysm. We identified causative variation in ten individuals with SGS in the proto-oncogene SKI, a known repressor of TGF-β activity. Cultured dermal fibroblasts from affected individuals showed enhanced activation of TGF-β signaling cascades and higher expression of TGF-β-responsive genes relative to control cells. Morpholino-induced silencing of SKI paralogs in zebrafish recapitulated abnormalities seen in humans with SGS. These data support the conclusions that increased TGF-β signaling is the mechanism underlying SGS and that high signaling contributes to multiple syndromic presentations of aortic aneurysm.

The TGF-β family of cytokines influences a diverse repertoire of cellular processes, including cell proliferation, differentiation, survival and synthetic activity. One of three TGF-β ligand isoforms (TGF-β1, -β2 or -β3) initiates signaling by binding to the TGF-β receptor complex, composed of type 1 and type 2 receptor subunits (TβRI and TβRII, respectively) [11]. The complex then transmits the signal through the canonical (SMAD dependent) pathway or noncanonical mitogen activated protein kinase (MAPK) cascades, prominently including extracellular signal-regulated kinase (ERK).

The ability of the nuclear SMAD proteins to precisely regulate gene transcription is further modulated by the recruitment of transcriptional co-activators such as p300 and CREB binding protein (CBP), and co-repressors such as the Sloan Kettering Institute proto-oncoprotein, SKI[12]. The SKI family of proteins, which also includes the SKI-like protein SKIL, negatively regulate SMAD-dependent TGF-β signaling by impeding SMAD2/3 activation, preventing nuclear translocation of activated R-SMAD/Co-SMAD complexes and inhibiting TGF-β target gene output by competing with p300/CBP for SMAD binding and recruiting transcriptional repressor proteins, such as mSin3A and HDAC1 (Supplementary Figure 1)[13-15]. A role for SKI in the regulation of MAPK cacades has not been described[16,17].

Dysregulated TGF-β signaling has been implicated in the pathogenesis of syndromic presentations of aneurysm. Excessive canonical and noncanonical TGF-β signaling is seen in the aortic wall and other diseased tissues in mouse models of MFS, a systemic connective tissue disorder caused by mutations in the gene (FBN1) encoding the extracellular matrix protein fibrillin-1[1,3,18-20]. Furthermore, pharmacologic antagonism of TGF-β and/or ERK signaling in mouse models of MFS has been shown to ameliorate multiple disease manifestations, including aortic aneurysm, suggesting that high TGF-β signaling drives disease progression[3,4].

An excess of canonical TGF-β signaling has also been demonstrated in the aortic wall of patients with LDS[2]. However, this has been difficult to reconcile with the character of the underlying mutations. LDS is predominantly caused by heterozygous missense substitutions in the kinase domains of either TβRI or TβRII (encoded by the TGFBR1 and TGFBR2 genes, respectively)[2,5]. When mutant proteins are expressed in cells naïve for the corresponding receptor subunit, there is failure to propagate SMAD-dependent TGF-β signaling[21]. More rarely, LDS-like phenotypes can be caused by haploinsufficiency for SMAD3 or the TGFβ2 ligand, positive effectors of TGF-β signaling[22,23]. Notably, heterozygous loss-of-function mutations in SMAD3 or TGFB2 also associate with high TGF-β signaling in the aortic wall of patients and mouse models[22,23]. Taken together, these apparently contradictory data have engendered considerable controversy regarding the precise role of TGF-β in the pathogenesis of aortic aneurysm.

Shprintzen-Goldberg syndrome is a systemic connective tissue disorder of unknown etiology that includes virtually all the craniofacial, skeletal, skin and cardiovascular manifestations of MFS and LDS, with the additional findings of mental retardation and severe skeletal muscle hypotonia[6,8]. We hypothesized that aberrant TGF-β activity also underlies SGS and that identification of the genetic basis would inform our understanding of other syndromic presentations of aneurysm[2]. We performed whole exome sequencing for a single affected SGS child-unaffected parent trio. An average of 6.9 gigabases of sequence was generated per individual as paired-end 75 base pair reads, out of which >98.9% mapped to the human reference genome (UCSC hg19). This revealed only one variant, a heterozygous missense change in exon 1 of the SKI gene (NM_003036; c.347G>A, p.Gly116Glu) that was not present in SNP databases, was predicted to be damaging (PolyPhen-2 score: 0.999, SIFT score: 0.05), was not present in either parent, and was a strong functional candidate based upon a described relationship to TGF-β signaling (Supplementary Figure 2)[24,25]. We subsequently sequenced SKI in 11 other sporadic SGS patients (Table 1 and Figure 1A) and identified heterozygous variants in 9 including 8 missense mutations and one 9 base pair deletion (Supplementary Figure 3). These mutations were absent from dbSNP 134, the 1000 genomes project and over 10,000 exomes reported on the NHLBI Exome Variant Server, and confirmed to be de novo when testing of the unaffected parents was possible (5/9). Collectively, 10 mutations in 10 patients with SGS were identified in SKI by a combination of whole exome and Sanger sequencing, including a recurrent mutation in 2 unrelated probands. Each mutation substituted or deleted residues that show complete evolutionary conservation among members of the SKI gene family (Supplementary Figure 4)[26]. No mutations were identified upon sequencing of SKIL (NM_005414) in the 2 remaining patients. The mutations in SKI clustered in two distinct regions at the N-terminus of the protein (Figure 1B). The first region is located in the SMAD2/3 binding domain of SKI (residues 17-45), while the second region localizes to a portion of the Dachshund-homology domain (DHD) of the SKI protein that mediates binding to SNW1 and N-CoR, proteins that are essential to the transforming activity of SKI and to the recruitment of transcriptional co-repressors such as histone deacetylases, respectively[15,27,28]. One of the mutations (p.Leu21Arg) substitutes an amino acid that was previously shown to be essential for SKI-SMAD3 interaction[29]. Two of the mutations (p.Gly116Glu and p.Gly117Arg) substitute immediately adjacent glycine residues at positions that contribute to an exposed β hairpin loop in the DHD[30]. In silico analysis suggests that both glycine substitutions maintain the β-turn; while there is variation in where the algorithm defines the margins of the adjacent β sheets, the overall structure of the DHD is fully overlapping for the reference and mutant sequences (Supplementary Figure 5). It is therefore possible that the principle consequence of these mutations is disruption of molecular interactions at the protein surface[31]. In this light, it is interesting to note that residues 125-131, which lie immediately distal to the β-turn, have been shown to influence binding of SMAD2 and SMAD3 to their more proximal N-terminal binding site[30,32].

Table 1

Clinical Manifestations in Shprintzen-Goldberg Syndrome (SGS) and Related Disorders

	SGS Ref.7	SGS Ref.8	Pt.1 43y F	Pt.2 6y M	Pt.3 16y M	Pt.4 12y M	Pt.5 22y F	Pt.6 21y M	Pt.7 2y F	Pt.8 6y M	Pt.9 5y M	Pt.10 4y F	SGS (SKI neg.)	Classic MFS	Classic LDS
Features typical of SGS
Craniofacial
Craniosynostosis	7/7	4/7	+	+	+	+	+	+	+	+	+	+	1/2	−	+
Dolichocephaly	16/17	11/14	+	+	+	+	+	+	−	+	+	+		+	+
Hypertelorism	15/16	8/13	+	+	+	+	+	+	+	+	+	+	0/2	−	+
Down-slanting eyes	14/16	11/14	+	+	+	+	+	+	+	+	+	+	2/2	+	+
Proptosis	13/17	11/12	+	−	+	+	+	−	+	+	+	+	1/2	−	−
Malar hypoplasia	15/15	7/8	+	+	+	+	+	+	+	+	+	+	2/2	+	+
High/narrow palate	17/17	12/12	+	+	+	+	+	+	+	+	+	+	2/2	+	+
Micrognathia	16/17	13/14	+	+	+	+	+	+	+	+	+	+	2/2	+	+
Low set ears	13/13	12/14	+	+	+	+	+	+	+	+	+	+	2/2	−	−
Skeletal
Arachnodactyly	15/17	14/14	+	+	+	+	+	+	+	+	+	+	2/2	+	+
Camptodactyly	10/17	6/14	+	+	+	+	+	+	−	−	+	+		+	+
Scoliosis	7/14	11/13	+	+	+	+	+	+	+	−	−	+	2/2	+	+
Pectus deformity	15/17	13/13	+	+	+	+	+	+	+	−	+	+	1/2	+	+
Joint hypermobility	9/14	10/10	+	+	+	+	+	+	+	+	+	+	2/2	+	+
Joint contracture	8/17	2/10	+	+	+	+	+	+	−	+	+	+	1/2	−	−
C1/2 spine malform.	4/16		−	−		+	+				+	+		−	+
Neuromuscular
Hypotonia	14/16	11/12	+	+	+	+	+	+	+	+	+	+	1/2	−	−
Developmental delay	15/17	12/14	+	+	+	+	+	+	+	+	+	+	2/2	−	−
Mitral valve prolapse	4/11	6/17	+	+	+	−	−			−	+	+	0/2	+	+
Aortic root dilatation1	5/17	3/11	+	+	+	+	+	+	+	−	−	+	0/2	+	+
Other features of MFS or LDS
Dural ectasia			+	+								+		+	+
Ectopia lentis	1/152		−	−	−	−	−	−	−	−	−	−	0/2	+	−
Cleft palate		2/11	−	−	−	+	−	−	+	+	−	−	0/2	−	+
Broad/Bifid uvula			−	+		−					−	+	1/2	−	+
Club foot deformity			−	−	+	−	−		−	+	−	−	0/2	−	+
Arterial tortuosity			−	−		−						+		−	+
Other aneurysm			+3	−		−				+4		−		−	+

Aortic root z-score ≥ 2.0 (range 2.1-5.7).

Atypical patient with FBN1 mutation.

Splenic artery aneurysm.

Splenic artery aneurysm with spontaneous rupture. Blank cells, not determined.

Figure 1

SKI mutations in patients with Shprintzen-Goldberg syndrome (SGS). A) Clinical features and mutations seen in SGS patients. Illustrated features involve the craniofacial (abnormal head shape due to craniosynostosis, widely-spaced eyes, small and receding chin and high-arched palate), skeletal (long fingers, joint contractures, chest wall deformity, spine curvature, foot deformity) and cardiovascular (aortic root aneurysm indicated by black arrowheads and mitral valve prolapse indicated by white arrowhead) systems. Pedigrees indicate that all cases were sporadic (affected, black symbols; unaffected, open symbols). Mutation status is indicated below each individual (−/−, mutation negative; +/−, heterozygous; blank, not available for testing). Permission to publish photographs was obtained from the affected individuals or their parents. B) Location of mutations in relation to known binding sites (highlighted in yellow) of SKI binding partners. The position of the Dachshund homology domain is indicated.

In order to assess the functional consequence of SKI mutations identified in SGS, we monitored TGF-β signaling in primary dermal fibroblasts derived from 2 SGS patients and 2 controls. Western blot analysis showed excessive SMAD2/3 and ERK1/2 phosphorylation (pSMAD2/3 and pERK1/2) in patient cells compared to controls, both at baseline and after acute (30 minute) stimulation with exogenous TGF-β2 (Figure 2A). This implies a loss of suppression of TGF-β dependent signaling cascades in SGS cells. In contrast, there was no difference in the activation of JNK or p38 between patient and control cells at baseline or in response to TGF-β2. These signaling alterations directly parallel previous findings in the aortas of both patients and mouse models of MFS and LDS[4,23]. We also observed a significant increase in mRNA expression for TGF-β dependent genes in SGS fibroblasts, including COL1A1, COL3A1, FN1, VIM and CDKN1A (encoding types I and III collagen, fibronectin, vimentin and p21, respectively) (Figure 2B). This is in keeping with prior work showing that transcription of these genes is normally suppressed by SKI[33]. TGF-β also induces transcription of genes encoding negative regulators of itself, including SKI, SKIL, and SMAD7, as part of an autoregulatory loop. We observed increased mRNA expression for each of these genes in SGS cells compared to controls, suggesting that SKI normally suppresses the expression of itself as well as other negative regulators of TGF-β signaling. Interestingly, the TGF-β-regulated target genes CTGF and SERPINE1 showed equal expression in SGS and control cells, suggesting either gene-specific differences in their sensitivity to attenuation of SKI function or to the regulatory influence of other transcriptional repressors such as SKIL or SMAD7.

Figure 2

Assessment of TGF-β signaling in dermal fibroblasts. A) Western blot analysis of pSMAD2, pSMAD3, pERK1/2, pJNK1/2, pp38 and β-Actin (loading control) at steady state (baseline) and in response to TGF-β2. Representative Shprintzen-Goldberg syndrome (SGS) and control blots are shown (for a single SGS patient and single control with 3 biological replicates for each). Graphs display Western blot quantification of 3 biological replicates in each of 2 SGS patients and 2 controls. All data are normalized to β-Actin. The control baseline value for each graph was set to 1.0; all other values represent fold change in comparison to this sample. B) Relative mRNA expression for TGF-β target genes (normalized to β-Actin), as assessed by quantitative polymerase chain reaction analysis. Quantification was performed on 3 biological replicates in each of 2 SGS patients and 2 controls. The upper and lower margins of the box define the 75th and 25th percentiles respectively; the internal line defines the median and the whiskers define the range. NS, not significant; * p<0.05; ** p<0.01; † p<0.001; †† p<0.0001.

Despite the near-complete phenotypic overlap between LDS and SGS, it is notable that the aneurysm phenotype in SGS is less penetrant, less diffuse (generally restricted to the aortic root) and less severe than that seen in LDS. We reasoned that this might relate to the temporal and regional expression pattern of SKI. To test this hypothesis, we performed a developmental survey of SKI expression in wild-type mice (Figure 3A). At embryonic day 13.5, SKI is robustly expressed throughout the vessel wall in the proximal ascending aorta, with less expression in the descending segment, and localizes to both the cytoplasm and nucleus. At birth, aortic expression is somewhat reduced, with predominant expression in the cytoplasm. In adult mice, SKI expression is further reduced in the medial layer of the aortic root, with exclusive cytoplasmic localization; in the more distal ascending aorta SKI expression is largely excluded from the central zone of the aortic media despite some residual expression in the intimal and adventitial layers. Taken together, these data support the hypothesis that the requirement for SKI for proper regulation of TGF-β signaling within the arterial media may predominate in the very proximal aorta at early stages of development.

Figure 3

Assessment of SKI in model systems. A) Murine SKI expression at embryonic day 13.5 (E13.5) is high in the mitral valve (MV), tricupsid valve (TV), pulmonary trunk (Pu), pulmonary valve (PuV), bronchus (Br) and proximal aorta (Ao; both nucleus and cytoplasm). Lower expression is seen in the descending aorta (dAo), esophagus (E), right and left ventricles (RV and LV) and left atrium (LA). Red, SKI; green, myosin heavy chain; blue, nuclei. At birth (P0), SKI is expressed in the aortic root (AoR), distal ascending aorta (ascAo) and endothelial surface of the mitral and aortic (AoV) valves, with less nuclear expression than at E13.5. In adulthood, SKI is expressed exclusively in the cytoplasm of the AoV and the media of the AoR; in the ascAo SKI is expressed in the intima (I) and adventitia (A), but excluded from the central media (M). B) Disruption of skia or skib expression in zebrafish. Top panels, assessment of cardiovascular anatomy at 3 days post-fertilization (dpf). Uninjected embryos display proper cardiac looping with normal relation between the atrium (A, dotted outline) and ventricle (V), and a distinct bulbous arteriosus (arrow) defining the outflow tract (OFT). skia morphants show incomplete looping with an irregular OFT, while skib morphants show failure of looping and an ill-defined OFT. Asterisk, atrioventricular cushion. Lower panels, alcian blue staining of cartilage at 6 dpf. skia and skib morphants display craniofacial cartilage deficits and prominent mandibular malformation. Meckel's cartilage, m; palatoquadrates, pq; ceratohyales, ch; ceratobranchial arches, cb. Scale bars, μm.

We show that the multisystem manifestations of SGS are caused by primary mutations in a prototypical repressor of TGF-β signaling, and that this associates with a cell-autonomous increase in SMAD2/3 activation at steady-state and in the acute-phase response to TGF-β ligand as well as with increased transcriptional output of TGF-β responsive genes. It remains to be determined whether the increased ERK1/2 activation seen in SGS cells manifests loss of a previously unrecognized primary function of SKI or secondary cellular events. The conclusion that release of the suppressive effects of SKI on TGF-β signaling underlies SGS is consistent with prior observation of CNS patterning alterations, skeletal muscle hypoplasia and craniofacial defects (all cardinal manifestations of SGS) in homozygous Ski-targeted mice[34]. It is further strengthened by our ability to replicate SGS-like defects after morpholino based knockdown of the 2 paralogs of mammalian SKI (skia and skib) in zebrafish. Morphant skia and skib zebrafish embryos showed significant craniofacial cartilage deficits including shortened and flat Meckel's cartilage, irregular lengths of palatoquadrates, shortened ceratohyales and depletion of ceratobranchial arches (Figure 3B, Supplementary Table 1). These deficits manifest in larval fish as maxillary hypoplasia, malformed ethmoid plate, micrognathia and microcephaly, and are frequently accompanied by ocular hypertelorism and spinal malformations. Furthermore, ski morphant embryos displayed severe cardiac anomalies, characterized by partial to complete failure in cardiac looping and malformations of the outflow tract (Figure 3B, Supplementary Table 1). Thus, in comparison to Ski-null mice, the zebrafish morphants more closely recapitulated the human SGS craniofacial phenotype and uniquely implicated a requirement for SKI in early cardiovascular development. Prior demonstration of a prominent role for SKI in both neurogenesis and myogenesis reconciles the highly penetrant developmental delay and hypotonia seen in SGS, respectively[35,36]. It is also notable that a recurrent mutation in SMAD4 that strongly impairs TGF-β transcriptional responses by the R-SMAD/Co-SMAD complex has recently been shown to cause Myhre syndrome, a condition characterized by short bones and aortic stenosis, a striking contrast to the bone overgrowth and aortic dilatation seen in SGS[37]. Patients with the recurrent 1p36 deletion syndrome are haploinsufficient for SKI and show some phenotypic overlap with SGS including hypotonia, developmental delay, craniofacial dysmorphism, vertebral abnormalities and structural heart disease[38,39]. Possible reasons for the features that distinguish these two conditions include the involvement of contiguous genes in the 1p36 deletion syndrome and/or the putative dominant-negative potential of mutant forms of SKI in SGS, that might retain their ability to form homodimeric complexes due to structural preservation of the SKI-interacting domain at their C-terminus, with functional deficits imposed by N-terminal mutations that selectively perturb R-SMAD and/or N-CoR interactions. Substitution of neighboring cysteine residues (p.Cys1223Tyr and p.Cys1221Tyr) in an epidermal growth factor-like domain of fibrillin-1 have been reported in association with an SGS-like phenotype including all features of MFS and SGS-specific features such as craniosynostosis and developmental delay[40,41]. While it remains to be determined if this tentative phenotype-genotype correlation manifests a predominant consequence of FBN1 genotype, modifier loci or chance, it seems notable that both fibrillin-1 and SKI are recognized regulators of TGF-β signaling.

In conclusion, this manuscript provides evidence that mutations predicted and observed to enhance TGF-β signaling are sufficient to cause human phenotypes that have variably been associated with low TGF-β signaling states, including craniosynostosis, altered palatogenesis and aortic aneurysm. The development of individualized treatment strategies will require further work to determine whether low and high signaling states achieve the same phenotypic consequence or, as we propose, if complex compensatory events during tissue morphogenesis and homeostasis culminate in a common pathogenic mechanism for apparently disparate disease etiologies[42].

Online Methods

Participants

Patients were recruited from the Connective Tissue Clinic at Johns Hopkins Hospital (H.C.D.), Radboud University Hospital/Antwerp University Hospital (B.L.L. and G.M.), The Nemours Children's Clinic (P.H.A.) and Genetic Medicine Central California (C.J.C.). All skin biopsies and research protocols were collected in compliance with the Institutional Review Board at each respective institution after informed consent was obtained. The diagnosis of SGS was made after exclusion of other known congenital syndromes on the basis of distinguishing phenotypic features. Echocardiograms were performed and interpreted as previously described[43]. Aortic root aneurysm was defined by a maximal aortic root z-score ≥ 2.0.

Whole Exome sequencing

Genomic DNA was extracted from peripheral blood lymphocytes using standard protocols. DNA fragmentation was performed using a Covaris S2 system (Covaris, Woburn, MA, USA), exon capture was performed using the Agilent SureSelect 38 Mb Human All Exon Target Enrichment system (Agilent Technologies, Inc., Santa Clara, CA, USA) and DNA sequencing was performed on an Illumina GAIIx instrument (Illumina, Inc, San Diego, CA, USA), all using standard protocols for a 75 bp paired-end run.

Bioinformatic analysis

Reads were mapped to the human reference genome (UCSC hg19) using the Burrows-Wheeler Aligner and a variant list created using SAMtools and annotated using ANNOVAR[44-46]. Local realignment and recalibration of base call quality scores was performed using the genome analysis toolkit[47,48]. Duplicates were marked using Picard (http://picard.sourceforge.net). We selected for novel variants (absent from dbSNP 134 or 1000 genome project databases), focusing only on exonic non-synonymous, splice site and indel variants. Variants were viewed directly using the integrated genome viewer and excluded if reads were only present in one direction, if ambiguously mapped reads were present or if an indel occurred within 3 bp of the end of the read[49].

Mutation validation and Sanger sequencing of candidate genes

The polymerase chain reactions were performed using a DNA Engine Dyad thermal cycler (Bio-Rad, Hercules, CA, USA). Phusion Flash High Fidelity PCR Master Mix was used in accordance with the manufacturer's instructions for each primer set (Dharmacon, Inc, Lafayette, CO, USA). Primer sequences and reaction conditions are listed in Supplementary Table 2. Cycle sequencing was performed using BigDye Terminator v3.1 kit and an ABI 3730 ×l DNA Analyzer/sequencing machine in accordance with the manufacturer's instructions (Life Technologies, Carlsbad, CA, USA). The samples were purified using QIAquick PCR purification kit (Qiagen, Inc, Valencia, CA, USA).

Cell culture

Primary human dermal fibroblasts were derived from forearm skin biopsies from 2 control individuals and 2 patients with SGS. They were cultured in Dulbecco's modified eagle medium with 10% fetal bovine serum in the presence of antibiotics and passaged at confluence. All cell experiments were conducted in serum-starved media, and stimulation was performed using 5 ng/mL recombinant human TGF-β2 (302-B2, R&D Systems, Minneapolis, MN, USA). Cells were collected at baseline, 30 minutes and 24 hours after stimulation.

Protein and RNA analysis

All protein and RNA was extracted, processed and analyzed using reagents and protocols as described previously[4,23]. The following pre-validated probes were used for mRNA analysis: Hs01060665 (β-ACTIN), Hs00164004 (COL1A1), Hs00943809 (COL3A1), Hs00287359 (FN1), Hs00185584 (VIM), Hs00355782 (CDKN1A), Hs01026927 (CTGF), Hs01126606 (SERPINE1), Hs00161707 (SKI), Hs01045418 (SKIL) and Hs00998193 (SMAD7) (Life Technologies, Carlsbad, CA, USA). In the Western blot analysis separate membranes were used for the control and SGS patients due to the large number of samples. The membranes were processed together at all stages, and exposed to the same piece of imaging film simultaneously, allowing for valid comparisons.

Developmental survey of SKI expression in mouse

Immunohistochemical analysis was performed using protocols as described previously[50]. Samples were incubated with an anti-mouse SKI polyclonal antibody diluted to 1:250 (EMD Millipore Corporation, Billerica, MA, USA) and myocyte-specific mouse monoclonal MF20-c diluted to 1:50 (Developmental Studies Hybridoma Bank, Iowa City, IA, USA). Elastin was detected using autofluorescence at 488 nm. Nuclei were stained with Hoechst dye diluted to 1:10,000 (Life Technologies, Carlsbad, CA, USA).

Zebrafish maintenance, imaging and staining

Adult AB and Tg(kdrl:G-RCFP; fluorescently marked endothelium) zebrafish lines were maintained according to standard methods51,52. Embryos were injected with previously published skia (12 ng) or skib (14 ng) morpholino antisense oligonucleotides (MO) at the 1-2 cell stage and analyzed at 2, 3, 5, 6, and 8 days post fertilization (dpf) for G-RCFP expression via fluorescence microscopy53. Morpholino oligonucleotides are listed in Supplementary Table 2. Prior work had established the specificity of these MOs via concomitant use of random MO and mRNA rescue experiments53. Cartilage staining was performed as previously described on zebrafish embryos fixed in 4% paraformaldehyde at 5, 6 and 8 dpf54. All experiments were performed in accordance with ethical permits by the Johns Hopkins Animal Care and Use Committee.

Statistical Analysis

All quantitative data are shown as box and whisker plots produced using the R statistical package. The upper and lower margins of the box define the 75th and 25th percentiles respectively; the internal line defines the median and the whiskers define the range. Statistical analysis was performed using two tailed t tests in Excel (Microsoft, Redmond, WA, USA). A p value <0.05 was considered statistically significant.