Transcription control reprogramming in genetic backup circuits

Altered interactions between cis-regulatory elements partially resolve BLADE-ON-PETIOLE genetic redundancy in Capsella rubella

Duplicated genes are thought to follow one of three evolutionary trajectories that resolve their redundancy: neofunctionalization, subfunctionalization, or pseudogenization. Differences in expression patterns have been documented for many duplicated gene pairs and interpreted as evidence of subfunctionalization and a loss of redundancy. However, little is known about the functional impact of such differences and about their molecular basis. Here, we investigate the genetic and molecular basis for the partial loss of redundancy between the two BLADE-ON-PETIOLE genes BOP1 and BOP2 in red shepherd's purse (Capsella rubella) compared to Arabidopsis (Arabidopsis thaliana). While both genes remain almost fully redundant in A. thaliana, BOP1 in C. rubella can no longer ensure wild-type floral organ numbers and suppress bract formation, due to an altered expression pattern in the region of the cryptic bract primordium. We use two complementary approaches, transgenic rescue of A. thaliana atbop1 atbop2 double mutants and deletions in the endogenous AtBOP1 promoter, to demonstrate that several BOP1 promoter regions containing conserved noncoding sequences interact in a nonadditive manner to control BOP1 expression in the bract primordium and that changes in these interactions underlie the evolutionary divergence between C. rubella and A. thaliana BOP1 expression and activity. Similarly, altered interactions between cis-regulatory regions underlie the divergence in functional promoter architecture related to the control of floral organ abscission by BOP1. These findings highlight the complexity of promoter architecture in plants and suggest that changes in the interactions between cis-regulatory elements are key drivers for evolutionary divergence in gene expression and the loss of redundancy.

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.

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.

Tackling redundancy: genetic mechanisms underlying paralog compensation in plants

Gene duplication is a powerful source of biological innovation giving rise to paralogous genes that undergo diverse fates. Redundancy between paralogous genes is an intriguing outcome of duplicate gene evolution, and its maintenance over evolutionary time has long been considered a paradox. Redundancy can also be dubbed 'a geneticist's nightmare': It hinders the predictability of genome editing outcomes and limits our ability to link genotypes to phenotypes. Genetic studies in yeast and plants have suggested that the ability of ancient redundant duplicates to compensate for dosage perturbations resulting from a loss of function depends on the reprogramming of gene expression, a phenomenon known as active compensation. Starting from considerations on the stoichiometric constraints that drive the evolutionary stability of redundancy, this review aims to provide insights into the mechanisms of active compensation between duplicates that could be targeted for breaking paralog dependencies - the next frontier in plant functional studies.

Genetic Channelization Mechanism of Four Chalcone Isomerase Homologous Genes for Synergistic Resistance to Fusarium wilt in Gossypium barbadense L

Duplication events occur very frequently during plant evolution. The genes in the duplicated pathway or network can evolve new functions through neofunctionalization and subfunctionalization. Flavonoids are secondary metabolites involved in plant development and defense. Our previous transcriptomic analysis of F6 recombinant inbred lines (RILs) and the parent lines after Fusarium oxysporum f. sp. vasinfectum (Fov) infection showed that CHI genes have important functions in cotton. However, there are few reports on the possible neofunctionalization differences of CHI family paralogous genes involved in Fusarium wilt resistance in cotton. In this study, the resistance to Fusarium wilt, expression of metabolic pathway-related genes, metabolite content, endogenous hormone content, reactive oxygen species (ROS) content and subcellular localization of four paralogous CHI family genes in cotton were investigated. The results show that the four paralogous CHI family genes may play a synergistic role in Fusarium wilt resistance. These results revealed a genetic channelization mechanism that can regulate the metabolic flux homeostasis of flavonoids under the mediation of endogenous salicylic acid (SA) and methyl jasmonate (MeJA) via the four paralogous CHI genes, thereby achieving disease resistance. Our study provides a theoretical basis for studying the evolutionary patterns of homologous plant genes and using homologous genes for molecular breeding.

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.

A Mouse Model with a Frameshift Mutation in the Nuclear Factor I/X (NFIX) Gene Has Phenotypic Features of Marshall-Smith Syndrome

The nuclear factor I/X (NFIX) gene encodes a ubiquitously expressed transcription factor whose mutations lead to two allelic disorders characterized by developmental, skeletal, and neural abnormalities, namely, Malan syndrome (MAL) and Marshall-Smith syndrome (MSS). NFIX mutations associated with MAL mainly cluster in exon 2 and are cleared by nonsense-mediated decay (NMD) leading to NFIX haploinsufficiency, whereas NFIX mutations associated with MSS are clustered in exons 6-10 and escape NMD and result in the production of dominant-negative mutant NFIX proteins. Thus, different NFIX mutations have distinct consequences on NFIX expression. To elucidate the in vivo effects of MSS-associated NFIX exon 7 mutations, we used CRISPR-Cas9 to generate mouse models with exon 7 deletions that comprised: a frameshift deletion of two nucleotides (Nfix Del2); in-frame deletion of 24 nucleotides (Nfix Del24); and deletion of 140 nucleotides (Nfix Del140). Nfix +/Del2, Nfix +/Del24, Nfix +/Del140, Nfix Del24/Del24, and Nfix Del140/Del140 mice were viable, normal, and fertile, with no skeletal abnormalities, but Nfix Del2/Del2 mice had significantly reduced viability (p < 0.002) and died at 2-3 weeks of age. Nfix Del2 was not cleared by NMD, and NfixDel2/Del2 mice, when compared to Nfix +/+ and Nfix +/Del2 mice, had: growth retardation; short stature with kyphosis; reduced skull length; marked porosity of the vertebrae with decreased vertebral and femoral bone mineral content; and reduced caudal vertebrae height and femur length. Plasma biochemistry analysis revealed Nfix Del2/Del2 mice to have increased total alkaline phosphatase activity but decreased C-terminal telopeptide and procollagen-type-1-N-terminal propeptide concentrations compared to Nfix +/+ and Nfix +/Del2 mice. Nfix Del2/Del2 mice were also found to have enlarged cerebral cortices and ventricular areas but smaller dentate gyrus compared to Nfix +/+ mice. Thus, Nfix Del2/Del2 mice provide a model for studying the in vivo effects of NFIX mutants that escape NMD and result in developmental abnormalities of the skeletal and neural tissues that are associated with MSS. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

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.

The Fate of Duplicated Enzymes in Prokaryotes: The Case of Isomerases

The isomerases are a unique enzymatic class of enzymes that carry out a great diversity of chemical reactions at the intramolecular level. This class comprises about 300 members, most of which are involved in carbohydrate and terpenoid/polyketide metabolism. Along with oxidoreductases and translocases, isomerases are one of the classes with the highest ratio of paralogous enzymes. Due to its relatively small number of members, it is plausible to explore it in greater detail to identify specific cases of gene duplication. Here, we present an analysis at the level of individual isomerases and identify different members that seem to be involved in duplication events in prokaryotes. As was suggested in a previous study, there is no homogeneous distribution of paralogs, but rather they accumulate into a few subcategories, some of which differ between Archaea and Bacteria. As expected, the metabolic processes with more paralogous isomerases have to do with carbohydrate metabolism but also with RNA modification (a particular case involving an rRNA-modifying isomerase is thoroughly discussed and analyzed in detail). Overall, our findings suggest that the most common fate for paralogous enzymes is the retention of the original enzymatic function, either associated with a dosage effect or with differential expression in response to changing environments, followed by subfunctionalization and, to a much lesser degree, neofunctionalization, which is consistent with what has been reported elsewhere.

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.

The HAPSTR2 retrogene buffers stress signaling and resilience in mammals

We recently identified HAPSTR1 (C16orf72) as a key component in a novel pathway which regulates the cellular response to molecular stressors, such as DNA damage, nutrient scarcity, and protein misfolding. Here, we identify a functional paralog to HAPSTR1: HAPSTR2. HAPSTR2 formed early in mammalian evolution, via genomic integration of a reverse transcribed HAPSTR1 transcript, and has since been preserved under purifying selection. HAPSTR2, expressed primarily in neural and germline tissues and a subset of cancers, retains established biochemical features of HAPSTR1 to achieve two functions. In normal physiology, HAPSTR2 directly interacts with HAPSTR1, markedly augmenting HAPSTR1 protein stability in a manner independent from HAPSTR1's canonical E3 ligase, HUWE1. Alternatively, in the context of HAPSTR1 loss, HAPSTR2 expression is sufficient to buffer stress signaling and resilience. Thus, we discover a mammalian retrogene which safeguards fitness.

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.

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.

Dynamic evolution of small signalling peptide compensation in plant stem cell control

Gene duplications are a hallmark of plant genome evolution and a foundation for genetic interactions that shape phenotypic diversity^1-5. Compensation is a major form of paralogue interaction^6-8 but how compensation relationships change as allelic variation accumulates is unknown. Here we leveraged genomics and genome editing across the Solanaceae family to capture the evolution of compensating paralogues. Mutations in the stem cell regulator CLV3 cause floral organs to overproliferate in many plants^9-11. In tomato, this phenotype is partially suppressed by transcriptional upregulation of a closely related paralogue¹². Tobacco lost this paralogue, resulting in no compensation and extreme clv3 phenotypes. Strikingly, the paralogues of petunia and groundcherry nearly completely suppress clv3, indicating a potent ancestral state of compensation. Cross-species transgenic complementation analyses show that this potent compensation partially degenerated in tomato due to a single amino acid change in the paralogue and cis-regulatory variation that limits its transcriptional upregulation. Our findings show how genetic interactions are remodelled following duplications and suggest that dynamic paralogue evolution is widespread over short time scales and impacts phenotypic variation from natural and engineered mutations.

Functional buffering via cell-specific gene expression promotes tissue homeostasis and cancer robustness

Functional buffering that ensures biological robustness is critical for maintaining tissue homeostasis, organismal survival, and evolution of novelty. However, the mechanism underlying functional buffering, particularly in multicellular organisms, remains largely elusive. Here, we proposed that functional buffering can be mediated via expression of buffering genes in specific cells and tissues, by which we named Cell-specific Expression-BUffering (CEBU). We developed an inference index (C-score) for CEBU by computing C-scores across 684 human cell lines using genome-wide CRISPR screens and transcriptomic RNA-seq. We report that C-score-identified putative buffering gene pairs are enriched for members of the same duplicated gene family, pathway, and protein complex. Furthermore, CEBU is especially prevalent in tissues of low regenerative capacity (e.g., bone and neuronal tissues) and is weakest in highly regenerative blood cells, linking functional buffering to tissue regeneration. Clinically, the buffering capacity enabled by CEBU can help predict patient survival for multiple cancers. Our results suggest CEBU as a potential buffering mechanism contributing to tissue homeostasis and cancer robustness in humans.

In rice splice variants that restore the reading frame after frameshifting indel introduction are common, often induced by the indels and sometimes lead to organism-level rescue

The introduction of frameshifting non-3n indels enables the identification of gene-trait associations. However, it has been hypothesised that recovery of the original reading frame owing to usage of non-canonical splice forms could cause rescue. To date there is very little evidence for organism-level rescue by such a mechanism and it is unknown how commonly indels induce, or are otherwise associated with, frame-restoring splice forms. We perform CRISPR/Cas9 editing of randomly selected loci in rice to investigate these issues. We find that the majority of loci have a frame-restoring isoform. Importantly, three quarters of these isoforms are not seen in the absence of the indels, consistent with indels commonly inducing novel isoforms. This is supported by analysis in the context of NMD knockdowns. We consider in detail the two top rescue candidates, in wax deficient anther 1 (wda1) and brittle culm (bc10), finding that organismal-level rescue in both cases is strong but owing to different splice modification routes. More generally, however, as frame-restoring isoforms are low abundance and possibly too disruptive, such rescue we suggest to be the rare exception, not the rule. Nonetheless, assuming that indels commonly induce frame-restoring isoforms, these results emphasize the need to examine RNA level effects of non-3n indels and suggest that multiple non-3n indels in any given gene are advisable to probe a gene’s trait associations.

As protein coding genes are read in units of three (codons), insertions or deletions (indels) that are not a multiple of three long (non 3n) are expected to be especially harmful. Whether they are is important both for interpreting the results of non-3n indel experiments to probe a gene’s functional importance and for diagnostics. Particularly enigmatic are incidences where some non-3n changes in a gene compromise phenotypes while other seemingly comparable ones do not. One explanation for the latter is that a non-3n indel might be rescued via a frame-restoring splice form. Here we examine this hypothesis by inducing non-3n indels in many genes in rice and find that many non-3n indels are associated with a splice form that restores the reading frame. In the majority of these cases the indel appears to induce the potential rescuing splice form. We examine two top hit cases in detail and show functional rescue by splice modification. More generally, the frame-restoring forms are, however, low abundance and probably result in compromised proteins. We conclude then that splice mediated rescue is possible, but probably uncommon. Nonetheless it should not be overlooked in experimental design and interpretation.

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.

Repurposing of the enhancer-promoter communication underlies the compensation of Mesp2 by Mesp1

Organisms are inherently equipped with buffering systems against genetic perturbations. Genetic compensation, the compensatory response by upregulating another gene or genes, is one such buffering mechanism. Recently, a well-conserved compensatory mechanism was proposed: transcriptional adaptation of homologs under the nonsense-mediated mRNA decay pathways. However, this model cannot explain the onset of all compensatory events. We report a novel genetic compensation mechanism operating over the Mesp gene locus. Mesp1 and Mesp2 are paralogs located adjacently in the genome. Mesp2 loss is partially rescued by Mesp1 upregulation in the presomitic mesoderm (PSM). Using a cultured PSM induction system, we reproduced the compensatory response in vitro and found that the Mesp2-enhancer is required to promote Mesp1. We revealed that the Mesp2-enhancer directly interacts with the Mesp1 promoter, thereby upregulating Mesp1 expression upon the loss of Mesp2. Of note, this interaction is established by genomic arrangement upon PSM development independently of Mesp2 disruption. We propose that the repurposing of this established enhancer-promoter communication is the mechanism underlying this compensatory response for the upregulation of the adjacent gene.

Genetic compensation, the compensatory response by upregulating another gene or genes, is one of the inherent mechanisms against gene disruption to confer cellular fitness. However, the regulatory mechanisms are largely unknown. Nonsense-mediated mutant mRNA degradation was recently proposed as a conserved mechanism across species to upregulate homologous genes to compensate for a disrupted gene, but this cannot explain compensation events with no mutant mRNA. This study investigated the compensation mechanism operating over adjacent paralogs, Mesp1 and Mesp2, in the genome. Mesp genes encode essential transcription factors in the presomitic mesoderm for development. In general, an enhancer is considered to activate a target gene when it physically interacts with the target. The communication of the Mesp2-enhancer with the Mesp1 promoter is established upon differentiation of the presomitic mesoderm, but this communication activates Mesp1 only when Mesp2 is disrupted, leading to compensation. We revealed a novel compensation mechanism depending on the repurposing of this enhancer-promoter communication by gene disruption. Our study also provides new insight into transcriptional regulation by providing the concept that an enhancer changes its target even among its physically interacting genes in a context-dependent manner.

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.

RNA chaperone activates Salmonella virulence program during infection

Organisms often harbor seemingly redundant proteins. In the bacterium Salmonella enterica serovar Typhimurium (S. Typhimurium), the RNA chaperones CspC and CspE appear to play redundant virulence roles because a mutant lacking both chaperones is attenuated, whereas mutants lacking only one exhibit wild-type virulence. We now report that CspC—but not CspE—is necessary to activate the master virulence regulator PhoP when S. Typhimurium experiences mildly acidic pH, such as inside macrophages. This CspC-dependent PhoP activation is specific to mildly acidic pH because a cspC mutant behaves like wild-type S. Typhimurium under other PhoP-activating conditions. Moreover, it is mediated by ugtL, a virulence gene required for PhoP activation inside macrophages. Purified CspC promotes ugtL translation by disrupting a secondary structure in the ugtL mRNA that occludes ugtL’s ribosome binding site. Our findings demonstrate that proteins that are seemingly redundant actually confer distinct and critical functions to the lifestyle of an organism.

Dissecting the Fitness Costs of Complex Mutations

The fitness cost of complex pleiotropic mutations is generally difficult to assess. On the one hand, it is necessary to identify which molecular properties are directly altered by the mutation. On the other, this alteration modifies the activity of many genetic targets with uncertain consequences. Here, we examine the possibility of addressing these challenges by identifying unique predictors of these costs. To this aim, we consider mutations in the RNA polymerase (RNAP) in Escherichia coli as a model of complex mutations. Changes in RNAP modify the global program of transcriptional regulation, with many consequences. Among others is the difficulty to decouple the direct effect of the mutation from the response of the whole system to such mutation. A problem that we solve quantitatively with data of a set of constitutive genes, those on which the global program acts most directly. We provide a statistical framework that incorporates the direct effects and other molecular variables linked to this program as predictors, which leads to the identification that some genes are more suitable to determine costs than others. Therefore, we not only identified which molecular properties best anticipate fitness, but we also present the paradoxical result that, despite pleiotropy, specific genes serve as more solid predictors. These results have connotations for the understanding of the architecture of robustness in biological systems.

Functional Roles of Homologous Recombination and Non-Homologous End Joining in DNA Damage Response and Microevolution in Cryptococcus neoformans

DNA double-strand breaks (DSBs) are the most deleterious type of DNA lesions because they cause loss of genetic information if not properly repaired. In eukaryotes, homologous recombination (HR) and non-homologous end joining (NHEJ) are required for DSB repair. However, the relationship of HR and NHEJ in DNA damage stress is unknown in the radiation-resistant fungus Cryptococcus neoformans. In this study, we found that the expression levels of HR- and NHEJ-related genes were highly induced in a Rad53-Bdr1 pathway-dependent manner under genotoxic stress. Deletion of RAD51, which is one of the main components in the HR, resulted in growth under diverse types of DNA damage stress, whereas perturbations of KU70 and KU80, which belong to the NHEJ system, did not affect the genotoxic stresses except when bleomycin was used for treatment. Furthermore, deletion of both RAD51 and KU70/80 renders cells susceptible to oxidative stress. Notably, we found that deletion of RAD51 induced a hypermutator phenotype in the fluctuation assay. In contrast to the fluctuation assay, perturbation of KU70 or KU80 induced rapid microevolution similar to that induced by the deletion of RAD51. Collectively, Rad51-mediated HR and Ku70/Ku80-mediated NHEJ regulate the DNA damage response and maintain genome stability.

Suboptimal Global Transcriptional Response Increases the Harmful Effects of Loss-of-Function Mutations

The fitness impact of loss-of-function mutations is generally assumed to reflect the loss of specific molecular functions associated with the perturbed gene. Here, we propose that rewiring of the transcriptome upon deleterious gene inactivation is frequently nonspecific and mimics stereotypic responses to external environmental change. Consequently, transcriptional response to gene deletion could be suboptimal and incur an extra fitness cost. Analysis of the transcriptomes of ∼1,500 single-gene deletion Saccharomyces cerevisiae strains supported this scenario. First, most transcriptomic changes are not specific to the deleted gene but are rather triggered by perturbations in functionally diverse genes. Second, gene deletions that alter the expression of dosage-sensitive genes are especially harmful. Third, by elevating the expression level of downregulated genes, we could experimentally mitigate the fitness defect of gene deletions. Our work shows that rewiring of genomic expression upon gene inactivation shapes the harmful effects of mutations.

Deletion of bglC triggers a genetic compensation response by awakening the expression of alternative beta-glucosidase

In the plant pathogen Streptomyces scabies, the gene bglC encodes a GH1 family cellobiose beta-glucosidase that is both required for primary metabolism and for inducing virulence of the bacterium. Deletion of bglC (strain ΔbglC) surprisingly resulted in the augmentation of the global beta-glucosidase activity of S. scabies. This paradoxical phenotype is highly robust as it has been observed in all bglC deletion mutants independently generated, thereby highlighting a phenomenon of genetic compensation. Comparative proteomics allowed to identify two glycosyl hydrolases - named BcpE1 and BcpE2 - of which peptide levels were significantly increased in strain ΔbglC. Quantitative RT-PCR revealed that the higher abundance of BcpE1 and BcpE2 is triggered at the transcriptional level, the expression of their respective gene being 100 and 15 times upregulated. Enzymatic studies with pure BcpE proteins showed that they both possess beta-glucosidase activity thereby explaining the genotypic-phenotypic discrepancy of the bglC deletion mutant. The GH1 family BcpE1 could hydrolyze cellobiose and generate glucose similarly to BglC itself thereby mainly contributing to the survival of strain ΔbglC when cellobiose is provided as sole nutrient source. The low affinity of BcpE2 for cellobiose suggests that this GH3 family beta-glucosidase would instead primarily target another and yet unknown glucose-beta-1,4-linked substrate. These results make S. scabies a new model system to study genetic compensation. Discovering how, either the bglC DNA locus, its mRNA, the BglC protein, or either its enzymatic activity controls bcpE genes' expression, will unveil new mechanisms directing transcriptional repression.

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.

Interactions among ryanodine receptor isotypes contribute to muscle fiber type development and function

Mutations affecting ryanodine receptor (RyR) calcium release channels commonly underlie congenital myopathies. Although these channels are known principally for their essential roles in muscle contractility, mutations in the human RYR1 gene result in a broad spectrum of phenotypes, including muscle weakness, altered proportions of fiber types, anomalous muscle fibers with cores or centrally placed nuclei, and dysmorphic craniofacial features. Currently, it is unknown which phenotypes directly reflect requirements for RyRs and which result secondarily to aberrant muscle function. To identify biological processes requiring RyR function, skeletal muscle development was analyzed in zebrafish embryos harboring protein-null mutations. RyR channels contribute to both muscle fiber development and function. Loss of some RyRs had modest effects, altering muscle fiber-type specification in the embryo without compromising viability. In addition, each RyR-encoding gene contributed to normal swimming behavior and muscle function. The RyR channels do not function in a simple additive manner. For example, although isoform RyR1a is sufficient for muscle contraction in the absence of RyR1b, RyR1a normally attenuates the activity of the co-expressed RyR1b channel in slow muscle. RyR3 also acts to modify the functions of other RyR channels. Furthermore, diminished RyR-dependent contractility affects both muscle fiber maturation and craniofacial development. These findings help to explain some of the heterogeneity of phenotypes that accompany RyR1 mutations in humans.

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.

Evolution of buffering in a genetic circuit controlling plant stem cell proliferation

Precise control of plant stem cell proliferation is necessary for the continuous and reproducible development of plant organs^1,2. The peptide ligand CLAVATA3 (CLV3) and its receptor protein kinase CLAVATA1 (CLV1) maintain stem cell homeostasis within a deeply conserved negative feedback circuit^1,2. In Arabidopsis, CLV1 paralogs also contribute to homeostasis, by compensating for the loss of CLV1 through transcriptional upregulation³. Here, we show that compensation^4,5 operates in diverse lineages for both ligands and receptors, but while the core CLV signaling module is conserved, compensation mechanisms have diversified. Transcriptional compensation between ligand paralogs operates in tomato, facilitated by an ancient gene duplication that impacted the domestication of fruit size. In contrast, we found little evidence for transcriptional compensation between ligands in Arabidopsis and maize, and receptor compensation differs between tomato and Arabidopsis. Our findings show that compensation among ligand and receptor paralogs is critical for stem cell homeostasis, but that diverse genetic mechanisms buffer conserved developmental programs.

Genome-wide prediction of synthetic rescue mediators of resistance to targeted and immunotherapy

Most patients with advanced cancer eventually acquire resistance to targeted therapies, spurring extensive efforts to identify molecular events mediating therapy resistance. Many of these events involve synthetic rescue (SR) interactions, where the reduction in cancer cell viability caused by targeted gene inactivation is rescued by an adaptive alteration of another gene (the rescuer). Here, we perform a genome-wide in silico prediction of SR rescuer genes by analyzing tumor transcriptomics and survival data of 10,000 TCGA cancer patients. Predicted SR interactions are validated in new experimental screens. We show that SR interactions can successfully predict cancer patients' response and emerging resistance. Inhibiting predicted rescuer genes sensitizes resistant cancer cells to therapies synergistically, providing initial leads for developing combinatorial approaches to overcome resistance proactively. Finally, we show that the SR analysis of melanoma patients successfully identifies known mediators of resistance to immunotherapy and predicts novel rescuers.

Minor Isozymes Tailor Yeast Metabolism to Carbon Availability

Gene duplication is one of the main evolutionary paths to new protein function. Typically, duplicated genes either accumulate mutations and degrade into pseudogenes or are retained and diverge in function. Some duplicated genes, however, show long-term persistence without apparently acquiring new function. An important class of isozymes consists of those that catalyze the same reaction in the same compartment, where knockout of one isozyme causes no known functional defect. Here we present an approach to assigning specific functional roles to seemingly redundant isozymes. First, gene expression data are analyzed computationally to identify conditions under which isozyme expression diverges. Then, knockouts are compared under those conditions. This approach revealed that the expression of many yeast isozymes diverges in response to carbon availability and that carbon source manipulations can induce fitness phenotypes for seemingly redundant isozymes. A driver of these fitness phenotypes is differential allosteric enzyme regulation, indicating isozyme divergence to achieve more-optimal control of metabolism.

Isozymes are enzymes that differ in sequence but catalyze the same chemical reactions. Despite their apparent redundancy, isozymes are often retained over evolutionary time, suggesting that they contribute to fitness. We developed an unsupervised computational method for identifying environmental conditions under which isozymes are likely to make fitness contributions. This method analyzes published gene expression data to find specific experimental perturbations that induce differential isozyme expression. In yeast, we found that isozymes are strongly enriched in the pathways of central carbon metabolism and that many isozyme pairs show anticorrelated expression during the respirofermentative shift. Building on these observations, we assigned function to two minor central carbon isozymes, aconitase 2 (ACO2) and pyruvate kinase 2 (PYK2). ACO2 is expressed during fermentation and proves advantageous when glucose is limiting. PYK2 is expressed during respiration and proves advantageous for growth on three-carbon substrates. PYK2’s deletion can be rescued by expressing the major pyruvate kinase only if that enzyme carries mutations mirroring PYK2’s allosteric regulation. Thus, central carbon isozymes help to optimize allosteric metabolic regulation under a broad range of potential nutrient conditions while requiring only a small number of transcriptional states.

IMPORTANCE Gene duplication is one of the main evolutionary paths to new protein function. Typically, duplicated genes either accumulate mutations and degrade into pseudogenes or are retained and diverge in function. Some duplicated genes, however, show long-term persistence without apparently acquiring new function. An important class of isozymes consists of those that catalyze the same reaction in the same compartment, where knockout of one isozyme causes no known functional defect. Here we present an approach to assigning specific functional roles to seemingly redundant isozymes. First, gene expression data are analyzed computationally to identify conditions under which isozyme expression diverges. Then, knockouts are compared under those conditions. This approach revealed that the expression of many yeast isozymes diverges in response to carbon availability and that carbon source manipulations can induce fitness phenotypes for seemingly redundant isozymes. A driver of these fitness phenotypes is differential allosteric enzyme regulation, indicating isozyme divergence to achieve more-optimal control of metabolism.

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.

Spatio-temporal expression dynamics differ between homologues of flowering time genes in the allopolyploid Brassica napus

Polyploidy is a recurrent feature of eukaryotic evolution and has been linked to increases in complexity, adaptive radiation and speciation. Within angiosperms such events have occurred repeatedly in many plant lineages. Here we investigate the retention and spatio-temporal expression dynamics of duplicated genes predicted to regulate the floral transition in Brassica napus (oilseed rape, OSR). We show that flowering time genes are preferentially retained relative to other genes in the OSR genome. Using a transcriptome time series in two tissues (leaf and shoot apex) across development we show that 67% of these retained flowering time genes are expressed. Furthermore, between 64% (leaf) and 74% (shoot apex) of the retained gene homologues show diverged expression patterns relative to each other across development, suggesting neo- or subfunctionalization. A case study of homologues of the shoot meristem identity gene TFL1 reveals differences in cis-regulatory elements that could explain this divergence. Such differences in the expression dynamics of duplicated genes highlight the challenges involved in translating gene regulatory networks from diploid model systems to more complex polyploid crop species.

Genetic Interaction Network as an Important Determinant of Gene Order in Genome Evolution

Although it is generally accepted that eukaryotic gene order is not random, the basic principles of gene arrangement on a chromosome remain poorly understood. Here, we extended existing population genetics theories that were based on two-locus models and proposed a hypothesis that genetic interaction networks drive the evolution of eukaryotic gene order. We predicted that genes with positive epistasis would move toward each other in evolution, during which a negative correlation between epistasis and gene distance formed. We tested and confirmed our prediction with computational simulations and empirical data analyses. Importantly, we demonstrated that gene order in the budding yeast could be successfully predicted from the genetic interaction network. Taken together, our study reveals the role of the genetic interaction network in the evolution of gene order, extends our understanding of the encoding principles in genomes, and potentially offers new strategies to improve synthetic biology.

Gene Duplication and Transference of Function in the paleoAP3 Lineage of Floral Organ Identity Genes

The floral organ identity gene APETALA3 (AP3) is a MADS-box transcription factor involved in stamen and petal identity that belongs to the B-class of the ABC model of flower development. Thalictrum (Ranunculaceae), an emerging model in the non-core eudicots, has AP3 homologs derived from both ancient and recent gene duplications. Prior work has shown that petals have been lost repeatedly and independently in Ranunculaceae in correlation with the loss of a specific AP3 paralog, and Thalictrum represents one of these instances. The main goal of this study was to conduct a functional analysis of the three AP3 orthologs present in Thalictrum thalictroides, representing the paleoAP3 gene lineage, to determine the degree of redundancy versus divergence after gene duplication. Because Thalictrum lacks petals, and has lost the petal-specific AP3, we also asked whether heterotopic expression of the remaining AP3 genes contributes to the partial transference of petal function to the first whorl found in insect-pollinated species. To address these questions, we undertook functional characterization by virus-induced gene silencing (VIGS), protein-protein interaction and binding site analyses. Our results illustrate partial redundancy among Thalictrum AP3s, with deep conservation of B-class function in stamen identity and a novel role in ectopic petaloidy of sepals. Certain aspects of petal function of the lost AP3 locus have apparently been transferred to the other paralogs. A novel result is that the protein products interact not only with each other, but also as homodimers. Evidence presented here also suggests that expression of the different ThtAP3 paralogs is tightly integrated, with an apparent disruption of B function homeostasis upon silencing of one of the paralogs that codes for a truncated protein. To explain this result, we propose two testable alternative scenarios: that the truncated protein is a dominant negative mutant or that there is a compensational response as part of a back-up circuit. The evidence for promiscuous protein-protein interactions via yeast two-hybrid combined with the detection of AP3 specific binding motifs in all B-class gene promoters provide partial support for these hypotheses.

Genome Mutational and Transcriptional Hotspots Are Traps for Duplicated Genes and Sources of Adaptations

Gene duplication generates new genetic material, which has been shown to lead to major innovations in unicellular and multicellular organisms. A whole-genome duplication occurred in the ancestor of Saccharomyces yeast species but 92% of duplicates returned to single-copy genes shortly after duplication. The persisting duplicated genes in Saccharomyces led to the origin of major metabolic innovations, which have been the source of the unique biotechnological capabilities in the Baker’s yeast Saccharomyces cerevisiae. What factors have determined the fate of duplicated genes remains unknown. Here, we report the first demonstration that the local genome mutation and transcription rates determine the fate of duplicates. We show, for the first time, a preferential location of duplicated genes in the mutational and transcriptional hotspots of S. cerevisiae genome. The mechanism of duplication matters, with whole-genome duplicates exhibiting different preservation trends compared to small-scale duplicates. Genome mutational and transcriptional hotspots are rich in duplicates with large repetitive promoter elements. Saccharomyces cerevisiae shows more tolerance to deleterious mutations in duplicates with repetitive promoter elements, which in turn exhibit higher transcriptional plasticity against environmental perturbations. Our data demonstrate that the genome traps duplicates through the accelerated regulatory and functional divergence of their gene copies providing a source of novel adaptations in yeast.

CRISPR Technology Reveals RAD(51)-ical Mechanisms of Repair in Roundworms: An Educational Primer for Use with "Promotion of Homologous Recombination by SWS-1 in Complex with RAD-51 Paralogs in Caenorhabditis elegans"

The mechanisms cells use to maintain genetic fidelity via DNA repair and the accuracy of these processes have garnered interest from scientists engaged in basic research to clinicians seeking improved treatment for cancer patients. Despite the continued advances, many details of DNA repair are still incompletely understood. In addition, the inherent complexity of DNA repair processes, even at the most fundamental level, makes it a challenging topic. This primer is meant to assist both educators and students in using a recent paper, "Promotion of homologous recombination by SWS-1 in complex with RAD-51 paralogs in Caenorhabditis elegans," to understand mechanisms of DNA repair. The goals of this primer are to highlight and clarify several key techniques utilized, with special emphasis on the clustered, regularly interspaced, short palindromic repeats technique and the ways in which it has revolutionized genetics research, as well as to provide questions for deeper in-class discussion.

Arabidopsis CP12 mutants have reduced levels of phosphoribulokinase and impaired function of the Calvin-Benson cycle

CP12 is a small, redox-sensitive protein, the most detailed understanding of which is the thioredoxin-mediated regulation of the Calvin-Benson cycle, where it facilitates the formation of a complex between glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and phosphoribulokinase (PRK) in response to changes in light intensity. In most organisms, CP12 proteins are encoded by small multigene families, where the importance of each individual CP12 gene in vivo has not yet been reported. We used Arabidopsis thaliana T-DNA mutants and RNAi transgenic lines with reduced levels of CP12 transcript to determine the relative importance of each of the CP12 genes. We found that single cp12-1, cp12-2, and cp12-3 mutants do not develop a severe photosynthetic or growth phenotype. In contrast, reductions of both CP12-1 and CP12-2 transcripts lead to reductions in photosynthetic capacity and to slower growth and reduced seed yield. No clear phenotype for CP12-3 was evident. Additionally, the levels of PRK protein are reduced in the cp12-1, cp12-1/2, and multiple mutants. Our results suggest that there is functional redundancy between CP12-1 and CP12-2 in Arabidopsis where these proteins have a role in determining the level of PRK in mature leaves and hence photosynthetic capacity.

Upon Accounting for the Impact of Isoenzyme Loss, Gene Deletion Costs Anticorrelate with Their Evolutionary Rates

System-level metabolic network models enable the computation of growth and metabolic phenotypes from an organism's genome. In particular, flux balance approaches have been used to estimate the contribution of individual metabolic genes to organismal fitness, offering the opportunity to test whether such contributions carry information about the evolutionary pressure on the corresponding genes. Previous failure to identify the expected negative correlation between such computed gene-loss cost and sequence-derived evolutionary rates in Saccharomyces cerevisiae has been ascribed to a real biological gap between a gene's fitness contribution to an organism "here and now" and the same gene's historical importance as evidenced by its accumulated mutations over millions of years of evolution. Here we show that this negative correlation does exist, and can be exposed by revisiting a broadly employed assumption of flux balance models. In particular, we introduce a new metric that we call "function-loss cost", which estimates the cost of a gene loss event as the total potential functional impairment caused by that loss. This new metric displays significant negative correlation with evolutionary rate, across several thousand minimal environments. We demonstrate that the improvement gained using function-loss cost over gene-loss cost is explained by replacing the base assumption that isoenzymes provide unlimited capacity for backup with the assumption that isoenzymes are completely non-redundant. We further show that this change of the assumption regarding isoenzymes increases the recall of epistatic interactions predicted by the flux balance model at the cost of a reduction in the precision of the predictions. In addition to suggesting that the gene-to-reaction mapping in genome-scale flux balance models should be used with caution, our analysis provides new evidence that evolutionary gene importance captures much more than strict essentiality.

Mitochondrial-Nuclear Interactions Mediate Sex-Specific Transcriptional Profiles in Drosophila

The assembly and function of mitochondria require coordinated expression from two distinct genomes, the mitochondrial DNA (mtDNA) and nuclear DNA (nDNA). Mutations in either genome can be a source of phenotypic variation, yet their coexpression has been largely overlooked as a source of variation, particularly in the emerging paradigm of mitochondrial replacement therapy. Here we tested how the transcriptome responds to mtDNA and nDNA variation, along with mitonuclear interactions (mtDNA × nDNA) in Drosophila melanogaster. We used two mtDNA haplotypes that differ in a substantial number of single nucleotide polymorphisms, with >100 amino acid differences. We placed each haplotype on each of two D. melanogaster nuclear backgrounds and tested for transcription differences in both sexes. We found that large numbers of transcripts were differentially expressed between nuclear backgrounds, and that mtDNA type altered the expression of nDNA genes, suggesting a retrograde, trans effect of mitochondrial genotype. Females were generally more sensitive to genetic perturbation than males, and males demonstrated an asymmetrical effect of mtDNA in each nuclear background; mtDNA effects were nuclear-background specific. mtDNA-sensitive genes were not enriched in male- or female-limited expression space in either sex. Using a variety of differential expression analyses, we show the responses to mitonuclear covariation to be substantially different between the sexes, yet the mtDNA genes were consistently differentially expressed across nuclear backgrounds and sexes. Our results provide evidence that the main mtDNA effects can be consistent across nuclear backgrounds, but the interactions between mtDNA and nDNA can lead to sex-specific global transcript responses.

Complementary Roles of GADD34- and CReP-Containing Eukaryotic Initiation Factor 2α Phosphatases during the Unfolded Protein Response

Phosphorylation of eukaryotic initiation factor 2α (eIF2α) controls transcriptome-wide changes in mRNA translation in stressed cells. While phosphorylated eIF2α (P-eIF2α) attenuates global protein synthesis, mRNAs encoding stress proteins are more efficiently translated. Two eIF2α phosphatases, containing GADD34 and CReP, catalyze P-eIF2α dephosphorylation. The current view of GADD34, whose transcription is stress induced, is that it functions in a feedback loop to resolve cell stress. In contrast, CReP, which is constitutively expressed, controls basal P-eIF2α levels in unstressed cells. Our studies show that GADD34 drives substantial changes in mRNA translation in unstressed cells, particularly targeting the secretome. Following activation of the unfolded protein response (UPR), rapid translation of GADD34 mRNA occurs and GADD34 is essential for UPR progression. In the absence of GADD34, eIF2α phosphorylation is persistently enhanced and the UPR translational program is significantly attenuated. This "stalled" UPR is relieved by the subsequent activation of compensatory mechanisms that include AKT-mediated suppression of PKR-like kinase (PERK) and increased expression of CReP mRNA, partially restoring protein synthesis. Our studies highlight the coordinate regulation of UPR by the GADD34- and CReP-containing eIF2α phosphatases to control cell viability.

Biological causal links on physiological and evolutionary time scales

Correlation does not imply causation. If two variables, say A and B, are correlated, it could be because A causes B, or that B causes A, or because a third factor affects them both. We suggest that in many cases in biology, the causal link might be bi-directional: A causes B through a fast-acting physiological process, while B causes A through a slowly accumulating evolutionary process. Furthermore, many trained biologists tend to consistently focus at first on the fast-acting direction, and overlook the slower process in the opposite direction. We analyse several examples from modern biology that demonstrate this bias (codon usage optimality and gene expression, gene duplication and genetic dispensability, stem cell division and cancer risk, and the microbiome and host metabolism) and also discuss an example from linguistics. These examples demonstrate mutual effects between the fast physiological processes and the slow evolutionary ones. We believe that building awareness of inference biases among biologists who tend to prefer one causal direction over another could improve scientific reasoning.

A tale of two paralogs: human Transformer2 proteins with differential RNA-binding affinities

The Transformer2 (Tra2) proteins in humans are homologues of the Drosophila Tra2 protein. One of the two RNA-binding paralogs, Tra2β, has been very well-studied over the past decade, but not much is known about Tra2α. It was very recently shown that the two proteins demonstrate the phenomenon of paralog compensation. Here, we provide a structural basis for this genetic backup circuit, using molecular modelling and dynamics studies. We show that the two proteins display similar binding specificities, but differential affinities to a short GAA-rich RNA stretch. Starting from the 6-nucleotide RNA in the solution structure, close to 4000 virtual mutations were modelled on RNA and the domain-RNA interactions were studied after energy minimisation to convergence. Separately, another known 13-nucleotide stretch was docked and the domain-RNA interactions were observed through a 100-ns dynamics trajectory. We have also demonstrated the 'compensatory' mechanism at the level of domains in one of the domain repeat-containing RNA-binding proteins.

Mechanisms of mutational robustness in transcriptional regulation

Robustness is the invariance of a phenotype in the face of environmental or genetic change. The phenotypes produced by transcriptional regulatory circuits are gene expression patterns that are to some extent robust to mutations. Here we review several causes of this robustness. They include robustness of individual transcription factor binding sites, homotypic clusters of such sites, redundant enhancers, transcription factors, redundant transcription factors, and the wiring of transcriptional regulatory circuits. Such robustness can either be an adaptation by itself, a byproduct of other adaptations, or the result of biophysical principles and non-adaptive forces of genome evolution. The potential consequences of such robustness include complex regulatory network topologies that arise through neutral evolution, as well as cryptic variation, i.e., genotypic divergence without phenotypic divergence. On the longest evolutionary timescales, the robustness of transcriptional regulation has helped shape life as we know it, by facilitating evolutionary innovations that helped organisms such as flowering plants and vertebrates diversify.

The origins of mutational robustness

Biological systems are resistant to genetic changes; a property known as mutational robustness, the origin of which remains an open question. In recent years, researchers have explored emergent properties of biological systems and mechanisms of genetic redundancy to reveal how mutational robustness emerges and persists. Several mechanisms have been proposed to explain the origin of mutational robustness, including molecular chaperones and gene duplication. The latter has received much attention, but its role in robustness remains controversial. Here, I examine recent findings linking genetic redundancy through gene duplication and mutational robustness. Experimental evolution and genome resequencing have made it possible to test the role of gene duplication in tolerating mutations at both the coding and regulatory levels. This evidence as well as previous findings on regulatory reprogramming of duplicates support the role of gene duplication in the origin of robustness.

Predicting human genetic interactions from cancer genome evolution

Synthetic Lethal (SL) genetic interactions play a key role in various types of biological research, ranging from understanding genotype-phenotype relationships to identifying drug-targets against cancer. Despite recent advances in empirical measuring SL interactions in human cells, the human genetic interaction map is far from complete. Here, we present a novel approach to predict this map by exploiting patterns in cancer genome evolution. First, we show that empirically determined SL interactions are reflected in various gene presence, absence, and duplication patterns in hundreds of cancer genomes. The most evident pattern that we discovered is that when one member of an SL interaction gene pair is lost, the other gene tends not to be lost, i.e. the absence of co-loss. This observation is in line with expectation, because the loss of an SL interacting pair will be lethal to the cancer cell. SL interactions are also reflected in gene expression profiles, such as an under representation of cases where the genes in an SL pair are both under expressed, and an over representation of cases where one gene of an SL pair is under expressed, while the other one is over expressed. We integrated the various previously unknown cancer genome patterns and the gene expression patterns into a computational model to identify SL pairs. This simple, genome-wide model achieves a high prediction power (AUC = 0.75) for known genetic interactions. It allows us to present for the first time a comprehensive genome-wide list of SL interactions with a high estimated prediction precision, covering up to 591,000 gene pairs. This unique list can potentially be used in various application areas ranging from biotechnology to medical genetics.

Somatic mutations favorable to patient survival are predominant in ovarian carcinomas

Somatic mutation accumulation is a major cause of abnormal cell growth. However, some mutations in cancer cells may be deleterious to the survival and proliferation of the cancer cells, thus offering a protective effect to the patients. We investigated this hypothesis via a unique analysis of the clinical and somatic mutation datasets of ovarian carcinomas published by the Cancer Genome Atlas. We defined and screened 562 macro mutation signatures (MMSs) for their associations with the overall survival of 320 ovarian cancer patients. Each MMS measures the number of mutations present on the member genes (except for TP53) covered by a specific Gene Ontology (GO) term in each tumor. We found that somatic mutations favorable to the patient survival are predominant in ovarian carcinomas compared to those indicating poor clinical outcomes. Specially, we identified 19 (3) predictive MMSs that are, usually by a nonlinear dose-dependent effect, associated with good (poor) patient survival. The false discovery rate for the 19 "positive" predictors is at the level of 0.15. The GO terms corresponding to these MMSs include "lysosomal membrane" and "response to hypoxia", each of which is relevant to the progression and therapy of cancer. Using these MMSs as features, we established a classification tree model which can effectively partition the training samples into three prognosis groups regarding the survival time. We validated this model on an independent dataset of the same disease (Log-rank p-value < 2.3 × 10(-4)) and a dataset of breast cancer (Log-rank p-value < 9.3 × 10(-3)). We compared the GO terms corresponding to these MMSs and those enriched with expression-based predictive genes. The analysis showed that the GO term pairs with large similarity are mainly pertinent to the proteins located on the cell organelles responsible for material transport and waste disposal, suggesting the crucial role of these proteins in cancer mortality.

A pair of homoeolog ClpP5 genes underlies a virescent yellow-like mutant and its modifier in maize

Gene-background interaction is a commonly observed phenomenon in many species, but the molecular mechanisms of such an interaction is less well understood. Here we report the cloning of a maize mutant gene and its modifier. A recessive mutant with a virescent yellow-like (vyl) phenotype was identified in an ethyl methanesulfonate-mutagenized population derived from the maize inbred line B73. Homozygous mutant maize plants exhibited a yellow leaf phenotype after emergence but gradually recovered and became indistinguishable from wild-type plants after approximately 2 weeks. Taking the positional cloning approach, the Chr.9_ClpP5 gene, one of the proteolytic subunits of the chloroplast Clp protease complex, was identified and validated as the candidate gene for vyl. When introgressed by backcross into the maize inbred line PH09B, the mutant phenotype of vyl lasted much longer in the greenhouse and was lethal in the field, implying the presence of a modifier(s) for vyl. A major modifier locus was identified on chromosome 1, and a paralogous ClpP5 gene was isolated and confirmed as the candidate for the vyl-modifier. Expression of Chr.1_ClpP5 is induced significantly in B73 by the vyl mutation, while the expression of Chr.1_ClpP5 in PH09B is not responsive to the vyl mutation. Moreover, expression and sequence analysis suggests that the PH09B Chr.1_ClpP5 allele is functionally weaker than the B73 allele. We propose that functional redundancy between duplicated paralogous genes is the molecular mechanism for the interaction between vyl and its modifier.

Genomic evidence for adaptation by gene duplication

Gene duplication is widely believed to facilitate adaptation, but unambiguous evidence for this hypothesis has been found in only a small number of cases. Although gene duplication may increase the fitness of the involved organisms by doubling gene dosage or neofunctionalization, it may also result in a simple division of ancestral functions into daughter genes, which need not promote adaptation. Hence, the general validity of the adaptation by gene duplication hypothesis remains uncertain. Indeed, a genome-scale experiment found similar fitness effects of deleting pairs of duplicate genes and deleting individual singleton genes from the yeast genome, leading to the conclusion that duplication rarely results in adaptation. Here we contend that the above comparison is unfair because of a known duplication bias among genes with different fitness contributions. To rectify this problem, we compare homologous genes from the budding yeast Saccharomyces cerevisiae and the fission yeast Schizosaccharomyces pombe. We discover that simultaneously deleting a duplicate gene pair in S. cerevisiae reduces fitness significantly more than deleting their singleton counterpart in S. pombe, revealing post-duplication adaptation. The duplicates-singleton difference in fitness effect is not attributable to a potential increase in gene dose after duplication, suggesting that the adaptation is owing to neofunctionalization, which we find to be explicable by acquisitions of binary protein-protein interactions rather than gene expression changes. These results provide genomic evidence for the role of gene duplication in organismal adaptation and are important for understanding the genetic mechanisms of evolutionary innovation.