Essential gene disruptions reveal complex relationships between phenotypic robustness,pleiotropy, and fitness

The concept of robustness in biology has gained much attention recently, but a mechanistic understanding of how genetic networks regulate phenotypic variation has remained elusive. One approach to understand the genetic architecture of variability has been to analyze dispensable gene deletions in model organisms; however, the most important genes cannot be deleted. Here, we have utilized two systems in yeast whereby essential genes have been altered to reduce expression. Using high-throughput microscopy and image analysis, we have characterized a large number of morphological phenotypes, and their associated variation, for the majority of essential genes in yeast. Our results indicate that phenotypic robustness is more highly dependent upon the expression of essential genes than on the presence of dispensable genes. Morphological robustness appears to be a general property of a genotype that is closely related to pleiotropy. While the fitness profile across a range of expression levels is idiosyncratic to each gene, the global pattern indicates that there is a window in which phenotypic variation can be released before fitness effects are observable.     

Large-scale analysis of gene function in Caenorhabditis elegans by high-throughput RNAi

Genome-wide analysis of gene function is essential for the post-genome era, and development of efficient and economical technology suitable for it has been in demand. Here we report a large-scale inactivation of the expressed genes in the nematode Caenorhabditis elegans. For this purpose, we have established a high-throughput "RNAi-by-soaking" methodology by modifying the conventional RNAi method [1, 2]. A set of tag-sequenced, nonredundant cDNAs corresponding to approximately 10,000 genes [3] (representing half of the predicted genes [4]) was used for the systematic RNAi analysis. We have processed approximately 2500 genes to date. In development, 27% of them showed detectable phenotypes, such as embryonic lethality, post-embryonic lethality, sterility, and morphological abnormality. Of these, we analyzed the phenotypes of F1 sterility in detail, and we have identified 24 genes that might play important roles in germline development. Combined with the ongoing analysis of expression patterns of these cDNAs [3, 5], the functional information obtained in this work will provide a starting point for the further analysis of each gene. Another finding from this screening is that the incidence of essential genes is significantly lower in the X chromosome than in the autosomes.     

Towards an Informative Mutant Phenotype for Every Bacterial Gene

Mutant phenotypes provide strong clues to the functions of the underlying genes and could allow annotation of the millions of sequenced yet uncharacterized bacterial genes. However, it is not known how many genes have a phenotype under laboratory conditions, how many phenotypes are biologically interpretable for predicting gene function, and what experimental conditions are optimal to maximize the number of genes with a phenotype. To address these issues, we measured the mutant fitness of 1,586 genes of the ethanol-producing bacterium Zymomonas mobilis ZM4 across 492 diverse experiments and found statistically significant phenotypes for 89% of all assayed genes. Thus, in Z. mobilis, most genes have a functional consequence under laboratory conditions. We demonstrate that 41% of Z. mobilis genes have both a strong phenotype and a similar fitness pattern (cofitness) to another gene, and are therefore good candidates for functional annotation using mutant fitness. Among 502 poorly characterized Z. mobilis genes, we identified a significant cofitness relationship for 174. For 57 of these genes without a specific functional annotation, we found additional evidence to support the biological significance of these gene-gene associations, and in 33 instances, we were able to predict specific physiological or biochemical roles for the poorly characterized genes. Last, we identified a set of 79 diverse mutant fitness experiments in Z. mobilis that are nearly as biologically informative as the entire set of 492 experiments. Therefore, our work provides a blueprint for the functional annotation of diverse bacteria using mutant fitness.     

Fitness effects of mutation: testing genetic redundancy in Arabidopsis thaliana

Screens of organisms with disruptive mutations in a single gene often fail to detect phenotypic consequences for the majority of mutants. One explanation for this phenomenon is that the presence of paralogous loci provides genetic redundancy. However, it is also possible that the assayed traits are affected by few loci, that effects could be subtle or that phenotypic effects are restricted to certain environments. We assayed a set of T‐DNA insertion mutant lines of Arabidopsis thaliana to determine the frequency with which mutation affected fitness‐related phenotypes. We found that between 8% and 42% of the assayed lines had altered fitness from the wild type. Furthermore, many of these lines exhibited fitness greater than the wild type. In a second experiment, we grew a subset of the lines in multiple environments and found whether a T‐DNA insert increased or decreased fitness traits depended on the assay environment. Overall, our evidence contradicts the hypothesis that genetic redundancy is a common phenomenon in A. thaliana for fitness traits. Evidence for redundancy from prior screens of knockout mutants may often be an artefact of the design of the phenotypic assays which have focused on less complex phenotypes than fitness and have used single environments. Finally, our study adds to evidence that beneficial mutations may represent a significant component of the mutational spectrum of A. thaliana.     

Metabolic network analysis revealed distinct routes of deletion effects between essential and non-essential genes

Interest in essential genes has arisen recently given their importance in antimicrobial drug development. Although knockouts of essential genes are commonly known to cause lethal phenotypes, there is insufficient understanding on the intermediate changes followed by genetic perturbation and to what extent essential genes correlate to other genes. Here, we characterized the gene knockout effects by using a list of affected genes, termed as 'damage lists'. These damage lists were identified through a refined cascading failure approach that was based on a previous topological flux balance analysis. Using an Escherichia coli metabolic network, we incorporated essentiality information into damage lists and revealed that the knockout of an essential gene mainly affects a large range of other essential genes whereas knockout of a non-essential gene only interrupts other non-essential genes. Also, genes sharing common damage lists tend to have the same essentiality. We extracted 72 core functional modules from the common damage lists of essential genes and demonstrated their ability to halt essential metabolites production. Overall, our network analysis revealed that essential and non-essential genes propagated their deletion effects via distinct routes, conferring mechanistic explanation to the observed lethality phenotypes of essential genes.     

Characteristics of Plant Essential Genes Allow for within- and between-Species Prediction of Lethal Mutant Phenotypes

Essential genes represent critical cellular components whose disruption results in lethality. Characteristics shared among essential genes have been uncovered in fungal and metazoan model systems. However, features associated with plant essential genes are largely unknown and the full set of essential genes remains to be discovered in any plant species. Here, we show that essential genes in Arabidopsis thaliana have distinct features useful for constructing within- and cross-species prediction models. Essential genes in A. thaliana are often single copy or derived from older duplications, highly and broadly expressed, slow evolving, and highly connected within molecular networks compared with genes with nonlethal mutant phenotypes. These gene features allowed the application of machine learning methods that predicted known lethal genes as well as an additional 1970 likely essential genes without documented phenotypes. Prediction models from A. thaliana could also be applied to predict Oryza sativa and Saccharomyces cerevisiae essential genes. Importantly, successful predictions drew upon many features, while any single feature was not sufficient. Our findings show that essential genes can be distinguished from genes with nonlethal phenotypes using features that are similar across kingdoms and indicate the possibility for translational application of our approach to species without extensive functional genomic and phenomic resources.     

Evolution in Saccharomyces cerevisiae: identification of mutations increasing fitness in laboratory populations

Since the publication of the complete sequence of the genome of Saccharomyces cerevisiae, a number of comprehensive investigations have been initiated to gain insight into cellular function. The focus of these studies has been to identify genes essential for survival in specific environments or those that when mutated cause gross phenotypic defects in growth. Here we describe Ty1-based mutational approaches designed to identify genes, which when mutated generate evolutionarily significant phenotypes causing small but positive increments on fitness. As expected, Ty1 mutations with a positive fitness effect were in the minority. However, mutations in two loci, one inactivating FAR3 and one upstream of CYR1, identified in evolving populations, were shown to have small but significantly positive fitness effects.     

Effects of juvenile hormone III on Reticulitermes flavipes: changes in hemolymph protein composition and gene expression

Termites express polyphenism during caste differentiation that is mostly undefined at the molecular level. Using the eastern subterranean termite, Reticulitermes flavipes Kollar, we wanted (1) to test juvenile hormone (JH) model assays for their ability to induce detectable molecular changes in worker termites and (2) to investigate hemolymph proteins and their corresponding genes during JH-induced soldier caste differentiation. Our results illustrate pronounced changes in two hemolymph proteins after JH treatment, as well as differences among several caste phenotypes. Significant increases in the expression of four genes encoding hemolymph proteins, including two vitellogenins and two hexamerins, were observed after JH exposure. These findings are the first to demonstrate such protein and gene expression changes during termite caste differentiation. These results also validate the utility of JH model assays for inducing detectable molecular changes in worker termites that have begun presoldier differentiation.     

Evidence-based annotation of gene function in Shewanella oneidensis MR-1 using genome-wide fitness profiling across 121 conditions

Most genes in bacteria are experimentally uncharacterized and cannot be annotated with a specific function. Given the great diversity of bacteria and the ease of genome sequencing, high-throughput approaches to identify gene function experimentally are needed. Here, we use pools of tagged transposon mutants in the metal-reducing bacterium Shewanella oneidensis MR-1 to probe the mutant fitness of 3,355 genes in 121 diverse conditions including different growth substrates, alternative electron acceptors, stresses, and motility. We find that 2,350 genes have a pattern of fitness that is significantly different from random and 1,230 of these genes (37% of our total assayed genes) have enough signal to show strong biological correlations. We find that genes in all functional categories have phenotypes, including hundreds of hypotheticals, and that potentially redundant genes (over 50% amino acid identity to another gene in the genome) are also likely to have distinct phenotypes. Using fitness patterns, we were able to propose specific molecular functions for 40 genes or operons that lacked specific annotations or had incomplete annotations. In one example, we demonstrate that the previously hypothetical gene SO_3749 encodes a functional acetylornithine deacetylase, thus filling a missing step in S. oneidensis metabolism. Additionally, we demonstrate that the orphan histidine kinase SO_2742 and orphan response regulator SO_2648 form a signal transduction pathway that activates expression of acetyl-CoA synthase and is required for S. oneidensis to grow on acetate as a carbon source. Lastly, we demonstrate that gene expression and mutant fitness are poorly correlated and that mutant fitness generates more confident predictions of gene function than does gene expression. The approach described here can be applied generally to create large-scale gene-phenotype maps for evidence-based annotation of gene function in prokaryotes.     

Diverse phenotypic and genetic responses to short‐term selection in evolving Escherichia coli populations

Beneficial mutations fuel adaptation by altering phenotypes that enhance the fit of organisms to their environment. However, the phenotypic effects of mutations often depend on ecological context, making the distribution of effects across multiple environments essential to understanding the true nature of beneficial mutations. Studies that address both the genetic basis and ecological consequences of adaptive mutations remain rare. Here, we characterize the direct and pleiotropic fitness effects of a collection of 21 first‐step beneficial mutants derived from naïve and adapted genotypes used in a long‐term experimental evolution of Escherichia coli. Whole‐genome sequencing was able to identify the majority of beneficial mutations. In contrast to previous studies, we find diverse fitness effects of mutations selected in a simple environment and few cases of genetic parallelism. The pleiotropic effects of these mutations were predominantly positive but some mutants were highly antagonistic in alternative environments. Further, the fitness effects of mutations derived from the adapted genotypes were dramatically reduced in nearly all environments. These findings suggest that many beneficial variants are accessible from a single point on the fitness landscape, and the fixation of alternative beneficial mutations may have dramatic consequences for niche breadth reduction via metabolic erosion.     

Higher duplicability of less important genes in yeast genomes

Gene duplication plays an important role in evolution because it is the primary source of new genes. Many recent studies showed that gene duplicability varies considerably among genes. Several considerations led us to hypothesize that less important genes have higher rates of successful duplications, where gene importance is measured by the fitness reduction caused by the deletion of the gene. Here, we test this hypothesis by comparing the importance of two groups of singleton genes in the yeast Saccharomyces cerevisiae (Sce). Group S genes did not duplicate in four other yeast species examined, whereas group D experienced duplication in these species. Consistent with our hypothesis, we found group D genes to be less important than group S genes. Specifically, 17% of group D genes are essential in Sce, compared to 28% for group S. Furthermore, deleting a group D gene in Sce reduces the fitness by 24% on average, compared to 38% for group S. Our subsequent analysis showed that less important genes have more cis-regulatory motifs, which could lead to a higher chance of subfunctionalization of duplicate genes and result in an enhanced rate of gene retention. Less important genes may also have weaker dosage imbalance effects and cause fewer genetic perturbations when duplicated. Regardless of the cause, our observation indicates that the previous finding of a less severe fitness consequence of deleting a duplicate gene than deleting a singleton gene is at least in part due to the fact that duplicate genes are intrinsically less important than singleton genes and suggests that the contribution of duplicate genes to genetic robustness has been overestimated.     

Dominance Genetic Variance for Traits Under Directional Selection in Drosophila serrata

In contrast to our growing understanding of patterns of additive genetic variance in single- and multi-trait combinations, the relative contribution of nonadditive genetic variance, particularly dominance variance, to multivariate phenotypes is largely unknown. While mechanisms for the evolution of dominance genetic variance have been, and to some degree remain, subject to debate, the pervasiveness of dominance is widely recognized and may play a key role in several evolutionary processes. Theoretical and empirical evidence suggests that the contribution of dominance variance to phenotypic variance may increase with the correlation between a trait and fitness; however, direct tests of this hypothesis are few. Using a multigenerational breeding design in an unmanipulated population of Drosophila serrata, we estimated additive and dominance genetic covariance matrices for multivariate wing-shape phenotypes, together with a comprehensive measure of fitness, to determine whether there is an association between directional selection and dominance variance. Fitness, a trait unequivocally under directional selection, had no detectable additive genetic variance, but significant dominance genetic variance contributing 32% of the phenotypic variance. For single and multivariate morphological traits, however, no relationship was observed between trait–fitness correlations and dominance variance. A similar proportion of additive and dominance variance was found to contribute to phenotypic variance for single traits, and double the amount of additive compared to dominance variance was found for the multivariate trait combination under directional selection. These data suggest that for many fitness components a positive association between directional selection and dominance genetic variance may not be expected.     

A Gene Family Derived from Transposable Elements during Early Angiosperm Evolution Has Reproductive Fitness Benefits in Arabidopsis thaliana

The benefits of ever-growing numbers of sequenced eukaryotic genomes will not be fully realized until we learn to decipher vast stretches of noncoding DNA, largely composed of transposable elements. Transposable elements persist through self-replication, but some genes once encoded by transposable elements have, through a process called molecular domestication, evolved new functions that increase fitness. Although they have conferred numerous adaptations, the number of such domesticated transposable element genes remains unknown, so their evolutionary and functional impact cannot be fully assessed. Systematic searches that exploit genomic signatures of natural selection have been employed to identify potential domesticated genes, but their predictions have yet to be experimentally verified. To this end, we investigated a family of domesticated genes called MUSTANG (MUG), identified in a previous bioinformatic search of plant genomes. We show that MUG genes are functional. Mutants of Arabidopsis thaliana MUG genes yield phenotypes with severely reduced plant fitness through decreased plant size, delayed flowering, abnormal development of floral organs, and markedly reduced fertility. MUG genes are present in all flowering plants, but not in any non-flowering plant lineages, such as gymnosperms, suggesting that the molecular domestication of MUG may have been an integral part of early angiosperm evolution. This study shows that systematic searches can be successful at identifying functional genetic elements in noncoding regions and demonstrates how to combine systematic searches with reverse genetics in a fruitful way to decipher eukaryotic genomes.     

The potential for sexually antagonistic polymorphism in different genome regions

Sex differences in the fitness effects of alleles at a single locus (intralocus sexual antagonism, or SA) have several evolutionary consequences. Among the consequences of SA, polymorphisms at genes partially linked to the sex-determining region of the sex chromosome pair potentially drive the evolution of suppressed recombination between the sex chromosomes. Understanding the conditions under which SA polymorphism can exist at such pseudo-autosomal (or PAR) loci should increase understanding of the evolution of recombination between sex chromosome pairs, and can help predict when we may expect potentially empirically detectable allele frequency differences between the sexes. Models so far published have concluded that PAR genes can maintain SA polymorphisms over a wider range of selection coefficients than autosomal ones, but have used restrictive assumptions. We expand the modeling of SA alleles at a single locus with the full range of degrees of linkage to the male-specific region, to include strong or weak selection and the possibility of different dominance coefficients in the two sexes. We confirm the previous major conclusion that SA polymorphisms are generally maintained in a larger region of parameter space if the locus is in the PAR than if it is autosomal.     

Quantitative target display: a method to screen yeast mutants conferring quantitative phenotypes by 'mutant DNA fingerprints'

Whole genome sequencing of several microbes has revealed thousands of genes of unknown function. A large proportion of these genes seem to confer subtle quantitative phenotypes or phenotypes that do not have a plate screen. We report a novel method to monitor such phenotypes, where the fitness of mutants is assessed in mixed cultures under competitive growth conditions, and the abundance of any individual mutant in the pool is followed by means of its unique feature, namely the mutation itself. A mixed population of yeast mutants, obtained through transposon mutagenesis, was subjected to selection. The DNA regions (targets) flanking the transposon, until nearby restriction sites, are then quantitatively amplified by means of a ligation-mediated PCR method, using transposon-specific and adapter-specific primers. The amplified PCR products correspond to mutated regions of the genome and serve as 'mutant DNA fingerprints' that can be displayed on a sequencing gel. The relative intensity of the amplified DNA fragments before and after selection match with the relative abundance of corresponding mutants, thereby revealing the fate of the mutants during selection. Using this method we demonstrate that UBI4, YDJ1 and HSP26 are essential for stress tolerance of yeast during ethanol production. We anticipate that this method will be useful for functional analysis of genes of any microbe amenable to insertional mutagenesis.     

Modelling the maintenance of male-fertile cytoplasm in a gynodioecious population

Gynodioecy is the co-occurrence of females and hermaphrodites in populations. It is usually due to the combined action of cytoplasmic male sterility (CMS) genes and nuclear genes that restore male fertility. According to previous theoretical studies, it is very difficult to explain the maintenance of gynodioecy with CMS and male-fertile cytotypes, although it has been observed in some species. However, only very specific situations have been investigated so far. We present a model to investigate the conditions that promote the maintenance of this breeding system in the case of an outcrossed species when CMS and male-fertile (non-CMS) cytotypes are present in an infinite panmictic population. We show that the type of cost of restoration strongly affects the conditions for stable maintenance of gynodioecy. Stable nuclear-cytoplasmic gynodioecy requires a female advantage, which is a classical condition for gynodioecy, but also a cost of CMS for female fitness, which had been rarely investigated. A cost of restoration is also needed, which could affect either pollen or seeds. Finally, we found that gynodioecy was attainable for a large set of parameter values, including low differences in fitness among genotypes and phenotypes. Our theoretical predictions are compared with previous theoretical work and with results of empirical studies on various gynodioecious species.     

Neutrality of foreign complex subunits in an experimental model of lateral gene transfer

Lateral gene transfer (LGT) is a powerful force in microbial evolution. However, the barriers that restrict this evolutionary phenomenon are not fully understood. It has long been observed that genes that encode subunits of complexes exhibit relatively compatible phylogenies, implying mostly vertical evolution. This may be explained by the failure of a new gene product to effectively interact with preexisting protein subunits, making its acquisition neutral--a theory termed the "complexity hypothesis." On the other hand, such genes may reduce the fitness of the host by disturbing the stoichiometric balance between complex subunits, resulting in purifying selection against gene retention. To examine these 2 alternative scenarios, we designed an experimental system that mimics the transfer of genes encoding homologs of essential complex subunits into the model bacterium Escherichia coli. In addition, we overexpressed the native E. coli gene in order to examine the contribution of gene dosage effects. We show that accumulation of native or foreign complex subunits in the cell does not result in loss of fitness, except for a minor fitness reduction observed for a single foreign homolog. Indeed, a series of genetic and biochemical assays failed to detect any interaction between the foreign subunits and the native polypeptides of the complex, implying an inability of such transfer events to generate positive selection for gene retention. We conclude that LGT of complex subunits may be mostly neutral and that forces operating against gene retention appear to be moderate.     

Natural variants suppress mutations in hundreds of essential genes

The consequence of a mutation can be influenced by the context in which it operates. For example, loss of gene function may be tolerated in one genetic background, and lethal in another. The extent to which mutant phenotypes are malleable, the architecture of modifiers and the identities of causal genes remain largely unknown. Here, we measure the fitness effects of ~ 1,100 temperature‐sensitive alleles of yeast essential genes in the context of variation from ten different natural genetic backgrounds and map the modifiers for 19 combinations. Altogether, fitness defects for 149 of the 580 tested genes (26%) could be suppressed by genetic variation in at least one yeast strain. Suppression was generally driven by gain‐of‐function of a single, strong modifier gene, and involved both genes encoding complex or pathway partners suppressing specific temperature‐sensitive alleles, as well as general modifiers altering the effect of many alleles. The emerging frequency of suppression and range of possible mechanisms suggest that a substantial fraction of monogenic diseases could be managed by modulating other gene products.     

Insertional mutagenesis in zebrafish: genes for development, genes for disease.

In order to rapidly identify a substantial fraction of the genes with a unique and essential role in vertebrate development, the laboratory of Nancy Hopkins at MIT has performed a large insertional mutagenesis screen in zebrafish using a pseudotyped retroviral vector as the mutagen. We have recovered mutations in about one-quarter of the embryonic essential genes in this organism, and have identified the mutated genes in nearly all of these (333). As the ease of gene identification allowed us to clone the mutated genes for nearly all of the mutants rather than prioritizing based upon the initially observed phenotypes, this has provided an unbiased view of the diversity of genes required for vertebrate development as well as a large collection of mutants to be screened for more specific phenotypes. In collaboration with other labs, we have screened the insertional mutant for the development of a variety of organs and cell types, as well as phenotypes that could represent disease models, such as cystic kidney and hepatomegaly. Furthermore, while all of these mutants are embryonic lethal in their homozygous state, we are investigating the heterozygous adults for additional phenotypes, such as cancer predisposition.