Post-translational buffering leads to convergent protein expression levels between primates

Macroevolutionary divergence of gene expression driven by selection on protein abundance

The regulation of messenger RNA (mRNA) and protein abundances is well-studied, but less is known about the evolutionary processes shaping their relationship. To address this, we derived a new phylogenetic model and applied it to multispecies mammalian data. Our analyses reveal (i) strong stabilizing selection on protein abundances over macroevolutionary time, (ii) mutations affecting mRNA abundances minimally impact protein abundances, (iii) mRNA abundances evolve under selection to align with protein abundances, and (iv) mRNA abundances adapt faster than protein abundances owing to greater mutational opportunity. These conclusions are supported by comparisons of model parameters with independent functional genomic data. By decomposing mutational and selective influences on mRNA-protein dynamics, our approach provides a framework for discovering the evolutionary rules that drive divergence in gene expression.

Multilevel Gene Expression Changes in Lineages Containing Adaptive Copy Number Variants

Copy number variants (CNVs) are an important class of genetic variation that can mediate rapid adaptive evolution. Whereas, CNVs can increase the relative fitness of the organism, they can also incur a cost due to the associated increased gene expression and repetitive DNA. We previously evolved populations of Saccharomyces cerevisiae over hundreds of generations in glutamine-limited (Gln-) chemostats and observed the recurrent evolution of CNVs at the GAP1 locus. To understand the role that gene expression plays in adaptation, both in relation to the adaptation of the organism to the selective condition and as a consequence of the CNV, we measured the transcriptome, translatome, and proteome of 4 strains of evolved yeast, each with a unique CNV, and their ancestor in Gln- chemostats. We find CNV-amplified genes correlate with higher mRNA abundance; however, this effect is reduced at the level of the proteome, consistent with post-transcriptional dosage compensation. By normalizing each level of gene expression by the abundance of the preceding step we were able to identify widespread differences in the efficiency of each level of gene expression. Genes with significantly different translational efficiency were enriched for potential regulatory mechanisms including either upstream open reading frames, RNA-binding sites for Ssd1, or both. Genes with lower protein expression efficiency were enriched for genes encoding proteins in protein complexes. Taken together, our study reveals widespread changes in gene expression at multiple regulatory levels in lineages containing adaptive CNVs highlighting the diverse ways in which genome evolution shapes gene expression.

Evolution of translational control and the emergence of genes and open reading frames in human and non-human primate hearts

Evolutionary innovations can be driven by changes in the rates of RNA translation and the emergence of new genes and small open reading frames (sORFs). In this study, we characterized the transcriptional and translational landscape of the hearts of four primate and two rodent species through integrative ribosome and transcriptomic profiling, including adult left ventricle tissues and induced pluripotent stem cell-derived cardiomyocyte cell cultures. We show here that the translational efficiencies of subunits of the mitochondrial oxidative phosphorylation chain complexes IV and V evolved rapidly across mammalian evolution. Moreover, we discovered hundreds of species-specific and lineage-specific genomic innovations that emerged during primate evolution in the heart, including 551 genes, 504 sORFs and 76 evolutionarily conserved genes displaying human-specific cardiac-enriched expression. Overall, our work describes the evolutionary processes and mechanisms that have shaped cardiac transcription and translation in recent primate evolution and sheds light on how these can contribute to cardiac development and disease.

Switch-like Gene Expression Modulates Disease Susceptibility

A fundamental challenge in biomedicine is understanding the mechanisms predisposing individuals to disease. While previous research has suggested that switch-like gene expression is crucial in driving biological variation and disease susceptibility, a systematic analysis across multiple tissues is still lacking. By analyzing transcriptomes from 943 individuals across 27 tissues, we identified 1,013 switch-like genes. We found that only 31 (3.1%) of these genes exhibit switch-like behavior across all tissues. These universally switch-like genes appear to be genetically driven, with large exonic genomic structural variants explaining five (~18%) of them. The remaining switch-like genes exhibit tissue-specific expression patterns. Notably, tissue-specific switch-like genes tend to be switched on or off in unison within individuals, likely under the influence of tissue-specific master regulators, including hormonal signals. Among our most significant findings, we identified hundreds of concordantly switched-off genes in the stomach and vagina that are linked to gastric cancer (41-fold, p<10-4) and vaginal atrophy (44-fold, p<10-4), respectively. Experimental analysis of vaginal tissues revealed that low systemic levels of estrogen lead to a significant reduction in both the epithelial thickness and the expression of the switch-like gene ALOX12. We propose a model wherein the switching off of driver genes in basal and parabasal epithelium suppresses cell proliferation therein, leading to epithelial thinning and, therefore, vaginal atrophy. Our findings underscore the significant biomedical implications of switch-like gene expression and lay the groundwork for potential diagnostic and therapeutic applications.

Switch-like Gene Expression Modulates Disease Susceptibility

A fundamental challenge in biomedicine is understanding the mechanisms predisposing individuals to disease. While previous research has suggested that switch-like gene expression is crucial in driving biological variation and disease susceptibility, a systematic analysis across multiple tissues is still lacking. By analyzing transcriptomes from 943 individuals across 27 tissues, we identified 1,013 switch-like genes. We found that only 31 (3.1%) of these genes exhibit switch-like behavior across all tissues. These universally switch-like genes appear to be genetically driven, with large exonic genomic structural variants explaining five (~18%) of them. The remaining switch-like genes exhibit tissue-specific expression patterns. Notably, tissue-specific switch-like genes tend to be switched on or off in unison within individuals, likely under the influence of tissue-specific master regulators, including hormonal signals. Among our most significant findings, we identified hundreds of concordantly switched-off genes in the stomach and vagina that are linked to gastric cancer (41-fold, p<10-4) and vaginal atrophy (44-fold, p<10-4), respectively. Experimental analysis of vaginal tissues revealed that low systemic levels of estrogen lead to a significant reduction in both the epithelial thickness and the expression of the switch-like gene ALOX12. We propose a model wherein the switching off of driver genes in basal and parabasal epithelium suppresses cell proliferation therein, leading to epithelial thinning and, therefore, vaginal atrophy. Our findings underscore the significant biomedical implications of switch-like gene expression and lay the groundwork for potential diagnostic and therapeutic applications.

Comprehensive translational profiling and STE AI uncover rapid control of protein biosynthesis during cell stress

Translational control is important in all life, but it remains a challenge to accurately quantify. When ribosomes translate messenger (m)RNA into proteins, they attach to the mRNA in series, forming poly(ribo)somes, and can co-localize. Here, we computationally model new types of co-localized ribosomal complexes on mRNA and identify them using enhanced translation complex profile sequencing (eTCP-seq) based on rapid in vivo crosslinking. We detect long disome footprints outside regions of non-random elongation stalls and show these are linked to translation initiation and protein biosynthesis rates. We subject footprints of disomes and other translation complexes to artificial intelligence (AI) analysis and construct a new, accurate and self-normalized measure of translation, termed stochastic translation efficiency (STE). We then apply STE to investigate rapid changes to mRNA translation in yeast undergoing glucose depletion. Importantly, we show that, well beyond tagging elongation stalls, footprints of co-localized ribosomes provide rich insight into translational mechanisms, polysome dynamics and topology. STE AI ranks cellular mRNAs by absolute translation rates under given conditions, can assist in identifying its control elements and will facilitate the development of next-generation synthetic biology designs and mRNA-based therapeutics.

Multilevel gene expression changes in lineages containing adaptive copy number variants

Copy-number variants (CNVs) are an important class of recurrent variants that mediate adaptive evolution. While CNVs can increase the relative fitness of the organism, they can also incur a cost. We previously evolved populations of Saccharomyces cerevisiae over hundreds of generations in glutamine-limited (Gln-) chemostats and observed the recurrent evolution of CNVs at the GAP1 locus. To understand the role that expression plays in adaptation, both in relation to the adaptation of the organism to the selective condition, and as a consequence of the CNV, we measured the transcriptome, translatome, and proteome of 4 strains of evolved yeast, each with a unique CNV, and their ancestor in Gln- conditions. We find CNV-amplified genes correlate with higher RNA abundance; however, this effect is reduced at the level of the proteome, consistent with post-transcriptional dosage compensation. By normalizing each level of expression by the abundance of the preceding step we were able to identify widespread divergence in the efficiency of each step in the gene in the efficiency of each step in gene expression. Genes with significantly different translational efficiency were enriched for potential regulatory mechanisms including either upstream open reading frames, RNA binding sites for SSD1, or both. Genes with lower protein expression efficiency were enriched for genes encoding proteins in protein complexes. Taken together, our study reveals widespread changes in gene expression at multiple regulatory levels in lineages containing adaptive CNVs highlighting the diverse ways in which adaptive evolution shapes gene expression.

Translatomics to explore dynamic differences in immunocytes in the tumor microenvironment

Protein is an essential component of all living organisms and is primarily responsible for life activities; furthermore, its synthesis depends on a highly complex and accurate translation system. For proteins, the regulation at the translation level exceeds the sum of that during transcription, mRNA degradation, and protein degradation. Therefore, it is necessary to study regulation at the translation level. Imbalance in the translation process may change the cellular landscape, which not only leads to the occurrence, maintenance, progression, invasion, and metastasis of cancer but also affects the function of immune cells and changes the tumor microenvironment. Detailed analysis of transcriptional and protein atlases is needed to better understand how gene translation occurs. However, a more rigorous direct correlation between mRNA and protein levels is needed, which somewhat limits further studies. Translatomics is a technique for capturing and sequencing ribosome-related mRNAs that can effectively identify translation changes caused by ribosome stagnation and local translation abnormalities during cancer occurrence to further understand the changes in the translation landscape of cancer cells themselves and immune cells in the tumor microenvironment, which can provide new strategies and directions for tumor treatment.

On the Decoupling of Evolutionary Changes in mRNA and Protein Levels

Variation in gene expression across lineages is thought to explain much of the observed phenotypic variation and adaptation. The protein is closer to the target of natural selection but gene expression is typically measured as the amount of mRNA. The broad assumption that mRNA levels are good proxies for protein levels has been undermined by a number of studies reporting moderate or weak correlations between the two measures across species. One biological explanation for this discrepancy is that there has been compensatory evolution between the mRNA level and regulation of translation. However, we do not understand the evolutionary conditions necessary for this to occur nor the expected strength of the correlation between mRNA and protein levels. Here, we develop a theoretical model for the coevolution of mRNA and protein levels and investigate the dynamics of the model over time. We find that compensatory evolution is widespread when there is stabilizing selection on the protein level; this observation held true across a variety of regulatory pathways. When the protein level is under directional selection, the mRNA level of a gene and the translation rate of the same gene were negatively correlated across lineages but positively correlated across genes. These findings help explain results from comparative studies of gene expression and potentially enable researchers to disentangle biological and statistical hypotheses for the mismatch between transcriptomic and proteomic data.

Stable maternal proteins underlie distinct transcriptome, translatome, and proteome reprogramming during mouse oocyte-to-embryo transition

BACKGROUND

The oocyte-to-embryo transition (OET) converts terminally differentiated gametes into a totipotent embryo and is critically controlled by maternal mRNAs and proteins, while the genome is silent until zygotic genome activation. How the transcriptome, translatome, and proteome are coordinated during this critical developmental window remains poorly understood.

RESULTS

Utilizing a highly sensitive and quantitative mass spectrometry approach, we obtain high-quality proteome data spanning seven mouse stages, from full-grown oocyte (FGO) to blastocyst, using 100 oocytes/embryos at each stage. Integrative analyses reveal distinct proteome reprogramming compared to that of the transcriptome or translatome. FGO to 8-cell proteomes are dominated by FGO-stockpiled proteins, while the transcriptome and translatome are more dynamic. FGO-originated proteins frequently persist to blastocyst while corresponding transcripts are already downregulated or decayed. Improved concordance between protein and translation or transcription is observed for genes starting translation upon meiotic resumption, as well as those transcribed and translated only in embryos. Concordance between protein and transcription/translation is also observed for proteins with short half-lives. We built a kinetic model that predicts protein dynamics by incorporating both initial protein abundance in FGOs and translation kinetics across developmental stages.

CONCLUSIONS

Through integrative analyses of datasets generated by ultrasensitive methods, our study reveals that the proteome shows distinct dynamics compared to the translatome and transcriptome during mouse OET. We propose that the remarkably stable oocyte-originated proteome may help save resources to accommodate the demanding needs of growing embryos. This study will advance our understanding of mammalian OET and the fundamental principles governing gene expression.

Novel start codons introduce novel coding sequences in the human genomes

Start-gain mutations can introduce novel start codons and generate novel coding sequences that may affect the function of genes. In this study, we systematically investigated the novel start codons that were either polymorphic or fixed in the human genomes. 829 polymorphic start-gain SNVs were identified in the human populations, and the novel start codons introduced by these SNVs have significantly higher activity in translation initiation. Some of these start-gain SNVs were reported to be associated with phenotypes and diseases in previous studies. By comparative genomic analysis, we found 26 human-specific start codons that were fixed after the divergence between the human and chimpanzee, and high-level translation initiation activity was observed on them. The negative selection signal was detected in the novel coding sequences introduced by these human-specific start codons, indicating the important function of these novel coding sequences.

Evolutionary trade-off and mutational bias could favor transcriptional over translational divergence within paralog pairs

How changes in the different steps of protein synthesis-transcription, translation and degradation-contribute to differences of protein abundance among genes is not fully understood. There is however accumulating evidence that transcriptional divergence might have a prominent role. Here, we show that yeast paralogous genes are more divergent in transcription than in translation. We explore two causal mechanisms for this predominance of transcriptional divergence: an evolutionary trade-off between the precision and economy of gene expression and a larger mutational target size for transcription. Performing simulations within a minimal model of post-duplication evolution, we find that both mechanisms are consistent with the observed divergence patterns. We also investigate how additional properties of the effects of mutations on gene expression, such as their asymmetry and correlation across levels of regulation, can shape the evolution of paralogs. Our results highlight the importance of fully characterizing the distributions of mutational effects on transcription and translation. They also show how general trade-offs in cellular processes and mutation bias can have far-reaching evolutionary impacts.

Genetic architecture of human brain evolution

Comparative studies of hominids have long sought to identify mutational events that shaped the evolution of the human nervous system. However, functional genetic differences are outnumbered by millions of nearly neutral mutations, and the developmental mechanisms underlying human nervous system specializations are difficult to model and incompletely understood. Candidate-gene studies have attempted to map select human-specific genetic differences to neurodevelopmental functions, but it remains unclear how to contextualize the relative effects of genes that are investigated independently. Considering these limitations, we discuss scalable approaches for probing the functional contributions of human-specific genetic differences. We propose that a systems-level view will enable a more quantitative and integrative understanding of the genetic, molecular and cellular underpinnings of human nervous system evolution.

Current advances in primate genomics: novel approaches for understanding evolution and disease

Primate genomics holds the key to understanding fundamental aspects of human evolution and disease. However, genetic diversity and functional genomics data sets are currently available for only a few of the more than 500 extant primate species. Concerted efforts are under way to characterize primate genomes, genetic polymorphism and divergence, and functional landscapes across the primate phylogeny. The resulting data sets will enable the connection of genotypes to phenotypes and provide new insight into aspects of the genetics of primate traits, including human diseases. In this Review, we describe the existing genome assemblies as well as genetic variation and functional genomic data sets. We highlight some of the challenges with sample acquisition. Finally, we explore how technological advances in single-cell functional genomics and induced pluripotent stem cell-derived organoids will facilitate our understanding of the molecular foundations of primate biology.

riboviz 2: a flexible and robust ribosome profiling data analysis and visualization workflow

MOTIVATION

Ribosome profiling, or Ribo-seq, is the state-of-the-art method for quantifying protein synthesis in living cells. Computational analysis of Ribo-seq data remains challenging due to the complexity of the procedure, as well as variations introduced for specific organisms or specialized analyses.

RESULTS

We present riboviz 2, an updated riboviz package, for the comprehensive transcript-centric analysis and visualization of Ribo-seq data. riboviz 2 includes an analysis workflow built on the Nextflow workflow management system for end-to-end processing of Ribo-seq data. riboviz 2 has been extensively tested on diverse species and library preparation strategies, including multiplexed samples. riboviz 2 is flexible and uses open, documented file formats, allowing users to integrate new analyses with the pipeline.

AVAILABILITY AND IMPLEMENTATION

riboviz 2 is freely available at github.com/riboviz/riboviz.

Selection of Cashmere Fineness Functional Genes by Translatomics

Cashmere fineness is an important index to evaluate cashmere quality. Liaoning Cashmere Goat (LCG) has a large cashmere production and long cashmere fiber, but its fineness is not ideal. Therefore, it is important to find genes involved in cashmere fineness that can be used in future endeavors aiming to improve this phenotype. With the continuous advancement of research, the regulation of cashmere fineness has made new developments through high-throughput sequencing and genome-wide association analysis. It has been found that translatomics can identify genes associated with phenotypic traits. Through translatomic analysis, the skin tissue of LCG sample groups differing in cashmere fineness was sequenced by Ribo-seq. With these data, we identified 529 differentially expressed genes between the sample groups among the 27197 expressed genes. From these, 343 genes were upregulated in the fine LCG group in relation to the coarse LCG group, and 186 were downregulated in the same relationship. Through GO enrichment analysis and KEGG enrichment analysis of differential genes, the biological functions and pathways of differential genes can be found. In the GO enrichment analysis, 491 genes were significantly enriched, and the functional region was mainly in the extracellular region. In the KEGG enrichment analysis, the enrichment of the human papillomavirus infection pathway was seen the most. We found that the COL6A5 gene may affect cashmere fineness.

CRYPTOCHROMES promote daily protein homeostasis

The daily organisation of most mammalian cellular functions is attributed to circadian regulation of clock-controlled protein expression, driven by daily cycles of CRYPTOCHROME-dependent transcriptional feedback repression. To test this, we used quantitative mass spectrometry to compare wild-type and CRY-deficient fibroblasts under constant conditions. In CRY-deficient cells, we found that temporal variation in protein, phosphopeptide, and K+ abundance was at least as great as wild-type controls. Most strikingly, the extent of temporal variation within either genotype was much smaller than overall differences in proteome composition between WT and CRY-deficient cells. This proteome imbalance in CRY-deficient cells and tissues was associated with increased susceptibility to proteotoxic stress, which impairs circadian robustness, and may contribute to the wide-ranging phenotypes of CRY-deficient mice. Rather than generating large-scale daily variation in proteome composition, we suggest it is plausible that the various transcriptional and post-translational functions of CRY proteins ultimately act to maintain protein and osmotic homeostasis against daily perturbation.

CRISPR-Cas Gene Perturbation and Editing in Human Induced Pluripotent Stem Cells

Directing the fates of human pluripotent stem cells (hPSC) to generate a multitude of differentiated cell types allows the study of the genetic regulation of human development and disease. The translational potential of hPSC is maximized by exploiting CRISPR to silence or activate genes with spatial and temporal precision permanently or reversibly. Here, we summarize the increasingly refined and diverse CRISPR toolkit for the latter forms of gene perturbation in hPSC and their downstream applications. We discuss newer methods to install edits efficiently with single nucleotide resolution and describe pooled CRISPR screens as a powerful means of unbiased discovery of genes associated with a phenotype of interest. Last, we discuss the potential of these combined technologies in the treatment of hitherto intractable human diseases and the challenges to their implementation in the clinic.

Assessing the functional relevance of splice isoforms

Alternative splicing of messenger RNA can generate an array of mature transcripts, but it is not clear how many go on to produce functionally relevant protein isoforms. There is only limited evidence for alternative proteins in proteomics analyses and data from population genetic variation studies indicate that most alternative exons are evolving neutrally. Determining which transcripts produce biologically important isoforms is key to understanding isoform function and to interpreting the real impact of somatic mutations and germline variations. Here we have developed a method, TRIFID, to classify the functional importance of splice isoforms. TRIFID was trained on isoforms detected in large-scale proteomics analyses and distinguishes these biologically important splice isoforms with high confidence. Isoforms predicted as functionally important by the algorithm had measurable cross species conservation and significantly fewer broken functional domains. Additionally, exons that code for these functionally important protein isoforms are under purifying selection, while exons from low scoring transcripts largely appear to be evolving neutrally. TRIFID has been developed for the human genome, but it could in principle be applied to other well-annotated species. We believe that this method will generate valuable insights into the cellular importance of alternative splicing.

Human-chimpanzee fused cells reveal cis-regulatory divergence underlying skeletal evolution

Gene regulatory divergence is thought to play a central role in determining human-specific traits. However, our ability to link divergent regulation to divergent phenotypes is limited. Here, we utilized human-chimpanzee hybrid induced pluripotent stem cells to study gene expression separating these species. The tetraploid hybrid cells allowed us to separate cis- from trans-regulatory effects, and to control for non-genetic confounding factors. We differentiated these cells into cranial neural crest cells (CNCCs), the primary cell type giving rise to the face. We discovered evidence of lineage-specific selection on the hedgehog signaling pathway, including a human-specific 6-fold down-regulation of EVC2 (LIMBIN), a key hedgehog gene. Inducing a similar down-regulation of EVC2 substantially reduced hedgehog signaling output. Mice and humans lacking functional EVC2 show striking phenotypic parallels to human-chimpanzee craniofacial differences, suggesting that the regulatory divergence of hedgehog signaling may have contributed to the unique craniofacial morphology of humans.

Divergence in alternative polyadenylation contributes to gene regulatory differences between humans and chimpanzees

While comparative functional genomic studies have shown that inter-species differences in gene expression can be explained by corresponding inter-species differences in genetic and epigenetic regulatory mechanisms, co-transcriptional mechanisms, such as alternative polyadenylation (APA), have received little attention. We characterized APA in lymphoblastoid cell lines from six humans and six chimpanzees by identifying and estimating the usage for 44,432 polyadenylation sites (PAS) in 9518 genes. Although APA is largely conserved, 1705 genes showed significantly different PAS usage (FDR 0.05) between species. Genes with divergent APA also tend to be differentially expressed, are enriched among genes showing differences in protein translation, and can explain a subset of observed inter-species protein expression differences that do not differ at the transcript level. Finally, we found that genes with a dominant PAS, which is used more often than other PAS, are particularly enriched for differentially expressed genes.

Determinants of genome-wide distribution and evolution of uORFs in eukaryotes

Upstream open reading frames (uORFs) play widespread regulatory functions in modulating mRNA translation in eukaryotes, but the principles underlying the genomic distribution and evolution of uORFs remain poorly understood. Here, we analyze ~17 million putative canonical uORFs in 478 eukaryotic species that span most of the extant taxa of eukaryotes. We demonstrate how positive and purifying selection, coupled with differences in effective population size (Ne), has shaped the contents of uORFs in eukaryotes. Besides, gene expression level is important in influencing uORF occurrences across genes in a species. Our analyses suggest that most uORFs might play regulatory roles rather than encode functional peptides. We also show that the Kozak sequence context of uORFs has evolved across eukaryotic clades, and that noncanonical uORFs tend to have weaker suppressive effects than canonical uORFs in translation regulation. This study provides insights into the driving forces underlying uORF evolution in eukaryotes.

Isoform-resolved correlation analysis between mRNA abundance regulation and protein level degradation

Profiling of biological relationships between different molecular layers dissects regulatory mechanisms that ultimately determine cellular function. To thoroughly assess the role of protein post‐translational turnover, we devised a strategy combining pulse stable isotope‐labeled amino acids in cells (pSILAC), data‐independent acquisition mass spectrometry (DIA‐MS), and a novel data analysis framework that resolves protein degradation rate on the level of mRNA alternative splicing isoforms and isoform groups. We demonstrated our approach by the genome‐wide correlation analysis between mRNA amounts and protein degradation across different strains of HeLa cells that harbor a high grade of gene dosage variation. The dataset revealed that specific biological processes, cellular organelles, spatial compartments of organelles, and individual protein isoforms of the same genes could have distinctive degradation rate. The protein degradation diversity thus dissects the corresponding buffering or concerting protein turnover control across cancer cell lines. The data further indicate that specific mRNA splicing events such as intron retention significantly impact the protein abundance levels. Our findings support the tight association between transcriptome variability and proteostasis and provide a methodological foundation for studying functional protein degradation.

Transcriptome and translatome co-evolution in mammals

Gene expression programs define shared and species-specific phenotypes, but their evolution remains largely uncharacterized beyond the transcriptome layer¹. Here we report an analysis of the co-evolution of translatomes and transcriptomes using ribosome-profling and matched RNA-sequencing data for three organs (brain, liver and testis) in fve mammals (human, macaque, mouse, opossum and platypus) and a bird (chicken). Our within-species analyses reveal that translational regulation is widespread in the diferent organs, in particular across the spermatogenic cell types of the testis. The between-species divergence in gene expression is around 20% lower at the translatome layer than at the transcriptome layer owing to extensive buffering between the expression layers, which especially preserved old, essential and housekeeping genes. Translational upregulation specifcally counterbalanced global dosage reductions during the evolution of sex chromosomes and the efects of meiotic sex-chromosome inactivation during spermatogenesis. Despite the overall prevalence of bufering, some genes evolved faster at the translatome layer—potentially indicating adaptive changes in expression; testis tissue shows the highest fraction of such genes. Further analyses incorporating mass spectrometry proteomics data establish that the co-evolution of transcriptomes and translatomes is refected at the proteome layer. Together, our work uncovers co-evolutionary patterns and associated selective forces across the expression layers, and provides a resource for understanding their interplay in mammalian organs.

Human and Nonhuman Primate Lineage-Specific Footprints in the Salivary Proteome

Proteins in saliva are needed for preprocessing food in the mouth, maintenance of tooth mineralization, and protection from microbial pathogens. Novel insights into human lineage-specific functions of salivary proteins and clues to their involvement in human disease can be gained through evolutionary studies, as recently shown for salivary amylase AMY1 and salivary agglutinin DMBT1/gp340. However, the entirety of proteins in saliva, the salivary proteome, has not yet been investigated from an evolutionary perspective. Here, we compared the proteomes of human saliva and the saliva of our closest extant evolutionary relatives, chimpanzees and gorillas, using macaques as an outgroup, with the aim to uncover features in saliva protein composition that are unique to each species. We found that humans produce a waterier saliva, containing less than half total protein than great apes and Old World monkeys. For all major salivary proteins in humans, we could identify counterparts in chimpanzee and gorilla saliva. However, we discovered unique protein profiles in saliva of humans that were distinct from those of nonhuman primates. These findings open up the possibility that dietary differences and pathogenic pressures may have shaped a distinct salivary proteome in the human lineage.

An analysis of tissue-specific alternative splicing at the protein level

The role of alternative splicing is one of the great unanswered questions in cellular biology. There is strong evidence for alternative splicing at the transcript level, and transcriptomics experiments show that many splice events are tissue specific. It has been suggested that alternative splicing evolved in order to remodel tissue-specific protein-protein networks. Here we investigated the evidence for tissue-specific splicing among splice isoforms detected in a large-scale proteomics analysis. Although the data supporting alternative splicing is limited at the protein level, clear patterns emerged among the small numbers of alternative splice events that we could detect in the proteomics data. More than a third of these splice events were tissue-specific and most were ancient: over 95% of splice events that were tissue-specific in both proteomics and RNAseq analyses evolved prior to the ancestors of lobe-finned fish, at least 400 million years ago. By way of contrast, three in four alternative exons in the human gene set arose in the primate lineage, so our results cannot be extrapolated to the whole genome. Tissue-specific alternative protein forms in the proteomics analysis were particularly abundant in nervous and muscle tissues and their genes had roles related to the cytoskeleton and either the structure of muscle fibres or cell-cell connections. Our results suggest that this conserved tissue-specific alternative splicing may have played a role in the development of the vertebrate brain and heart.

We manually curated a set of 255 splice events detected in a large-scale tissue-based proteomics experiment and found that more than a third had evidence of significant tissue-specific differences. Events that were significantly tissue-specific at the protein level were highly conserved; almost 75% evolved over 400 million years ago. The tissues in which we found most evidence for tissue-specific splicing were nervous tissues and cardiac tissues. Genes with tissue-specific events in these two tissues had functions related to important cellular structures in brain and heart tissues. These splice events may have been essential for the development of vertebrate heart and muscle. However, our data set may not be representative of alternative exons as a whole. We found that most tissue specific splicing was strongly conserved, but just 5% of annotated alternative exons in the human gene set are ancient. More than three quarters of alternative exons are primate-derived. Although the analysis does not provide a definitive answer to the question of the functional role of alternative splicing, our results do indicate that alternative splice variants may have played a significant part in the evolution of brain and heart tissues in vertebrates.

A hierarchical Bayesian mixture model for inferring the expression state of genes in transcriptomes

How do the cells of an organism—each with an identical genome—give rise to tissues of incredible phenotypic diversity? Key to answering this question is the transcriptome: the set of genes expressed in a given tissue. We would clearly benefit from the ability to identify qualitative differences in expression (whether a gene is active or inactive in a given tissue/species). Inferring the expression state of genes is surprisingly difficult, owing to the complex biological processes that give rise to transcriptomes and to the vagaries of techniques used to generate transcriptomic datasets. We develop a hierarchical Bayesian mixture model that—by describing those biological and technical processes—allows us to infer the expression state of genes from replicate transcriptomic datasets.

Transcriptomes are key to understanding the relationship between genotype and phenotype. The ability to infer the expression state (active or inactive) of genes in the transcriptome offers unique benefits for addressing this issue. For example, qualitative changes in gene expression may underly the origin of novel phenotypes, and expression states are readily comparable between tissues and species. However, inferring the expression state of genes is a surprisingly difficult problem, owing to the complex biological and technical processes that give rise to observed transcriptomic datasets. Here, we develop a hierarchical Bayesian mixture model that describes this complex process and allows us to infer expression state of genes from replicate transcriptomic libraries. We explore the statistical behavior of this method with analyses of simulated datasets—where we demonstrate its ability to correctly infer true (known) expression states—and empirical-benchmark datasets, where we demonstrate that the expression states inferred from RNA-sequencing (RNA-seq) datasets using our method are consistent with those based on independent evidence. The power of our method to correctly infer expression states is generally high and remarkably, approaches the maximum possible power for this inference problem. We present an empirical analysis of primate-brain transcriptomes, which identifies genes that have a unique expression state in humans. Our method is implemented in the freely available R package zigzag.

The translational landscape of ground state pluripotency

Translational control plays a central role in regulation of gene expression and can lead to significant divergence between mRNA- and protein-abundance. Here, we used genome-wide approaches combined with time-course analysis to measure the mRNA-abundance, mRNA-translation rate and protein expression during the transition of naïve-to-primed mouse embryonic stem cells (ESCs). We find that the ground state ESCs cultured with GSK3-, MEK-inhibitors and LIF (2iL) display higher ribosome density on a selective set of mRNAs. This set of mRNAs undergo strong translational buffering to maintain stable protein expression levels in 2iL-ESCs. Importantly, we show that the global alteration of cellular proteome during the transition of naïve-to-primed pluripotency is largely accompanied by transcriptional rewiring. Thus, we provide a comprehensive and detailed overview of the global changes in gene expression in different states of ESCs and dissect the relative contributions of mRNA-transcription, translation and regulation of protein stability in controlling protein abundance.

Translational control of gene expression can lead to significant divergence between mRNA and protein abundance. Here, the authors describe transcriptional rewiring and translational buffering during transition from naïve to primed pluripotency through quantitation of mRNA-abundance, translation rate and protein expression.

Neurology-related protein biomarkers are associated with cognitive ability and brain volume in older age

Identifying biological correlates of late life cognitive function is important if we are to ascertain biomarkers for, and develop treatments to help reduce, age-related cognitive decline. Here, we investigated the associations between plasma levels of 90 neurology-related proteins (Olink® Proteomics) and general fluid cognitive ability in the Lothian Birth Cohort 1936 (LBC1936, N = 798), Lothian Birth Cohort 1921 (LBC1921, N = 165), and the INTERVAL BioResource (N = 4451). In the LBC1936, 22 of the proteins were significantly associated with general fluid cognitive ability (β between -0.11 and -0.17). MRI-assessed total brain volume partially mediated the association between 10 of these proteins and general fluid cognitive ability. In an age-matched subsample of INTERVAL, effect sizes for the 22 proteins, although smaller, were all in the same direction as in LBC1936. Plasma levels of a number of neurology-related proteins are associated with general fluid cognitive ability in later life, mediated by brain volume in some cases.