The substrate quality of CK2 target sites has a determinant role on their function and evolution.

Most biological processes are regulated by signaling modules that bind to short linear motifs. For protein kinases, substrates may have full or only partial matches to the kinase recognition motif, a property known as "substrate quality." However, it is not clear whether differences in substrate quality represent neutral variation or if they have functional consequences. We examine this question for the kinase CK2, which has many fundamental functions. We show that optimal CK2 sites are phosphorylated at maximal stoichiometries and found in many conditions, whereas minimal substrates are more weakly phosphorylated and have regulatory functions. Optimal CK2 sites tend to be more conserved, and substrate quality is often tuned by selection. For intermediate sites, increases or decreases in substrate quality may be deleterious, as we demonstrate for a CK2 substrate at the kinetochore. The results together suggest a strong role for substrate quality in phosphosite function and evolution. A record of this paper's transparent peer review process is included in the supplemental information.

Deep mutational scanning of Pneumocystis jirovecii dihydrofolate reductase reveals allosteric mechanism of resistance to an antifolate.

Pneumocystis jirovecii is a fungal pathogen that causes pneumocystis pneumonia, a disease that mainly affects immunocompromised individuals. This fungus has historically been hard to study because of our inability to grow it in vitro. One of the main drug targets in P. jirovecii is its dihydrofolate reductase (PjDHFR). Here, by using functional complementation of the baker's yeast ortholog, we show that PjDHFR can be inhibited by the antifolate methotrexate in a dose-dependent manner. Using deep mutational scanning of PjDHFR, we identify mutations conferring resistance to methotrexate. Thirty-one sites spanning the protein have at least one mutation that leads to resistance, for a total of 355 high-confidence resistance mutations. Most resistance-inducing mutations are found inside the active site, and many are structurally equivalent to mutations known to lead to resistance to different antifolates in other organisms. Some sites show specific resistance mutations, where only a single substitution confers resistance, whereas others are more permissive, as several substitutions at these sites confer resistance. Surprisingly, one of the permissive sites (F199) is without direct contact to either ligand or cofactor, suggesting that it acts through an allosteric mechanism. Modeling changes in binding energy between F199 mutants and drug shows that most mutations destabilize interactions between the protein and the drug. This evidence points towards a more important role of this position in resistance than previously estimated and highlights potential unknown allosteric mechanisms of resistance to antifolate in DHFRs. Our results offer unprecedented resources for the interpretation of mutation effects in the main drug target of an uncultivable fungal pathogen.

Parallel Nonfunctionalization of CK1δ/ε Kinase Ohnologs Following a Whole-Genome Duplication Event.

Whole-genome duplication (WGD) followed by speciation allows us to examine the parallel evolution of ohnolog pairs. In the yeast family Saccharomycetaceae, HRR25 is a rare case of repeated ohnolog maintenance. This gene has reverted to a single copy in Saccharomyces cerevisiae where it is now essential, but has been maintained as pairs in at least 7 species post-WGD. In S. cerevisiae, HRR25 encodes the casein kinase 1δ/ε and plays a role in a variety of functions through its kinase activity and protein-protein interactions (PPIs). We hypothesized that the maintenance of duplicated HRR25 ohnologs could be a result of repeated subfunctionalization. We tested this hypothesis through a functional complementation assay in S. cerevisiae, testing all pairwise combinations of 25 orthologs (including 7 ohnolog pairs). Contrary to our expectations, we observed no cases of pair-dependent complementation, which would have supported the subfunctionalization hypothesis. Instead, most post-WGD species have one ohnolog that failed to complement, suggesting their nonfunctionalization or neofunctionalization. The ohnologs incapable of complementation have undergone more rapid protein evolution, lost most PPIs that were observed for their functional counterparts and singletons from post-WGD and non-WGD species, and have nonconserved cellular localization, consistent with their ongoing loss of function. The analysis in Naumovozyma castellii shows that the noncomplementing ohnolog is expressed at a lower level and has become nonessential. Taken together, our results indicate that HRR25 orthologs are undergoing gradual nonfunctionalization.

Parallel nonfunctionalization of CK1δ/ε kinase ohnologs following a whole-genome duplication event.

Whole genome duplication (WGD) followed by speciation allows us to examine the parallel evolution of ohnolog pairs. In the yeast family Saccharomycetaceae, HRR25 is a rare case of repeated ohnolog maintenance. This gene has reverted to a single copy in S. cerevisiae where it is now essential, but has been maintained as pairs in at least 7 species post WGD. In S. cerevisiae, HRR25 encodes the casein kinase (CK) 1δ/ε and plays a role in a variety of functions through its kinase activity and protein-protein interactions (PPIs). We hypothesized that the maintenance of duplicated HRR25 ohnologs could be a result of repeated subfunctionalization. We tested this hypothesis through a functional complementation assay in S. cerevisiae, testing all pairwise combinations of 25 orthologs (including 7 ohnolog pairs). Contrary to our expectations, we observed no cases of pair-dependent complementation, which would have supported the subfunctionalization hypothesis. Instead, most post-WGD species have one ohnolog that failed to complement, suggesting their nonfunctionalization or neofunctionalization. The ohnologs incapable of complementation have undergone more rapid protein evolution, lost most PPIs that were observed for their functional counterparts and singletons from post and non-WGD species, and have non-conserved cellular localization, consistent with their ongoing loss of function. The analysis in N. castelli shows that the non-complementing ohnolog is expressed at a lower level and has become non-essential. Taken together, our results indicate that HRR25 orthologs are undergoing gradual nonfunctionalization.

Evolutionary Trajectories are Contingent on Mitonuclear Interactions.

Critical mitochondrial functions, including cellular respiration, rely on frequently interacting components expressed from both the mitochondrial and nuclear genomes. The fitness of eukaryotic organisms depends on a tight collaboration between both genomes. In the face of an elevated rate of evolution in mtDNA, current models predict that the maintenance of mitonuclear compatibility relies on compensatory evolution of the nuclear genome. Mitonuclear interactions would therefore exert an influence on evolutionary trajectories. One prediction from this model is that the same nuclear genome evolving with different mitochondrial haplotypes would follow distinct molecular paths toward higher fitness. To test this prediction, we submitted 1,344 populations derived from 7 mitonuclear genotypes of Saccharomyces cerevisiae to >300 generations of experimental evolution in conditions that either select for a mitochondrial function or do not strictly require respiration for survival. Performing high-throughput phenotyping and whole-genome sequencing on independently evolved individuals, we identified numerous examples of gene-level evolutionary convergence among populations with the same mitonuclear background. Phenotypic and genotypic data on strains derived from this evolution experiment identify the nuclear genome and the environment as the main determinants of evolutionary divergence, but also show a modulating role for the mitochondrial genome exerted both directly and via interactions with the two other components. We finally recapitulated a subset of prominent loss-of-function alleles in the ancestral backgrounds and confirmed a generalized pattern of mitonuclear-specific and highly epistatic fitness effects. Together, these results demonstrate how mitonuclear interactions can dictate evolutionary divergence of populations with identical starting nuclear genotypes.

Epistasis between promoter activity and coding mutations shapes gene evolvability.

The evolution of protein-coding genes proceeds as mutations act on two main dimensions: regulation of transcription level and the coding sequence. The extent and impact of the connection between these two dimensions are largely unknown because they have generally been studied independently. By measuring the fitness effects of all possible mutations on a protein complex at various levels of promoter activity, we show that promoter activity at the optimal level for the wild-type protein masks the effects of both deleterious and beneficial coding mutations. Mutations that are deleterious at low activity but masked at optimal activity are slightly destabilizing for individual subunits and binding interfaces. Coding mutations that increase protein abundance are beneficial at low expression but could potentially incur a cost at high promoter activity. We thereby demonstrate that promoter activity in interaction with protein properties can dictate which coding mutations are beneficial, neutral, or deleterious.

Expression attenuation as a mechanism of robustness against gene duplication.

Gene duplication is ubiquitous and a major driver of phenotypic diversity across the tree of life, but its immediate consequences are not fully understood. Deleterious effects would decrease the probability of retention of duplicates and prevent their contribution to long-term evolution. One possible detrimental effect of duplication is the perturbation of the stoichiometry of protein complexes. Here, we measured the fitness effects of the duplication of 899 essential genes in the budding yeast using high-resolution competition assays. At least 10% of genes caused a fitness disadvantage when duplicated. Intriguingly, the duplication of most protein complex subunits had small to nondetectable effects on fitness, with few exceptions. We selected four complexes with subunits that had an impact on fitness when duplicated and measured the impact of individual gene duplications on their protein-protein interactions. We found that very few duplications affect both fitness and interactions. Furthermore, large complexes such as the 26S proteasome are protected from gene duplication by attenuation of protein abundance. Regulatory mechanisms that maintain the stoichiometric balance of protein complexes may protect from the immediate effects of gene duplication. Our results show that a better understanding of protein regulation and assembly in complexes is required for the refinement of current models of gene duplication.

Frequent Assembly of Chimeric Complexes in the Protein Interaction Network of an Interspecies Yeast Hybrid.

Hybrids between species often show extreme phenotypes, including some that take place at the molecular level. In this study, we investigated the phenotypes of an interspecies diploid hybrid in terms of protein-protein interactions inferred from protein correlation profiling. We used two yeast species, Saccharomyces cerevisiae and Saccharomyces uvarum, which are interfertile, but yet have proteins diverged enough to be differentiated using mass spectrometry. Most of the protein-protein interactions are similar between hybrid and parents, and are consistent with the assembly of chimeric complexes, which we validated using an orthogonal approach for the prefoldin complex. We also identified instances of altered protein-protein interactions in the hybrid, for instance, in complexes related to proteostasis and in mitochondrial protein complexes. Overall, this study uncovers the likely frequent occurrence of chimeric protein complexes with few exceptions, which may result from incompatibilities or imbalances between the parental proteomes.

The role of structural pleiotropy and regulatory evolution in the retention of heteromers of paralogs.

Gene duplication is a driver of the evolution of new functions. The duplication of genes encoding homomeric proteins leads to the formation of homomers and heteromers of paralogs, creating new complexes after a single duplication event. The loss of these heteromers may be required for the two paralogs to evolve independent functions. Using yeast as a model, we find that heteromerization is frequent among duplicated homomers and correlates with functional similarity between paralogs. Using in silico evolution, we show that for homomers and heteromers sharing binding interfaces, mutations in one paralog can have structural pleiotropic effects on both interactions, resulting in highly correlated responses of the complexes to selection. Therefore, heteromerization could be preserved indirectly due to selection for the maintenance of homomers, thus slowing down functional divergence between paralogs. We suggest that paralogs can overcome the obstacle of structural pleiotropy by regulatory evolution at the transcriptional and post-translational levels.

Turnover of ribosome-associated transcripts from de novo ORFs produces gene-like characteristics available for de novo gene emergence in wild yeast populations.

Little is known about the rate of emergence of de novo genes, what their initial properties are, and how they spread in populations. We examined wild yeast populations (Saccharomyces paradoxus) to characterize the diversity and turnover of intergenic ORFs over short evolutionary timescales. We find that hundreds of intergenic ORFs show translation signatures similar to canonical genes, and we experimentally confirmed the translation of many of these ORFs in laboratory conditions using a reporter assay. Compared with canonical genes, intergenic ORFs have lower translation efficiency, which could imply a lack of optimization for translation or a mechanism to reduce their production cost. Translated intergenic ORFs also tend to have sequence properties that are generally close to those of random intergenic sequences. However, some of the very recent translated intergenic ORFs, which appeared <110 kya, already show gene-like characteristics, suggesting that the raw material for functional innovations could appear over short evolutionary timescales.

A collection of barcoded natural isolates of Saccharomyces paradoxus to study microbial evolutionary ecology.

While the use of barcoded collections of laboratory microorganisms and the development of barcode-based cell tracking are rapidly developing in genetics and genomics research, tools to track natural populations are still lacking. The yeast Saccharomyces paradoxus is an emergent microbial model in ecology and evolution. More than five allopatric and sympatric lineages have been identified and hundreds of strains have been isolated for this species, allowing to assess the impact of natural diversity on complex traits. We constructed a collection of 550 barcoded and traceable strains of S. paradoxus, including all three North American lineages SpB, SpC, and SpC*. These strains are diploid, many have their genome fully sequenced and are barcoded with a unique 20 bp sequence that allows their identification and quantification. This yeast collection is functional for competitive experiments in pools as the barcodes allow to measure each lineage's and individual strains' fitness in common conditions. We used this tool to demonstrate that in the tested conditions, there are extensive genotype-by-environment interactions for fitness among S. paradoxus strains, which reveals complex evolutionary potential in variable environments. This barcoded collection provides a valuable resource for ecological genomics studies that will allow gaining a better understanding of S. paradoxus evolution and fitness-related traits.

Extended Linkers Improve the Detection of Protein-protein Interactions (PPIs) by Dihydrofolate Reductase Protein-fragment Complementation Assay (DHFR PCA) in Living Cells.

Understanding the function of cellular systems requires describing how proteins assemble with each other into transient and stable complexes and to determine their spatial relationships. Among the tools available to perform these analyses on a large scale is Protein-fragment Complementation Assay based on the dihydrofolate reductase (DHFR PCA). Here we test how longer linkers between the fusion proteins and the reporter fragments affect the performance of this assay. We investigate the architecture of the RNA polymerases, the proteasome and the conserved oligomeric Golgi (COG) complexes in living cells and performed large-scale screens with these extended linkers. We show that longer linkers significantly improve the detection of protein-protein interactions and allow to measure interactions further in space than the standard ones. We identify new interactions, for instance between the retromer complex and proteins related to autophagy and endocytosis. Longer linkers thus contribute an enhanced additional tool to the existing toolsets for the detection and measurements of protein-protein interactions and protein proximity in living cells.

Deep transcriptome annotation enables the discovery and functional characterization of cryptic small proteins.

Recent functional, proteomic and ribosome profiling studies in eukaryotes have concurrently demonstrated the translation of alternative open-reading frames (altORFs) in addition to annotated protein coding sequences (CDSs). We show that a large number of small proteins could in fact be coded by these altORFs. The putative alternative proteins translated from altORFs have orthologs in many species and contain functional domains. Evolutionary analyses indicate that altORFs often show more extreme conservation patterns than their CDSs. Thousands of alternative proteins are detected in proteomic datasets by reanalysis using a database containing predicted alternative proteins. This is illustrated with specific examples, including altMiD51, a 70 amino acid mitochondrial fission-promoting protein encoded in MiD51/Mief1/SMCR7L, a gene encoding an annotated protein promoting mitochondrial fission. Our results suggest that many genes are multicoding genes and code for a large protein and one or several small proteins.

Gene duplication can impart fragility, not robustness, in the yeast protein interaction network.

The maintenance of duplicated genes is thought to protect cells from genetic perturbations, but the molecular basis of this robustness is largely unknown. By measuring the interaction of yeast proteins with their partners in wild-type cells and in cells lacking a paralog, we found that 22 out of 56 paralog pairs compensate for the lost interactions. An equivalent number of pairs exhibit the opposite behavior and require each other's presence for maintaining their interactions. These dependent paralogs generally interact physically, regulate each other's abundance, and derive from ancestral self-interacting proteins. This reveals that gene duplication may actually increase mutational fragility instead of robustness in a large number of cases.

Evolutionary rescue by compensatory mutations is constrained by genomic and environmental backgrounds.

Since deleterious mutations may be rescued by secondary mutations during evolution, compensatory evolution could identify genetic solutions leading to therapeutic targets. Here, we tested this hypothesis and examined whether these solutions would be universal or would need to be adapted to one's genetic and environmental make-ups. We performed experimental evolutionary rescue in a yeast disease model for the Wiskott-Aldrich syndrome in two genetic backgrounds and carbon sources. We found that multiple aspects of the evolutionary rescue outcome depend on the genotype, the environment, or a combination thereof. Specifically, the compensatory mutation rate and type, the molecular rescue mechanism, the genetic target, and the associated fitness cost varied across contexts. The course of compensatory evolution is therefore highly contingent on the initial conditions in which the deleterious mutation occurs. In addition, these results reveal biologically favored therapeutic targets for the Wiskott-Aldrich syndrome, including the target of an unrelated clinically approved drug. Our results experimentally illustrate the importance of epistasis and environmental evolutionary constraints that shape the adaptive landscape and evolutionary rate of molecular networks.

Systematic identification of signal integration by protein kinase A.

Cellular processes and homeostasis control in eukaryotic cells is achieved by the action of regulatory proteins such as protein kinase A (PKA). Although the outbound signals from PKA directed to processes such as metabolism, growth, and aging have been well charted, what regulates this conserved regulator remains to be systematically identified to understand how it coordinates biological processes. Using a yeast PKA reporter assay, we identified genes that influence PKA activity by measuring protein-protein interactions between the regulatory and the two catalytic subunits of the PKA complex in 3,726 yeast genetic-deletion backgrounds grown on two carbon sources. Overall, nearly 500 genes were found to be connected directly or indirectly to PKA regulation, including 80 core regulators, denoting a wide diversity of signals regulating PKA, within and beyond the described upstream linear pathways. PKA regulators span multiple processes, including the antagonistic autophagy and methionine biosynthesis pathways. Our results converge toward mechanisms of PKA posttranslational regulation by lysine acetylation, which is conserved between yeast and humans and that, we show, regulates protein complex formation in mammals and carbohydrate storage and aging in yeast. Taken together, these results show that the extent of PKA input matches with its output, because this kinase receives information from upstream and downstream processes, and highlight how biological processes are interconnected and coordinated by PKA.

N-terminal entrance loop of yeast Yps1 and O-glycosylation of substrates are determinant factors controlling the shedding activity of this GPI-anchored endopeptidase.

BACKGROUND: S. cerevisiae Yps1 is the prototypical aspartic endopeptidase of the fungal yapsin family. This glycosylphosphatidylinositol (GPI) anchored enzyme was recently shown to be involved in the shedding of the GPI proteins Utr2, Gas1 and itself. It was also proposed to be part of a novel quality control mechanism that eliminates excess and/or misfolded GPI proteins. What regulates its shedding activity at the cell surface is however poorly understood. Yps1 is initially synthesized as a zymogen requiring proteolytic activation to remove a pro-peptide and further processing within a large insertion loop (N-entrance loop) generates a two-subunit endopeptidase. To investigate the role of this loop on its shedding activity, which typically takes place within Ser/Thr-rich domains, it was replaced with the short peptide found at the analogous position in Yps3. We also tested whether O-glycosylation might protect against proteolytic processing by Yps1. RESULTS: We show here that replacement of the N-entrance loop (N-ent loop) of Yps1 generates a single chain endopeptidase that undergoes partial (pH 6.0) or complete (pH 3.0) pro-peptide removal. At both pH, the shedding activity of the chimeric endopeptidase (Yps1-DL) toward Gas1 and itself is strongly and drastically increased, respectively. A direct correlation between endoproteolytic cleavage of this loop in native Yps1 and its shedding is observed. The Yps1-dependent shedding of two model GPI proteins (Gas1 and Yps1) is also stimulated by the absence of the O-mannosyltransferases, Pmt4 and Pmt2 respectively, involved in O-glycosylation of their Ser/Thr-rich domains. Under these conditions, some Yps1-independent shedding is also observed. CONCLUSIONS: Partial pro-peptide removal is essential to produce a functional Yps1 endopeptidase. The Yps1 N-ent loop plays a major role in regulating the shedding activity of the endopeptidase, most likely by limiting access to the active site, and its cleavage in native Yps1 is associated with its shedding. O-glycosylation protects against Yps1-dependent and -independent shedding of GPI proteins. It is postulated that hypoglycosylation of cell surface proteins, which may occur for misfolded proteins that escaped the ER-associated degradation, might target their elimination through shedding by Yps1 and possibly other yapsin members.

Detecting functional divergence after gene duplication through evolutionary changes in posttranslational regulatory sequences.

Gene duplication is an important evolutionary mechanism that can result in functional divergence in paralogs due to neo-functionalization or sub-functionalization. Consistent with functional divergence after gene duplication, recent studies have shown accelerated evolution in retained paralogs. However, little is known in general about the impact of this accelerated evolution on the molecular functions of retained paralogs. For example, do new functions typically involve changes in enzymatic activities, or changes in protein regulation? Here we study the evolution of posttranslational regulation by examining the evolution of important regulatory sequences (short linear motifs) in retained duplicates created by the whole-genome duplication in budding yeast. To do so, we identified short linear motifs whose evolutionary constraint has relaxed after gene duplication with a likelihood-ratio test that can account for heterogeneity in the evolutionary process by using a non-central chi-squared null distribution. We find that short linear motifs are more likely to show changes in evolutionary constraints in retained duplicates compared to single-copy genes. We examine changes in constraints on known regulatory sequences and show that for the Rck1/Rck2, Fkh1/Fkh2, Ace2/Swi5 paralogs, they are associated with previously characterized differences in posttranslational regulation. Finally, we experimentally confirm our prediction that for the Ace2/Swi5 paralogs, Cbk1 regulated localization was lost along the lineage leading to SWI5 after gene duplication. Our analysis suggests that changes in posttranslational regulation mediated by short regulatory motifs systematically contribute to functional divergence after gene duplication.

Modulation of the yeast protein interactome in response to DNA damage.

Cells deploy diverse mechanisms to physiologically adapt to potentially detrimental perturbations. These mechanisms include changes in the organization of protein-protein interaction networks (PINs). Most PINs characterized to date are portrayed in a single environmental condition and are thus likely to miss important connections among biological processes. In this report, we show that the yeast DHFR-PCA on high-density arrays allows to detects modulations of protein-protein interactions (PPIs) in different conditions by testing more than 1000 PPIs in standard and in a drug-inducing DNA damage conditions. We identify 156 PPIs that show significant modulation in response to DNA damage. We provide evidence that modulated PPIs involve essential genes (NOP7, EXO84 and LAS17) playing critical roles in response to DNA damage. Additionally, we show that a significant proportion of PPI changes are likely explained by changes in protein localization and, to a lesser extent, protein abundance. The protein interaction modules affected by changing PPIs support the role of mRNA stability and translation, protein degradation and ubiquitylation and the regulation of the actin cytoskeleton in response to DNA damage. Overall, we provide a valuable tool and dataset for the study of the rewiring of PINs in response to environmental perturbations. BIOLOGICAL SIGNIFICANCE: We show that the DHFR-PCA is a high-throughput method that allows the detection of changes in PPIs associated with different environmental conditions using DNA damage response as a testbed. We provide a valuable resource for the study of DNA damage in eukaryotic cells. This article is part of a Special Issue: Can Proteomics Fill the Gap Between Genomics and Phenotypes?