Exposing the fitness contribution of duplicated genes

Frustration can Limit the Adaptation of Promiscuous Enzymes Through Gene Duplication and Specialisation

Virtually all enzymes catalyse more than one reaction, a phenomenon known as enzyme promiscuity. It is unclear whether promiscuous enzymes are more often generalists that catalyse multiple reactions at similar rates or specialists that catalyse one reaction much more efficiently than other reactions. In addition, the factors that shape whether an enzyme evolves to be a generalist or a specialist are poorly understood. To address these questions, we follow a three-pronged approach. First, we examine the distribution of promiscuity in empirical enzymes reported in the BRENDA database. We find that the promiscuity distribution of empirical enzymes is bimodal. In other words, a large fraction of promiscuous enzymes are either generalists or specialists, with few intermediates. Second, we demonstrate that enzyme biophysics is not sufficient to explain this bimodal distribution. Third, we devise a constraint-based model of promiscuous enzymes undergoing duplication and facing selection pressures favouring subfunctionalization. The model posits the existence of constraints between the catalytic efficiencies of an enzyme for different reactions and is inspired by empirical case studies. The promiscuity distribution predicted by our constraint-based model is consistent with the empirical bimodal distribution. Our results suggest that subfunctionalization is possible and beneficial only in certain enzymes. Furthermore, the model predicts that conflicting constraints and selection pressures can cause promiscuous enzymes to enter a 'frustrated' state, in which competing interactions limit the specialisation of enzymes. We find that frustration can be both a driver and an inhibitor of enzyme evolution by duplication and subfunctionalization. In addition, our model predicts that frustration becomes more likely as enzymes catalyse more reactions, implying that natural selection may prefer catalytically simple enzymes. In sum, our results suggest that frustration may play an important role in enzyme evolution.

RNA cis-regulators are important for Streptococcus pneumoniae in vivo success

Bacteria have evolved complex transcriptional regulatory networks, as well as many diverse regulatory strategies at the RNA level, to enable more efficient use of metabolic resources and a rapid response to changing conditions. However, most RNA-based regulatory mechanisms are not well conserved across different bacterial species despite controlling genes important for virulence or essential biosynthetic processes. Here, we characterize the activity of, and assess the fitness benefit conferred by, twelve cis-acting regulatory RNAs (including several riboswitches and a T-box), in the opportunistic pathogen Streptococcus pneumoniae TIGR4. By evaluating native locus mutants of each regulator that result in constitutively active or repressed expression, we establish that growth defects in planktonic culture are associated with constitutive repression of gene expression, while constitutive activation of gene expression is rarely deleterious. In contrast, in mouse nasal carriage and pneumonia models, strains with either constitutively active and repressed gene expression are significantly less fit than matched control strains. Furthermore, two RNA-regulated pathways, FMN synthesis/transport and pyrimidine synthesis/transport display exceptional sensitivity to mis-regulation or constitutive gene repression in both planktonic culture and in vivo environments. Thus, despite lack of obvious phenotypes associated with constitutive gene expression in vitro, the fitness benefit conferred on bacteria via fine-tuned metabolic regulation through cis-acting regulatory RNAs is substantial in vivo, and therefore easily sufficient to drive the evolution and maintenance of diverse RNA regulatory mechanisms.

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.

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.

Fit-Seq2.0: An Improved Software for High-Throughput Fitness Measurements Using Pooled Competition Assays

The fitness of a genotype is defined as its lifetime reproductive success, with fitness itself being a composite trait likely dependent on many underlying phenotypes. Measuring fitness is important for understanding how alteration of different cellular components affects a cell's ability to reproduce. Here, we describe an improved approach, implemented in Python, for estimating fitness in high throughput via pooled competition assays.

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.

POInTbrowse: orthology prediction and synteny exploration for paleopolyploid genomes

We describe POInT_browse, a web portal that gives access to the orthology inferences made for polyploid genomes with POInT, the Polyploidy Orthology Inference Tool. Ancient, or paleo-, polyploidy events are widely distributed across the eukaryotic phylogeny, and the combination of duplicated and lost duplicated genes that these polyploidies produce can confound the identification of orthologous genes between genomes. POInT uses conserved synteny and phylogenetic models to infer orthologous genes between genomes with a shared polyploidy. It also gives confidence estimates for those orthology inferences. POInT_browse gives both graphical and query-based access to these inferences from 12 different polyploidy events, allowing users to visualize genomic regions produced by polyploidies and perform batch queries for each polyploidy event, downloading genes trees and coding sequences for orthologous genes meeting user-specified criteria. POInT_browse and the associated data are online at https://wgd.statgen.ncsu.edu .

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 3,500 instances where only one of two paralogous proteins undergoes post-translational 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 post-translational 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.

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.

Using colony size to measure fitness in Saccharomyces cerevisiae

Competitive fitness assays in liquid culture have been a mainstay for characterizing experimental evolution of microbial populations. Growth of microbial strains has also been extensively characterized by colony size and could serve as a useful alternative if translated to per generation measurements of relative fitness. To examine fitness based on colony size, we established a relationship between cell number and colony size for strains of Saccharomyces cerevisiae robotically pinned onto solid agar plates in a high-density format. This was used to measure growth rates and estimate relative fitness differences between evolved strains and their ancestors. After controlling for edge effects through both normalization and agar-trimming, we found that colony size is a sensitive measure of fitness, capable of detecting 1% differences. While fitnesses determined from liquid and solid mediums were not equivalent, our results demonstrate that colony size provides a sensitive means of measuring fitness that is particularly well suited to measurements across many environments.

Synthetic negative genome screen of the GPN-loop GTPase NPA3 in Saccharomyces cerevisiae

The GPN-loop GTPase Npa3 is encoded by an essential gene in the yeast Saccharomyces cerevisiae. Npa3 plays a critical role in the assembly and nuclear accumulation of RNA polymerase II (RNAPII), a function that may explain its essentiality. Genetic interactions describe the extent to which a mutation in a particular gene affects a specific phenotype when co-occurring with an alteration in a second gene. Discovering synthetic negative genetic interactions has long been used as a tool to delineate the functional relatedness between pairs of genes participating in common or compensatory biological pathways. Previously, our group showed that nuclear targeting and transcriptional activity of RNAPII were unaffected in cells expressing exclusively a C-terminal truncated mutant version of Npa3 (npa3∆C) lacking the last 106 residues naturally absent from the single GPN protein in Archaea, but universally conserved in all Npa3 orthologs of eukaryotes. To gain insight into novel cellular functions for Npa3, we performed here a genome-wide Synthetic Genetic Array (SGA) study coupled to bulk fluorescence monitoring to identify negative genetic interactions of NPA3 by crossing an npa3∆C strain with a 4,389 nonessential gene-deletion collection. This genetic screen revealed previously unknown synthetic negative interactions between NPA3 and 15 genes. Our results revealed that the Npa3 C-terminal tail extension regulates the participation of this essential GTPase in previously unknown biological processes related to mitochondrial homeostasis and ribosome biogenesis.

Evolution of binding preferences among whole-genome duplicated transcription factors

Throughout evolution, new transcription factors (TFs) emerge by gene duplication, promoting growth and rewiring of transcriptional networks. How TF duplicates diverge was studied in a few cases only. To provide a genome-scale view, we considered the set of budding yeast TFs classified as whole-genome duplication (WGD)-retained paralogs (~35% of all specific TFs). Using high-resolution profiling, we find that ~60% of paralogs evolved differential binding preferences. We show that this divergence results primarily from variations outside the DNA-binding domains (DBDs), while DBD preferences remain largely conserved. Analysis of non-WGD orthologs revealed uneven splitting of ancestral preferences between duplicates, and the preferential acquiring of new targets by the least conserved paralog (biased neo/sub-functionalization). Interactions between paralogs were rare, and, when present, occurred through weak competition for DNA-binding or dependency between dimer-forming paralogs. We discuss the implications of our findings for the evolutionary design of transcriptional networks.

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.

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.

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.

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.

Degree of Functional Divergence in Duplicates Is Associated with Distinct Roles in Plant Evolution

Gene duplication is a major mechanism to create new genes. After gene duplication, some duplicated genes undergo functionalization, whereas others largely maintain redundant functions. Duplicated genes comprise various degrees of functional diversification in plants. However, the evolutionary fate of high and low diversified duplicates is unclear at genomic scale. To infer high and low diversified duplicates in Arabidopsis thaliana genome, we generated a prediction method for predicting whether a pair of duplicate genes was subjected to high or low diversification based on the phenotypes of knock-out mutants. Among 4,017 pairs of recently duplicated A. thaliana genes, 1,052 and 600 are high and low diversified duplicate pairs, respectively. The predictions were validated based on the phenotypes of generated knock-down transgenic plants. We determined that the high diversified duplicates resulting from tandem duplications tend to have lineage-specific functions, whereas the low diversified duplicates produced by whole-genome duplications are related to essential signaling pathways. To assess the evolutionary impact of high and low diversified duplicates in closely related species, we compared the retention rates and selection pressures on the orthologs of A. thaliana duplicates in two closely related species. Interestingly, high diversified duplicates resulting from tandem duplications tend to be retained in multiple lineages under positive selection. Low diversified duplicates by whole-genome duplications tend to be retained in multiple lineages under purifying selection. Taken together, the functional diversities determined by different duplication mechanisms had distinct effects on plant evolution.

Genome-Wide Identification, Structure Characterization, and Expression Pattern Profiling of the Aquaporin Gene Family in Betula pendula

Aquaporin water channels (AQPs) constitute a large family of transmembrane proteins present throughout all kingdoms of life. They play key roles in the flux of water and many solutes across the membranes. The AQP diversity, protein features, and biological functions of silver birch are still unknown. A genome analysis of Betula pendula identified 33 putative genes encoding full-length AQP sequences (BpeAQPs). They are grouped into five subfamilies, representing ten plasma membrane intrinsic proteins (PIPs), eight tonoplast intrinsic proteins (TIPs), eight NOD26-like intrinsic proteins (NIPs), four X intrinsic proteins (XIPs), and three small basic intrinsic proteins (SIPs). The BpeAQP gene structure is conserved within each subfamily, with exon numbers ranging from one to five. The predictions of the aromatic/arginine selectivity filter (ar/R), Froger's positions, specificity-determining positions, and 2D and 3D biochemical properties indicate noticeable transport specificities to various non-aqueous substrates between members and/or subfamilies. Nevertheless, overall, the BpePIPs display mostly hydrophilic ar/R selective filter and lining-pore residues, whereas the BpeTIP, BpeNIP, BpeSIP, and BpeXIP subfamilies mostly contain hydrophobic permeation signatures. Transcriptional expression analyses indicate that 23 BpeAQP genes are transcribed, including five organ-related expressions. Surprisingly, no significant transcriptional expression is monitored in leaves in response to cold stress (6 °C), although interesting trends can be distinguished and will be discussed, notably in relation to the plasticity of this pioneer species, B. pendula. The current study presents the first detailed genome-wide analysis of the AQP gene family in a Betulaceae species, and our results lay a foundation for a better understanding of the specific functions of the BpeAQP genes in the responses of the silver birch trees to cold stress.

Expression attenuation as a mechanism of robustness against gene duplication

Many studies have focused on the mechanisms of long-term retention of gene duplicates, such as the gain of functions or reciprocal losses. However, such changes are more likely to occur if the duplicates are maintained for a long period. This time span will be short if duplication is immediately deleterious. We measured the distribution of fitness effects of gene duplication for 899 genes in budding yeast. We find that gene duplication is more likely to be deleterious than beneficial. However, contrary to previous models, in general, gene duplication does not affect fitness by altering the organization of protein complexes. We show that expression attenuation may protect complexes from the effects of gene duplication.

Gene duplication is ubiquitous and a major driver of phenotypic diversity across the tree of life, but its immediate consequences are not fully understood. Deleterious effects would decrease the probability of retention of duplicates and prevent their contribution to long-term evolution. One possible detrimental effect of duplication is the perturbation of the stoichiometry of protein complexes. Here, we measured the fitness effects of the duplication of 899 essential genes in the budding yeast using high-resolution competition assays. At least 10% of genes caused a fitness disadvantage when duplicated. Intriguingly, the duplication of most protein complex subunits had small to nondetectable effects on fitness, with few exceptions. We selected four complexes with subunits that had an impact on fitness when duplicated and measured the impact of individual gene duplications on their protein–protein interactions. We found that very few duplications affect both fitness and interactions. Furthermore, large complexes such as the 26S proteasome are protected from gene duplication by attenuation of protein abundance. Regulatory mechanisms that maintain the stoichiometric balance of protein complexes may protect from the immediate effects of gene duplication. Our results show that a better understanding of protein regulation and assembly in complexes is required for the refinement of current models of gene duplication.

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.

The Pathophysiology of Rett Syndrome With a Focus on Breathing Dysfunctions

Rett syndrome (RTT), an X-chromosome-linked neurological disorder, is characterized by serious pathophysiology, including breathing and feeding dysfunctions, and alteration of cardiorespiratory coupling, a consequence of multiple interrelated disturbances in the genetic and homeostatic regulation of central and peripheral neuronal networks, redox state, and control of inflammation. Characteristic breath-holds, obstructive sleep apnea, and aerophagia result in intermittent hypoxia, which, combined with mitochondrial dysfunction, causes oxidative stress-an important driver of the clinical presentation of RTT.

Multiplex enCas12a screens detect functional buffering among paralogs otherwise masked in monogenic Cas9 knockout screens

BACKGROUND

Pooled library CRISPR/Cas9 knockout screening across hundreds of cell lines has identified genes whose disruption leads to fitness defects, a critical step in identifying candidate cancer targets. However, the number of essential genes detected from these monogenic knockout screens is low compared to the number of constitutively expressed genes in a cell.

RESULTS

Through a systematic analysis of screen data in cancer cell lines generated by the Cancer Dependency Map, we observe that half of all constitutively expressed genes are never detected in any CRISPR screen and that these never-essentials are highly enriched for paralogs. We investigated functional buffering among approximately 400 candidate paralog pairs using CRISPR/enCas12a dual-gene knockout screening in three cell lines. We observe 24 synthetic lethal paralog pairs that have escaped detection by monogenic knockout screens at stringent thresholds. Nineteen of 24 (79%) synthetic lethal interactions are present in at least two out of three cell lines and 14 of 24 (58%) are present in all three cell lines tested, including alternate subunits of stable protein complexes as well as functionally redundant enzymes.

CONCLUSIONS

Together, these observations strongly suggest that functionally redundant paralogs represent a targetable set of genetic dependencies that are systematically under-represented among cell-essential genes in monogenic CRISPR-based loss of function screens.

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.

Phylogeny, evolution, and potential ecological relationship of cytochrome CYP52 enzymes in Saccharomycetales yeasts

Cytochrome P450s from the CYP52 family participate in the assimilation of alkanes and fatty acids in fungi. In this work, the evolutionary history of a set of orthologous and paralogous CYP52 proteins from Saccharomycetales yeasts was inferred. Further, the phenotypic assimilation profiles were related with the distribution of cytochrome CYP52 members among species. The maximum likelihood phylogeny of CYP52 inferred proteins reveled a frequent ancient and modern duplication and loss events that generated orthologous and paralogous groups. Phylogeny and assimilation profiles of alkanes and fatty acids showed a family expansion in yeast isolated from hydrophobic-rich environments. Docking analysis of deduced ancient CYP52 proteins suggests that the most ancient function was the oxidation of C4-C11 alkanes, while the oxidation of >10 carbon alkanes and fatty acids is a derived character. The ancient CYP52 paralogs displayed partial specialization and promiscuous interaction with hydrophobic substrates. Additionally, functional optimization was not evident. Changes in the interaction of ancient CYP52 with different alkanes and fatty acids could be associated with modifications in spatial orientations of the amino acid residues that comprise the active site. The extended family of CYP52 proteins is likely evolving toward functional specialization, and certain redundancy for substrates is being maintained.

The lasting after-effects of an ancient polyploidy on the genomes of teleosts

The ancestor of most teleost fishes underwent a whole-genome duplication event three hundred million years ago. Despite its antiquity, the effects of this event are evident both in the structure of teleost genomes and in how the surviving duplicated genes still operate to drive form and function. I inferred a set of shared syntenic regions that survive from the teleost genome duplication (TGD) using eight teleost genomes and the outgroup gar genome (which lacks the TGD). I then phylogenetically modeled the TGD's resolution via shared and independent gene losses and applied a new simulation-based statistical test for the presence of bias toward the preservation of genes from one parental subgenome. On the basis of that test, I argue that the TGD was likely an allopolyploidy. I find that duplicate genes surviving from this duplication in zebrafish are less likely to function in early embryo development than are genes that have returned to single copy at some point in this species' history. The tissues these ohnologs are expressed in, as well as their biological functions, lend support to recent suggestions that the TGD was the source of a morphological innovation in the structure of the teleost retina. Surviving duplicates also appear less likely to be essential than singletons, despite the fact that their single-copy orthologs in mouse are no less essential than other genes.

Similarities in biological processes can be used to bridge ecology and molecular biology

Much of the research in biology aims to understand the origin of diversity. Naturally, ecological diversity was the first object of study, but we now have the necessary tools to probe diversity at molecular scales. The inherent differences in how we study diversity at different scales caused the disciplines of biology to be organized around these levels, from molecular biology to ecology. Here, we illustrate that there are key properties of each scale that emerge from the interactions of simpler components and that these properties are often shared across different levels of organization. This means that ideas from one level of organization can be an inspiration for novel hypotheses to study phenomena at another level. We illustrate this concept with examples of events at the molecular level that have analogs at the organismal or ecological level and vice versa. Through these examples, we illustrate that biological processes at different organization levels are governed by general rules. The study of the same phenomena at different scales could enrich our work through a multidisciplinary approach, which should be a staple in the training of future scientists.

Gene family information facilitates variant interpretation and identification of disease-associated genes in neurodevelopmental disorders

BACKGROUND

Classifying pathogenicity of missense variants represents a major challenge in clinical practice during the diagnoses of rare and genetic heterogeneous neurodevelopmental disorders (NDDs). While orthologous gene conservation is commonly employed in variant annotation, approximately 80% of known disease-associated genes belong to gene families. The use of gene family information for disease gene discovery and variant interpretation has not yet been investigated on a genome-wide scale. We empirically evaluate whether paralog-conserved or non-conserved sites in human gene families are important in NDDs.

METHODS

Gene family information was collected from Ensembl. Paralog-conserved sites were defined based on paralog sequence alignments; 10,068 NDD patients and 2078 controls were statistically evaluated for de novo variant burden in gene families.

RESULTS

We demonstrate that disease-associated missense variants are enriched at paralog-conserved sites across all disease groups and inheritance models tested. We developed a gene family de novo enrichment framework that identified 43 exome-wide enriched gene families including 98 de novo variant carrying genes in NDD patients of which 28 represent novel candidate genes for NDD which are brain expressed and under evolutionary constraint.

CONCLUSION

This study represents the first method to incorporate gene family information into a statistical framework to interpret variant data for NDDs and to discover new NDD-associated genes.

Global Genetic Networks and the Genotype-to-Phenotype Relationship

Genetic interactions identify combinations of genetic variants that impinge on phenotype. With whole-genome sequence information available for thousands of individuals within a species, a major outstanding issue concerns the interpretation of allelic combinations of genes underlying inherited traits. In this Review, we discuss how large-scale analyses in model systems have illuminated the general principles and phenotypic impact of genetic interactions. We focus on studies in budding yeast, including the mapping of a global genetic network. We emphasize how information gained from work in yeast translates to other systems, and how a global genetic network not only annotates gene function but also provides new insights into the genotype-to-phenotype relationship.

A Single, Shared Triploidy in Three Species of Parasitic Nematodes

The root-knot nematodes of the genus Meloidogyne are important and damaging parasites capable of infecting most flowering plants. Within this genus, several species of the Meloidogyne incognita group show evidence of paleopolyploidy in their genomes. We used our software tool POInT, the Polyploidy Orthology Inference Tool, to phylogenetically model the gene losses that followed that polyploidy. These models, and simulations based on them, show that three of these species (M. incognita, M. arenaria and M. javanica) descend from a single common hybridization event that yielded triplicated genomes with three distinguishable subgenomes. While one of the three subgenomes shows elevated gene loss rates relative to the other two, this subgenome does not show elevated sequence divergence. In all three species, ancestral loci where two of the three gene copies have been lost are less likely to have orthologs in Caenorhabditis elegans that are lethal when knocked down than are ancestral loci with surviving duplicate copies.

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.

Characterisation of a rice vacuolar invertase isoform, OsINV2, for growth and yield-related traits

OsINV2, a rice vacuolar invertase isoform, was assessed for its functional roles in plant growth and development with key focus on its agronomic traits such as grain weight, grain filling percentage, grain number and dry weights at various stages until harvest. Lack of differences between the wild-type and the mutants with respect to any of the aforementioned traits tested revealed a possibility of functional compensation of OsINV2 in the mutants conceivably by its isoform. This was confirmed by OsINV2 promoter::GUS studies, where its spatial and temporal expression in the panicle elongation stages showed that although OsINV2 expression was observed from the stage with young panicles ~1 cm in length to the flag leaf stage, significant differences with respect to panicle and spikelet phenotypes between the wild-type and the mutant were not present. However, complement lines displaying an overexpression phenotype of OsINV2 possessed a higher stem non-structural carbohydrate content under both monoculm and normal tillering conditions. A trade-off between the spikelet number and grain weight in the complement lines grown under monoculm conditions was also observed, pointing towards the necessity of OsINV2 regulation for grain yield-related traits.

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.

Redundancy, Feedback, and Robustness in the Arabidopsis thaliana BZR/BEH Gene Family

Organismal development is remarkably robust, tolerating stochastic errors to produce consistent, so-called canalized adult phenotypes. The mechanistic underpinnings of developmental robustness are poorly understood, but recent studies implicate certain features of genetic networks such as functional redundancy, connectivity, and feedback. Here, we examine the BZR/BEH gene family, whose function contributes to embryonic stem development in the plant Arabidopsis thaliana, to test current assumptions on functional redundancy and trait robustness. Our analyses of BZR/BEH gene mutants and mutant combinations revealed that functional redundancy among these gene family members is not necessary for trait robustness. Connectivity is another commonly cited determinant of robustness; however, we found no correlation between connectivity among gene family members or their connectivity with other transcription factors and effects on developmental robustness. Instead, our data suggest that BEH4, the earliest diverged family member, modulates developmental robustness. We present evidence indicating that regulatory cross-talk among gene family members is integrated by BEH4 to promote wild-type levels of developmental robustness. Further, the chaperone HSP90, a known determinant of developmental robustness, appears to act via BEH4 in maintaining robustness of embryonic stem length. In summary, we demonstrate that even among closely related transcription factors, trait robustness can arise through the activity of a single gene family member, challenging common assumptions about the molecular underpinnings of robustness.

The DEAD-Box RNA Helicases of Bacillus subtilis as a Model to Evaluate Genetic Compensation Among Duplicate Genes

The presence of duplicated genes in organisms is well documented. There is increasing interest in understanding how these genes subfunctionalize and whether functional overlap can explain the fact that some of these genes are dispensable. Bacillus subtilis possesses four DEAD-box RNA helicases (DBRH) genes, cshA, cshB, deaD/yxiN, and yfmL that make a good case to study to what extent they can complement each other despite their subfunctionalization. They possess the highly conserved N-terminal catalytic domain core common to RNA helicases, but different carboxy-terminal ends. All four genes have been shown to have independent functions although all participate in rRNA assembly. None of the B. subtilis DBRH is essential for growth at 37°C, and all single deletion mutants exhibit defective growth at 18°C except for ΔdeaD/yxiN. Evaluation of double mutants did not reveal negative epistasis, suggesting that they do not have overlapping functions. The absence of any one gene distorts the expression pattern of the others, but not in a specific pattern suggestive of compensation. Overexpression of these paralogous genes in the different mutant backgrounds did not result in cross-complementation, further confirming their lack of buffering capability. Since no complementation could be observed among full sized proteins, we evaluated to what extent the superfamily 2 (SF2) helicase core of the smallest DBRH, YfmL, could be functional when hooked to each of the C-terminal end of CshA, CshB, and DeaD/YxiN. None of the different chimeras complemented the different mutants, and instead, all chimeras inhibited the growth of the ΔyfmL mutant, and other combinations were also deleterious. Our findings suggest that the long time divergence between DEAD-box RNA helicase genes has resulted in specialized activities in RNA metabolism and shows that these duplicated genes cannot buffer one another.

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.

Mapping DNA damage-dependent genetic interactions in yeast via party mating and barcode fusion genetics

Condition‐dependent genetic interactions can reveal functional relationships between genes that are not evident under standard culture conditions. State‐of‐the‐art yeast genetic interaction mapping, which relies on robotic manipulation of arrays of double‐mutant strains, does not scale readily to multi‐condition studies. Here, we describe barcode fusion genetics to map genetic interactions (BFG‐GI), by which double‐mutant strains generated via en masse “party” mating can also be monitored en masse for growth to detect genetic interactions. By using site‐specific recombination to fuse two DNA barcodes, each representing a specific gene deletion, BFG‐GI enables multiplexed quantitative tracking of double mutants via next‐generation sequencing. We applied BFG‐GI to a matrix of DNA repair genes under nine different conditions, including methyl methanesulfonate (MMS), 4‐nitroquinoline 1‐oxide (4NQO), bleomycin, zeocin, and three other DNA‐damaging environments. BFG‐GI recapitulated known genetic interactions and yielded new condition‐dependent genetic interactions. We validated and further explored a subnetwork of condition‐dependent genetic interactions involving MAG1,SLX4, and genes encoding the Shu complex, and inferred that loss of the Shu complex leads to an increase in the activation of the checkpoint protein kinase Rad53.

Whole-Genome Analysis of Three Yeast Strains Used for Production of Sherry-Like Wines Revealed Genetic Traits Specific to Flor Yeasts

Flor yeast strains represent a specialized group of Saccharomyces cerevisiae yeasts used for biological wine aging. We have sequenced the genomes of three flor strains originated from different geographic regions and used for production of sherry-like wines in Russia. According to the obtained phylogeny of 118 yeast strains, flor strains form very tight cluster adjacent to the main wine clade. SNP analysis versus available genomes of wine and flor strains revealed 2,270 genetic variants in 1,337 loci specific to flor strains. Gene ontology analysis in combination with gene content evaluation revealed a complex landscape of possibly adaptive genetic changes in flor yeast, related to genes associated with cell morphology, mitotic cell cycle, ion homeostasis, DNA repair, carbohydrate metabolism, lipid metabolism, and cell wall biogenesis. Pangenomic analysis discovered the presence of several well-known "non-reference" loci of potential industrial importance. Events of gene loss included deletions of asparaginase genes, maltose utilization locus, and FRE-FIT locus involved in iron transport. The latter in combination with a flor-yeast-specific mutation in the Aft1 transcription factor gene is likely to be responsible for the discovered phenotype of increased iron sensitivity and improved iron uptake of analyzed strains. Expansion of the coding region of the FLO11 flocullin gene and alteration of the balance between members of the FLO gene family are likely to positively affect the well-known propensity of flor strains for velum formation. Our study provides new insights in the nature of genetic variation in flor yeast strains and demonstrates that different adaptive properties of flor yeast strains could have evolved through different mechanisms of genetic variation.

New oenological practice to promote non-Saccharomyces species of interest: saturating grape juice with carbon dioxide

Non-Saccharomyces yeast species, naturally found in grape must, may impact wine quality positively or negatively. In this study, a mixture of five non-Saccharomyces species (Torulaspora delbrueckii, Metschnikowia spp., Starmerella bacillaris (formerly called Candida zemplinina), Hanseniaspora uvarum, Pichia kluyveri), mimicking the composition of the natural non-Saccharomyces community found in grape must, was used for alcoholic fermentation. The impact of CO₂ saturation of the grape juice was studied first on this mixture alone, and then in the presence of Saccharomyces cerevisiae. Two isogenic strains of this species were used: the first with a short and the second a long fermentation lag phase. This study demonstrated that saturating grape juice with CO₂ had interesting potential as an oenological technique, inhibiting undesirable species (S. bacillaris and H. uvarum) and stimulating non-Saccharomyces of interest (T. delbrueckii and P. kluyveri). This stimulating effect was particularly marked when CO₂ saturation was associated with the presence of S. cerevisiae with long fermentation lag phase. The direct consequence of this association was an enhancement of 3-SH levels in the resulting wine.

Ploidy Variation in Fungi: Polyploidy, Aneuploidy, and Genome Evolution

The ability of an organism to replicate and segregate its genome with high fidelity is vital to its survival and for the production of future generations. Errors in either of these steps (replication or segregation) can lead to a change in ploidy or chromosome number. While these drastic genome changes can be detrimental to the organism, resulting in decreased fitness, they can also provide increased fitness during periods of stress. A change in ploidy or chromosome number can fundamentally change how a cell senses and responds to its environment. Here, we discuss current ideas in fungal biology that illuminate how eukaryotic genome size variation can impact the organism at a cellular and evolutionary level. One of the most fascinating observations from the past 2 decades of research is that some fungi have evolved the ability to tolerate large genome size changes and generate vast genomic heterogeneity without undergoing canonical meiosis.

Intrinsic adaptive value and early fate of gene duplication revealed by a bottom-up approach

The population genetic mechanisms governing the preservation of gene duplicates, especially in the critical very initial phase, have remained largely unknown. Here, we demonstrate that gene duplication confers per se a weak selective advantage in scenarios of fitness trade-offs. Through a precise quantitative description of a model system, we show that a second gene copy serves to reduce gene expression inaccuracies derived from pervasive molecular noise and suboptimal gene regulation. We then reveal that such an accuracy in the phenotype yields a selective advantage in the order of 0.1% on average, which would allow the positive selection of gene duplication in populations with moderate/large sizes. This advantage is greater at higher noise levels and intermediate concentrations of the environmental molecule, when fitness trade-offs become more evident. Moreover, we discuss how the genome rearrangement rates greatly condition the eventual fixation of duplicates. Overall, our theoretical results highlight an original adaptive value for cells carrying new-born duplicates, broadly analyze the selective conditions that determine their early fates in different organisms, and reconcile population genetics with evolution by gene duplication.

Protein Moonlighting Revealed by Noncatalytic Phenotypes of Yeast Enzymes

A single gene can partake in several biological processes, and therefore gene deletions can lead to different-sometimes unexpected-phenotypes. However, it is not always clear whether such pleiotropy reflects the loss of a unique molecular activity involved in different processes or the loss of a multifunctional protein. Here, using Saccharomyces cerevisiae metabolism as a model, we systematically test the null hypothesis that enzyme phenotypes depend on a single annotated molecular function, namely their catalysis. We screened a set of carefully selected genes by quantifying the contribution of catalysis to gene deletion phenotypes under different environmental conditions. While most phenotypes were explained by loss of catalysis, slow growth was readily rescued by a catalytically inactive protein in about one-third of the enzymes tested. Such noncatalytic phenotypes were frequent in the Alt1 and Bat2 transaminases and in the isoleucine/valine biosynthetic enzymes Ilv1 and Ilv2, suggesting novel "moonlighting" activities in these proteins. Furthermore, differential genetic interaction profiles of gene deletion and catalytic mutants indicated that ILV1 is functionally associated with regulatory processes, specifically to chromatin modification. Our systematic study shows that gene loss phenotypes and their genetic interactions are frequently not driven by the loss of an annotated catalytic function, underscoring the moonlighting nature of cellular metabolism.

Gene duplication can impart fragility, not robustness, in the yeast protein interaction network

The maintenance of duplicated genes is thought to protect cells from genetic perturbations, but the molecular basis of this robustness is largely unknown. By measuring the interaction of yeast proteins with their partners in wild-type cells and in cells lacking a paralog, we found that 22 out of 56 paralog pairs compensate for the lost interactions. An equivalent number of pairs exhibit the opposite behavior and require each other's presence for maintaining their interactions. These dependent paralogs generally interact physically, regulate each other's abundance, and derive from ancestral self-interacting proteins. This reveals that gene duplication may actually increase mutational fragility instead of robustness in a large number of cases.