Selective Disruption of Synaptic BMP Signaling by a Smad Mutation Adjacent to the Highly Conserved H2 Helix.

Bone morphogenetic proteins (BMPs) shape normal development and function via canonical and noncanonical signaling pathways. BMPs initiate canonical signaling by binding to transmembrane receptors that phosphorylate Smad proteins and induce their translocation into the nucleus and regulation of target genes. Phosphorylated Smads also accumulate at cellular junctions, but this noncanonical, local BMP signaling modality remains less defined. We have recently reported that phosphorylated Smad (pMad in Drosophila) accumulates at synaptic junctions in protein complexes with genetically distinct composition and regulation. Here, we examined a wide collection of Drosophila Mad alleles and searched for molecular features relevant to pMad accumulation at synaptic junctions. We found that strong Mad alleles generally disrupt both synaptic and nuclear pMad, whereas moderate Mad alleles have a wider range of phenotypes and can selectively impact different BMP signaling pathways. Interestingly, regulatory Mad mutations reveal that synaptic pMad appears to be more sensitive to a net reduction in Mad levels than nuclear pMad. Importantly, a previously uncharacterized allele, Mad 8 , showed markedly reduced synaptic pMad but only moderately diminished nuclear pMad. The postsynaptic composition and electrophysiological properties of Mad 8 neuromuscular junctions (NMJs) were also altered. Using biochemical approaches, we examined how a single point mutation in Mad 8 could influence the Mad-receptor interface and identified a key motif, the H2 helix. Our study highlights the biological relevance of Smad-dependent, synaptic BMP signaling and uncovers a highly conserved structural feature of Smads, critical for normal development and function.

BONE morphogenetic proteins (BMPs) modulate a wide variety of cellular processes via canonical and noncanonical signaling pathways (Massague 1990; Hogan 1996; Derynck and Zhang 2003). Misregulation of BMP signaling is associated with developmental abnormalities and disease states, highlighting the need for tight regulation of various BMP pathways. Like all members of the TGF-β superfamily, BMPs form biologically active dimers that initiate signaling by binding to a hetero-tetrameric complex of Ser/Thr kinases, known as Type I and Type II BMP receptors (BMPRs). Type II receptors are constitutive kinases that phosphorylate Type I receptors within a regulatory glycine-serine-rich (GS) domain to activate them. Activated Type I receptors bind to and phosphorylate the intracellular transducers of the BMP pathway, the R-Smads (Smad-1, -5 or -8 in vertebrates and Mad in Drosophila) (Feng and Derynck 2005; Schmierer and Hill 2007). Phosphorylated R-Smads (pSmads) have a propensity to trimerize that favors their dissociation from the receptors (Kawabata ). In the canonical pathway, pSmads associate with co-Smads, translocate into the nucleus, and, in conjunction with other transcription factors, modulate expression of target genes. Activated BMPRs can also signal independently of Smads through noncanonical pathways that include mitogen-activated protein kinase (MAPK), LIM kinase, phosphatidylinositol 3-kinase/Akt (PI3K/Akt), and Rho-like small GTPases (Derynck and Zhang 2003; Moustakas and Heldin 2005; Zhang 2009). More recently, pSmad accumulation at cell membranes has been reported in at least two instances: (i) at tight junctions during neural tube closure (Eom ), and (ii) at the Drosophila neuromuscular junction (NMJ) (Dudu ; Smith ). During neural tube closure, pSmad1/5/8 binds to apical polarity complexes and mediates stabilization of BMP/BMPR complexes at tight junctions (Eom ); prolonged BMP blockade disrupts tight junctions and affects epithelial organization (Eom ). In both cases, these Smad-dependent functions do not require transcriptional regulation. In epithelial cells, nuclear pSmad signals tend to be very strong, obscuring junctional pSmad signals. In contrast, at the fly NMJ, junctional pMad localizes at synaptic terminals [the contacts between motor neurons (MNs) and body-wall muscles], whereas nuclear pMad accumulates in MN nuclei. The spatial separation between nuclear and synaptic pMad was instrumental in the initial characterization of the Smad-dependent local BMP signaling pathway (Sulkowski , 2016).

Flies rely on BMP signaling for NMJ growth and neurotransmitter release (Marqués and Zhang 2006). Many of these functions are triggered by Gbb—a BMP7 homolog that binds to presynaptic BMPRII, Wishful thinking (Wit), and the BMPRIs, Thickveins (Tkv) and Saxophone (Sax). During canonical signaling, this high-order BMP/BMPR complex is endocytosed and retrogradely transported to the MN soma, where Mad is phosphorylated and regulates transcriptional programs with distinct roles in the development of the NMJ. Gbb and Wit also signal noncanonically through the effector protein LIM kinase 1 (LIMK1) to regulate synapse stability (Eaton and Davis 2005). A third BMP signaling modality, Smad-dependent synaptic BMP signaling, does not require Gbb but does require Wit, Tkv, Sax, and the postsynaptic type-A glutamate receptors (DiAntonio 2006; Sulkowski ). We have previously demonstrated that synaptic pMad is involved in a positive feedback loop across the synaptic cleft. Active, postsynaptic type-A (GluRIIA-containing) glutamate receptors trigger presynaptic accumulation of pMad, which, in turn, functions to stabilize type-A receptors at postsynaptic densities (Sulkowski , 2016). Genetic and cell biology studies suggest that synaptic Mad is phosphorylated locally by activated BMPRs confined to the active zone, a region with massive (synaptic vesicles) exocytosis but no endocytosis. This region effectively traps BMPRs, precluding them from endocytosis and participation in canonical BMP signaling. But how does pMad, a product of an enzymatic reaction, remain associated with its own kinase?

To learn more about the structural features of Mad that may influence its association with the BMPRs, we have collected most of the existing Drosophila alleles and compared them side-by side for the ability to sustain the two Smad-dependent signaling modalities: canonical BMP signaling, marked by pMad accumulation in MN nuclei, and synaptic BMP signaling, marked by pMad accumulation at synaptic terminals. Mad is a modular, highly conserved protein that contains an N-terminal MH1 (Mad homology 1) DNA binding domain, a C-terminal MH2 protein interaction domain, and a linker that has been implicated in the crosstalk with other signaling pathways (Hoodless ). Changes in the structure and the oligomeric state of the MH2 domain contribute in two ways to the directionality of the signaling process (Wu ). First, the MH2 mediates Mad binding to BMPRI and thus phosphorylation of its C-terminal SSXS motif, the site of BMP-dependent phosphorylation (Hoodless ; Macías-Silva ; Zhang ). Second, the MH2 is critical to formation of Smad trimers that dissociate from BMPRs (Kawabata ). In particular, the L3 loop has been implicated in mutually exclusive interactions with the BMPRI and the phosphorylated C-terminal SSXS motif (Wu ). Most molecular lesions in our collection map to the MH2 domain, including in the L3 loop, making this allelic series particularly suitable for studies on Mad–Tkv interactions.

Within this comprehensive collection, we found that strong alleles generally disrupt both synaptic and nuclear pMad accumulation, whereas moderate alleles have a broader range of phenotypes and can differentially impact different BMP signaling modalities. In particular, Mad showed drastically reduced synaptic pMad levels, but only moderately diminished nuclear pMad. The postsynaptic composition and electrophysiological properties of Mad NMJs were likewise altered. Using biochemical assays and structural modeling, we have examined how point mutations such as S359L, present in Mad, could influence the Mad–Tkv interface. Our study identified a new molecular determinant for this Mad–Tkv interaction, the well-conserved H2 helix. Several genetic variants identified in human patients map to H2, underscoring the relevance of this motif for normal development and function.

Materials and Methods

Fly stocks

Drosophila stocks used in this study are as follows: w, Mad, Mad, Mad, Mad, Mad, Mad, Mad, Mad, Mad, Mad, Mad, Mad (Sekelsky ), Mad (here Mad) (Dworkin and Gibson 2006), Mad (Bellen ), Mad (Wiersdorff ), Df(2L)C28 (Raftery ), and twitMI06552 (Venken ). The flies were reared at 25° on Jazz-Mix food (Fisher Scientific). To control for larval crowding, eight to ten females were crossed with five to seven males per vial and passed to fresh vials every 3 days.

Molecular constructs

Flag-Mad plasmid was previously described (Shimmi and O’Connor 2003). Flag-Mad8 was generated by PCR followed by Gibson assembly (NEB) and the S359L substitution was verified by sequencing. The following PCR primers were utilized:

Mad8-For: 5′-CCAGTACTCGTTCCTCGCCACTCGGAATTCGCGCCCGGTCACTCG;

Mad8-Rev: 5′-GCGGTATTTTGCACACGGTTAGCGGATGGAATCCGTGGTGG.

To generate the Tac–Tkv construct, Tac extracellular and transmembrane sequences (HindIII–EcoRI fragment) (Ren ) were joined with the PCR amplified wild-type or activated Tkv cytoplasmic domain in the pIB/V5-His vector (ThermoFisher). The activated Tkv chimera (Tac–TkvA) had a single residue substitution Q199D at the end of the GS box (Wieser ).

The following Tkv primers were utilized:

Tkv-For: 5′-GCGTCCTCCTCCTGAGTGGGCTCTGTTTCACCTACAAGCGACGCGAGAAGC;

Tkv-Rev: 5′-GGCTTACCTTCGAACCGCGGGCCCTCTAGACAATCTTAATGGGCACATCG.

Immunohistochemistry

Larvae were dissected as described previously in ice-cooled Ca2+-free HL-3 solution (Stewart ; Budnik ). The samples were fixed in either 4% formaldehyde (Polysciences) for 30 min or in Bouin’s fixative (Bio-Rad) for 3 min and washed in phosphate-buffered saline (PBS) containing 0.5% Triton X-100. Primary antibodies from the Developmental Studies Hybridoma Bank were used at the following dilutions: mouse anti-GluRIIA (MH2B), 1:200; rat anti-Elav (7E8A10), 1:200. Other primary antibodies were as follows: rabbit anti-phosphorylated Mad (pMad), 1:500, (a gift from C. H. Heldin), rabbit anti-GluRIIC, 1:1000 (Ramos ), Cy5-conjugated goat anti-horseradish peroxidase (HRP), 1:1000 (Jackson ImmunoResearch Laboratories, Inc.). Alexa Fluor 488-, Alexa Fluor 568-, and Alexa Fluor 647- conjugated secondary antibodies (Molecular Probes) were used at 1:400. Larval filets and brains were mounted in ProLong Gold (Invitrogen). Samples of different genotypes were processed simultaneously and imaged under identical confocal settings using laser scanning confocal microscopes (CarlZeiss LSM780).

Fluorescence intensity measurements

Maximum intensity projections were used for quantitative image analysis in Fiji. To measure pMad at the synaptic terminal, synapse surface area was calculated by creating a mask around the HRP signal that labels the neuronal membrane. For each sample, a threshold was applied manually to the HRP and pMad channels to remove irrelevant low intensity pixels outside the NMJ area. The mean puncta intensity was calculated as the total fluorescence intensity signal of the puncta divided by their area. All intensities were normalized to control values within an experimental set. The same method was used to quantify synaptic GluRIIA and GluRIIC levels. For nuclear pMad, MNs nuclei were labeled with anti-Elav, then thresholds were applied to the Elav and pMad signals as above. Mean puncta intensity was calculated as the total fluorescence intensity signal of the puncta divided by the area of the puncta. Bouton area was calculated by separating the HRP-positive regions of interest surrounding synaptic pMad immunoreactivities. The GluRIIC-positive puncta were counted with Imaris 9.5.1 (Bitplane) by first converting the confocal stacks into 3D reconstructed NMJs, then applying built-in spot detection algorithms. Briefly, the point of maximum intensity in each spot was defined as the center of GluRIIC puncta. The intensity threshold was adjusted to detect all GluRIIC-positive spots. False positive and false negative spots were inspected manually and corrected via the edit and select functions.

Student’s t-test was performed using R software (version 3.6.2) to evaluate statistical significance. All graphs represent mean value of all animals of indicated genotypes ± SEM.

Cell-based assay

A cell-based assay for BMP signaling was described previously (Shimmi and O’Connor 2003; Serpe and O’Connor 2006). In brief, S2 cells were transfected with Flag-Mad or Flag-Mad8. At 3 days after transfection, cells were incubated with 1nM Dpp (R&D Systems) for 1 hr, and cell extracts were analyzed by Western blotting. The pMad levels were revealed with anti-phospho-Smad 1/5 (Ser463/465) 41D10 (Cell Signaling) 1:200 and anti-Tubulin (Invitrogen) 1:500, and quantified relative to the Flag signal detected with anti-Flag M2 (Sigma) at 1:2000 using IRDye secondary antibodies for simultaneous detection on an Odyssey Infrared Imaging System (Li-Cor Biosciences).

For immunocytochemistry experiments, eight-well chambers (Fisher Scientific) were coated with 0.1 mg/ml Concanavalin A (Sigma) for 1 hr at 25°, then with 1 ng/ml anti-Tac antibodies M-A251, (BioLegend) for 1 hr at 25°. Transiently transfected S2 cells were grown for 3 days before spreading on ConA/anti-Tac- coated chamber for 1 hr at 25°. The surface-attached cells were fixed in 4% formaldehyde (Polysciences) for 15 min, then stained for pMad (anti-phospho-Smad 1/5, 41D10, 1:200, Cell Signaling), and Flag (anti-Flag, M2, 1:500, Sigma). FITC-conjugated goat anti-phalloidin, 1:500 (Invitrogen), Alexa Fluor 568-, and Alexa Fluor 647- conjugated secondary antibodies (Molecular Probes) were used at 1:400. Cells were mounted in DAPI-containing Prolong Gold (Invitrogen).

Electrophysiology

Recordings were performed on muscle 6, segment A3 of third instar larvae as previously reported (Qin ). Briefly, wandering third-instar larvae were dissected in ice-cold, calcium-free physiological HL-3 saline (Stewart ) and immersed in HL-3 containing Ca2+ before being shifted to the recording chamber. The calcium-free HL-3 saline contains (in mM): 70 NaCl, 5 KCl, 20 MgCl2, 10 HCO3, 5 trehalose, 115 sucrose, 5 HEPES, pH adjusted to 7.2 at room temperature. The recording solution was HL-3 with 0.5 mM CaCl2. Intracellular electrodes (borosilicate glass capillaries of 1 mm diameter) were filled with 3 M KCl and resistances ranged from 12 to 25 MΩ. Recordings were done at room temperature from muscle cells with an initial membrane potential between −50 and −70 mV, and input resistances of ≥ 4 MΩ. For miniature excitatory junction potentials (mEJC) recordings, the muscle cells were clamped to −80 mV. To calculate mean amplitudes and frequency of mEJCs, 100–150 events from each muscle were measured and averaged using the Mini Analysis program (Synaptosoft). Minis with a slow rise and falling time arising from neighboring electrically coupled muscle cells were excluded from analysis (Gho 1994; Zhang ). To compare decay time constant of mEJCs between genotypes, 50 clear representative events from each recording were averaged and fitted using a double exponential equation of the formwhere I is the peak current amplitude of a decay component and τx is the corresponding decay time constant. To allow for easier comparison of decay times between genotypes weighted τ (ms) were calculated using the formulaElectrical signals were recorded with an Axoclamp 2B amplifier (Axon Instruments). The signals were filtered at 1 kHz and digitized at 10 kHz by using an analog-digital converter (Digidata 1440A) and pCLAMP software (version 10.0, Axon Instruments). Data are presented as mean ± SEM. Two-tailed unpaired Student’s t-test was used to assess statistically significant differences among genotypes. Differences were considered significant at P < 0.05.

Data availability

All reagents are available on request. The authors affirm that all data necessary for confirming the conclusions of the article are present within the article, main figures and tables.

Results

Analyses of pMad in larval MNs in distinct Mad alleles

To test the effect of mutations on the distribution and levels of synaptic and nuclear pMad, 15 available alleles were examined (Figure 1, A and B). Most of these alleles were isolated as dominant maternal enhancers of recessive mutations (Raftery ; Sekelsky ; Chen ). We crossed each of these alleles with Df(2L)C28, a small deficiency covering the locus, counted the resulted progenies, and calculated the percent of expected adult progeny that were trans-heterozygous for the allele and the deletion (Figure 1C). As expected, we observed only a few Mad/Df(2L)C28 adult escapers, with the largest escaper percentage from the weak Mad and Mad alleles, previously classified as hemizygous viable (Sekelsky ). In all crosses, we observed larval and pupal lethality as -deficient animals were unable to enter and/or complete metamorphosis. The Mad/Df(2L)C28 larvae, henceforth referred to as “Mad mutants,” have developmental delays and appear transparent due to reduced fat body. Since BMP signaling, primarily through Gbb/BMP7, is a central player in the energy homeostasis (Ballard ), we used fat body accumulation as an additional metric for the severity of mutant phenotypes (below).

Figure 1

alleles analyzed. (A) Diagram of the locus and transposable element alleles. (B) Mad domain organization and point mutations. (C) List of alleles analyzed with viability data.

To quantify the accumulation of pMad in MN nuclei, we costained third-instar larval ventral ganglia for pMad and Elav, a neuronal nuclear factor. For synaptic pMad, we quantified the pMad NMJ signals relative to anti-HRP, which labels neuronal membranes (Jan and Jan 1982). We used as a negative control the well-characterized Mad mutant to examine any residual pMad staining at the larval body wall muscle NMJ and in MNs (Figure 2, A and B). The protein encoded by Mad allele has a C-terminal truncation (Q417* stop) preceding the BMP-dependent phosphorylation site, S451SXS, that abolishes any zygotic pMad signals (Sekelsky ). Indeed, in Mad mutants, synaptic, and nuclear pMad were reduced by 80 ± 5% and 91 ± 5%, respectively (Table 1). The residual signals may be due to (i) maternal wild-type Mad protein that may persist at very low levels during third-instar larval stages; (ii) cross reaction of the pMad antibodies with pSmox, the R-Smad of the Activin signaling pathway in Drosophila, that shares high homology with Mad over the C-terminus, or (iii) background staining.

Figure 2

Strong alleles disrupt both synaptic and nuclear pMad. (A–D) Representative confocal images of NMJ 6/7 boutons from third-instar larvae of indicated genotypes labeled for pMad (red) and HRP (green). Synaptic pMad is shown in Fire-lut representation in the upper panels; on the intensity scale, white represents peak intensity (∼6000 arbitrary units, A.U.). (E–H) Confocal images of ventral nerve cord of larvae of indicated genotypes labeled with pMad (red in the merged panels) and Elav (green). In comparison with control animals (E), nuclear pMad is strongly reduced in Mad (F), Mad (G), and Mad (H) mutants. (I–L) Third instar control larvae (I) have an opaque appearance, while mutants (J–L) are almost transparent due to severely reduced fat body. (M and N) Quantification of synaptic (M) and nuclear pMad (N) levels. See also Table 1. Bar, 10 μm (A–H) and 1 μm (details); 100 μm (I–L). Error bars indicate SEM. ***P < 0.0001.

Table 1

Summary of synaptic and nuclear pMad levels with the corresponding bouton area

Figure	Genotype	Synaptic pMad				Nuclear pMad				Bouton area
		Mean	SEM	P	N	Mean	SEM	P	N	Mean	SEM	P	N
2	Control	100	2.47	1	30	100	2.08	1	29	100	2.08	1	16
2	Mad12/Df	19.70	5.25	4.243E-11	18	9.36	5.19	<2.2E-16	17	22.37	3.73	1.82E-07	18
2	Mad7/Df	23.55	4.44	1.698E-10	21	13.26	4.85	<2.2E-16	14	15.23	3.25	3.148E-08	14
2	Mad10/Df	22.88	4.51	1.331E-11	23	20.49	3.60	<2.2E-16	15	19.54	1.05	4.878E-09	14
2	Mad4/Df	24.07	1.51	<2.2E-16	17	7.82	2.33	<2.2E-16	12	8.63	2.17	1.595E-08	12
Additional data	Mad5/Df	99.05	5.38	0.83	14	97.56	4.64	0.6252	17	97.67	4.71	0.8487	17
Additional data	Mad6/Df	88.27	2.18	0.87	14	93.67	4.43	0.7717	13	84.73	5.61	0.1168	13
3	Mad237/Df	63.99	3.01	1.585E-11	41	75.31	3.36	1.31E-06	19	59.30	3.65	4.766E-06	17
Additional data	Mad237/Mad237	58.58	3.02	1.47E-13	36	60.03	4.20	9.42E-09	17
Additional data	Mad8−2/Df	50.23	4.40	2.42E-11	21	57.35	4.11	3.1E-09	16	64.05	4.64	0.007	16
3	Mad9/Df	62.92	4.43	1.516E-07	17	79.94	2.01	1.92E-05	16	65.34	4.38	0.0004	15
3	Mad11/Df	61.75	1.48	1.398E-11	19	74.02	4.97	0.000199	13	68.79	4.33	0.001	15
4	Mad1/Df	45.98	2.86	5.63E-10	19	31.62	5.06	1.02E-08	15	36.76	2.32	8.854E-08	18
	Mad2/Df	N/A				N/A				N/A
	Mad3/Df	N/A				N/A				N/A
	MadKG00581/Df	N/A				N/A				N/A
5	Mad8/Df	53.52	1.67	3.44E-09	22	75.42	3.26	1.42E-12	15	57.18	3.31	3.00E-07	26

Many alleles have diminished HRP-labeled NMJ terminals and a significantly decreased bouton area, prompting us to find additional metrics to assess differences in staining at synaptic terminals. For example, in maximum intensity projections, the HRP-positive NMJ area in Mad mutants was reduced to 41.68 ± 6.15% (n = 18, P < 0.0001), compared to control (n = 16), whereas the bouton area, measured as HRP- and pMad-positive, was reduced to 22.37 ± 3.73% (Table 1 and Materials and Methods). For this study, we have selected the bouton area as the metric for NMJ growth. This parameter consolidates the number and size of synaptic boutons, and appears less prone to technical variability than the HRP-positive NMJ area. For consistency, we performed all NMJ analyses at NMJ 6/7 in the A3 abdominal segment.

Strong Mad alleles drastically reduce both synaptic and nuclear pMad signals

Mad and Mad are strong alleles with single residue substitutions within the MH2 domain: V419M in Mad and G409S in Mad, respectively (Sekelsky ; Takaesu ). Both residues map to the L3 loop of the MH2 domain, which has been implicated in mutually exclusive interactions with either BMPRIs or with the phosphorylated pS-X-pS tail of adjacent Smads in trimeric complexes (Wu ). As expected, both Mad and Mad exhibited a drastic reduction of synaptic and nuclear pMad (Figure 2, C and D, quantified in M,N). Compared to the w control, synaptic pMad levels were reduced to 24 ± 4% in Mad and to 23 ± 4% in Mad. Such levels were very close to those observed in the negative control, Mad (Figure 2B and Table 1). Nuclear pMad signals were similarly reduced to 13 ± 5% and 20 ± 5% in Mad and Mad, respectively. The Mad, Mad, and Mad mutant larvae appear translucent, with strongly reduced fat body, indicating severely impaired canonical BMP signaling (Figure 2, I–L).

Mad, Mad, and Mad were previously classified as null alleles based on their maternal effect enhancement of (Sekelsky ). Our analyses indicate that these alleles also behave as very strong hypomorphs during larval stages. In contrast, Mad was initially classified as a moderate allele; however, we found that both synaptic and nuclear pMad levels were drastically reduced in the Mad mutant larvae to 24 ± 2% and 8 ± 2% respectively, compared to control (Table 1). This result could not be explained by additional lesions on the Mad chromosome, because we analyzed Mad/Df(2L)C28 trans-heterozygous animals. Rather, we speculate that the unidentified molecular lesion in the Mad allele affects critical function(s) during larval development.

In contrast, the weak Mad and Mad mutants showed no significant changes in either synaptic or nuclear pMad levels during larval stages (Table 1) and had normal larval fat body (not shown). These alleles were isolated independently but have the same single residue change, R272H (Sekelsky ; Takaesu ). Nonetheless, these mutants yielded few adult escapers (Figure 1B), indicating that the substitution of this residue triggers critical deficits during later developmental stages.

Moderate Mad alleles exhibit differential effects on synaptic and nuclear pMad

While single point mutations disrupt functional domains within the Mad protein, the regulatory Mad allele (or Mad) has reduced overall Mad expression due to the insertion of a transposable element within the 5′UTR of the Mad mRNA (Dworkin and Gibson 2006). Previous studies reported significant NMJ defects in Mad homozygous and trans-heterozygous combinations (Merino ). We also found disruptions of pMad levels in Mad third-instar larvae (Figure 3, A and B and Table 1). Interestingly, the pMad levels were differentially disrupted in Mad mutants, with the synaptic pMad reduced to 64 ± 3% and the nuclear pMad to 75 ± 3% compared to control. A separate regulatory allele, Mad, with a lesion in the first intron (Wiersdorff ), showed a very similar phenotype (Table 1). Thus, synaptic pMad appears to be more sensitive to suboptimal Mad levels than nuclear pMad.

Figure 3

Moderate alleles exhibit differential effects on synaptic and nuclear pMad. (A–H) Confocal images of NMJ 6/7 boutons and ventral nerve cord from third-instar larvae of indicated genotypes labeled for pMad (red) and HRP (green) (A–D) or Elav (green) (E–H). Synaptic pMad is showed in Fire-lut in the upper panels; on the intensity scale, white represents peak intensity (∼6000 arbitrary units, A.U.). (E–H) Confocal images of ventral nerve cord of larvae of indicated genotypes labeled with pMad (red in the merged panels) and Elav (green). (I–L) Images of third instar larvae. mutants (J–L) appear slightly translucent in comparison with control (I) due to mildly reduced fat body. (M and N) Quantification of synaptic (M) and nuclear pMad (N) levels. See also Table 1. Bar, 10 μm (A–H) and 1 μm (details); 100 μm (I–L). Error bars indicate SEM. ***P < 0.0001.

A similar trend was observed for two other moderate alleles, Mad, and Mad, which code for single residue substitutions within the L3 loop of the MH2 domain (Sekelsky ; Takaesu ) (Figure 1). In these mutants, synaptic pMad was reduced to 63 ± 4% and 62 ± 4%, respectively, while nuclear pMad was reduced to 80 ± 4% and 74 ± 4% compared to control (Figure 3, C and D, quantified in M and N). A relatively small reduction in canonical BMP signaling was also evident when examining the larval fat body: these mutants had modestly reduced fat body in comparison with the control (Figure 3, I–L). The greater reduction in synaptic pMad observed in moderate alleles suggests that synaptic BMP signaling is sensitive to Mad net levels and to the integrity of the Mad-Type-I receptor interface, including the L3 loop.

Mad is the only allele with a predominant effect on nuclear pMad

Mad is the only known Drosophila allele with a mutation in the DNA binding domain of Mad, (Q90L within the MH1) (Takaesu ). Mad has been classified as a gain-of-function allele with dominant-negative activity because it exhibited an enhanced effect on dpp/dpp vein phenotypes in comparison with a deletion (Takaesu ). Interestingly, Mad was one the few alleles in our study with relative reduction of nuclear pMad signals exceeding that of synaptic pMad: in Mad mutants, the pMad levels were reduced to 46 ± 3% and 32 ± 5% compared to control for synaptic and nuclear pMad, respectively (Figure 4, A–D). This likely reflects the molecular nature of this mutation that may disrupt Mad binding to DNA more than Mad interacting with other proteins. The fat body was also strongly reduced in Mad mutant, consistent with severe defects for canonical BMP signaling (Figure 4, E and F).

Figure 4

Mad is the only allele with a predominant effect on nuclear pMad. (A–D) Confocal images of third instar NMJ 6/7 boutons and ventral nerve cord from control and Mad mutants labeled for pMad (red) and HRP (green) (A and B) or Elav (green) (C and D). (E and F) The level of fat body is strongly reduced in Mad mutants (F) in comparison to control (E). pMad is shown in Fire-lut in the upper panels (A and B) and is quantified in (G–H). Bar, 10 μm (A–D) and 1 μm (details); 100 μm (E and F). Error bars indicate SEM. ***P < 0.0001.

In our analyses, Mad and Mad alleles appeared more severe than the strong alleles noted above, Mad, Mad, Mad, and Mad. This is because Mad (or Mad)/Df(2L)C28 animals died before reaching the third instar stage, whereas all other mutants exhibited partial lethality as third-instar larvae. A similar result was obtained for Mad, a strong allele generated by the insertion of a transposon within the gene (Bellen ). We did not characterize these strong alleles further. Instead, in our search for molecular determinants relevant for pMad accumulation at synaptic junctions, we turned our attention to alleles that differentially modulate different BMP signaling pathways.

Mad exhibits a disproportionate reduction of synaptic pMad

Aside from Mad and the loop L3 mutants, most alleles showed a relatively proportional reduction in both synaptic and nuclear pMad levels (Table 1). A notable exception was Mad, which codes for a protein with a single residue change, S359L, outside of any known functional motif (Takaesu ). Mad exhibited a very significant reduction in synaptic pMad (to 54 ± 5% compared to control), but a more modest reduction in nuclear pMad (to 75 ± 5%) (Figure 5, A–D, Table 1). Moreover, there was a mild reduction of the fat body in Mad mutants in comparison to control third-instar larvae (Figure 5, E and F). We further confirmed the small effect of the Mad allele on nuclear pMad signals by examining MN expression of twit (), a gene regulated by BMP signaling (Kim and Marques 2010). Using a twit::GFP insertion (Venken ; Sulkowski ), we found that Twit::GFP levels were reduced to 70 ± 4% in the Mad mutant MNs in comparison to the control (Figure 5, I–K). These data are consistent with a relatively moderate reduction of nuclear pMad in Mad larval MNs.

Figure 5

Disproportionate reduction of synaptic pMad in Mad mutants (A–D) Confocal images of NMJ 6/7 boutons and ventral ganglia in control and Mad third-instar larvae labeled for pMad (red) and HRP (green) (A and B) or Elav (green) (C and D). (E and F) Compared with control (E), fat body levels appear modestly reduced in Mad mutant larva (F). (G–I) Quantitative analysis of synaptic pMad (G), nuclear pMad (H) and Twit::GFP (I). (J and K) Confocal images of ventral nerve cord of indicated genotypes showing MN nuclei labeled with Twit::GFP (green) and Elav (red). Twit::GFP, a MiMIC-generated chimera, provides an additional read-out for canonical BMP signaling. Genotypes: control (Mi(MIC)twit/+), Mad (Mad/Df(2L)C28, Mi(MIC)twit). Bar, 10 μm (A–D) and (J and K), 1 μm (details), 100 μm (E and F). Error bars indicate SEM. ***P < 0.0001.

Since synaptic pMad mirrors the postsynaptic type-A glutamate receptors, we next examined the GluRIIA levels relative to total glutamate receptors, stained for GluRIIC, a subunit common for both type-A and type-B receptors (Marrus ). GluRIIA levels were significantly reduced at Mad NMJs (to 67 ± 5% compared to control, n = 30) (Figure 6, A and B, quantified in 6C). In contrast, synaptic GluRIIC levels appeared only modestly diminished at Mad NMJs comparing to control (to 89 ± 6%, n = 30). This predominant GluRIIA reduction correlates with the disproportionate reduction of synaptic pMad at Mad NMJs.

Figure 6

Mad mutants have altered postsynaptic iGluRs composition. (A and B) Confocal images of NMJ6/7 boutons from third instar larvae of control and Mad mutants labeled for HRP (blue), GluRIIA (green), and GluRIIC (red or black). The Mad NMJs show severe reduction of GluRIIA synaptic signals but modest decrease of GluRIIC levels (quantified in C). The bouton area and the number of synaptic contacts (GluRIIC-positive puncta) are similarly reduced in Mad mutants compared to control (quantified in D). Bar, 10 μm, 1 μm (details). Error bars indicate SEM. ***P < 0.0001; *P < 0.05; ns, P > 0.05.

To parse the effects of canonical and synaptic BMP signaling, we further examined bouton area and the number of glutamate receptor fields at Mad terminals (Table 1 and Figure 6A’’’–B’’’, quantified in 6D). Both metrics were significantly reduced in Mad mutants, to 64 ± 4% of control for bouton area and to 70 ± 3% of control for the number of GluRIIC-positive puncta per NMJ. This follows very closely the observed reduction in canonical BMP signaling, as captured by diminished nuclear pMad and Twit reporter intensity (Figure 5). The slightly more severe reduction in bouton area is probably caused by technical differences, since the bouton area was measured in collapsed [two-dimensional (2D)] NMJ stacks and the GluRIIC puncta in three-dimensional (3D) reconstituted NMJs. Together, our data indicate that the Mad NMJ deficits include (i) ∼30% reduction in synaptic contacts and bouton area, reflecting the reduction in the canonical BMP signaling, and (ii) > 45% reduction in the synaptic GluRIIA levels, mirroring the more profound disruption of synaptic pMad.

Electrophysiology recordings confirm strong deficits of local BMP signaling at Mad NMJs

GluRIIA-containing (type-A) receptors are determinants of the quantal size (the postsynaptic response to the release of single presynaptic vesicles) (DiAntonio ). The reduced GluRIIA/GluRIIC ratio at Mad NMJs predicts that this mutant will have diminished quantal size. We tested this prediction by recording spontaneous junction currents from muscle 6 of control and Mad mutant third-instar larvae. The mEJC amplitude was significantly reduced at Mad mutant NMJs (0.59 ± 0.03 nA in Mad, n = 11, vs. 0.84 ± 0.04 nA in control, n = 11, P = 0.0003) (Figure 7, A–C). Furthermore, the decay time constant was decreased in Mad mutants (5.7 ± 0.23 ms for Mad vs. 8.51 ± 0.40 ms for control, P < 0.0001), indicating a switch to faster desensitizing receptors at Mad NMJs (Figure 7D). Since type-B receptors desensitize much faster than type-A (DiAntonio ), these data are consistent with the reduced relative levels of type-A receptors (reduced GluRIIA/GluRIIC ratio) observed at these synapses (Figure 6). In addition, the mEJC frequency was significantly reduced in Mad mutants (0.38 ± 0.08 Hz in Mad vs. 1.01 ± 0.14 Hz in control; P = 0.001) (Figure 7E). The reduction in mEJC frequency could result from defects in the presynaptic release probability and the number and density of postsynaptic glutamate receptors consistent with the observed reduction in the number of postsynaptic receptors at Mad NMJs. Also, Twit, a target of canonical BMP signaling, has been implicated in presynaptic neurotransmitter release (Kim and Marques 2012). Together, our data suggest that Mad electrophysiological deficits represent a compound phenotype that includes (i) canonical BMP signaling deficits, such as reduced postsynaptic receptors and impaired neurotransmitter release, and (ii) synaptic BMP signaling defects, such as a disproportionate loss of GluRIIA-containing receptors.

Figure 7

Electrophysiological deficits at Mad NMJs. (A and B) Representative traces of spontaneous junction currents (left) and average mEJC traces (right) recorded at 0.5 mM Ca2+ from muscle 6, segment A3, of control (A) and Mad NMJs (B). Summary graphs showing the mean amplitude (C), decay time constant (D) and mean frequency (E) of mEJCs. Error bars indicate SEM. ***P < 0.0001.

Synaptic and nuclear pMad capture different NMJ phenotypes

To test whether the reduction of synaptic pMad at Mad NMJs reflects the unique molecular nature of this mutation or is an enhanced version of a common feature in other alleles, we have expanded our in-depth analysis to include (i) Mad, the allele with a predominant effect on nuclear pMad; (ii) Mad, a strong allele; and (iii) Mad, a moderate allele with a mutation in the L3 loop (Figure 1). We examined the synaptic levels of GluRIIA and GluRIIC, the IIA/IIC ratio, bouton area, the number of synaptic contacts (GluRIIC-positive puncta) per NMJ, and the Twit::GFP levels; and compared these results to control and measurements (Tables 1, 2 and 3 and below).

Table 2

Summary of GluRIIA and GluRIIC synaptic levels

Genotype	GluRIIA			GluRIIC			IIA/IIC			N
	Mean	SEM	P	Mean	SEM	P	Mean	SEM	P
Control	100	5.31	1	100	4.50	1	100	2.09	1	19
Mad10/Df	37.65	2.15	<2.2E-16	45.08	4.37	2.3E-12	83.51	4.24	0.023	15
Mad1/Df	37.45	1.56	<2.2E-16	47.20	3.08	<2.2E-16	79.34	6.24	0.0487	26
Mad9/Df	55.65	2.28	1.088E-11	79.71	4.09	6.194E-09	69.82	2.76	6.05E-09	15
Mad8/Df	60.21	1.38	2.775E-08	88.72	1.91	0.0002	67.87	1.64	1.756E-11	30

Table 3

Summary of Twit::GFP levels and the number of synaptic contacts (IIC puncta)

Genotype	Twit::GFP				GluRIIC puncta
	Mean	SEM	P	N	Mean	SEM	P	N
Control	100	2.08	1	29	100	2.08	1	19
Mad10/Df	27.45	2.91	<2.2E-16	18	31.04	1.89	<2.2E-16	13
Mad1/Df	38.69	1.52	1.602E-15	13	38.25	1.72	<2.2E-16	26
Mad9/Df	77.49	2.48	9.144E-09	24	69.97	2.62	4.018E-12	15
Mad8/Df	69.86	3.87	2.174E-09	23	69.99	2.94	1.84E-11	30
Mad237/Df	71.48	2.56	4.921E-12	33

We found that the strong alleles, Mad and Mad, have faint GluRIIA and GluRIIC synaptic signals; and drastically reduced number of synaptic contacts (Figure 8, A–C, quantified in I and J). Interestingly, both GluRIIA and GluRIIC immunoreactivities were similarly reduced at Mad and Mad NMJs, to 37% and 47% in Mad and 31% and 37% in Mad compared to control. Consequently, the IIA/IIC ratios were diminished by only ∼20%, (Mad, to 79 ± 6% of control, n = 26; Mad, to 83 ± 4%, n = 13). In contrast, at Mad NMJs, the GluRIIA levels were reduced to 56 ± 2%, while the GluRIIC signals remained relatively robust (at 80 ± 4% of control). This unequal reduction of GluRIIA diminished the IIA/IIC ratio to 70 ± 3%, a value comparable to that observed for Mad (68 ± 2%, compare Figures 6D and 8D, quantified in 8I).

Figure 8

Synaptic and nuclear pMad capture different NMJ phenotypes. (A–D) Confocal images of NMJ6/7 boutons from third-instar larvae of indicated genotypes labeled for HRP (blue), GluRIIA (green), and GluRIIC (red or black). (E–H) Confocal images of ventral nerve cord of indicated genotypes showing MN nuclei labeled with Twit::GFP (green) and Elav (red). (I) Quantification of GluRIIA and GluRIIC synaptic signals for the indicated alleles and comparison with the synaptic pMad. The correspondence between synaptic pMad and GluRIIA signals is particularly evident in moderate Mad and Mad alleles. (J) Quantification of bouton area, number of synaptic contacts (GluRIIC-positive puncta), Twit::GFP and nuclear pMad signal intensity for the indicated alleles. Note the tight correlation of these metrics in these alleles. Bar, 10 μm, 1 μm (details). Error bars indicate SEM. ***P < 0.0001; ns, P > 0.05.

This analysis also showed that strong Mad alleles such as Mad and Mad have diminished bouton area and drastically reduced GluRIIC-positive puncta per NMJ, a direct consequence of their decreased nuclear pMad (Figure 8J and Table 1). In our hands, Mad larvae were extremely frail and showed significant lethality, suggesting our measurements underestimate the severity of Mad phenotypes. The Twit::GFP levels were likewise reduced in Mad and Mad mutants and remained apparent only in selective MNs (Figure 8, F and G, quantified in 8J). In contrast, Mad showed more robust NMJ growth, as illustrated by bouton area, the number of synaptic contacts, and modestly reduced Twit::GFP levels (Figure 8, H and J).

Notably, for each of the alleles examined in detail, the metrics describing canonical BMP signaling (nuclear pMad and Twit::GFP levels, bouton area, and the number of synapses per NMJ) changed in synchrony, in tight correlation with each other (Figure 8J). This set of metrics therefore provides a good diagnostic for the status of canonical BMP signaling. In contrast, synaptic pMad segregates away from the canonical BMP signaling metrics and follows the GluRIIA levels instead. This is particularly evident for the moderate alleles, Mad and Mad, which showed a more severe drop in GluRIIA levels than GluRIIC (Figure 8I). We suspect the same is true for strong alleles, including Mad and Mad; however, these alleles had extremely low levels of both synaptic pMad and postsynaptic glutamate receptors, impairing our ability to compare and correlate synaptic pMad with GluRIIA levels. Taking into account the shrunken NMJ areas, we estimated that the net levels of synaptic pMad, over an entire synaptic terminal, dropped to 16.8% and 4.4%, and the GluRIIC levels to 13.9% and 8.1% in Mad and Mad, respectively, compared with controls.

We have previously shown that synaptic pMad accumulates at active zones in puncta juxtaposing the postsynaptic densities (Sulkowski ). This distribution near the presynaptic membrane together with the short life-time of pMad prompted us to propose that synaptic pMad (i) is generated locally by active BMP/BMPR signaling complexes, and (ii) remains associated with its own kinase, the activated Type I BMPR Tkv, in complexes anchored at the active zone by trans-synaptic interactions. Such signaling complexes, comprised of ligands bound to tetrameric (Type I and II) BMPRs and pSmads, have also been observed at vertebrate cellular junctions (Eom ). At the fly NMJ, Mad mutations that impact the Mad-BMPR interface should disrupt not only the binding and phosphorylation of Mad, but also the accumulation of pMad-BMPR complexes at synaptic sites. The Mad allele codes for a single residue substitution (S421L) in the L3 loop at the interface between Mad and Tkv. This mutant is predicted to (i) bind to Tkv suboptimally, which can account for its reduced nuclear pMad levels, and to (ii) dissociate faster, which can further diminish its synaptic pMad. Indeed, our analysis of Mad NMJs captures both phenomena. In addition, our comparison of multiple alleles confirms that synaptic and nuclear pMad specify distinct NMJ features. In particular, the correlation between synaptic pMad and GluRIIA is consistent with our previous findings that synaptic pMad is part of a positive feedback loop that functions to stabilize GluRIIA-containing receptors at synaptic sites (Sulkowski , 2016).

Biochemical analyses reveal faster dissociation of Mad8-Tkv complexes

To understand how Mad can have disproportionate influence on synaptic pMad accumulation, we examined its phosphorylation levels using a tissue culture-based signaling assay (Serpe and O’Connor 2006). We transfected S2 insect cells with either wild-type or S359L Mad constructs and presented them with 0.1 nM Dpp, the fly ortholog of BMP2/4. Dpp triggered a robust increase in pMad levels compared with the nontreated control (Figure 9, A and B); however, pMad accumulation in cells transfected with Flag-Mad-S359L, called here Mad8, was only 50% when compared to control. This result could not be explained by a disruption of Mad levels; western blot analysis indicated similar levels of the internal control, the Flag epitope. Instead, the difference in pMad signals suggests that the S359L substitution impacts the ability of Mad8 to bind to, and/or be phosphorylated by, activated Tkv.

Figure 9

Biochemical analysis of Mad-Tkv interaction (A) Western blot analysis of whole extracts from S2 cells expressing Flag-Mad and Flag-Mad8 and treated with Dpp as indicated. Signaling is measured by the relative levels of pMad (red) and Flag (green) and is quantified in the histogram (A’). Compared to control, Mad8 has significantly reduced pMad levels upon Dpp exposure. (B) Confocal images of S2 cells transfected with Flag-Mad variants and Tac-Tkv chimeras as indicated, spread on anti-Tac coated surfaces and labeled for pMad (red), Flag (green), actin (phalloidin - magenta), and DNA (DAPI - blue). The activated Tac-TkvA chimera induces nuclear pMad accumulation when cotransfected with Flag-Mad, and, to a lesser extent, with Flag-Mad8, as quantified in (B’). pMad signals localize to cell surfaces only in Tac-TkvA/Flag-Mad cotransfected cells. (C) Structure of the Type I receptor (PDM code 3TZM). The L45 loop (magenta) interacts specifically with R-Smads. The N-lobe of the receptor, including the GS box and L45, forms a docking surface for MH2 and positions the S-V-S C-tail of Mad in the catalytic pocket. (D and E) Structure of the MH2 domain of Drosophila Mad (PDM code 3DIT) shown as monomer (D) and trimer (E). The L3 loop (yellow) is engaged in exclusive interactions with either the L45 loop of the receptor or the phosphorylated C-tail of another MH2 domain. (F) Map of MH2 Mad residues mutated in various alleles. The two views of the structure are related by a 90° rotation around a vertical axis. (G) S359L in silico mutagenesis. S359 and its adjacent peptide backbone form hydrogen bonds with H357 and residues on the H2 helix (purple); the S359L substitution breaks the hydrogen bonds with H357 (red asterisk) and introduces a bulky moiety, shifting the H2. (H) Lateral view of the Mad MH2 trimer showing the close proximity of the H2 helix (purple) to the charged L3 surface (colored by atoms). (I) Alignment of R-Smad sequences indicating class specific residues in the H2 region, including the S359 (yellow). Bar, 10 μm (B). Error bars indicate SEM. ***P < 0.0001; **P < 0.001; *P < 0.05.

Since Mad predominantly disrupts synaptic pMad, we hypothesize that Mad8 also prematurely dissociates from BMP/BMPR signaling complexes at the cell membrane. To test this possibility, we generated a Tac–Tkv chimera, where the Tkv intracellular domain was fused with the N-terminal and transmembrane domain of human Tac, the low affinity IL-2 receptor alpha chain (Ren ). We also generated a Tac-activated Tkv chimera, referred here as Tac–TkvA, by introducing the Q199D substitution at the end of the GS box of Tkv; this substitution induces activation of the Ser/Thr kinase domains within TGF-β Type I receptors (Wieser ).

In the presence of Tac–Tkv, both Mad and Mad8 were detectable throughout cells, including at cell membranes (Figure 9B); however, no pMad immunoreactivity was visible, indicative of an inactive Ser/Thr kinase. In the presence of activated Tac–TkvA, cotransfection of Mad induced massive pMad accumulation throughout the cell. In contrast, cotransfection of Mad8 elicited only nuclear pMad accumulation. Quantification of the pMad levels showed a drastic reduction of nuclear pMad for Mad8 (to 50% compared to control, Mad transfected cells), and an even more profound reduction of pMad at cell membranes (to ∼30% of control levels). Similar to signaling assays (Figure 9A), these reduced pMad levels could reflect (i) a disruption of the Tkv–Mad8 enzyme-substrate (E-S) interaction by reducing the E-S association (via reduced binding and/or increased dissociation) or increasing the energy of activation, or (ii) an incorrect, premature dissociation of the enzyme-product (E-P) complexes. Since pMad immunoreactivity at the cell membrane is diminished beyond the level in nuclei, the S359L substitution likely induces a faster dissociation of the Tkv-pMad8 (E-P) complexes. This interpretation is consistent with our in vivo findings that Mad mutants disproportionally affect pMad accumulation at synaptic terminals.

At this time, structures for Smad-Type I receptor complexes are not available, but site-directed mutagenesis and studies on related structures indicate that the L3 loops of R-Smads bind to the L45 loops of the receptors (Huse ; Durocher ) (Figure 9, C and D); residues on both L3 and L45 confer class-specificity during signaling, i.e., BMP vs. Activin pathway. Also, the GS domain and N-terminal lobe of the receptor kinase appear to form an MH2 docking interface that positions the Smad SSXS C-tail adjacent to the catalytic pocket of the receptor kinase. An opposing concave interface has been observed for all Smad MH2 domains, including the Drosophila Mad MH2. Once phosphorylated, pSmads have a high propensity to trimerize that favors their dissociation from receptors (Kawabata ). Within a trimeric complex, the L3 loop engages the phosphorylated pS-X-pS tail of an adjacent Smad (Figure 9E). In addition, the L3 loop of R-Smads has been implicated in mutually exclusive interactions with the Type I receptors or with the phosphorylated pS-X-pS tail of other Smads in trimeric complexes (Wu ). This model is also consistent with our experimental observations: Strong L3 mutations could effectively impair BMP signaling by disrupting both complexes, Mad-Tkv and Mad hetero-trimer complexes, while moderate alleles may cause more limited disruptions of Mad-Tkv interaction.

S359, the residue mutated in Mad8, maps within the same concave cavity of MH2 as the L3 loop and follows a highly conserved helix, H2 (Figure 9, F and I). Inspection of this region indicates that S359 stabilizes the highly conserved H357 via hydrogen bonds. Additional hydrogen bonds connect negative charges of the peptide backbone with N349 and Q346. Together these interactions appear to anchor and stabilize the H2 helix. Similar to L3, the H2 helix contains class-specific residues: Y352HH in Smads of the BMP pathway, and QRY in the equivalent position in Smads of the Activin pathway (Figure 9I). The class-specific features of the H2 helix and its spatial orientation and proximity to the L3 loop suggest that H2 is critical for the Smad-Type I receptor interaction. In this scenario, the S359L substitution in Mad8 disrupts the Mad-Tkv interaction by (1) replacing a polar interface residue with a bulky hydrophobic one, (2) destroying the hydrogen bond S359-H357, and (3) likely shifting the position of H2, and the YHH class-specific residues, relative to the Mad-Tkv interface. This is consistent with our in vivo findings that the S359L substitution induces a predominant reduction of synaptic pMad.

Discussion

Our comprehensive analysis of numerous Drosophila alleles revealed a rich phenotypic complexity for this important gene during different BMP signaling modalities and uncovered a new molecular determinant that confers class-specificity during TGF-β signal transduction.

Using the Drosophila NMJ system, we have quantified the effect of alleles on both canonical and synaptic BMP signaling pathways. Drosophila requires canonical BMP signaling for NMJ growth. This pathway is triggered by muscle-derived BMPs that bind their receptors on MNs and form active BMP signaling complexes that are endocytosed and trafficked retrogradely to MN soma where they phosphorylate Mad (Marqués and Zhang 2006). pMad accumulation in MN nuclei quantitatively captures this canonical BMP response (Kim and Marques 2010). In addition, pMad accumulates at synaptic terminals as a result of a local, Smad-dependent BMP pathway that sculpts postsynaptic composition as a function of synapse activity (Sulkowski , 2016). We have previously shown that synaptic pMad is generated locally by BMP signaling complexes confined to the active zone, a region protected from endocytosis (Sulkowski ). This suggests that MNs can monitor synapse status and then deploy BMP receptors to (i) the active zone, where they strengthen the synapse via local BMP signaling, or (ii) outside the active zone, where they engage in canonical signaling to expand the NMJ.

As expected, strong alleles drastically reduce both nuclear and synaptic pMad levels (Figure 2 and Table 1). In contrast, moderate alleles reveal differential requirements for different signaling modalities. First, uniform reduction of Mad levels preferentially reduced local pMad, indicating that pMad accumulation at synaptic junctions is more sensitive than pMad accumulation in MN nuclei (note Mad and Mad, Figure 3 and Table 1). This reduction may reflect (i) the distribution of Mad protein within the MNs, presumably more abundant in the soma than in neurites, or (ii) the nature of the Mad-BMPRs complexes, more transient in the soma but more stable at the synaptic junctions. Second, the prominent reduction of nuclear pMad in Mad mutants reinforces the idea that the MH1, the DNA binding domain of Mad affected in this mutant, is critical for nuclear accumulation of pMad and proper transcriptional control (Figure 4 and Table 1). Third, mutations in the L3 loop of the MH2 domain produced a set of moderate alleles (Mad and Mad) with predominant deficits in synaptic pMad, but also produced two of the strongest Mad alleles (Mad and Mad) (Figures 2 and 3). Since the L3 loop of R-Smads engages in a mutually exclusive interaction with the Type I receptors or with the phosphorylated C-tail of other Smads, we reason that strong L3 mutations may disrupt both sets of L3-mediated interactions, while moderate alleles may cause primarily Mad-Tkv dissociation. Finally, the disproportionate deficits of synaptic pMad in Mad revealed a role for the H2 helix in the modulation of the Mad-Tkv interaction. This highly conserved helix resides adjacent to the L3 loop and includes additional class-specific residues. These H2 features and its relative position to the MH2 concave cavity suggest that the H2 region is critical for the Smad-Type I receptor interaction. H2 maps outside the Smad trimer interface (Figure 7H), but it may further interfere with Smad-dependent transcription via Smad/cofactor interactions. Our in vivo findings favor a role for H2 in the modulation of Mad-Tkv interaction. The H2 contribution may be direct, by shaping the Mad-Tkv interface, or indirect, via recruiting other protein(s) that may stabilize Mad-Tkv complexes at cell junctions.

Indeed, other proteins appear to contribute to the stabilization of BMP signaling complexes at specialized cell junctions. Previous studies in the chick neural tube have shown that BMP signaling controls apicobasal polarity partly by enabling the pSmad1/5/8-dependent association of BMP signaling complexes with the PAR3-PAR5-aPKC complex at the tight junctions (Eom ). Reduced Smad phosphorylation destabilizes the PAR complex and disrupts the tight junctions. In flies, loss-of-function disruptions of Bazooka(Par-3)-Par-6-aPKC complexes produced NMJs with significantly reduced number of boutons and increased levels of postsynaptic GluRIIA receptors (Ruiz-Canada ). These phenotypes have been attributed to severe disruptions of microtubule stability, which may obscure a role for these complexes in modulation of BMP signaling.

In addition, posttranslational modifications of Mad and Tkv may strengthen or weaken their interaction and selectively modify nuclear and synaptic pMad. For example, Mad phosphorylation at S25 by Nemo (i) disrupts Mad association with Tkv at synaptic terminals, and (ii) favors nuclear export of pMad (Zeng ; Merino ). In both flies and mammals, posttranslational modifications such as ubiquitination limit receptor activity and trigger degradation or deactivation of both Smads and receptors, keeping these signaling components in check (Zhu ; Dupont , Stinchfield ). Interestingly, disruption of a number of components of the endocytic machinery elevate synaptic pMad and diminish nuclear pMad presumably by shifting receptor allocation from canonical BMP signaling to local, noncanonical signaling (O’Connor-Giles ; Vanlandingham ; Heo ).

Misregulation of BMP signaling is associated with many developmental abnormalities and disease states. Our study indicates that the H2 helix region is critical to local BMP signaling and suggests that mutations in this region may disrupt specialized tight junctions throughout the animal kingdom. We searched for human genetic variants in the H2 region using the MARRVEL database (Model Organism Aggregated Resources for Rare Variant ExpLoration) (Wang ). Interestingly, a Smad1 point mutation, Y362C that changes a class-specific residue within H2, was reported in the relative of a patient with colonic atresia, a congenital intestinal malformation that results in failure to pass meconium in newborns. Also, N361D substitution in Smad5 was found in a patient with malformation of the heart and great vessels. Finally, a N361-G365 deletion in Smad5 was reported in a patient with epileptic encephalopathy. While further studies are required to elucidate the functional impact of these genetic alterations on synaptic BMP signaling, these variants underscore the relevance of the H2 helix for normal development and function.