HAM-5 Functions As a MAP Kinase Scaffold during Cell Fusion in
Cell fusion between genetically identical cells of the fungus Neurospora crassa occurs when germinating asexual cells (conidia) sense each other's proximity and redirect their growth. Chemotropic growth is dependent upon the assembly of a MAPK cascade (NRC-1/MEK-2/MAK-2) at the cell cortex (conidial anastomosis tubes; CATs), followed by disassembly over an ∼8 min cycle. A second protein required for fusion, SO, also assembles and disassembles at CAT tips during chemotropic growth, but with perfectly opposite dynamics to the MAK-2 complex. This process of germling chemotropism, oscillation and cell fusion is regulated by many genes and is poorly understood. Via a phosphoproteomics approach, we identify HAM-5, which functions as a scaffold for the MAK-2 signal transduction complex. HAM-5 is required for assembly/disassembly and oscillation of the MAK-2 complex during chemotropic growth. Our data supports a model whereby regulated modification of HAM-5 controls the disassembly of the MAK-2 MAPK complex and is essential for modulating the tempo of oscillation during chemotropic interactions.
Published in the journal:
. PLoS Genet 10(11): e32767. doi:10.1371/journal.pgen.1004783
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pgen.1004783
Summary
Cell fusion between genetically identical cells of the fungus Neurospora crassa occurs when germinating asexual cells (conidia) sense each other's proximity and redirect their growth. Chemotropic growth is dependent upon the assembly of a MAPK cascade (NRC-1/MEK-2/MAK-2) at the cell cortex (conidial anastomosis tubes; CATs), followed by disassembly over an ∼8 min cycle. A second protein required for fusion, SO, also assembles and disassembles at CAT tips during chemotropic growth, but with perfectly opposite dynamics to the MAK-2 complex. This process of germling chemotropism, oscillation and cell fusion is regulated by many genes and is poorly understood. Via a phosphoproteomics approach, we identify HAM-5, which functions as a scaffold for the MAK-2 signal transduction complex. HAM-5 is required for assembly/disassembly and oscillation of the MAK-2 complex during chemotropic growth. Our data supports a model whereby regulated modification of HAM-5 controls the disassembly of the MAK-2 MAPK complex and is essential for modulating the tempo of oscillation during chemotropic interactions.
Introduction
Fusion between genetically identical cells occurs in many different organisms and plays pivotal roles in different developmental processes, such as myoblast fusion during muscle formation, macrophage fusion involved in tissue remodeling and fusion of trophoblasts during placental development [1], [2]. Cell fusion is also important for the formation of the interconnected mycelial network that is the hallmark of filamentous fungal growth [3], [4], [5]. In addition to hyphal fusion, fusion can also occur between genetically identical germinating asexual spores (conidia) of filamentous fungi [6], [7], [8]. Both hyphal and germling fusion are integral to the formation of an interconnected hyphal network and impart fitness benefits, as well as mediating genetic mixing and the sharing of resources [3], [4], [9], [10], [11], [12], [13].
In the filamentous ascomycete fungus, Neurospora crassa, germinated conidia (germlings) in close proximity fuse via specialized conidial anastomosis tubes (CATs) that form at germ tube tips or between conidia [14]. In a fungal colony, hyphal fusion is observed in the central parts of the colony, in contrast to the peripheral parts where hyphae avoid each other. A large number of genes have been identified in N. crassa that are important for the process of sensing, chemotropic interactions and CAT fusion and have contributed to an understanding of this complex developmental system in filamentous ascomycete fungi [7], [15], [16], [17]. An essential part of chemotropic interactions in N. crassa is the oscillatory recruitment of three kinases of a MAPK cascade (NRC-1, MEK-2 and MAK-2) and of a protein of unknown function, SOFT (SO), to CAT tips [18], [19]. In strains carrying loss-of-function mutations in these genes, oscillatory recruitment of NRC-1/MEK-2/MAK-2 or SO, chemotropic interactions and fusion do not occur [8], [18], [19], [20]. It was proposed that the alternating oscillation of MAK-2 and SO to CAT tips may function to establish two distinct physiological states in interacting germlings to enable chemotropism to persist, avoid self-stimulation and assure a rapid and efficient cell fusion [19], [21]. Recently, it has been shown that an ortholog of SOFT in the related filamentous ascomycete species, Sordaria macrospora, PRO40, is a scaffold protein for the cell wall integrity MAP kinase pathway (MAK-1) [22]; Δmak-1 mutants in both N. crassa and S. macrospora are fusion mutants [22], [23].
Previously, it was shown that kinase activity of MAK-2 is required to maintain oscillatory recruitment of both MAK-2 and SO; addition of ATP-analog 1NM-PP1 to a mutant containing an inhibitable MAK-2 protein encoded by a mak-2Q100G allele, disrupted the oscillation of both MAK-2 and SO in communicating germlings and stalled chemotropic interactions and the fusion process [19]. The use of mak-2Q100G strain also contributed to the identification of downstream genes whose expression levels depend on functional MAK-2 [24]. However, MAK-2 kinase targets involved in the oscillation process have not been identified and the characterization of such potential targets could help unravel molecular mechanisms associated with this highly regulated and complex process.
In recent years, highly sensitive liquid chromatography-mass spectrometry (LC-MS) based quantitative phosphoproteomic techniques have contributed to our understanding of kinase pathway function in eukaryotic cells [25], [26], [27], [28]. To identify MAK-2 kinase targets in N. crassa, we took a global approach by identifying phosphopeptides in mak-2Q100G germlings, treated or not with 1NM-PP1. From the phosphoproteomic screen, a number of candidate MAK-2 target proteins were identified, one of which encodes a protein previously identified as being essential for germling/hyphal fusion, HAM-5 [23], [29]. We show that HAM-5 oscillates during chemotropic interactions with components of the MAK-2 pathway, physically interacts with NRC-1, MEK-2 and MAK-2 and was required for localization of MAK-2 and MEK-2 to puncta. Our data supports the hypothesis that HAM-5 functions as a scaffold protein by binding to and co-localizing to CAT tips with all three kinases in the MAK-2 cascade during chemotropic growth in germlings as well as in hyphae undergoing fusion events. These studies shed new light on the mechanisms of oscillation during communication and chemotropic interactions between genetically identical cells and which may be important for function of this conserved MAPK pathway in other filamentous ascomycete fungi.
Results
Identification of MAK-2 targets in Δmak-2Q100G germlings using a phosphoproteomic approach
To better understand the role of MAK-2 during chemotropic interactions, we set out to identify putative kinase targets using a global quantitative phosphoproteomics approach using a strain carrying an inhibitable mak-2Q100G allele. Altering a specific amino acid in the ATP binding site (glutamine for glycine) renders MAK-2 sensitive to inhibition to the ATP analogue 1NM-PP1, but does not affect MAK-2 kinase activity in the absence of inhibitor [19]. In the absence of the inhibitor 1NM-PP1, strains and germlings containing the mak-2Q100G allele (his3::mak-2Q100G; Δmak-2) showed wild-type growth and fused normally, while in the presence of inhibitor, the mak-2Q100G strain showed a mutant mak-2 phenotype; chemotropic interactions and fusion were not observed, consistent with inactivation of MAK-2 kinase activity [19]. For identifying MAK-2-dependent phosphorylation events, 4.5-hr old mak-2Q100G germlings were treated with DMSO (two samples and two biological replicates) or 1NM-PP1 (dissolved in DMSO; two samples and two biological replicates) for 10 min. Proteins were extracted, digested, and enriched phosphopeptides identified and quantified using isobaric peptide tags for relative and absolute quantification (iTRAQ)-based LC-MS/MS analyses. Although the same peptides across experimental conditions are labeled with different iTRAQ reagents and indistinguishable by mass, different masses will be generated in the tandem MS by releasing the reporter ions for the 4-plex iTRAQ method (Figure S1). We performed the full experiment twice, resulting in a total of eight samples, which were analyzed for phosphopeptide identity and abundance.
From these experiments a total of 3200 unique phosphopeptides were identified. These 3200 unique phosphopeptides originated from 1164 proteins (Dataset S1). A small percentage of the identified peptides have multiple phosphorylation sites (12%, Figure 1B), as compared to peptides where only a single phosphorylation site was identified (88%, Figure 1B). Phosphorylation sites were predominantly identified on serine residues (75%), and to a lesser extent on threonine (22%) and tyrosine (3%) residues (Figure 1C). FunCat analysis [30] showed the set of identified phosphorylated proteins in germlings originated from a wide spectrum of functional categories (Figure 1A), including, not unexpectedly, proteins involved in metabolism, energy, cell cycle and DNA processing, protein synthesis and transcription. However, proteins within the functional categories of cellular communication/signal transduction, cell defense, and interaction with the environment were also identified, suggesting germlings are poised to respond to variations in their environment (Dataset S1).
Of the 3200 phosphopeptides, 96 unique phosphopeptides from 67 proteins were>1.5 times less abundant (p<0.05) in at least one replica experiment in the 1NM-PP1 treated mak-2Q100G germlings relative to the DMSO-treated control cells (Table S1). Functional category analyses [30] of this set of proteins showed enrichment for genes involved in metabolism, energy, cell cycle and DNA processing, transcription, protein fate, regulation of metabolism and protein function, cellular communication/signal transduction, interaction with the environment, cell fate and cell type differentiation (p<0.01) (Figure 1A; Dataset S1). Three proteins in this group had previously been shown to be required for hyphal fusion, including HAM-9, HAM-11 and MAK-1 [23], [24], [31], [32]. In addition to the MAP kinase involved in the osmotic response signaling (OS2) [32], CUT-1, which is implicated in the osmotic stress response [33], [34] as well as the transcription factor that is a target of the OS-2 pathway (ASL-1/ATF-1, NCU01345) [35] were identified (Table 1; Table S1). Of these 3200 unique phosphopeptides, 33 phosphopeptides (from 27 proteins) were identified with an increased abundance of at least 1.5 fold (p<0.05) in the 1NM-PP1 treated germlings in at least one replica experiment. Functional category analyses [30] of this set of 27 proteins showed over-representation for genes/proteins involved in phosphate metabolism (perhaps as a response to 1NM-PP1 exposure), protein synthesis, protein fate and protein with binding function or cofactor requirements (Dataset S1).
To identify potential direct MAK-2 targets, we inspected the 96 phosphopeptides that showed reduced abundance in the 1NM-PP1 treated mak-2Q100G germlings for MAPK consensus phosphorylation sites (P-X-S*/T*-P) [36], [37]. Nine proteins were identified, of which five were annotated as hypothetical proteins. One of the proteins was ASL-1, the transcription factor that is a target of the OS-2 pathway and which shows an ascospore-lethal phenotype [35]; Δmak-2 mutants also show an ascospore-lethal phenotype [8], [38]. The remaining proteins included a predicted 6-phosphofructo-2-kinase (NCU01728), a predicted trehalose phosphatase (NCU05041) and a protein previously reported to be required for hyphal fusion, HAM-5 [23], [29] (Table 1). In addition to a predicted MAPK phosphorylation site, a predicted MAPK docking motif (R/K-R/K-(X)1-5-I/L-X-I/L) [39] was predicted in NCU07868 (hypothetical protein) and HAM-5 (Table 1). Of the nine genes encoding potential MAK-2 phosphorylation targets, five showed a reduction in expression levels in either a Δpp-1 mutant (transcription factor that is a target of the MAK-2 pathway) or in a mak-2Q100G germlings treated with 1NM-PP1 [24], including ham-5 (Table 1).
To determine whether the putative direct or indirect MAK-2 phosphorylation targets were involved in germling fusion, strains carrying individual deletions of 8 out of the 9 genes whose proteins contained a MAPK consensus phosphorylation site and showed decreased abundance in the 1NM-PP1-treated mak-2Q100G germlings (including strains carrying deletions in all five hypothetical proteins) (Table 1), plus an additional 31 deletion mutants selected from a subset of proteins that showed decreased abundance in 1NM-PP1-treated mak-2Q100G germlings, but which lacked a predicted MAPK consensus phosphorylation site (Table S1), were evaluated for the ability to undergo chemotropic interactions and cell fusion. Of these, Δham-5 in the first set, and Δham-9, Δham-11, and Δmak-1 strains in the second set, were fusion defective.
Putative MAK-2 target HAM-5-GFP localizes in an oscillatory manner to cell tips during chemotrophic interactions
To further characterize putative MAK-2 targets that are required for germling fusion, we GFP-tagged proteins where localization during chemotropic interactions/cell fusion had not been determined. This effort included HAM-5. The ham-5 mutant was identified in a forward screen for mutants that failed to form a heterokaryon [29] and encodes a large protein of 1686 amino acids (aa) with seven putative WD40 repeats at the N-terminus (aa 14-313, grey boxes; Figure 2A). At the C-terminus, an unstructured region of low complexity was identified (shaded white boxes; Figure 2A) that contains stretches of predominately proline and glutamine residues. Two coiled coil domains were also predicted in HAM-5 (aa 1168-1190 and 1257-1286, red boxes; Figure 2A). In total, 16 phosphorylation sites were identified, located mainly at the middle section of HAM-5 (Figure 2A). Three phosphorylation sites were identified from our phosphoproteomics analysis, one of which was a putative MAPK site at amino acid residue 506 (serine) (Figure 2A; Table 1). Thirteen additional HAM-5 phosphopeptide sites were identified in a recent phosphoproteomics study of N. crassa hyphal cultures exposed to different carbon sources [40]. A putative MAPK docking site was also identified between residues 1128-1136 (RRKPPALDL) in the C-terminus of HAM-5 (yellow bar; Figure 2A).
In conidia, germlings and in hyphae, HAM-5-GFP fluorescence was observed in small, cytoplasmically localized puncta (white arrows, Figure S2). In mature hyphae, localization to septa was also observed (red arrow, Figure S2A). During chemotropic interactions between germlings, HAM-5-GFP driven by the tef-1 promoter was observed as small puncta in interacting germlings and also at the tip (Figure 2B). Importantly, HAM-5 oscillated to CAT tips in germling pairs during chemotropic interactions, with localization dynamics similar to MAK-2 or SO: HAM-5-GFP appeared at the CAT tip of one germling, but was absent from the CAT tip in the partner germling, while 4 min later, HAM-5-GFP was present at the CAT tip of the second partner germling, but was absent from the CAT tip of the first germling (Movie S1). During the fusion process, a HAM-5-GFP signal was also observed in puncta that localized either to the cell periphery or to points close to nuclear compartments, which were devoid of HAM-5-GFP (Figure 2B, asterisk). When HAM-5-GFP localized to the CAT tips, the number of cytoplasmically localized puncta in the partner germling was reduced and the GFP signal was more dispersed in the cytoplasm (Figure 2B). HAM-5-GFP was detected at the site of germling contact and remained there until the cytoplasm of the two germlings mixed (Figure 2C).
To assess whether the HAM-5-GFP-tagged versions were biologically functional, we crossed wild type (WT) strains carrying ham-5-gfp driven by either the native, tef-1 or ccg-1 promoter into the Δham-5 mutant. Progeny bearing both ham-5-gfp and Δham-5 showed a fusion phenotype and frequency similar to WT, indicating that ham-5-gfp driven either by the native, ccg-1 or tef-1 promoter was functional. Localization of HAM-5-GFP driven by either ccg-1 or tef-1 in the Δham-5 strain was similar to HAM-5-GFP in the WT strain, and showed localization to puncta and oscillation to CAT tips during chemotropic interactions. ham-5-gfp driven by its native promoter in the Δham-5 strain also showed localization to puncta (Figure S3). However, due to low expression levels of HAM-5-GFP in this strain, localization during chemotropic interactions could not be fully assessed. In the following sections, HAM-5-GFP is shown driven by the tef-1 promoter unless stated otherwise. A heterokaryon of a strain carrying ham-5-gfp driven by tef-1 and a strain carrying histone H1-dsRED (a nuclear marker), showed non-overlapping fluorescence (Figure S2), indicating that HAM-5-GFP was excluded from the nucleus. This localization pattern is similar to SO-GFP [19], but different than MAK-2-GFP, which localizes to the cytoplasm and to nuclei (Figure S2) [19].
MAK-2 regulates HAM-5 dynamics
HAM-5 is predicted to be a highly phosphorylated protein, with 16 predicted phosphorylation sites (Figure 2A). To assess whether HAM-5 is a phosphorylation target of MAK-2, as indicated by the phosphoproteomics data, we introduced ham-5-GFP into a Δmak-2 strain and determined its phosphorylation status during germling fusion. HAM-5-GFP was immunoprecipitated from WT (ham-5-gfp) and Δmak-2 (ham-5-gfp) 5 hr-old germlings using anti-GFP antibodies and assayed for phosphorylation status using Western blot analysis with anti-phosphoserine/threonine antibodies that specifically recognize phosphoserine or phosphothreonine sites followed by a proline residue (MAPK phosphorylation sites). The results showed that HAM-5-GFP from both WT and Δmak-2 cells was specifically phosphorylated (Figure 2D). The HAM-5-GFP protein band from WT germlings showed a slight smear upwards and was slightly larger than that observed in Δmak-2 (ham-5-gfp) germlings (Figure 2D). These data support the phosphoproteomics data, which suggested a MAK-2-dependent modification of HAM-5 during germling fusion (Table 1). However, it is clear that HAM-5-GFP was also phosphorylated in the Δmak-2 mutant, suggested possible additional regulatory inputs into HAM-5 via phosphorylation by other proteins during conidial germination/germling fusion.
To assess whether the identified phosphosite in HAM-5 (serine 506) was required for HAM-5 function during chemotropic interactions and germling fusion, strains carrying site-directed mutations whereby serine 506 was altered to an alanine (phosphorylation impaired) or a glutamate (phosphorylation mimic) residue were evaluated: chemotropic behavior and cell fusion in germlings carrying the ham-5S506A or ham-5S506E mutations driven by the ccg-1 promoter were indistinguishable from wild type germlings and fully complemented the growth phenotype of the Δham-5 strain (Figure S4A). To assess whether mutations in the predicted MAPK docking site affected HAM-5 function, a strain was constructed where the first three amino acids of the MAPK docking site were changed to alanine (RRKPPALDL to AAAPPALDL). However, germlings bearing this ham-5 allele (ham-5RRK1128AAA) also showed a similar communication and fusion phenotype to WT germlings resulting in similar growth phenotype to WT and complemented fully the growth phenotype of the Δham-5 strain (Figure S4A). Thus, the predicted MAPK docking site in HAM-5 might not be functional or additional MAPK docking sites are present making this site redundant for function.
We predicted that HAM-5-GFP would show altered localization when introduced into Δmak-2 germlings, due to the inability of these cells to undergo chemotropic interactions. However, HAM-5-GFP localized to puncta in Δmak-2; ham-5-gfp germlings, as observed in wild type cells (Figure 2E). However, no oscillation of HAM-5-GFP puncta to cell tips was observed, consistent with the lack of fusion and chemotropic interactions in Δmak-2 germlings. HAM-5-GFP puncta were localized randomly within the cell, with some puncta close to the nuclear periphery and membrane, but cells also contained at least one HAM-5-GFP puncta localized to the cell tip (Figure 2E). The stable, tip-anchored localization of HAM-5-GFP (Figure 2E, arrows) was not observed in WT germlings in the absence of chemotropic interactions (in isolated germlings), and was unique to Δmak-2; ham-5-gfp cells.
HAM-5-GFP oscillates with the MAP kinase cascade members to CAT tips
The observation that HAM-5-GFP showed oscillation during germling fusion, but still localized to puncta in Δmak-2 germlings, suggested that either HAM-5 interacted with SO or with other proteins in the MAK-2 signal transduction pathway. Previously, it was shown that the components of the MAK-2 pathway, the MAPKKK (NRC-1) and the MAPKK (MEK-2), oscillate with MAK-2 during chemotropic interactions in germlings [18]. To determine if HAM-5 oscillated with the components of the MAK-2 complex or with SO, we used heterokaryons of strains carrying ham-5-gfp and either mCherry-tagged mak-2, mek-2, nrc-1 or so alleles and examined co-localization of these proteins during chemotropic interactions in germlings. As shown in Figure 3, HAM-5-GFP + MAK-2-mCherry, or HAM-5-GFP + MEK-2-mCherry or HAM-5-GFP + NRC-1-mCherry were co-recruited to the CAT tips during chemotropic interactions. In homokaryotic strains carrying both HAM-5-GFP and MAK-2-mCherry, the dynamics of the two proteins were identical and showed simultaneous oscillatory recruitment to CAT tips (Movie S2). HAM-5-GFP + MAK-2-mCherry and HAM-5-GFP + MEK-2-mCherry also co-localized to cytoplasmic puncta during chemotropic interactions (Figure 3). As observed in [18], the NRC-1-mCherry signal was weak, and recruitment of HAM-5-GFP + NRC-1-mCherry to cytoplasmic puncta could not be assessed. In contrast to MAK-2/MEK-2/NRC-1, SO-mCherry + HAM-5-GFP always appeared at opposite CAT tips (Figure 3D) and with exactly opposite dynamics during chemotropic interactions (Movie S3).
The localization of HAM-5 during chemotropic interactions suggested that HAM-5 physically interacts with MAK-2, MEK-2 and/or NRC-1. To test this hypothesis, we performed co-immunoprecipitation experiments using strains carrying HAM-5-GFP + MAK-2-mCherry, HAM-5-GFP + MEK-2-mCherry or HAM-5-GFP + NRC-1-mCherry. As controls, we used strains carrying MAK-2-mCherry, MEK-2-mCherry or NRC-1-mCherry-tagged proteins with GFP driven by the ccg-1 promoter; GFP in these strains showed only cytoplasmic localization and was never observed in puncta. A specific interaction between HAM-5-GFP and MEK-2-mCherry and HAM-5-GFP and NRC-1-mCherry was detected when HAM-5-GFP was immunoprecipitated using anti-GFP antibodies from 5 hr-old germlings and subsequently re-probed using anti-mCherry antibodies (Figure 3E), consistent with the co-localization of these proteins observed by confocal microscopy (Figure 3A–C). No interaction between the mCherry-tagged proteins and cytoplasmic GFP was observed (Figure 3E; Figure S4C). An interaction between MAK-2-mCherry and HAM-5 could not be assessed, as MAK-2-mCherry showed non-specific binding. However, using phospho-specific anti-P42/P44 (Erk1/Erk2) antibodies that recognize phosphorylated MAK-2 [8], an interaction between HAM-5-GFP and MAK-2 was detected via co-immunoprecipitation (Figure 4D). By contrast, no interaction was detected between HAM-5-GFP and SO-mCherry (Figure 3E).
The HAM-5 WD40 domain is required for function and stability and binds to MAK-2 but not to MEK-2
HAM-5 is a large protein with predicted protein-protein interaction domains including seven WD40 repeats, which are predicted to form β-propeller structures that have been implicated in coordinating protein assemblages in other systems [41]. To test the hypothesis that the WD40 domain is involved in HAM-5-MAK/MEK-2/NRC-1 interactions, we constructed two mutant ham-5 alleles: one in which the WD40 motifs (aa 67-348) were removed (HAM-5Δ67-348) and one in which only the first 351 aa including the WD40 motifs of HAM-5 were retained (HAM-51-351). Both alleles were tagged with gfp, and function and localization were assessed in both WT and Δham-5 mutant strains.
The ham-51-351 construct failed to complement the growth or fusion defects of the Δham-5 mutant (Figure S4A). When observed microscopically, the localization of the HAM-51-351-GFP in Δham-5 germlings was cytoplasmic and nuclear; no puncta were observed (Figure 4B; Figure S5A). This result is in contrast to the full length HAM-5-GFP, which localized to puncta and was excluded from the nucleus (Figure 2; Figure S2). However, in a WT background, a portion of the HAM-51-351-GFP localized to puncta and showed oscillation during chemotropic interactions between germlings (Figure 4A). These observations suggest that in WT germlings, HAM-51-351-GFP may bind the native untagged HAM-5, resulting in localization to puncta when the complex oscillates to CAT tips during chemotropic interactions. To determine whether HAM-51-351-GFP (in a Δham-5 mutant) physically interacted with MAK-2 or MEK-2, we performed co-immunoprecipitation experiments using either anti-mCherry antibodies (for MEK-2-mCherry) or anti-p42/44 antibodies for MAK-2 [8]. HAM-51-351 specifically immunoprecipitated phosphorylated MAK-2 but not MEK-2-mCherry (Figure 4C, D; Figure S4D).
The ham-5Δ67-348 construct also failed to complement the growth or fusion defects of the Δham-5 mutant (Figure S4A), consistent with an essential role for the HAM-5 WD40 domain. Cellular fluorescence of HAM-5Δ67-348-GFP in germlings or hyphae was not observed and less protein was produced than other GFP-tagged proteins (Figure S4B), suggesting that the WD40 domain is required for HAM-5-GFP stability. However, co-immunoprecipitation experiments revealed a specific interaction between HAM-5Δ67-348 and MEK-2-mCherry, but not with MAK-2 (Figure 4C; Figure S4F). These biochemical interaction studies indicated that the WD40 domain of HAM-5 is important for interactions with MAK-2 and the C-terminus is important for interactions with MEK-2. The predicted MAPK docking site, which was not required for function by mutational analyses (see above), is located in the C-terminus of HAM-5, indicating that this site is not essential for MAK-2-HAM-5 interactions.
HAM-5 is required for assembly of the MAPK-cascade members MAK-2 and MEK-2 into puncta
Scaffold proteins such as Ste5 in Saccharomyces cerevisiae [42], which assembles the pheromone response MAPK pathway, regulates spatial functionality of this pathway via nuclear/plasma membrane shuttling of Fus3 during mating. We hypothesized that HAM-5 was also required for the assembly of the MAK-2 cascade members in complexes and subsequent recruitment of these complexes to their correct cellular location during chemotropic interactions (e.g. the CAT tip). To test this hypothesis, we introduced MAK-2-GFP or MEK-2-mCherry into a Δham-5 strain. In contrast to WT and Δmak-2 germlings where HAM-5-GFP was localized to puncta (Figure 2), MAK-2-GFP showed cytoplasmic and nuclear localization in the Δham-5 mutant; no puncta were observed (Figure 5A, B). In Δham-5 hyphae, MEK-2-mCherry localized to the cytoplasm and to septa, but puncta were not observed as in WT hyphae (Figure 5C, D; Movie S4). We then tested whether HAM-5 was required to establish a stable interaction between MEK-2 and MAK-2. When MEK-2-mCherry was immunoprecipitated from WT germlings, phosphorylated MAK-2 was also detected, while in Δham-5 germlings, co-immunoprecipitation of phosphorylated MAK-2 with MEK-2-mCherry was not detectable (Figure 5E; Figure S4E). SO-GFP was also cytoplasmically localized in Δham-5 (so-gfp) germlings (Figure S5), a localization pattern identical to that observed in WT (so-gfp) germlings not undergoing chemotropic interactions.
Phosphorylation of MAK-2 by the upstream kinase, MEK-2, is required for fusion [18], but is not fully dependent on functional HAM-5. In a Δham-5 mutant, phosphorylated MAK-2 is still detectable in hyphal preparations [29]. We confirmed this result in germlings, where phosphorylated MAK-2 was also observed in the Δham-5 samples (Figure 6A). MAK-2 phosphorylation is also observed in two other fusion mutant strains: Δham-7 and Δham-11 [24], [43] (Figure 6A). To investigate whether the localization of MAK-2 complexes to puncta was dependent on HAM-5 or the ability to undergo chemotropic interactions, we expressed HAM-5-GFP and MAK-2-mCherry in Δham-7 and Δham-11 cells. Although Δham-7 and Δham-11 germlings are unable to undergo chemotropic interactions and cell fusion [23], [24], [43], co-localization of HAM-5-GFP and MAK-2-mCherry to puncta was still observed (Figure 6B-D). Interestingly, as seen in the Δmak-2; ham-5-gfp strain, tip-anchored HAM-5-GFP and MAK-2-mCherry were present in Δham-7 and Δham-11 germlings (Figure 6, arrows).
HAM-5-GFP oscillates with MAK-2 during hyphal fusion
Most mutants affected in germling fusion are also deficient in hyphal fusion [7], although localization of MAK-2 or SO during hyphal interactions has not been previously reported. Whether hyphal fusion is similarly coordinated as germling fusion is unknown, as different avoidance and fusion signals may be present at the periphery and older parts of a colony [44]. Another difference between germlings and hyphae is the presence of cytoplasmic flow in hyphae [10] that may influence the oscillation of proteins to sites of fusion. We therefore evaluated the localization of HAM-5-GFP during hyphal fusion in a mature colony. As shown in Figure 7, oscillation of HAM-5-GFP to the tips of hyphae undergoing chemotropic interactions was observed with dynamics very similar to that during germling fusion (where MAK-2 has a cycling time of ∼8 min at a single hyphal tip). Similar to germlings, MAK-2-mCherry also showed co-localization with HAM-5-GFP and oscillated with identical dynamics. MAK-2-mCherry was also observed in nuclei, while HAM-5-GFP was excluded (Figure S6A and Movie S5).
In addition to localization to sites at fusion tips of hyphae, puncta containing HAM-5-GFP and MAK-2-mCherry within adjoining hyphal compartments also showed oscillation (Figure 7 and Figure S6). Interestingly, upon membrane merger (t = 22 min), oscillation in both fusion hyphae was completely coordinated for an additional 30 minutes (Figure 7C,D, Figure S6 and Movie S5). We further assessed how far oscillation of HAM-5-GFP puncta extended in hyphal compartments that surrounded a fusion point. In filamentous fungi like N. crassa, hyphal compartments are delineated by septa, but septa contain a pore through which organelles, including nuclei, can move [45]. We observed coordinated oscillation of HAM-5 in over six hyphal compartments that were ∼100 µm from the point of fusion for a total distance of ∼200 µm (Figure S7). Hyphal compartments that showed different or no oscillation of HAM-5 distant from the point of fusion were delineated by septa (Movie S6). These data show that the oscillation of HAM-5-GFP in the hyphal network was coordinated over large distances surrounding a fusion point, but could be restricted by septa. Cytoplasmic flow, septal plugging and fusion events in nearby hyphae may also affect the oscillation of HAM-5 in compartments surrounding hyphal fusion points.
Discussion
In this study, we show that HAM-5 functions as a scaffold protein for the MAK-2 MAP kinase complex and is required for oscillation of this complex during chemotropic interactions during germling and hyphal fusion in N. crassa. Our findings are complemented by the accompanying study of Dettmann et al., [46] that show physical interaction between HAM-5 and MAK-2/MEK-2/NRC-1 via mass spectrometry and yeast two hybrid, assessing both indirect and direct physical interactions of HAM-5 with the MAK-2 kinase complex. In other filamentous ascomycete species, mutations in nrc-1, mek-2, mak-2 orthologs results in strains unable to undergo vegetative cell fusion [47], [48], [49], [50] as well as defects in growth, reproduction, virulence and host colonization phenotypes, indicating expanded functions for this MAPK pathway in filamentous fungi as compared to yeast [51], [52], [53], [54], [55], [56]. ham-5 is highly conserved in the genomes of filamentous ascomycete species [51]. We predict that these ham-5 homologs will function as a scaffold in these species for mak-2/mek-2/nrc-1 orthologs, and which may be important for mediating growth, reproduction and virulence functions of this important and conserved signal transduction pathway.
The MAK-2 MAPK signal transduction pathway (MAK-2, MEK-2 and NRC-1) in filamentous fungi is orthologous to the pheromone response pathway in S. cerevisiae (Fus3, Ste7 and Ste11). Previously, it was shown that a FUS3/KSS1 ortholog in the filamentous ascomycete species Magnaporthe grisea, (PMK1) as well as its ortholog in Aspergillus nidulans (mpkB) complements the pheromone response/mating defect of a S. cerevisiae fus3Δ/kss1Δ mutant [49], [57]. In S. cerevisiae, the Ste5 scaffold protein allosterically facilitates Ste7 phosphorylation of N. crassa MAK-2 (called N. cra mpkB) [58] and A. nidulans MpkB [57], indicating conservation of regulation of these kinases by allosteric motifs within Ste5. However, although components of these two MAPK pathways are highly homologous in fungi, an ortholog of STE5 is absent in the genomes of filamentous ascomycete species. Future experiments comparing the function of these non-homologous scaffold proteins (STE5 and HAM-5) in regulating conserved signal transduction pathways will reveal how selection and evolution has shaped convergent evolution of these processes.
In S. cerevisiae, Gβγ is involved in the recruitment of Ste5 to the plasma membrane upon pheromone exposure, thereby recruiting the Fus3 MAPK cascade to the membrane. Membrane binding of Ste5 likely concentrates the bound MAP kinases spatially, promoting amplification of the signal [59], [60]. In N. crassa, how HAM-5 and MAK-2 are recruited to the membrane is still elusive, but is dissimilar from Ste5 since the Gβγ ortholog in N. crassa is not involved in germling fusion [61], [62]. Other upstream factors shared between S. cerevisiae and N. crassa might regulate the activation of the MAPK pathway. One is STE50, a component in yeast that helps to activate Ste11. In the accompanying article, Dettmann et al., [46] identified a role for N. crassa STE-50 as an activator of NRC-1; Δste-50 mutants were fusion deficient. Three other proteins acting upstream of NRC-1 and STE-50 were also identified: the MAP4 kinase STE-20, the small GTPase RAS-2/SMCO-7 and the capping protein of the adenylate cyclase (AC) complex, CAP-1/NCU08008. Strains carrying a deletion of any of these three genes still showed residual germling fusion, suggesting multiple and redundant inputs into the MAK-2 signal transduction pathway. In S. cerevisiae, Bem1 interacts with Ste20, Ste5 and actin [63]. In N. crassa, strains carrying either a deletion of bem-1, encoding a predicted scaffold for NADPH oxidase (NOX), or its regulator (NOXR), are germling/hyphal fusion deficient [23]. Activated RAC-1 also localizes to CAT tips during chemotropic interactions and may also function as an upstream activator [64]; Δrac-1 mutants are also are germling fusion defective [23].
During chemotropic interactions, HAM-5/MAK-2 complex assembles in puncta at the CAT tip and in the cytoplasm, followed by disassembly of HAM-5/MAK-2 complex from puncta, not just at the CAT tip, but from puncta observed throughout the germling and fusion hyphae, a cycle that repeats itself during chemotropic interactions every ∼8 min. In the absence of HAM-5, localization of MAK-2 kinase complex to puncta was impaired, while in Δmak-2 mutants, HAM-5 puncta were still observed. These observations indicate that MAK-2 kinase activity is essential for disassembly of the HAM-5/MAK-2 complex (Figure 8) during chemotropic interactions. The formation of HAM-5/MAK-2 complexes in puncta was not disrupted in other fusion mutants, such as Δham-7 and Δham-11; cortical localization of the HAM-5/MAK-2 complexes were observed at cell tips, although oscillation was not. These data suggest that the HAM-5/MAK-2/MEK-2/NRC-1 complexes are poised to signal for chemotropic interactions, but that the absence of HAM-7 or HAM-11 disrupts signaling that results in disassembly of the HAM-5/MAK-2 complexes, both within the cell and at the cell cortex.
Few downstream targets of MAK-2 have been identified in filamentous fungi. The PP-1 protein, a transcription factor similar to Ste12 from yeast, is a likely downstream factor that is required for the activation of genes that play a role during the cell fusion and membrane merger [24], [46], [65]. Another target of MAK-2 is MOB-3, a protein of the STRIPAK complex involved in cell fusion that assures correct nuclear localization of MAK-1 [66]. Among the phosphorylated proteins in addition to HAM-5 identified in this study are members of the osmosensing (OS-2, CUT-1 and ASL-1) and cell wall integrity pathways (MAK-1) and other proteins of unknown functions but which are required for fusion (HAM-9 and HAM-11) (Table S1). The accompanying study [46] also identified MAK-1, OS-2 and CUT-1 in MAK-2 complexes via mass spectrometry. Both studies also identified other proteins that interacted with MAK-2/MEK-2/NRC-1 complex [46] or as potential phosphorylation targets of MAK-2 (this study; Table S1), including a glucokinase (NCU00575), SUC (pyruvate decarboxylase), an aminotransferase (NCU03500), a trehalose-phosphatase (NCU05041), CAMK-4 calcium/calmodulin-dependent kinase-4, and four hypothetical proteins (NCU00627, NCU00935, NCU006247 and NCU08330) [46]. The identification of MAK-2 phosphorylation targets provides new clues to the interconnectivity of signaling pathways in N. crassa; these two combined datasets will be a rich resource for further studies on the MAPK pathway function in fungi.
It is challenging to identify kinase targets that are often present in low abundance and have low phosphorylation stoichiometry, from a complex whole cell lysate with limited sample size. The multiplexed iTRAQ quantitation strategy, high specificity of phosphopeptide enrichment and high resolution nano-flow LC separation coupled to MS, as used in this study, have together contributed to the success of identifying and quantifying thousands of phosphopeptides from a small size sample (∼200 µg protein per sample condition). Our phosphoproteomics dataset from 5 hr old germlings provides information on stage-specific phosphorylation events on over ∼1100 proteins (Dataset S1). This dataset can be further compared to a recently published study on ∼3500 phosphorylated proteins identified under hyphal conditions and different carbon sources [40]. Both of these studies provide rich datasets for the filamentous fungal research community to interrogate the identity and function of phosphorylation sites on a large fraction of proteins in the N. crassa proteome. For example, three phosphorylation sites on HAM-5 from germlings were identified from this study, but an additional 13 HAM-5 phosphopeptides were identified in a sample from a 20 hr-old hyphal culture [40]. For chemotropic interactions, further studies on the additional phosphorylation sites and further dissection of the protein domains in the C-terminus of HAM-5 may explain how this scaffold protein itself is recruited to puncta, and may reveal additional binding partners for HAM-5. Such studies will elucidate the molecular mechanism and function of oscillation of the HAM-5/MAK-2/MEK-2/NRC-1 complex during chemotropic interactions. Understanding the molecular basis of germling/hyphal fusion in filamentous fungi provides a window into fungal language and communication and provides a paradigm for self-signaling mechanisms in multicellular eukaryotic species.
Materials and Methods
Molecular techniques and strain construction
Deletion strains used to screen for fusion mutants and strains constructed for this study are listed in Table S2. Strains were grown on Vogel's minimal medium (VMM) [67] (with supplements as required) and were crossed on Westergaard's medium [68]. Transformations and other N. crassa molecular techniques were performed as described [69] or using protocols available at the Neurospora home page at the FGSC (http://www.fgsc.net/Neurospora/NeurosporaProtocolGuide.htm).
To construct the ham-5 alleles, PCR was performed with the restriction enzyme linkers included in the primer region. We amplified ham-5 alleles using primers 1-5: ham-5FXbaI tttttctagaATGTCGGTCCCCGGACACA; ham-5RpacI aaaattaattaaGATCATCTCACTATGATGCAAC; ham-5 WD40onlyR tttttaattaaGTTAGCAGGATGTTGAACGTTG; ham-5RWD40 tttatgcatatttaaaTCATGGTGGCAGCATACAATC; ham-5FWD40 tttatttaaatCCTGCTAACATGTTACCTCC; and cloned the fragments into pCR-Blunt vector (Invitrogen). For constructing the point mutations at the predicted phosphorylation site we used fusion PCR strategy with primers CGTCATCGGGGGCGCGCGGCACCTCATGGCG and GGTGCCGCGCGCCCCCGATGACGCGAAAGTTGT (S → A) and CGTCATCGGGCTCGCGCGGCACCTCATGGCG and GGTGCCGCGCGAGCCCGATGACGCGAAAGTTGT (S → E) together with the primers TTTGATGCATCACAATGCTGACC and TTAAGGGCCGAATTCTTCGC. The mutated constructs were ligated into the NsiI and EcoRI sites of the HAM-5 gene. For constructing the mutation at the predicted docking site, we used primers ggtcgctgacaaactcgaat and tttgccggctgcCTCGGAACTGCGCGCGCGG that has the restriction site NaeI in the linker and TTTCGGCCAGCATCATGAGA and tttagcgctCCTCCAGCACTCGACCTTCGC that has the restriction site AfeI in the linker. The two respective products were digested with NsiI and NaeI and AfeI and Tth111I, respectively and ligated using a three-point ligation in the vector with HAM-5 cut open with NsiI and Tth111I.
We sequenced and digested the constructs from the pCR-Blunt vector with the appropriate restriction enzymes. The different fragments were ligated into plasmid pMF272 (AY598428) [70], [71]. For the MAK-2-mCherry, SO-mCherry, MEK-2-mCherry and NRC-1-mCherry strains, plasmid TSL84C was used. To generate plasmid pTSL84C, sGFP, from plasmid pMF272, was removed by digestion with PacI and EcoRI restriction enzymes and replaced by a version of mCherry that is codon-optimized for N. crassa and includes a C-terminal 6x-His-tag (mCherryNc-6xHis). Plasmid pMFP26 [72] contains untagged codon-optimized mCherryNc and was used as a template for PCR using forward primer OTS177 (AAATTAATTAACGTGAGCAAGGGCGAGGAGGATAAC) and reverse primer OTS202 (AAAGAATTCCTAGTGGTGGTGGTGGTGGTGGCTGCCCTTGTACAGCTCGTCCATGCCGCCG), which contained information for the 6xHis tag and a stop codon. Plasmid derivatives with the tef-1 promoter instead of ccg-1 promoter were obtained by swapping the ccg-1 promoter for the tef-1 promoter using the restriction enzymes NotI and XbaI. A tandem construct of tef-1-ham-5-gfp and tef-1-mak-2-mCherry was created by digesting the tef-1-mak-2-mCherry pMF272 plasmid with restriction enzymes PspOMI and BstBI and the tef-1-ham-5-gfp construct with NotI and BstBI. The latter fragment was ligated into the tef-1-mak-2-mCherry pMF272 plasmid to create tef-1-mak-2-mCherry and tef-1-ham-5-gfp. All constructs were transformed into the WT his-3 strain with selection for His+ prototrophy. Homokaryotic strain was obtained via microconidial purification [73]. A strain bearing cytoplasmic GFP was obtained by transformation of the empty pMF272 plasmid into the WT his-3 strain. All micrographs with HAM-5-GFP are with strains bearing the tef-1-ham-5-gfp constructs unless stated otherwise.
Deletion strains were obtained from the FGSC [74] that were generated as part of the N. crassa functional genomics project [69], [75]. For each deletion strain, both the mating type A and mating type a strains were analyzed, if available.
Phenotypic analyses
To assess the ability of fusion between conidia of a deletion strain as compared to wild type, slant tubes containing the strains were grown for 4–6 days or until significant conidiation occurred. Conidia were harvested by vortexing the slant tube with 2 ml ddH2O and subsequently filtered by pouring over cheesecloth to remove hyphal fragments. Conidia were diluted to a concentration of 3.3×107 conidia/ml. For each sample, 300 µl of spore suspension was spread on a 9 cm solid VMM plate. The plates were dried in a fume hood for 20-30 minutes and incubated for 3–4 hours at 30°C. Squares of 1 cm were excised and observed with a Zeiss Axioskop 2 using a 40× Plan-Neofluor oil immersion objective. The ability of germlings to communicate was determined by evaluating whether germlings displayed homing behavior when germinated conidia were within ∼15 um of each other.
Phosphoproteomics sample preparation
Samples for protein extraction were grown for 4.5 h at 30°C and 200 rpm. Shaking was stopped and samples were grown for 20 minutes longer to encourage cell-cell interaction. The mak-2Q100G strain was treated with 10 µM 1NM-PP1 final concentration in DMSO or with DMSO alone for 10 minutes. Cells were harvested by filtration and frozen before protein extraction. Protein was extracted using the Trizol procedure (according to manufacture's protocols) and was kept in a solution of 6 M guanidine, 50 mM ammonium bicarbonate at pH 7.4. 200 µg for each sample was used for guanidine digestion: 1) pH was adjusted to pH 7.4 with the iTRAQ resuspending buffer (500 mM), 2) 1 hr reduction using 5 mM DTT at 56°C, 3) 1 h alkylation using 10 mM iodoacetamide at RT in the dark, 4) samples were diluted 10 fold with 25 mM NH4HCO3, pH 7.8, 5) 2 mM CaCl2 and trypsin was added at a 1∶50 (trypsin-to-protein) ratio and incubated for 3 hr at 37°C with gentle shaking, 6) trypsin was added again in the same ratio and samples were incubated over-night at 37°C with gentle shaking, 7) subsequently, a standard C18 Solid phase extraction (SPE) was performed with 80% ACN and no TFA, 8) samples were dried in a Speed-Vac and a bicinchoninic acid (BCA) assay was performed.
For 4-plex iTRAQ labeling, 100 µg lyophilized sample per iTRAQ label was used: 1) samples were reconstituted with 30.0 µL of dissolution buffer (500 mM triethylammonium bicarbonate), sonicated and vortexed to resuspend the peptides, 2) sample pH was checked (∼pH 8.5), 3) each vial of iTRAQ reagent (114, 115, 116 and 117) was brought to room temperature and 70 µL of ethanol was added to each iTRAQ reagent vial, 4) samples were vortexed for 1 min and spun down, 5) each labeled reaction mix was added to one separate sample, 6) samples were vortexed, spun down and incubated for 1 hr at RT. Samples were subsequently hydrolyzed by adding 300 uL of 0.05% TFA (3 times the volume), vortexed, spun down and incubated at room temperature for another 30 min, 7) samples were concentrated to 40 µL using a Speed-Vac, 8) samples were pooled into a fresh 2 mL silanized tube and concentrated to ∼100 µL using a Speed-Vac before a desalting (SPE C18) step was performed.
For phosphopeptide enrichment, magnetic Ni-NTA-agarose beads were obtained from Qiagen (Valencia, CA Part N#36111): 1) 50 µL of the 5% suspension metal ion activated NTA was used for 100 ug peptides, 2) beads were first prepared by washing 3× with nano-pure water (800.0 µL of water per 1.0 mL of bead suspension), 3) beads were then treated with 100 mM EDTA, pH 8.0 (800.0 µL of 100 mM EDTA per 1.0 mL of bead suspension) for 30 min with end-over-end rotation, 4) EDTA solution was removed, and beads were washed 3× with nano-pure water (at the same ratio). Subsequently, the beads were treated with 10 mM aqueous metal ion solution (800.0 uL of 100 mM FeCl3 per 1.0 mL of bead suspension) for 30 min with end-over-end rotation, 5) after removing excess metal ions, beads were washed 3× with water (at the same ratio), and resuspended in 1∶1∶1 acetonitrile/methanol/0.01% acetic acid for aliquoting into microcentrifuge tubes, 6) peptide samples were resuspended in 200.0 µl wash/resuspension buffer (80% acetonitrile, 0.1% TFA), 7) beads were washed with of 80% acetonitrile with 0.1% TFA (200 µL of 80% acetonitrile per 50.0 µL of beads) and precipitated using the magnetic stand; the supernatant was discarded. The resuspended samples (100 µg peptides in 200 µl of 80% acetonitrile, 0.1% TFA) were added to the activated beads and incubated for 30 min with end-over-end rotation, 8) beads were precipitated using the magnetic stand and washed for 1 min with 80% acetonitrile, 0.1% TFA. This step was repeated three more times, 9) phosphopeptides were eluted from the beads using an appropriate amount of elution buffer (50.0 µL of the elution buffer per every 50.0 uL of beads/100.0 ug of peptides) after incubating for 5 min, 10) the samples were then acidified to pH 4.0 by concentrating the samples down to 5−10 µl in a Speed-Vac and reconstituted in 30 µL with 0.1% TFA.
Mass-spectrometry based analysis
All peptide samples were analyzed using an automated home-built constant flow nano LC system (Agilent) coupled to an LTQ Orbitrap Velos mass spectrometer (Thermo Fisher Scientific) operating in data-dependent mode [76]. Electrospray emitters were custom made using either 360 µm o.d. ×20 µm i.d. chemically etched fused silica. The nano LC system for phosphoproteomics analysis has an online 4-cm×360 µm o.d. ×150 µm i.d. C18 SPE column (5- µm Jupiter C18, Phenomenex, Torrence, CA) to desalt each phosphopeptide sample (20 µL), which is connected to a home-made 60-cm×360 µm o.d. ×50 µm i.d. capillary column (3- µm Jupiter C18, Phenomenex, Torrence, CA). Mobile phase flow rate was 100 nL/min and consisted of 0.1 M acetic acid in water and 0.1 M acetic acid in 70∶30 (v/v) acetonitrile:water. For each sample, three technical replicates of LC-MS analyses were performed as shown in Figure S1. These included (i) an LC gradient of 300 min with the LTQ Orbitrap Velos mass spectrometer acquiring higher-energy collisional dissociation (HCD) scans; (ii) an LC gradient of 300 min with the LTQ Orbitrap Velos mass spectrometer acquiring alternating collision-induced dissociation (CID), ETD (electron transfer dissociation), and higher-energy collisional dissociation (HCD) scans; (iii) an LC gradient of 180 min with the LTQ Orbitrap Velos mass spectrometer acquiring alternating collision-induced dissociation (CID), ETD (electron transfer dissociation), and higher-energy collisional dissociation (HCD) scans.
Peptide identification and quantification
For peptide identification, MS/MS spectra were searched against a decoy Neurospora protein sequence database using SEQUEST [77]. Search parameters included: trypsin enzyme specificity with a maximum of two missed cleavages, +/- 50 ppm precursor mass tolerance, +/- 0.05 Da product mass tolerance, and carbamidomethylation of cysteines and iTRAQ labeling of lysines and peptide N-termini as fixed modifications. Allowed variable modifications were phosphorylation of serine, threonine or tyrosine residues. MSGF spectra probability value [78] was also calculated for peptides identified from SEQUEST search. Measured mass accuracy and MSGF spectra probability were used to filter identified peptides to <1% false discovery rate (FDR) at spectrum level.
iTRAQ reporter ions were extracted using the MASIC software [79] within 10 ppm mass tolerance of each expected iTRAQ reporter ion from each MS/MS spectrum. The sum of the individual iTRAQ reporter ion values from all MS/MS spectra for a given peptide was used for calculating their relative abundance across different conditions. To correct any systematic error due to pipetting, data were normalized by the median of iTRAQ reporter ion of the individual sample. Fold change of each phosphopeptide was calculated by dividing the data points from the two different conditions (i.e. control and 1NM-PP1-treated) and transformed into Log2 scale. Statistical analyses were performed using Students t-test. Only data with changes exceeding 1.5× greater in control versus 1NM-PP1-treated (p<0.05) were considered differential.
Fluorescence microscopy
The strain used to cross so-gfp into a Δham-5 strain was AF-SoT8 and the strain used to cross mak-2-gfp into a Δham-5 strain was AF-M512 (Table S2).
Oscillation studies performed with HAM-5-GFP and mCherry tagged strains were prepared as described above with modifications from [20]. Images were taken using a Leica SD6000 microscope with a 100×1.4 NA oil-immersion objective equipped with a Yokogawa CSU-X1 spinning disk head and a 488-nm or 561-nm laser controlled by Metamorph software. Multiple pairs of interacting germlings were analyzed per experiment and representative pairs are shown for each strain. The ImageJ software was used for image analysis.
Immunoprecipitations and western blotting
Harvested conidia (1×106/ml) were inoculated in 100 ml VMM in flasks and incubated for 2.5 hrs at 30°C with shaking at 200 rpm, then an additional 2.5 hrs at 30°C without shaking. Germlings from 3 flasks were harvested by vacuum filtration over a nitrocellulose membrane and frozen in liquid nitrogen. Protein extraction from ground mycelium was performed using 1 ml lysis buffer described in [8] containing complete protease inhibitors, phosphatase inhibitor and Triton X-100. 20 µl supernatant was used for western blotting and the remaining fraction was used for immunoprecipitation using Protein G Dynabeads (Invitrogen), according to manufacturer's instructions, with the following exceptions: mouse or rabbit anti-GFP antibody (Roche or Life Technologies, respectively) or rabbit anti-mCherry antibody (Bio-vision) was covalently bound to the beads using BS3 (Sulfo-DSS, Fisher scientific) or DMP (dimethylpimelimidate) according to manufacturer's instructions.
Supernatant samples were incubated with the beads overnight at 4°C. Beads were washed with standard PBS for three times before protein was removed from the beads by heating at 70°C for 10 min in 1× loading buffer, and samples were run on a 4–12% Nu-Page Bis-Tris GelGel (NOVEX, Life Technologies). Protein gels were subjected to Western blot analysis using standard methods. Samples for the MAK-1 and MAK-2 phosphorylation western blots were treated similarly, except, after protein extraction with 1 ml lysis buffer, 25 µl of protein sample was directly loaded on a 7% NuPage Bis-Tris GelGel (NOVEX, Life Technologies) protein gel. Gels were subjected to Western blot analysis using standard methods and detection of phosphorylated MAK-1 and MAK-2 was carried out using anti-phospho p44/42 MAP kinase antibodies (1∶3000 dilution) (PhosphoPlus antibody kit; Cell Signaling Technology) as described [8]. Detection of phosphorylated HAM-5-GFP was performed using anti-phosphothreonine-proline/phosphoserine-proline antibodies (Abcam).
Supporting Information
Zdroje
1. AguilarPS, BayliesMK, FleissnerA, HelmingL, InoueN, et al. (2013) Genetic basis of cell-cell fusion mechanisms. Trends Genet 29: 427–437.
2. ChenEH, OlsonEN (2005) Unveiling the mechanisms of cell-cell fusion. Science 308: 369–373.
3. Fricker MD, Boddy L, Bebber D (2007) Network organization of filamentous fungi. In: Howard, R and Gow, N. A R, Eds. Biology of the Fungal Cell. Berlin: Springer-Verlag. pp.309-330.
4. Rayner ADM (1996) Interconnectedness and individualism in fungal mycelia; Sutton BC, editor. Cambridge: Cambridge University Press. pp.193-232.
5. GlassNL, RasmussenC, RocaMG, ReadND (2004) Hyphal homing, fusion and mycelial interconnectedness. Trends Microbiol 12: 135–141.
6. Roca GM, Read ND, Wheals AE (2005) Conidial anastomosis tubes in filamentous fungi. FEMS Microbio Letts 249: : 191-198.
7. Read ND, Fleißner A, Roca MG, Glass NL (2010) Hyphal Fusion. In: Borkovich, K.A and Ebbole, D., Eds. Cellular and Molecular Biology of Filamentous Fungi. Washington, D.C. ASM Press. pp.260-273.
8. PandeyA, RocaMG, ReadND, GlassNL (2004) Role of a mitogen-activated protein kinase pathway during conidial germination and hyphal fusion in Neurospora crassa. Eukaryot Cell 3: 348–358.
9. SimoninA, Palma-GuerreroJ, FrickerM, GlassNL (2012) Physiological significance of network organization in fungi. Eukaryot Cell 11: 1345–1352.
10. RoperM, SimoninA, HickeyPC, LeederA, GlassNL (2013) Nuclear dynamics in a fungal chimera. Proc Natl Acad Sci U S 110: 12875–12880.
11. RichardF, GlassNL, PringleA (2012) Cooperation among germinating spores facilitates the growth of the fungus, Neurospora crassa. Biol Lett 8: 419–422.
12. CharltonND, ShojiJY, GhimireSR, NakashimaJ, CravenKD (2012) Deletion of the fungal gene soft disrupts mutualistic symbiosis between the grass endophyte Epichloe festucae and the host plant. Eukaryot Cell 11: 1463–1471.
13. TanakaA, CartwrightGM, SaikiaS, KayanoY, TakemotoD, et al. (2013) ProA, a transcriptional regulator of fungal fruiting body development, regulates leaf hyphal network development in the Epichloe festucae-Lolium perenne symbiosis. Mol Microbiol 90: 551–568.
14. RocaM, ArltJ, JeffreeC, ReadN (2005) Cell biology of conidial anastomosis tubes in Neurospora crassa. Eukaryot Cell 4: 911–919.
15. ReadND, LichiusA, ShojiJY, GoryachevAB (2009) Self-signalling and self-fusion in filamentous fungi. Curr Opin Microbiol 12: 608–615.
16. IshikawaFH, SouzaEA, ShojiJY, ConnollyL, FreitagM, et al. (2012) Heterokaryon incompatibility is suppressed following conidial anastomosis tube fusion in a fungal plant pathogen. PLoS One 7: e31175.
17. FleissnerA, SimoninAR, GlassNL (2008) Cell fusion in the filamentous fungus, Neurospora crassa. Methods Mol Biol 475: 21–38.
18. DettmannA, IllgenJ, MarzS, SchurgT, FleissnerA, et al. (2012) The NDR kinase scaffold HYM1/MO25 is essential for MAK-2 map kinase signaling in Neurospora crassa. PLoS Genet 8: e1002950.
19. FleissnerA, LeederAC, RocaMG, ReadND, GlassNL (2009) Oscillatory recruitment of signaling proteins to cell tips promotes coordinated behavior during cell fusion. Proc Natl Acad Sci USA 106: 19387–19392.
20. FleissnerA, GlassNL (2007) SO, a protein involved in hyphal fusion in Neurospora crassa, localizes to septal plugs. Eukaryot Cell 6: 84–94.
21. GoryachevAB, LichiusA, WrightGD, ReadND (2012) Excitable behavior can explain the "ping-pong" mode of communication between cells using the same chemoattractant. Bioessays 34: 259–266.
22. TeichertI, SteffensEK, SchnabN, FranzelB, KrispC, et al. (2014) PRO40 is a scaffold protein of the cell wall integrity pathway, linking the MAP kinase module to the upstream activator protein kinase C. PLoS Genet 10: e1004582.
23. FuC, IyerP, HerkalA, AbdullahJ, StoutA, et al. (2011) Identification and characterization of genes required for cell-to-cell fusion in Neurospora crassa. Eukaryot Cell 10: 1100–1109.
24. LeederAC, JonkersW, LiJ, GlassNL (2013) Germination and early colony establishment in Neurospora crassa requires a MAP kinase regulatory network. Genetics 195: 883–898.
25. KosakoH, NaganoK (2011) Quantitative phosphoproteomics strategies for understanding protein kinase-mediated signal transduction pathways. Expert Rev Proteomics 8: 81–94.
26. WangY, DingSJ, WangW, JacobsJM, QianWJ, et al. (2007) Profiling signaling polarity in chemotactic cells. Proc Natl Acad Sci USA 104: 8328–8333.
27. KirkpatrickDS, BustosDJ, DoganT, ChanJ, PhuL, et al. (2013) Phosphoproteomic characterization of DNA damage response in melanoma cells following MEK/PI3K dual inhibition. Proc Natl Acad Sci USA 110: 19426–19431.
28. MertinsP, YangF, LiuT, ManiDR, PetyukVA, et al. (2014) Ischemia in tumors induces early and sustained phosphorylation changes in stress kinase pathways but does not affect global protein levels. Mol Cell Proteomics 13: 1690–1704.
29. AldabbousMS, RocaMG, StoutA, HuangIC, ReadND, et al. (2010) The ham-5, rcm-1 and rco-1 genes regulate hyphal fusion in Neurospora crassa. Microbiol 156: 2621–2629.
30. RueppA, ZollnerA, MaierD, AlbermannK, HaniJ, et al. (2004) The FunCat, a functional annotation scheme for systematic classification of proteins from whole genomes. Nucl Acids Res 32: 5539–5545.
31. MaerzS, ZivC, VogtN, HelmstaedtK, CohenN, et al. (2008) The nuclear Dbf2-related kinase COT1 and the mitogen-activated protein kinases MAK-1 and MAK-2 genetically interact to regulate filamentous growth, hyphal fusion and sexual development in Neurospora crassa. Genetics 179: 1313–1325.
32. ZhangY, LammR, PillonelC, LamS, XuJR (2002) Osmoregulation and fungicide resistance: the Neurospora crassa os-2 gene encodes a HOG1 mitogen-activated protein kinase homologue. Appl Environ Microbiol 68: 532–538.
33. YoussarL, SchmidhauserTJ, AvalosJ (2005) The Neurospora crassa gene responsible for the cut and ovc phenotypes encodes a protein of the haloacid dehalogenase family. Mol Microbiol 55: 828–838.
34. FujimuraM, OchiaiN, OshimaM, MotoyamaT, IchiishiA, et al. (2003) Putative homologs of SSK22 MAPKK kinase and PBS2 MAPK kinase of Saccharomyces cerevisiae encoded by os-4 and os-5 genes for osmotic sensitivity and fungicide resistance in Neurospora crassa. Biosci Biotechnol Biochem 67: 186–191.
35. YamashitaK, ShiozawaA, WatanabeS, FukumoriF, KimuraM, et al. (2008) ATF-1 transcription factor regulates the expression of ccg-1 and cat-1 genes in response to fludioxonil under OS-2 MAP kinase in Neurospora crassa. Fungal Genet Biol 45: 1562–1569.
36. MokJ, KimPM, LamHY, PiccirilloS, ZhouX, et al. (2010) Deciphering protein kinase specificity through large-scale analysis of yeast phosphorylation site motifs. Science Signaling 3: ra12.
37. ParnellSC, MarottiLAJr, KiangL, TorresMP, BorchersCH, et al. (2005) Phosphorylation of the RGS protein Sst2 by the MAP kinase Fus3 and use of Sst2 as a model to analyze determinants of substrate sequence specificity. Biochem 44: 8159–8166.
38. LiD, BobrowiczP, WilkinsonHH, EbboleDJ (2005) A mitogen-activated protein kinase pathway essential for mating and contributing to vegetative growth in Neurospora crassa. Genetics 170: 1091–1104.
39. RemenyiA, GoodMC, BhattacharyyaRP, LimWA (2005) The role of docking interactions in mediating signaling input, output, and discrimination in the yeast MAPK network. Mol Cell 20: 951–962.
40. XiongY, CoradettiST, LiX, GritsenkoMA, ClaussT, et al. (2014) The proteome and phosphoproteome of Neurospora crassa in response to cellulose, sucrose and carbon starvation. Fungal Genet Biol S1087-1845(14)00083-8.
41. XuC, MinJ (2011) Structure and function of WD40 domain proteins. Protein & Cell 2: 202–214.
42. GoodMC, ZalatanJG, LimWA (2011) Scaffold proteins: hubs for controlling the flow of cellular information. Science 332: 680–686.
43. MaddiA, DettmanA, FuC, SeilerS, FreeSJ (2012) WSC-1 and HAM-7 are MAK-1 MAP kinase pathway sensors required for cell wall integrity and hyphal fusion in Neurospora crassa. PLoS ONE 7: e42374.
44. LeederAC, Palma-GuerreroJ, GlassNL (2011) The social network: deciphering fungal language. Nat Rev Microbiol 9: 440–451.
45. HickeyPC, JacobsonD, ReadND, Louise GlassNL (2002) Live-cell imaging of vegetative hyphal fusion in Neurospora crassa. Fungal Genet Biol 37: 109–119.
46. DettmannA, HeiligY, ValeriusO, LudwigS, SeilerS (2014) Fungal communication requires the MAK-2 pathway elements STE-20, RAS-2, the NRC-1 adapter STE-50 and the MAP kinase scaffold HAM-5. PLoS Genet 10: e1004762.
47. HouZ, XueC, PengY, KatanT, KistlerHC, et al. (2002) A mitogen-activated protein kinase gene (MGV1) in Fusarium graminearum is required for female fertility, heterokaryon formation, and plant infection. Mol Plant Microbe Interact 15: 1119–1127.
48. Di PietroA, Garcia-MacEiraFI, MegleczE, RonceroMIG (2001) A MAP kinase of the vascular wilt fungus Fusarium oxysporum is essential for root penetration and pathogenesis. Mol Microbiol 39: 1140–1152.
49. BayramO, BayramOS, AhmedYL, MaruyamaJ, ValeriusO, et al. (2012) The Aspergillus nidulans MAPK module AnSte11-Ste50-Ste7-Fus3 controls development and secondary metabolism. PLoS Genet 8: e1002816.
50. Ruiz-RoldanMC, KohliM, RonceroMI, PhilippsenP, Di PietroA, et al. (2010) Nuclear dynamics during germination, conidiation, and hyphal fusion of Fusarium oxysporum. Eukaryot Cell 9: 1216–1224.
51. Jamet-ViernyC, DebuchyR, PrigentM, SilarP (2007) IDC1, a pezizomycotina-specific gene that belongs to the PaMpk1 MAP kinase transduction cascade of the filamentous fungus Podospora anserina. Fungal Genet Biol 44: 1219–1230.
52. XuJR, StaigerCJ, HamerJE (1998) Inactivation of the mitogen-activated protein kinase Mps1 from the rice blast fungus prevents penetration of host cells but allows activation of plant defense responses. Proc Natl Acad Sci USA 95: 12713–12718.
53. LevS, SharonA, HadarR, MaH, HorwitzBA (1999) A mitogen-activated protein kinase of the corn leaf pathogen Cochliobolus heterostrophus is involved in conidiation, appressorium formation, and pathogenicity: diverse roles for mitogen-activated protein kinase homologs in foliar pathogens. Proc Natl Acad Sci USA 96: 13542–13547.
54. ZhengL, CampbellM, MurphyJ, LamS, XuJR (2000) The BMP1 gene is essential for pathogenicity in the gray mold fungus Botrytis cinerea. Mol Plant-Microbe Interact 13: 724–732.
55. KojimaK, KikuchiT, TakanoY, OshiroE, OkunoT (2002) The mitogen-activated protein kinase gene MAF1 is essential for the early differentiation phase of appressorium formation in Colletotrichum lagenarium. Mol Plant Microbe Interact 15: 1268–1276.
56. LengelerKB, DavidsonRC, D'SouzaC, HarashimaT, ShenW-C, et al. (2000) Signal transduction cascades regulating fungal development and virulence. Microbiol Mol Biol Rev 64: 746–785.
57. XuJR, HamerJE (1996) MAP kinase and cAMP signaling regulate infection structure formation and pathogenic growth in the rice blast fungus Magnaporthe grisea. Genes Devel 10: 2696–2706.
58. CoyleSM, FloresJ, LimWA (2013) Exploitation of latent allostery enables the evolution of new modes of MAP kinase regulation. Cell 154: 875–887.
59. BhattacharyyaRP, RemenyiA, GoodMC, BashorCJ, FalickAM, et al. (2006) The Ste5 scaffold allosterically modulates signaling output of the yeast mating pathway. Science 311: 822–826.
60. StrickfadenSC, WintersMJ, Ben-AriG, LamsonRE, TyersM, et al. (2007) A mechanism for cell-cycle regulation of MAP kinase signaling in a yeast differentiation pathway. Cell 128: 519–531.
61. WonS, MichkovAV, KrystofovaS, GarudAV, BorkovichKA (2012) Genetic and physical interactions between Galpha subunits and components of the Gbetagamma dimer of heterotrimeric G proteins in Neurospora crassa. Eukaryot Cell 11: 1239–1248.
62. EatonCJ, CabreraIE, ServinJA, WrightSJ, CoxMP, et al. (2012) The guanine nucleotide exchange factor RIC8 regulates conidial germination through Galpha proteins in Neurospora crassa. PLoS One 7: e48026.
63. LeeuwT, Fourest-LieuvinA, WuC, ChenevertJ, ClarkK, et al. (1995) Pheromone response in yeast: association of Bem1p with proteins of the MAP kinase cascade and actin. Science 270: 1210–1213.
64. LichiusA, GoryachevAB, FrickerMD, ObaraB, Castro-LongoriaE, et al. (2014) CDC-42 and RAC-1 regulate opposite chemotropisms in Neurospora crassa. J Cell Sci 127: 1953–1965.
65. Palma-GuerreroJ, LeederAC, WelchJ, GlassNL (2014) Identification and characterization of LFD1, a novel protein involved in membrane merger during cell fusion in Neurospora crassa. Mol Microbiol 92: 164–182.
66. DettmannA, HeiligY, LudwigS, SchmittK, IllgenJ, et al. (2013) HAM-2 and HAM-3 are central for the assembly of the Neurospora STRIPAK complex at the nuclear envelope and regulate nuclear accumulation of the MAP kinase MAK-1 in a MAK-2-dependent manner. Mol Microbiol 90: 796–812.
67. VogelHJ (1956) A convenient growth medium for Neurospora. Microbiol Genet Bull 13: 42–46.
68. WestergaardM, MitchellHK (1947) Neurospora V. A synthetic medium favoring sexual reproduction. Amer J Bot 34: 573–577.
69. ColotHV, ParkG, TurnerGE, RingelbergC, CrewCM, et al. (2006) A high-throughput gene knockout procedure for Neurospora reveals functions for multiple transcription factors. Proc Natl Acad Sci USA 103: 10352–10357.
70. MargolinBS, FreitagM, SelkerEU (1997) Improved plasmids for gene targeting at the his-3 locus of Neurospora crassa by electroporation. Fungal Genet Newslett 44: 34–36.
71. FreitagM, HickeyPC, RajuNB, SelkerEU, ReadND (2004) GFP as a tool to analyze the organization, dynamics and function of nuclei and microtubules in Neurospora crassa. Fungal Genet Biol 41: 897–910.
72. Castro-LongoriaE, FerryM, Bartnicki-GarciaS, HastyJ, BrodyS (2010) Circadian rhythms in Neurospora crassa: dynamics of the clock component frequency visualized using a fluorescent reporter. Fungal Genet Biol 47: 332–341.
73. PanditA, MaheshwariR (1994) A simple method of obtaining pure microconidia in Neurospora crassa. Fungal Genet Newsl 41: 64–65.
74. McCluskeyK (2003) The Fungal Genetics Stock Center: from molds to molecules. Adv Appl Microbiol 52: 245–262.
75. DunlapJC, BorkovichKA, HennMR, TurnerGE, SachsMS, et al. (2007) Enabling a community to dissect an organism: Overview of the Neurospora functional genomics project. Adv Genet 57: 49–96.
76. NguyenTH, BrechenmacherL, AldrichJT, ClaussTR, GritsenkoMA, et al. (2012) Quantitative phosphoproteomic analysis of soybean root hairs inoculated with Bradyrhizobium japonicum. Mol Cell Proteomics 11: 1140–1155.
77. EngJK, McCormackAL, YatesJR (1994) An approach to correlate tandem mass spectral data of peptides with amino acid sequences in a protein database. J Amer Soc Mass Spectrometry 5: 976–989.
78. KimS, GuptaN, PevznerPA (2008) Spectral probabilities and generating functions of tandem mass spectra: a strike against decoy databases. J Proteome Res 7: 3354–3363.
79. MonroeME, ShawJL, DalyDS, AdkinsJN, SmithRD (2008) MASIC: a software program for fast quantitation and flexible visualization of chromatographic profiles from detected LC-MS(/MS) features. Comp Biol Chem 32: 215–217.
Štítky
Genetika Reprodukční medicínaČlánek vyšel v časopise
PLOS Genetics
2014 Číslo 11
Nejčtenější v tomto čísle
- An RNA-Seq Screen of the Antenna Identifies a Transporter Necessary for Ammonia Detection
- Systematic Comparison of the Effects of Alpha-synuclein Mutations on Its Oligomerization and Aggregation
- Functional Diversity of Carbohydrate-Active Enzymes Enabling a Bacterium to Ferment Plant Biomass
- Regularized Machine Learning in the Genetic Prediction of Complex Traits