Patterns of evolutionary conservation of essential genes correlate with their compensability

Evolutionary adaptation after crippling cell polarization follows reproducible trajectories

Cells are organized by functional modules, which typically contain components whose removal severely compromises the module's function. Despite their importance, these components are not absolutely conserved between parts of the tree of life, suggesting that cells can evolve to perform the same biological functions with different proteins. We evolved Saccharomyces cerevisiae for 1000 generations without the important polarity gene BEM1. Initially the bem1∆ lineages rapidly increase in fitness and then slowly reach >90% of the fitness of their BEM1 ancestors at the end of the evolution. Sequencing their genomes and monitoring polarization reveals a common evolutionary trajectory, with a fixed sequence of adaptive mutations, each improving cell polarization by inactivating proteins. Our results show that organisms can be evolutionarily robust to physiologically destructive perturbations and suggest that recovery by gene inactivation can lead to rapid divergence in the parts list for cell biologically important functions.

DOI:http://dx.doi.org/10.7554/eLife.09638.001

Cells use the genetic instructions provided by genes in particular combinations called ‘modules’ to perform particular jobs. Very different organisms can share many of the same modules because certain abilities are fundamental to the survival of all cells and so they have been retained over the course of evolution. That said, these modules may not necessarily involve the same genes because it is often possible to achieve the same result using different components.

One way to study how those modules can diversify is to deliberately disrupt one of the genes in a module, and observe how the organism and its descendants respond over many generations. Other genes in these organisms may acquire genetic mutations that enable the genes to take on the role of the missing protein. However, the removal of a single component can be detrimental to the survival of the organisms or may affect many different processes. This can make it difficult to understand what is going on.

A gene called BEM1 is crucial for yeast cells to establish polarity, that is, to allow the different sides of a cell to become distinct from one another. This activity is essential for the yeast to replicate itself. Previous studies have shown that the BEM1 gene had a different role in other species of fungi, which suggests that yeast may have other genes that previously assumed the role that BEM1 does now. In this study, Laan et al. removed BEM1 from yeast and allowed the population of mutant cells to evolve for a thousand generations. The approach differs from previous studies because Laan et al. deliberately selected for yeast that had acquired multiple genetic mutations that can together almost fully compensate for the loss of BEM1.

Initially, the mutant cells grew very slowly, were abnormal in shape and likely to burst open. However, by the end of the experiment, the cells were able to grow almost as well as the original yeast cells had before the gene deletion. Genetic analysis revealed that the deletion of BEM1 triggers the inactivation of other genes that are also involved in the regulation of polarity, which largely restored the ability of the disrupted polarity module to work. This restoration follows a ‘reproducible trajectory’, as the same genes were switched off in the same order in different populations of yeast that were studied at the same time.

The work is an example of reproducible evolution, whereby a specific order of changes to gene activity repeatedly enables cells with severe defects in important processes to adapt and restore a gene module, using whatever components they have left. The next challenge will be to understand how the particular roles of important modules affect their adaptability.

DOI:http://dx.doi.org/10.7554/eLife.09638.002

Genomic Adaptations to the Loss of a Conserved Bacterial DNA Methyltransferase

CcrM is an orphan DNA methyltransferase nearly universally conserved in a vast group of Alphaproteobacteria. In Caulobacter crescentus, it controls the expression of key genes involved in the regulation of the cell cycle and cell division. Here, we demonstrate, using an experimental evolution approach, that C. crescentus can significantly compensate, through easily accessible genetic changes like point mutations, the severe loss in fitness due to the absence of CcrM, quickly improving its growth rate and cell morphology in rich medium. By analyzing the compensatory mutations genome-wide in 12 clones sampled from independent ΔccrM populations evolved for ~300 generations, we demonstrated that each of the twelve clones carried at least one mutation that potentially stimulated ftsZ expression, suggesting that the low intracellular levels of FtsZ are the major burden of ΔccrM mutants. In addition, we demonstrate that the phosphoenolpyruvate-carbohydrate phosphotransfer system (PTS) actually modulates ftsZ and mipZ transcription, uncovering a previously unsuspected link between metabolic regulation and cell division in Alphaproteobacteria. We present evidence that point mutations found in genes encoding proteins of the PTS provide the strongest fitness advantage to ΔccrM cells cultivated in rich medium despite being disadvantageous in minimal medium. This environmental sign epistasis might prevent such mutations from getting fixed under changing natural conditions, adding a plausible explanation for the broad conservation of CcrM.

In bacteria, DNA methylation has a variety of functions, including the control of DNA replication and/or gene expression. The cell cycle-regulated DNA methyltransferase CcrM modulates the transcription of many genes and is critical for fitness in Caulobacter crescentus. Here, we used an original experimental evolution approach to determine which of its many targets make CcrM so important physiologically. We show that populations lacking CcrM evolve quickly, accumulating an excess of mutations affecting, directly or indirectly, the expression of the ftsZ cell division gene. This finding suggests that the most critical function of CcrM in C. crescentus is to promote cell division by enhancing FtsZ intracellular levels. During this work, we also discovered an unexpected link between metabolic regulation and cell division that might extend to other Alphaproteobacteria.

Exploring the utility of cross-laboratory RAD-sequencing datasets for phylogenetic analysis

BACKGROUND

Restriction site-Associated DNA sequencing (RAD-Seq) is widely applied to generate genome-wide sequence and genetic marker datasets. RAD-Seq has been extensively utilised, both at the population level and across species, for example in the construction of phylogenetic trees. However, the consistency of RAD-Seq data generated in different laboratories, and the potential use of cross-species orthologous RAD loci in the estimation of genetic relationships, have not been widely investigated. This study describes the use of SbfI RAD-Seq data for the estimation of evolutionary relationships amongst ten teleost fish species, using previously established phylogeny as a benchmark.

RESULTS

The number of orthologous SbfI RAD loci identified decreased with increasing evolutionary distance between the species, with several thousand loci conserved across five salmonid species (divergence ~50 MY), and several hundred conserved across the more distantly related teleost species (divergence ~100-360 MY). The majority (>70%) of loci identified between the more distantly related species were genic in origin, suggesting that the bias of SbfI towards genic regions is useful for identifying distant orthologs. Interspecific single nucleotide variants at each orthologous RAD locus were identified. Evolutionary relationships estimated using concatenated sequences of interspecific variants were congruent with previously published phylogenies, even for distantly (divergence up to ~360 MY) related species.

CONCLUSION

Overall, this study has demonstrated that orthologous SbfI RAD loci can be identified across closely and distantly related species. This has positive implications for the repeatability of SbfI RAD-Seq and its potential to address research questions beyond the scope of the original studies. Furthermore, the concordance in tree topologies and relationships estimated in this study with published teleost phylogenies suggests that similar meta-datasets could be utilised in the prediction of evolutionary relationships across populations and species with readily available RAD-Seq datasets, but for which relationships remain uncharacterised.

Statistical Analysis of Hurst Exponents of Essential/Nonessential Genes in 33 Bacterial Genomes

Methods for identifying essential genes currently depend predominantly on biochemical experiments. However, there is demand for improved computational methods for determining gene essentiality. In this study, we used the Hurst exponent, a characteristic parameter to describe long-range correlation in DNA, and analyzed its distribution in 33 bacterial genomes. In most genomes (31 out of 33) the significance levels of the Hurst exponents of the essential genes were significantly higher than for the corresponding full-gene-set, whereas the significance levels of the Hurst exponents of the nonessential genes remained unchanged or increased only slightly. All of the Hurst exponents of essential genes followed a normal distribution, with one exception. We therefore propose that the distribution feature of Hurst exponents of essential genes can be used as a classification index for essential gene prediction in bacteria. For computer-aided design in the field of synthetic biology, this feature can build a restraint for pre- or post-design checking of bacterial essential genes. Moreover, considering the relationship between gene essentiality and evolution, the Hurst exponents could be used as a descriptive parameter related to evolutionary level, or be added to the annotation of each gene.

Insights into the miRNA regulations in human disease genes

BACKGROUND

MicroRNAs are a class of short non-coding RNAs derived from either cellular or viral transcripts that act post-transcriptionally to regulate mRNA stability and translation. In recent days, increasing numbers of miRNAs have been shown to be involved in the development and progression of a variety of diseases. We, therefore, intend to enumerate miRNA targets in several known disease classes to explore the degree of miRNA regulations on them which is unexplored till date.

RESULTS

Here, we noticed that miRNA hits in cancer genes are remarkably higher than other diseases in human. Our observation suggests that UTRs and the transcript length of cancer related genes have a significant contribution in higher susceptibility to miRNA regulation. Moreover, gene duplication, mRNA stability, AREScores and evolutionary rate were likely to have implications for more miRNA targeting on cancer genes. Consequently, the regression analysis have confirmed that the AREScores plays most important role in detecting miRNA targets on disease genes. Interestingly, we observed that epigenetic modifications like CpG methylation and histone modification are less effective than miRNA regulations in controlling the gene expression of cancer genes.

CONCLUSIONS

The intrinsic properties of cancer genes studied here, for higher miRNA targeting will enhance the knowledge on cancer gene regulation.

Bacillus subtilis and Escherichia coli essential genes and minimal cell factories after one decade of genome engineering

Investigation of essential genes, besides contributing to understanding the fundamental principles of life, has numerous practical applications. Essential genes can be exploited as building blocks of a tightly controlled cell ‘chassis’. Bacillus subtilis and Escherichia coli K-12 are both well-characterized model bacteria used as hosts for a plethora of biotechnological applications. Determination of the essential genes that constitute the B. subtilis and E. coli minimal genomes is therefore of the highest importance. Recent advances have led to the modification of the original B. subtilis and E. coli essential gene sets identified 10 years ago. Furthermore, significant progress has been made in the area of genome minimization of both model bacteria. This review provides an update, with particular emphasis on the current essential gene sets and their comparison with the original gene sets identified 10 years ago. Special attention is focused on the genome reduction analyses in B. subtilis and E. coli and the construction of minimal cell factories for industrial applications.

Integrons: Vehicles and pathways for horizontal dissemination in bacteria

Integrons are genetic elements first described at the end of the 1980s. Although most integrons were initially described in human clinical isolates, they have now been identified in many non-clinical environments, such as water and soil. Integrons are present in ≈10% of the sequenced bacterial genomes and are frequently linked to mobile genetic elements (MGEs); particularly the class 1 integrons. Genetic linkage to a diverse set of MGEs facilitates horizontal transfer of class 1 integrons within and between bacterial populations and species. The mechanistic aspects limiting transfer of MGEs will therefore limit the transfer of class 1 integrons. However, horizontal movement due to genes provided in trans and homologous recombination can result in class 1 integron dynamics independent of MGEs. A key determinant for continued dissemination of class 1 integrons is the probability that transferred MGEs will be vertically inherited in the recipient bacterial population. Heritability depends both on genetic stability as well as the fitness costs conferred to the host. Here we review the factors known to govern the dissemination of class 1 integrons in bacteria.

The genomic landscape of compensatory evolution

The Genomic Landscape of Compensatory Evolution Laboratory selection experiment explains how organisms compensate for the loss of genes during evolution, and reveals the deleterious side-effects of this process when adapting to novel environments.

Adaptive evolution is generally assumed to progress through the accumulation of beneficial mutations. However, as deleterious mutations are common in natural populations, they generate a strong selection pressure to mitigate their detrimental effects through compensatory genetic changes. This process can potentially influence directions of adaptive evolution by enabling evolutionary routes that are otherwise inaccessible. Therefore, the extent to which compensatory mutations shape genomic evolution is of central importance. Here, we studied the capacity of the baker's yeast genome to compensate the complete loss of genes during evolution, and explored the long-term consequences of this process. We initiated laboratory evolutionary experiments with over 180 haploid baker's yeast genotypes, all of which initially displayed slow growth owing to the deletion of a single gene. Compensatory evolution following gene loss was rapid and pervasive: 68% of the genotypes reached near wild-type fitness through accumulation of adaptive mutations elsewhere in the genome. As compensatory mutations have associated fitness costs, genotypes with especially low fitnesses were more likely to be subjects of compensatory evolution. Genomic analysis revealed that as compensatory mutations were generally specific to the functional defect incurred, convergent evolution at the molecular level was extremely rare. Moreover, the majority of the gene expression changes due to gene deletion remained unrestored. Accordingly, compensatory evolution promoted genomic divergence of parallel evolving populations. However, these different evolutionary outcomes are not phenotypically equivalent, as they generated diverse growth phenotypes across environments. Taken together, these results indicate that gene loss initiates adaptive genomic changes that rapidly restores fitness, but this process has substantial pleiotropic effects on cellular physiology and evolvability upon environmental change. Our work also implies that gene content variation across species could be partly due to the action of compensatory evolution rather than the passive loss of genes.

While core cellular processes are generally conserved during evolution, the constituent genes differ somewhat between related species with similar lifestyles. Why should this be so? In this work, we propose that gene loss may initially be deleterious, but organisms can recover fitness by the accumulation of compensatory mutations elsewhere in the genome. To investigate this process in the laboratory, we investigated 180 haploid yeast strains, each of which initially displayed slow growth owing to the deletion of a single gene. Laboratory evolutionary experiments revealed that defects in a broad range of molecular processes can readily be compensated during evolution. Genomic analyses and functional assays demonstrated that compensatory evolution generates hidden genetic and physiological variation across parallel evolving lines, which can be revealed when the environment changes. Strikingly, despite nearly full recovery of fitness, the wild-type genomic expression pattern is generally not restored. Based on these results, we argue that genomes undergo major changes not simply to adapt to external conditions but also to compensate for previously accumulated deleterious mutations.

Geptop: a gene essentiality prediction tool for sequenced bacterial genomes based on orthology and phylogeny

Integrative genomics predictors, which score highly in predicting bacterial essential genes, would be unfeasible in most species because the data sources are limited. We developed a universal approach and tool designated Geptop, based on orthology and phylogeny, to offer gene essentiality annotations. In a series of tests, our Geptop method yielded higher area under curve (AUC) scores in the receiver operating curves than the integrative approaches. In the ten-fold cross-validations among randomly upset samples, Geptop yielded an AUC of 0.918, and in the cross-organism predictions for 19 organisms Geptop yielded AUC scores between 0.569 and 0.959. A test applied to the very recently determined essential gene dataset from the Porphyromonas gingivalis, which belongs to a phylum different with all of the above 19 bacterial genomes, gave an AUC of 0.77. Therefore, Geptop can be applied to any bacterial species whose genome has been sequenced. Compared with the essential genes uniquely identified by the lethal screening, the essential genes predicted only by Gepop are associated with more protein-protein interactions, especially in the three bacteria with lower AUC scores (<0.7). This may further illustrate the reliability and feasibility of our method in some sense. The web server and standalone version of Geptop are available at http://cefg.uestc.edu.cn/geptop/ free of charge. The tool has been run on 968 bacterial genomes and the results are accessible at the website.

Conserved expression of natural antisense transcripts in mammals

Background

Recent studies had found thousands of natural antisense transcripts originating from the same genomic loci of protein coding genes but from the opposite strand. It is unclear whether the majority of antisense transcripts are functional or merely transcriptional noise.

Results

Using the Affymetrix Exon array with a modified cDNA synthesis protocol that enables genome-wide detection of antisense transcription, we conducted large-scale expression analysis of antisense transcripts in nine corresponding tissues from human, mouse and rat. We detected thousands of antisense transcripts, some of which show tissue-specific expression that could be subjected to further study for their potential function in the corresponding tissues/organs. The expression patterns of many antisense transcripts are conserved across species, suggesting selective pressure on these transcripts. When compared to protein-coding genes, antisense transcripts show a lesser degree of expression conservation. We also found a positive correlation between the sense and antisense expression across tissues.

Conclusion

Our results suggest that natural antisense transcripts are subjected to selective pressure but to a lesser degree compared to sense transcripts in mammals.

Essential genes in Bacillus subtilis: a re-evaluation after ten years

In 2003, an initial study on essential genes in the Gram-positive model bacterium described 271 genes as essential. In the past decade, the functions of many unknown genes and their encoded proteins have been elucidated. Moreover, detailed analyses have revealed that 31 genes that were thought to be essential are in fact non-essential whereas 20 novel essential genes have been described. Thus, 261 genes coding for 259 proteins and two functional RNAs are regarded essential as of January 2013. Among the essential proteins, the largest group is involved in protein synthesis, secretion and protein quality control. Other large sets of essential proteins are involved in lipid biosynthesis, cell wall metabolism and cell division, and DNA replication. Another interesting group of essential proteins protects the cell against endogenous toxic proteins, metabolites, or other intermediates. There are only six essential proteins in B. subtilis, for which no function is known. The functional analysis of these important proteins is predicted to be a key issue in the research on this model organism in the coming years.

Three independent determinants of protein evolutionary rate

One of the most widely accepted ideas related to the evolutionary rates of proteins is that functionally important residues or regions evolve slower than other regions, a reasonable outcome of which should be a slower evolutionary rate of the proteins with a higher density of functionally important sites. Oddly, the role of functional importance, mainly measured by essentiality, in determining evolutionary rate has been challenged in recent studies. Several variables other than protein essentiality, such as expression level, gene compactness, protein-protein interactions, etc., have been suggested to affect protein evolutionary rate. In the present review, we try to refine the concept of functional importance of a gene, and consider three factors-functional importance, expression level, and gene compactness, as independent determinants of evolutionary rate of a protein, based not only on their known correlation with evolutionary rate but also on a reasonable mechanistic model. We suggest a framework based on these mechanistic models to correctly interpret the correlations between evolutionary rates and the various variables as well as the interrelationships among the variables.