Pan-cancer analyses reveal cancer-type-specific fungal ecologies and bacteriome interactions.

Cancer-microbe associations have been explored for centuries, but cancer-associated fungi have rarely been examined. Here, we comprehensively characterize the cancer mycobiome within 17,401 patient tissue, blood, and plasma samples across 35 cancer types in four independent cohorts. We report fungal DNA and cells at low abundances across many major human cancers, with differences in community compositions that differ among cancer types, even when accounting for technical background. Fungal histological staining of tissue microarrays supported intratumoral presence and frequent spatial association with cancer cells and macrophages. Comparing intratumoral fungal communities with matched bacteriomes and immunomes revealed co-occurring bi-domain ecologies, often with permissive, rather than competitive, microenvironments and distinct immune responses. Clinically focused assessments suggested prognostic and diagnostic capacities of the tissue and plasma mycobiomes, even in stage I cancers, and synergistic predictive performance with bacteriomes.

Meiotic Recombination: Genetics' Good Old Scalpel.

In the era of genome engineering, a new study returns to classical genetics to decipher genotype-phenotype relationships in unprecedented throughput and with unprecedented accuracy. Capitalizing on natural variation in yeast strains and frequent meiotic recombination, She and Jarosz (2018) dissect and map to nucleotide resolution, simple and complex determinants of diverse phenotypic traits.

Repertoires of tRNAs: The Couplers of Genomics and Proteomics.

The pool of transfer RNA (tRNA) molecules in cells allows the ribosome to decode genetic information. This repertoire of molecular decoders is positioned in the crossroad of the genome, the transcriptome, and the proteome. Omics and systems biology now allow scientists to explore the entire repertoire of tRNAs of many organisms, revealing basic exciting biology. The tRNA gene set of hundreds of species is now characterized, in addition to the tRNA genes of organelles and viruses. Genes encoding tRNAs for certain anticodon types appear in dozens of copies in a genome, while others are universally absent from any genome. Transcriptome measurement of tRNAs is challenging, but in recent years new technologies have allowed researchers to determine the dynamic expression patterns of tRNAs. These advances reveal that availability of ready-to-translate tRNA molecules is highly controlled by several transcriptional and posttranscriptional regulatory processes. This regulation shapes the proteome according to the cellular state. The tRNA pool profoundly impacts many aspects of cellular and organismal life, including protein expression level, translation accuracy, adequacy of folding, and even mRNA stability. As a result, the shape of the tRNA pool affects organismal health and may participate in causing conditions such as cancer and neurological conditions.

A relay race on the evolutionary adaptation spectrum.

Adaptation is the process in which organisms improve their fitness by changing their phenotype using genetic or non-genetic mechanisms. The adaptation toolbox consists of varied molecular and genetic means that we posit span an almost continuous "adaptation spectrum." Different adaptations are characterized by the time needed for organisms to attain them and by their duration. We suggest that organisms often adapt by progressing the adaptation spectrum, starting with rapidly attained physiological and epigenetic adaptations and culminating with slower long-lasting genetic ones. A tantalizing possibility is that earlier adaptations facilitate realization of later ones.

A dual program for translation regulation in cellular proliferation and differentiation.

A dichotomous choice for metazoan cells is between proliferation and differentiation. Measuring tRNA pools in various cell types, we found two distinct subsets, one that is induced in proliferating cells, and repressed otherwise, and another with the opposite signature. Correspondingly, we found that genes serving cell-autonomous functions and genes involved in multicellularity obey distinct codon usage. Proliferation-induced and differentiation-induced tRNAs often carry anticodons that correspond to the codons enriched among the cell-autonomous and the multicellularity genes, respectively. Because mRNAs of cell-autonomous genes are induced in proliferation and cancer in particular, the concomitant induction of their codon-enriched tRNAs suggests coordination between transcription and translation. Histone modifications indeed change similarly in the vicinity of cell-autonomous genes and their corresponding tRNAs, and in multicellularity genes and their tRNAs, suggesting the existence of transcriptional programs coordinating tRNA supply and demand. Hence, we describe the existence of two distinct translation programs that operate during proliferation and differentiation.

The evolutionary safety of mutagenic drugs should be assessed before drug approval.

Some drugs increase the mutation rate of their target pathogen, a potentially concerning mechanism as the pathogen might evolve faster toward an undesired phenotype. We suggest a four-step assessment of evolutionary safety for the approval of such treatments.

Systematic Detection of Amino Acid Substitutions in Proteomes Reveals Mechanistic Basis of Ribosome Errors and Selection for Translation Fidelity.

The translation machinery and the genes it decodes co-evolved to achieve production throughput and accuracy. Nonetheless, translation errors are frequent, and they affect physiology and protein evolution. Mapping translation errors in proteomes and understanding their causes is hindered by lack of a proteome-wide experimental methodology. We present the first methodology for systematic detection and quantification of errors in entire proteomes. Following proteome mass spectrometry, we identify, in E. coli and yeast, peptides whose mass indicates specific amino acid substitutions. Most substitutions result from codon-anticodon mispairing. Errors occur at sites that evolve rapidly and that minimally affect energetic stability, indicating selection for high translation fidelity. Ribosome density data show that errors occur at sites where ribosome velocity is higher, demonstrating a trade-off between speed and accuracy. Treating bacteria with an aminoglycoside antibiotic or deprivation of specific amino acids resulted in particular patterns of errors. These results reveal a mechanistic and evolutionary basis for translation fidelity.

Evolutionary safety of lethal mutagenesis driven by antiviral treatment.

Nucleoside analogs are a major class of antiviral drugs. Some act by increasing the viral mutation rate causing lethal mutagenesis of the virus. Their mutagenic capacity, however, may lead to an evolutionary safety concern. We define evolutionary safety as a probabilistic assurance that the treatment will not generate an increased number of mutants. We develop a mathematical framework to estimate the total mutant load produced with and without mutagenic treatment. We predict rates of appearance of such virus mutants as a function of the timing of treatment and the immune competence of patients, employing realistic assumptions about the vulnerability of the viral genome and its potential to generate viable mutants. We focus on the case study of Molnupiravir, which is an FDA-approved treatment against Coronavirus Disease-2019 (COVID-19). We estimate that Molnupiravir is narrowly evolutionarily safe, subject to the current estimate of parameters. Evolutionary safety can be improved by restricting treatment with this drug to individuals with a low immunological clearance rate and, in future, by designing treatments that lead to a greater increase in mutation rate. We report a simple mathematical rule to determine the fold increase in mutation rate required to obtain evolutionary safety that is also applicable to other pathogen-treatment combinations.

Aneuploidy Can Be an Evolutionary Diversion on the Path to Adaptation.

Aneuploidy is common in eukaryotes, often leading to decreased fitness. However, evidence from fungi and human tumur cells suggests that specific aneuploidies can be beneficial under stressful conditions and facilitate adaptation. In a previous evolutionary experiment with yeast, populations evolving under heat stress became aneuploid, only to later revert to euploidy after beneficial mutations accumulated. It was therefore suggested that aneuploidy is a "stepping stone" on the path to adaptation. Here, we test this hypothesis. We use Bayesian inference to fit an evolutionary model with both aneuploidy and mutation to the experimental results. We then predict the genotype frequency dynamics during the experiment, demonstrating that most of the evolved euploid population likely did not descend from aneuploid cells, but rather from the euploid wild-type population. Our model shows how the beneficial mutation supply-the product of population size and beneficial mutation rate-determines the evolutionary dynamics: with low supply, much of the evolved population descends from aneuploid cells; but with high supply, beneficial mutations are generated fast enough to outcompete aneuploidy due to its inherent fitness cost. Our results suggest that despite its potential fitness benefits under stress, aneuploidy can be an evolutionary "diversion" rather than a "stepping stone": it can delay, rather than facilitate, the adaptation of the population, and cells that become aneuploid may leave less descendants compared to cells that remain diploid.

An evolutionarily conserved mechanism for controlling the efficiency of protein translation.

Recent years have seen intensive progress in measuring protein translation. However, the contributions of coding sequences to the efficiency of the process remain unclear. Here, we identify a universally conserved profile of translation efficiency along mRNAs computed based on adaptation between coding sequences and the tRNA pool. In this profile, the first approximately 30-50 codons are, on average, translated with a low efficiency. Additionally, in eukaryotes, the last approximately 50 codons show the highest efficiency over the full coding sequence. The profile accurately predicts position-dependent ribosomal density along yeast genes. These data suggest that translation speed and, as a consequence, ribosomal density are encoded by coding sequences and the tRNA pool. We suggest that the slow "ramp" at the beginning of mRNAs serves as a late stage of translation initiation, forming an optimal and robust means to reduce ribosomal traffic jams, thus minimizing the cost of protein expression.

Gene Architectures that Minimize Cost of Gene Expression.

Gene expression burdens cells by consuming resources and energy. While numerous studies have investigated regulation of expression level, little is known about gene design elements that govern expression costs. Here, we ask how cells minimize production costs while maintaining a given protein expression level and whether there are gene architectures that optimize this process. We measured fitness of ∼14,000 E. coli strains, each expressing a reporter gene with a unique 5' architecture. By comparing cost-effective and ineffective architectures, we found that cost per protein molecule could be minimized by lowering transcription levels, regulating translation speeds, and utilizing amino acids that are cheap to synthesize and that are less hydrophobic. We then examined natural E. coli genes and found that highly expressed genes have evolved more forcefully to minimize costs associated with their expression. Our study thus elucidates gene design elements that improve the economy of protein expression in natural and heterologous systems.

Correction: A spatial vaccination strategy to reduce the risk of vaccine-resistant variants.

[This corrects the article DOI: 10.1371/journal.pcbi.1010391.].

Genetic redundancy: new tricks for old genes.

Many crucial components of signal transduction, developmental, and metabolic pathways have functionally redundant copies. Further, these redundancies show surprising evolutionary stability over prolonged time scales. We propose that redundancies are not just archeological leftovers of ancient gene duplications, but rather that synergy arising from feedback between redundant copies may serve as an information processing element that facilitates signal transduction and the control of gene expression.

A broad analysis of splicing regulation in yeast using a large library of synthetic introns.

RNA splicing is a key process in eukaryotic gene expression, in which an intron is spliced out of a pre-mRNA molecule to eventually produce a mature mRNA. Most intron-containing genes are constitutively spliced, hence efficient splicing of an intron is crucial for efficient regulation of gene expression. Here we use a large synthetic oligo library of ~20,000 variants to explore how different intronic sequence features affect splicing efficiency and mRNA expression levels in S. cerevisiae. Introns are defined by three functional sites, the 5' donor site, the branch site, and the 3' acceptor site. Using a combinatorial design of synthetic introns, we demonstrate how non-consensus splice site sequences in each of these sites affect splicing efficiency. We then show that S. cerevisiae splicing machinery tends to select alternative 3' splice sites downstream of the original site, and we suggest that this tendency created a selective pressure, leading to the avoidance of cryptic splice site motifs near introns' 3' ends. We further use natural intronic sequences from other yeast species, whose splicing machineries have diverged to various extents, to show how intron architectures in the various species have been adapted to the organism's splicing machinery. We suggest that the observed tendency for cryptic splicing is a result of a loss of a specific splicing factor, U2AF1. Lastly, we show that synthetic sequences containing two introns give rise to alternative RNA isoforms in S. cerevisiae, demonstrating that merely a synthetic fusion of two introns might be suffice to facilitate alternative splicing in yeast. Our study reveals novel mechanisms by which introns are shaped in evolution to allow cells to regulate their transcriptome. In addition, it provides a valuable resource to study the regulation of constitutive and alternative splicing in a model organism.

Dynamic changes in tRNA modifications and abundance during T cell activation.

The tRNA pool determines the efficiency, throughput, and accuracy of translation. Previous studies have identified dynamic changes in the tRNA (transfer RNA) supply and mRNA (messenger RNA) demand during cancerous proliferation. Yet dynamic changes may also occur during physiologically normal proliferation, and these are less well characterized. We examined the tRNA and mRNA pools of T cells during their vigorous proliferation and differentiation upon triggering their antigen receptor. We observed a global signature of switch in demand for codons at the early proliferation phase of the response, accompanied by corresponding changes in tRNA expression levels. In the later phase, upon differentiation, the response of the tRNA pool relaxed back to the basal level, potentially restraining excessive proliferation. Sequencing of tRNAs allowed us to evaluate their diverse base-modifications. We found that two types of tRNA modifications, wybutosine and ms2t6A, are reduced dramatically during T cell activation. These modifications occur in the anticodon loops of two tRNAs that decode "slippery codons," which are prone to ribosomal frameshifting. Attenuation of these frameshift-protective modifications is expected to increase the potential for proteome-wide frameshifting during T cell proliferation. Indeed, human cell lines deleted of a wybutosine writer showed increased ribosomal frameshifting, as detected with an HIV gag-pol frameshifting site reporter. These results may explain HIV's specific tropism toward proliferating T cells since it requires ribosomal frameshift exactly on the corresponding codon for infection. The changes in tRNA expression and modifications uncover a layer of translation regulation during T cell proliferation and expose a potential tradeoff between cellular growth and translation fidelity.

Identification and functional implications of pseudouridine RNA modification on small noncoding RNAs in the mammalian pathogen Trypanosoma brucei.

Trypanosoma brucei, the parasite that causes sleeping sickness, cycles between an insect and a mammalian host. However, the effect of RNA modifications such as pseudouridinylation on its ability to survive in these two different host environments is unclear. Here, two genome-wide approaches were applied for mapping pseudouridinylation sites (Ψs) on small nucleolar RNA (snoRNA), 7SL RNA, vault RNA, and tRNAs from T. brucei. We show using HydraPsiSeq and RiboMeth-seq that the Ψ on C/D snoRNA guiding 2'-O-methylation increased the efficiency of the guided modification on its target, rRNA. We found differential levels of Ψs on these noncoding RNAs in the two life stages (insect host and mammalian host) of the parasite. Furthermore, tRNA isoform abundance and Ψ modifications were characterized in these two life stages demonstrating stage-specific regulation. We conclude that the differential Ψ modifications identified here may contribute to modulating the function of noncoding RNAs involved in rRNA processing, rRNA modification, protein synthesis, and protein translocation during cycling of the parasite between its two hosts.

Fitness Landscape Analysis of a tRNA Gene Reveals that the Wild Type Allele is Sub-optimal, Yet Mutationally Robust.

Fitness landscape mapping and the prediction of evolutionary trajectories on these landscapes are major tasks in evolutionary biology research. Evolutionary dynamics is tightly linked to the landscape topography, but this relation is not straightforward. Here, we analyze a fitness landscape of a yeast tRNA gene, previously measured under four different conditions. We find that the wild type allele is sub-optimal, and 8-10% of its variants are fitter. We rule out the possibilities that the wild type is fittest on average on these four conditions or located on a local fitness maximum. Notwithstanding, we cannot exclude the possibility that the wild type might be fittest in some of the many conditions in the complex ecology that yeast lives at. Instead, we find that the wild type is mutationally robust ("flat"), while more fit variants are typically mutationally fragile. Similar observations of mutational robustness or flatness have been so far made in very few cases, predominantly in viral genomes.

A spatial vaccination strategy to reduce the risk of vaccine-resistant variants.

The COVID-19 pandemic demonstrated that the process of global vaccination against a novel virus can be a prolonged one. Social distancing measures, that are initially adopted to control the pandemic, are gradually relaxed as vaccination progresses and population immunity increases. The result is a prolonged period of high disease prevalence combined with a fitness advantage for vaccine-resistant variants, which together lead to a considerably increased probability for vaccine escape. A spatial vaccination strategy is proposed that has the potential to dramatically reduce this risk. Rather than dispersing the vaccination effort evenly throughout a country, distinct geographic regions of the country are sequentially vaccinated, quickly bringing each to effective herd immunity. Regions with high vaccination rates will then have low infection rates and vice versa. Since people primarily interact within their own region, spatial vaccination reduces the number of encounters between infected individuals (the source of mutations) and vaccinated individuals (who facilitate the spread of vaccine-resistant strains). Thus, spatial vaccination may help mitigate the global risk of vaccine-resistant variants.

Evolution of intron splicing towards optimized gene expression is based on various Cis- and Trans-molecular mechanisms.

Splicing expands, reshapes, and regulates the transcriptome of eukaryotic organisms. Despite its importance, key questions remain unanswered, including the following: Can splicing evolve when organisms adapt to new challenges? How does evolution optimize inefficiency of introns' splicing and of the splicing machinery? To explore these questions, we evolved yeast cells that were engineered to contain an inefficiently spliced intron inside a gene whose protein product was under selection for an increased expression level. We identified a combination of mutations in Cis (within the gene of interest) and in Trans (in mRNA-maturation machinery). Surprisingly, the mutations in Cis resided outside of known intronic functional sites and improved the intron's splicing efficiency potentially by easing tight mRNA structures. One of these mutations hampered a protein's domain that was not under selection, demonstrating the evolutionary flexibility of multi-domain proteins as one domain functionality was improved at the expense of the other domain. The Trans adaptations resided in two proteins, Npl3 and Gbp2, that bind pre-mRNAs and are central to their maturation. Interestingly, these mutations either increased or decreased the affinity of these proteins to mRNA, presumably allowing faster spliceosome recruitment or increased time before degradation of the pre-mRNAs, respectively. Altogether, our work reveals various mechanistic pathways toward optimizations of intron splicing to ultimately adapt gene expression patterns to novel demands.