Genome-Wide Identification,Classification, and Expression Analysis of 14-3-3 Gene Family in Populus

Background

In plants, 14-3-3 proteins are encoded by a large multigene family and are involved in signaling pathways to regulate plant development and protection from stress. Although twelve Populus 14-3-3s were identified based on the Populus trichocarpa genome V1.1 in a previous study, no systematic analysis including genome organization, gene structure, duplication relationship, evolutionary analysis and expression compendium has been conducted in Populus based on the latest P. trichocarpa genome V3.0.

Principal Findings

Here, a comprehensive analysis of Populus 14-3-3 family is presented. Two new 14-3-3 genes were identified based on the latest P. trichocarpa genome. In P. trichocarpa, fourteen 14-3-3 genes were grouped into ε and non-ε group. Exon-intron organizations of Populus 14-3-3s are highly conserved within the same group. Genomic organization analysis indicated that purifying selection plays a pivotal role in the retention and maintenance of Populus 14-3-3 family. Protein conformational analysis indicated that Populus 14-3-3 consists of a bundle of nine α-helices (α1-α9); the first four are essential for formation of the dimer, while α3, α5, α7, and α9 form a conserved peptide-binding groove. In addition, α1, α3, α5, α7, and α9 were evolving at a lower rate, while α2, α4, and α6 were evolving at a relatively faster rate. Microarray analyses showed that most Populus 14-3-3s are differentially expressed across tissues and upon exposure to various stresses.

Conclusions

The gene structures and their coding protein structures of Populus 14-3-3s are highly conserved among group members, suggesting that members of the same group might also have conserved functions. Microarray and qRT-PCR analyses showed that most Populus 14-3-3s were differentially expressed in various tissues and were induced by various stresses. Our investigation provided a better understanding of the complexity of the 14-3-3 gene family in poplars.     

Dissection of Binding between a Phosphorylated Tyrosine Hydroxylase Peptide and 14-3-3ζ: A Complex Story Elucidated by NMR

Human tyrosine hydroxylase activity is regulated by phosphorylation of its N-terminus and by an interaction with the modulator 14-3-3 proteins. We investigated the binding of singly or doubly phosphorylated and thiophosphorylated peptides, comprising the first 50 amino acids of human tyrosine hydroxylase, isoform 1 (hTH1), that contain the critical interaction domain, to 14-3-3ζ, by ³¹P NMR. Single phosphorylation at S19 generates a high affinity 14-3-3ζ binding epitope, whereas singly S40-phosphorylated peptide interacts with 14-3-3ζ one order-of-magnitude weaker than the S19-phosphorylated peptide. Analysis of the binding data revealed that the 14-3-3ζ dimer and the S19- and S40-doubly phosphorylated peptide interact in multiple ways, with three major complexes formed: 1), a single peptide bound to a 14-3-3ζ dimer via the S19 phosphate with the S40 phosphate occupying the other binding site; 2), a single peptide bound to a 14-3-3ζ dimer via the S19 phosphorous with the S40 free in solution; or 3), a 14-3-3ζ dimer with two peptides bound via the S19 phosphorous to each binding site. Our system and data provide information as to the possible mechanisms by which 14-3-3 can engage binding partners that possess two phosphorylation sites on flexible tails. Whether these will be realized in any particular interacting pair will naturally depend on the details of each system.     

Regulation of the Water Channel Aquaporin-2 via 14-3-3θ and -ζ

The 14-3-3 family of proteins are multifunctional proteins that interact with many of their cellular targets in a phosphorylation-dependent manner. Here, we determined that 14-3-3 proteins interact with phosphorylated forms of the water channel aquaporin-2 (AQP2) and modulate its function. With the exception of σ, all 14-3-3 isoforms were abundantly expressed in mouse kidney and mouse kidney collecting duct cells (mpkCCD₁₄). Long-term treatment of mpkCCD₁₄ cells with the type 2 vasopressin receptor agonist dDAVP increased mRNA and protein levels of AQP2 alongside 14-3-3β and -ζ, whereas levels of 14-3-3η and -θ were decreased. Co-immunoprecipitation (co-IP) studies in mpkCCD₁₄ cells uncovered an AQP2/14-3-3 interaction that was modulated by acute dDAVP treatment. Additional co-IP studies in HEK293 cells determined that AQP2 interacts selectively with 14-3-3ζ and -θ. Use of phosphatase inhibitors in mpkCCD₁₄ cells, co-IP with phosphorylation deficient forms of AQP2 expressed in HEK293 cells, or surface plasmon resonance studies determined that the AQP2/14-3-3 interaction was modulated by phosphorylation of AQP2 at various sites in its carboxyl terminus, with Ser-256 phosphorylation critical for the interactions. shRNA-mediated knockdown of 14-3-3ζ in mpkCCD₁₄ cells resulted in increased AQP2 ubiquitylation, decreased AQP2 protein half-life, and reduced AQP2 levels. In contrast, knockdown of 14-3-3θ resulted in increased AQP2 half-life and increased AQP2 levels. In conclusion, this study demonstrates phosphorylation-dependent interactions of AQP2 with 14-3-3θ and -ζ. These interactions play divergent roles in modulating AQP2 trafficking, phosphorylation, ubiquitylation, and degradation.     

Identification of Novel 14-3-3 Residues That Are Critical for Isoform-specific Interaction with GluN2C to Regulate N-Methyl-d-aspartate (NMDA) Receptor Trafficking

The 14-3-3 family of proteins is widely distributed in the CNS where they are major regulators of essential neuronal functions. There are seven known mammalian 14-3-3 isoforms (ζ,, τ, ϵ, η, β, and σ), which generally function as adaptor proteins. Previously, we have demonstrated that 14-3-3ϵ isoform dynamically regulates forward trafficking of GluN2C-containing NMDA receptors (NMDARs) in cerebellar granule neurons, that when expressed on the surface, promotes neuronal survival following NMDA-induced excitotoxicity. Here, we report 14-3-3 isoform-specific binding and functional regulation of GluN2C. In particular, we show that GluN2C C-terminal domain (CTD) binds to all 14-3-3 isoforms except 14-3-3σ, and binding is dependent on GluN2C serine 1096 phosphorylation. Co-expression of 14-3-3 (ζ and ϵ) and GluN1/GluN2C promotes the forward delivery of receptors to the cell surface. We further identify novel residues serine 145, tyrosine 178, and cysteine 189 on α-helices 6, 7, and 8, respectively, within ζ-isoform as part of the GluN2C binding motif and independent of the canonical peptide binding groove. Mutation of these conserved residues abolishes GluN2C binding and has no functional effect on GluN2C trafficking. Reciprocal mutation of alanine 145, histidine 180, and isoleucine 191 on 14-3-3σ isoform promotes GluN2C binding and surface expression. Moreover, inhibiting endogenous 14-3-3 using a high-affinity peptide inhibitor, difopein, greatly diminishes GluN2C surface expression. Together, these findings highlight the isoform-specific structural and functional differences within the 14-3-3 family of proteins, which determine GluN2C binding and its essential role in targeting the receptor to the cell surface to facilitate glutamatergic neurotransmission.     

Identification of 14-3-3 Proteins Phosphopeptide-Binding Specificity Using an Affinity-Based Computational Approach

The 14-3-3 proteins are a highly conserved family of homodimeric and heterodimeric molecules, expressed in all eukaryotic cells. In human cells, this family consists of seven distinct but highly homologous 14-3-3 isoforms. 14-3-3σ is the only isoform directly linked to cancer in epithelial cells, which is regulated by major tumor suppressor genes. For each 14-3-3 isoform, we have 1,000 peptide motifs with experimental binding affinity values. In this paper, we present a novel method for identifying peptide motifs binding to 14-3-3σ isoform. First, we propose a sampling criteria to build a predictor for each new peptide sequence. Then, we select nine physicochemical properties of amino acids to describe each peptide motif. We also use auto-cross covariance to extract correlative properties of amino acids in any two positions. Finally, we consider elastic net to predict affinity values of peptide motifs, based on ridge regression and least absolute shrinkage and selection operator (LASSO). Our method tests on the 1,000 known peptide motifs binding to seven 14-3-3 isoforms. On the 14-3-3σ isoform, our method has overall pearson-product-moment correlation coefficient (PCC) and root mean squared error (RMSE) values of 0.84 and 252.31 for N–terminal sublibrary, and 0.77 and 269.13 for C–terminal sublibrary. We predict affinity values of 16,000 peptide sequences and relative binding ability across six permutated positions similar with experimental values. We identify phosphopeptides that preferentially bind to 14-3-3σ over other isoforms. Several positions on peptide motifs are in the same amino acid category with experimental substrate specificity of phosphopeptides binding to 14-3-3σ. Our method is fast and reliable and is a general computational method that can be used in peptide-protein binding identification in proteomics research.     

Activation of Phosphatidylinositol 3-Kinase (PI3K) and Mitogen-activated Protein Kinase (MAPK) Signaling and the Consequent Induction of Transformation by Overexpressed 14-3-3γ Protein Require Specific Amino Acids within 14-3-3γ N-terminal Variable Region II

Members of the 14-3-3 superfamily regulate numerous cellular functions by binding phosphoproteins. The seven human isoforms (and the myriad of other eukaryotic 14-3-3 proteins) are highly conserved in amino acid sequence and secondary structure, yet there is abundant evidence that the various isoforms manifest disparate as well as common functions. Several of the human 14-3-3 isoforms are dysregulated in certain cancers and thus have been implicated in oncogenesis; experimentally, 14-3-3γ behaves as an oncogene, whereas 14-3-3σ acts as a tumor suppressor. In this study, we sought to localize these opposing phenotypes to specific regions of the two isoforms and then to individual amino acids therein. Using a bioinformatics approach, six variable regions (VRI–VRVI) were identified. Using this information, two sets of constructs were created in which N-terminal portions (including either VRI–IV or only VRI and VRII) of 14-3-3γ and 14-3-3σ were swapped; NIH3T3 cells overexpressing the four chimeric proteins were tested for transformation activity (focus formation, growth in soft agar) and activation of PI3K and MAPK signaling. We found that the specific phenotypes of 14-3-3γ are associated with the N-terminal 40 amino acids (VRI and VRII); in like fashion, VRI and VRII of 14-3-3σ dictated its tumor suppressor function. Using individual amino acid substitutions within the 14-3-3γ VRII, we identified two residues required for and two contributing to the γ-specific phenotypes. Our observations suggest that isoform-specific phenotypes are dictated by a relatively few amino acids within variable regions.     

Modulation of 14-3-3/Phosphotarget Interaction by Physiological Concentrations of Phosphate and Glycerophosphates

Molecular mechanisms governing selective binding of a huge number of various phosphorylated protein partners to 14-3-3 remain obscure. Phosphate can bind to 14-3-3 and therefore being present at high intracellular concentration, which undergoes significant changes under physiological conditions, phosphate can theoretically regulate interaction of 14-3-3 with phosphorylated targets. In order to check this hypothesis we analyzed effect of phosphate and other natural abundant anions on interaction of 14-3-3 with phosphorylated human small heat shock protein HspB6 (Hsp20) participating in regulation of different intracellular processes. Inorganic phosphate, glycerol-1-phosphate and glycerol-2-phosphate at physiologically relevant concentrations (5-15 mM) significantly destabilized complexes formed by 14-3-3ζ and phosphorylated HspB6 (pHspB6), presumably, via direct interaction with the substrate-binding site of 14-3-3. Phosphate also destabilized complexes between pHspB6 and 14-3-3γ or the monomeric mutant form of 14-3-3ζ. Inorganic sulfate and pyrophosphate were less effective in modulation of 14-3-3 interaction with its target protein. The inhibitory effect of all anions on pHspB6/14-3-3 interaction was concentration-dependent. It is hypothesized that physiological changes in phosphate anions concentration can modulate affinity and specificity of interaction of 14-3-3 with its multiple targets and therefore the actual phosphointeractome of 14-3-3.     

Virus-induced gene silencing of <Emphasis Type="Italic">14-3-3</Emphasis> genes abrogates dark repression of nitrate reductase activity in <Emphasis Type="Italic">Nicotiana benthamiana</Emphasis>

In order to study the effect of repression of 14-3-3 genes on actual activity of the nitrate reductase (NR) in Nicotiana benthamiana leaves, Nb14-3-3a gene was silenced by virus-induced gene silencing (VIGS) method using potato virus X (PVX). Expression of Nb14-3-3a as well as Nb14-3-3b genes was altogether repressed in the leaves of PVX-14-3a-infected plants. Furthermore, two-dimensional gel electrophoresis and immunoblot analysis with anti-14-3-3 antiserum suggested that the expressions of Nb14-3-3a and Nb14-3-3b proteins are accordingly repressed in PVX-14-3a-infected plants. It is well known that binding of 14-3-3 proteins to phosphorylated NR leads to substantial decrease in NR activity of leaves under darkness. Therefore, we studied the changes in NR activity in response to light/dark transitions in the leaves of PVX-14-3a-infected plants. NR activation state was kept at a high level under darkness in PVX-14-3a-infected plants, but not in PVX-green fluorescent protein (GFP)-infected and control plants. This result suggests that Nb14-3-3a and/or Nb14-3-3b proteins are indeed involved in the inactivation of NR activity under darkness in N. benthamiana.     

Aberrant Upregulation of 14-3-3σ and EZH2 Expression Serves as an Inferior Prognostic Biomarker for Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is the fifth most common malignancy in the world. It is of important significance to find biomarkers for the prognostic monitoring of HCC. The 14-3-3σ and EZH2 proteins are involved in cell cycle regulation and epigenetic silencing. We herein examined the significance of 14-3-3 σ and EZH2 in HCC (n = 167) by immunohistochemistry, RT-PCR and qRT-PCR. The correlation between 14-3-3σ and EZH2 expression and patients'' clinicopathologic features were examined, as was the correlation between 14-3-3σ and EZH2 expression and the prognosis of HCC patients. We found that 14-3-3σ and EZH2 were highly expressed in HCC (71% and 90%), the expression of EZH2, but not 14-3-3σ, is associated with vascular invasion and tumor differentiation (p<0.01). The coexistence of 14-3-3σ and EZH2 overexpression is associated with a relatively unfavorable prognosis (p<0.01), suggesting that aberrant upregulation of 14-3-3σ and EZH2 expression serves as an inferior prognostic biomarker for HCC.     

Biophysical Characterization of Essential Phosphorylation at the Flexible C-Terminal Region of C-Raf with 14-3-3ζ Protein

Phosphorylation at the C-terminal flexible region of the C-Raf protein plays an important role in regulating its biological activity. Auto-phosphorylation at serine 621 (S621) in this region maintains C-Raf stability and activity. This phosphorylation mediates the interaction between C-Raf and scaffold protein 14-3-3ζ to activate the downstream MEK kinase pathway. In this study, we have defined the interaction of C-terminal peptide sequence of C-Raf with 14-3-3ζ protein and determined the possible structural adaptation of this region. Biophysical elucidation of the interaction was carried out using phosphopeptide (residue number 615–630) in the presence of 14-3-3ζ protein. Using isothermal titration calorimetry (ITC), a high binding affinity with micro-molar range was found to exist between the peptide and 14-3-3ζ protein, whereas the non-phosphorylated peptide did not show any appreciable binding affinity. Further interaction details were investigated using several biophysical techniques such as circular dichroism (CD), fluorescence, and nuclear magnetic resonance (NMR) spectroscopy, in addition to molecular modeling. This study provides the molecular basis for C-Raf C-terminal-derived phosphopeptide interaction with 14-3-3ζ protein as well as structural insights responsible for phosphorylated S621-mediated 14-3-3ζ binding at an atomic resolution.     

Antibodies to assess phosphorylation of spinach leaf nitrate reductase on serine 543 and its binding to 14-3-3 proteins.

To monitor site-specific phosphorylation of spinach leaf nitrate reductase (NR) and binding of the enzyme to 14-3-3 proteins, serum antibodies were raised that select for either serine 543 phospho- or dephospho-NR. The dephospho-specific antibodies blocked NR phosphorylation on serine 543. The phospho-specific antibodies prevented NR binding to 14-3-3s, NR inhibition by 14-3-3s, NR dephosphorylation on serine 543, and did not precipitate 14-3-3s together with NR. Together, this confirms that 14-3-3s bind to NR at hinge 1 after it has been phosphorylated on serine 543. The amounts of individual NR forms were determined in leaf extracts by immunoblotting and immunoprecipitation. The phosphorylation state of NR on serine 543 increased 2-3-fold in leaves upon a light/ dark transition. Before the transition, one-third of NR was already phosphorylated on serine 543 but was not bound to 14-3-3s. Phosphorylation of serine 543 seems not to be enough to bind to 14-3-3s in leaves.     

Inhibition of the Arabidopsis Salt Overly Sensitive Pathway by 14-3-3 Proteins

The Salt Overly Sensitive (SOS) pathway regulates intracellular sodium ion (Na⁺) homeostasis and salt tolerance in plants. Until recently, little was known about the mechanisms that inhibit the SOS pathway when plants are grown in the absence of salt stress. In this study, we report that the Arabidopsis thaliana 14-3-3 proteins λ and κ interact with SOS2 and repress its kinase activity. Growth in the presence of salt decreases the interaction between SOS2 and the 14-3-3 proteins, leading to kinase activation in planta. 14-3-3 λ interacts with the SOS2 junction domain, which is important for its kinase activity. A phosphorylation site (Ser-294) is identified within this domain by mass spectrometry. Mutation of Ser-294 to Ala or Asp does not affect SOS2 kinase activity in the absence of the 14-3-3 proteins. However, in the presence of 14-3-3 proteins, the inhibition of SOS2 activity is decreased by the Ser-to-Ala mutation and enhanced by the Ser-to-Asp exchange. These results identify 14-3-3 λ and κ as important regulators of salt tolerance. The inhibition of SOS2 mediated by the binding of 14-3-3 proteins represents a novel mechanism that confers basal repression of the SOS pathway in the absence of salt stress.     

A single Arabidopsis GF14 isoform possesses biochemical characteristics of diverse 14-3-3 homologues

Arabidopsis cDNA clones of GF14 proteins originally were isolated on the basis of their association with the G-box DNA/protein complex by a monoclonal antibody screening approach. GF14 proteins are homologous to the 14-3-3 family of mammalian proteins. Here we demonstrate that recombinant GF14 , one member of the Arabidopsis GF14 protein family, is a dimeric protein that possesses many of the attributes of diverse mammalian 14-3-3 homologues. GF14 activates rat brain tryptophan hydroxylase and protein kinase C in a manner similar to the bovine 14-3-3 protein. It also activates exoenzyme S of Pseudomonas aeruginosa as does bovine brain factor activating exoenzyme S (FAS), which is itself a member of 14-3-3 proteins. In addition, GF14 binds calcium, as does the human 14-3-3 homologue reported to be a phospholipase A2. These results indicate that a single isoform of this plant protein family can have multiple functions and that individual GF14 isoforms may have multiple roles in mediating signal transductions in plants. However, GF14 does not regulate growth in an in vivo test for functional similarity to the yeast 14-3-3 homologue, BMH1. Thus, while a single plant GF14 isoform can exhibit many of the biochemical attributes of diverse mammalian 14-3-3 homologues, open questions remain regarding the physiological functions of GF14/14-3-3 proteins.