MAPK- and glycogen synthase kinase 3–mediated phosphorylation regulates the DEAD-box protein modulator Gle1 for control of stress granule dynamics

Rapid expression of critical stress response factors is a key survival strategy for diseased or stressed cells. During cell stress, translation is inhibited, and a pre-existing pool of cytoplasmic mRNA–protein complexes reversibly assembles into cytoplasmic stress granules (SGs). Gle1 is a conserved modulator of RNA-dependent DEAD-box proteins required for mRNA export, translation, and stress responses. Proper Gle1 function is critical as reflected by some human diseases such as developmental and neurodegenerative disorders and some cancers linked to gle1 mutations. However, the mechanism by which Gle1 controls SG formation is incompletely understood. Here, we show that human Gle1 is regulated by phosphorylation during heat shock stress. In HeLa cells, stress-induced Gle1 hyperphosphorylation was dynamic, primarily in the cytoplasmic pool, and followed changes in translation factors. MS analysis identified 14 phosphorylation sites in the Gle1A isoform, six of which clustered in an intrinsically disordered, low-complexity N-terminal region flanking the coil-coiled self-association domain. Of note, two mitogen-activated protein kinases (MAPKs), extracellular signal–regulated kinase (ERK) and c-Jun N-terminal kinase (JNK), phosphorylated the Gle1A N-terminal domain, priming it for phosphorylation by glycogen synthase kinase 3 (GSK3). A phosphomimetic gle1A6D variant (in which six putative Ser/Thr phosphorylation sites were substituted with Asp) perturbed self-association and inhibited DEAD-box helicase 3 (X-linked) (DDX3) ATPase activity. Expression of alanine-substituted, phosphodeficient GFP-gle1A6A promoted SG assembly, whereas GFP-gle1A6D enhanced SG disassembly. We propose that MAPKs and GSK3 phosphorylate Gle1A and thereby coordinate SG dynamics by altering DDX3 function.

A key survival pathway for cells impacted by disease states and environmental stresses is the rapid expression of critical stress response factors. To achieve this, nuclear export of nonessential transcripts is blocked, whereas temporally induced mRNA transcripts in protein complexes (mRNPs) 3 are selectively exported through nuclear pore complexes (NPCs) and then translated. Concurrently, general translation is inhibited, and a pre-existing pool of cytoplasmic mRNPs reversibly assembles into cytoplasmic stress granule (SG) complexes of RNA processing machinery, translation factors, and signaling proteins (1). In addition to the changes in gene expression profiles, the dynamics of mRNP occupancy in SGs affect rates of mRNA decay, global translational activity, and cell signaling cascades (2)(3)(4)(5). Increasing evidence suggests the involvement of altered SG assembly and disassembly in the pathogenesis of neurodegenerative diseases; however, the underlying means by which changes in SG dynamics impact diseases is not well understood (6 -9).
Several recent reports propose mechanisms for controlling the highly dynamic process of SG assembly and disassembly. SG components have liquid droplet properties and assemble through phase transitions in response to changing temperature, pH, or ionic strength (7,10,11). For some components, post-translational modifications (PTMs) mediate their SG assembly (12,13). For example, phosphorylation of the SG component G3BP by casein kinase 2 promotes dissociation of G3BP from SGs and triggers SG disassembly (14). Notably, aberrant activity of mitogen-activated kinases (MAPKs) is linked to amyotrophic lateral sclerosis (ALS), a neurodegenerative disease characterized by SG accumulation (15). Finally, SG dynamics are altered by the actions of RNA-dependent DEADbox proteins (Dbps in yeast/DDXs in humans) and their cofactors, possibly through remodeling of RNA-protein complexes or unwinding RNA-RNA duplexes (16).
The essential and highly conserved Dbp/DDX modulator Gle1 plays key roles in nuclear mRNA export, translation initiation and termination, and SG dynamics (17)(18)(19)(20)(21)(22)(23)(24). In each of these roles, Gle1 binds to and modulates the ATPase activity of distinct Dbps/DDXs bound to RNA. In humans, the GLE1 gene is alternatively spliced to generate at least two isoforms, Gle1A and Gle1B (25). Both human isoforms share high sequence similarity and common functional motifs: an amino (N)-terminal region that interacts with the human NPC component Nup155; a coiled-coil region that is involved in Gle1 self-association; a carboxyl (C)-terminal domain that interacts with DDX19B, DDX3, and inositol hexakisphosphate (IP 6 ); and a shuttling domain that mediates its translocation between the nucleus and cytoplasm (21,(25)(26)(27). Gle1B exhibits pancellular localization with pronounced steady-state enrichment at the NPCs that is partially dependent on a unique 39-amino acid C-terminal extension that mediates Nup42 binding (25,26). Our work further shows that Nup42 interaction and IP 6 are independently required to activate Dbp5/DDX19B for proper mRNA export (26). In contrast, Gle1A lacks the Nup42-binding domain, and it is not functional in mRNA export at NPCs (17). Instead, Gle1A localizes predominantly in the cytoplasm (25) where it interacts with DDX3 to modulate SG dynamics and translational repression in response to stress (17,28). Thus, Gle1A and Gle1B reside in distinct subcellular pools and perform nonoverlapping functions. These specificities provide human Gle1A and Gle1B with the capacity to regulate multiple steps of gene expression simultaneously, a critical aspect of the cellular stress response.
Several mutations in GLE1 are linked with human diseases, including ALS, multiple forms of arthrogryposis multiplex congenita, a broad array of developmental defects, and cancer (29 -33). Our prior studies of human gle1 variants linked to ALS and the arthrogryposis multiplex congenita disease lethal congenital contracture syndrome 1 (LCCS1) suggest that proper subcellular localization and the separable functions of each Gle1 isoform are central to both mechanisms of pathogenesis (17,21,28). Consistent with this premise is the observation that Gle1 subcellular localization is disrupted in mouse models of Huntington's disease, which exhibit nuclear localization of Gle1 in the brain cortex (34). Thus, controlled Gle1 subcellular localization and segregation of isoform-specific functions are critical for normal cell physiology.
Based on Gle1's intracellular dynamics and roles in stress responses, we speculated that mechanisms might exist to control Gle1 function in response to stress or disease. Here, we show that, under stress conditions, human Gle1 is hyperphosphorylated by MAPK and GSK3 in an N-terminal, low-complexity region. The basally phosphorylated pool of Gle1A promotes SG assembly and stimulates DDX3 activity, whereas Gle1A hyperphosphorylation promotes SG disassembly, inhibits DDX3 activity, and is disrupted in its capacity to oligomerize in vitro. We propose that the Gle1A hyperphosphorylation cycle effectively alters Gle1A self-association and DDX3 functions. These insights allow a further understanding of how perturbing Gle1 might contribute to different disease states and be a key connection between reported neurodegeneration and SG dysfunction.

Gle1 is phosphorylated during stress
We previously found a critical role for Gle1A in the SG response (17). However, the molecular events directing Gle1A activity during stress were not defined. To test whether PTMs might be involved, HeLa cell lysates from cells grown under stress conditions were analyzed by SDS-PAGE and anti-Gle1 immunoblotting. As shown in Fig. 1A, Gle1 in either heat shock-or arsenite-stressed cell lysates exhibited reduced electrophoretic mobility (lanes 4 and 7) compared with that of the untreated sample (lane 1). Similar shifts in mobility were detected for Gle1 in RPE-1 and HEK-293 cells following heat stress (data not shown). To test whether this slower migration was due to phosphorylation, cell lysates from untreated control, heat-shocked, or sodium arsenite-treated cells were incubated with -phosphatase (-PPase) in the presence or absence of phosphatase inhibitors. The -PPase treatment resulted in collapse of the slower migrating Gle1 band in lysates from stresstreated cells (Fig. 1A, lanes 5 and 8), whereas no change was observed in the presence of both -PPase and phosphatase inhibitors (Fig. 1A, lanes 6 and 9). Notably, a minor fraction of basally phosphorylated Gle1 was reproducibly observed in control cell lysates, whereas lysates from stress-treated cells contained a significantly greater proportion of phosphorylated Gle1 (Fig. 1A, compare lane 1 with lanes 4 and 7). To examine Gle1 phosphorylation by an alternative method, we utilized Phos-tag gel electrophoresis followed by immunoblot analysis (35). Three distinct bands for Gle1 were observed in heatshocked cell lysates (Fig. 1B, lane 2), indicating varying degrees of Gle1 phosphorylation. Interestingly, treatment with the broad-spectrum kinase inhibitor K252a prior to heat shock abolished the altered Gle1 mobility ( Fig. S1A and Fig. 1B, lane  3). Thus, upon stress treatment, phosphorylation of endogenous Gle1 was enhanced.
To test whether Gle1 phosphorylation occurred independently in the C-terminal (residues 400 -698) and/or N-terminal (residues 1-400) regions of the protein, plasmids expressing GFP-GLE1  or GFP-GLE1 400 -698 were transfected into HeLa cells. Of note (as shown in Fig. S1 and data not shown), heat shock treatment consistently reduced the steady-state level of all tested, exogenously expressed, GFP-tagged Gle1 proteins through a mechanism that remains to be determined. Analysis of respective cell lysates showed that only the GFP-Gle1 1-400 protein exhibited two distinct electrophoretic mobilities on a Phos-tag gel with an increased proportion of the slower migrating band present following stress (Fig. S1B, lanes  2 and 3). As observed for endogenous Gle1, K252a pretreatment decreased the amount of phosphorylated Gle1 present (Fig. S1B, compare lane 3 with lane 4). Because in silico phosphorylation analysis predicted a clustering of potential Gle1 phosphorylation sites within the first N-terminal 120 amino acids (data not shown), we also examined GFP-Gle1 120 -698 protein and found no electrophoretic shift in response to stress (Fig. S1B, lanes 10 and 11). Taken together, we concluded that residues in the N-terminal region of Gle1 were phosphorylated under these stress conditions.

Phosphorylation of Gle1A regulates DDX3 and stress response The Gle1 phosphorylation state correlates with preceding cycles of MAPK and translation factor phosphorylation
To investigate the temporal dynamics of Gle1 phosphorylation during stress and recovery, a time course analysis of heat shock-induced phosphorylation was conducted. Cell lysates were collected at time 0 and every 30 min throughout the 45°C heat stress and subsequent 37°C recovery phase. A dramatic increase in Gle1 phosphorylation was observed from the 60-min time point during heat stress to the 120-min time point in the recovery phase ( Fig. 1, C and D, lanes 1-4). This level of Gle1 phosphorylation persisted in the recovery phase for approximately 3 h (240 min experimental time point; lane 5) after which it slowly diminished and was undetectable 9 h poststress treatment (lane 8).
As both yeast Gle1 and human Gle1A function to regulate translation (17,19,20), the time course of Gle1 phosphorylation was further assessed by evaluating the concomitant phosphorylation state of various translation regulators (Fig. 1D). During stress, MAPKs and mechanistic targets of rapamycin (mTORs) are activated to repress global translation initiation and stimulate translation of stress response proteins (36). Changes in eIF2␣ and 4EBP1 phosphorylation states are well characterized and effectively halt translation initiation by, respectively, blocking formation of the ternary complex and eIF4F complex (37)(38)(39). As cells recover from stress, the phosphorylation states of these proteins are reversed, and global translation resumes. Immunoblotting was used to assess kinase activation (Fig. 1C) and to analyze translation factor phosphorylation (Fig. 1D) in cell lysates from the heat shock time course. As shown in Fig. 1C, the stress-induced phosphorylation cycles of ERK and JNK slightly preceded Gle1 phosphorylation and dephosphorylation (Fig. 1C). Strong phosphorylation of JNK was detected within 30 min (Fig. 1C, lane 2) of heat stress, whereas, consistent with other reports and quite reproducibly (40,41), a rapid early burst of ERK phosphorylation was observed at 30 min of stress treatment (lane 2) followed by a later peak of phosphorylation at 120 min (lane 4). JNK and ERK phosphorylation was reduced by the 240-min time point (lane 5) and continued to decrease throughout the monitored recovery period. Consistent with their roles in triggering the onset of translational reprogramming (42), the cycles of eIF2␣ phosphorylation and 4EBP dephosphorylation also preceded the timing of Gle1 phosphorylation changes (Fig. 1D). Overall, stress-induced Gle1 phosphorylation correlated with a response to Figure 1. Cytoplasmic Gle1 is dynamically phosphorylated in response to stress. A and B, Gle1 is phosphorylated during stress. A, HeLa cells were either left untreated, exposed to heat shock at 45°C, or treated with 0.5 mM sodium arsenite for 60 min. Cell lysates were prepared and incubated with either phosphatase buffer alone, -PPase, or -PPase and phosphatase inhibitors together for 30 min. Samples were resolved by SDS-PAGE and immunoblotted with anti-Gle1 antibodies. B, HeLa cells were either left untreated (UT) or preincubated with the broad-spectrum kinase inhibitor K252a (1 M) for 60 min at 37°C followed by heat shock at 45°C for 60 min. Cell lysates were resolved on a Phos-tag SDS-PAGE gel and immunoblotted with anti-Gle1 antibodies. C and D, Gle1 phosphorylation follows MAPK activation (C) and phosphorylation (p) of translation factors (D). HeLa cells subjected to heat shock at 45°C for 60 min followed by a recovery phase at 37°C were harvested at the indicated time points. Lysates were analyzed by immunoblotting using the indicated antibodies. E, phosphorylated Gle1 localizes to the cytoplasm. Following heat shock at 45°C for 60 min, nuclear (N) and cytosolic (C) HeLa cell fractions were prepared, resolved by SDS-PAGE, and immunoblotted with anti-Gle1, anti-lamin B1, and anti-GAPDH antibodies. T, total cell lysate.

Phosphorylation of Gle1A regulates DDX3 and stress response
MAPK activation and translation repression cues, and Gle1 dephosphorylation correlated with MAPK deactivation and translation activation.
Based on the link to translation factor changes, we reasoned that cytoplasmic Gle1 might be the predominant target for phosphorylation compared with the nuclear pool. To test this, nuclear and cytosolic fractions were prepared from untreated and heat-shocked HeLa cells, and Gle1 phosphorylation was analyzed by immunoblotting (Fig. 1E). Lamin B1 and GAPDH were controls for nuclear and cytosolic fractions, respectively. The slower migrating, phosphorylated Gle1 was observed in the cytosolic fraction for both untreated and stressed cells (Fig. 1E, lanes 2 and 5). In contrast, Gle1 detected in the nuclear fraction migrated at the molecular mass of nonphosphorylated Gle1 (Fig. 1E, lanes 3 and 6). Thus, in both untreated and heat shock-stressed cells, phosphorylated Gle1 was predominantly cytoplasmic. The presence of phosphorylated Gle1 in the cytoplasm may reflect the targeting of only cytoplasmic Gle1 by kinases or, alternatively, rapid export of Gle1 out of the nucleus following phosphorylation.

Gle1A is phosphorylated at periodically repeating serine/threonine residues in the intrinsically disordered N-terminal region
Studies suggest that protein phosphorylation often occurs within intrinsically disordered regions (43,44). Based on our analysis ( Fig. 2A) and as recently reported (45), the first 120 residues of the Gle1 N-terminal region are predicted to comprise an intrinsically disordered, low-complexity domain. Together with the domain mapping analysis in Fig. S1B, we A, the N-terminal region of Gle1A is predicted to be a highly disordered, low-complexity domain. Disorder potential was predicted using the IUPred web server. Scores above 0.5 indicate a high degree of predicted highly disordered. B, schematic representation of phosphorylated Gle1A residues identified by MS. The diagram depicts the functional and structural domains of Gle1A (adapted from Folkmann et al. (21)) and location of phosphorylated residues therein. Gle1A residues shown in gray were phosphorylated in both untreated and stress-treated samples. Residues detected only in stress-treated samples are shown in black. Phosphorylation of Ser 41 was detected only in untreated samples. The asterisk indicates the position of the PFQ insertion generated by the gle1-Fin major mutation. C, a periodicity is evident in the pattern of phosphorylation events across the N-terminal cluster of phosphorylated amino acids. Dotted arrows across two consecutive SXXXS(P) motifs in green show the in silico prediction of GSK3 priming by MAPK phosphorylation. Purple residues illustrate the PXXP polyproline motif of amino acids 97-103. D, modification of the Ser 88 -Thr 102 cluster alters Gle1A electrophoretic mobility. HeLa cells transfected with plasmids expressing GFP-GLE1A, GFP-gle1A 6A , or GFP-gle1A 6D were either untreated or exposed to 0.5 mM sodium arsenite for 60 min. Cell lysates were resolved by Phos-tag SDS-PAGE and immunoblotted with anti-GFP antibodies. E, phosphorylation alters the b-isox precipitation of Gle1. Nuclear (N) and cytosolic (C) fractions prepared from HeLa cells exposed to 45°C heat shock for 60 min were incubated with 100 M b-isox, pelleted by centrifugation, and analyzed by immunoblotting with anti-Gle1 and anti-actin antibodies. Sup, supernatant.

Phosphorylation of Gle1A regulates DDX3 and stress response
speculated that the Gle1 residues modified by stress-induced phosphorylation reside in this predicted low-complexity domain. Given that we observed stress-induced phosphorylation specifically in the cytoplasmic fraction of endogenous Gle1 (Fig. 1E) and the predominantly cytoplasmic isoform Gle1A, but not Gle1B, is specifically required during stress to modulate translation (17), we selected Gle1A for a proteomics approach to identify the phosphorylated amino acid residues. HeLa cells ectopically expressing GFP-GLE1A were either left untreated or treated with 0.5 mM sodium arsenite for 60 min. GFP-Gle1A protein was then isolated by immunoprecipitation, excised from an SDS-PAGE gel (Fig. S2), trypsin-digested, and processed for MS. In total, MS/MS analysis identified 14 Gle1A phosphorylation sites in samples subjected to stress. The locations of modified residues are shown in the domain map of Fig.  2B. Phosphorylation of five residues (Ser 88 , Ser 93 , Ser 99 , Thr 427 , and Ser 465 ) was also detected in Gle1A isolated from untreated cells, consistent with the observed basal Gle1 phosphorylation (Fig. 2B, denoted in gray). The majority of stress-induced phosphorylation sites concentrated in two small clusters adjacent to either the N-or C-terminal boundary of the coiled-coil domain. With regard to the hyperphosphorylated N-terminal region (Ser 88 , Ser 92 , Ser 93 , Ser 99 , and Thr 102 ), we validated the modifications at residues Ser 88 , Ser 99 , and Thr 102 by manual inspection of the MS/MS scans (Fig. S3). Although the MS/MS scans obtained for the Ser 92 and Ser 93 sites were less conclusive, evidence for phosphorylation of all five residues in the N-terminal cluster as well as for Ser 96 is also reported in published large-scale data sets from high-throughput phosphoproteomic screens (46 -48).
The Ser 88 /Ser 92 /Ser 93 /Ser 96 /Ser 99 /Thr 102 cluster of phosphorylated Gle1A residues (hereafter termed the Ser 88 -Thr 102 cluster) lies within the N-terminal 120-amino acid region of Gle1A responsible for the phosphatase-sensitive mobility shift according to our domain mapping analysis (Fig. S1B). Although our MS data indicated that three of these residues (Ser 88 , Ser 93 , and Ser 99 ) were also phosphorylated in Gle1A isolated from untreated cells, a quantitative analysis of individual phosphorylation levels was not possible because the multiply phosphorylated peptides eluted together in our LC-MS experiments. Therefore, we could not definitively exclude increased phosphorylation at these residues as a possible contributor to this change in mass. Conversely, the C-terminal region of Gle1 did not demonstrate a mobility shift under stress conditions (Fig.  S1B), suggesting that the degree of stress-induced modification at the seven identified C-terminal phosphorylation sites is less significant. These data therefore singled out the N-terminal Ser 88 -Thr 102 cluster as a putative region of hyperphosphorylation during stress.
Remarkably, a striking periodicity was evident in the potential pattern of phosphorylation events across this N-terminal cluster, including modifications occurring at every third residue from Ser 93 to Thr 102 and within two sequential SXXXS(P) motifs (where S(P) represents phosphorylation of serine) (Fig.  2C). Two polyproline II helix peptide motifs (PXXP), which when phosphorylated could markedly alter Gle1 structure and interactions, are also present in this sequence adjacent to Ser 96 and containing Ser 99 and Thr 102 (PXXP; Fig. 2C) (49 -51). Phos-phoNet in silico kinase prediction analysis across this sequence suggested that Gle1 residues Ser 88 , Ser 99 , and Thr 102 are potential targets for phosphorylation by ERK, and that Ser 93 and Ser 96 are targeted by JNK. Furthermore, when primed by Ser 96 phosphorylation, GSK3 was predicted to sequentially phosphorylate Ser 92 and Ser 88 (Fig. 2C). Taken together, the Gle1 phosphorylation domain mapping and in silico analysis strongly indicated that Ser 88 -Thr 102 is phosphorylated en masse following the stress-induced activation of MAPK/ERK1/2 and Akt/GSK3 cascades with potentially important functional repercussions.
Given the preponderance of data implicating the entire Ser 88 -Thr 102 cluster as a likely target of massive stress-induced MAPK and GSK3 phosphorylation, we investigated the impact that global modification of these residues has on Gle1 function. To test whether modifications in the Ser 88 -Thr 102 cluster were responsible for the mobility shift of Gle1 from cells grown under stress conditions, mutations in the coding sequences were generated that result in nonphosphorylatable alanines (gle1A 6A ) or phosphomimetic aspartic acids (gle1A 6D ) in place of the six putative serine or threonine phosphorylation sites. Plasmids expressing GFP-GLE1A, GFP-gle1A 6A , or GFP-gle1A 6D were transfected into HeLa cells, and lysates from untreated or arsenite-treated samples were analyzed by Phostag electrophoresis and anti-GFP immunoblotting. Arsenite stress treatment resulted in reduced electrophoretic mobility of GFP-Gle1A compared with untreated protein (Fig. 2D, compare lane 4 with lane 1), whereas no mobility shift was observed for GFP-gle1A 6A (lanes 2 and 5). GFP-gle1A 6D exhibited slower mobility irrespective of treatment (Fig. 2E, lanes 3 and 6) as would be predicted for a phosphomimetic form. Thus, modification of residues in the Ser 88 -Thr 102 cluster was required to produce the gel mobility shift in response to stress.

Phosphorylation of Gle1 alters its biochemical properties
Others have reported that phosphorylation events within low-complexity domains alter the respective protein's polymerization characteristics and function (43,44,52,53). Because low-complexity domains can mediate the type of phase transitions that potentially drive SG dynamics (11), phosphorylation of SG components may serve to regulate that process. Others also have demonstrated that the isoxazole (b-isox) compound specifically interacts with highly disordered, low-complexity domains and results in protein precipitation. Indeed, many SG components precipitate with b-isox (52,54), and a recent study shows that phosphorylation of Cdc19's low-complexity domain abolishes its precipitation by b-isox (53,55). Thus, we established assays to measure the propensity for Gle1 precipitation by b-isox as a test of the effect of phosphorylation on its biochemical properties. The nuclear and cytoplasmic fractions isolated from heat-shocked HeLa cells were treated with b-isox, and the resulting pelleted and soluble portions were tested for Gle1 precipitation by immunoblotting (Fig. 2E). Consistent with our data in Fig. 1E, phosphorylated Gle1 was in the total cytoplasmic fraction (Fig. 2E, lane 5). In the presence of b-isox, the phosphodeficient Gle1 in the nuclear fraction was distributed between the soluble and pelleted fractions (Fig. 2E, lanes 1  and 3). However, phosphorylated Gle1 was only detected in the supernatant following b-isox incubation (Fig. 2E, lane 4). The Phosphorylation of Gle1A regulates DDX3 and stress response same distribution was observed in untreated cells with basally phosphorylated Gle1 remaining soluble when incubated with b-isox (data not shown). These data indicated that the N-terminal low-complexity domain of phosphodeficient Gle1 is highly disordered and prone to aggregation like that of many SG components. Conversely, Gle1 hyperphosphorylation appeared to change the protein's biochemical properties, reducing its b-isox precipitation propensity. Given the high degree of disorder inherently surrounding many of the Gle1 phosphorylation sites, it is likely that the altered behavior arises from charge-based disruption of intramolecular bonds or oligomerization interactions. This effect was previously described for FUS, which is also phosphorylated in a low-complexity domain. Modification of FUS does not alter the disordered structure of the domain, but rather the increased negative charge introduced by multiple phosphorylation events (or the phosphomimetics) decreases intramolecular and intermolecular interactions and thereby disrupts FUS aggregation (56).

MAPKs prime Gle1 for phosphorylation by GSK3
ERK and JNK belong to the MAPK family and play critical roles in stress responses (57,58). Translation inhibitors such as cycloheximide and anisomycin activate MAPKs (59) and prevent polysome disassembly, trapping the ribosomes in different states (60) and blocking SG formation to varying degrees as shown by anti-G3BP immunofluorescence (Fig. S4A). In agreement with the prediction of MAPK sequence motifs in the Ser 88 -Thr 102 cluster, treatment of HeLa cells with either cycloheximide or anisomycin (in the absence or presence of stress) dramatically enhanced the phosphorylation of Gle1 (Fig. S4B). Thus, increased Gle1 phosphorylation was potentially not dependent on SG assembly but might be linked to the regulation of translation.
To test whether MAPKs are required for Gle1 phosphorylation, specific kinase inhibitors of the three major MAPK subgroups (ERK1/2, JNK, and p38) were tested for inhibition of Gle1 phosphorylation in HeLa cells exposed to heat shock. The resulting cell lysates were analyzed by anti-Gle1 immunoblotting. Addition of p38 inhibitor had no impact upon Gle1 phosphorylation (Fig. 3A, lanes 3 and 7). Similarly, inhibition of either JNK or ERK alone did not alter the relative degree of Gle1 phosphorylation (Fig. 3A, lanes 4, 5, and 6) compared to control with no inhibitors (Fig. 3A, lane 2). However, in the presence of inhibitors for both JNK and ERK, a faster migrating Gle1 band indicative of reduced phosphorylation was observed (Fig. 3A,  lane 8). The intermediate mobility of this band might indicate less than 100% inhibition of ERK and JNK or represent the actions of another kinase on Gle1. This latter possibility was supported by the greater efficacy of the MEK1/2-targeted ERK inhibitor U0126 than the direct ERK inhibitor FR180204 at preventing Gle1 phosphorylation when combined with JNK inhibitor (Fig. 3A, lanes 8 and 9, respectively). In sum, these results suggested that actions of MAPKs, including the ERK and JNK, are involved in Gle1 phosphorylation.
To further investigate whether ERK and JNK can directly phosphorylate Gle1, recombinant purified N-terminal Gle1 1-360 or gle1 1-360 6A was incubated with or without purified active ERK and/or JNK in the presence of [␥-32 P]ATP. In the absence of exogenous substrate, autophosphorylation of JNK was observed (Fig. 3B, lane 2). Incubation of Gle1 1-360 with ERK and/or JNK resulted in phosphorylation of the N-terminal Gle1 protein (Fig. 3B, lanes 4, 5, and 6). Phosphorylation of gle1 1-360 6A was minimal compared with WT ( Fig. 3B; lanes 8, 9, and 10), and phosphorylation of the C-terminal Gle1 360 -698 protein was also negligible (Fig. S5). Together, these data demonstrated that ERK and JNK specifically phosphorylate Gle1 and that phosphorylation requires the N-terminal Ser 88 -Thr 102 cluster of residues.
Based on the sequence motifs surrounding the Ser 88 -Thr 102 cluster, we also tested whether Gle1 is phosphorylated by GSK3. GSK3 typically targets substrates primed by prior phosphorylation events at serine or threonine residues located four residues C-terminal to the GSK3 phosphorylation site (61). We reasoned that Ser 92 might be a GSK3 target site that is primed by ERK or JNK phosphorylation at Ser 96 (Fig. 2C). To determine whether the MAPKs prime Gle1 for GSK3 phosphorylation, we performed sequential in vitro kinase assays. Purified Gle1 1-360 was preincubated with ERK and JNK followed by the addition of active recombinant GSK3 in the presence of [␥-32 P]ATP (Fig. 3C, schematic). A minor band corresponding to GSK3 autophosphorylation was detected in the absence of any Gle1 protein (Fig. 3D, lane 15). As expected, incubation with GSK3 alone resulted in minimal phosphorylation of Gle1 1-360 (Fig. 3D, lane 4). However, when preincubated with ERK or JNK, robust phosphorylation of Gle1 1-360 by GSK3 was detected (Fig. 3D, lanes 5-7). To investigate whether this phosphorylation signal arises solely from primed GSK3 activity and not from residual MAPK activity, ERK and/or JNK specific inhibitors were added to the reactions after the respective MAPK preincubation step. In the absence of GSK3, addition of the MAPK inhibitors after pretreatment, but prior to incubation with [␥-32 P]ATP, completely abolished the incorporation of radiolabeled 32 P into Gle1 (Fig. 3D, lanes 8-10). In contrast, if GSK3 was included after addition of the MAPK inhibitors, dramatic Gle1 phosphorylation was observed (Fig. 3, D, lanes [12][13][14], and E, lanes [3][4][5]. Consistent with its lack of phosphorylation by ERK and JNK, GSK3 also did not phosphorylate gle1 1-360 6A (Fig. 3E, lanes 8 -10). Together, these results strongly implicated a kinase cascade of ERK/JNK priming events followed by GSK3 phosphorylation in regulating Gle1A function during stress.

The Gle1 phosphorylation cycle directs its role in SG assembly and disassembly
Given that Gle1A is required for proper SG assembly (17) and stress-induced Gle1A phosphorylation occurred concurrently with changes in translation factors (Fig. 1D), we speculated that Gle1A phosphorylation might play a role during SG assembly and/or disassembly. In ectopic expression experiments, both GFP-gle1A 6A and GFP-gle1A 6D localized to SGs upon heat shock (Fig. S6). To assay for SG perturbations resulting from specific Gle1 phosphorylation states, we performed knockdown/addback experiments using our previously established methods (17). Non-targeting (CTRL) or GLE1 siRNA-treated HeLa cells expressing either GFP, GFP-GLE1A, GFP-gle1A 6A , or GFP-gle1A 6D were heat-shocked for 60 min and processed Phosphorylation of Gle1A regulates DDX3 and stress response for anti-G3BP immunofluorescence analysis (Fig. 4A). In CTRL siRNA-treated cells after heat shock for 60 min (Fig. 4A), the expression of phosphomimetic GFP-gle1A 6D behaved similarly to that of GFP alone with no significant difference in the number of SGs per cell. In contrast, GFP-GLE1A and GFP-gle1A 6A expression resulted in substantially fewer SGs per cell in the heat-shocked, CTRL siRNA-treated samples (Fig. 4A). Thus, as we previously reported for GFP-GLE1A (17,28), overexpression of GFP-gle1A 6A also altered SG morphology. Consistent with our previous reports, the number of SGs also significantly increased in GLE1 siRNA-treated cells expressing GFP compared with CTRL siRNA-treated cells expressing GFP (p value Ͻ0.0001), and this phenotype was rescued by expressing either GFP-GLE1A, GFP-gle1A 6A , or GFP-gle1A 6D (Fig. 4A). Interest-ingly, in the heat-shocked GLE1 siRNA-treated cells, GFP-gle1A 6A expression resulted in greater rescue of the SG defect than GFP-GLE1A or GFP-gle1A 6D (Fig. 4A). This analysis indicated that Gle1A-mediated SG assembly is not dependent on phosphorylation of the N-terminal Ser 88 -Thr 102 cluster.
Because Gle1 dephosphorylation coincided with the release of translation factors from inhibitory phosphorylation states (Fig. 1D), we speculated that Gle1 dephosphorylation might be required for SG disassembly. As a measure of SG disassembly, the number of cells with no anti-G3BP-labeled SGs was quantified for each knockdown/addback condition wherein, after heat shock for 60 min, duplicate cell samples were allowed to recover at 37°C for 60 min (Fig. 4B). Compared with GFP alone, an increase in the number of cells without SGs would indicate

Phosphorylation of Gle1A regulates DDX3 and stress response
enhanced SG disassembly during the 60-min recovery period. As reported previously (28), GFP-GLE1A expression, in CTRL and GLE1 siRNA-treated cells, resulted in a greater percentage of cells with no SGs compared with GFP alone (Fig. 4B). Thus, GFP-Gle1A promoted SG disassembly. Expression of GFP-gle1A 6D in CTRL and GLE1 siRNA-treated cells resulted in similarly high levels of cells lacking SGs within 60 min of recovery from heat shock as with GFP-GLE1A expression. This suggested that GFP-gle1A 6D mediates SG disassembly to the same extent as WT GFP-Gle1A within a 60-min recovery phase. However, with GFP-gle1A 6A expression when endogenous Gle1 is depleted (GLE1 siRNA cells), fewer cells were without SGs compared with GFP-GLE1A or GFP-gle1A 6D (Fig. 4B). Thus, GFP-gle1A 6A alone was not as effective as WT or GFP-gle1A 6D in mediating SG disassembly during the 60-min recovery period. Collectively, the impact of Gle1 on SG dynamics was distinctly phosphodependent. These results suggested that, during stress, an unphosphorylated or basally phosphorylated pool of Gle1A promotes SG assembly. Then, during recovery from stress, increasing Gle1A hyperphosphorylation in the Ser 88 -Thr 102 cluster facilitates SG disassembly.

Phosphorylation of the Ser 88 -Thr 102 cluster in Gle1A inhibits DDX3 ATPase activity
Our prior studies in human cells linked the role of Gle1A function in SG dynamics to the regulation of DDX3 (17). DDX3 modulates translation through interaction with the preinitiation complex and is required for proper SG assembly (60). Because Gle1A modification at the Ser 88 -Thr 102 cluster promoted SG disassembly (Fig. 4B), we tested for the impact of Gle1A phosphorylation on DDX3 ATPase activity. In vitro ATPase assays were conducted with bacterially expressed, purified recombinant DDX3 and Gle1A, gle1 6A , or gle1A 6D at various molar ratios. As expected, DDX3 exhibited RNA-dependent ATPase activity (Fig. 5, bar 3) (62), which was further stimulated by addition of Gle1A (Fig. 5, bars 7 and 8). Addition of gle1A 6A also stimulated DDX3 activity to a level equivalent to WT Gle1A (Fig. 5, bars 9 and 10). However, addition of gle1A 6D inhibited DDX3 ATPase activity in a dose-dependent manner (Fig. 5, bars 11-13). IP 6 did not appear to play a role in Gle1 stimulation of DDX3 as its addition to the assay had no effect (Fig. 5, bars 8, 10, and 14). We concluded that the phosphorylation state of Gle1 alters its ability to stimulate DDX3. From these data, we further speculated that unphosphorylated or basally phosphorylated Gle1 promotes SG assembly by activat- Figure 4. Gle1A phosphorylation state alters SG assembly and disassembly. A, gle1A 6A promotes SG assembly. CTRL or GLE1 siRNA-treated HeLa cells transfected with GFP, GFP-GLE1A, GFP-gle1A 6A , or GFP-gle1A 6D and heat-shocked at 45°C for 60 min were processed for anti-G3BP immunofluorescence and analyzed for the number of SGs present. B, Gle1A phosphorylation promotes disassembly of SG granules. CTRL or GLE1 siRNA-treated cells transfected with GFP, GFP-GLE1A, GFP-gle1A 6A , or GFP-gle1A 6D were heat-shocked at 45°C for 60 min and processed for anti-G3BP immunofluorescence before and after recovery at 37°C for 60 min, and the number of cells without SGs postrecovery was quantified. Analysis was performed using a two-tailed unpaired t test, and error bars represent S.E.

Phosphorylation of Gle1A regulates DDX3 and stress response
ing DDX3. Likewise, once the stress is removed, inhibition of DDX3 activity by hyperphosphorylated Gle1 would reverse the process and support SG disassembly.

Gle1A is a stable component of SGs
The SG assembly factors typically designated as nucleators, such as G3BP and TIA-1, are largely mobile within SGs and exhibit rapid exchange rates of 2-3 s (63,64). In contrast, RNAbinding proteins such as YB1 and HuR (63) and the Fas-activated serine/threonine phosphoprotein (FAST) (64) are stable SG components. Although Gle1A function and phosphorylation state modulate several aspects of the stress response, including SG assembly, disassembly, and DDX3 activity, it was unclear whether Gle1A is transiently localized at SGs like a nucleator or whether it is a stable SG component. To address this question, the exchange rate of Gle1A between SGs and the cytoplasm was examined under heat shock conditions. HeLa cells ectopically expressing GFP-GLE1A and mCHERRY-G3BP were subjected to heat shock at 45°C for 45 min, and then fluorescence recovery after photobleaching (FRAP) was measured only for SGs that had colocalized mCherry-G3BP and GFP-Gle1A. Large GFP-Gle1 aggregates devoid of mCherry-G3BP were excluded from analysis. SGs were photobleached using a 561 nm laser for mCherry-G3BP FRAP (Fig. 6A, top  panel) or using a 488 nm laser for GFP-Gle1A FRAP (Fig. 6A,  bottom panel). Consistent with published studies (63, 64), mCherry-G3BP moved rapidly in and out of SGs (mobile fraction of 63.8%) and exhibited a fast half-time recovery of 4.4 s. In contrast, GFP-Gle1A was relatively immobile in stress granules (mobile fraction of 6.1%). No significant difference in mobility was observed for the phosphodeficient GFP-Gle1A 6A (mobile fraction of 7.8%) and phosphomimetic GFP-Gle1A 6D (mobile fraction of 10.2%) compared with GFP-Gle1A (Fig. S7). These results suggested that GFP-Gle1 is a stable component of SGs.

Gle1 self-association is required for SG assembly and perturbed in vitro with gle1A 6D
The full N-terminal domain of Gle1 from residues 1 to 360, containing both the intrinsically disordered, low-complexity region and a coiled-coil domain, self-associates in vitro to form higher-order disk-like oligomer structures (21). Because the two major clusters of phosphorylated Gle1 residues map to regions flanking the coiled-coil domain, we investigated whether phosphorylation altered the disk-like, oligomerized Gle1 structures. High-molecular-mass complexes of purified, recombinant Gle1 1-360 , gle1 1-360 6A , and gle1 1-360 6D were isolated and analyzed by negative-stain EM. As expected, Gle1 1-360 formed organized disk-like structures (Fig. 7A). The gle1 1-360 6A protein also formed higher-order symmetrical disks that were indistinguishable from Gle1 1-360 (Fig. 7B). However, gle1 1-360 6D assembled into less ordered, more elliptical disks (Fig. 7C). Quantification of the disk circularity, as shown by the outline of disk shape and major/minor axis cross-hairs depicted in Fig. 7, A-C, insets, documented that gle1 1-360 6D disks were more elongated compared with Gle1 1-360 or gle1 1-360 6A disks (Fig. 7D). Thus, the phosphomimetic variant slightly altered Gle1 self-association in vitro, further suggesting that Gle1 hyperphosphorylation might perturb Gle1 self-association in a cellular context.
To evaluate whether proper self-association of Gle1 plays a role in SG assembly, we also examined the LCCS1-causative gle1-Fin major mutation that perturbs Gle1 self-association and leads to a defect in mRNA export (21). The gle1-Fin major mutation results in a PFQ insertion in the coiled-coil domain of Gle1 (31), and in vitro analysis shows gle1-Fin major disk-like structures are deformed (21), although in a different manner than gle1 1-360 6D disks. GFP-gle1A-Fin major expression in heatshocked GLE1 siRNA cells rescued the SG defects, however not as robustly as compared with GFP-GLE1A (Fig. 7, E and F). The extent of SG rescue by GFP-gle1A-Fin major expression was similar to that observed for GFP-gle1A 6D (Fig. 4B). Thus, proper self-association was required for effective SG assembly. Moreover, the N-terminal domain of Gle1 was specifically implicated in regulation of SG biology, as either the PFQ insertion in the coiled-coil domain or phosphomimetic alterations in the intrinsically disordered, low-complexity region impacted SG dynamics.

Discussion
Although phosphorylation is a highly temporal and rapidly reversible PTM found on many SG regulatory proteins, it has not been fully resolved how such stress-induced phosphorylation events might control SG dynamics. Here, we demonstrate that human Gle1 is hyperphosphorylated in response to cellular

Phosphorylation of Gle1A regulates DDX3 and stress response
stress, with this phosphorylation persisting until modifications that alter translation factors are reversed. As represented in the cartoon model (Fig. 8), we show that Gle1A residues in the N-terminal Ser 88 -Thr 102 cluster are substrates for sequential phosphorylation by the MAPKs ERK and JNK and by GSK3. Furthermore, the effect of mutations that change the coding sequence to that for alanine or aspartic acid at these residues suggests that the phosphorylation status of the Ser 88 -Thr 102 cluster impacts SG assembly/disassembly dynamics as well as DDX3 activity. In the early stages of stress response, we propose that Gle1A in its basally modified state self-associates via its N-terminal domain and promotes SG assembly. We further posit that the accumulation of phosphorylation events in the Ser 88 -Thr 102 cluster over the course of the stress response changes the low-complexity character of Gle1A's N terminus, altering Gle1A self-association and thereby promoting SG disassembly. As such, the stress-induced phosphorylation of Gle1A and subsequent dephosphorylation during recovery provide a critical layer of control over the gene expression pathway and cell survival. It remains to be determined whether each of the six phosphosites in the Ser 88 -Thr 102 cluster or specific combinations thereof are necessary for these functional changes. This will be a focus of our future studies.
We found that a sequential cascade of stress-induced MAPK and GSK3 activity underlies Gle1 phosphorylation. During Insets depict disk outline (dashed lines) with major and minor axes indicated (solid lines). Scale bar, 50 nm. D, circularity was quantified from at least 171 disks per protein and analyzed by unpaired t test. Major axis/minor axis represents the ratio of longest to shortest diameter of the elliptical as shown by insets in A-C. A ratio of 1 indicates a circle. Error bars represent S.D. E and F, GFP-gle1A-Fin major exhibits reduced capacity for rescuing SG defects. E, CTRL or GLE1 siRNA-treated cells were transfected with GFP, GFP-GLE1A, or GFP-gle1A-Fin major plasmids, heat-shocked for 60 min at 45°C, and processed for immunofluorescence using anti-G3BP antibodies. Scale bar, 10 m. F, quantification from three independent experiments, performed in duplicate, is shown for the number of SGs in CTRL or GLE1 siRNA-treated cells expressing the indicated plasmids. Data were quantified for at least 65 cells per condition and analyzed using a two-tailed, unpaired t test. Error bars represent S.E.

Phosphorylation of Gle1A regulates DDX3 and stress response
stress conditions, multiple signaling pathways are triggered to mediate a multifaceted stress response. MAPKs and GSK3 are implicated in proper SG function and disease. For instance, aberrant MAPK p38 activity induces a gain of function for the ALS-linked SG component FUS (65), and, conversely, inhibition of MAPK signaling improves survival of SOD1 E100G mutant motor neurons derived from ALS patient induced pluripotent stem cells (66). Likewise, tight regulation of GSK3 activity appears to underlie proper SG responses, given that GSK3 inhibition dramatically reduces stress granule formation (12,67), and aberrant activation of GSK3 is associated with ALS (68). Moreover, several reports show reductions in ALS phenotypes using GSK3 inhibitors, positioning it as a potential therapeutic target (69,70). Given our prior studies of ALS-linked gle1 mutations that alter Gle1 function (28,30), phosphoregulation of Gle1 by both MAPK and GSK3 phosphorylation offers insights into new mechanisms by which stress responses in ALS cells might be modulated.
A model for how Gle1 phosphorylation regulates SG biology is based on several pieces of evidence. We previously reported that Gle1A regulates SGs by modulating DDX3's role in translation repression. In brief, heat shock stress of Gle1-depleted HeLa cells results in defective assembly of mature SGs and loss of translation inhibition, both of which are rescued by DDX3 overexpression (17). Here, we show mutagenesis data indicating that Gle1A phosphorylation states modulate the ATPase activity of DDX3 and regulate SG assembly and disassembly. The phosphodeficient gle1A 6A stimulates DDX3 ATPase activity and rescues SG assembly defects in Gle1-depleted cells, arguing that Gle1A stimulates SG assembly in the absence of stress-induced hyperphosphorylation. In contrast, the phosphomimetic gle1A 6D inhibits DDX3 ATPase activity and promotes SG disassembly. This suggests that Gle1A drives DDX3 ATPase activity to stimulate SG assembly. Although previous studies have argued that DDX3 does not require its ATPase activity for translational repression (62), acting instead by displacement of eIF4E in the cap-binding complex (71,72), interpretation of this prior work is complicated by experimental conditions based on overexpression. Furthermore, another study recently reported conflicting data, showing that a medulloblastoma-linked ATPase-dead variant of DDX3 is unable to repress translation following arsenite stress (73). Taken together with our findings, it is plausible that Gle1A regulation of DDX3 catalytic activity does play a role in translation inhibition and SG assembly, perhaps via modulation of RNA structure and proteins associated with translation machinery or by mRNP remodeling in SGs (models are reviewed in Ref. 74). Similarly, eIF4B and eIF4G modulate the ATPase activity of eIF4A, which is required for SG assembly (75,76). Moreover, this mode of regulating ribonucleoprotein bodies appears to be conserved, since recent work in budding yeast showed that the Not1 cofactor stimulates the ATPase activity of Dhh1, which controls P body assembly and disassembly (77). Our finding suggesting that Gle1A hyperphosphorylation alters DDX3 activity and SG dynamics adds an additional layer of complexity to this model. Future studies will be necessary to provide a better understanding of how Gle1A modulates DDX3 functions in vivo during SG assembly and translation inhibition.
If SG assembly is driven in large part by phase separation as proposed in current models (78), then Gle1A phosphorylation may directly alter this physical property and thereby diminish SG formation. There are numerous examples of SG components that, when phosphorylated, have altered SG assembly properties (10,11,36,54,79). For example, G3BP is subject to phosphorylation in its intrinsically disordered region, which abrogates its function as an SG nucleator (13). Likewise, the ALS-linked SG nucleator FUS is hyperphosphorylated in its low-complexity domain by DNA-dependent protein kinase following DNA damage, which reduces its phase separation properties and aggregation (56). Both examples mirror the effects of phosphorylation on Gle1A, wherein hyperphosphorylation of its intrinsically disordered N-terminal domain changes its functioninSGassemblycomparedwiththenonhyperphosphorylated form. Our data suggest that this does not result from a loss of Gle1A recruitment to SGs when it is hyperphosphorylated because both GFP-gle1A 6A and GFP-gle1A 6D localize to stress granules (Fig. S6). Indeed, recruitment to SGs does not appear to be a factor in dictating the effect of Gle1 on SGs, as Gle1B also localizes to SGs but is incapable of supporting SG dynamics (17).
Although the mechanism by which Gle1A's phosphorylation state dictates DDX3 activity is not resolved, several pieces of data point to Gle1 self-association as a potential switch. We know that Gle1 self-associates via its N-terminal 360 amino acids and that disruption of Gle1B self-association by the LCCS1-linked gle1-Fin major mutation (annotated in Fig. 2) impedes its ability to stimulate DDX19B-mediated mRNA export (21). It remains to be determined whether Gle1 selfassociation is a common determinant of its ability to stimulate other Dbps/DDXs. However, our analysis of the phosphodeficient gle1A 6A and phosphomimetic gle1A 6D suggests that phosphorylation may dictate the oligomerization state of Gle1A and During initiation of the stress response, Gle1A in its basally modified, self-associated state activates DDX3 functions and promotes SG formation (1). As eIF2␣ is phosphorylated and the integrated stress response proceeds, the activation of MAPK cascades stimulate ERK1/2 and JNK to phosphorylate a cluster of N-terminal Gle1A residues, which then primes GSK3 to phosphorylate additional residues in the Gle1A cluster (2). Accumulation of the multiple phosphate groups over the course of the stress response alters the low-complexity character of the N terminus and perturbs Gle1 self-association (3). In this state, hyperphosphorylated Gle1A inhibits DDX3 activity and thereby promotes SG disassembly as translation resumes (4).

Phosphorylation of Gle1A regulates DDX3 and stress response
thereby modulate its regulation of DDX3 and SG dynamics. Phosphorylation often alters protein function and interactions by changing the character of intrinsically disordered, low-complexity domains (51,80,81), and the Ser 88 -Thr 102 cluster lies within a low-complexity domain of Gle1. Intriguingly, multimerization is required for the RNA helicase activity of Ded1, the yeast homologue of human DDX3 (82). It follows that Gle1A oligomerization might also be required to stimulate a multimeric form of DDX3 during the stress response.
Altered DDX3 activation may be only one of the ways that a phosphosensitive Gle1 oligomerization state could modulate function. Alternatively, but not mutually exclusive, Gle1 oligomers might serve as scaffolds for higher-order organization of SGs by promoting multivalent interactions. In a similar way, the assembly of membraneless organelles such as germ granules and paraspeckles is initiated by self-association of coiled-coil domains in PGL1/3 and Drosophila behavior/human splicing proteins, respectively (83,84). Recent work also suggests that SGs are made up of dynamic shells surrounding stable core substructures, which are thought to assemble through stable interactions between low-complexity domains in the RNAbinding proteins that decorate stalled mRNPs (16). The FRAP data presented here show that Gle1A is a stable component of SGs. Moreover, SG maturation is promoted by WT Gle1A or gle1 6A overexpression (both of which are capable of uniform oligomerization in vitro) and not by gle1 6D (Figs. 4A and 7, A-D). Thus, it is possible that a scaffold of nonhyperphosphorylated Gle1A oligomers supports SG core formation through recruitment of SG nucleators and that hyperphosphorylation dismantles this Gle1A scaffold as part of SG disassembly. However, despite its disparate ability to oligomerize in vitro, gle1 6D localization is similarly immobile compared with WT Gle1A or gle1 6A . Additional studies are required to understand the role of Gle1A oligomerization in SG biology.
It is also intriguing that the PFQ insertion in the Gle1 coiledcoil domain from the LCCS1 disease-linked gle1-Fin major allele impacts Gle1 self-association and SG dynamics in a manner comparable with that caused by aspartic acid substitutions that mimic stress-induced hyperphosphorylation in the intrinsically disordered, low-complexity region (gle1A 6D ) (21) (Fig. 7, E and  F). This suggests that the pathologies underlying LCCS1 may be due to both altered mRNA export (21) and perturbed SG dynamics. As such, multiple coincident defects in the gene expression pathway and in recovery from stress responses might contribute to the disease pathology.
Although Gle1A is localized to the cytoplasm at steady state, it shuttles in and out of the nucleus. It is unclear whether Gle1 is hyperphosphorylated in the nucleus or in the cytoplasm, but our data demonstrate that it is localized in the cytoplasm when phosphorylated. Gle1 phosphorylation is not dependent upon formation of SGs because translation inhibitors that block SG assembly (Fig. S4A) but activate MAPKs are still able to stimulate Gle1 hyperphosphorylation (Fig. S4B). This suggests that Gle1 hyperphosphorylation might also occur under conditions that do not stimulate SG assembly, and therefore may regulate other functions of Gle1. Furthermore, the longer persistence of Gle1 hyperphosphorylation through the period when eIF2␣ and 4EBP phosphorylation states have begun to reverse (Fig.   1D) suggests that the Gle1 phosphorylation state may act as a homeostatic signal.
Gle1 serves as a critical lynchpin in multiple steps of gene expression from mRNA export to translation (17). In this light, the advantage of having multiple constraints on Gle1 function is clear from the expression of functionally distinct isoforms to their discrete subcellular compartmentalization and targeting for specific interaction partners. Dynamic phosphorylation of Gle1A adds another essential layer to this complex array of regulatory mechanisms. This also reveals further potential key connections between SG and Gle1 dysfunction in neurodegenerative diseases, including ALS. Future studies will focus on combining structural and mechanistic analyses to further define how Gle1 hyperphosphorylation regulates individual steps in gene expression.

Cell culture, treatments, and reagents
HeLa cells were cultured in complete DMEM (Gibco) supplemented with 10% fetal bovine serum (Atlanta Biologicals, Flowery Branch, GA) at 37°C incubator in 5% CO 2 . For heat shock stress, cells were incubated at 45°C in the absence of CO 2 . To induce oxidative stress, cells were treated with 0.5 mM sodium arsenite for 60 min at 37°C. Kinase inhibitors were obtained from the following sources: broad-spectrum kinase inhibitor K252a (Tocris, Minneapolis, MN), JNK inhibitor JNK-IN-8 (Calbiochem, Darmstadt, Germany), ERK inhibitor UO126 (Cell Signaling Technology, Danvers, MA) and MEK-ERK inhibitor FR 180204 (Tocris). AllStars negative control siRNA and GLE1 siRNAs were supplied by Qiagen (Germantown, MD).

Phosphorylation of Gle1A regulates DDX3 and stress response Immunoblotting, immunofluorescence, and SG quantification
Immunoblotting and indirect immunofluorescence were performed as described (17). Briefly, HeLa cells were plated on 1.5-mm round coverslips in a 24-well plate (Fisher Scientific) for immunofluorescence. Following transfection and 60-min heat shock treatment at 45°C, cells were fixed with 4% paraformaldehyde and permeabilized with 0.2% Triton X-100. Coverslips were blocked with 10% FBS in PBS prior to sequential incubations with anti-G3BP antibody (1:300; BD Transduction Laboratories) and Alexa Fluor-conjugated secondary antibodies (Life Technologies). Imaging was performed using a 63ϫ 1.4 numerical aperture oil-immersion objective on a confocal microscope (Leica TCS SP5), and measurements of SG number and SG recovery were performed as described (17). Briefly, postimage processing was done using the ImageJ plug-in 3D object counter to quantify the number of SGs per cell. The 3D objects counter measurement parameters were set to "surface," and the minimum size filter was set to 3. For SG recovery, the number of cells with SGs was determined as 100 ϫ (Number of cells with G3BP-positive foci)/(Total number of cells). Prism6 (GraphPad, La Jolla, CA) was used to graph data from at least three independent experiments as mean Ϯ S.E. and perform unpaired, two-tailed Student's t tests of statistical significance. For immunoblotting, the following antibodies were used: anti-Gle1 (ASW48.

Fluorescence recovery after photobleaching
HeLa cells expressing mCHERRY-G3BP and GFP-GLE1A were heat-shocked at 45°C for 45 min prior to FRAP assay and subsequently imaged for 45 min. A Zeiss LSM710 confocal microscope was prewarmed to 45°C with 5% CO 2 ; set to 1024 ϫ 1024 frame size, eight-bit pixel depth, optical zoom of 3, and 1 airy unit; and used with a 63ϫ objective. The region of interest (ROI) for photobleaching was defined using mCherry-G3BP as a marker of SGs. Three prebleach images were acquired, and then three bleaching iterations of the ROI were performed at a scan speed of 2 using a 561 nm laser (Fig. 6) or a 488 nm laser (Fig. S7), each at 100% intensity. Postbleach images were acquired at 200-ms intervals. Localization of mCherry-G3BP (Fig. 6) or GFP-Gle1A (Fig. S7) was used to confirm SG locations in postbleach images. Analyses of FRAP recovery were performed as described previously (21) using the following equation, where F bkg is background fluorescence, F(t) ROI is intensity of bleaching ROI at each time point, F(t) cell is intensity of the entire cell at each time point, F(i) cell is initial whole-cell intensity, and F(i) ROI is initial bleaching ROI intensity. Nonregression analysis was used to determine the recovery kinetics and fitted to the following one-phase association equation using GraphPad Prism software, where Y 0 is the Y value when X (time) is zero, Plateau is the Y value at infinite times, and K is the rate constant. Span is the difference between Y 0 and Plateau and represents the mobile fraction, which was multiplied by 100 and expressed as a percentage of total fluorescence intensity. An unpaired two-tailed t test was used to test for statistical differences on half-life and rate constant.

Phos-tag gel electrophoresis
HeLa cells were lysed with radioimmunoprecipitation assay buffer (Sigma-Aldrich) supplemented with cOmplete EDTAfree protease inhibitor mixture (Roche Applied Science), 50 mM NaF, and 100 M Zn 2ϩ . A 5.5% Phos-tag acrylamide gel was prepared according to the manufacturer's instructions, and lysates were resolved by electrophoresis at 25 mA for 5-6 h at room temperature. The gel was washed two times with Trisglycine transfer buffer containing 2 mM EDTA, once with transfer buffer alone, and then wet transferred to nitrocellulose membrane overnight at 25 V at 4°C. The membrane was blocked with 5% milk blocking buffer and processed for immunoblotting with the indicated antibodies.

-Phosphatase assay
Following the indicated treatments, equivalent cell numbers were lysed in a buffer containing 50 mM Tris-HCl (pH 7.5), 150 mM NaCl, 1.5 mM MgCl 2 , 1% NP-40, and protease inhibitor mixture. Equal volumes of cell lysate were incubated at 30°C for 30 min with phosphatase buffer alone (10ϫ PPase buffer and 10ϫ MnCl 2 ) or -phosphatase enzyme (500 units; New England Biolabs) in the presence or absence of 50 mM NaF and 1ϫ PhosSTOP (Roche Applied Science). Reactions were terminated by addition of 6ϫ Laemmli sample denaturation buffer, and samples were analyzed by immunoblotting.

Subcellular fractionation
Following the indicated treatments, HeLa cells were trypsinized and collected by centrifugation at 100 ϫ g for 3 min. Nuclear and cytosolic fractions were prepared using the NE-PER kit (Thermo Fisher) following the manufacturer's instructions. In brief, cell pellets were lysed in ice-cold CER I buffer on ice for 10 min. CER II buffer was added, and lysates were vortexed and further incubated for 1 min on ice. Samples were centrifuged at 16,000 ϫ g for 7 min after which supernatants were transferred to new tubes and labeled as cytoplasmic fractions. The pellets were washed once in CER I buffer and then resuspended in NER solution to lyse the nuclei. Samples were centrifuged for 10 min at ϳ16,000 ϫ g, and supernatants were labeled as nuclear fractions.

Phosphorylation of Gle1A regulates DDX3 and stress response B-isox precipitation of proteins
Nuclear and cytosolic fractions were incubated with 100 M b-isox (Sigma-Aldrich) for 1 h at 4°C with shaking. Following incubation, reactions were centrifuged at 16,000 ϫ g for 15 min at 4°C. Pellets were washed twice with 20 mM Tris-HCl (pH 7.5), 150 mM NaCl, 5 mM MgCl 2 , 20 mM ␤-mercaptoethanol, 0.5% NP-40, and protease inhibitors. Supernatant and pellets were resuspended in Laemmli sample denaturation buffer, heated at 95°C for 5 min, and analyzed by immunoblotting.

In vitro kinase assay
For in vitro kinase assays with MAPKs, recombinant bacterially produced proteins were incubated with either 100 units of ERK enzyme (New England Biolabs) or 100 ng of JNK enzyme (Thermo Scientific, Waltham, MA) for 45 min at 30°C in 30-l reactions with the supplied kinase reaction buffer (New England Biolabs) supplemented with 100 M ATP and 1 Ci of [␥-32 P]ATP. Reactions were stopped by addition of Laemmli sample denaturation buffer, and then samples were heated at 95°C for 5 min and resolved by 7.5% SDS-gel electrophoresis. Gels were dried and analyzed by autoradiography. 32 P incorporation was detected and quantified with a Typhoon FLA7000 biomolecular imager (GE Healthcare). For in vitro kinase assays with GSK3, the indicated proteins were incubated with ERK and/or JNK in 30-l reactions with kinase buffer supplemented with 100 M ATP for 60 min at 30°C. GSK-3 (125 units; New England Biolabs) and 1 Ci of [␥-32 P]ATP were then added to the reactions in the presence or absence of ERK and or JNK inhibitors, and all reaction volumes were adjusted to 40 l. After incubation for 45 min at 30°C, reactions were terminated by addition of Laemmli sample denaturation buffer. Samples were heated at 95°C for 5 min, resolved by SDS-PAGE, and analyzed by autoradiography.

Immunoprecipitation
24 h post-transfection, HeLa cells were either left untreated or treated with 0.5 mM sodium arsenite for 60 min at 37°C. Cells were then lysed in 500 l of lysis buffer containing 20 mM Tris-HCl (pH 7.5), 150 mM NaCl, 2 mM MgCl 2 , 1 mM EDTA, 0.5% NP-40, and protease and phosphatase inhibitors. Lysates were spun at 16,000 ϫ g for 10 min at 4°C. Agarose beads conjugated to anti-GFP V H H single domain antibody (GFP-TRAP, ChromoTek, Hauppauge, NY) were added to the supernatant, and the slurry was incubated with shaking for 2 h at 4°C. Bound complexes were pelleted at 2500 ϫ g for 1 min, washed five times with lysis buffer, and eluted with 2ϫ Laemmli sample denaturation buffer. Proteins were resolved by 7.5% SDS-gel electrophoresis, and the gel was stained with colloidal blue. Protein bands corresponding to GFP-Gle1A were excised from the gel and processed for MS.

Mass spectrometry
Acrylamide gel samples containing GFP-Gle1A were cut into 1-mm 3 pieces and treated with 45 mM DTT, and available Cys residues were carbamidomethylated with 100 mM iodoacetamide. The samples were destained with 50% MeCN in 25 mM ammonium bicarbonate, and proteins were digested with tryp-sin (10 ng/l) in 25 mM ammonium bicarbonate overnight at 37°C. Peptides were extracted by gel dehydration (60% MeCN, 0.1% TFA), the extract was dried by SpeedVac centrifugation, and peptides were reconstituted in 0.1% formic acid. The peptide solutions were then loaded onto a capillary reverse-phase analytical column (360-m outer diameter ϫ 100-m inner diameter) using an Eksigent NanoLC HPLC and autosampler. The analytical column was packed with C 18 reverse-phase resin (Jupiter, 3-m beads, 300 Å, Phenomenex, Torrance, CA) directly into a laser-pulled emitter tip. Mobile-phase solvents consisted of 0.1% formic acid, 99.9% water (solvent A) and 0.1% formic acid, 99.9% acetonitrile (solvent B). A 90-min gradient was performed, and eluting peptides were mass-analyzed on an LTQ Orbitrap Velos mass spectrometer (Thermo Scientific) equipped with a nanoelectrospray ionization source. The instrument was operated using a data-dependent method with dynamic exclusion enabled. Full-scan (m/z 300 -2000) spectra were acquired with the Orbitrap, and the top 16 most abundant ions in each MS scan were selected for fragmentation via collision-induced dissociation in the LTQ. An isolation width of 2 m/z, activation time of 10 ms, and 35% normalized collision energy were used to generate MS2 spectra. For additional LC-MS/MS analysis, the instrument was operated using a combination method of data-dependent (top 12) and targeted scan events. Targeted scans were acquired for m/z 985.43 and m/z 657.29, which correspond to phosphorylated Gle1 peptide SPDASSAFSPASPATPNGTK. For identification of peptides, tandem mass spectra were searched with Sequest (Thermo Scientific) against a human subset database created from the Uni-ProtKB protein database. Variable modifications of ϩ57.0214 on Cys (carbamidomethylation); ϩ15.9949 on Met (oxidation); and ϩ79.9663 on Ser, Thr, and Tyr (phosphorylation) were included for database searching. Search results were assembled using Scaffold 4.3.2 (Proteome Software, Portland, OR).

Protein purification
Plasmids were transformed into Escherichia coli Rosetta cells and cultured in Terrific Broth under kanamycin and chloramphenicol selection. Expression was induced with 0.2 mM isopropyl 1-thio-␤-D-galactopyranoside at an A 600 of 0.8 for 18 h at 18°C. Bacteria were lysed by sonication in 20 mM HEPES (pH 7.5), 400 mM NaCl, 0.5 mM TCEP, and 20% glycerol supplemented with cOmplete EDTA-free protease inhibitor mixture and 2 mM phenylmethylsulfonyl fluoride. MBP-Gle1A 1-360 -His was purified by nickel-immobilized metal affinity chromatography (IMAC). His-PPS was added to the pooled peak elution fractions to cleave the His-MBP tag, and fractions were dialyzed overnight to remove imidazole. Gle1A 1-360 was then separated from the tag and His-PPS by nickel-IMAC (Thermo Scientific) and applied to an S200 column equilibrated with buffer containing 20 mM HEPES (pH 7.5), 150 mM NaCl, 20% glycerol, and 0.5 mM TCEP for final purification. His-MBP-tagged fulllength Gle1A proteins were purified by amylose resin chromatography according to the manufacturer's recommendations followed by nickel-IMAC. His-tagged DDX3 was purified by nickel-IMAC and dialyzed overnight in buffer containing 20 mM HEPES (pH 7.5), 50 mM NaCl, 0.5 mM TCEP, 2 mM MgCl 2 , and 20% glycerol. The protein was further purified by heparin Phosphorylation of Gle1A regulates DDX3 and stress response chromatography and eluted with a salt gradient from 50 to 1000 mM NaCl.

Electron microscopy
Uranyl formate-stained samples were prepared as described (85). Samples were imaged on an FEI Morgagni electron microscope operated at an acceleration voltage of 100 kV. Images were recorded at a magnification of 28,000ϫ and collected using a 1000 ϫ 1000 charge-coupled device camera (ATM). Particle shape was quantified by taking the ratio of the longest axis (major axis) and the axis perpendicular to the major axis (minor axis).

ATPase assay
Colorimetric enzyme-coupled ATPase rate assays were performed as described (86).