Determinants of Protein Abundance and Translation Efficiency in S. cerevisiae

The translation efficiency of most Saccharomyces cerevisiae genes remains fairly constant across poor and rich growth media. This observation has led us to revisit the available data and to examine the potential utility of a protein abundance predictor in reinterpreting existing mRNA expression data. Our predictor is based on large-scale data of mRNA levels, the tRNA adaptation index, and the evolutionary rate. It attains a correlation of 0.76 with experimentally determined protein abundance levels on unseen data and successfully cross-predicts protein abundance levels in another yeast species (Schizosaccharomyces pombe). The predicted abundance levels of proteins in known S. cerevisiae complexes, and of interacting proteins, are significantly more coherent than their corresponding mRNA expression levels. Analysis of gene expression measurement experiments using the predicted protein abundance levels yields new insights that are not readily discernable when clustering the corresponding mRNA expression levels. Comparing protein abundance levels across poor and rich media, we find a general trend for homeostatic regulation where transcription and translation change in a reciprocal manner. This phenomenon is more prominent near origins of replications. Our analysis shows that in parallel to the adaptation occurring at the tRNA level via the codon bias, proteins do undergo a complementary adaptation at the amino acid level to further increase their abundance

Human cachexia induces changes in mitochondria, autophagy and apoptosis in the skeletal muscle

Cachexia is a wasting syndrome characterized by the continuous loss of skeletal muscle mass due to imbalance between protein synthesis and degradation, which is related with poor prognosis and compromised quality of life. Dysfunctional mitochondria are associated with lower muscle strength and muscle atrophy in cancer patients, yet poorly described in human cachexia. We herein investigated mitochondrial morphology, autophagy and apoptosis in the skeletal muscle of patients with gastrointestinal cancer-associated cachexia (CC), as compared with a weight-stable cancer group (WSC). CC showed prominent weight loss and increased circulating levels of serum C-reactive protein, lower body mass index and decreased circulating hemoglobin, when compared to WSC. Electron microscopy analysis revealed an increase in intermyofibrillar mitochondrial area in CC, as compared to WSC. Relative gene expression of Fission 1, a protein related to mitochondrial fission, was increased in CC, as compared to WSC. LC3 II, autophagy-related (ATG) 5 and 7 essential proteins for autophagosome formation, presented higher content in the cachectic group. Protein levels of phosphorylated p53 (Ser46), activated caspase 8 (Asp384) and 9 (Asp315) were also increased in the skeletal muscle of CC. Overall, our results demonstrate that human cancer-associated cachexia leads to exacerbated muscle-stress response that may culminate in muscle loss, which is in part due to disruption of mitochondrial morphology, dysfunctional autophagy and increased apoptosis. To the best of our knowledge, this is the first report showing quantitative morphological alterations in skeletal muscle mitochondria in cachectic patients

A signature of saliva-derived exosomal small RNAs as predicting biomarker for esophageal carcinoma:a multicenter prospective study

BACKGROUND: The tRNA-derived small RNAs (tsRNAs) are produced in a nuclease-dependent manner in responses to variety of stresses that are common in cancers. We focus on a cancer-enriched tsRNA signature to develop a salivary exosome-based non-invasive biomarker for human esophageal squamous cell carcinoma (ESCC). METHODS: Cancer-enriched small RNAs were identified by RNA sequencing of salivary exosomes obtained from ESCC patients (n = 3) and healthy controls (n = 3) in a pilot study and further validated in discovery cohort (n = 66). A multicenter prospective observational study was conducted in two ESCC high-incidence regions (n = 320 and 200, respectively) using the newly developed biomarker signature. RESULTS: The tsRNA (tRNA-GlyGCC-5) and a previously undocumented small RNA were specifically enriched in salivary exosomes of ESCC patients, ESCC tissues and ESCC cells. The bi-signature composed of these small RNAs was able to discriminate ESCC patients from the controls with high sensitivity (90.50%) and specificity (94.20%). Based on the bi-signature Risk Score for Prognosis (RSP), patients with high-RSP have both shorter overall survival (OS) (HR 4.95, 95%CI 2.90–8.46) and progression-free survival (PFS) (HR 3.69, 95%CI 2.24–6.10) than those with low-RSP. In addition, adjuvant therapy improved OS (HR 0.47, 95%CI 0.29–0.77) and PFS (HR 0.36, 95%CI 0.21–0.62) only for patients with high but not low RSP. These findings are consistent in both training and validation cohort. CONCLUSIONS: The tsRNA-based signature not only has the potential for diagnosis and prognosis but also may serve as a pre-operative biomarker to select patients who would benefit from adjuvant therapy. TRIAL REGISTRATION: A prospective study of diagnosis biomarkers of esophageal squamous cell carcinoma, ChiCTR2000031507. Registered 3 April 2016 - Retrospectively registered. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s12943-022-01499-8

A signature of saliva-derived exosomal small RNAs as predicting biomarker for esophageal carcinoma:a multicenter prospective study

BACKGROUND: The tRNA-derived small RNAs (tsRNAs) are produced in a nuclease-dependent manner in responses to variety of stresses that are common in cancers. We focus on a cancer-enriched tsRNA signature to develop a salivary exosome-based non-invasive biomarker for human esophageal squamous cell carcinoma (ESCC). METHODS: Cancer-enriched small RNAs were identified by RNA sequencing of salivary exosomes obtained from ESCC patients (n = 3) and healthy controls (n = 3) in a pilot study and further validated in discovery cohort (n = 66). A multicenter prospective observational study was conducted in two ESCC high-incidence regions (n = 320 and 200, respectively) using the newly developed biomarker signature. RESULTS: The tsRNA (tRNA-GlyGCC-5) and a previously undocumented small RNA were specifically enriched in salivary exosomes of ESCC patients, ESCC tissues and ESCC cells. The bi-signature composed of these small RNAs was able to discriminate ESCC patients from the controls with high sensitivity (90.50%) and specificity (94.20%). Based on the bi-signature Risk Score for Prognosis (RSP), patients with high-RSP have both shorter overall survival (OS) (HR 4.95, 95%CI 2.90–8.46) and progression-free survival (PFS) (HR 3.69, 95%CI 2.24–6.10) than those with low-RSP. In addition, adjuvant therapy improved OS (HR 0.47, 95%CI 0.29–0.77) and PFS (HR 0.36, 95%CI 0.21–0.62) only for patients with high but not low RSP. These findings are consistent in both training and validation cohort. CONCLUSIONS: The tsRNA-based signature not only has the potential for diagnosis and prognosis but also may serve as a pre-operative biomarker to select patients who would benefit from adjuvant therapy. TRIAL REGISTRATION: A prospective study of diagnosis biomarkers of esophageal squamous cell carcinoma, ChiCTR2000031507. Registered 3 April 2016 - Retrospectively registered. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s12943-022-01499-8

Single-stranded genomic architecture constrains optimal codon usage

Viral codon usage is shaped by the conflicting forces of mutational pressure and selection to match host patterns for optimal expression. We examined whether genomic architecture (single- or double-stranded DNA) influences the degree to which bacteriophage codon usage differ from their primary bacterial hosts and each other. While both correlated equally with their hosts' genomic nucleotide content, the coat genes of ssDNA phages were less well adapted than those of dsDNA phages to their hosts' codon usage profiles due to their preference for codons ending in thymine. No specific biases were detected in dsDNA phage genomes. In all nine of ten cases of codon redundancy in which a specific codon was overrepresented, ssDNA phages favored the NNT codon. A cytosine to thymine biased mutational pressure working in conjunction with strong selection against non-synonymous mutations appears be shaping codon usage bias in ssDNA viral genomes

The Evolutionary and Functional Roles of Synonymous Codon Usage in Eukaryotes

Most amino acids are encoded by multiple synonymous codons. Although alternative usage of synonymous codons does not affect the amino acid sequences of proteins, researchers have been reporting evidence for functional synonymous codon usage at the species- and gene-specific levels for over four decades. It has been shown that variations in synonymous codon usage can affect phenotypes through diverse mechanisms such as shaping translation efficiency and mRNA stability. On the other hand, the common view that cellular and organismal phenotypes are primarily determined by proteins whose functions are primarily determined by amino acid sequences, often drives the assumption that synonymous mutations are evolutionarily neutral. Consequently, this assumption has been used extensively in evolutionary biology, population genetics, and structural biology. One explanation of the apparent contradiction between the empirical findings, which indicate that synonymous mutations can affect related phenotypes, and the theoretical models, which stipulate that synonymous mutations are neutral, is that neutral synonymous mutations represent the general rule while non-neutral synonymous mutations represent the rare exceptions. In my thesis, I examined this explanation by applying computational and experimental approaches, which indicated that: 1) Non-neutral synonymous mutations significantly affect a considerable proportion of protein-coding genes; 2) Gene-specific codon usage patterns, such as the preference for a specific combination of rare codons, are possibly associated with specific gene functions, such as enhancing tissue-specific gene expression; 3) Some protein-coding genes include codon clusters whose codon usage patterns cannot be explained by selection-independent processes, and thus such codon clusters seem to serve as domains affecting protein functions. Together, these data suggest that synonymous mutations should not be a priori considered neutral. Furthermore, my studies suggest that the biochemical functions of at least some proteins are not only shaped by the constituent amino acid residues but also by codon usage biases at the gene-specific and sub-genic levels. In conclusion, my thesis work suggests that many of the commonly used approaches for analyzing the selection on protein-coding DNA sequences, which rely on the assumption that synonymous mutations are generally neutral, may generate biased results. Furthermore, my studies indicate that selection on gene-specific codon usage bias has evolved to serve diverse biological functions, which are still mostly uncharacterized

Multivariate Analysis of Prokaryotic Amino Acid Usage Bias: A Computational Method for Understanding Protein Building Block Selection in Primitive Organisms

Organisms expend a significant fraction of their overall energy budget in the creation of proteins, particularly for those that are produced in large quantities. Recent research has demonstrated that genes encoding these proteins are shaped by natural selection to produce the proteins with low cost building blocks (amino acids) whenever possible. The negative correlation between protein production rate and their energetic costs has been established for two bacterial genomes: Escherichia coli and Bacillus subtilis. This thesis provides scientific validation of this theory by automating the analysis and extending the research to additional genomes. Investigations into building block selection are highly computational in nature. Diverse methodologies, including principal component analysis, calculation of Mahalanobis distance, and the execution of Mantel-Haenszel and Bonferroni tests, are required in order to automate the process. In order to verify that the cause of the observed trend is energetic cost minimization it is necessary to eliminate as many alternative explanations as possible. This is accomplished through demonstration that the trend is not localized to any particular region of the protein’s primary structure and that the trend is consistent across all genes regardless of functionality. This investigation of the energetic cost of polypeptide synthesis provides valuable insights into protein building block selection. As an example, parasitic organisms appear to exhibit no correlation between protein production rate and amino acid cost. When the costs associated with building blocks that the parasite obtains from its host are removed,however, a trend once again becomes evident

Maintaining protein localization, structure, and functional interactions via codon usage and coevolution of gene expression: Combining evolutionary bioinformatics with omics-scale data to test hypotheses related to protein function

A major challenge of the omics-era is identifying how a protein functions, both in terms of its specific function and within the context of the various biological processes necessary for the cell\u27s survival. Key elements necessary for a protein to perform its function are efficient and accurate protein localization, protein folding, and interactions with other proteins. Previous work implicated codon usage as a means to modulate protein localization and folding. Using a mechanistic model rooted in population genetics, I examine potential selective differences in codon usage in signal peptides (localization) and protein secondary structures. Although previous work argued signal peptides were under selection for increased translation inefficiency, I find selection is generally consistent with the 5\u27-regions of non-secreted proteins. I also find that previous work was likely confounded by biases in signal peptide amino acid usage and gene expression. Although the direction of selection on codon usage is mostly consistent between protein secondary structures, the strength of this selection does vary for certain codons. After successful folding and localization of a protein, it must be able to function within the context of other proteins in the cell, often through protein-protein interactions of metabolic pathways. Previous work suggests proteins which are part of the same functional processes within a cell are co-expressed across time and environmental conditions. Using the concept of guilt-by-association, I combine empirical protein abundances (measured via mass spectrometry) with sequence homology based function prediction tools to identify potential functions of proteins of unknown function in \textit{C. thermocellum}. Building upon the concept that functionally-related genes are co-expressed within a species, I demonstrate how phylogenetic comparative methods can be used to detect signals of gene expression coevolution across species while accounting for the shared ancestry of the species in question

Amino Acid Usage Is Asymmetrically Biased in AT- and GC-Rich Microbial Genomes.

INTRODUCTION: Genomic base composition ranges from less than 25% AT to more than 85% AT in prokaryotes. Since only a small fraction of prokaryotic genomes is not protein coding even a minor change in genomic base composition will induce profound protein changes. We examined how amino acid and codon frequencies were distributed in over 2000 microbial genomes and how these distributions were affected by base compositional changes. In addition, we wanted to know how genome-wide amino acid usage was biased in the different genomes and how changes to base composition and mutations affected this bias. To carry this out, we used a Generalized Additive Mixed-effects Model (GAMM) to explore non-linear associations and strong data dependences in closely related microbes; principal component analysis (PCA) was used to examine genomic amino acid- and codon frequencies, while the concept of relative entropy was used to analyze genomic mutation rates. RESULTS: We found that genomic amino acid frequencies carried a stronger phylogenetic signal than codon frequencies, but that this signal was weak compared to that of genomic %AT. Further, in contrast to codon usage bias (CUB), amino acid usage bias (AAUB) was differently distributed in AT- and GC-rich genomes in the sense that AT-rich genomes did not prefer specific amino acids over others to the same extent as GC-rich genomes. AAUB was also associated with relative entropy; genomes with low AAUB contained more random mutations as a consequence of relaxed purifying selection than genomes with higher AAUB. CONCLUSION: Genomic base composition has a substantial effect on both amino acid- and codon frequencies in bacterial genomes. While phylogeny influenced amino acid usage more in GC-rich genomes, AT-content was driving amino acid usage in AT-rich genomes. We found the GAMM model to be an excellent tool to analyze the genomic data used in this study

Evolutionary dynamics promoting and accompanying rapid adaptive trait loss

Understanding the conditions which promote adaptation is a key goal of evolutionary biology, and a pressing issue across fields of biology. Addressing this involves investigating not just the genetic and developmental mechanisms through which adaptive phenotypes arise, but also the environmental and ecological conditions which promote their spread. A major challenge in addressing these aims is that contemporary examples of rapid adaptive evolution are difficult to study, owing to the difficulty of identifying traits under selection during the early stages of adaptation. In this thesis, I use a Hawaiian field cricket system which provides a useful exception; males of the species Teleogryllus oceanicus ordinarily sing to attract females, however adaptive male-silencing (‘flatwing’) phenotypes have recently emerged and spread on at least three islands, under selection against male song exerted by a parasitoid fly, Ormia ochracea, which is attracted to singing males. Prior work indicates at least two of these flatwing phenotypes, from islands of Kauai and Oahu, have evolved independently under this shared selection pressure. This example of rapid, convergent evolution provides an opportunity to identify conditions which have promoted and resulted from rapid adaptation in wild populations evolving under extreme selection pressure. I investigate features which have contributed to the ability of these populations to rapidly, and repeatedly, adapt under strong selection against male song. The results indicate convergent sexual trait loss has been promoted by sex-biased development pathways maintained by sexually antagonistic selection; that pleiotropic, or associated, effects of adaptive mutation(s) in both sexes have played an important role in their spread; that adaptive male song-loss phenotypes have evolved repeatedly, above and beyond flatwing morphology; and that silent males nevertheless invest as much energy in practicing wing movement patterns associated with song and, despite reduced sexual dimorphism, are just as likely to be involved in aggressive intrasexual contests."This work was supported by the Natural Environment Research Council [grant numbers NE/I027800/1, NE/G014906/1, NE/L011255/1]. This work was supported by the University of St Andrews [School of Biology]." -- Funding (p. 2