Profiling of Histidine Phosphoproteome in Danio rerio by TiO2 Enrichment

Histidine phosphorylation is a reversible post‐translational modification that is known to regulate signal transduction in prokaryotes. However, functional studies in eukaryotes have been largely neglected due to the labile nature of N‐linked phosphorylated amino acids. In an effort to help elucidate the heretofore hidden vertebrate phosphoproteome, this report presents a global phosphorylation analysis of Danio rerio (zebrafish) larvae. Phosphopeptide enrichment is performed using a TiO₂ affinity technique. A total of 68 unique phosphohistidine sites are detected on 63 proteins among 1076 unique phosphosites on 708 proteins. Data are available via ProteomeXchange with identifier PXD012735. This report provides the first phosphohistidine dataset obtained from zebrafish.     

Protein kinases phosphorylate long disordered regions in intrinsically disordered proteins

Phosphorylation is a major post‐translational modification that plays a central role in signaling pathways. Protein kinases phosphorylate substrates (phosphoproteins) by adding phosphate at Ser/Thr or Tyr residues (phosphosites). A large amount of data identifying and describing phosphosites in phosphoproteins has been reported but the specificity of phosphorylation is not fully resolved. In this report, data of kinase‐substrate pairs identified by the Kinase‐Interacting Substrate Screening (KISS) method were used to analyze phosphosites in intrinsically disordered regions (IDRs) of intrinsically disordered proteins. We compared phosphorylated and nonphosphorylated IDRs and found that the phosphorylated IDRs were significantly longer than nonphosphorylated IDRs. The phosphorylated IDR is often the longest IDR (71%) in a phosphoprotein when only a single phosphosite exists in the IDR, and when the phosphoprotein has multiple phosphosites in an IDR(s), the phosphosites are primarily localized in a single IDR (78%) and this IDR is usually the longest one (81%). We constructed a stochastic model of phosphorylation to estimate the effect of IDR length. The model that accounted for IDR length produced more realistic results when compared with a model that excluded the IDR length. We propose that the IDR length is a significant determinant for locating kinase phosphorylation sites in phosphoproteins.     

In‐depth and 3‐dimensional exploration of the budding yeast phosphoproteome

Phosphorylation is one of the most dynamic and widespread post‐translational modifications regulating virtually every aspect of eukaryotic cell biology. Here, we assemble a dataset from 75 independent phosphoproteomic experiments performed in our laboratory using Saccharomyces cerevisiae. We report 30,902 phosphosites identified from cells cultured in a range of DNA damage conditions and/or arrested in distinct cell cycle stages. To generate a comprehensive resource for the budding yeast community, we aggregate our dataset with the Saccharomyces Genome Database and another recently published study, resulting in over 46,000 budding yeast phosphosites. With the goal of enhancing the identification of functional phosphorylation events, we perform computational positioning of phosphorylation sites on available 3D protein structures and systematically identify events predicted to regulate protein complex architecture. Results reveal hundreds of phosphorylation sites mapping to or near protein interaction interfaces, many of which result in steric or electrostatic “clashes” predicted to disrupt the interaction. With the advancement of Cryo‐EM and the increasing number of available structures, our approach should help drive the functional and spatial exploration of the phosphoproteome.     

On the incidence of intron loss and gain in paralogous gene families

Understanding gene duplication and gene structure evolution are fundamental goals of molecular evolutionary biology. A previous study by Babenko et al. (2004. Prevalence of intron gain over intron loss in the evolution of paralogous gene families. Nucleic Acids Res. 32:3724-3733) employed Dollo parsimony to infer spliceosomal intron losses and gains in paralogous gene families and concluded that there was a general excess of gains over losses. This result contrasts with patterns in orthologous genes, in which most lineages show an excess of intron losses over gains, suggesting the possibility of fundamentally different modes of intron evolution between orthologous and paralogous genes. We further studied the data and found a low level of intron position conservation with outgroups, and this led to problems with using Dollo parsimony to analyze the data. Statistical reanalysis of the data suggests, instead, that intron losses have outnumbered intron gains in paralogous gene families.     

Phosphorylation of ubiquitin at Ser65 affects its polymerization,targets, and proteome‐wide turnover

Ubiquitylation is an essential post‐translational modification that regulates numerous cellular processes, most notably protein degradation. Ubiquitin can itself be phosphorylated at nearly every serine, threonine, and tyrosine residue. However, the effect of this modification on ubiquitin function is largely unknown. Here, we characterized the effects of phosphorylation of yeast ubiquitin at serine 65 in vivo and in vitro. We find this post‐translational modification to be regulated under oxidative stress, occurring concomitantly with the restructuring of the ubiquitin landscape into a highly polymeric state. Phosphomimetic mutation of S65 recapitulates the oxidative stress phenotype, causing a dramatic accumulation of ubiquitylated proteins and a proteome‐wide reduction of protein turnover rates. Importantly, this mutation impacts ubiquitin chain disassembly, chain linkage distribution, ubiquitin interactions, and substrate targeting. These results demonstrate that phosphorylation is an additional mode of ubiquitin regulation with broad implications in cellular physiology.     

A model for the evolution of paralog families in genomes

We introduce and analyse a simple probabilistic model of genome evolution. It is based on three fundamental evolutionary events: gene loss, duplication and accumulated change. This is motivated by previous works which consisted in fitting the available genomic data into, what is called paralog distributions. This formalism is described by a system of infinite number of linear equations. We show that this system generates a semigroup of linear operators on the space l ¹. We prove that size distribution of paralogous gene families in a genome converges to the equilibrium as time goes to infinity. Moreover we show that when probabilities of gene removal and duplication are close to each other, then the resulting distribution is close to logarithmic distribution. Some empirical results for yeast genomes are presented.     

Yeast ribosomes: variety is the spice of life

In the budding yeast Saccharomyces cerevisiae, 59 of the 79 cytoplasmic ribosomal proteins are encoded by two genes, stemming from an ancient genome duplication event. Komili et al. (2007) now report that these paralogous genes are not functionally equivalent, suggesting the possible existence of a "ribosome code."     

Evolution and regulation of cellular periodic processes: a role for paralogues

Several cyclic processes take place within a single organism. For example, the cell cycle is coordinated with the 24 h diurnal rhythm in animals and plants, and with the 40 min ultradian rhythm in budding yeast. To examine the evolution of periodic gene expression during these processes, we performed the first systematic comparison in three organisms (Homo sapiens, Arabidopsis thaliana and Saccharomyces cerevisiae) by using public microarray data. We observed that although diurnal‐regulated and ultradian‐regulated genes are not generally cell‐cycle‐regulated, they tend to have cell‐cycle‐regulated paralogues. Thus, diverged temporal expression of paralogues seems to facilitate cellular orchestration under different periodic stimuli. Lineage‐specific functional repertoires of periodic‐associated paralogues imply that this mode of regulation might have evolved independently in several organisms.     

Functional analysis of the single Est1/Ebs1 homologue in Kluyveromyces lactis reveals roles in both telomere maintenance and rapamycin resistance

Est1 and Ebs1 in Saccharomyces cerevisiae are paralogous proteins that arose through whole-genome duplication and that serve distinct functions in telomere maintenance and translational regulation. Here we present our functional analysis of the sole Est1/Ebs1 homologue in the related budding yeast Kluyveromyces lactis (named KlEst1). We show that similar to other Est1s, KlEst1 is required for normal telomere maintenance in vivo and full telomerase primer extension activity in vitro. KlEst1 also associates with telomerase RNA (Ter1) and an active telomerase complex in cell extracts. Both the telomere maintenance and the Ter1 association functions of KlEst1 require its N-terminal domain but not its C terminus. Analysis of clusters of point mutations revealed residues in both the N-terminal TPR subdomain and the downstream helical subdomain (DSH) that are important for telomere maintenance and Ter1 association. A UV cross-linking assay was used to establish a direct physical interaction between KlEst1 and a putative stem-loop in Ter1, which also requires both the TPR and DSH subdomains. Moreover, similar to S. cerevisiae Ebs1 (ScEbs1) (but not ScEst1), KlEst1 confers rapamycin sensitivity and may be involved in nonsense-mediated decay. Interestingly, unlike telomere regulation, this apparently separate function of KlEst1 requires its C-terminal domain. Our findings provide insights on the mechanisms and evolution of Est1/Ebs1 homologues in budding yeast and present an attractive model system for analyzing members of this multifunctional protein family.     

Dynamic Phosphoproteome Profiling of Zebrafish Embryonic Fibroblasts during Cold Acclimation

Temperature affects almost all aspects of the fish life. To cope with low temperature, fish have evolved the ability of cold acclimation for survival. However, intracellular signaling events underlying cold acclimation in fish remain largely unknown. Here, the formation of cold acclimation in zebrafish embryonic fibroblasts (ZF4) is monitored and the phosphorylation events during the process are investigated through a large‐scale quantitative phosphoproteomic approach. In total, 11 474 phosphorylation sites are identified on 4066 proteins and quantified 5772 phosphosites on 2519 proteins. Serine, threonine, and tyrosine (Ser/Thr/Tyr) phosphorylation accounted for 85.5%, 13.3%, and 1.2% of total phosphosites, respectively. Among all phosphosites, 702 phosphosites on 510 proteins show differential regulation during cold acclimation of ZF4 cells. These phosphosites are divided into six clusters according to their dynamic changes during cold exposure. Kinase–substrate prediction reveals that mitogen‐activated protein kinase (MAPK) among the kinase groups is predominantly responsible for phosphorylation of these phosphosites. The differentially regulated phosphoproteins are functionally associated with various cellular processes such as regulation of actin cytoskeleton and MAPK signaling pathway. These data enrich the database of protein phosphorylation sites in zebrafish and provide key clues for the elucidation of intracellular signaling networks during cold acclimation of fish.     

Analysis of lamprey and hagfish genes reveals a complex history of gene duplications during early vertebrate evolution

It has been proposed that two events of duplication of the entire genome occurred early in vertebrate history (2R hypothesis). Several phylogenetic studies with a few gene families (mostly Hox genes and proteins from the MHC) have tried to confirm these polyploidization events. However, data from a single locus cannot explain the evolutionary history of a complete genome. To study this 2R hypothesis, we have taken advantage of the phylogenetic position of the lamprey to study the history of gene duplications in vertebrates. We selected most gene families that contain several paralogous genes in vertebrates and for which lamprey genes and an out-group are known in databases. In addition, we isolated members of the nuclear receptor superfamily in lamprey. Hagfish genes were also analyzed and found to confirm the lamprey gene analysis. Consistent with the 2R hypothesis, the phylogenetic analysis of 33 selected gene families, dispersed through the whole genome, revealed that one period of gene duplication arose before the lamprey-gnathostome split and this was followed by a second period of gene duplication after the lamprey-gnathostome split. Nevertheless, our analysis suggests that numerous gene losses and other gene-genome duplications occurred during the evolution of the vertebrate genomes. Thus, the complexity of all the paralogy groups present in vertebrates should be explained by the contribution of genome duplications (2R hypothesis), extra gene duplications, and gene losses.     

A general tendency for conservation of protein length across eukaryotic kingdoms

Protein elongation can occur in many ways, such as domain duplication or insertion and as recruitment of a transposable element fragment into the coding region, and it is believed to be a general tendency in protein evolution. Indeed, a previous study showed that yeast proteins are, on average, longer than their orthologs in bacteria, and in this study, we found that proteins in yeast, nematode, Drosophila, human, and Arabidopsis are, on average, longer than their orthologs in Escherichia coli. Surprisingly, however, we found conservation of protein sequence length across eukaryotic kingdoms. We collected 1,252 orthologous proteins from yeast, nematode, Drosophila, human, and Arabidopsis and found that the total length of these proteins is very similar among the five species and that there is no general tendency for a protein to increase or decrease in length. Furthermore, although paralogous proteins tend to undergo more sequence-length changes, there is also no general tendency for length increase. However, proteins that are commonly shared by Drosophila and human but not by yeast are, on average, substantially longer than proteins that are shared by yeast, Drosophila, and human. This is a puzzle that begs for an answer.     

Has gene duplication impacted the evolution of Eutherian longevity?

One of the greatest unresolved questions in aging biology is determining the genetic basis of interspecies longevity variation. Gene duplication is often the key to understanding the origin and evolution of important Eutherian phenotypes. We systematically identified longevity‐associated genes in model organisms that duplicated throughout Eutherian evolution. Longevity‐associated gene families have a marginally significantly higher rate of duplication compared to non‐longevity‐associated gene families. Anti‐longevity‐associated gene families have significantly increased rate of duplication compared to pro‐longevity gene families and are enriched in neurodegenerative disease categories. Conversely, duplicated pro‐longevity‐associated gene families are enriched in cell cycle genes. There is a cluster of longevity‐associated gene families that expanded solely in long‐lived species that is significantly enriched in pathways relating to 3‐UTR‐mediated translational regulation, metabolism of proteins and gene expression, pathways that have the potential to affect longevity. The identification of a gene cluster that duplicated solely in long‐lived species involved in such fundamental processes provides a promising avenue for further exploration of Eutherian longevity evolution.     

Phosphoproteome dynamics during mitotic exit in budding yeast

The cell division cycle culminates in mitosis when two daughter cells are born. As cyclin‐dependent kinase (Cdk) activity reaches its peak, the anaphase‐promoting complex/cyclosome (APC/C) is activated to trigger sister chromatid separation and mitotic spindle elongation, followed by spindle disassembly and cytokinesis. Degradation of mitotic cyclins and activation of Cdk‐counteracting phosphatases are thought to cause protein dephosphorylation to control these sequential events. Here, we use budding yeast to analyze phosphorylation dynamics of 3,456 phosphosites on 1,101 proteins with high temporal resolution as cells progress synchronously through mitosis. This reveals that successive inactivation of S and M phase Cdks and of the mitotic kinase Polo contributes to order these dephosphorylation events. Unexpectedly, we detect as many new phosphorylation events as there are dephosphorylation events. These correlate with late mitotic kinase activation and identify numerous candidate targets of these kinases. These findings revise our view of mitotic exit and portray it as a dynamic process in which a range of mitotic kinases contribute to order both protein dephosphorylation and phosphorylation.     

Phosphorylation of the FUS low‐complexity domain disrupts phase separation,aggregation, and toxicity

Neuronal inclusions of aggregated RNA‐binding protein fused in sarcoma (FUS) are hallmarks of ALS and frontotemporal dementia subtypes. Intriguingly, FUS's nearly uncharged, aggregation‐prone, yeast prion‐like, low sequence‐complexity domain (LC) is known to be targeted for phosphorylation. Here we map in vitro and in‐cell phosphorylation sites across FUS LC. We show that both phosphorylation and phosphomimetic variants reduce its aggregation‐prone/prion‐like character, disrupting FUS phase separation in the presence of RNA or salt and reducing FUS propensity to aggregate. Nuclear magnetic resonance spectroscopy demonstrates the intrinsically disordered structure of FUS LC is preserved after phosphorylation; however, transient domain collapse and self‐interaction are reduced by phosphomimetics. Moreover, we show that phosphomimetic FUS reduces aggregation in human and yeast cell models, and can ameliorate FUS‐associated cytotoxicity. Hence, post‐translational modification may be a mechanism by which cells control physiological assembly and prevent pathological protein aggregation, suggesting a potential treatment pathway amenable to pharmacologic modulation.     

Phosphorylation of Argonaute proteins affects mRNA binding and is essential for microRNA‐guided gene silencing in vivo

Argonaute proteins associate with microRNAs and are key components of gene silencing pathways. With such a pivotal role, these proteins represent ideal targets for regulatory post‐translational modifications. Using quantitative mass spectrometry, we find that a C‐terminal serine/threonine cluster is phosphorylated at five different residues in human and Caenorhabditis elegans. In human, hyper‐phosphorylation does not affect microRNA binding, localization, or cleavage activity of Ago2. However, mRNA binding is strongly affected. Strikingly, on Ago2 mutants that cannot bind microRNAs or mRNAs, the cluster remains unphosphorylated indicating a role at late stages of gene silencing. In C. elegans, the phosphorylation of the conserved cluster of ALG‐1 is essential for microRNA function in vivo. Furthermore, a single point mutation within the cluster is sufficient to phenocopy the loss of its complete phosphorylation. Interestingly, this mutant retains its capacity to produce and bind microRNAs and represses expression when artificially tethered to an mRNA. Altogether, our data suggest that the phosphorylation state of the serine/threonine cluster is important for Argonaute–mRNA interactions.     

Co‐ and post‐translational modifications of the 26S proteasome in yeast

The yeast (Saccharomyces cerevisiae) 26S proteasome consists of the 19S regulatory particle (19S RP) and 20S proteasome subunits. We detected comprehensively co‐ and post‐translational modifications of these subunits using proteomic techniques. First, using MS/MS, we investigated the N‐terminal modifications of three 19S RP subunits, Rpt1, Rpn13, and Rpn15, which had been unclear, and found that the N‐terminus of Rpt1 is not modified, whereas that of Rpn13 and Rpn15 is acetylated. Second, we identified a total of 33 Ser/Thr phosphorylation sites in 15 subunits of the proteasome. The data obtained by us and other groups reveal that the 26S proteasome contains at least 88 phospho‐amino acids including 63 pSer, 23 pThr, and 2 pTyr residues. Dephosphorylation treatment of the 19S RP with λ phosphatase resulted in a 30% decrease in ATPase activity, demonstrating that phosphorylation is involved in the regulation of ATPase activity in the proteasome. Third, we tried to detect glycosylated subunits of the 26S proteasome. However, we identified neither N‐ and O‐linked oligosaccharides nor O‐linked β‐N‐acetylglucosamine in the 19S RP and 20S proteasome subunits. To date, a total of 110 co‐ and post‐translational modifications, including N^α‐acetylation, N^α‐myristoylation, and phosphorylation, in the yeast 26S proteasome have been identified.     

Linkage mapping with paralogs exposes regions of residual tetrasomic inheritance in chum salmon (Oncorhynchus keta)

Gene sequence similarity due to shared ancestry after a duplication event, that is paralogy, complicates the assessment of genetic variation, as sequences originating from paralogs can be difficult to distinguish. These confounded sequences are often removed prior to further analyses, leaving the underlying loci uncharacterized. Salmonids have only partially rediploidized subsequent to a whole‐genome duplication; residual tetrasomic inheritance has been observed in males. We present a maximum‐likelihood‐based method to resolve confounded paralogous loci by observing the segregation of alleles in gynogenetic haploid offspring and demonstrate its effectiveness by constructing two linkage maps for chum salmon (Oncorhynchus keta), with and without these newly resolved loci. We find that the resolved paralogous loci are not randomly distributed across the genome. A majority are clustered in expanded subtelomeric regions of 14 linkage groups, suggesting a significant fraction of the chum salmon genome may be missed by the exclusion of paralogous loci. Transposable elements have been proposed as drivers of genome evolution and, in salmonids, may have an important role in the rediploidization process by driving differentiation between homeologous chromosomes. Consistent with that hypothesis, we find a reduced fraction of transposable element annotations among paralogous loci, and these loci predominately occur in the genomic regions that lag in the rediploidization process.