Pervasive and persistent redundancy among duplicated genes in yeast

Evaluating sequence and structural similarity metrics for predicting shared paralog functions

Gene duplication is the primary source of new genes, resulting in most genes having identifiable paralogs. Over time, paralog pairs may diverge in some respects but many retain the ability to perform the same functional role. Protein sequence identity is often used as a proxy for functional similarity and can predict shared functions between paralogs as revealed by synthetic lethal experiments. However, the advent of alternative protein representations, including embeddings from protein language models (PLMs) and predicted structures from AlphaFold, raises the possibility that alternative similarity metrics could better capture functional similarity between paralogs. Here, using two species (budding yeast and human) and two different definitions of shared functionality (shared protein-protein interactions and synthetic lethality), we evaluated a variety of alternative similarity metrics. For some tasks, predicted structural similarity or PLM similarity outperform sequence identity, but more importantly these similarity metrics are not redundant with sequence identity, i.e. combining them with sequence identity leads to improved predictions of shared functionality. By adding contextual features, representing similarity to homologous proteins within and across species, we can significantly enhance our predictions of shared paralog functionality. Overall, our results suggest that alternative similarity metrics capture complementary aspects of functional similarity beyond sequence identity alone.

Cryptic genetic variation shapes the fate of gene duplicates in a protein interaction network

Paralogous genes are often functionally redundant for long periods of time. While their functions are preserved, paralogs accumulate cryptic changes in sequence and expression, which could modulate the impact of future mutations through epistasis. We examine the impact of mutations on redundant myosin proteins that have maintained the same binding preference despite having accumulated differences in expression levels and amino acid substitutions in the last 100 million years. By quantifying the impact of all single-amino acid substitutions in their SH3 domains on the physical interaction with their interaction partners, we show that the same mutations in the paralogous SH3s change binding in a paralog-specific and interaction partner-specific manner. This contingency is explained by the difference in promoter strength of the two paralogous myosin genes and epistatic interactions between the mutations introduced and cryptic divergent sites within the SH3s. One significant consequence of this contingency is that while some mutations would be sufficient to nonfunctionalize one paralog, they would have minimal impact on the other. Our results reveal how cryptic divergence, which accumulates while maintaining functional redundancy in cellular networks, could bias gene duplicates to specific fates.

Genome engineering on size reduction and complexity simplification: A review

BACKGROUND

Genome simplification is an important topic in the field of life sciences that has attracted attention from its conception to the present day. It can help uncover the essential components of the genome and, in turn, shed light on the underlying operating principles of complex biological systems. This has made it a central focus of both basic and applied research in the life sciences. With the recent advancements in related technologies and our increasing knowledge of the genome, now is an opportune time to delve into this topic.

AIM OF REVIEW

Our review investigates the progress of genome simplification from two perspectives: genome size reduction and complexity simplification. In addition, we provide insights into the future development trends of genome simplification.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Reducing genome size requires eliminating non-essential elements as much as possible. This process has been facilitated by advances in genome manipulation and synthesis techniques. However, we still need a better and clearer understanding of living systems to reduce genome complexity. As there is a lack of quantitative and clearly defined standards for this task, we have opted to approach the topic from various perspectives and present our findings accordingly.

Cryptic genetic variation shapes the fate of gene duplicates in a protein interaction network

Paralogous genes are often redundant for long periods of time before they diverge in function. While their functions are preserved, paralogous proteins can accumulate mutations that, through epistasis, could impact their fate in the future. By quantifying the impact of all single-amino acid substitutions on the binding of two myosin proteins to their interaction partners, we find that the future evolution of these proteins is highly contingent on their regulatory divergence and the mutations that have silently accumulated in their protein binding domains. Differences in the promoter strength of the two paralogs amplify the impact of mutations on binding in the lowly expressed one. While some mutations would be sufficient to non-functionalize one paralog, they would have minimal impact on the other. Our results reveal how functionally equivalent protein domains could be destined to specific fates by regulatory and cryptic coding sequence changes that currently have little to no functional impact.

LTP2 hypomorphs show genotype-by-environment interaction in early seedling traits in Arabidopsis thaliana

Isogenic individuals can display seemingly stochastic phenotypic differences, limiting the accuracy of genotype-to-phenotype predictions. The extent of this phenotypic variation depends in part on genetic background, raising questions about the genes involved in controlling stochastic phenotypic variation. Focusing on early seedling traits in Arabidopsis thaliana, we found that hypomorphs of the cuticle-related gene LIPID TRANSFER PROTEIN 2 (LTP2) greatly increased variation in seedling phenotypes, including hypocotyl length, gravitropism and cuticle permeability. Many ltp2 hypocotyls were significantly shorter than wild-type hypocotyls while others resembled the wild-type. Differences in epidermal properties and gene expression between ltp2 seedlings with long and short hypocotyls suggest a loss of cuticle integrity as the primary determinant of the observed phenotypic variation. We identified environmental conditions that reveal or mask the increased variation in ltp2 hypomorphs and found that increased expression of its closest paralog LTP1 is necessary for ltp2 phenotypes. Our results illustrate how decreased expression of a single gene can generate starkly increased phenotypic variation in isogenic individuals in response to an environmental challenge.

Paralog dispensability shapes homozygous deletion patterns in tumor genomes

Genomic instability is a hallmark of cancer, resulting in tumor genomes having large numbers of genetic aberrations, including homozygous deletions of protein coding genes. That tumor cells remain viable in the presence of such gene loss suggests high robustness to genetic perturbation. In model organisms and cancer cell lines, paralogs have been shown to contribute substantially to genetic robustness-they are generally more dispensable for growth than singletons. Here, by analyzing copy number profiles of > 10,000 tumors, we test the hypothesis that the increased dispensability of paralogs shapes tumor genome evolution. We find that genes with paralogs are more likely to be homozygously deleted and that this cannot be explained by other factors known to influence copy number variation. Furthermore, features that influence paralog dispensability in cancer cell lines correlate with paralog deletion frequency in tumors. Finally, paralogs that are broadly essential in cancer cell lines are less frequently deleted in tumors than non-essential paralogs. Overall, our results suggest that homozygous deletions of paralogs are more frequently observed in tumor genomes because paralogs are more dispensable.

Synthetic Lethality Screening with Recursive Feature Machines

Synthetic lethality refers to a genetic interaction where the simultaneous perturbation of gene pairs leads to cell death. Synthetically lethal gene pairs (SL pairs) provide a potential avenue for selectively targeting cancer cells based on genetic vulnerabilities. The rise of large-scale gene perturbation screens such as the Cancer Dependency Map (DepMap) offers the opportunity to identify SL pairs automatically using machine learning. We build on a recently developed class of feature learning kernel machines known as Recursive Feature Machines (RFMs) to develop a pipeline for identifying SL pairs based on CRISPR viability data from DepMap. In particular, we first train RFMs to predict viability scores for a given CRISPR gene knockout from cell line embeddings consisting of gene expression and mutation features. After training, RFMs use a statistical operator known as average gradient outer product to provide weights for each feature indicating the importance of each feature in predicting cellular viability. We subsequently apply correlation-based filters to re-weight RFM feature importances and identify those features that are most indicative of low cellular viability. Our resulting pipeline is computationally efficient, taking under 3 minutes for analyzing all 17, 453 knockouts from DepMap for candidate SL pairs. We show that our pipeline more accurately recovers experimentally verified SL pairs than prior approaches. Moreover, our pipeline finds new candidate SL pairs, thereby opening novel avenues for identifying genetic vulnerabilities in cancer.

Active compensation for changes in TDH3 expression mediated by direct regulators of TDH3 in Saccharomyces cerevisiae

Genetic networks are surprisingly robust to perturbations caused by new mutations. This robustness is conferred in part by compensation for loss of a gene's activity by genes with overlapping functions, such as paralogs. Compensation occurs passively when the normal activity of one paralog can compensate for the loss of the other, or actively when a change in one paralog's expression, localization, or activity is required to compensate for loss of the other. The mechanisms of active compensation remain poorly understood in most cases. Here we investigate active compensation for the loss or reduction in expression of the Saccharomyces cerevisiae gene TDH3 by its paralog TDH2. TDH2 is upregulated in a dose-dependent manner in response to reductions in TDH3 by a mechanism requiring the shared transcriptional regulators Gcr1p and Rap1p. TDH1, a second and more distantly related paralog of TDH3, has diverged in its regulation and is upregulated by another mechanism. Other glycolytic genes regulated by Rap1p and Gcr1p show changes in expression similar to TDH2, suggesting that the active compensation by TDH3 paralogs is part of a broader homeostatic response mediated by shared transcriptional regulators.

LTP2 hypomorphs show genotype-by-environment interaction in early seedling traits in Arabidopsis thaliana

Isogenic individuals can display seemingly stochastic phenotypic differences, limiting the accuracy of genotype-to-phenotype predictions. The extent of this phenotypic variation depends in part on genetic background, raising questions about the genes involved in controlling stochastic phenotypic variation. Focusing on early seedling traits in Arabidopsis thaliana, we found that hypomorphs of the cuticle-related gene LTP2 greatly increased variation in seedling phenotypes, including hypocotyl length, gravitropism and cuticle permeability. Many ltp2 hypocotyls were significantly shorter than wild-type hypocotyls while others resembled the wild type. Differences in epidermal properties and gene expression between ltp2 seedlings with long and short hypocotyls suggest a loss of cuticle integrity as the primary determinant of the observed phenotypic variation. We identified environmental conditions that reveal or mask the increased variation in ltp2 hypomorphs, and found that increased expression of its closest paralog LTP1 is necessary for ltp2 phenotypes. Our results illustrate how decreased expression of a single gene can generate starkly increased phenotypic variation in isogenic individuals in response to an environmental challenge.

Dynamic interplay between RPL3- and RPL3L-containing ribosomes modulates mitochondrial activity in the mammalian heart

The existence of naturally occurring ribosome heterogeneity is now a well-acknowledged phenomenon. However, whether this heterogeneity leads to functionally diverse 'specialized ribosomes' is still a controversial topic. Here, we explore the biological function of RPL3L (uL3L), a ribosomal protein (RP) paralogue of RPL3 (uL3) that is exclusively expressed in skeletal muscle and heart tissues, by generating a viable homozygous Rpl3l knockout mouse strain. We identify a rescue mechanism in which, upon RPL3L depletion, RPL3 becomes up-regulated, yielding RPL3-containing ribosomes instead of RPL3L-containing ribosomes that are typically found in cardiomyocytes. Using both ribosome profiling (Ribo-seq) and a novel orthogonal approach consisting of ribosome pulldown coupled to nanopore sequencing (Nano-TRAP), we find that RPL3L modulates neither translational efficiency nor ribosome affinity towards a specific subset of transcripts. In contrast, we show that depletion of RPL3L leads to increased ribosome-mitochondria interactions in cardiomyocytes, which is accompanied by a significant increase in ATP levels, potentially as a result of fine-tuning of mitochondrial activity. Our results demonstrate that the existence of tissue-specific RP paralogues does not necessarily lead to enhanced translation of specific transcripts or modulation of translational output. Instead, we reveal a complex cellular scenario in which RPL3L modulates the expression of RPL3, which in turn affects ribosomal subcellular localization and, ultimately, mitochondrial activity.

Evolutionary conservation of sequence motifs at sites of protein modification

Gene duplications are common in biology and are likely to be an important source of functional diversification and specialization. The yeast Saccharomyces cerevisiae underwent a whole-genome duplication event early in evolution, and a substantial number of duplicated genes have been retained. We identified more than 3500 instances where only one of two paralogous proteins undergoes posttranslational modification despite having retained the same amino acid residue in both. We also developed a web-based search algorithm (CoSMoS.c.) that scores conservation of amino acid sequences based on 1011 wild and domesticated yeast isolates and used it to compare differentially modified pairs of paralogous proteins. We found that the most common modifications-phosphorylation, ubiquitylation, and acylation but not N-glycosylation-occur in regions of high sequence conservation. Such conservation is evident even for ubiquitylation and succinylation, where there is no established 'consensus site' for modification. Differences in phosphorylation were not associated with predicted secondary structure or solvent accessibility but did mirror known differences in kinase-substrate interactions. Thus, differences in posttranslational modification likely result from differences in adjoining amino acids and their interactions with modifying enzymes. By integrating data from large-scale proteomics and genomics analysis, in a system with such substantial genetic diversity, we obtained a more comprehensive understanding of the functional basis for genetic redundancies that have persisted for 100 million years.

Targeting synthetic lethal paralogs in cancer

Synthetic lethal interactions, where mutation of one gene renders cells sensitive to inhibition of another gene, can be exploited for the development of targeted therapeutics in cancer. Pairs of duplicate genes (paralogs) often share common functionality and hence are a potentially rich source of synthetic lethal interactions. Because the majority of human genes have paralogs, exploiting such interactions could be a widely applicable approach for targeting gene loss in cancer. Moreover, existing small-molecule drugs may exploit synthetic lethal interactions by inhibiting multiple paralogs simultaneously. Consequently, the identification of synthetic lethal interactions between paralogs could be extremely informative for drug development. Here we review approaches to identify such interactions and discuss some of the challenges of exploiting them.

Imaginal disk growth factors are Drosophila chitinase-like proteins with roles in morphogenesis and CO2 response

Chitinase-like proteins (CLPs) are members of the family 18 glycosyl hydrolases, which include chitinases and the enzymatically inactive CLPs. A mutation in the enzyme's catalytic site, conserved in vertebrates and invertebrates, allowed CLPs to evolve independently with functions that do not require chitinase activity. CLPs normally function during inflammatory responses, wound healing, and host defense, but when they persist at excessive levels at sites of chronic inflammation and in tissue-remodeling disorders, they correlate positively with disease progression and poor prognosis. Little is known, however, about their physiological function. Drosophila melanogaster has 6 CLPs, termed Imaginal disk growth factors (Idgfs), encoded by Idgf1, Idgf2, Idgf3, Idgf4, Idgf5, and Idgf6. In this study, we developed tools to facilitate characterization of the physiological roles of the Idgfs by deleting each of the Idgf genes using the CRISPR/Cas9 system and assessing loss-of-function phenotypes. Using null lines, we showed that loss of function for all 6 Idgf proteins significantly lowers viability and fertility. We also showed that Idgfs play roles in epithelial morphogenesis, maintaining proper epithelial architecture and cell shape, regulating E-cadherin and cortical actin, and remarkably, protecting these tissues against CO2 exposure. Defining the normal molecular mechanisms of CLPs is a key to understanding how deviations tip the balance from a physiological to a pathological state.

Active compensation for changes in TDH3 expression mediated by direct regulators of TDH3 in Saccharomyces cerevisiae

Genetic networks are surprisingly robust to perturbations caused by new mutations. This robustness is conferred in part by compensation for loss of a gene's activity by genes with overlapping functions, such as paralogs. Compensation occurs passively when the normal activity of one paralog can compensate for the loss of the other, or actively when a change in one paralog's expression, localization, or activity is required to compensate for loss of the other. The mechanisms of active compensation remain poorly understood in most cases. Here we investigate active compensation for the loss or reduction in expression of the Saccharomyces cerevisiae gene TDH3 by its paralogs TDH1 and TDH2. TDH1 and TDH2 are upregulated in a dose-dependent manner in response to reductions in TDH3 by a mechanism requiring the shared transcriptional regulators Gcr1p and Rap1p. Other glycolytic genes regulated by Rap1p and Gcr1p show changes in expression similar to TDH2, suggesting that the active compensation by TDH3 paralogs is part of a broader homeostatic response mediated by shared transcriptional regulators.

Variable paralog expression underlies phenotype variation

Human faces are variable; we look different from one another. Craniofacial disorders further increase facial variation. To understand craniofacial variation and how it can be buffered, we analyzed the zebrafish mef2ca mutant. When this transcription factor encoding gene is mutated, zebrafish develop dramatically variable craniofacial phenotypes. Years of selective breeding for low and high penetrance of mutant phenotypes produced strains that are either resilient or sensitive to the mef2ca mutation. Here, we compared gene expression between these strains, which revealed that selective breeding enriched for high and low mef2ca paralog expression in the low- and high-penetrance strains, respectively. We found that mef2ca paralog expression is variable in unselected wild-type zebrafish, motivating the hypothesis that heritable variation in paralog expression underlies mutant phenotype severity and variation. In support, mutagenizing the mef2ca paralogs, mef2aa, mef2b, mef2cb, and mef2d demonstrated modular buffering by paralogs. Specifically, some paralogs buffer severity while others buffer variability. We present a novel, mechanistic model for phenotypic variation where variable, vestigial paralog expression buffers development. These studies are a major step forward in understanding the mechanisms of facial variation, including how some genetically resilient individuals can overcome a deleterious mutation.

Functional Divergence in a Multi-gene Family Is a Key Evolutionary Innovation for Anaerobic Growth in Saccharomyces cerevisiae

The amplification and diversification of genes into large multi-gene families often mark key evolutionary innovations, but this process often creates genetic redundancy that hinders functional investigations. When the model budding yeast Saccharomyces cerevisiae transitions to anaerobic growth conditions, the cell massively induces the expression of seven serine/threonine-rich anaerobically-induced cell wall mannoproteins (anCWMPs): TIP1, TIR1, TIR2, TIR3, TIR4, DAN1, and DAN4. Here, we show that these genes likely derive evolutionarily from a single ancestral anCWMP locus, which was duplicated and translocated to new genomic contexts several times both prior to and following the budding yeast whole genome duplication (WGD) event. Based on synteny and their phylogeny, we separate the anCWMPs into four gene subfamilies. To resolve prior inconclusive genetic investigations of these genes, we constructed a set of combinatorial deletion mutants to determine their contributions toward anaerobic growth in S. cerevisiae. We found that two genes, TIR1 and TIR3, were together necessary and sufficient for the anCWMP contribution to anaerobic growth. Overexpressing either gene alone was insufficient for anaerobic growth, implying that they encode non-overlapping functional roles in the cell during anaerobic growth. We infer from the phylogeny of the anCWMP genes that these two important genes derive from an ancient duplication that predates the WGD event, whereas the TIR1 subfamily experienced gene family amplification after the WGD event. Taken together, the genetic and molecular evidence suggests that one key anCWMP gene duplication event, several auxiliary gene duplication events, and functional divergence underpin the evolution of anaerobic growth in budding yeasts.

Individual Genetic Heterogeneity

Genetic variation has been widely covered in literature, however, not from the perspective of an individual in any species. Here, a synthesis of genetic concepts and variations relevant for individual genetic constitution is provided. All the different levels of genetic information and variation are covered, ranging from whether an organism is unmixed or hybrid, has variations in genome, chromosomes, and more locally in DNA regions, to epigenetic variants or alterations in selfish genetic elements. Genetic constitution and heterogeneity of microbiota are highly relevant for health and wellbeing of an individual. Mutation rates vary widely for variation types, e.g., due to the sequence context. Genetic information guides numerous aspects in organisms. Types of inheritance, whether Mendelian or non-Mendelian, zygosity, sexual reproduction, and sex determination are covered. Functions of DNA and functional effects of variations are introduced, along with mechanism that reduce and modulate functional effects, including TARAR countermeasures and intraindividual genetic conflict. TARAR countermeasures for tolerance, avoidance, repair, attenuation, and resistance are essential for life, integrity of genetic information, and gene expression. The genetic composition, effects of variations, and their expression are considered also in diseases and personalized medicine. The text synthesizes knowledge and insight on individual genetic heterogeneity and organizes and systematizes the central concepts.

A Simple Evolutionary Model of Genetic Robustness After Gene Duplication

When a dispensable gene is duplicated (referred to the ancestral dispensability denoted by O⁺), genetic buffering and duplicate compensation together maintain the duplicate redundancy, whereas duplicate compensation is the only mechanism when an essential gene is duplicated (referred to the ancestral essentiality denoted by O^-). To investigate these evolutionary scenarios of genetic robustness, I formulated a simple mixture model for analyzing duplicate pairs with one of the following states: double dispensable (DD), semi-dispensable (one dispensable one essential, DE), or double essential (EE). This model was applied to the yeast duplicate pairs from a whole-genome duplication (WGD) occurred about 100 million years ago (mya), and the mouse duplicate pairs from a WGD occurred about more than 500 mya. Both case studies revealed that the proportion of essentiality for those duplicates with ancestral essentiality [P_E(O^-)] was much higher than that for those with ancestral dispensability [P_E(O⁺)]. While it was negligible in the yeast duplicate pairs, P_E(O⁺) (about 20%) was shown statistically significant in the mouse duplicate pairs. These findings, together, support the hypothesis that both sub-functionalization and neo-functionalization may play some roles after gene duplication, though the former may be much faster than the later.

MULTI-SCULPT: Multiplex Integration via Selective, CRISPR-Mediated, Ultralong Pathway Transformation in Yeast for Plant Natural Product Synthesis

Yeast has been a versatile model host for complex and valuable natural product biosynthesis via the reconstruction of heterologous biosynthetic pathways. Recent advances in natural product pathway elucidation have uncovered many large and complicated plant pathways that contain 10-30 genes for the biosynthesis of structurally complex, valuable natural products. However, the ability to reconstruct ultralong pathways efficiently in yeast does not match the increasing demand for valuable plant natural product biomanufacturing. Here, we developed a one-pot, multigene pathway integration method in yeast, named MULTI-SCULPT for multiplex integration via selective, CRISPR-mediated, ultralong pathway transformation. Leveraging multilocus genomic disruption via CRISPR/Cas9, newly developed native and synthetic genetic parts, and fine-tuned gene integration and characterization methods, we managed to integrate 21 DNA inserts that contain a 12-gene plant isoflavone biosynthetic pathway into yeast with a 90-100% success rate in 12 days. This method enables fast and efficient ultralong biosynthetic pathway integration and can allow for the fast iterative integration of even longer pathways in the future. Ultimately, this method will accelerate combinatorial optimization of elucidated plant natural product pathways and accelerate putative pathway characterization heterologously.

Genome-Wide Analysis of the MADS-Box Gene Family in Holoparasitic Plants (Balanophora subcupularis and Balanophora fungosa var. globosa)

MADS-box is an important transcription factor family that is involved in the regulation of various stages of plant growth and development, especially flowering regulation and flower development. Being a holoparasitic plant, the body structure of Balanophoraceae has changed dramatically over time, and its vegetative and reproductive organs have been extensively modified, with rudimentary flower organs. Meanwhile, extraordinary gene losses have been identified in holoparasitic plants compared with autotrophs. Our study reveals that the MADS-box gene family contracted sharply in Balanophora subcupularis and Balanophora fungosa var. globosa, and some subfamilies were lost, exhibiting reduced redundancy in both. The genes that functioned in the transition from the vegetative to floral production stages suffered a significant loss, but the ABCE model genes remained intact. We further investigated genes related to flowering regulation in B. subcupularis and B. fungosa var. globosa, vernalization and autonomous ways of regulating flowering time remained comparatively integrated, while genes in photoperiod and circadian clock pathways were almost lost. Convergent gene loss in flowering regulation occurred in Balanophora and another holoparasitic plant Sapria himalayana (Rafflesiaceae). The genome-wide analysis of the MADS-box gene family in Balanophora species provides valuable information for understanding the classification, gene loss pattern, and flowering regulation mechanism of MADS-box gene family in parasitic plants.

Beyond Mendelian Inheritance: Genetic Buffering and Phenotype Variability

Understanding the way genes work amongst individuals and across generations to shape form and function is a common theme for many genetic studies. The recent advances in genetics, genome engineering and DNA sequencing reinforced the notion that genes are not the only players that determine a phenotype. Due to physiological or pathological fluctuations in gene expression, even genetically identical cells can behave and manifest different phenotypes under the same conditions. Here, we discuss mechanisms that can influence or even disrupt the axis between genotype and phenotype; the role of modifier genes, the general concept of genetic redundancy, genetic compensation, the recently described transcriptional adaptation, environmental stressors, and phenotypic plasticity. We furthermore highlight the usage of induced pluripotent stem cells (iPSCs), the generation of isogenic lines through genome engineering, and sequencing technologies can help extract new genetic and epigenetic mechanisms from what is hitherto considered 'noise'.

Near Chromosome-Level Genome Assembly and Annotation of Rhodotorula babjevae Strains Reveals High Intraspecific Divergence

The genus Rhodotorula includes basidiomycetous oleaginous yeast species. Rhodotorula babjevae can produce compounds of biotechnological interest such as lipids, carotenoids, and biosurfactants from low value substrates such as lignocellulose hydrolysate. High-quality genome assemblies are needed to develop genetic tools and to understand fungal evolution and genetics. Here, we combined short- and long-read sequencing to resolve the genomes of two R. babjevae strains, CBS 7808 (type strain) and DBVPG 8058, at chromosomal level. Both genomes are 21 Mbp in size and have a GC content of 68.2%. Allele frequency analysis indicates that both strains are tetraploid. The genomes consist of a maximum of 21 chromosomes with a size of 0.4 to 2.4 Mbp. In both assemblies, the mitochondrial genome was recovered in a single contig, that shared 97% pairwise identity. Pairwise identity between most chromosomes ranges from 82 to 87%. We also found indications for strain-specific extrachromosomal endogenous DNA. A total of 7591 and 7481 protein-coding genes were annotated in CBS 7808 and DBVPG 8058, respectively. CBS 7808 accumulated a higher number of tandem duplications than DBVPG 8058. We identified large translocation events between putative chromosomes. Genome divergence values between the two strains indicate that they may belong to different species.

Multi-Omics Analysis of Multiple Glucose-Sensing Receptor Systems in Yeast

The yeast Saccharomyces cerevisiae has long been used to produce alcohol from glucose and other sugars. While much is known about glucose metabolism, relatively little is known about the receptors and signaling pathways that indicate glucose availability. Here, we compare the two glucose receptor systems in S. cerevisiae. The first is a heterodimer of transporter-like proteins (transceptors), while the second is a seven-transmembrane receptor coupled to a large G protein (Gpa2) that acts in coordination with two small G proteins (Ras1 and Ras2). Through comprehensive measurements of glucose-dependent transcription and metabolism, we demonstrate that the two receptor systems have distinct roles in glucose signaling: the G-protein-coupled receptor directs carbohydrate and energy metabolism, while the transceptors regulate ancillary processes such as ribosome, amino acids, cofactor and vitamin metabolism. The large G-protein transmits the signal from its cognate receptor, while the small G-protein Ras2 (but not Ras1) integrates responses from both receptor pathways. Collectively, our analysis reveals the molecular basis for glucose detection and the earliest events of glucose-dependent signal transduction in yeast.

Yeast Crf1p: An activator in need is an activator indeed

Ribosome biogenesis is an energetically costly process, and tight regulation is required for stoichiometric balance between components. This requires coordination of RNA polymerases I, II, and III. Lack of nutrients or the presence of stress leads to downregulation of ribosome biogenesis, a process for which mechanistic target of rapamycin complex I (mTORC1) is key. mTORC1 activity is communicated by means of specific transcription factors, and in yeast, which is a primary model system in which transcriptional coordination has been delineated, transcription factors involved in regulation of ribosomal protein genes include Fhl1p and its cofactors, Ifh1p and Crf1p. Ifh1p is an activator, whereas Crf1p has been implicated in maintaining the repressed state upon mTORC1 inhibition. Computational analyses of evolutionary relationships have indicated that Ifh1p and Crf1p descend from a common ancestor. Here, we discuss recent evidence, which suggests that Crf1p also functions as an activator. We propose a model that consolidates available experimental evidence, which posits that Crf1p functions as an alternate activator to prevent the stronger activator Ifh1p from re-binding gene promoters upon mTORC1 inhibition. The correlation between retention of Crf1p in related yeast strains and duplication of ribosomal protein genes suggests that this backup activation may be important to ensure gene expression when Ifh1p is limiting. With ribosome biogenesis as a hallmark of cell growth, failure to control assembly of ribosomal components leads to several human pathologies. A comprehensive understanding of mechanisms underlying this process is therefore of the essence.

Retention of duplicated genes in evolution

Gene duplication is a prevalent phenomenon across the tree of life. The processes that lead to the retention of duplicated genes are not well understood. Functional genomics approaches in model organisms, such as yeast, provide useful tools to test the mechanisms underlying retention with functional redundancy and divergence of duplicated genes, including fates associated with neofunctionalization, subfunctionalization, back-up compensation, and dosage amplification. Duplicated genes may also be retained as a consequence of structural and functional entanglement. Advances in human gene editing have enabled the interrogation of duplicated genes in the human genome, providing new tools to evaluate the relative contributions of each of these factors to duplicate gene retention and the evolution of genome structure.

Impact of Genome Reduction in Microsporidia

Microsporidia represent an evolutionary outlier in the tree of life and occupy the extreme edge of the eukaryotic domain with some of their biological features. Many of these unicellular fungi-like organisms have reduced their genomic content to potentially the lowest limit. With some of the most compacted eukaryotic genomes, microsporidia are excellent model organisms to study reductive evolution and its functional consequences. While the growing number of sequenced microsporidian genomes have elucidated genome composition and organization, a recent increase in complementary post-genomic studies has started to shed light on the impacts of genome reduction in these unique pathogens. This chapter will discuss the biological framework enabling genome minimization and will use one of the most ancient and essential macromolecular complexes, the ribosome, to illustrate the effects of extreme genome reduction on a structural, molecular, and cellular level. We outline how reductive evolution in microsporidia has shaped DNA organization, the composition and function of the ribosome, and the complexity of the ribosome biogenesis process. Studying compacted mechanisms, processes, or macromolecular machines in microsporidia illuminates their unique lifestyle and provides valuable insights for comparative eukaryotic structural biology.

Ancient rapid functional differentiation and fixation of the duplicated members in rice Dof genes after whole genome duplication

Whole genome duplication (WGD) in plants is typically followed by genomic downsizing, where large portions of the new genome are lost. Whether this downsizing is accompanied by increased or decreased evolutionary rates of the remaining genes is poorly known, not least because homeolog pairings are often obscured by chromosomal rearrangement. Here, we use the newly published genome from a sedge, namely Kobresia littledalei, and CRISPR/Cas-9 editing to investigate how the Rho WGD event 70 million years ago (MYA) affected transcription factor evolutionary rates, fates, and function in rice (Oryza sativa) and sorghum (Sorghum bicolor). We focus on the 30-member DNA-binding with one zinc finger (Dof) transcription factor family in both crops due to their agronomic importance. Using the known speciation dates of rice from Kobresia (97 MYA) and sorghum (50 MYA), we find that rates of amino acid substitution in the critical Dof domain region were over twofold higher during the 20-million-year period following the WGD than before or afterward. Through comparison of synteny blocks, we report that at least 11% of Dof genes were purged from 70 to 50 MYA, while only 6% have been lost in the most recent 50-million-year interval. CRISPR/Cas9 editing revealed widespread fitness-related defects in flowering and lack of redundancy of paired members, as well as significant differences in expression between gene pairs. Together these findings demonstrate the strength of Dof genes as a model for deep evolutionary study and offer one of the most detailed portraits yet of the Rho WGD impact on a gene lineage.

Deletion of ribosomal paralogs Rpl39 and Rpl39l compromises cell proliferation via protein synthesis and mitochondrial activity

Accumulating evidences suggest that the composition and functional roles of ribosomes are heterogeneous in cells, partly due to the temporal-spatial expression of paralogous ribosomal proteins (RPs), of which functional relationships remain largely unexplored. In mouse, the X chromosome-linked RPL39 and its male germline specific paralog RPL39L are thought to express mutually exclusively due to the meiotic sex chromosome inactivation, hinders the understanding of their functional relationships. In the present study, we investigated the expression and functional relations of Rpl39 and Rpl39l in a proliferative mouse cell line, in which both genes are expressed simultaneously, with the expression level of Rpl39 higher than that of Rpl39l. Disruption of Rpl39 via CRISPR/Cas9 method caused decreased cell proliferation, nascent protein synthesis and altered mitochondrial functions, whereas double mutations of Rpl39 and Rpl39l augmented these phenotypes, suggesting that both RPs contribute to the cellular physiology. Consistently, overexpression of Rpl39, Rpl39l or an alanine mutant of RPL39, rescued cell proliferation similarly in Rpl39-/-::Rpl39l-/- dual gene null cells. Deletion of Rpl39l induced compensatory expression of Rpl39, rendering the deleterious effects of Rpl39l mutation. Supporting this, Rpl39l mutation was more detrimental to cells under a low serum condition, under which the compensatory expression of Rpl39 was inhibited. Moreover, the low serum condition induced expression of both genes, suggesting that they possess stress responsive roles. Taken together, these data indicate that both RPL39 and RPL39L influence cell proliferation via protein synthesis and mitochondrial functions, suggesting a link between protein translation and cellular metabolism through these ribosomal protein paralogs.

Comprehensive prediction of robust synthetic lethality between paralog pairs in cancer cell lines

Pairs of paralogs may share common functionality and, hence, display synthetic lethal interactions. As the majority of human genes have an identifiable paralog, exploiting synthetic lethality between paralogs may be a broadly applicable approach for targeting gene loss in cancer. However, only a biased subset of human paralog pairs has been tested for synthetic lethality to date. Here, by analyzing genome-wide CRISPR screens and molecular profiles of over 700 cancer cell lines, we identify features predictive of synthetic lethality between paralogs, including shared protein-protein interactions and evolutionary conservation. We develop a machine-learning classifier based on these features to predict which paralog pairs are most likely to be synthetic lethal and to explain why. We show that our classifier accurately predicts the results of combinatorial CRISPR screens in cancer cell lines and furthermore can distinguish pairs that are synthetic lethal in multiple cell lines from those that are cell-line specific. A record of this paper's transparent peer review process is included in the supplemental information.

Discovery of synthetic lethal and tumor suppressor paralog pairs in the human genome

CRISPR screens have accelerated the discovery of important cancer vulnerabilities. However, single-gene knockout phenotypes can be masked by redundancy among related genes. Paralogs constitute two-thirds of the human protein-coding genome, so existing methods are likely inadequate for assaying a large portion of gene function. Here, we develop paired guide RNAs for paralog genetic interaction mapping (pgPEN), a pooled CRISPR-Cas9 single- and double-knockout approach targeting more than 2,000 human paralogs. We apply pgPEN to two cell types and discover that 12% of human paralogs exhibit synthetic lethality in at least one context. We recover known synthetic lethal paralogs MEK1/MEK2, important drug targets CDK4/CDK6, and other synthetic lethal pairs including CCNL1/CCNL2. Additionally, we identify ten tumor suppressor paralog pairs whose compound loss promotes cell proliferation. These findings nominate drug targets and suggest that paralog genetic interactions could shape the landscape of positive and negative selection in cancer.

Trade-offs Among Resilience, Robustness, Stability, and Performance and How We Might Study Them

Biological systems are likely to be constrained by trade-offs among robustness, resilience, and performance. A better understanding of these trade-offs is important for basic biology, as well as applications where biological systems can be designed for different goals. We focus on redundancy and plasticity as mechanisms governing some types of trade-offs, but mention others as well. Whether trade-offs are due to resource constraints or "design" constraints (i.e., structure of nodes and links within a network) will also affect the types of trade-offs that are important. Identifying common themes across scales of biological organization will require that researchers use similar approaches to quantifying robustness, resilience, and performance, using units that can be compared across systems.

An insight on the impact of teleost whole genome duplication on the regulation of the molecular networks controlling skeletal muscle growth

Fish muscle growth is a complex process regulated by multiple pathways, resulting on the net accumulation of proteins and the activation of myogenic progenitor cells. Around 350–320 million years ago, teleost fish went through a specific whole genome duplication (WGD) that expanded the existent gene repertoire. Duplicated genes can be retained by different molecular mechanisms such as subfunctionalization, neofunctionalization or redundancy, each one with different functional implications. While the great majority of ohnolog genes have been identified in the teleost genomes, the effect of gene duplication in the fish physiology is still not well characterized. In the present study we studied the effect of WGD on the transcription of the duplicated components controlling muscle growth. We compared the expression of lineage-specific ohnologs related to myogenesis and protein balance in the fast-skeletal muscle of pacus (Piaractus mesopotamicus—Ostariophysi) and Nile tilapias (Oreochromis niloticus—Acanthopterygii) fasted for 4 days and refed for 3 days. We studied the expression of 20 ohnologs and found that in the great majority of cases, duplicated genes had similar expression profiles in response to fasting and refeeding, indicating that their functions during growth have been conserved during the period after the WGD. Our results suggest that redundancy might play a more important role in the retention of ohnologs of regulatory pathways than initially thought. Also, comparison to non-duplicated orthologs showed that it might not be uncommon for the duplicated genes to gain or loss new regulatory elements simultaneously. Overall, several of duplicated ohnologs have similar transcription profiles in response to pro-growth signals suggesting that evolution tends to conserve ohnolog regulation during muscle development and that in the majority of ohnologs related to muscle growth their functions might be very similar.

Yeast cell fate control by temporal redundancy modulation of transcription factor paralogs

Recent single-cell studies have revealed that yeast stress response involves transcription factors that are activated in pulses. However, it remains unclear whether and how these dynamic transcription factors temporally interact to regulate stress survival. Here we show that budding yeast cells can exploit the temporal relationship between paralogous general stress regulators, Msn2 and Msn4, during stress response. We find that individual pulses of Msn2 and Msn4 are largely redundant, and cells can enhance the expression of their shared targets by increasing their temporal divergence. Thus, functional redundancy between these two paralogs is modulated in a dynamic manner to confer fitness advantages for yeast cells, which might feed back to promote the preservation of their redundancy. This evolutionary implication is supported by evidence from Msn2/Msn4 orthologs and analyses of other transcription factor paralogs. Together, we show a cell fate control mechanism through temporal redundancy modulation in yeast, which may represent an evolutionarily important strategy for maintaining functional redundancy between gene duplicates.

Extensive non-redundancy in a recently duplicated developmental gene family

BACKGROUND

It has been proposed that recently duplicated genes are more likely to be redundant with one another compared to ancient paralogues. The evolutionary logic underpinning this idea is simple, as the assumption is that recently derived paralogous genes are more similar in sequence compared to members of ancient gene families. We set out to test this idea by using molecular phylogenetics and exploiting the genetic tractability of the model nematode, Caenorhabditis elegans, in studying the nematode-specific family of Hedgehog-related genes, the Warthogs. Hedgehog is one of a handful of signal transduction pathways that underpins the development of bilaterian animals. While having lost a bona fide Hedgehog gene, most nematodes have evolved an expanded repertoire of Hedgehog-related genes, ten of which reside within the Warthog family.

RESULTS

We have characterised their evolutionary origin and their roles in C. elegans and found that these genes have adopted new functions in aspects of post-embryonic development, including left-right asymmetry and cell fate determination, akin to the functions of their vertebrate counterparts. Analysis of various double and triple mutants of the Warthog family reveals that more recently derived paralogues are not redundant with one another, while a pair of divergent Warthogs do display redundancy with respect to their function in cuticle biosynthesis.

CONCLUSIONS

We have shown that newer members of taxon-restricted gene families are not always functionally redundant despite their recent inception, whereas much older paralogues can be, which is considered paradoxical according to the current framework in gene evolution.

Knockout of RSN1, TVP18 or CSC1-2 causes perturbation of Golgi cisternae in Pichia pastoris

The structural organization of the Golgi stacks in mammalian cells is intrinsically linked to function, including glycosylation, but the role of morphology is less clear in lower eukaryotes. Here we investigated the link between the structural organization of the Golgi and secretory pathway function using Pichia pastoris as a model system. To unstack the Golgi cisternae, we disrupted 18 genes encoding proteins in the secretory pathway without loss of viability. Using biosensors, confocal microscopy and transmission electron microscopy we identified three strains with irreversible perturbations in the stacking of the Golgi cisternae, all of which had disruption in genes that encode proteins with annotated function as or homology to calcium/calcium permeable ion channels. Despite this, no variation in the secretory pathway for ER size, whole cell glycomics or recombinant protein glycans was observed. Our investigations showed the robust nature of the secretory pathway in P. pastoris and suggest that Ca²⁺ concentration, homeostasis or signalling may play a significant role for Golgi stacking in this organism and should be investigated in other organisms.

Genomic signatures of globally enhanced gene duplicate accumulation in the megadiverse higher Diptera fueling intralocus sexual conflict resolution

Gene duplication is an important source of evolutionary innovation. To explore the relative impact of gene duplication during the diversification of major insect model system lineages, we performed a comparative analysis of lineage-specific gene duplications in the fruit fly Drosophila melanogaster (Diptera: Brachycera), the mosquito Anopheles gambiae (Diptera: Culicomorpha), the red flour beetle Tribolium castaneum (Coleoptera), and the honeybee Apis mellifera (Hymenoptera). Focusing on close to 6,000 insect core gene families containing maximally six paralogs, we detected a conspicuously higher number of lineage-specific duplications in Drosophila (689) compared to Anopheles (315), Tribolium (386), and Apis (223). Based on analyses of sequence divergence, phylogenetic distribution, and gene ontology information, we present evidence that an increased background rate of gene duplicate accumulation played an exceptional role during the diversification of the higher Diptera (Brachycera), in part by providing enriched opportunities for intralocus sexual conflict resolution, which may have boosted speciation rates during the early radiation of the megadiverse brachyceran subclade Schizophora.

Regulation of ribosomal protein genes: An ordered anarchy

Ribosomal protein genes are among the most highly expressed genes in most cell types. Their products are generally essential for ribosome synthesis, which is the cornerstone for cell growth and proliferation. Many cellular resources are dedicated to producing ribosomal proteins and thus this process needs to be regulated in ways that carefully balance the supply of nascent ribosomal proteins with the demand for new ribosomes. Ribosomal protein genes have classically been viewed as a uniform interconnected regulon regulated in eukaryotic cells by target of rapamycin and protein kinase A pathway in response to changes in growth conditions and/or cellular status. However, recent literature depicts a more complex picture in which the amount of ribosomal proteins produced varies between genes in response to two overlapping regulatory circuits. The first includes the classical general ribosome‐producing program and the second is a gene‐specific feature responsible for fine‐tuning the amount of ribosomal proteins produced from each individual ribosomal gene. Unlike the general pathway that is mainly controlled at the level of transcription and translation, this specific regulation of ribosomal protein genes is largely achieved through changes in pre‐mRNA splicing efficiency and mRNA stability. By combining general and specific regulation, the cell can coordinate ribosome production, while allowing functional specialization and diversity. Here we review the many ways ribosomal protein genes are regulated, with special focus on the emerging role of posttranscriptional regulatory events in fine‐tuning the expression of ribosomal protein genes and its role in controlling the potential variation in ribosome functions.

This article is categorized under:

Translation > Ribosome Biogenesis

Translation > Ribosome Structure/Function

Translation > Translation Regulation

Arabidopsis paralogous genes RPL23aA and RPL23aB encode functionally equivalent proteins

BACKGROUND

In plants, each ribosomal protein (RP) is encoded by a small gene family but it is largely unknown whether the family members are functionally diversified. There are two RPL23a paralogous genes (RPL23aA and RPL23aB) encoding cytoplasmic ribosomal proteins in Arabidopsis thaliana. Knock-down of RPL23aA using RNAi impeded growth and led to morphological abnormalities, whereas knock-out of RPL23aB had no observable phenotype, thus these two RPL23a paralogous proteins have been used as examples of ribosomal protein paralogues with functional divergence in many published papers.

RESULTS

In this study, we characterized T-DNA insertion mutants of RPL23aA and RPL23aB. A rare non-allelic non-complementation phenomenon was found in the F1 progeny of the rpl23aa X rpl23ab cross, which revealed a dosage effect of these two genes. Both RPL23aA and RPL23aB were found to be expressed almost in all examined tissues as revealed by GUS reporter analysis. Expression of RPL23aB driven by the RPL23aA promoter can rescue the phenotype of rpl23aa, indicating these two proteins are actually equivalent in function. Interestingly, based on the publicly available RNA-seq data, we found that these two RPL23a paralogues were expressed in a concerted manner and the expression level of RPL23aA was much higher than that of RPL23aB at different developmental stages and in different tissues.

CONCLUSIONS

Our findings suggest that the two RPL23a paralogous proteins are functionally equivalent but the two genes are not. RPL23aA plays a predominant role due to its higher expression levels. RPL23aB plays a lesser role due to its lower expression. The presence of paralogous genes for the RPL23a protein in plants might be necessary to maintain its adequate dosage.

An Overview of Duplicated Gene Detection Methods: Why the Duplication Mechanism Has to Be Accounted for in Their Choice

Gene duplication is an important evolutionary mechanism allowing to provide new genetic material and thus opportunities to acquire new gene functions for an organism, with major implications such as speciation events. Various processes are known to allow a gene to be duplicated and different models explain how duplicated genes can be maintained in genomes. Due to their particular importance, the identification of duplicated genes is essential when studying genome evolution but it can still be a challenge due to the various fates duplicated genes can encounter. In this review, we first describe the evolutionary processes allowing the formation of duplicated genes but also describe the various bioinformatic approaches that can be used to identify them in genome sequences. Indeed, these bioinformatic approaches differ according to the underlying duplication mechanism. Hence, understanding the specificity of the duplicated genes of interest is a great asset for tool selection and should be taken into account when exploring a biological question.

Ribosomal Protein uL11 as a Regulator of Metabolic Circuits Related to Aging and Cell Cycle

Aging is a biological phenomenon common to all living organisms. It is thought that the rate of aging is influenced by diverse factors, in many cases related to the control of energy metabolism, i.e., the so-called pro-longevity effects of starvation. Translation, regarded as the main energy consumption process, lies at the center of interest, as it has a significant impact on the longevity phenomenon. It has been shown that perturbations in the translational apparatus may lead to a lower rate of aging. Therefore, the main aim of this study was to investigate aging in relation to the protein biosynthesis circuit, taking into account the uL11 ribosomal protein as a vital ribosomal element. To this end, we used set of yeast mutants with deleted single uL11A or uL11B genes and a double disruptant uL11AB mutant. We applied an integrated approach analyzing a broad range of biological parameters of yeast mutant cells, especially the longevity phenomenon, supplemented with biochemical and high throughput transcriptomic and metobolomic approaches. The analysis showed that the longevity phenomenon is not fully related to the commonly considered energy restriction effect, thus the slow-down of translation does not represent the sole source of aging. Additionally, we showed that uL11 can be classified as a moonlighting protein with extra-ribosomal function having cell-cycle regulatory potential.

Exploring whole-genome duplicate gene retention with complex genetic interaction analysis

Whole-genome duplication has played a central role in the genome evolution of many organisms, including the human genome. Most duplicated genes are eliminated, and factors that influence the retention of persisting duplicates remain poorly understood. We describe a systematic complex genetic interaction analysis with yeast paralogs derived from the whole-genome duplication event. Mapping of digenic interactions for a deletion mutant of each paralog, and of trigenic interactions for the double mutant, provides insight into their roles and a quantitative measure of their functional redundancy. Trigenic interaction analysis distinguishes two classes of paralogs: a more functionally divergent subset and another that retained more functional overlap. Gene feature analysis and modeling suggest that evolutionary trajectories of duplicated genes are dictated by combined functional and structural entanglement factors.

Functional characterization and expression patterns of PnATX genes under different abiotic stress treatments in Populus

The copper chaperone ATX1 has been investigated previously in the herbaceous plants Arabidopsis and rice. However, the molecular mechanisms of ATX1 underlying copper transport and functional characteristics in the woody plant Populus are poorly understood. In this study, PnATX1 and PnATX2 of Populus simonii × P. nigra were identified and characterized. Sequence analysis showed that PnATXs contained the metal-binding motif MXCXXC in the N-terminus and a lysine-rich region. Phylogenetic analysis of ATX protein sequences revealed that PnATXs were clustered in the same group as AtATX1. PnATX proteins were localized in the cytoplasm and nucleus. Tissue-specific expression analysis showed that PnATX1 and PnATX2 were expressed in all analyzed tissues and, in particular, expressed to a higher relative expression level in young leaves. Quantitative real-time PCR analysis indicated that each PnATX gene was differentially expressed in different tissues under treatments with copper, zinc, iron, jasmonate and salicylic acid (SA). The copper-response element GTAC, methyl jasmonate and salicylic acid responsiveness elements and other cis-acting elements were identified in the PnATX1 and PnATX2 promoters. Expression of β-glucuronidase driven by the PnATX1 promoter was observed in the apical meristem of 7-day-old Arabidopsis transgenic seedlings, and the signal strength was not influenced by deficient or excessive copper conditions. Both PnATX1 and PnATX2 functionally rescued the defective phenotypes of yeast atx1Δ and sod1Δ strains. Under copper excess and deficiency conditions, transgenic Arabidopsis atx1 mutants harboring 35S::PnATX constructs exhibited root length and fresh weight similar to those of the wild type and higher than those of Arabidopsis atx1 mutants. Superoxide dismutase activity decreased in transgenic lines compared with that of atx1 mutants, whereas peroxidase and catalase activities increased significantly under excess copper. The results provide a basis for elucidating the role of Populus PnATX genes in copper homeostasis.

Location, location! cellular relocalization primes specialized metabolic diversification

Specialized metabolites are structurally diverse and cell- or tissue-specific molecules produced in restricted plant lineages. In contrast, primary metabolic pathways are highly conserved in plants and produce metabolites essential for all of life, such as amino acids and nucleotides. Substrate promiscuity - the capacity to accept non-native substrates - is a common characteristic of enzymes, and its impact is especially apparent in generating specialized metabolite variation. However, promiscuity only leads to metabolic diversity when alternative substrates are available; thus, enzyme cellular and subcellular localization directly influence chemical phenotypes. We review a variety of mechanisms that modulate substrate availability for promiscuous plant enzymes. We focus on examples where evolution led to modification of the 'cellular context' through changes in cell-type expression, subcellular relocalization, pathway sequestration, and cellular mixing via tissue damage. These varied mechanisms contributed to the emergence of structurally diverse plant specialized metabolites and inform future metabolic engineering approaches.

Paralog buffering contributes to the variable essentiality of genes in cancer cell lines

What makes a gene essential for cellular survival? In model organisms, such as budding yeast, systematic gene deletion studies have revealed that paralog genes are less likely to be essential than singleton genes and that this can partially be attributed to the ability of paralogs to buffer each other's loss. However, the essentiality of a gene is not a fixed property and can vary significantly across different genetic backgrounds. It is unclear to what extent paralogs contribute to this variation, as most studies have analyzed genes identified as essential in a single genetic background. Here, using gene essentiality profiles of 558 genetically heterogeneous tumor cell lines, we analyze the contribution of paralogy to variable essentiality. We find that, compared to singleton genes, paralogs are less frequently essential and that this is more evident when considering genes with multiple paralogs or with highly sequence-similar paralogs. In addition, we find that paralogs derived from whole genome duplication exhibit more variable essentiality than those derived from small-scale duplications. We provide evidence that in 13–17% of cases the variable essentiality of paralogs can be attributed to buffering relationships between paralog pairs, as evidenced by synthetic lethality. Paralog pairs derived from whole genome duplication and pairs that function in protein complexes are significantly more likely to display such synthetic lethal relationships. Overall we find that many of the observations made using a single strain of budding yeast can be extended to understand patterns of essentiality in genetically heterogeneous cancer cell lines.

Somewhat surprisingly, the majority of human genes can be mutated or deleted in individual cell lines without killing the cells. This observation raises a number of questions—which genes can be lost and why? Here we address these questions by analyzing data on which genes are essential for the growth of over 500 cancer cell lines. In general we find that paralog genes are essential in fewer cell lines than genes that are not paralogs. Paralogs are genes that have been duplicated at some point in evolutionary history, resulting in our genome having two copies of the same gene—a paralog pair. These paralog pairs are a potential source of redundancy, similar to a car having a spare tire. If this is the case, one might expect that losing one gene from a paralog pair could be tolerated by cells, due to the existence of a 'backup gene', but losing both members would cause cells to die. By analyzing the cancer cell lines we estimate this to be the case for 13–17% of paralog pairs, and that this provides an explanation for why some genes are essential in some cell lines but not others.

Paralog dependency indirectly affects the robustness of human cells

The protective redundancy of paralogous genes partly relies on the fact that they carry their functions independently. However, a significant fraction of paralogous proteins may form functionally dependent pairs, for instance, through heteromerization. As a consequence, one could expect these heteromeric paralogs to be less protective against deleterious mutations. To test this hypothesis, we examined the robustness landscape of gene loss-of-function by CRISPR-Cas9 in more than 450 human cell lines. This landscape shows regions of greater deleteriousness to gene inactivation as a function of key paralog properties. Heteromeric paralogs are more likely to occupy such regions owing to their high expression and large number of protein-protein interaction partners. Further investigation revealed that heteromers may also be under stricter dosage balance, which may also contribute to the higher deleteriousness upon gene inactivation. Finally, we suggest that physical dependency may contribute to the deleteriousness upon loss-of-function as revealed by the correlation between the strength of interactions between paralogs and their higher deleteriousness upon loss of function.

Functional Analysis of the Ribosomal uL6 Protein of Saccharomyces cerevisiae

The genome-wide duplication event observed in eukaryotes represents an interesting biological phenomenon, extending the biological capacity of the genome at the expense of the same genetic material. For example, most ribosomal proteins in Saccharomyces cerevisiae are encoded by a pair of paralogous genes. It is thought that gene duplication may contribute to heterogeneity of the translational machinery; however, the exact biological function of this event has not been clarified. In this study, we have investigated the functional impact of one of the duplicated ribosomal proteins, uL6, on the translational apparatus together with its consequences for aging of yeast cells. Our data show that uL6 is not required for cell survival, although lack of this protein decreases the rate of growth and inhibits budding. The uL6 protein is critical for the efficient assembly of the ribosome 60S subunit, and the two uL6 isoforms most likely serve the same function, playing an important role in the adaptation of translational machinery performance to the metabolic needs of the cell. The deletion of a single uL6 gene significantly extends the lifespan but only in cells with a high metabolic rate. We conclude that the maintenance of two copies of the uL6 gene enables the cell to cope with the high demands for effective ribosome synthesis.

Long-Term Conservation of Ohnologs Through Partial Tetrasomy Following Whole-Genome Duplication in Salmonidae

Whole-genome duplications (WGDs) have occurred repeatedly and broadly throughout the evolutionary history of eukaryotes. However, the effects of WGD on genome function and evolution remain unclear. The salmonid WGD that occurred approximately 88 million years ago presents an excellent opportunity for studying the effects of WGD as ∼10-15% of each salmonid genome still exhibits tetrasomic inheritance. Herein, we utilized the rainbow trout (Oncorhynchus mykiss) genome assembly and brain transcriptome data to examine the fate of gene pairs (ohnologs) following the salmonid whole-genome duplication. We find higher sequence identity between ohnologs located within known tetrasomic regions than between ohnologs found in disomic regions, and that tetrasomically inherited ohnologs showed greater similarity in patterns of gene expression and per ohnolog were lower expressed, than disomically inherited ohnologs. Enrichment testing for Gene Ontology terms identified 49 over-represented terms in tetrasomically inherited ohnologs compared to disomic ohnologs. However, why these ohnologs are retained as tetrasomic is difficult to answer. It could be that we have identified salmonid specific "dangerous duplicates", that is, genes that cannot take on new roles following WGD. Alternatively, there may be adaptive advantages for retaining genes as functional duplicates in tetrasomic regions, as presumably, movement of these genes into disomic regions would affect both their sequence identity and their gene expression patterns.

Plasticity and transient binding are key ingredients of the periplasmic chaperone network

SurA, Skp, FkpA, and DegP constitute a chaperone network that ensures biogenesis of outer membrane proteins (OMPs) in Gram-negative bacteria. Both Skp and FkpA are holdases that prevent the self-aggregation of unfolded OMPs, whereas SurA accelerates folding and DegP is a protease. None of these chaperones is essential, and we address here how functional plasticity is manifested in nine known null strains. Using a comprehensive computational model of this network termed OMPBioM, our results suggest that a threshold level of steady state holdase occupancy by chaperones is required, but the cell is agnostic to the specific holdase molecule fulfilling this function. In addition to its foldase activity, SurA moonlights as a holdase when there is no expression of Skp and FkpA. We further interrogate the importance of chaperone-client complex lifetime by conducting simulations using lifetime values for Skp complexes that range in length by six orders of magnitude. This analysis suggests that transient occupancy of durations much shorter than the Escherichia coli doubling time is required. We suggest that fleeting chaperone occupancy facilitates rapid sampling of the periplasmic conditions, which ensures that the cell can be adept at responding to environmental changes. Finally, we calculated the network effects of adding multivalency by computing populations that include two Skp trimers per unfolded OMP. We observe only modest perturbations to the system. Overall, this quantitative framework of chaperone-protein interactions in the periplasm demonstrates robust plasticity due to its dynamic binding and unbinding behavior.

Fine-Tuning the Expression of Duplicate Genes by Translational Regulation in Arabidopsis and Maize

Plant genomes are extensively shaped by various types of gene duplication. However, in this active area of investigation, the vast majority of studies focus on the sequence and transcription of duplicate genes, leaving open the question of how translational regulation impacts the expression and evolution of duplicate genes. We explored this issue by analyzing the ribo- and mRNA-seq data sets across six tissue types and stress conditions in Arabidopsis thaliana and maize (Zea mays). We dissected the relative contributions of transcriptional and translational regulation to the divergence in the abundance of ribosome footprint (RF) for different types of duplicate genes. We found that the divergence in RF abundance was largely programmed at the transcription level and that translational regulation plays more of a modulatory role. Intriguingly, translational regulation is characterized by its strong directionality, with the divergence in translational efficiency (TE) globally counteracting the divergence in mRNA abundance, indicating partial buffering of the transcriptional divergence between paralogs by translational regulation. Divergence in TE was associated with several sequence features. The faster-evolving copy in a duplicate pair was more likely to show lower RF abundance, which possibly results from relaxed purifying selection compared with its paralog. A considerable proportion of duplicates displayed differential TE across tissue types and stress conditions, most of which were enriched in photosynthesis, energy production, and translation-related processes. Additionally, we constructed a database TDPDG-DB (http://www.plantdupribo.tk), providing an online platform for data exploration. Overall, our study illustrates the roles of translational regulation in fine-tuning duplicate gene expression in plants.

Systemic signalling and local effectors in developmental stability, body symmetry, and size

Symmetric growth and the origins of fluctuating asymmetry are unresolved phenomena of biology. Small, and sometimes noticeable, deviations from perfect bilateral symmetry reflect the vulnerability of development to perturbations. The degree of asymmetry is related to the magnitude of the perturbations and the ability of an individual to cope with them. As the left and right sides of an individual were presumed to be genetically identical, deviations of symmetry were traditionally attributed to non-genetic effects such as environmental and developmental noise. In this review, we draw attention to other possible sources of variability, especially to somatic mutations and transposons. Mutations are a major source of phenotypic variability and recent genomic data have highlighted somatic mutations as ubiquitous, even in phenotypically normal individuals. We discuss the importance of factors that are responsible for buffering and stabilizing the genome and for maintaining size robustness and quality through elimination of less-fit or damaged cells. However, the important question that arises from these studies is whether this self-correcting capacity and intrinsic organ size controls are sufficient to explain how symmetric structures can reach an identical size and shape. Indeed, recent discoveries in the fruit fly have uncovered a conserved hormone of the insulin/IGF/relaxin family, Dilp8, that is responsible for stabilizing body size and symmetry in the face of growth perturbations. Dilp8 alarm signals periphery growth status to the brain, where it acts on its receptor Lgr3. Loss of Dilp8-Lgr3 signaling renders flies incapable of detecting growth perturbations and thus maintaining a stable size and symmetry. These findings help to understand how size and symmetry of somatic tissues remain undeterred in noisy environments, after injury or illnesses, and in the presence of accumulated somatic mutations.