|Home | About | Journals | Submit | Contact Us | Français|
Communication between pre- and post-synaptic cells is a key process in the development and modulation of synapses. Reciprocal induction between pre- and postsynaptic cells involves regulation of gene transcription, yet the underlying genetic program remains largely unknown. To investigate how innervation-dependent gene expression in postsynaptic cells supports synaptic differentiation, we performed comparative microarray analysis of Drosophila muscles before and after innervation, and of prospero mutants, which show a delay in motor axon outgrowth. We identified 84 candidate genes that are potentially up- or downregulated in response to innervation. By systematic functional analysis, we found that one of the downregulated genes, longitudinals lacking (lola), which encodes a BTB-Zn-finger transcription factor, is required for proper expression of glutamate receptors. When the function of lola was knocked down in muscles by RNAi, the abundance of glutamate receptors (GluRs), GluRIIA, GluRIIB and GluRIII, as well as that of p-21 activated kinase (PAK), was greatly reduced at the neuromuscular junctions (NMJs). Recordings of the synaptic response revealed a decrease in postsynaptic quantal size, consistent with the reduction in GluR levels. Lola appears to regulate the expression of GluRs and PAK at the level of transcription, because the amount of mRNAs encoding these molecules was also reduced in the mutants. The transcriptional level of lola, in turn, is downregulated by increased neural activity. We propose that Lola coordinates expression of multiple postsynaptic components by transcriptional regulation.
During initial synapse formation, reciprocal interaction between innervating neurons and their targets are essential for assembly of synaptic components, cytoskeletal organization and activation of gene expression (Goda and Davis, 2003; Li and Sheng, 2003; McAllister, 2007). Similarly, mutual trans-synaptic signaling is important for activity-dependent refinement of differentiated synapses (Kandel, 2001; Flavell and Greenberg, 2008). While short-term changes in synaptic properties can be induced by modulation of pre-existing proteins and/or mRNAs, such as trafficking, modification and local translation, long-term changes require transcriptional control of gene expression. During the formation of vertebrate neuromuscular junctions (NMJs), signals from presynaptic motor neurons are necessary for the regulation of gene expression of postsynaptic transmitter receptors, which are acetylcholine receptors (AChRs) (Sanes and Lichtman, 2001; Schaeffer et al., 2001; Kummer et al., 2006). Many studies have focused on the role of immediate-early genes (IEGs) and CREB-mediated transcriptional regulation in long-term synaptic plasticity and memory formation (Flavell and Greenberg, 2008; Cohen and Greenberg, 2008). Nonetheless, a large gap still exists in our knowledge about how multiple molecular pathways integrate and orchestrate the development and plasticity of synapses. In particular, despite extensive work on activity-induced genes, very few studies have established functional links between these activity-induced genes and the downstream target genes that ultimately regulate synaptic properties.
In this study, we used the Drosophila NMJ as a model to study gene expression changes in postsynaptic muscle cells in response to presynaptic innervation. The Drosophila NMJ is a glutamatergic synapse expressing ionotropic glutamate receptors (GluRs) and contains a number of synaptic components commonly found in mammalian synapses, such as the postsynaptic density protein Discs-Large/PSD-95 (Keshishian et al., 1996; Griffith and Budnik, 2006). Previous studies showed that immediate–early transcription factors such as CREB and AP-1 regulate the strength and/or morphology of this synapse (Davis et al., 1996; Sanyal et al., 2002) (reviewed in Sanyal and Ramaswami, 2006). Signaling pathways mediated by secreted factors, such as Wnts and Bmps, are known to regulate anterograde and/or retrograde interaction between the motor neurons and muscles that are important for synaptic development (McCabe et al., 2003; Ataman et al., 2008; Korkut et al., 2009) (reviewed in Griffith and Budnik, 2006). However, the final targets of these signaling cascades—the molecules that directly regulate the changes in synaptic structure and function—remain largely unknown.
In this study, we performed genome-wide microarray analyses of specific muscle cells and identified 84 candidate genes whose expression changed in response to innervation. By systematic functional analyses of the candidate genes, we found that longitudinals lacking (lola), a gene downregulated by innervation, plays a prominent role in the transcriptional control of a number of postsynaptic components. lola encodes a BTB-Zn-finger transcription factor with a number of different isoforms (Goeke et al., 2003; Horiuchi et al., 2003). This transcription factor, Lola, has been implicated in a wide range of developmental and cellular processes including axon guidance, neural specification and tumorigenesis (Madden et al., 1999; Crowner et al., 2002; Goeke et al., 2003; Ferres-Marco et al., 2006; Spletter et al., 2007). Previous studies suggest that Lola may execute its function by directly binding to DNA and regulating the expression of the target genes. Here we show that postsynaptic Lola transcriptionally regulates the expression level of the glutamate receptors GluRIIA, GluRIIB and GluRIII, as well as p-21 activated kinase (PAK). We also show that the transcriptional level of lola is downregulated by increased neural activity. We propose that postsynaptic Lola functions as a transcription factor that controls synapse formation and/or maturation by regulating the expression of multiple synaptic components.
For microarray analysis, we used an allele of prospero (prosM4) (Broadie and Bate, 1993b) and y, w. Forced expression and RNA interference (RNAi) analyses were performed using the GAL4-UAS system (Brand and Perrimon, 1993). Elav-gal4 (Luo et al., 1994), 24B-Gal4 (Brand and Perrimon, 1993) or G14-Gal4 (Shishido et al., 1998), or 5053A-Gal4 (Ritzenthaler et al., 2000) were used to induce expression in all neurons, all muscles, or in M12, respectively. UAS-RNAi lines were obtained from the Vienna Drosophila RNAi Center (VDRC) and Fly Stocks of the National Institute of Genetics (NIG) (Dietzl et al., 2007). Lines and alleles used for the systematic functional analyses are listed in supplementary material Table S2. Animals were raised at 29°C for RNAi analyses and at 25°C for other analyses. Alleles of lola, lolaORC46 (Crowner et al., 2002) and lolaORE119 (Goeke et al., 2003), were used.
The collection of embryonic somatic muscles was performed as previously described (Inaki et al., 2007) with the following modifications. For collection of muscles at 18 hr after egg laying (AEL), the preparation was treated with 1 mg / ml collagenase (Sigma, St. Louis, Missouri) for ~30 sec to weaken intersegmental muscle-muscle attachments. For chip analysis, we prepared three samples for each developmental stage or genotype: one sample was prepared from 200 collected cells and two from 50 collected cells. We performed two or three rounds of cRNA amplification before biotin labeling for samples prepared from 200 cells or those prepared from 50 cells, respectively. Synthesis of biotin-labeled cRNA was performed with the IVT Labeling Kit (Affymetrix, Santa Clara, CA). Hybridization of the fragmented and labeled cRNA to Affymetrix Drosophila Genome 2.0 Genechip arrays were performed according to the Genechip Expression Analysis technical manual (Affymetrix). Gene expression data analyses were performed using Affymetrix GeneChip Operating Software 1.4 (Affymetrix). By comparing gene expression signals between the samples prepared from the same number of cells, we obtained the Change Call (I, increase; D, decrease; NC, no change). After three pairs of comparisons, we selected the genes that displayed Change Call ‘I’ in all three pairs. We carried out Gene Ontology analyses by batch query using NetAffx (Affymetrix).
Real-time reverse-transcription PCR (RT-PCR) assays were performed with ABI Prism 7000 SDS or Applied Biosystems StepOne Real Time PCR System (Applied Biosystems, Foster City, CA) with SYBR Green fluorescence according to the manufacturer protocol. We used as templates the second-round amplified cDNAs of wild type muscles and pros mutant muscles, prepared as described for microarray analyses. The gene expression values were normalized using the Myosin heavy chain (Mhc) gene as a reference. To examine gene expression in larval muscles, total RNA was extracted from the body wall of wall-climbing third-instar larvae using Sepasol RNA I Super (Nacalai, Kyoto, Japan) according to the manufacturer instructions. Approximately 2 µg of total RNA was reverse-transcribed into cDNA using oligo-dT primer and SuperScript III Reverse Transcriptase (Invitrogen, Carlsbad, CA). The obtained values were normalized using the Gapdh1 gene as a reference.
Immunohistochemical staining of dissected larvae and embryos was performed as described previously (Inaki et al., 2007). The following primary antibodies were used at the indicated concentrations: goat anti-horse radish peroxidase (HRP; Jackson, West Grove, PA, 14000); rabbit anti-Lola polyclonal antibodies (1150) (Giniger et al., 1994); mouse monoclonal anti-GluRIIA (8B4D2) and mouse monoclonal anti-Dlg (Developmental Studies Hybridoma Bank, University of Iowa, Iowa City, IA, 110 and 150, respectively); rabbit anti-GluRIIB and rabbit anti-GluRIII (12500, and 15000, respectively) (Marrus et al., 2004); mouse Nc82 (1100) (Wagh et al., 2006); rabbit anti-Pak antibody (1500) (Sone et al., 2000).
We quantified the morphology of the NMJ with Imaris software (Bitplane, Zurich, Switzerland). The number of boutons and branches and terminal length were scored based on the morphology visualized with anti-HRP staining as previously described (Coyle et al., 2004). The terminal length was normalized to muscle area. The bouton and branch number were also normalized to muscle area in the analyses of pst RNAi mutants because there was a slight change in muscle area in the mutants (surface area, 84.3±2.8 [n=10] in pst-RNAi-8588R-4 [p<0.05] and 66.4±2.7 [n=17] in pst-RNAi-107243 [p<0.05] compared to 75.0±3.0 [n=32] in control, × 103 µm2). To assess the content of particular molecules at synapses, confocal images were projected to create a two-dimensional image and analyzed with IPLab software (Scanalytics, Fairfax, VA, USA). To control for relative intensity, experimental and control samples were stained in the same dish and imaged under identical conditions. We defined the synaptic area as delimited by HRP immunoreactivity on the muscle surface and the synaptic intensity of GluR, PAK or BRP as the average immunoreactive signals within the synaptic area. The total area of GluRIIA-, PAK- and BRP-positive clusters were defined by the pixels with an intensity above a threshold set arbitrarily, and were normalized by the total area of muscles. The GluRIIA expression level of the embryonic NMJ was defined as the total immunoreactive signals within the area of GluRIIA-positive clusters. For analysis of Dlg levels, we measured net intensity of immunoreactive signals normalized by the stained area. The extra-junctional GluR expression level was defined as the average immunoreactive signal within arbitrarily chosen regions covering ~1% of muscle surface area. All statistical analyses were done by Student's t-test. All quantitative data in the larvae were obtained from muscle 6/7 NMJs in abdominal segment 2.
Two hundred first-instar larvae were ground in RIPA buffer containing protease inhibitors (Complete mini; Roche Diagnostics, Mannheim, Germany) as described by Thomas et al. (Thomas et al., 1997) and loaded onto 4%–12% SDS-PAGE gels (TEFCO, Tokyo, Japan). Immunoblotting was performed using mouse monoclonal anti-Tubulin (12G10; Developmental Studies Hybridoma Bank, 12000), anti-GluRIIA (MAb 8B4D2 described above, 150), goat anti-mouse IgG HRP conjugate (Jackson, 11000) and donkey anti-goat IgG IRDye 680(LI-COR Biosciences, Lincoln, NE, USA, 12500) and imaged and quantified by Odyssey infrared imaging system (LI-OCR). To control for loading, the level of GluRIIA was normalized to the level of Tubulin.
Intracellular recordings were made from muscle 6 in abdominal segments 3 and 4, and excitatory junctional potentials (EJPs) and miniature EJPs were recorded as described previously (Morimoto et al., 2010). HL3-containing medium at varying concentrations of extracellular Ca2+ were used (0.375 mM for mEJP recordings and 1.5 mM for EJP recordings). The peak amplitudes of ten EJPs were measured and averaged in each recording. All statistical analyses were done by Student's t-test.
Embryos of elav-Gal4/UAS-dTRPA1 (Pulver et al., 2009) were raised at 22°C, collected at 18 hr AEL and transferred to a PCR tube containing 25 µl of Halocarbon 700 oil (Sigma). Repetitive, spaced stimuli, consisting of 2-min intervals at a higher temperature (27°C) spaced by 5 min of rest (22°C) were applied in a PCR machine for 12 hrs. The larvae were then transferred to an agar plate and reared at 22°C until third instar. Control embryos (elav-Gal4/+) underwent the same temperature cycling.
In the Drosophila NMJ, synaptogenesis starts at ~13 hr AEL, when motor neurons contact the target muscles for the first time (Broadie and Bate, 1993a). During the process of intimate interaction between the presynaptic motor neurons and the postsynaptic muscles, the synapses gradually mature and become fully functional by 17 hr AEL (Broadie and Bate, 1993a; Kohsaka et al., 2007). Using Affymetrix DNA microarrays, we performed two series of comparative gene expression analysis of muscles (Fig. 1A). First, we compared gene expression profiling of muscles before (12.75 hr AEL) and after (18 hr AEL) innervation (“before versus after” comparison). Second, we compared gene expression in muscles of wild-type versus prospero (pros) mutants; ~80% of the muscles in pros mutants lack innervation but otherwise differentiate normally (“with versus without” comparison, at 18 hr AEL) (Broadie and Bate, 1993b). Genes that are differentially expressed in the first comparison would include genes that are up- or downregulated by innervation. However, these would also include genes that are involved in other aspects of muscle differentiation and whose expression is regulated at this stage, such as genes involved in the maturation of the muscle contraction apparatus. On the other hand, the second comparison would identify genes that are specifically regulated by innervation because the differentiation of muscle itself is normal in pros mutants (Broadie and Bate, 1993b). By conducting comparative analyses based on these two categories, we aimed to identify genes that function in synapse formation and whose expression is regulated by innervation.
Using the single-cell collection technique that we developed previously (Inaki et al., 2007), we performed microarray analyses of four neighboring ventral muscles, M6, M7, M12 and M13, using Affymetrix Drosophila Genome chips (see Materials and Methods for detail). Because these muscles show similar morphology, run in parallel and attach to adjacent muscle insertion sites, they likely share common genetic programs for synapse formation. We therefore used a combination of the four muscles as a source for the microarray analyses. We prepared three independent samples each for the two comparative expression analyses described above (“before versus after” and “with versus without” innervation) and focused on genes that displayed differential expression reproducibly. We identified 72 genes that showed increased signals in muscles “after” and “with” innervation compared respectively to muscles “before” and “without” innervation. We defined these as “genes upregulated by innervation” (Fig. 1B; Table 1). We also identified 12 genes that showed decreased signals in the same comparisons and defined them as “genes downregulated by innervation” (Fig. 1C; Table 1). To validate the microarray results, we analyzed the expression of randomly-chosen subsets of the identified genes by quantitative RT-PCR (qPCR). Differential expression (“with vs. without”) was confirmed for 22 of the 26 genes assayed, indicating 85% concordance (Table 2).
The 72 upregulated and 12 downregulated candidate genes described above were categorized into multiple functional classes of proteins (supplementary material Table S1). Among the upregulated genes was cactus (cact), one of the few genes known to function in muscles to regulate synaptic development and/or function (Cantera et al., 1999; Heckscher et al., 2007). To study whether other candidate genes also have roles in synaptic development, we performed loss-of-function (LOF) analyses on 37 of the candidate genes whose predicted molecular structure and/or function are related to synaptogenesis (supplementary material Table S2), using existing LOF mutant alleles and RNA interference (RNAi)-mediated gene knock-down. We visualized the NMJ by staining the presynaptic membrane (anti-HRP) and postsynaptic GluRIIA, and looked for changes in the terminal morphology and/or receptor clustering in the mutants. This screening revealed a function for pastrel (pst) in proper synaptic growth and lola in proper GluRIIA expression. No dramatic phenotypes were observed in the LOF mutants of the remaining genes examined.
RNAi knock-down of pastrel (pst) with two independent constructs that target different portions of the mRNA resulted in decreased bouton number, branch number and NMJ length (Fig. 2) (number of boutons per muscle area, 1.1±0.094 in pst-RNAi-8588R-4 [p<0.0005] versus 1.7±0.072 in control, × 10−3/µm2; number of branches per muscle area, 0.24±0.021 in pst-RNAi-8588R-4 [p<0.05] and 0.24±0.012 in pst-RNAi-107243 [p<0.0005] versus 0.32±0.014 in control, × 10−3/µm2; NMJ length per muscle area, 6.9±0.3 in pst-RNAi-8588R-4 [p<0.0005] and 8.9±0.3 in pst-RNAi-107243 [p<0.05] versus 9.7±0.3 in control, × 10−3 µm−1; n=10, 17, 32).
pst has been implicated in the processes of learning and memory (Dubnau et al., 2003). The expression of pst appears to be upregulated with innervation. Therefore, pst may be part of the genetic program that is induced by innervation to promote synaptic growth.
The systematic functional analyses described above revealed that lola plays a critical role in the regulation of postsynaptic GluRIIA (Fig. 3). lola encodes a BTB-Zn-finger transcription factor with at least 20 unique protein isoforms, some of which are expressed in somatic muscles (Goeke et al., 2003). In our microarray screen, lola was identified as a gene downregulated by innervation; the expression level was 5-fold lower in 18-hr AEL compared to 12.75 hr AEL and 4-fold higher in the pros mutants compared to wild type (Table 1). qPCR analyses further confirmed differential expression: lola expression was 1.6-fold higher in pros mutants compared to the wild type (Table 2). Thus, lola expression appears to be downregulated by innervation. Lola protein was expressed in the nuclei of muscles in 12.75 hr and 18 hr embryos and in third instar larvae (supplementary material Fig. S1); see also Giniger et al. (Giniger et al., 1994).
We knocked down lola expression in muscles using RNAi constructs. The level of lola mRNA in the body wall was reduced to 37±6.5% (n=9) of control levels in the RNAi mutants, as assessed by qPCR analysis. We observed a dramatic decrease in GluRIIA abundance in NMJs in the mutants (Fig. 3A,B). This phenotype was seen when two independent RNAi constructs that target different regions of the mRNA were used. In contrast, there were no abnormalities in the morphology of the NMJs as assessed by anti-HRP staining (Fig. 3A′,B′, Fig. 6H). Thus, lola knock-down specifically alters GluRIIA abundance without affecting the area of the NMJs. We quantified the change in the abundance of GluRIIA by measuring the fluorescence intensity of GluRIIA at the NMJ as described previously (Albin and Davis, 2004). The abundance of GluRIIA decreased to ~40% of control levels (Fig. 3G) (38±1.9% in lola-RNAi-12573 [n=15] and 38±2.3% in lola-RNAi-101925 [n=17] compared to 100±2.6% in 24B-Gal4/+ [n=27] and 98±9.3% in UAS-lola-RNAi-12573/+ [n=8]). Similar changes in GluRIIA abundance in the RNAi mutants were also observed at the late embryonic NMJ (20 hr AEL) (supplementary material Fig. S2). Furthermore, even more dramatic changes in GluRIIA abundance were observed in the lolaORC46 embryos (Fig. 4). lolaORC46 is a null allele and is embryonic lethal. In lolaORC46 embryos, motor axon projection and muscle patterning are partly deranged with some muscles entirely lacking nerve innervation (Madden et al., 1999). We therefore focused on muscles that received normally-sized nerve terminals (as visualized with a presynaptic marker, anti-HRP) and found a dramatic decrease in the density of synaptic GluRIIA (Fig. 4) (39±3.7 in lolaORC46 [n=27] compared to 100±3.0% in yw [n=101], p<0.0005). These data indicate that Lola is necessary for normal GluRIIA expression at the NMJ.
The decrease in GluRIIA abundance at the NMJ could be due to the decrease in total GluRIIA expression in muscles or to defects in proper accumulation of GluRIIA at the NMJs. If the latter is the case, then an increase in the level of GluRIIA in the extra-junctional region would be expected. However, we did not detect any such increase in the RNAi mutants. Instead, fluorescence intensity of GluRIIA in the extra-junctional region was significantly decreased in lola RNAi mutants (Fig. 3H) (70±6.3 in 24B-Gal4/UAS-lola-RNAi-12573 [n=15] and 67±6.3 in 24B-Gal4/UAS-lola-RNAi-101925 [n=17] compared to 100±3.6% in 24B-Gal4/+ [n=27] and 110±30 in UAS-lola-RNAi-12573/+ [n=8]). Thus, postsynaptic lola knock-down appears to impair GluRIIA levels at the NMJ without altering the overall distribution of the receptors. We also examined the total amount of GluRIIA expression by immunoblotting. The total amount of GluRIIA protein in lola RNAi larvae was significantly reduced compared to control (Fig. 5) (expression in lola-RNAi [24B-Gal4/UAS-lola-RNAi-101925] was reduced to 35±10% of the control [24B-Gal4/+] level [n=2]). These results strongly suggest that Lola is required for the control of total GluRIIA expression in muscles rather than its proper accumulation at the NMJs. Overexpression of lola in muscles driven by 24B-Gal4 did not induce obvious changes in NMJ morphology and level of GluRIIA expression (supplementary material Table S2).
We next asked if Lola is required for proper expression of other GluR subunits, GluRIIB and GluRIII (Featherstone et al., 2005; Marrus et al., 2004; Petersen et al., 1997; Qin et al., 2005). As is the case for GluRIIA, the fluorescence intensity of GluRIIB and GluRIII at the NMJ was dramatically decreased in the mutants (Fig. 3C–G) (GluRIIB, 60±2.1[n=15] compared to 100±3.3 [n=12], p<0.0005; GluRIII, 32±0.8 [n=12] compared to 100±3.4 [n=12], p<0.0005; data represent % value in lola-RNAi-12573 compared to 24B-Gal4/+). The abundance of GluRIIB and GluRIII at extra-junctional muscle areas was not also increased (Fig. 3H) (GluRIIB, 89±4.8 [n=15] compared to 100±2.5 [n=12]; GluRIII, 92±7.0 [n=12] compared to 100±4.7 [n=12]). Thus, Lola appears to be required for proper expression of at least three GluR subunits—GluRIIA, GluRIIB and GlulRIII—at the NMJs.
We also examined the expression of other postsynaptic components. Pak kinase (PAK) colocalizes with GluRIIA at the post-synaptic density (PSD) and is involved in the PSD assembly (Albin and Davis, 2004; Parnas et al., 2001; Rasse et al., 2005). In lola RNAi mutants, the level of synaptic PAK was significantly reduced (Fig. 6A,B,G) (synaptic intensity, 87±3.7 [n=17] compared to 100±3.2 [n=16], p<0.05; expression field area, 48±8.2 [n=17] compared to 100±13 [n=16], p<0.05). We also found a significant reduction in synaptic levels of Dlg, a scaffolding protein highly enriched at the perisynaptic muscle membrane (Fig. 6E,F) (postsynaptic intensity, 85±2.0 [n=22] compared to 100±2.4 [n=24], p<0.0005). These data indicate that postsynaptic Lola is required for proper expression of a number of postsynaptic components.
The lola RNAi mutants developed into normally sized third-instar larvae with no obvious behavioral abnormalities. No apparent abnormality was found in the morphology or size of muscles in the mutants (surface area, 82.7±2.9 [n=19] in UAS-lola-RNAi compared to 86.8±2.9 [n=22] in control, × 103 µm2 [p>0.05]). The morphology of the presynaptic terminals also appeared normal (Fig. 3A′,B′, Fig. 6H). There was no significant change in the bouton number, branch number or NMJ length (bouton number, 92±3.6 [n=17] compared to 100±4.2% [n=17]; branch number, 26±1.5 [n=19] compared to 26±1.0 [n=22]; NMJ length per muscle area, 8.1±0.4 [n=19] compared to 9.0±0.4 [n=22], × 10−3 µm−1; data represent values in 24B-Gal4/UAS-lola-RNAi-12573 compared to control [24B-Gal4/+]). We also examined the expression and distribution of an active zone marker Bruchpilot (BRP) (Wagh et al., 2006), and found no abnormality (Fig. 6C,D,G) (synaptic intensity, 92±2.6 [n=10] compared to 100±3.4% [n=16]; expression field area, 92±10 [n=10] compared to 100±6.2% [n=16]). Thus, despite the severe reduction in the expression of a number of postsynaptic components, the morphology and molecular expression of the apposing active zones appeared normal.
To investigate the physiological effects of the postsynaptic lola knock-down, we performed intracellular recordings from muscle 6 in the third-instar larvae. We tested the postsynaptic response to both spontaneous and evoked transmitter release. The mean amplitude of spontaneous miniature excitatory junctional potentials (mEJPs), or quantal size, is a measure of postsynaptic sensitivity to transmitters. A significant reduction in quantal size was observed in lola RNAi mutants (Fig. 7A,B) (0.91±0.068 mV in 24B-GAL4/UAS-lola-RNAi-12573 [n=11] compared to 1.33±0.10 mV in 24B-Gal4/+ [n=8], p<0.05). Previous studies suggested that the density of channels is an important determinant of the quantal size (Petersen et al., 1997). Therefore, our data are consistent with the decreased GluR expression observed in the imaging analysis. However, there was no significant change in the amplitude of excitatory junctional potential (EJP) in lola RNAi mutants, despite the decrease in the quantal size (Fig. 7C,D) (36.0±2.07 mV in lola-RNAi [n=11] compared to 35.8±4.58 mV in control [n=8]). Similar compensation for the decrease in quantal size has been reported for several mutants, including GluRIIA (Davis et al., 1998; Paradis et al., 2001; Petersen et al., 1997). These previous studies suggested that retrograde signals from the postsynaptic cells increased the presynaptic transmitter release to compensate for postsynaptic defects. The mechanism for the compensation in lola mutants remains to be determined.
The imaging and physiological analyses described above show that postsynaptic Lola is necessary for normal expression of GluRs and other synaptic components. Since lola encodes a BTB-Zn-finger transcription factor (Goeke et al., 2003; Cavarec et al., 1997), Lola might regulate the transcriptional level of GluRs. To assess this possibility, we examined the level of GluR mRNAs in the mutants. We extracted total RNA from the larval body wall and performed real-time quantitative PCR. We found that the levels of GluRIIA and GluRIIB mRNAs in lola RNAi mutants are dramatically reduced (Fig. 8) (GluRIIA, 42±7.9 [n=9] compared to 100±2.9 [n=9], p<0.0005; GluRIIB, 56±14.4 [n=9] compared to 100±6.4 [n=9], p<0.05; data represent % values in lola-RNAi-12573 compared to control [24B-Gal4/+]). We also found a smaller but significant decrease in the level of GluRIII (75±8.7 [n=9] compared to 100±3.5 [n=9], p<0.05). Furthermore, we found a significant reduction of PAK mRNAs (59±6.4 [n=9] compared to 100±2.9 [n=9], p<0.0005). On the other hand, the mRNA transcript level of Dlg and Myosin heavy chain (Mhc), a protein involved in muscle contraction, was not affected (Dlg, 87±11.3 [n=9] compared to 100±9.6 [n=9]; Mhc,116±8.4 [n=9] compared to 100±4.5 [n=9]). These data indicate that Lola positively regulates the transcription level of glutamate receptors and PAK but not Dlg.
Having shown that Lola is involved in the transcriptional regulation of a number of synaptic components, we next asked how the expression of lola itself is regulated. Our microarray analyses showed that the level of lola transcript is changed by innervation. Therefore, one possibility is that neural activity modulates lola expression. To address this, we used the Drosophila warmth-activated cation channel TRPA1 (dTRPA1) (Pulver et al., 2009) to activate all neurons during late embryonic and early 1st instar larval stages and studied the effect on the expression of lola in the third instar larvae. We applied repetitive 2-min intervals at a higher temperature (27°C) spaced by 5 min of rest (22°C), for 12 hrs from AEL18. Such spaced induction of increased activity with Channelrhodopsin or high K+ has previously been shown to induce changes in synaptic structure and function in larval NMJs (Ataman et al., 2008). The spaced stimuli with dTRPA1 increased the number of boutons (Fig. 9A–C), indicating that our stimulation protocol successfully induced plastic changes in the NMJs. The stimulation protocol also decreased the level of lola mRNAs (Fig. 9E), indicating that neural activity negatively regulates lola expression. The levels of GluRIIA mRNA and protein were also decreased upon the activity elevation (Fig. 9D,F). These results suggest that increased neural activity downregulates the transcription level of lola and downstream GluRIIA.
We used the highly accessible and well-characterized Drosophila NMJ to study gene expression changes in the postsynaptic cells upon presynaptic innervation. We performed comparative genome-wide gene expression profiling of muscles by individually picking up the cells at specific developmental stages (before and after innervation) and in pros mutants that lack nerve innervation. Through two series of comparative expression profiling (“before versus after” and “with versus without” innervation), we identified genes whose expression is potentially up- or downregulated upon innervation (72 and 12 genes, respectively). Subsequent systematic functional analyses of the candidate genes identified two genes, pst and lola, that play roles in synaptic development. Muscle-specific knock-down of pst decreased the number of boutons and branches. Since manipulating the expression of pst in muscles affected the growth of the presynaptic terminals, some retrograde signals derived from postsynaptic muscles are likely to be involved in causing the phenotypes. Previous studies describe several molecules that act as muscle-secreted retrograde signals (e.g., the BMP ortholog Gbb) and coordinate synaptic growth (McCabe et al., 2003; McCabe et al., 2004). pst may function in the molecular pathway(s) that regulate the activity of such retrograde signals.
We found that another candidate gene identified by the microarray analyses, lola, plays a prominent role in the regulation of GluR abundance. lola encodes a transcriptional factor with a BTB-Zn-finger motif. Previous studies on the Drosophila NMJ identified several molecular pathways that regulate GluR abundance through post-transcriptional mechanisms. Subsynaptic protein synthesis, mediated by aggregates of protein synthesis machinery, affects the level of postsynaptic GluRIIA (Sigrist et al., 2000; Sigrist et al., 2003). NF-κB/Dorsal, IκB/Cactus and IRAK/Pelle colocalize postsynaptically around GluR clusters and regulate postsynaptic GluR density through unknown post-transcriptional mechanism(s) (Cantera et al., 1999; Heckscher et al., 2007). A serine threonine kinase, PAK, which colocalizes with GluRs, regulates GluR abundance via its interaction with the adaptor protein Dreadlocks (Dock) (Albin and Davis, 2004; Parnas et al., 2001). Unlike these molecules, Lola regulates the abundance of GluRs at the level of mRNA. When the function of lola was knocked down in muscles, the mRNA levels of both of the specific subunits, GluRIIA and the GluRIIB, as well as the common subunit GluRIII were decreased. The abundance of GluR proteins at the synapses was also decreased.
In addition to the three GluR genes, lola also regulates the transcriptional level of another postsynaptic component, PAK. Thus, lola is a potential candidate gene that may control synaptic function by coordinating transcription of multiple post-synaptic components. A similar role has been proposed for Lola during axon guidance. During midline crossing of commissural axons, Lola regulates the expression of the guidance cue Slit and its receptor Robo (Crowner et al., 2002). We also observed that synaptic accumulation of Dlg protein but not the level of its mRNA is abnormal in lola mutants. This is likely a secondary consequence of the change in the expression level of other synaptic regulators whose mRNA level is directly controlled by Lola. A strong candidate is PAK since PAK is known to regulate synaptic abundance of Dlg (Albin and Davis, 2004). Thus, Lola appears to regulate the expression of synaptic components both directly and indirectly.
It remains to be determined how Lola controls the transcriptional level of GluRs. Lola is a transcription factor that can activate gene expression by binding to specific DNA sequences in cultured cells (Cavarec et al., 1997). Thus, Lola may directly modulate the transcription level of GluRs by binding to the regulatory region of the genes. Alternatively, Lola may indirectly modulate GluR levels by regulating the transcription and/or activity of other molecules. lola is a highly complex genetic locus and encodes at least 20 different splice isoforms (Goeke et al., 2003; Horiuchi et al., 2003). Most isoforms have four common exons that contain a BTB or POZ domain and distinct exons that encode one or a pair of zinc-finger domains. The lola sequence that was used in the present RNAi analyses targets all Lola isoforms. While some RNAi lines that target specific Lola isoforms are available, none of them targets the isoforms that are expressed in muscles. Isoform-specific conventional mutants are available for two Lola isoforms (lolaORE4 and lolaORE50 for isoform K, and lolaORE119 for isoform L) (Goeke et al., 2003). Of these, isoform L is expressed in muscle cells, so we tested whether lolaORE119 mutants have defects in GluR abundance but found no abnormality. It is thus currently difficult to determine isoform-specific functions of Lola at the NMJ.
The results of the microarray and qPCR analyses showed that the expression level of lola is downregulated by innervation. lola expression was downregulated after innervation and upregulated in pros mutants, which lack innervation. Similarly, increasing neural activity decreased the transcriptional level of lola. Thus, lola mRNA expression appears to be negatively regulated by innervation and neural activity. Since the expression of GluRs and other synaptic components are in turn regulated by lola, lola may link changes in neural activity to changes in postsynaptic gene expression. Consistent with this view, increased neural activity decreased the expression of the GluRIIA transcript and protein, concomitant with the decrease in lola expression. It is however currently unknown whether lola is solely responsible for the changes in GluRIIA expression in the activity misregulated NMJs. We also could not detect obvious changes in the expression of Lola protein in the aneural muscles (in pros mutants) or in activity-misregulated muscles (in dTRPA1 expressing larvae). Role of lola in coordinating activity-dependent changes in postsynaptic components remains to be explored in future studies.
How lola expression is regulated by neural innervation and/or activity also remains to be determined. Signaling mediated by PKA and CaMKII has been previously implicated in activity-dependent regulation of GluRIIA expression (Davis et al., 1998; Haghighi et al., 2003; Kazama et al., 2003; Morimoto-Tanifuji et al., 2004). Whether lola is regulated by these or other activity-dependent signaling mechanisms will be an interesting question to address in the future.
The negative regulation of lola by innervation and increased activity is somewhat counterintuitive in light of the fact that lola positively regulates the level of GluR and other postsynaptic components. The synaptic expression level of GluRs is known to increase after innervation (Broadie and Bate, 1993c; Chen and Featherstone, 2005). Some protocols that increase larval activity are known to induce an increase in GluRIIA expression (Schuster, 2006). Increasing motoneuron activity also increases their terminal size (Ataman et al., 2008; Budnik et al., 1990). One possibility is that Lola is part of the homeostatic mechanism(s) that maintain a proper level of synaptic strength. It has been shown that a decrease in postsynaptic function is compensated by an increase in presynaptic synaptic release at the Drosophila NMJs (Davis, 2006; Dickman and Davis, 2009; Frank et al., 2009). Similarly, Lola may control the appropriate level of postsynaptic components according to the level of presynaptic activity. For example, when presynaptic activity is too strong, expression of lola may be downregulated to reduce the level of GluRs at the postsynaptic cells. Conversely, when presynaptic innervation drops to an insufficient level in pros mutants, lola may be upregulated to compensate for the decreased transmission by increasing the level of GluRs.
Another possibility is involvement in the pruning of ectopic synapses. During the formation of neuromuscular connectivity in the embryos, reduction in electric activity in motor neurons results in an increase in ectopic synapses on non-target muscles (Jarecki and Keshishian, 1995; Carrillo et al., 2010). While involvement of presynaptic mechanisms has been shown by a recent study (Carrillo et al., 2010), postsynaptic Lola may also contribute to the regulation of ectopic synapses. For example, lola expression may be downregulated after innervation by proper motor neurons, to prevent the formation of ectopic synapses. Possible roles of lola in these processes will be determined in future studies.
Several studies suggest that Lola can regulate the chromatin structure and function in an epigenetic manner (Ferres-Marco et al., 2006; Zhang et al., 2003). Thus, Lola may be involved in the epigenetic control of synaptic plasticity, as shown in other systems (Guan et al., 2004; Levenson and Sweatt, 2005). Identification of Lola as a regulator of multiple postsynaptic components will provide a unique entry point to the understanding of transcriptional regulation of synapse formation and plasticity at the Drosophila NMJ.
We are grateful to E. Giniger, T. Aigaki, A. DiAntonio, M. Saitoe, A. Hofbauer, C. Hama, Bloomington and Kyoto Stock Center, Vienna Drosophila RNAi Center (VDRC), Fly Stocks of the National Institute of Genetics (NIG), The Exelixis Collection at the Harvard Medical School, and the Developmental Studies Hybridoma Bank for fly stocks and reagents. We acknowledge H. Miyakawa and members of the Nose lab for helpful discussion, H. Kohsaka for help in staining of the embryos, Y. Fujioka, H. Kato and M. Sotodate for technical assistance. This study was supported by a Grant-in-Aid for Scientific Research B and for Scientific Research on Priority Areas-Molecular Brain Science-from the MEXT to A.N.
Competing interests: We declare that there are no competing interests.