Identification of novel phosphorylation modification sites in human proteins that originated after the human-chimpanzee divergence

Reorganization of 3D genome structure may contribute to gene regulatory evolution in primates

A growing body of evidence supports the notion that variation in gene regulation plays a crucial role in both speciation and adaptation. However, a comprehensive functional understanding of the mechanisms underlying regulatory evolution remains elusive. In primates, one of the crucial missing pieces of information towards a better understanding of regulatory evolution is a comparative annotation of interactions between distal regulatory elements and promoters. Chromatin conformation capture technologies have enabled genome-wide quantifications of such distal 3D interactions. However, relatively little comparative research in primates has been done using such technologies. To address this gap, we used Hi-C to characterize 3D chromatin interactions in induced pluripotent stem cells (iPSCs) from humans and chimpanzees. We also used RNA-seq to collect gene expression data from the same lines. We generally observed that lower-order, pairwise 3D genomic interactions are conserved in humans and chimpanzees, but higher order genomic structures, such as topologically associating domains (TADs), are not as conserved. Inter-species differences in 3D genomic interactions are often associated with gene expression differences between the species. To provide additional functional context to our observations, we considered previously published chromatin data from human stem cells. We found that inter-species differences in 3D genomic interactions, which are also associated with gene expression differences between the species, are enriched for both active and repressive marks. Overall, our data demonstrate that, as expected, an understanding of 3D genome reorganization is key to explaining regulatory evolution.

The way in which a genome folds affects the regulation of gene expression. This is often due to loops in the three-dimensional structure that bring linearly distant genes and regulatory elements into close proximity. Most studies examining three-dimensional structure genome-wide are limited to a single species. In this study, we compared three-dimensional structure in the genomes of induced pluripotent stem cells from humans and chimpanzees. We collected gene expression data from the same samples, which allowed us to assess the contribution of three-dimensional chromatin conformation to gene regulatory evolution in primates. Our results demonstrate that gene expression differences between the species may often be mediated by differences in three-dimensional genomic interactions. Our data also suggest that large-scale chromatin structures (i.e. topologically associating domains, TADs) are not well conserved in their placement across species. We hope the analytical paradigms we present here could serve as a basis for future comparative studies of three-dimensional genome organization, elucidating the putative functional regulatory loci driving speciation.

Multiple Innovations in Genetic and Epigenetic Mechanisms Cooperate to Underpin Human Brain Evolution

Our knowledge of how the human brain differs from those of other species in terms of evolutionary adaptations and functionality is limited. Comparative genomics reveal valuable insight, especially the expansion of human-specific noncoding regulatory and repeat-containing regions. Recent studies add to our knowledge of evolving brain function by investigating cellular mechanisms such as protein emergence, extensive sequence editing, retrotransposon activity, dynamic epigenetic modifications, and multiple noncoding RNA functions. These findings present an opportunity to combine newly discovered genetic and epigenetic mechanisms with more established concepts into a more comprehensive picture to better understand the uniquely evolved human brain.

Structural-functional interactions of NS1-BP protein with the splicing and mRNA export machineries for viral and host gene expression

The influenza virulence factor NS1 protein interacts with the cellular NS1-BP protein to promote splicing and nuclear export of the viral M mRNAs. The viral M1 mRNA encodes the M1 matrix protein and is alternatively spliced into the M2 mRNA, which is translated into the M2 ion channel. These proteins have key functions in viral trafficking and budding. To uncover the NS1-BP structural and functional activities in splicing and nuclear export, we performed proteomics analysis of nuclear NS1-BP binding partners and showed its interaction with constituents of the splicing and mRNA export machineries. NS1-BP BTB domains form dimers in the crystal. Full-length NS1-BP is a dimer in solution and forms at least a dimer in cells. Mutations suggest that dimerization is important for splicing. The central BACK domain of NS1-BP interacts directly with splicing factors such as hnRNP K and PTBP1 and with the viral NS1 protein. The BACK domain is also the site for interactions with mRNA export factor Aly/REF and is required for viral M mRNA nuclear export. The crystal structure of the C-terminal Kelch domain shows that it forms a β-propeller fold, which is required for the splicing function of NS1-BP. This domain interacts with the polymerase II C-terminal domain and SART1, which are involved in recruitment of splicing factors and spliceosome assembly, respectively. NS1-BP functions are not only critical for processing a subset of viral mRNAs but also impact levels and nuclear export of a subset of cellular mRNAs encoding factors involved in metastasis and immunity.

Loss of conserved ubiquitylation sites in conserved proteins during human evolution

Ubiquitylation of lysine residues in proteins serves a pivotal role in the efficient removal of misfolded or unused proteins and in the control of various regulatory pathways by monitoring protein activity that may lead to protein degradation. The loss of ubiquitylated lysines may affect the ubiquitin‑mediated regulatory network and result in the emergence of novel phenotypes. The present study analyzed mouse ubiquitylation data and orthologous proteins from 62 mammals to identify 193 conserved ubiquitylation sites from 169 proteins that were lost in the Euarchonta lineage leading to humans. A total of 8 proteins, including betaine homocysteine S‑methyltransferase, clin and CBS domain divalent metal cation transport mediator 3, ribosome‑binding protein 1 and solute carrier family 37 member 4, lost 1 conserved lysine residue, which was ubiquitylated in the mouse ortholog, following the human‑chimpanzee divergence. A total of 17 of the lost ubiquitylated lysines are also known to be modified by acetylation and/or succinylation in mice. In 8 cases, a novel lysine evolved at positions flanking the lost conserved lysine residues, potentially as a method of compensation. We hypothesize that the loss of ubiquitylation sites during evolution may lead to the development of advantageous phenotypes, which are then fixed by selection. The ancestral ubiquitylation sites identified in the present study may be a useful resource for investigating the association between loss of ubiquitylation sites and the emergence of novel phenotypes during evolution towards modern humans.

Natural variability of minimotifs in 1092 people indicates that minimotifs are targets of evolution

Since the function of a short contiguous peptide minimotif can be introduced or eliminated by a single point mutation, these functional elements may be a source of human variation and a target of selection. We analyzed the variability of ∼300 000 minimotifs in 1092 human genomes from the 1000 Genomes Project. Most minimotifs have been purified by selection, with a 94% invariance, which supports important functional roles for minimotifs. Minimotifs are generally under negative selection, possessing high genomic evolutionary rate profiling (GERP) and sitewise likelihood-ratio (SLR) scores. Some are subject to neutral drift or positive selection, similar to coding regions. Most SNPs in minimotif were common variants, but with minor allele frequencies generally <10%. This was supported by low substation rates and few newly derived minimotifs. Several minimotif alleles showed different intercontinental and regional geographic distributions, strongly suggesting a role for minimotifs in adaptive evolution. We also note that 4% of PTM minimotif sites in histone tails were common variants, which has the potential to differentially affect DNA packaging among individuals. In conclusion, minimotifs are a source of functional genetic variation in the human population; thus, they are likely to be an important target of selection and evolution.

RIP1 Cleavage in the Kinase Domain Regulates TRAIL-Induced NF-κB Activation and Lymphoma Survival

Although TRAIL is considered a potential anticancer agent, it enhances tumor progression by activating NF-κB in apoptosis-resistant cells. Cellular FLICE-like inhibitory protein (cFLIP) overexpression and caspase-8 activation have been implicated in TRAIL-induced NF-κB activation; however, the underlying mechanisms are unknown. Here, we report that caspase-8-dependent cleavage of RIP1 in the kinase domain (KD) and intermediate domain (ID) determines the activation state of the NF-κB pathway in response to tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) treatment. In apoptosis-sensitive cells, caspase-8 cleaves RIP1 in the KD and ID immediately after the recruitment of RIP1 to the receptor complex, impairing IκB kinase (IKK) recruitment and NF-κB activation. In apoptosis-resistant cells, cFLIP restricts caspase-8 activity, resulting in limited RIP1 cleavage and generation of a KD-cleaved fragment capable of activating NF-κB but not apoptosis. Notably, depletion of the cytoplasmic pool of TRAF2 and cIAP1 in lymphomas by CD40 ligation inhibits basal RIP1 ubiquitination but does not prompt cell death, due to CD40L-induced cFLIP expression and limited RIP1 cleavage. Inhibition of RIP1 cleavage at the KD suppresses NF-κB activation and cell survival even in cFLIP-overexpressing lymphomas. Importantly, RIP1 is constitutively cleaved in human and mouse lymphomas, suggesting that cFLIP-mediated and caspase-8-dependent limited cleavage of RIP1 is a new layer of mechanism that promotes NF-κB activation and lymphoma survival.

Loss of gene function and evolution of human phenotypes

Humans have acquired many distinct evolutionary traits after the human-chimpanzee divergence. These phenotypes have resulted from genetic changes that occurred in the human genome and were retained by natural selection. Comparative primate genome analyses reveal that loss-of-function mutations are common in the human genome. Some of these gene inactivation events were revealed to be associated with the emergence of advantageous phenotypes and were therefore positively selected and fixed in modern humans (the "less-ismore" hypothesis). Representative cases of human gene inactivation and their functional implications are presented in this review. Functional studies of additional inactive genes will provide insight into the molecular mechanisms underlying acquisition of various human-specific traits.

Evolution of transcription factor function as a mechanism for changing metazoan developmental gene regulatory networks

The form that an animal takes during development is directed by gene regulatory networks (GRNs). Developmental GRNs interpret maternally deposited molecules and externally supplied signals to direct cell-fate decisions, which ultimately leads to the arrangements of organs and tissues in the organism. Genetically encoded modifications to these networks have generated the wide range of metazoan diversity that exists today. Most studies of GRN evolution focus on changes to cis-regulatory DNA, and it was historically theorized that changes to the transcription factors that bind to these cis-regulatory modules (CRMs) contribute to this process only rarely. A growing body of evidence suggests that changes to the coding regions of transcription factors play a much larger role in the evolution of developmental gene regulatory networks than originally imagined. Just as cis-regulatory changes make use of modular binding site composition and tissue-specific modules to avoid pleiotropy, transcription factor coding regions also predominantly evolve in ways that limit the context of functional effects. Here, we review the recent works that have led to this unexpected change in the field of Evolution and Development (Evo-Devo) and consider the implications these studies have had on our understanding of the evolution of developmental processes.

The acquisition of novel N-glycosylation sites in conserved proteins during human evolution

Background

N-linked protein glycosylation plays an important role in various biological processes, including protein folding and trafficking, and cell adhesion and signaling. The acquisition of a novel N-glycosylation site may have significant effect on protein structure and function, and therefore, on the phenotype.

Results

We analyzed the human glycoproteome data set (2,534 N-glycosylation sites in 1,027 proteins) and identified 112 novel N-glycosylation sites in 91 proteins that arose in the human lineage since the last common ancestor of Euarchonta (primates and treeshrews). Three of them, Asn-196 in adipocyte plasma membrane-associated protein (APMAP), Asn-91 in cluster of differentiation 166 (CD166/ALCAM), and Asn-76 in thyroglobulin, are human-specific. Molecular evolutionary analysis suggested that these sites were under positive selection during human evolution. Notably, the Asn-76 of thyroglobulin might be involved in the increased production of thyroid hormones in humans, especially thyroxine (T4), because the removal of the glycan moiety from this site was reported to result in a significant decrease in T4 production.

Conclusions

We propose that the novel N-glycosylation sites described in this study may be useful candidates for functional analyses to identify innovative genetic modifications for beneficial phenotypes acquired in the human lineage.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0468-5) contains supplementary material, which is available to authorized users.

Turnover of protein phosphorylation evolving under stabilizing selection

Most proteins are regulated by posttranslational modifications and changes in these modifications contribute to evolutionary changes as well as to human diseases. Phosphorylation of serines, threonines, and tyrosines are the most common modifications identified to date in eukaryotic proteomes. While the mode of action and the function of most phosphorylation sites remain unknown, functional studies have shown that phosphorylation affects protein stability, localization and ability to interact. Two broad modes of action have been described for protein phosphorylation. The first mode corresponds to the canonical and qualitative view whereby single phosphorylation sites act as molecular switches that either turn on or off specific protein functions through direct or allosteric effects. The second mode is more akin to a rheostat than a switch. In this case, a group of phosphorylation sites in a given protein region contributes collectively to the modification of the protein, irrespective of the precise position of individual sites, through an aggregate property. Here we discuss these two types of regulation and examine how they affect the rate and patterns of protein phosphorylation evolution. We describe how the evolution of clusters of phosphorylation sites can be studied under the framework of complex traits evolution and stabilizing selection.

Functional divergence and evolutionary turnover in mammalian phosphoproteomes

Protein phosphorylation is a key mechanism to regulate protein functions. However, the contribution of this protein modification to species divergence is still largely unknown. Here, we studied the evolution of mammalian phosphoregulation by comparing the human and mouse phosphoproteomes. We found that 84% of the positions that are phosphorylated in one species or the other are conserved at the residue level. Twenty percent of these conserved sites are phosphorylated in both species. This proportion is 2.5 times more than expected by chance alone, suggesting that purifying selection is preserving phosphoregulation. However, we show that the majority of the sites that are conserved at the residue level are differentially phosphorylated between species. These sites likely result from false-negative identifications due to incomplete experimental coverage, false-positive identifications and non-functional sites. In addition, our results suggest that at least 5% of them are likely to be true differentially phosphorylated sites and may thus contribute to the divergence in phosphorylation networks between mouse and humans and this, despite residue conservation between orthologous proteins. We also showed that evolutionary turnover of phosphosites at adjacent positions (in a distance range of up to 40 amino acids) in human or mouse leads to an over estimation of the divergence in phosphoregulation between these two species. These sites tend to be phosphorylated by the same kinases, supporting the hypothesis that they are functionally redundant. Our results support the hypothesis that the evolutionary turnover of phosphorylation sites contributes to the divergence in phosphorylation profiles while preserving phosphoregulation. Overall, our study provides advanced analyses of mammalian phosphoproteomes and a framework for the study of their contribution to phenotypic evolution.

Understanding how differences in cellular regulation lead to phenotypic differences between species remains an open challenge in evolutionary genetics. The extensive phosphorylation data currently available allows to compare the human and mouse phosphoproteomes and to measure changes in their phosphoregulation. We found a general conservation of phosphorylation sites between these two species. However, a fraction of sites are conserved at the sequence level (the same amino acid is present in both species) but differ in their phosphorylation status. These sites represent candidate sites that have the potential to explain differences between human and mouse signalling networks that do not depend on the divergence of orthologous residues. Furthermore, we identified several sites where to a phosphorylation site in one species corresponds a non-phosphorylatable residue in the other one. These cases represent clear differences in protein regulation. Recent studies suggest that phosphorylation sites can shift position during evolution, leading to configurations in which pairs of divergent phosphorylation sites are functionally redundant. We identified more than 100 putative such cases, suggesting that divergence in amino acid does not necessarily imply functional divergence when comparing phosphoproteomes. Overall, our study provides new key concepts and data for the study of how regulatory differences may be linked to phenotypic ones at the network level.

Gains of ubiquitylation sites in highly conserved proteins in the human lineage

Background

Post-translational modification of lysine residues of specific proteins by ubiquitin modulates the degradation, localization, and activity of these target proteins. Here, we identified gains of ubiquitylation sites in highly conserved regions of human proteins that occurred during human evolution.

Results

We analyzed human ubiquitylation site data and multiple alignments of orthologous mammalian proteins including those from humans, primates, other placental mammals, opossum, and platypus. In our analysis, we identified 281 ubiquitylation sites in 252 proteins that first appeared along the human lineage during primate evolution: one protein had four novel sites; four proteins had three sites each; 18 proteins had two sites each; and the remaining 229 proteins had one site each. PML, which is involved in neurodevelopment and neurodegeneration, acquired three sites, two of which have been reported to be involved in the degradation of PML. Thirteen human proteins, including ERCC2 (also known as XPD) and NBR1, gained human-specific ubiquitylated lysines after the human-chimpanzee divergence. ERCC2 has a Lys/Gln polymorphism, the derived (major) allele of which confers enhanced DNA repair capacity and reduced cancer risk compared with the ancestral (minor) allele. NBR1 and eight other proteins that are involved in the human autophagy protein interaction network gained a novel ubiquitylation site.

Conclusions

The gain of novel ubiquitylation sites could be involved in the evolution of protein degradation and other regulatory networks. Although gains of ubiquitylation sites do not necessarily equate to adaptive evolution, they are useful candidates for molecular functional analyses to identify novel advantageous genetic modifications and innovative phenotypes acquired during human evolution.

Human-specific protein isoforms produced by novel splice sites in the human genome after the human-chimpanzee divergence

Background

Evolution of splice sites is a well-known phenomenon that results in transcript diversity during human evolution. Many novel splice sites are derived from repetitive elements and may not contribute to protein products. Here, we analyzed annotated human protein-coding exons and identified human-specific splice sites that arose after the human-chimpanzee divergence.

Results

We analyzed multiple alignments of the annotated human protein-coding exons and their respective orthologous mammalian genome sequences to identify 85 novel splice sites (50 splice acceptors and 35 donors) in the human genome. The novel protein-coding exons, which are expressed either constitutively or alternatively, produce novel protein isoforms by insertion, deletion, or frameshift. We found three cases in which the human-specific isoform conferred novel molecular function in the human cells: the human-specific IMUP protein isoform induces apoptosis of the trophoblast and is implicated in pre-eclampsia; the intronization of a part of SMOX gene exon produces inactive spermine oxidase; the human-specific NUB1 isoform shows reduced interaction with ubiquitin-like proteins, possibly affecting ubiquitin pathways.

Conclusions

Although the generation of novel protein isoforms does not equate to adaptive evolution, we propose that these cases are useful candidates for a molecular functional study to identify proteomic changes that might bring about novel phenotypes during human evolution.

Kinetochore KMN network gene CASC5 mutated in primary microcephaly

Several genes expressed at the centrosome or spindle pole have been reported to underlie autosomal recessive primary microcephaly (MCPH), a neurodevelopmental disorder consisting of an important brain size reduction present since birth, associated with mild-to-moderate mental handicap and no other neurological feature nor associated malformation. Here, we report a mutation of CASC5 (aka Blinkin, or KNL1, or hSPC105) in MCPH patients from three consanguineous families, in one of which we initially reported the MCPH4 locus. The combined logarithm of odds score of the three families was >6. All patients shared a very rare homozygous mutation of CASC5. The mutation induced skipping of exon 18 with subsequent frameshift and truncation of the predicted protein. CASC5 is part of the KMN network of the kinetochore and is required for proper microtubule attachment to the chromosome centromere and for spindle-assembly checkpoint (SAC) activation during mitosis. Like MCPH gene ASPM, CASC5 is upregulated in the ventricular zone (VZ) of the human fetal brain. CASC5 binds BUB1, BUBR1, ZWINT-1 and interestingly it binds to MIS12 through a protein domain which is truncated by the mutation. CASC5 localized at the equatorial plate like ZWINT-1 and BUBR1, while ASPM, CEP152 and PCTN localized at the spindle poles in our patients and in controls. Comparison of primate and rodent lineages indicates accelerated evolution of CASC5 in the human lineage. Our data provide strong evidence for CASC5 as a novel MCPH gene, and underscore the role of kinetochore integrity in proper volumetric development of the human brain.

Rewiring of PDZ domain-ligand interaction network contributed to eukaryotic evolution

PDZ domain-mediated interactions have greatly expanded during metazoan evolution, becoming important for controlling signal flow via the assembly of multiple signaling components. The evolutionary history of PDZ domain-mediated interactions has never been explored at the molecular level. It is of great interest to understand how PDZ domain-ligand interactions emerged and how they become rewired during evolution. Here, we constructed the first human PDZ domain-ligand interaction network (PDZNet) together with binding motif sequences and interaction strengths of ligands. PDZNet includes 1,213 interactions between 97 human PDZ proteins and 591 ligands that connect most PDZ protein-mediated interactions (98%) in a large single network via shared ligands. We examined the rewiring of PDZ domain-ligand interactions throughout eukaryotic evolution by tracing changes in the C-terminal binding motif sequences of the PDZ ligands. We found that interaction rewiring by sequence mutation frequently occurred throughout evolution, largely contributing to the growth of PDZNet. The rewiring of PDZ domain-ligand interactions provided an effective means of functional innovations in nervous system development. Our findings provide empirical evidence for a network evolution model that highlights the rewiring of interactions as a mechanism for the development of new protein functions. PDZNet will be a valuable resource to further characterize the organization of the PDZ domain-mediated signaling proteome.

Rewiring of interactions is a powerful tool for the evolution of organism complexity. Rewiring among preexisting proteins provides a simple mechanism for the development of new signaling circuits by redirecting information flows without a gain or loss of genes. Particularly, interactions mediated by short linear motifs can be easily changed by mutations during evolution, resulting in a rewiring of interactions. However, how interaction rewiring of linear motif interactions facilitates the emergence of new protein function during evolution is poorly understood. Here, we systematically investigated the rewiring of interactions mediated by PDZ domains, which are one of the most commonly found peptide recognition modules. We found that PDZ domain-ligand interactions are frequently rewired by C-terminal sequence mutations in PDZ ligands during evolution. Especially, rewiring of PDZ domain-ligand interactions was involved in neuronal function development, occurring concurrently with the emergence of vertebrates and suggesting that reorganization of signaling pathways by rewiring PDZ domain-ligand interactions significantly contributed to the evolution of nervous systems in vertebrates. Our findings highlight the rewiring of interactions as an effective means for functional innovation, providing new insight into eukaryotic evolution, which has not been fully explained by only the expansion of protein families.