FEZ2 Has Acquired Additional Protein Interaction Partners Relative to FEZ1: Functional and Evolutionary Implications

Background The FEZ (fasciculation and elongation protein zeta) family designation was purposed by Bloom and Horvitz by genetic analysis of C. elegans unc-76. Similar human sequences were identified in the expressed sequence tag database as FEZ1 and FEZ2. The unc-76 function is necessary for normal axon fasciculation and is required for axon-axon interactions. Indeed, the loss of UNC-76 function results in defects in axonal transport. The human FEZ1 protein has been shown to rescue defects caused by unc-76 mutations in nematodes, indicating that both UNC-76 and FEZ1 are evolutionarily conserved in their function. Until today, little is known about FEZ2 protein function. Methodology/Principal Findings Using the yeast two-hybrid system we demonstrate here conserved evolutionary features among orthologs and non-conserved features between paralogs of the FEZ family of proteins, by comparing the interactome profiles of the C-terminals of human FEZ1, FEZ2 and UNC-76 from C. elegans. Furthermore, we correlate our data with an analysis of the molecular evolution of the FEZ protein family in the animal kingdom. Conclusions/Significance We found that FEZ2 interacted with 59 proteins and that of these only 40 interacted with FEZ1. Of the 40 FEZ1 interacting proteins, 36 (90%), also interacted with UNC-76 and none of the 19 FEZ2 specific proteins interacted with FEZ1 or UNC-76. This together with the duplication of unc-76 gene in the ancestral line of chordates suggests that FEZ2 is in the process of acquiring new additional functions. The results provide also an explanation for the dramatic difference between C. elegans and D. melanogaster unc-76 mutants on one hand, which cause serious defects in the nervous system, and the mouse FEZ1 -/- knockout mice on the other, which show no morphological and no strong behavioural phenotype. Likely, the ubiquitously expressed FEZ2 can completely compensate the lack of neuronal FEZ1, since it can interact with all FEZ1 interacting proteins and additional 19 proteins.


Introduction
The FEZ (fasciculation and elongation protein zeta) family designation was first purposed by Bloom and colleagues in 1997 by genetic analysis of C. elegans unc-76. The 376 and 385 amino acid containing isoforms of the UNC-76 protein arise by alternative splicing of the same gene and showed no strong similarity to any previously characterized proteins, but similar human sequences were identified in the expressed sequence tag database, dbEST [1,2] and Bloom & Horvitz named the two identified human genes FEZ1 and FEZ2.
The unc-76 function is necessary for normal fascicle structure and is required specifically for axon-axon interactions in C. elegans. The human FEZ1 gene was able to restore partial locomotion and axonal fasciculation in the C. elegans unc-76 mutants in germ-line transformation experiments, indicating that both the function and the structure of the FEZ proteins have been conserved in evolution [1]. Loss of Drosophila Unc-76 function results in locomotion and axonal transport defects reminiscent of phenotypes observed in kinesin mutants and thereby suggesting that UNC-76 is required for kinesin-dependent axonal transport [3]. The FEZ1-deficient mice however, did not exhibit any obvious abnormal brain architecture, although they manifest slight behavioral abnormalities, including a hyperlocomotion phenotype and enhanced responsiveness to psychostimulants [4].
The homologous proteins UNC-76, FEZ1 and FEZ2, share a conserved predicted coiled-coil region at their C-terminal regions [1,5,6]. Coiled-coils are autonomous folding units consisting of two to four a-helices that wrap around each other with a slight lefthanded super-helical twist [7] and thereby mediate sub-unit oligomerization [8] as well as protein-protein interactions [9]. Up to now 4 interactions were reported for C. elegans UNC-76, 45 interactions for human FEZ1 and 17 interactions for human FEZ2 according to BIND, HPRD, BioGRID data and previous research published by our group [10]. All of these interactions occur at the C-terminus with the predicted coiled-coil region. Some of these interactions implicate FEZ1 protein in the context of schizophrenia (DISC1) [11,12], aspects of viral infection such as post-entry block of retroviral infection [13] and inhibition of the release of progeny virions (agnoprotein of the human polyomavirus JC virus) [14], and naturally in neuritogenesis (NBR1, PKCf, DISC1) [5,11,15]. Furthermore, the association between FEZ1 and E4B was enhanced by co-expression of a constitutively active form of protein kinase Cf, and phosphorylation of FEZ1 by this kinase and its subsequent ubiquitylation by E4B resulted in neurite extension in PC12 cells [4,16]. FEZ1 also contributes to the polarization of hippocampal neurons by controlling mitochondrial motility [4,17].
Here, we report that nucleotide sequences belonging to the FEZ family can be found in many eukaryotic species, but that most of the putative homologs are not only functionally uncharacterized but even remain unrecognized as FEZ family members in the databases. Our phylogenetic analysis of the FEZ family of proteins indicates that the ancestral gene of the FEZ family is found only in the animal kingdom. The gene duplication of unc-76 to FEZ1 and FEZ2 is predicted to have occurred during the two rounds of whole-genome duplication in the chordate ancestral line, after the divergence of cephalochordates but before the splitting of the teleosts and tetrapods taxa.
The most conserved region in FEZ proteins is its C-terminus, and we and others showed previously that this region is involved in the association with other proteins [10,18]. Here, yeast two-hybrid assays with C. elegans UNC-76 (248-372), human FEZ1 (221-392) and human FEZ2 (207-353) as baits and human prey clones from cDNA libraries demonstrated that the pattern of protein-protein interactions (PPIs) is highly conserved between UNC-76 and FEZ1. These orthologues did however not interact with 19 proteins that interacted specifically with FEZ2, although FEZ2 interacted with all 40 proteins that interact with FEZ1/UNC-76. In summary, we found largely overlapping PPI patterns for FEZ1 and UNC-76 and extended PPIs for FEZ2. Interestingly, our data provide an explanation for the ability of FEZ1 to rescue the defects caused by unc-76 mutations in nematodes and to the lack of strong defects in FEZ1 deficient mice, since the continuous presence of the ubiquitous expressed paralogue FEZ2 may be able to compensate lack of FEZ1 because its PPI pattern contains all the interactions found for FEZ1.

Protein sequence analysis and multiple sequence alignments
Protein Psi-Blast [19] searches with the full length C. elegans UNC-76 sequence were performed at the NCBI Web site http://www. ncbi.nlm.nih.gov/BLAST/ using the non-redundant protein sequence database available at August 16, 2009. After six rounds of iteration, all UNC-76, FEZ1 and FEZ2 orthologs with an E-value of 0.005 or below were selected and all redundant sequences were excluded (except to B. floridae). All sequences were collected in FASTA format for further analysis as shown in Table S1. The identification and posterior naming of the protein sequences as either UNC-76, FEZ1 or FEZ2 is based on the phylogenetical analyses shown in Figure 1 and Figure S1. The sequences were aligned using the ClustalW2 Web Site http://www.ebi.ac.uk/ Tools/clustalw2/index.html at default settings. The alignments were then shaded using the multiple sequence alignment editor GENEDOC http://www.nrbsc.org/gfx/genedoc/index.html. The identity and similarity were calculated using default settings by NPS@: Network Protein Sequence Analysis trough the alignments of human proteins FEZ1 and FEZ2 and C. elegans UNC-76 protein (http://npsa-pbil.ibcp.fr/cgi-bin/align_clustalw.pl) [20].

Phylogenetical analyses
PHYLIP version 3.5c [21] was used for the phylogenetical analyses at the web site http://mobyle.pasteur.fr/cgi-bin/portal.py?jobs = http://mobyle.pasteur.fr/data/jobs/protpars. Parsimony analyses were performed using the protein alignment as input. Bootstrap values were obtained by using SEQBOOT and creating 1000 Bootstrap data sets. Analysis was then performed using PROTPARS with standard parameters. The ''M'' option was invoked for the analysis of the multiple data sets generated by SEQBOOT.
Yeast two-hybrid screen and DNA sequence analysis Using FEZ1 as a bait in a screening, Assmann et al. (2006) [10] identified 16 proteins interacting with FEZ1 and all of them interacted with FEZ2 protein ( Fig. 2A). In this work, yeast twohybrid screens [22] of human fetal brain and bone marrow cDNA libraries (Clontech) were performed as described previously by using the yeast strain L40 (trp1-901, his3D200, leu2-3, ade2 LYS2::(lex-Aop)4-HIS3 URA3::(lexAop)8-lac GAL4) and human FEZ2 (207-353) as a bait in fusion with the LexA protein as encoded in the recombinant vector pBTM116 [23]. This fragment of FEZ2 does not auto-activate the yeast reporter genes (as FEZ proteins full length), neither did so the C-terminal constructs of FEZ1 and UNC-76 used in one-to-one confirmation. The autonomous activation test for HIS3 was performed in minimal medium plates without tryptophane and histidine but containing 10, 20, 30, or 50 mM of 3amino-1,2,4-triazole (3-AT). Furthermore, the autonomous activation of LacZ was measured by the beta-galactosidase filter assay as described below. Yeast cells were transformed according to the protocols supplied by the cDNA library manufacturer (Clontech). The screening with FEZ2 (207-353) as bait was performed in minimal medium plates without tryptophan, leucine, and histidine and with addition of 10 mM of 3-AT to repress unspecific background growth.Recombinant pACT2 plasmids of positive clones were isolated and their insert DNAs were sequenced with a DNA sequencer model 377S (Applied Biosystems, Foster City, CA). The obtained DNA sequence data were compared with sequences in the NCBI data bank using the BLASTX 2.2.12 program [19]. The corresponding Accession numbers of the DNA sequences identified are given in the Table S1. As no new sequences have been obtained no new sequence data have been deposited in the GenBank. After identification, one representative clone of each prey protein was selected to yeast two-hybrid one-to-one confirmation.
Yeast Two-hybrid assays one-to-one confirmation with the homologues FEZ1, FEZ2 and UNC-76 Yeast two-hybrid assays were performed with human FEZ1 (221-392) and C. elegans UNC-76 (248-372) as additional baits  to LexA in vector pBTM116, as described above and previously [10]. The recombinant pACT2 plasmids of positives clones isolated from the library screens with bait FEZ2 (207-353) were used as preys and the interactions with FEZ1 and UNC-76 (and control FEZ2) were analyzed by growth in minimal medium plates and activation of the LacZ using the standard betagalactosidase filter assay described in the following section (Fig. 2B).
Assay for beta-galactosidase activity in yeast cells beta-Galactosidase activity in yeast cells was determined by the filter assay method. Yeast transformants (Leu + , Trp + , and His + ) were transferred onto nylon membranes, permeabilized in liquid nitrogen, and placed on Whatman 3 MM paper previously soaked in Z buffer (60 mM Na 2 HPO 4 , 40 mM NaH 2 PO 4 , 10 mM MgCl 2 , 50 mM 2-mercaptoethanol, pH 7.0) containing 1 mg/ml 5-bromo-4-chloro-3-indolyl-b-D-galactoside (X-gal). After incubation at 37uC for 1 h, the yeast cells forming dark blue colonies were taken from replica plates for further analysis (primary screen).

Results
Identification of members of the FEZ family and of conserved regions in their amino acid sequences A database of sequences judged to be members of the FEZ protein family was compiled (Table S1, see Material and Methods section). Altogether 47 members for the FEZ family were identified by Psi-Blast searches with the C. elegans UNC-76 sequence in the non-redundant protein sequence database. FEZ1 family sequences were found in a variety of species ranging from Nematodes, via Arthropods to Mammals but not in Plants, Fungi or Protists. Multiple Clustal W alignments of sequences of the FEZ family identified a highly conserved region both at the Nterminus (ca. aa 160-180 in FEZ1) and at the C-terminus ( Figure  S1). The latter consists mainly in coiled-coil regions [10,24,25] located at the C-terminus which mediate the majority of the protein-protein interactions (ca. aa 230-307 in FEZ1) [10]. The alignment revealed, that the nuclear localization signal (NLS) (KKRRK, ca. aa 290-294 in FEZ1) [26], which consists of basic amino acids, is quite conserved in all FEZ1 genes that resulted Ciona intestinalis, Da -Drosophila ananassae, Dg -Drosophila grimshawi, Dm -Drosophila melanogaster, Dmj -Drosophila mojavensis, Dp -Drosophila persimilis, Dpo -Drosophila pseudoobscura, Dr -Danio rerio, Ds -Drosophila sechellia, Dv -Drosophila virilis, Dw -Drosophila willistoni, Dy -Drosophila yakuba, Ec -Equus caballus, Gg -Gallus gallus, Hs -Homo sapiens, Md -Monodelphis domestica, Mf -Macaca fascicularis, Mm -Mus musculus, Mmt -Macaca mulatta, Nv -Nasonia vitripennis, Ph -Pediculus humanus, Pt -Pan troglodytes, Rn -Rattus norvegicus, Ss -Salmo salar, Tc -Tribolium castaneum, Tg -Taeniopygia guttata, Tn -Tetraodon nigroviridis, Xl -Xenopus laevis, Xt -Xenopus tropicalis. doi:10.1371/journal.pone.0017426.g001 from the putative unc-76 gene duplication (see details below), suggesting that new nuclear functions may have been acquired in FEZ1 afterwards. Although there is no striking NLS in FEZ2 sequences, they contain a conserved polybasic region (e.g. KKKKK) with conserved amino acid substitutions in the FEZ1 related NLS. The NLS is less well conserved in all UNC-76 sequences.

Phylogenetic analysis of the FEZ gene duplication
A phylogenetic tree ( Fig. 1) was generated from the Clustal W alignment ( Figure S1) of the all identified FEZ protein sequences (Table S1) using the parsimony method. Bootstrap options to resampling methods were defined to contain 1000 replications. The presence of FEZ1 and FEZ2 in all analyzed Actinopterygii (ray-finned fishes) genomes suggests that the FEZ1 and FEZ2 genes originated from the ancestral unc-76 by gene or whole genome duplications before the radiation of the Actinopterygii but after the divergence of the amphioxus branch and hence concomitant origin of the chordata. In fact, according to an amphioxus-human synteny analysis by Putnam [27] two rounds of whole-genome duplication seem to have occurred in the chordate stemline, after the divergence of the cephalochordates but before the split into teleosts and tetrapods (Fig. 1). Hence, the unc-76 gene seems to have duplicated in the chordate stemline into FEZ1 and FEZ2 genes.

Identification of proteins that interact with FEZ2
So far, only interactions for C-terminal of FEZ proteins are described, as by our screenings using C termini of FEZ proteins (due to auto-activation of full length and N-terminal constructs of the FEZ proteins as baits in these assays), as well as by other publications describing FEZ proteins recovered as prey in other screens (in this case, the interaction region was mapped to be the C-terminal region) [5,14]. In a previous study, we had identified proteins that interact with human FEZ1 using its C-terminal region (221-392) as a bait in a yeast two-hybrid screen of a human fetal brain cDNA library [10]. We further found that all of the 16 interacting proteins identified were also able to interact with FEZ2 (207-353). The interacting proteins are functionally involved in transcriptional regulation (6 proteins), neuronal cell development (2 proteins), intracellular transport processes (3 proteins), apoptosis (2 proteins), or were of unknown function (3 proteins). Although these data suggested that FEZ1 and FEZ2 have largely overlapping protein interaction profiles we wanted to deepen our insights by performing another two-hybrid screen with the human paralogue FEZ2 as well as by introducing an evolutionary component by comparing the interaction profiles of both proteins with the orthologous C. elegans protein UNC-76, which may have preserved ancient features.
To this end we started by employing the yeast two-hybrid system [23] to screen humans fetal brain and bone marrow cDNA libraries. The truncated FEZ2 (207-353), which showed no autoactivation (see Materials and Methods), was used as a bait (Fig. 2). All grown colonies that showed a strong blue color in the subsequent b-galactosidase filter assay had their plasmid DNA extracted and sequenced. A total of 166 plasmid DNAs (69 from human bone marrow and 97 from human fetal brain cDNA libraries) from clones positive for both HIS3 and LacZ reporters were sequenced. 59 different proteins were identified using FEZ2 (207-353) as a bait (Table 1, Fig. 3). These can be organized into the following groups according the major described GO (Gene-Ontology) process attributed to them: transcription (10 proteins), translation (6 proteins), apoptosis (5 proteins), signal transduction (5 proteins), neuronal cell development (5 proteins), cytoskeleton/ centrosome (3 proteins), unknown function (10 proteins) and other functions (15).
Most interestingly, of these 59 FEZ2 interactors only 40 interacted also with FEZ1 (Fig. 3, Table 1) and indeed 8 of the 40 were identical with proteins that had been identified in the screen with FEZ1 [10]. This result may suggest that FEZ2 gained in respect to FEZ1 additional 19 new interactors, within the set of preys analyzed, which are specific to FEZ2. Such an evolutionary gain of function is compatible with the histological expression patterns of both proteins. Whereas FEZ1 (as its ortholog UNC-76) is expressed only in the nervous system, FEZ2 shows a ubiquitous expression pattern. Thinking of the role of FEZ proteins as transport adaptor proteins, FEZ2 may have evolved to transport additional cargo proteins, which may be expressed in nonneuronal tissues.
UNC-76 and FEZ1 share a conserved protein-protein interaction pattern Next we were interested to involve UNC-76 as a third component in our analysis, to try to understand how the protein-protein interaction pattern may have evolved during evolution of the FEZ family proteins. To this end we employed a C-terminal construct of UNC-76 (242-378) in a direct one-on-one analysis against all 59 FEZ2 interacting proteins (Fig. 3, Table 1).
The cDNA for C. elegans UNC-76 was a kind gift from Dr. Hengartner.
To our great surprise, the pattern of interaction of UNC-76, in relation to the proteins identified in the screening assay of twohybrid and FEZ2 as bait, was almost identical to that observed for FEZ1 (Table 1, Fig. 3). This suggests that the protein-protein interaction profile of the orthologs FEZ1 and UNC-76 is highly conserved. Only 4 of the 40 proteins that interacted with FEZ2 and FEZ1 did not interact with UNC-76. In conclusion the protein interaction pattern for both human FEZ1 and C. elegans UNC-76 seemed to be mostly conserved since the divergence of the deuterostomes (represented by the human FEZ1) from the ''protostomes'' (represented by C. elegans UNC-76), which occurred around 650 million years ago.
In this context it is worthwhile to differentiate between nuclear and cytoplasmic proteins that interact with the FEZ family proteins. A significant proportion of the interacting proteins are nuclear or predicted to be nuclear (18 of 59, 30.51%). Among these 18 nuclear proteins identified in the screen with FEZ2, a high proportion of 15 (ca. 83%) also interacted with FEZ1 (Fig. 4). This may be of potential relevance because we have shown in previous cell fractionation experiments that GFP-FEZ1 can be found in the nuclear fraction [26]. The exact role of FEZ1 in the nucleus remains however to be established. Interestingly, of the 41 cytoplasmic FEZ2 interacting proteins identified, a significant fraction of 16 (ca. 40%) interact exclusively with FEZ2. This again may suggest the acquisition of new interaction partners and thereby possibly the aggregation of new ''cytoplasmic functions'' for FEZ2 relative to FEZ1/UNC-76.
In 2000, Walhout and collaborators [28] introduced the concept of an ''interologue'', referring to a protein-protein interaction pattern that is conserved between pairs of orthologs. This notion is intuitive in molecular biology and has been documented by many interactions and may be conserved even in distantly related organisms, such as bacteria and animals. Our results show not only a high degree of preservation of the interactions between FEZ1, FEZ2 and UNC-76 but also the acquisition of new PPIs for FEZ2. Makino and coworkers [29] suggested that the rate of evolution of a protein that has more partners PPI is much smaller than the one that has few partners  PPI. This is in good agreement with the high degree of sequence similarity among FEZ1/FEZ2/UNC-76, which is especially pronounced in the C-terminal coiled-coil region, which is essential for most protein-protein interactions ( Figure S2). Only 12.8% of the amino acids in the coiled-coil regions of FEZ1 and FEZ2 are different. Three different fates have been proposed after the duplication of a genes [29]. First, one of the copies can be silenced by the accumulation of deleterious mutations and eventually become indistinguishable from genomic non-coding neighboring regions, while the other copy retains the original function. Second, while one copy maintains the original function, the other acquires a new function possibly by aggregation of advantageous mutations. Third, both copies accumulate mutations that alter the original function, but the original function is retained and compensated cooperatively. Our results seem to strongly suggest that for the FEZ family the second pathway is in progression. All interactors identified for FEZ1 in our laboratory have been confirmed to interact with FEZ2, too [10], but only 60% of the interactions identified for FEZ2 were confirmed for FEZ1 and UNC-76. Thus, FEZ2 seems to have acquired additional new functions but still preserved the majority of interactions found for FEZ1/UNC-76.
In this context it is also interesting to analyze the origin of the proteins fished with the FEZ2 bait from either bone marrow or fetal brain cDNA libraries and to correlate this with the supposed expression patterns (FEZ1: brain, FEZ2: ubiquitous) and function (FEZ1: neuron specific, FEZ2: ubiquitous) ( Table 1). Of the 17 proteins that were of origin from either bone marrow or brain (not from both, like 2 additonal proteins) and that were specific to interact with FEZ2, 13 were from bone marrow and only 4 from brain. Of the 52 proteins that were of origin from either bone marrow or brain only (7 additional proteins were fished from both libraries!) and that interacted with both FEZ1 and FEZ2, 32 were from the bone marrow library and 20 from the brain. Thus, FEZ1 interactions are relatively enriched in clones from the brain, while FEZ2 interactions are of more ubiquitous origin. In fact the ratio bone marrow proteins/brain protein for FEZ2 is almost three times higher (3.25; 13/4) than for FEZ1 (1.19; 19/16). This tendency may already indicate the acquisition of additional interacting partners and functions in the more ubiquitously expressed FEZ2 relative to the brain specific FEZ1 since the candidate cargo proteins for transport of FEZ1 can be mostly found in brain tissue, whereas those of FEZ2 preferentially in other tissues, including bone marrow.

Discussion
The large number of proteins that were identified to interact FEZ1, led to its classification as a hub protein [10]. Most proteins initially characterized as hubs tend to maintain this status after additional studies [30]. Indeed, here we were able to confirm the hub status for human FEZ1 and even extend it to the entire family of FEZ proteins, since both human FEZ2 as well as C. elegans UNC-76 interacted also with a large number of proteins. Genomewide studies have shown that deletions of hub protein encoding genes are three times more likely to be lethal than deletions of nonhubs, a phenomenon known as the centrality-lethality rule [31]. This is confirmed in C. elegans mutations in the unc-76 gene, which  It is depicted the mininum length of the retrieved sequences which could be visualized by forward DNA sequencing only, 3 Interaction confirmed (+) by yeast two hybrid system with UNC-76, FEZ1 or FEZ2 proteins; 4 Biological process based on the GO database (other functions may be known). doi:10.1371/journal.pone.0017426.t001 although not lethal, cause deleterious abnormalities in the elongation of neuronal axons along other axonal surfaces (but not over non-neuronal surfaces) during the animals development, which results in a ''paralyzed'' phenotype, where the worm has severly reduced body movements (therefore the designation unc for uncoordinated movements). These abnormalities can be complemented by the expression of human FEZ1 [1]. This result already shows the high functional conservation which we confirmed and detailed here: the protein interaction profiles of FEZ1 and UNC-76 are highly overlapping, since 36 of 40 FEZ1 interactors also interact with UNC-76 (90%). In this context it is also interesting to note that although the knockout of UNC-76 had severe and deleterious effects in the worm a recently described knockout mice for FEZ1 showed no dramatic changes in the brains anatomy or, the morphology of neuronal cell bodies, dendrites or axons, nor of the general locomotion or behavior of the animals [32]. Only after new social interactions or under stress conditions a slight hyper locomotion behavior was observed as well as a greater response to certain psycho-stimulants [32]. Our comparative analysis of the FEZ1 and FEZ2 interactomes becomes interesting at this point, because in mammals, and chordates in general, we have two FEZ family gene members in contrast to the C. elegans genome, which contains only unc-76. The duplication of the unc-76 gene occurred most likely in the stemline of the chordates (Table S1, Fig. 1) [6]. Our current interpretation of the lack of any greater neurophysiological or behavioural abnormalities in the FEZ1-/mice, may be due to a compensatory expression of the ubiquitously expressed FEZ2, which -as we demonstrated here -not only acquired additional new interactors in respect to FEZ1, but at the same time also maintained all of the interactors of FEZ1. From our point of view it is therefore not at all surprising that the FEZ1-/-mice had such a mild phenotype in comparison to the UNC-76 -/-worm [32]. It may even be that expression of FEZ2 in the FEZ1 k.o. mice, by hitherto unknown mechanisms, may be up-regulated in comparison to wild-type mice. It is interesting to note that in Figure S1A of the article of Sakae and co-workers, we can observe a new band at the height of the FEZ1 protein in extracts from the stomach of FEZ1-/-but not wild type mice (anti-FEZ1 Western blot). Since FEZ2 has almost the same MW as FEZ1, this band may represent up-regulated FEZ2 protein, in response to the lack of FEZ1. The observed relatively mild phenotype in the FEZ -/-mice may be as well attributed to additional proteins interacting with FEZ2, which may interfere in the way of a dominant negative in FEZ1 specific functions, as to specific but not yet identified PPI for FEZ1.
Although FEZ2 has more interactors than FEZ1 there is not a 100% functional redundancy. This would explain why both genes are still maintained in the gene pool (from a population genetics point of view). The reasons why both are maintained may be manifold, including dosage effects or a promoter specific regulation that must be maintained in the neurons (FEZ1). Data from the literature suggest that FEZ1 expression is increased upon NGF treatment, a neuronal tissue specific factor [6,33]. It may be even speculated that FEZ1 is somehow involved in neurogenesis per se. FEZ1 mRNA expression shows a peak that coincides with neurogenesis in rats, and decreases with the progress of adulthood. The mRNA expression FEZ2, although with higher levels in embryogenesis, is relatively constant, always less than that of FEZ1. Additional data on FEZ2 however are not available so far but one could speculate that its more ubiquitous expression depends on internal genetic programs rather than on external stimuli, like in the case of FEZ1.
Many of the proteins identified by our screen with FEZ2 interact with proteins that also interact with each other, increasing therefore the likelihood that these are true interactions of biological relevance. We like to highlight here the interactions of pleiotrophin (PTN), protein tyrosine phosphatase, receptor type S (PTPRS) ( [18], and retinoic acid receptor alpha (RARA). In 2005, Stelzl et al. [18] performed a yeast two-hybrid screen on a large scale, reporting more than 3000 interactions between human proteins. The human protein FEZ1 interacted with 21 proteins. Of these 21 proteins, 4 interact with pleiotrophin (general transcription factor IIF, protocollagen hydroxylase, neurokinin beta, translocase of outer mitochondrial membrane 20 homolog). FEZ1 interacted further with a receptor protein tyrosine phosphatase (PTP) (sigma) of the same family of receptors of pleiotrophin (beta and zeta) [34]. Pleiotrophin is also known as neurite growth-promoting factor 1 (NEGF1), has mitotic activity and influences the growth of neuritis. It acts in the developing nervous system [35] and signals via the inactivation of a receptordependent ligand of RPTP beta/zeta (or PTPRS). It increases the level of tyrosine phosphorylated signaling proteins in the cell, promoting thereby downstream signaling events [34]. In the nervous system, PTN is expressed in both neurons and glia cells [36] as well as in neuronal stem cells [37]. SAGE analysis indicated a 10-fold elevation of PTN transcripts in response to long-term NGF treatment of PC12 cells. It is of interest, that in P19 neuroprogenitor cells, PTN was also induced by retinoic acid during neuronal differentiation [38,39]. Retinoic acid is essential for the differentiation and morphogenesis of various tissues, including the nervous system [40]. Silencing of FEZ1 by RNA interference (RNAi) strongly inhibited the NGF-induced differentiation and efficiently reduced the anterograde transport of mitochondria in PC12 cells, suggesting that FEZ1 is essential for NGF-induced neuronal differentiation of PC12 cells [33,41].
Retinoic acid receptors (RARs) are transcriptional activators which are activated by binding of their ligands and act in concert with a combination retinoid X type receptors (RXR alpha, beta and gamma). Even in the absence of ligand, RAR-RXR heterodimers bind to DNA sequences known as response elements for RA, but in this state they recruit co-repressors. These corepressors mediate the negative modulation of transcription by recruitment of the histone deacetylase complex and the transfer of methyl radicals to the DNA bound histones, thereby stabilizing the nucleosome. The binding of the ligand retinoic acid causes however conformational changes in the binding domain of RAR, which in return promote the release of the co-repressors and the recruitment of co-activators. While some co-activators interact with factors of the basal transcription machinery, others induce chromatin remodeling and specific transcriptional activation. Data from our previous two-hybrid and pull-down assays confirmed the interaction of FEZ1 and FEZ2 proteins related to activation of transcription (BAF60a) and the repression of transcription (SAP30L) [10]. Furthermore, RXR/RAR interacts with N-Cor, thereby recruiting the complex of repression of transcription [42]. It is known that SAP30 also interacts with N-Cor [43] promoting the compaction of chromatin. These data gain potential relevance by the observation that both RAR and FEZ1/FEZ2 interact with protein complexes of chromatin remodeling, and with each other. Thus, FEZ1 may be also involved in the regulation of gene transcription mediated by transcription factors.
In summary, we have described here the molecular expansion of the FEZ protein family in the chordate stemline, most likely after two rounds of whole-genome duplication. We have identified many new members of the FEZ family from different species and a highly conserved C-terminal domain, which is responsible to mediate the majority of the protein-protein interactions. Our analyses indicate that the ancestral UNC-76 protein function and protein interaction profile is surprisingly conserved in animals, evolutionary distinct as vertebrates and nematodes. The analysis of the phylogeny and the detailed protein interaction profiles uncovered a likely functional divergence between FEZ1 and FEZ2, since FEZ2 acquired in respect to FEZ1/UNC-76 additional new interactors and is ubiquitously expressed and not restricted to neuronal cells, as is FEZ1/UNC-76. Our data provide an explanation for the ability of human FEZ1 to rescue the defects caused by unc-76 mutations in nematodes, since the PPI pattern of FEZ1 and UNC-76 are highly similar. Furthermore, no strong defects are caused in mice lacking the FEZ1, probably due to the compensatory presence of its paralogue FEZ2, which interacted with all FEZ1 interacting proteins, and additional new interactors that are FEZ2 specific. Our new data will hopefully stimulate and facilitate further studies on the functional role of proteins of the FEZ family in the development and function of both neuronal and non-neuronal cells. Furthermore, studies such as ours may contribute to our ever growing understanding of how protein networks function and specially how they came about and were modified and adapted during evolution. In the future we may devise new and more complete ways to visualize and comprehend how these networks were shaped throughout time. Figure S1 Amino acid sequence alignment of the members of the FEZ family. Names of the sequences are given in the Table S1. The residues in the alignment are shaded light grey, grey, or black to indicate shared identity in 40%, 70% or 100% of the analyzed sequences, respectively. The bars indicate regions predicted to form coiled-coil. (PDF) Figure S2 Identities and similarities between human proteins FEZ1 and FEZ2 and C. elegans UNC-76. A general scheme of the FEZ family proteins is shown at the top. The identity and similarity comparisons were made of two-by-two proteins both for the complete protein alignment as well as for local alignment of FEZ fragments by NPS@ (http://npsa-pbil.ibcp.fr/cgi-bin/npsa_automat.pl?page = npsa_clustalw.html). I = identity, S = strongly similar, W = weakly similar, D = different. (TIF)