Mutations in SKI in Shprintzen–Goldberg syndrome lead to attenuated TGF-β responses through SKI stabilization
Abstract
Shprintzen–Goldberg syndrome (SGS) is a multisystemic connective tissue disorder, with considerable clinical overlap with Marfan and Loeys–Dietz syndromes. These syndromes have commonly been associated with enhanced TGF-β signaling. In SGS patients, heterozygous point mutations have been mapped to the transcriptional co-repressor SKI, which is a negative regulator of TGF-β signaling that is rapidly degraded upon ligand stimulation. The molecular consequences of these mutations, however, are not understood. Here we use a combination of structural biology, genome editing, and biochemistry to show that SGS mutations in SKI abolish its binding to phosphorylated SMAD2 and SMAD3. This results in stabilization of SKI and consequently attenuation of TGF-β responses, both in knockin cells expressing an SGS mutation and in fibroblasts from SGS patients. Thus, we reveal that SGS is associated with an attenuation of TGF-β-induced transcriptional responses, and not enhancement, which has important implications for other Marfan-related syndromes.
Introduction
Shprintzen–Goldberg syndrome (SGS) is a multisystemic connective tissue disorder. Common features observed in SGS patients include craniofacial, skeletal, and cardiovascular anomalies, ranging from heart valve defects to thoracic aortic aneurysms, all of which are also characteristic of Marfan syndrome (MFS) and Loeys–Dietz syndrome (LDS) (Cook et al., 2015a; Verstraeten et al., 2016; Loeys et al., 2005; Williams et al., 2007). In addition, SGS patients present with craniosynostosis, intellectual disability, and skeletal muscle hypotonia (Shprintzen and Goldberg, 1982; Glesby and Pyeritz, 1989; Greally et al., 1998). All three syndromes have been linked to deregulation of the transforming growth factor β (TGF-β) signaling pathway (Cannaerts et al., 2015).
The TGF-β family of ligands comprises the TGF-βs themselves, Activins, Nodal, bone morphogenetic proteins (BMPs), and growth differentiation factors (GDFs) and they play pleiotropic roles in embryonic development and tissue homeostasis. In addition, their signaling is deregulated in diverse pathologies (Miller and Hill, 2016). They exert their action by binding to type I and type II serine/threonine kinase receptors at the cell surface (TGFBR1 and TGFBR2, respectively, for the TGF-βs) (Massagué, 2012). In the resulting ligand-bound heterotetrameric receptor complex, the type II receptor phosphorylates and activates the type I receptor, which in turn phosphorylates the intracellular mediators, the receptor-regulated SMADs (R-SMADs). Once phosphorylated, the R-SMADs (SMAD2 and SMAD3 in the case of TGF-β, Activin, and Nodal) associate with the common mediator of the pathway, SMAD4. The resulting heterotrimeric complexes accumulate in the nucleus where they interact with other transcriptional regulators to activate or repress target gene expression (Massagué, 2012). Two highly related co-repressors, SKI and SKIL (formerly known as SnoN), act as negative regulators in the pathway (Deheuninck and Luo, 2009; see below).
The role of deregulated TGF-β signaling in Marfan-related syndromes is controversial. MFS is caused by loss-of-function mutations in the extracellular matrix protein, Fibrillin 1 (FBN1) (Dietz et al., 1991). These mutations are thought to increase the bioavailability of TGF-β ligands, as FBN1 binds the latent form of the TGF-βs (Neptune et al., 2003; Kaartinen and Warburton, 2003). Supporting the idea that excessive TGF-β signaling contributes to the manifestations of MFS, a TGF-β neutralizing antibody significantly improved the lung phenotype in a mouse model of MFS (homozygous Fbn1mgΔ) (Neptune et al., 2003; Cannaerts et al., 2015) and reduced the occurrence of aortic aneurysms in the Fbn1C1039G/+ mouse model of MFS (Habashi et al., 2006). Contradicting these results, others have shown that the aortopathy in the Fbn1C1039G/+ mouse model is not mediated by excessive TGF-β signaling and in fact is exacerbated by loss of TGF-β signaling in smooth muscle cells (Wei et al., 2017). Furthermore, TGF-β signaling protects against abdominal aortic aneurysms in angiotensin II-infused mice (Angelov et al., 2017). This controversy emphasizes the importance of understanding exactly how TGF-β signaling is impacted in MFS. Furthermore, the related syndrome LDS is caused by pathogenic mutations in several different components of the TGF-β pathway, TGFBR1, TGFBR2, SMAD2, SMAD3, and the ligands, TGFB2 and TGFB3. These mutations all cause missense amino acid substitutions that have been either verified in vitro, or are predicted to be loss of function, implying that LDS is caused by attenuated TGF-β signaling (Horbelt et al., 2010; Cardoso et al., 2012; Schepers et al., 2018). However, paradoxically, histological and biochemical studies of aortic tissue derived from LDS patients reveal an apparent high TGF-β signaling signature (van de Laar et al., 2012; Gallo et al., 2014; Lindsay et al., 2012). SGS is caused by mutations in SKI, and both SMAD-mediated and non-SMAD-mediated TGF-β signaling has been reported to be increased in primary dermal fibroblasts from SGS patients (Doyle et al., 2012).
The co-repressors SKI and SKIL play important roles in a number of different cellular processes including proliferation, differentiation, transformation, and tumor progression (Bonnon and Atanasoski, 2012). They are dimeric proteins that interact with both phosphorylated SMAD2 and SMAD3 (PSMAD2 or PSMAD3) via short motifs at their N-termini, and with SMAD4 via a SAND domain (named after Sp100, AIRE-1, NucP41/75, DEAF-1) in the middle of both proteins (Deheuninck and Luo, 2009). Between these two domains lies a Dachshund homology domain (DHD), which is thought to also be important for R-SMAD binding (Wilson et al., 2004; Ueki and Hayman, 2003). SKI and SKIL both contain a leucine zipper domain in their C-termini, through which they dimerize (Deheuninck and Luo, 2009). They are negative regulators of TGF-β/Activin signaling, with two distinct mechanisms of regulation having been proposed. In one model, SKI and SKIL bind with SMAD4 to SMAD binding elements (SBEs) of TGF-β/Activin target genes, and recruit co-repressors such as NCOR1 or SIN3A (Tokitou et al., 1999; Nomura et al., 1999; Stroschein et al., 1999; Deheuninck and Luo, 2009). They thus maintain the transcription of these target genes suppressed in the absence of signal. Upon TGF-β/Activin signaling, SKI and SKIL are rapidly degraded by the E3 ubiquitin ligase, RNF111 (formerly known as Arkadia), a process that requires SKI/SKIL binding to PSMAD2 or PSMAD3 (Le Scolan et al., 2008; Levy et al., 2007; Nagano et al., 2007). This then allows the activated SMAD3–SMAD4 complexes to bind the exposed SBEs and activate target gene transcription (Levy et al., 2007; Stroschein et al., 1999). In the competing model, SKI and SKIL act as repressors of active signaling simply by binding to PSMAD2 or PSMAD3 and SMAD4 in such a way as to disrupt the activated PSMAD2/PSMAD3–SMAD4 complexes (Luo, 2004; Ueki and Hayman, 2003; Wu et al., 2002). The heterozygous missense mutations that cause SGS have been mapped in SKI to the N-terminal R-SMAD-binding domain, with some small deletions and point mutations also found in the DHD, which is also necessary for R-SMAD binding (Carmignac et al., 2012; Doyle et al., 2012; Schepers et al., 2015). Thus, depending on the mechanism whereby SKI inhibits TGF-β/Activin signaling, loss of the interaction with PSMAD2/PSMAD3 would be predicted to have opposite effects on signaling output. If the PSMAD2/PSMAD3 interaction is required for SKI degradation, its loss would inhibit TGF-β signaling. However, if SKI binding to PSMAD2/PSMAD3 disrupts active SMAD complexes, then its loss would promote TGF-β signaling.
Here we use a combination of genome editing, structural biology, biochemistry, and analysis of patient samples to elucidate the molecular mechanism underlying SGS and to resolve the paradox surrounding the role of TGF-β signaling in Marfan-related syndromes. We first determine at the molecular level how SKI/SKIL function in the TGF-β/Activin signaling pathways and show that an intact ternary phosphorylated R-SMAD–SMAD4 complex is required for ligand-induced SKI/SKIL degradation. We demonstrate that the SGS mutations in SKI abolish interaction with PSMAD2 and PSMAD3 and this results in an inability of SKI to be degraded in response to TGF-β/Activin signaling. We go on to show that SKI stabilization results in an attenuation of the TGF-β transcriptional response in both knockin HEK293T cells and fibroblasts from SGS patients. Our work unequivocally establishes that SGS mutations lead to an attenuated TGF-β response, which has major implications for all the Marfan-related syndromes.
Results
A PSMAD2/3–SMAD4 ternary complex is essential for TGF-β/Activin-induced degradation of SKI/SKIL
To understand the consequences of SKI mutations in SGS and to resolve the paradox surrounding the function of TGF-β signaling in Marfan-related syndromes, we first set out to determine exactly how SKI and SKIL act as negative regulators of TGF-β and Activin signaling. We and others have previously demonstrated that SKI and SKIL are rapidly degraded upon TGF-β/Activin stimulation by the E3 ubiquitin ligase RNF111, and this requires PSMAD2 or PSMAD3 (Le Scolan et al., 2008; Levy et al., 2007; Nagano et al., 2007). Knockdown experiments suggested that SMAD4 was not necessary (Levy et al., 2007), but we subsequently showed that tumor cells deleted for SMAD4 or containing mutations in SMAD4 that abolish interactions with activated R-SMADs, abrogated TGF-β-induced degradation of SKI/SKIL (Briones-Orta et al., 2013). Whether the requirement for SMAD4 was direct or indirect was not clear.
To define the role of SMAD4 in SKI/SKIL degradation, we used CRISPR/Cas9 technology to delete SMAD4 in transformed embryonic kidney cells HEK293T, which express both SKI and SKIL, and in the human keratinocyte cell line, HaCaT, which predominantly express SKIL (Levy et al., 2007; Figure 1—source data 1). In wild-type (WT) cells, TGF-β/Activin induced rapid SKI and SKIL degradation, compared to cells treated with the TGFBR1 inhibitor, SB-431542 (Inman et al., 2002; Figure 1A,B). Deletion of SMAD4 in multiple clones of both cell types abolished ligand-induced SKI/SKIL degradation (Figure 1A,B). We validated these SMAD4-null cell lines by demonstrating that transient expression of SMAD4 could rescue TGF-β/Activin-induction of the SMAD3–SMAD4 reporter, CAGA12-Luciferase (Figure 1—figure supplement 1A,B). Furthermore, we could show that loss of SMAD4 inhibited the ligand-induced expression of a number of endogenous TGF-β and BMP target genes (Figure 1—figure supplement 1C). By knocking out SMAD2 or SMAD3 individually or together, we also confirmed that these R-SMADs are absolutely required for TGF-β/Activin-induced degradation of SKI and SKIL and act redundantly (Figure 1C; Figure 1—source data 1). Thus, R-SMADs and SMAD4 are all essential for TGF-β/Activin-dependent SKI/SKIL degradation.
In addition to forming a ternary complex with PSMAD2 or PSMAD3, SMAD4 has also been shown to interact directly with SKI and SKIL through their SAND domains (Walldén et al., 2017; Wu et al., 2002). To determine which of these SMAD4 interactions were important for TGF-β/Activin-induced SKI/SKIL degradation, we stably reintroduced enhanced GFP (EGFP) fusions of WT or mutated SMAD4 into HaCaT SMAD4-null cells. We selected two missense mutations on opposite faces of the C-terminal Mad homology 2 (MH2) domain of SMAD4: Asp351->His (D351H) and Asp537-->Tyr (D537Y) (Shi et al., 1997). These have been shown to occur naturally in the human colorectal cancer cell lines CACO-2 and SW948, and have lost the ability to bind phosphorylated R-SMADs (De Bosscher et al., 2004). In addition, we used the crystal structure of the MH2 domain of SMAD4 and the SAND domain of SKI, to design two mutations Ala433-->Glu (A433E) and Ile435-->Tyr (I435Y), that would be expected to abolish SMAD4 binding to SKI and SKIL (Wu et al., 2002).
We confirmed that these SMAD4 mutants behaved as expected in the rescue cell lines by testing their interaction with SKIL and R-SMADs by immunoprecipitation. As endogenous RNF111 triggers SKIL degradation in TGF-β/Activin-dependent manner, the stable SMAD4-expressing HaCaT rescue cell lines were incubated with the proteasome inhibitor, MG-132 for 3 hr prior to TGF-β stimulation, to block SKIL degradation. As predicted, the D351H and D537Y SMAD4 mutants had lost their ability to bind SMAD2 upon TGF-β induction, but retained the interaction with SKIL. By contrast, A433E and I435Y SMAD4 mutants were unable to bind SKIL, but could interact with SMAD2 upon TGF-β stimulation (Figure 2A). Furthermore, as expected, D351H and D537Y SMAD4 mutants failed to rescue the ability of TGF-β to induce expression of CAGA12-Luciferase in HaCaT SMAD4-null cells or rescue TGF-β-induced transcription of target genes, but the A433E and I435Y SMAD4 mutants rescued these responses almost as well as WT SMAD4 (Figure 2—figure supplement 1A,B).
Having demonstrated that these mutants behaved as designed, we asked which were able to mediate TGF-β-induced SKIL degradation, using three different assays. In a Western blot assay using nuclear extract, we found that reintroduction of WT SMAD4 in SMAD4-null cells caused a 50% reduction in SKIL levels in TGF-β-induced cells compared to those treated with SB-431542 (Figure 2B). However, none of the four SMAD4 mutants could rescue TGF-β-induced SKIL degradation (Figure 2B). We then established a flow cytometry assay to quantify SKIL protein stability in EGFP/EGFP-SMAD4-expressing cells (Figure 2C; Figure 2—figure supplement 1C). Treatment with TGF-β for 1 hr caused a 52% reduction in the relative median fluorescence intensity in the EGFP-SMAD4 WT-expressing cells, reflecting SKIL levels, compared to cells treated with SB-431542 (Figure 2C). However, for all four SMAD4 mutants tested, the median fluorescence was not decreased by TGF-β treatment (Figure 2C). Finally, we used an immunofluorescence analysis to monitor SKIL protein stability following TGF-β exposure. SMAD4-null cells showed strong nuclear staining of SKIL in the non-signaling condition (SB-431542), which remained unchanged by TGF-β treatment (Figure 3). Reintroduction of WT EGFP-SMAD4 conferred the ability to degrade SKIL upon TGF-β treatment, whereas none of the mutant SMAD4s were able to rescue SKIL degradation (Figure 3, arrows). Thus, all three assays demonstrate that a ternary R-SMAD–SMAD4 complex is absolutely necessary for TGF-β-induced SKIL degradation, as is the ability of SMAD4 to interact with SKIL itself. This suggests that within a canonical activated ternary SMAD complex, the R-SMAD component binds to the N-terminal region of SKIL/SKI, whilst SMAD4 binds the SAND domain, and both interactions are absolutely required for SKIL/SKI degradation.
SGS mutations inhibit the interaction of SKI with phosphorylated R-SMADs
We next investigated the consequences of the SGS mutations on SKI and SKIL’s ability to interact with the R-SMADs. SKI and SKIL share a highly conserved region at their N-terminus comprising the domain known to be important for R-SMAD binding (Deheuninck and Luo, 2009; Figure 4—figure supplement 1A). We first determined the minimal region of SKI required for R-SMAD binding using peptide pulldown assays with biotinylated SKI peptides and whole-cell extract from uninduced and TGF-β-treated HaCaT cells. This revealed that amino acids 11–45 of SKI are sufficient for binding to PSMAD2 and PSMAD3 upon TGF-β stimulation, whilst the unphosphorylated SMADs did not bind to any of the SKI peptides (Figure 4—figure supplement 1B). SMAD4 is also pulled down in these assays in a ligand-induced manner, by virtue of its interaction with the phosphorylated R-SMADs.
The SGS mutations discovered so far mostly cluster within this 11–45 region of SKI and a few deletions and point mutations have additionally been mapped in the DHD domain (Carmignac et al., 2012; Doyle et al., 2012; Schepers et al., 2015). The residues mutated are completely conserved, both between species, and also in the related protein, SKIL (Figure 4—figure supplement 1A; Carmignac et al., 2012). To determine the effect of these mutations on R-SMAD interaction, we introduced six different SGS mutations into the SKI peptide 11–45 and showed that they all prevented binding of PSMAD2 and PSMAD3, and as a result, also SMAD4 (Figure 4A). These results were also confirmed with the equivalent mutations in SKIL (Figure 4B). We proved that the interaction with SMAD2 was mediated via its MH2 domain using a mouse embryonic fibroblast cell line that expresses a truncated SMAD2 protein comprising just the MH2 domain (Piek et al., 2001; Das et al., 2009; Figure 4C). We confirmed this using recombinant human phosphorylated SMAD2 MH2 domain produced in insect cells by co-expressing the SMAD2 MH2 domain with the kinase domain of TGFBR1 (Figure 4D). In both cases, the SGS mutations prevented interaction of the SKI peptide with the SMAD2 MH2 domain.
We next used a peptide array to gain a better understanding of which amino acids can be tolerated at the positions found to be mutated in SGS and to determine which other amino acids in this region of SKI are essential for the R-SMAD interaction. The SKI peptide corresponding to amino acids 11–45 was synthesized as an array on a cellulose sheet such that each residue in the sequence between residues 19 and 35 was substituted with all 19 alternative amino acids (Figure 4E; Figure 4—source data 1). The array was probed with a recombinant PSMAD3–SMAD4 trimer, generated by co-expressing SMAD3 and SMAD4 with the TGFBR1 kinase domain in insect cells. The PSMAD3–SMAD4 complex was then detected using a fluorescently-labeled SMAD2/3 antibody. Eight residues are intolerant to almost any amino acid substitution (Thr20, Leu21, Phe24, Ser28, Ser31, Leu32, Gly34, and Pro35). Strikingly, six of these residues are the amino acids known to be mutated in SGS patients, and the array results readily explain why these residues are mutated to a number of different amino acids in SGS (Figure 4E; for quantification, see Figure 4—figure supplement 1C and Figure 4—source data 2). In addition, Thr20 and Phe24 are also crucial residues for binding the PSMAD3–SMAD4 complex, but have not yet been reported as disease mutations. Mutations in the other nine amino acids do not impair the binding, and almost any other amino acid apart from proline can be tolerated at these positions.
Crystal structure of the SKI peptide with the phosphorylated SMAD2 MH2 domain
To discover why these eight amino acids were so crucial for R-SMAD binding, and also to understand why SKI and SKIL only recognize phosphorylated R-SMADs, we solved the crystal structure of the SKI peptide (amino acids 11–45) with a phosphorylated homotrimer of the SMAD2 MH2 domain, produced in insect cells as described above. We confirmed using SEC-MALLS that the phosphorylated SMAD2 MH2 domain was indeed trimeric in solution (Figure 5—figure supplement 1A). Analysis of the binding affinity of the SKI peptide to the SMAD2 MH2 domain trimer indicated that the dissociation constant (Kd) was in the low nanomolar range (Figure 5—figure supplement 1B). The structure was determined by molecular replacement and refined at 2 Å resolution and readily explained why the crucial amino acids identified in the peptide were required for SMAD2 binding (Figure 5; Figure 5—figure supplement 1C).
The SKI peptide binds on the outside face of the MH2 domain at the so-called three helix bundle, comprising helices 3, 4, and 5 (Wu et al., 2001; Figure 5A). The N-terminal helix of SKI packs against helix 3 of SMAD2, and the C-terminal portion of the SKI peptide, which contains the critical Gly34 and Pro35, forms a sharp turn that is stabilized by pi-stacking coordination between Phe24 of SKI, Trp448 of SMAD2, and Pro35 of SKI (Figure 5B). Moreover, the NE1 of the Trp448 side chain forms a H-bond to the main chain carbonyl group of Gly33, which in turn positions Pro35 for the interaction with Trp448 (Figure 5B). Furthermore, Glu270 in SMAD2 provides a pocket, which has a negatively charged base that ties down SKI Gly34 through hydrogen-bonding to its main chain amides. Other key interactions involving amino acids identified above as crucial for binding include the main chain carbonyl of Ser31, which forms a hydrogen bond to the ND1 of Asn387 in helix 3 (Figure 5C), and the hydroxyl group of SKI Thr20, which forms a hydrogen bond with the Gln455 at the end of helix 5 of SMAD2, and is nearly completely buried in the interface (Figure 5D). The two leucine residues (Leu21 and Leu32) that are mutated in SGS are both buried in the structure (Figure 5E,F). The structure we obtained is consistent with a SKI–SMAD2 MH2 domain structure that was published by others, while this work was in progress (Miyazono et al., 2018). In that case, a pseudo-phosphorylated SMAD2 MH2 domain produced in Escherichia coli was used, complexed with a SKI peptide containing a C-terminal acidic tag (Ser-Asp-Glu-Asp).
Since our structure was generated with phosphorylated SMAD2, we were able to explore why SKI only binds phosphorylated R-SMADs and not monomeric unphosphorylated R-SMADs. To do this we compared the structure of the unphosphorylated SMAD2 MH2 domain bound to a region of ZFYVE9 (formerly called SARA) (Wu et al., 2000) with our current structure of phosphorylated SMAD2 MH2 domain complexed with SKI. It was clear that in the unphosphorylated SMAD2 structure, Tyr268 in the so-called β1’ strand (amino acids 261–274) is locked in a stable conformation in a hydrophobic pocket, and also forms a number of hydrogen bonds (Figure 5G). Crucially, this conformation forces Trp448 into flattened orientation, which is incompatible with SKI binding through the pi-stacking involving SKI Phe24, SMAD2 Trp448, and SKI Pro35 (Figure 5G). MH2 domain trimerization generates a new binding site for the β1’ strand on the adjacent MH2 domain subunit (Figure 5H; Video 1). The central residue driving this is Tyr268. In the trimer, the hydroxy group of Tyr268 makes hydrogen bond contact with the carbonyl group of Asp450 and the main chain of Lys451 on the adjacent MH2 domain subunit. As a consequence, Trp448 moves into an upright position in the trimer, allowing engagement with SKI. Thus, SKI can only bind SMAD2 in its phosphorylated trimeric state.
Knockin of an SGS mutation into HEK293T cells inhibits Activin-induced SKI degradation and attenuates PSMAD3–SMAD4-mediated transcriptional activity
We have shown that the presence of SGS mutations prevent the interaction of SKI/SKIL with phosphorylated R-SMADs and have demonstrated that SKI/SKIL degradation requires an activated R-SMAD–SMAD4 complex. We therefore went on to investigate the functional effect of the SGS mutations on SKI degradation and TGF-β/Activin-induced transcriptional responses. To do this, we chose to focus on Pro35 because of its crucial role in forming the stacking interaction with Trp448 in SMAD2. In SGS patients, Pro35 is mutated to Ser or Gln (Carmignac et al., 2012; Doyle et al., 2012; Schepers et al., 2015), substitutions not tolerated in activated R-SMAD–SMAD4 binding (Figure 4E).
We used CRISPR/Cas9 technology with a single-stranded template oligonucleotide to knock in the Pro35 → Ser (P35S) mutation into HEK293T cells, and we efficiently generated a number of homozygous clones (Figure 6—figure supplement 1A). In three independent clonal cell lines carrying the P35S SKI mutation, the binding to endogenous phosphorylated SMAD2 was severely compromised, compared with WT SKI (Figure 6A). The binding to SMAD4, however, was unchanged in the mutant cell lines, as the SGS mutations do not affect the SKI SAND domain, which is responsible for SMAD4 binding (Figure 6A). To assess the impact of the P35S SKI mutation on the SKI and SKIL degradation, cells were treated with Activin for 1 or 2 hr and SKI/SKIL levels determined by immunoblotting. At both time points, we clearly demonstrated that P35S SKI levels remained stable, whilst in the parental cell lines, SKI protein is almost entirely degraded after 1 hr of Activin treatment (Figure 6B). The presence of mutated SKI had no effect on the Activin-induced degradation of SKIL in these lines (Figure 6B). Thus, the P35S mutation renders SKI completely resistant to ligand-induced degradation.
To determine whether SGS mutations had the same effect in SKIL, we introduced a G103V mutation into SKIL, corresponding to the SGS mutation G34V in SKI (referred to as SKIL ΔS2/3). Transfection of G103V SKIL in HEK293T cells led to reduction of SMAD2 binding in parental cells (Figure 6—figure supplement 1B). The residual binding was mediated via SMAD2’s interaction with SMAD4, as it was lost in the SMAD4 knockout cells (Figure 6—figure supplement 1B). Binding of SMAD4 in the absence or presence of signal was unaffected by the mutation. As observed above for SKI, the SGS mutation in SKIL led to resistance to Activin-induced degradation (Figure 6—figure supplement 1C), indicating that the R-SMAD interaction was essential. In addition, we made a version of SKIL with mutations in the SAND domain (R314A, T315A, H317A, and W318E) that rendered it unable to interact with SMAD4 (referred to as SKIL ΔS4). This mutant was also not degraded upon Activin stimulation (Figure 6—figure supplement 1C), demonstrating an essential requirement for SMAD4 binding.
SKI and SKIL bind DNA in conjunction with SMAD4 at SBEs of TGF-β/Activin target genes in the absence of ligand stimulation. The ligand-induced degradation of SKI and SKIL then allows the activated R-SMAD–SMAD4 complexes access to the SBEs to activate transcription of target genes (Levy et al., 2007; Stroschein et al., 1999). We hypothesized that if the SGS mutations render SKI resistant to ligand-induced degradation, then mutant SKI and SMAD4 would remain bound to the DNA. To test this, we used a DNA pulldown assay with an oligonucleotide corresponding to SBEs from the JUN promoter (Levy et al., 2007). Consistent with our prediction, both WT and P35S SKI bound the SBEs with SMAD4 in the absence of signal, but after Activin stimulation, the binding of WT SKI is lost, whilst the binding of P35S SKI is retained (Figure 6C).
The SKI–SMAD4 complex bound to SBEs in the absence of signal are transcriptionally repressive (Levy et al., 2007; Stroschein et al., 1999). Since P35S SKI remains bound with SMAD4 in Activin-stimulated cells, we reasoned that this would inhibit Activin-induced gene expression. To address this, we stably expressed luciferase reporters (CAGA12-Luciferase and BRE-Luciferase), together with TK-Renilla as an internal control (Dennler et al., 1998; Korchynskyi and ten Dijke, 2002), in parental HEK293T cells and in two independent clones of the knockin P35S SKI cells. The CAGA12-Luciferase reporter responds to TGF-β and Activin, is induced by PSMAD3–SMAD4 complexes, and is sensitive to SKI and SKIL levels, whilst the BRE-Luciferase reporter is induced by SMAD1/5–SMAD4 complexes in response to BMPs and is not affected by SKI and SKIL (Levy et al., 2007). Strikingly, we found a significant reduction in Activin-induced CAGA12-Luciferase activity in the P35S SKI cells compared to the parental cell line (Figure 6D), while BMP4-induced BRE-Luciferase activity was similar in all cell lines (Figure 6E). The results indicate that SGS mutations in SKI lead to inhibition of TGF-β/Activin-induced transcription mediated by PSMAD3–SMAD4 complexes.
Dermal fibroblasts from SGS patients exhibit an attenuated transcriptional TGF-β response
To gain further insights into the functional consequences of the SGS mutations, we obtained dermal fibroblasts from two SGS patients: a female carrying the heterozygous point mutation L32V and a male carrying a heterozygous deletion of 12 base pairs corresponding to codons 94–97 (ΔS94-97) (Carmignac et al., 2012). Both patients present the classical features of SGS such as marfanoid habitus and intellectual disability, whilst the patient-carrying L32V mutation also manifests craniosynostosis. In addition, we obtained dermal fibroblast from a healthy male subject as a control.
We investigated whether the SGS mutations rendered SKI resistant to TGF-β-induced degradation in the dermal fibroblasts, as demonstrated in the knockin HEK293T cells. Indeed, in control fibroblasts, SKI expression is abrogated upon TGF-β stimulation after 1 hr, while in the SGS-derived fibroblasts, SKI protein remains relatively stable (Figure 7A). Note that these cells are heterozygous for the mutation, while the HEK293T P35S SKI knockins used above were homozygous, thus accounting for the incomplete SKI stabilization exhibited by the SGS fibroblasts, compared with the knockin cells.
To determine the effect of SKI stabilization on global TGF-β-induced transcription we performed genome-wide RNA-sequencing (RNA-seq) in three different conditions: SB-431542-treated cells (non-signaling condition), 1 hr TGF-β-treated and 8 hr TGF-β-treated, and compared the L32V and ΔS94-97 SKI fibroblasts to control fibroblasts. The samples separated in a principal component analysis according to the cell line used and the treatment performed (Figure 7—figure supplement 1A), and we confirmed that the differentially enriched genes after TGF-β treatment in the control fibroblasts were characteristic of pathways related to TGF-β signaling (Figure 7—figure supplement 1B). We then performed a pairwise comparison between the TGF-β-treated and SB-431542-treated samples for each of the cell lines individually. We found that of the 339 genes that were differentially expressed in normal fibroblasts after 1 hr of TGF-β treatment, 60% (202 genes) were induced or repressed less efficiently in the L32V mutant fibroblasts, and of these, 97 genes were not significantly differentially expressed in the mutant cells at all (Figure 7—source data 1, Figure 7B,C). After an 8 hr TGF-β induction, we found that 4769 genes were differentially expressed in the normal fibroblasts, and of these 75% (3556 genes) were induced or repressed less efficiently in the L32V mutant cells, and 880 genes were not significantly differentially expressed in the mutant cells at all (Figure 7—source data 1, Figure 7B,C). We observed similar results when comparing the TGF-β responses in the normal fibroblasts versus ΔS94-97 SKI fibroblasts, although the effects were less dramatic (Figure 7—source data 1, Figure 7—figure supplement 1C,D). To illustrate the magnitude of these effects we validated six gene expression profiles (ISLR2, CALB2, SOX11, ITGB6, HEY1, COL7A1) in the normal versus mutant fibroblasts by qPCR (Figure 7—figure supplement 2).
Thus, we conclude that the presence of SGS point mutations in SKI that render it resistant to ligand-induced degradation, result in attenuated TGF-β responses for a substantial subset of target genes.
Discussion
SGS mutations in SKI lead to stabilization of SKI and attenuated TGF-β/Activin transcriptional responses
In this study, we have resolved the mechanism of action of SKI and the related protein SKIL, which has allowed us to elucidate the molecular consequences of SKI mutations in SGS. Our proposed mechanism of how SKI acts as a transcriptional repressor of TGF-β/Activin signaling in health and disease is illustrated in Figure 7D, and our results suggest that the mechanism of action of SKIL is equivalent. In the absence of ligand stimulation, in both healthy and diseased cells, SKI and SKIL bind in conjunction with SMAD4 at SBEs of TGF-β/Activin target genes. Here, they repress transcription by recruiting corepressors such as NCOR1 or SIN3A (Tokitou et al., 1999; Nomura et al., 1999; Stroschein et al., 1999; Deheuninck and Luo, 2009). We and others have also shown that in unstimulated cells, SKI and SKIL interact with the E3 ubiquitin ligase, RNF111, although this binding per se does not lead to SKI/SKIL degradation (Le Scolan et al., 2008; Levy et al., 2007; Nagano et al., 2007). In healthy cells, upon TGF-β/Activin stimulation, SKI/SKIL form a complex with a canonical PSMAD2/PSMAD3–SMAD4 trimer, which induces rapid degradation of SKI/SKIL via RNF111 (Le Scolan et al., 2008; Levy et al., 2007; Nagano et al., 2007). A possible mechanism explaining this would be that the binding of the PSMAD2/PSMAD3–SMAD4 complex induces an activating conformational change in RNF111, although this has not yet been demonstrated. Degradation of SKI and SKIL removes the repressors from the SBEs, allowing access of activated PSMAD3–SMAD4 complexes to the SBEs to regulate transcription of target genes. In SGS cells, mutated SKI can no longer interact with PSMAD2/PSMAD3, and it is therefore not degraded upon TGF-β/Activin signaling. It thus remains bound with SMAD4 to SBEs, resulting in an attenuation of TGF-β/Activin transcriptional responses. To demonstrate the functional consequences of the SGS mutations, we used genome-wide RNA-seq analysis of fibroblasts derived from SGS patients and have shown that SGS mutations indeed lead to a reduction in the magnitude of TGF-β transcriptional responses.
The mechanism underlying SKI and SKIL function
Since the discovery that SKI and SKIL interact with SMAD2 and SMAD3, and act as negative regulators of TGF-β/Activin pathways (Luo et al., 1999; Stroschein et al., 1999; Sun et al., 1999), two different mechanisms of action have been proposed. One mechanism is as described in the paragraph above – an initial version of which was first proposed in 1999 (Stroschein et al., 1999). The second mechanism was based on the crystal structure of the SKI SAND domain with the MH2 domain of SMAD4 (Wu et al., 2002). In this crystal structure, the binding of SKI with SMAD4, which is mediated via the I-loop of the SKI SAND domain with the L3 loop of SMAD4, was mutually exclusive with the binding of SMAD4 to activated R-SMADs, which also requires the L3 loop of SMAD4. Thus, the authors concluded that the mechanism whereby SKI (and by analogy, SKIL) inhibited TGF-β/Activin signaling was by binding to the activated R-SMADs and SMAD4 in such a way as to disrupt the phosphorylated R-SMAD–SMAD4 complexes required for transcriptional activation (Wu et al., 2002). This mechanism has been supported by the observation that overexpression of SKI and SKIL inhibits TGF-β-induced functional responses (Luo, 2004). The two mechanisms are fundamentally different. In the first, SKI and SKIL are constitutive repressors that need to be degraded to allow pathway activation. In the second, SKI and SKIL act as inducible repressors, as they repress only upon ligand induction, by virtue of their ability to disrupt activated SMAD complexes.
Our biochemical analysis of the role of SMAD4 in SKI/SKIL function now resolves the controversy between the two models. First, the efficient SKIL and SKI degradation that we and others have observed upon ligand stimulation in effect rules out the second model, at least in the first hours after ligand induction, as there would be little or no nuclear SKI or SKIL to disrupt activated SMAD complexes. Second, we clearly demonstrate for the first time that an intact functional phosphorylated R-SMAD–SMAD4 trimer is required to bind to SKIL to induce its ligand-dependent degradation. The key piece of evidence for this comes from our analysis of the inability of SMAD4 point mutants to restore ligand-induced SKIL degradation in SMAD4-null HaCaTs. Critically, we show that SMAD4 mutants that cannot form a canonical activated R-SMAD–SMAD4 trimer cannot rescue ligand-induced SKIL degradation, neither can SMAD4 mutants that do not interact with the SAND domain of SKIL. Thus, we conclude that TGF-β/Activin-induced SKIL degradation occurs only when SKIL interacts simultaneously with phosphorylated R-SMADs and with SMAD4, which in turn must interact with each other in a transcriptionally active trimer. This strongly indicates that SKIL binding to the R-SMAD–SMAD4 complex does not disrupt it.
This conclusion is supported by a recent crystal structure of the SAND domain of SKIL with the SMAD4 MH2 domain (Walldén et al., 2017). This structure revealed that SKIL interacts with SMAD4 in two states: an ‘open’ and a ‘closed’ conformation. In the open conformation, the authors showed that SKIL can bind the R-SMAD–SMAD4 complex without intermolecular clashes or further structural readjustment, whereas in the closed state, structural reorganization within the SMAD heterotrimer is required to allow binding of SKIL, as has been observed in the previous structure of SKI with SMAD4 (Wu et al., 2002). Molecular modelling has subsequently confirmed that SKIL in the open conformation forms a stable ternary SKIL–SMAD3–SMAD4 complex (Ji et al., 2019). Furthermore, surface plasmon resonance indicated only one dominant binding mode for SKIL and SMAD4, leading to the conclusion that the open conformation is the biologically and functionally relevant mode and that the closed conformation may be the result of crystal packing forces (Walldén et al., 2017). Note that the residues that allow the binding in the open conformation are highly conserved between SKI and SKIL, suggesting that both repressors bind to an intact activated R-SMAD–SMAD4 complex, which is required for their degradation. This would exclude the disruption model and thus favors the degradation model. It will now be important to solve the structure of an activated R-SMAD–SMAD4 trimer with the N-terminal half of SKI or SKIL that contains the R-SMAD binding motif, the DHD domain and the SAND domain that contacts the SMAD4 moiety. The role of the DHD domain is particularly intriguing as SGS mutations also occur in this domain (Carmignac et al., 2012; Doyle et al., 2012; Schepers et al., 2015) and appear from our patient sample analysis to have a similar effect on inhibiting ligand-induced SKI degradation.
Our structural data also elegantly explain why SKI and SKIL only bind to phosphorylated SMAD2 and SMAD3 in the context of an activated SMAD trimer, and not to monomeric SMAD2 or SMAD3. The key residue for this discrimination is Trp448 in SMAD2, the equivalent of which would be Trp406 in SMAD3. In the trimer, this residue is in a conformation compatible with stacking between Phe24 and Pro35 of SKI. In the monomer, however, it is rotated approximately 90°, prohibiting SKI binding. The binding mode of SKI to SMAD2 is distinct from that of other SMAD2-binding partners, for example, the transcription factor FOXH1, which contains two binding motifs (Miyazono et al., 2018). One of which, the so-called SIM, binds SMAD2 in both monomeric and trimeric forms, whilst the so-called FM only binds phosphorylated trimeric SMAD2 because it recognizes the interface of the SMAD trimer (Miyazono et al., 2018; Randall et al., 2002; Randall et al., 2004). It will be interesting in the future to discover whether any other SMAD2/3-binding partner uses the same mode of interaction as SKI and SKIL.
The molecular mechanism of SGS
As discussed in Introduction, SGS is a Marfan-related syndrome, with patients exhibiting many of the same features characteristic of MFS and LDS. These syndromes have been considered as TGF-β signalopathies, as the causal mutations are either direct components of the TGF-β signaling pathway or, as in the case of MFS, a component of the microfibrils in the ECM, FBN1, that is known to bind latent TGF-β in complex with latent TGF-β binding proteins (LTBPs) (Cannaerts et al., 2015; Ramirez et al., 2004; Robertson et al., 2015). There has been controversy over whether the manifestations of these syndromes result from too little TGF-β signaling, or too much. This is obviously a crucial issue to resolve, as it is influencing the types of treatments being developed for patients with these syndromes.
The first suggestion that MFS resulted from excessive TGF-β came from mouse models, where key phenotypes could be rescued by a TGF-β neutralizing antibody (Habashi et al., 2006; Neptune et al., 2003). However, later studies using a potent murine anti-TGF-β antibody, or genetic methods for reducing TGF-β signaling have not corroborated these findings, and have worsened, rather than improved disease in MFS mouse models (Cook et al., 2015b; Holm et al., 2011; Lindsay et al., 2012; Wei et al., 2017). Furthermore, administration of small-molecule inhibitors of the TGF-β type I receptor, or a pan-TGF-β neutralizing antibody, has been associated with serious adverse cardiovascular toxicities, such as valve defects, similar to those found in MFS (Anderton et al., 2011; Stauber et al., 2014; Mitra et al., 2020). The finding that mutations in FBN1 that prevent binding of LTBPs might result in lower levels of TGF-β signaling, rather than excessive signaling, is not so surprising given the more recent understanding of how TGF-β is activated. In order for mature TGF-β ligands to be released from the latent complex, either force has to be applied via integrins, to partially unfold the cleaved TGF-β pro-domain allowing release of the mature domain, or the pro-domain must be degraded by proteases (Dong et al., 2017; Rifkin et al., 2018; Robertson and Rifkin, 2016). For the traction mechanism to occur, the integrin must be anchored to the actin cytoskeleton, and the LTBP must be tethered to ECM, via FBN1 and fibronectin. Thus, release of latent TGF-β alone is not therefore sufficient to produce active TGF-β ligands.
Consistent with the view that lower levels of TGF-β signaling might be responsible for MFS, the mutations that give rise to LDS are all loss-of-function mutations in TGF-β pathway components (Schepers et al., 2018). Paradoxically though, signatures of higher TGF-β signaling were observed over time in mouse models of LDS (Gallo et al., 2014; MacFarlane et al., 2019). However, in this case, the pathology could not be rescued by neutralizing TGF-β activity. One possibility with both the mouse models of MFS and LDS is that the mutations do initially lead to lowered TGF-β signaling, but over time cells compensate by up-regulating either TGF-β ligands themselves or other TGF-β family ligands that signal through PSMAD2/PSMAD3, ultimately leading to the enhanced signaling signatures observed.
With respect to SGS, there has been much less research into the consequences of the SGS mutations on TGF-β signaling responses. Based on the SMAD complex disruption model of SKI/SKIL action, it has been assumed that loss-of-function mutations in a negative regulator would lead to an increase in TGF-β signaling, and in fact, this has been used to support the idea that the Marfan-related syndromes are caused by excessive TGF-β signaling (Doyle et al., 2012; Gallo et al., 2014). Here we unequivocally show that the opposite is true. We demonstrate that these mutations lead to loss of ligand-induced SKI degradation. As a result, the stabilized SKI remains bound to SBEs with SMAD4 as a repressive complex, and hence, a subset of TGF-β/Activin transcriptional responses are attenuated. We have proven this in both HEK293T knockin cells and patient-derived fibroblasts. Moreover, we also find no evidence for increased PSMAD2 or PSMAD3 signaling. Indeed, neither model of SKI function would actually predict that the SKI mutations would affect levels of phosphorylated R-SMADs, since SKI acts downstream of R-SMAD phosphorylation. Finally, our finding that SGS mutations in SKI lead to its stabilization and are not equivalent to loss of SKI function also explains why patients with 1p36 deletion syndrome, who are haploinsufficient for SKI, do not have SGS. However, unsurprisingly, many of the same organs are affected in both syndromes (Colmenares et al., 2002; Zhu et al., 2013).
It will now be important to use animal models to explore how attenuation of transcription of specific TGF-β/Activin target genes leads to the manifestations of SGS, and to understand why SGS patients exhibit additional defects compared with LDS and MFS patients. We anticipate that our new understanding that the SKI mutations lead to attenuation of TGF-β responses will resolve the paradoxes surrounding the role of aberrant TGF-β signaling in the other Marfan-related disorders and will help inform the development of new therapeutic approaches.
Materials and methods
Cell lines
Request a detailed protocolHEK293T and HaCaT cells were obtained from the Francis Crick Institute Cell Services and cultured in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% fetal bovine serum (FBS) and 1% Penicillin/Streptomycin (Pen/Strep). All CRISPR-Cas9 edited cell lines were cultured in the same media. Dermal fibroblasts from healthy subjects were kindly provided by David Abraham (UCL-Medical School Royal Free Campus) under the ethics of the Health Research Authority, NRES Committee London – Hampstead, Research Ethics Committee (REC) reference, 6398. L32V and ΔS94-97 SKI dermal fibroblasts were obtained from Laurence Faivre and Virginie Carmignac (Université de Bourgogne UMR1231 GAD, Dijon, France) under the ethics of the GAD collection, number DC2011-1332 (Carmignac et al., 2012). The mutations were confirmed by Sanger sequencing and RNA sequencing. The fibroblasts were all cultured in DMEM supplemented with 10% FBS, 1% Pen/Strep, and 1% insulin–transferrin–selenium (Thermo Fisher). Mouse embryo-derived fibroblasts harboring the homozygous null allele Smad2ex2 (MEF SMAD2ex2) (Piek et al., 2001) were maintained in DMEM supplemented with 10% FBS and 1% Pen/Strep. All cell lines have been banked by the Francis Crick Institute Cell Services and certified negative for mycoplasma. The identity of all cell lines was also authenticated by confirming that their responses to ligands and their phenotypes were consistent with published history. All the cell lines are listed in Key Resources Table.
Ligands, chemicals, and cell treatments
Request a detailed protocolLigands and inhibitors were used at the following concentrations: TGF-β (PeproTech), 2 ng/ml; Activin A (PeproTech), 20 ng/ml; BMP4 (PeproTech), 20 ng/ml; SB-431542 (Tocris), 10 µM; MG-132 (Tocris), 25 µM. All treatments were performed in full serum or, where required, in serum-starved (0.5% FBS) DMEM. Unless otherwise stated, cells were incubated with 10 µM SB-431542 overnight to inhibit autocrine signalling, then were washed three times with warm media, and stimulated with either Activin A or TGF-β. For proteasome inhibition, cells were treated for 3 hr with 25 µM MG-132 prior stimulation with Activin A or TGF-β.
Plasmids
Plasmids are listed in Key Resources Table. CAGA12-Luciferase, BRE-Luciferase, TK-Renilla, pEGFP-C1, and pEGFP-SMAD4 were as described previously (Dennler et al., 1998; Korchynskyi and ten Dijke, 2002; Levy et al., 2007; Nicolás et al., 2004). The pEGFP-SMAD4 mutants (D351H and D537Y) were generated by swapping the mutated SMAD4 coding region from the EF-HA vector into pEGFP-C1 (De Bosscher et al., 2004). The pEGFP-SMAD4 SMAD4 mutants (A433E and I435Y) were generated by PCR using oligonucleotides listed in Key Resources Table. EF-Flag-SKIL G103V was generated from pEF-FLAG-SKIL by PCR using oligonucleotides listed in Key Resources Table, while the mutant containing the R314A, T315A, H317A, and W318E mutations was generated by synthesizing the SKIL region between BSTEII and AVRII sites containing the mutations and cloning that fragment into pEF-FLAG-SKIL. For plasmids used for generating recombinant proteins, see below.
Transfections, generation of stable cell lines, and reporter assays
Request a detailed protocolCells were transfected with the appropriate plasmids using Fugene 6 (Roche) according to the manufacturer’s instructions. Luciferase reporter assays were performed as previously described, using the Dual-Glo assay system (Promega) following the manufacturer's instructions (Levy et al., 2007).
HaCaT SMAD4 KO lines stably expressing either EGFP or EGFP-SMAD4 WT or mutants were generated by transfecting the cells with the appropriate plasmids. Transfected cells were selected with 500 µg/ml of G418 (Invitrogen), then FACS sorted for EGFP-positive cells, and expanded. EGFP expression was confirmed by microscopy. To generate stable HEK293T cell lines expressing either CAGA12-Luciferase or BRE-Luciferase together with TK-Renilla, cells were transfected with the appropriate plasmids together with a plasmid carrying the puromycin resistance gene (pSUPER-retro-puro; OligoEngine). Cells were then selected with 2 μg/ml puromycin (Sigma).
CRISPR/Cas9-mediated knockout of SMADs in HEK293T and HaCaT cells
Request a detailed protocolFor the generation of knockin or knockout HEK293T cells, a parental clone was selected, which was a representative clone from the HEK293T pool that showed a robust Activin-induced SKI degradation and responded to TGF-β family ligands in the same way as the starting pool. For HaCaTs, a pool of cells was used as starting material for knockouts.
For SMAD2 and SMAD3 knockouts, a guide RNA in the MH2 domain of the protein was selected, whereas for SMAD4, two guide RNAs were picked, one targeting the MH1 domain and the other targeting at the end of the MH2 domain. The guide RNAs are shown in Key Resources Table. The guide RNAs were expressed from the plasmid pSpCas9(BB)−2A-GFP (pX458) (Addgene, #48138) (Ran et al., 2013). HEK293T and HaCaT were transfected with the appropriate plasmid, and for the double knockout SMAD2 and SMAD3, the two plasmids were transfected simultaneously. Forty-eight hours after transfection, cells were sorted for EGFP expression, plated as single cells in 96-well plates, and screened by Western blot to assess the loss of the protein. For HEK293Ts, two knockout clones for SMAD2, SMAD3, SMAD2/SMAD3, and SMAD4 were used in these studies. For HaCaTs, four independent SMAD4 knockout clones were used. The sequences of the knockout alleles are shown in Figure 1—source data 1.
Knockin of P35S SKI at the endogenous locus
Request a detailed protocolTo introduce the P35S mutation into SKI, a gRNA was selected immediately downstream of codon 35. A 120 bp ssODN, where the codon CCG (P35) was mutated to TCC (S35) and codons 33 and 34 where silently mutated from GGC to GGA, was made and purified by Sigma (Key Resources Table). The ssODN contained phosphorothioate bonds between the first two and last two nucleotides at the 5′ and 3′ ends, respectively, to avoid ssODN degradation by endogenous nucleases. The silent mutations at codons 33 and 34 were introduced to increase the specificity of the downstream screening primer. The mutation at codon 35 also disrupts the PAM sequence.
Cells were cotransfected with the pX458 plasmid expressing the gRNA and 10 μM ssODN using Fugene 6. After 48 hr, cells were FACS sorted for GFP expression and plated as single cells in 96-well plates. Subsequently, clones were consolidated, and from replicate plates, genomic DNA was extracted using Quickextract DNA extraction solution (Lucigen) according to manufacturer’s instructions. PCR was performed using a universal reverse primer and two different forward primers: a primer which allows detection of WT SKI and one that detects the P35S SKI (Key Resources Table). Clones positive for P35S SKI knockin were selected, verified by Sanger Sequencing, and used for further analysis.
Western blotting, immunoprecipitations, DNA pulldowns, and immunofluorescence
Request a detailed protocolWhole-cell extracts were prepared as previously described (Inman et al., 2002), while nuclear lysates were prepared according to Wong et al., 1999. Western blots were carried out using standard methods. The list of the antibodies used is shown in Key Resources Table. Immunoprecipitations using GFP-Trap beads (Cromotek) were performed according to the manufacturer’s instructions. Immunoprecipitations using antibodies coupled to protein G-Sepharose beads (Sigma) were as described previously (Levy et al., 2007).
DNA pulldown assays were performed as previously described with some modifications (Levy et al., 2007). Nuclear lysates were extracted using buffer containing 360 mM NaCl, and the DNA pulldowns were performed in the presence of a 40 μg of non-biotinylated mutant oligonucleotide to reduce non-specific binding. The oligonucleotides corresponding to WT and mutated SBE of the JUN promoter are shown in Key Resources Table. Immunofluorescence was performed as previously described (Pierreux et al., 2000), except that cells were washed and fixed for 5 min in methanol at −20 C (Levy et al., 2007). Nuclei were counter stained with DAPI (0.1 µ/ml). Imaging was performed on a Zeiss Upright 780 confocal microscope. Z-stacks were acquired for all channels, and maximum intensity projection images are shown.
For all of these techniques, a representative experiment of at least two biological repeats is shown.
Flow cytometry
Request a detailed protocolAfter 1 hr of TGF-β induction, HaCaT SMAD4 KO rescue cells, expressing either EGFP or EGFP-SMAD4 fusions, were washed, trypsinized, and pelleted. Cell pellets were fixed with methanol for 5 min at −20 C. Fixed cells were incubated with primary antibody against SKIL. Cells were then washed three times in phosphate-buffered saline (PBS) and incubated with secondary antibody conjugated with anti-rabbit Alexa 647. As a negative control, we used cells incubated with secondary antibody only. Antibodies used are listed in Key Resources Table. Subsequently, cells were washed three times with PBS, and pellets were resuspended with 500 µl of PBS and filtered to achieve a single-cell suspension. Cells were then analyzed for EGFP fluorescence on an LSRII flow cytometer (BD Biosciences), gated for viable, single cells. We then quantified the fluorescence emitted by Alexa 647 in cells expressing EGFP as a measure of SKIL protein. The FlowJo program was used to analyze the results.
Peptide pulldown assays and peptide array
Request a detailed protocolFor peptide pulldowns, N-terminal biotinylated peptides were synthesized by the Peptide Chemistry Facility at the Francis Crick Institute using standard procedures. The peptide pulldown assays were performed as described previously (Randall et al., 2002). Where recombinant protein was used, it was dissolved in buffer Y (50 mM Tris–HCl, pH 7.5, 150 mM NaCl, 1 mM EDTA, 1% [vol/vol] NP-40) and used at the concentrations given in the legend to Figure 4. Peptides sequences are given in Key Resources Table.
A peptide array was generated using peptides corresponding to SKI amino acids 11–45, of which the amino acids 19–35 were mutated one by one to every other amino acid. Arrays were synthesized on an Intavis ResPepSL Automated Peptide Synthesiser (Intavis Bioanalytical Instruments, Germany) on a cellulose membrane by cycles of N(a)-Fmoc amino acids coupling via activation of carboxylic acid groups with diisopropylcarbodiimide in the presence of hydroxybenzotriazole (HOBt), followed by removal of the temporary α-amino protecting group by piperidine treatment. Subsequent to chain assembly, side chain protection groups were removed by treatment of membranes with a deprotection cocktail (20 ml 95% trifluoroacetic acid, 3% triisopropylsilane, 2% water for 4 hr at room temperature) and then washed (4 x dichloromethane, 4 x ethanol, 2 x water, 1 x ethanol), prior to being air dried.
Peptide array membranes were blocked with 5% Milk in 0.01% Tween 20 in PBS and then incubated with purified PSMAD3–SMAD4 complex (see below) overnight. Subsequently, the membranes were washed and incubated with an antibody against SMAD2/3 (BD Biosciences) conjugated to Alexa 488 using a Zenon Mouse IgG Labeling Kit (Life Technology) according to the manufacturer’s instructions. Fluorescence was detected where binding occurred between SKI and PSMAD3–SMAD4 complex and measured with a Typhoon FLA 9500 biomolecular imager (GE Healthcare).
In all peptide experiments, a representative of at least two biological repeats is shown.
Generation of phosphorylated SMAD2 MH2 domain in insect cells
Request a detailed protocolThe cDNA encoding a fusion protein consisting of GST followed by a 3C cleavage site and the human SMAD2 MH2 domain (residues 241–465) was inserted into the MCSI of the pFastBac Dual vector (Thermo Fisher Scientific) in the Sal1 and Spe1 restriction sites. The cDNA encoding a constitutively active version of the human TGFBR1 kinase domain (residues 175–503 with a T204D mutation) was cloned into the MCSII using the SmaI and SphI restriction sites. A high-titre baculovirus (>108 pfu/ml) was generated using standard published protocols (Fitzgerald et al., 2006). Expression of the trimeric phosphorylated SMAD2 MH2 domain was performed by infecting Sf21 cells at a density of 1.5 × 106 cells/ml at a MOI: 1 and incubating for 72 hr at 27°C with rotation at 110 r.p.m. Cells were harvested by centrifugation at 1000 × g for 10 min and stored at −80°C until required.
Expression of the human phosphorylated SMAD3–SMAD4 complex
Request a detailed protocolAn expression construct was generated encoding GST-fused full-length human SMAD3 (with 3C protease site) and inserted into the MCSI of the pFastBac dual. As with the SMAD2 MH2 domain construct, the constitutively active human TGFBR1 kinase domain was cloned into MCSII. Another vector encoding full-length human SMAD4 alone was also constructed, where SMAD4 was inserted into the pBacPAK plasmid. The resulting vectors were used to generate high-titre virus (>108 pfu/ml) using a standard published protocol. For expression of the phosphorylated SMAD3–SMAD4 complex both viruses were used to infect cultures of Sf21 insect cells at a density of 1.5 × 106 cells/ml at a MOI: 1. Infected cultures were allowed to grow for 72 hr at 27°C with rotation at 110 r.p.m. Cells were harvested by centrifugation at 1000 x g for 10 min and stored at −80°C until required.
Purification of trimeric phosphorylated SMAD2 MH2 domain and phosphorylated SMAD3–SMAD4 complexes
Request a detailed protocolThe same procedure was used to purify both phosphorylated SMAD2 MH2 domain trimers and phosphorylated SMAD3–SMAD4 complexes. Typically, 500 ml of infected Sf21 cells were lysed in 30 ml of a lysis buffer consisting of 50 mM HEPES (pH 8.0), 250 mM NaCl, 10% (vol/vol) glycerol, 1% (vol/vol) Triton X-100, 10 mM β-glycerophosphate, 1 mM NaF, 10 mM benzamidine, and 1 mM dithiothreitol (DTT) supplemented with 5 µl BaseMuncher (2500 U/µl), 5 mM MgCl2, and phosphatase (phosphatase inhibitor cocktail 3, Sigma) and protease inhibitors (EDTA-free cOmplete protease inhibitors, Roche). After incubation for 20 min, the suspension was sonicated to ensure complete lysis. The insoluble fraction was pelleted by centrifugation (100,000 × g at 4°C for 30 min). The soluble fraction was incubated with 500 µl bed volume of Glutathione 4B Sepharose (Cytiva) for 2 hr at 4°C with gentle agitation. The resin was washed extensively with buffer containing 50 mM HEPES (pH 7.5), 200 mM NaCl, and 1 mM DTT. The phosphorylated SMAD2 MH2 domain or phosphorylated SMAD3–SMAD4 complexes were eluted from the GSH resin by cleavage with GST-3C protease (20 µg) in 5 ml wash buffer overnight at 4°C with gentle agitation. The proteins were concentrated to 0.5 ml and applied to a S200 10/300 Increase (Cytiva) size exclusion column equilibrated with 50 mM HEPES (pH 7.5), 200 mM NaCl, 5% (vol/vol) glycerol, and 1 mM DTT. Fractions were analyzed by SDS–PAGE and concentrated to 2 mg/ml and snap frozen as 50 µl aliquots and stored at −80°C until required. Purified phosphorylated SMAD2 MH2 domain was quantified using a molar extinction coefficient value of 39,420 M−1cm−1. The molar extinction coefficients used for SMAD3 and SMAD4 were 68,870 and 70,820 M−1cm−1, respectively.
Size exclusion chromatography with multiangle laser light scattering
Request a detailed protocolThe trimeric arrangement of the phosphorylated SMAD2 MH2 domain was confirmed by SEC-MALLS. In brief, an S200 10/300 Increase column was attached to an AKTA Micro FPLC system (GE Healthcare), which was connected to a Heleos Dawn 8+ followed by an Optilab TRex (Wyatt). For data collection, 100 µl of a 2 mg/ml stock of phosphorylated SMAD2 MH2 domain was injected and data collected at 0.5 ml/min for 60 min. The data were analyzed by ASTRA6.1 software.
Biolayer interferometry
Request a detailed protocolBiolayer interferometry was carried out using an Octet RED96 instrument (ForteBio). Biotinylated N-SKI peptide (residues 11–45) was immobilized on streptavidin-coated biosensors (ForteBio) at a concentration of 1 μg/ml in buffer containing 50 mM HEPES pH 7.5, 200 mM NaCl, 1 mg/ml bovine serum albumin, and 0.1% Tween-20 for 100 s. The immobilization typically reached a response level of 2 nm. Association and dissociation curves were obtained through addition of a dilution series of trimeric phosphorylated SMAD2 MH2 domain complex (15.6–1.95 μM) for 100 s followed by dissociation in buffer for 350 s using the Octet acquisition software. The binding data were fitted using the Octet analysis software.
Crystallization of the N-SKI–SMAD2 MH2 domain complex
Request a detailed protocolThe SKI peptide (amino acids 11–45) was synthesized by the Peptide Chemistry Group and added in a 2:1 ratio to the phosphorylated SMAD2 MH2 domain trimer. The complex was concentrated to 6 mg/ml and subject to crystallization trials. Initial screening gave rise to fine needles, in several conditions, and these were used as seedstock for rescreening. This gave rise to 25 µm crystals, with a cubic morphology, in the Crystal Screen Cryo screen (Hampton Research); the condition being 1.5 M ammonium sulphate, 0.15 M K Na tartrate, 0.08 M Na3citrate, and 25% (vol/vol) glycerol. The crystals were flash-frozen in liquid nitrogen, and data were collected on the I24 beamline at Diamond Light Source (DLS). Data was processed automatically using the DLS Xia2/XDS pipeline (Winter, 2010; Kabsch, 2010). The crystals belong to the I213 space group and diffracted to a resolution of 2.0 Å.
Structure determination
Request a detailed protocolMolecular replacement was undertaken with the CCP4 program Phaser (McCoy et al., 2007), utilising pdbfile 1khx (with the C-terminal tail removed), as the search model (Winn et al., 2011). Initial structure refinement was undertaken with Refmac (Murshudov et al., 2011), with manual model building in Coot (Emsley et al., 2010), before switching to Phenix.Refine to finalize the model (Liebschner et al., 2019; Afonine et al., 2012). Coordinates and data are available from the Protein Data Bank, with accession code: 6ZVQ.
RNA extraction, qRT-PCR, and RNA-sequencing
Request a detailed protocolTotal RNA was extracted using Trizol (Thermo Fisher Scientific) according to the manufacturer’s instructions. cDNA synthesis and qPCRs were performed as described (Grönroos et al., 2012). Primer sequences are listed in Key Resources Table. All qPCRs were performed with the PowerUp SYBR Green Master Mix (Thermo Fisher Scientific) with 300 nM of each primer and 2 µl of diluted cDNA. Fluorescence acquisition was performed on either a 7500 FAST machine or QuantStudio 12 Flex (Thermo Fisher Scientific). Calculations were performed using the ΔΔCt method, and levels of mRNA are expressed as fold change relative to untreated or SB-431542-treated control cells. Means ± SEM from at least three independent experiments are shown. Results were analyzed using GraphPad Prism 8 software and statistics were performed on these data using one-way or two-way analysis of variance (ANOVA) as stated in the figure legends.
For the RNA-sequencing experiment four biological replicates were used. Total RNA was extracted, and the quality of the RNA was assessed using a bioanalyzer (Agilent). Libraries were prepared using the KAPA mRNA HyperPrep kit (Roche), and single-end reads were generated using an Illumina HiSeq 4000.
RNA sequencing analysis
Request a detailed protocolRaw reads were quality and adapter trimmed using cutadapt-1.9.1 (Martin, 2011) prior to alignment. Reads were then aligned and quantified using RSEM-1.3.0/STAR-2.5.2 (Dobin et al., 2013; Li and Dewey, 2011) against the human genome GRCh38 and annotation release 89, both from Ensembl. TPM (transcripts per kilobase million) values were also generated using RSEM/STAR. Differential gene expression analysis was performed in R-3.6.1 (R Development Core Team, 2009) using the DESeq2 package (version 1.24.0) (Love et al., 2014). Differential genes were selected using a 0.05 false discovery rate threshold using the pairwise comparisons between time points within each cell line individually. Normalization and variance-stabilizing transformation was applied on raw counts before performing principal component analysis and Euclidean distance-based clustering.
For the heatmaps, we selected those genes that were significantly differentially expressed in the control cell line with an absolute log2 (fold change) of at least 0.5, not significantly differentially expressed in the SKI-mutated cell lines, and detected with at least a TPM value of 2 in the control cell lines in any condition. We then visualized the log2 (fold change) for each time point against the control condition. Heatmaps were generated in R-3.6.1 using the ComplexHeatmap (Gu et al., 2016) package (version 2.0.0). The raw data files have been uploaded to the European Genome-phenome Archive (EGA), accession number EGAS00001004908.
Quantifications
Request a detailed protocolQuantification of Western blots was performed by densitometry measurements of each lane using ImageJ software. The measurements were normalized to the loading control in the same blot. In each case, quantifications were normalized to the SB-431542-treated samples. Quantification for the flow cytometry was performed by measurement of the fluorescence emitted by Alexa 647 in EGFP-expressing single cells as a measure of SKIL protein, and levels were normalized to the SB-431542-treated sample in the SMAD4 knockout cells. The program FlowJo was used to analyze the results. For the peptide arrays, two independent experiments were quantified. The staining intensity of each array was quantified by systematically moving a circular selection (diameter 20 pixels) across the array and measuring the 8-bit greyscale intensity for each spot. Each intensity measure was normalized to the average intensity of 60 positive controls of the WT peptide after subtracting the background, measured from the average intensity of 60 negative controls (truncated SKI peptide C as indicated in Figure 4—figure supplement 1B).
Statistical analysis
Request a detailed protocolStatistical analysis was performed in Prism 8 (GraphPad). At least three independent experiments were performed for statistical analysis unless otherwise specified in the figure legends. Normalized values were log transformed for the statistical analysis. For comparison between more than two groups with one variable, one-way ANOVA was used followed by the Sidak’s correction test. For comparison between groups that have been split on two independent variables, two-way ANOVA was performed followed by Tukey’s multiple comparison tests.
Appendix 1
Data availability
Sequencing data have been uploaded to the European Genome-phenome Archive (EGA), accession number EGAS00001004908. Diffraction data have been deposited in PDB under the accession code 6ZVQ. All data generated or analysed during this study are included in the manuscript and supporting files. Source data files have been provided for Figures 1, 2, 4, 5, 6, 7, Figure 1 Supplement 1, Figure 2 Supplement 1, Figure 7 Supplement 2.
-
European Genome-Phenome ArchiveID EGAS00001004908. Mutations in SKI in Shprintzen-Goldberg syndrome lead to attenuated TGF-β responses through SKI stabilization.
-
RCSB Protein Data BankID 6ZVQ. Complex betweenSMAD2 MH2 domain and peptide from Ski corepressor.
References
-
Towards automated crystallographic structure refinement with phenix.refineActa Crystallographica Section D Biological Crystallography 68:352–367.https://doi.org/10.1107/S0907444912001308
-
Induction of heart valve lesions by small-molecule ALK5 inhibitorsToxicologic Pathology 39:916–924.https://doi.org/10.1177/0192623311416259
-
TGF-β (Transforming growth Factor-β) Signaling protects the thoracic and abdominal aorta from angiotensin II-Induced pathology by distinct mechanismsArteriosclerosis, Thrombosis, and Vascular Biology 37:2102–2113.https://doi.org/10.1161/ATVBAHA.117.309401
-
c-Ski in health and diseaseCell and Tissue Research 347:51–64.https://doi.org/10.1007/s00441-011-1180-z
-
TGF-β signalopathies as a paradigm for translational medicineEuropean Journal of Medical Genetics 58:695–703.https://doi.org/10.1016/j.ejmg.2015.10.010
-
TGFBR1 mutations associated with Loeys-Dietz syndrome are inactivatingJournal of Receptors and Signal Transduction 32:150–155.https://doi.org/10.3109/10799893.2012.664553
-
In-frame mutations in exon 1 of SKI cause dominant Shprintzen-Goldberg syndromeThe American Journal of Human Genetics 91:950–957.https://doi.org/10.1016/j.ajhg.2012.10.002
-
Dimorphic effects of transforming growth factor-β signaling during aortic aneurysm progression in mice suggest a combinatorial therapy for marfan syndromeArteriosclerosis, Thrombosis, and Vascular Biology 35:911–917.https://doi.org/10.1161/ATVBAHA.114.305150
-
STAR: ultrafast universal RNA-seq alignerBioinformatics 29:15–21.https://doi.org/10.1093/bioinformatics/bts635
-
Features and development of cootActa Crystallographica. Section D, Biological Crystallography 66:486–501.https://doi.org/10.1107/S0907444910007493
-
Protein complex expression by using multigene baculoviral vectorsNature Methods 3:1021–1032.https://doi.org/10.1038/nmeth983
-
Angiotensin II-dependent TGF-β signaling contributes to Loeys-Dietz syndrome vascular pathogenesisJournal of Clinical Investigation 124:448–460.https://doi.org/10.1172/JCI69666
-
Shprintzen-Goldberg syndrome: a clinical analysisAmerican Journal of Medical Genetics 76:202–212.https://doi.org/10.1002/(SICI)1096-8628(19980319)76:3<202::AID-AJMG2>3.0.CO;2-S
-
Computational investigation of a ternary model of SnoN-SMAD3-SMAD4 complexComputational Biology and Chemistry 83:107159.https://doi.org/10.1016/j.compbiolchem.2019.107159
-
XDSActa Crystallographica. Section D, Biological Crystallography 66:125–132.https://doi.org/10.1107/S0907444909047337
-
Arkadia activates Smad3/Smad4-dependent transcription by triggering signal-induced SnoN degradationMolecular and Cellular Biology 27:6068–6083.https://doi.org/10.1128/MCB.00664-07
-
Macromolecular structure determination using X-rays, neutrons and electrons: recent developments in phenixActa Crystallographica Section D Structural Biology 75:861–877.https://doi.org/10.1107/S2059798319011471
-
The Ski oncoprotein interacts with the Smad proteins to repress TGFβ signalingGenes & Development 13:2196–2206.https://doi.org/10.1101/gad.13.17.2196
-
Ski and SnoN: negative regulators of TGF-β signalingCurrent Opinion in Genetics & Development 14:65–70.https://doi.org/10.1016/j.gde.2003.11.003
-
Lineage-specific events underlie aortic root aneurysm pathogenesis in Loeys-Dietz syndromeJournal of Clinical Investigation 129:659–675.https://doi.org/10.1172/JCI123547
-
TGFβ signalling in contextNature Reviews Molecular Cell Biology 13:616–630.https://doi.org/10.1038/nrm3434
-
Phaser crystallographic softwareJournal of Applied Crystallography 40:658–674.https://doi.org/10.1107/S0021889807021206
-
BookTGF-β superfamily signallingIn: Bradshaw R. A, Stahl P. D, editors. Encyclopedia of Cell Biology. Elsevier. pp. 1–308.
-
REFMAC5 for the refinement of macromolecular crystal structuresActa Crystallographica Section D Biological Crystallography 67:355–367.https://doi.org/10.1107/S0907444911001314
-
Arkadia induces degradation of SnoN and c-Ski to enhance transforming growth factor-beta signalingJournal of Biological Chemistry 282:20492–20501.https://doi.org/10.1074/jbc.M701294200
-
Analysis of smad nucleocytoplasmic shuttling in living cellsJournal of Cell Science 117:4113–4125.https://doi.org/10.1242/jcs.01289
-
Functional characterization of transforming growth factor beta signaling in Smad2- and Smad3-deficient fibroblastsJournal of Biological Chemistry 276:19945–19953.https://doi.org/10.1074/jbc.M102382200
-
Transforming growth factor beta-independent shuttling of Smad4 between the cytoplasm and nucleusMolecular and Cellular Biology 20:9041–9054.https://doi.org/10.1128/MCB.20.23.9041-9054.2000
-
SoftwareR: A language and environment for statistical computingR Foundation for Statistical Computing, Vienna, Austria.
-
Genome engineering using the CRISPR-Cas9 systemNature Protocols 8:2281–2308.https://doi.org/10.1038/nprot.2013.143
-
Recognition of phosphorylated-Smad2-containing complexes by a novel smad interaction motifMolecular and Cellular Biology 24:1106–1121.https://doi.org/10.1128/MCB.24.3.1106-1121.2004
-
LTBPs in biology and medicine: ltbp diseasesMatrix Biology 71-72:90–99.https://doi.org/10.1016/j.matbio.2017.11.014
-
Regulation of the bioavailability of TGF-β and TGF-β-Related proteinsCold Spring Harbor Perspectives in Biology 8:a021907.https://doi.org/10.1101/cshperspect.a021907
-
The SMAD-binding domain of SKI: a hotspot for de novo mutations causing Shprintzen-Goldberg syndromeEuropean Journal of Human Genetics 23:224–228.https://doi.org/10.1038/ejhg.2014.61
-
A mutation update on the LDS-associated genes TGFB2/3 and SMAD2/3Human Mutation 39:621–634.https://doi.org/10.1002/humu.23407
-
A recurrent pattern syndrome of craniosynostosis associated with arachnodactyly and abdominal herniasJournal of Craniofacial Genetics and Developmental Biology 2:65–74.
-
Viral ski inhibits retinoblastoma protein (Rb)-mediated transcriptional repression in a dominant negative fashionJournal of Biological Chemistry 274:4485–4488.https://doi.org/10.1074/jbc.274.8.4485
-
Phenotypic spectrum of the SMAD3-related aneurysms-osteoarthritis syndromeJournal of Medical Genetics 49:47–57.https://doi.org/10.1136/jmedgenet-2011-100382
-
Marfan syndrome and related disorders: 25 years of gene discoveryHuman Mutation 37:524–531.https://doi.org/10.1002/humu.22977
-
SnoN stabilizes the SMAD3/SMAD4 protein complexScientific Reports 7:46370.https://doi.org/10.1038/srep46370
-
Aortopathy in a mouse model of marfan syndrome is not mediated by altered transforming growth factor β signalingJournal of the American Heart Association 6:e004968.https://doi.org/10.1161/JAHA.116.004968
-
Early surgical experience with Loeys-Dietz: a new syndrome of aggressive thoracic aortic aneurysm diseaseThe Annals of Thoracic Surgery 83:S757–S763.https://doi.org/10.1016/j.athoracsur.2006.10.091
-
Overview of the CCP4 suite and current developmentsActa Crystallographica. Section D, Biological Crystallography 67:235–242.https://doi.org/10.1107/S0907444910045749
-
xia2 : an expert system for macromolecular crystallography data reductionJournal of Applied Crystallography 43:186–190.https://doi.org/10.1107/S0021889809045701
-
Smad3-Smad4 and AP-1 complexes synergize in transcriptional activation of the c-Jun promoter by transforming growth factor betaMolecular and Cellular Biology 19:1821–1830.https://doi.org/10.1128/MCB.19.3.1821
Article and author information
Author details
Funding
Francis Crick Institute (FC10095)
- Ilaria Gori
- Roger George
- Andrew G Purkiss
- Stephanie Strohbuecker
- Rebecca A Randall
- Roksana Ogrodowicz
- Dhira Joshi
- Svend Kjaer
- Caroline S Hill
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We would like to thank David Abraham for the normal dermal fibroblasts. We are very grateful to the Francis Crick Institute Advanced Sequencing Facility, the Light Microscopy Facility, the Flow Cytometry Facility, Cell Services, the Genomics Equipment Park and to the Peptide Chemistry Facility. We thank Andrew Economou for quantification of the peptide arrays. We thank all the members of the Hill lab and Davide Coda, Daniel Miller, and Anassuya Ramachandran for helpful discussions and very useful comments on the manuscript. This work was supported by the Francis Crick Institute which receives its core funding from Cancer Research UK (FC001095), the UK Medical Research Council (FC001095), and the Wellcome Trust (FC001095).
Ethics
Human subjects: Dermal fibroblasts from healthy subjects were kindly provided by David Abraham (UCL-Medical School Royal Free Campus) under the ethics of the Health Research Authority, NRES Committee London - Hampstead, Research Ethics Committee (REC) reference, 6398. L32V and ΛS94-97 SKI dermal fibroblasts were obtained from Laurence Faivre and Virginie Carmignac (Université de Bourgogne UMR1231 GAD, Dijon, France) under the ethics of the GAD collection, number DC2011-1332 (Carmignac et al., 2012).
Copyright
© 2021, Gori et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,276
- views
-
- 255
- downloads
-
- 22
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
Pre-mRNA splicing is catalyzed in two steps: 5ʹ splice site (SS) cleavage and exon ligation. A number of proteins transiently associate with spliceosomes to specifically impact these steps (first and second step factors). We recently identified Fyv6 (FAM192A in humans) as a second step factor in Saccharomyces cerevisiae; however, we did not determine how widespread Fyv6’s impact is on the transcriptome. To answer this question, we have used RNA sequencing (RNA-seq) to analyze changes in splicing. These results show that loss of Fyv6 results in activation of non-consensus, branch point (BP) proximal 3ʹ SS transcriptome-wide. To identify the molecular basis of these observations, we determined a high-resolution cryo-electron microscopy (cryo-EM) structure of a yeast product complex spliceosome containing Fyv6 at 2.3 Å. The structure reveals that Fyv6 is the only second step factor that contacts the Prp22 ATPase and that Fyv6 binding is mutually exclusive with that of the first step factor Yju2. We then use this structure to dissect Fyv6 functional domains and interpret results of a genetic screen for fyv6Δ suppressor mutations. The combined transcriptomic, structural, and genetic studies allow us to propose a model in which Yju2/Fyv6 exchange facilitates exon ligation and Fyv6 promotes usage of consensus, BP distal 3ʹ SS.
-
- Biochemistry and Chemical Biology
- Neuroscience
TIMM50, an essential TIM23 complex subunit, is suggested to facilitate the import of ~60% of the mitochondrial proteome. In this study, we characterized a TIMM50 disease-causing mutation in human fibroblasts and noted significant decreases in TIM23 core protein levels (TIMM50, TIMM17A/B, and TIMM23). Strikingly, TIMM50 deficiency had no impact on the steady-state levels of most of its putative substrates, suggesting that even low levels of a functional TIM23 complex are sufficient to maintain the majority of TIM23 complex-dependent mitochondrial proteome. As TIMM50 mutations have been linked to severe neurological phenotypes, we aimed to characterize TIMM50 defects in manipulated mammalian neurons. TIMM50 knockdown in mouse neurons had a minor effect on the steady state level of most of the mitochondrial proteome, supporting the results observed in patient fibroblasts. Amongst the few affected TIM23 substrates, a decrease in the steady state level of components of the intricate oxidative phosphorylation and mitochondrial ribosome complexes was evident. This led to declined respiration rates in fibroblasts and neurons, reduced cellular ATP levels, and defective mitochondrial trafficking in neuronal processes, possibly contributing to the developmental defects observed in patients with TIMM50 disease. Finally, increased electrical activity was observed in TIMM50 deficient mice neuronal cells, which correlated with reduced levels of KCNJ10 and KCNA2 plasma membrane potassium channels, likely underlying the patients’ epileptic phenotype.