Predicting genome-wide redundancy using machine learning

Background  

Gene duplication can lead to genetic redundancy, which masks the function of mutated genes in genetic analyses. Methods to increase sensitivity in identifying genetic redundancy can improve the efficiency of reverse genetics and lend insights into the evolutionary outcomes of gene duplication. Machine learning techniques are well suited to classifying gene family members into redundant and non-redundant gene pairs in model species where sufficient genetic and genomic data is available, such as Arabidopsis thaliana, the test case used here.     

Unequal genetic redundancies in Arabidopsis--a neglected phenomenon?

Genetic redundancy is a common phenomenon in Arabidopsis and is thought to be responsible for the absence of phenotypes in the majority of single loss-of-function mutants. In this review, we highlight an increasing number of examples in which redundancy between homologous genes is limited or absent despite functional equivalence of the respective proteins. In particular, we focus on cases of unequal redundancy, where the absence of a mutant phenotype in loss-of-function mutants of one gene contrasts with a strong phenotype in mutants of its homolog. In the double mutants, this phenotype is strongly enhanced. Possible explanations for such scenarios are discussed. We propose that the study of unequally redundant gene pairs offers a unique opportunity to understand global patterns of functional genome evolution.     

The cellular robustness by genetic redundancy in budding yeast

The frequent dispensability of duplicated genes in budding yeast is heralded as a hallmark of genetic robustness contributed by genetic redundancy. However, theoretical predictions suggest such backup by redundancy is evolutionarily unstable, and the extent of genetic robustness contributed from redundancy remains controversial. It is anticipated that, to achieve mutual buffering, the duplicated paralogs must at least share some functional overlap. However, counter-intuitively, several recent studies reported little functional redundancy between these buffering duplicates. The large yeast genetic interactions released recently allowed us to address these issues on a genome-wide scale. We herein characterized the synthetic genetic interactions for ～500 pairs of yeast duplicated genes originated from either whole-genome duplication (WGD) or small-scale duplication (SSD) events. We established that functional redundancy between duplicates is a pre-requisite and thus is highly predictive of their backup capacity. This observation was particularly pronounced with the use of a newly introduced metric in scoring functional overlap between paralogs on the basis of gene ontology annotations. Even though mutual buffering was observed to be prevalent among duplicated genes, we showed that the observed backup capacity is largely an evolutionarily transient state. The loss of backup capacity generally follows a neutral mode, with the buffering strength decreasing in proportion to divergence time, and the vast majority of the paralogs have already lost their backup capacity. These observations validated previous theoretic predictions about instability of genetic redundancy. However, departing from the general neutral mode, intriguingly, our analysis revealed the presence of natural selection in stabilizing functional overlap between SSD pairs. These selected pairs, both WGD and SSD, tend to have decelerated functional evolution, have higher propensities of co-clustering into the same protein complexes, and share common interacting partners. Our study revealed the general principles for the long-term retention of genetic redundancy.     

Gene duplication in the evolution of sexual dimorphism

Males and females share most of the same genes, so selection in one sex will typically produce a correlated response in the other sex. Yet, the sexes have evolved to differ in a multitude of behavioral, morphological, and physiological traits. How did this sexual dimorphism evolve despite the presence of a common underlying genome? We investigated the potential role of gene duplication in the evolution of sexual dimorphism. Because duplication events provide extra genetic material, the sexes each might use this redundancy to facilitate sex‐specific gene expression, permitting the evolution of dimorphism. We investigated this hypothesis at the genome‐wide level in Drosophila melanogaster, using the presence of sex‐biased expression as a proxy for the sex‐specific specialization of gene function. We expected that if sexually antagonistic selection is a potent force acting upon individual genes, duplication will result in paralog families whose members differ in sex‐biased expression. Gene members of the same duplicate family can have different expression patterns in males versus females. In particular, duplicate pairs containing a male‐biased gene are found more frequently than expected, in agreement with previous studies. Furthermore, when the singleton ortholog is unbiased, duplication appears to allow one of the paralog copies to acquire male‐biased expression. Conversely, female‐biased expression is not common among duplicates; fewer duplicate genes are expressed in the female‐soma and ovaries than in the male‐soma and testes. Expression divergence exists more in older than in younger duplicates pairs, but expression divergence does not correlate with protein sequence divergence. Finally, genomic proximity may have an effect on whether paralogs differ in sex‐biased expression. We conclude that the data are consistent with a role of gene duplication in fostering male‐biased, but not female‐biased, gene expression, thereby aiding the evolution of sexual dimorphism.     

Evolutionary Analysis of the Two-Component Systems in <Emphasis Type="Italic">Pseudomonas aeruginosa</Emphasis> PAO1

Gene organization and functional motif analyses of the 123 two-component system (2CS) genes in Pseudomonas aeruginosa PAO1 were carried out. In addition, NJ and ML trees for the sensor kinases and the response regulators were constructed, and the distances measured and comparatively analyzed. It was apparent that more than half of the sensor-regulator gene pairs, especially the 2CSs with OmpR-like regulators, are derivatives of a common ancestor and have most likely co-evolved through gene pair duplication. Several of the 2CS pairs, especially those with NarL-like regulators, however, appeared to be relatively divergent. This is supportive of the recruitment model, in which a sensor gene and regulator gene with different phylogenetic history are assembled to form a 2CS. Correlation of the classification of sensor kinases and response regulators provides further support for these models. Upon comparison of the phylogenetic trees comprised of sensors and regulators, we have identified six congruent clades, which represent the group of the most recently duplicated 2CS gene pairs. Analyses of the congruent 2CS pairs of each of the clades revealed that certain paralogous 2CS pairs may carry a redundant function even after a gene duplication event. Nevertheless, comparative analysis of the putative promoter regions of the paralogs suggested that functional redundancy could be prevented by a differential control. Both codon usage and G+C content of these 2CS genes were found to be comparable with those of the P. aeruginosa genome, suggesting that they are not newly acquired genes.Reviewing Editor: Dr. Martin Kreitman     

A Minimal Set of Glycolytic Genes Reveals Strong Redundancies in Saccharomyces cerevisiae Central Metabolism

As a result of ancestral whole-genome and small-scale duplication events, the genomes of Saccharomyces cerevisiae and many eukaryotes still contain a substantial fraction of duplicated genes. In all investigated organisms, metabolic pathways, and more particularly glycolysis, are specifically enriched for functionally redundant paralogs. In ancestors of the Saccharomyces lineage, the duplication of glycolytic genes is purported to have played an important role leading to S. cerevisiae''s current lifestyle favoring fermentative metabolism even in the presence of oxygen and characterized by a high glycolytic capacity. In modern S. cerevisiae strains, the 12 glycolytic reactions leading to the biochemical conversion from glucose to ethanol are encoded by 27 paralogs. In order to experimentally explore the physiological role of this genetic redundancy, a yeast strain with a minimal set of 14 paralogs was constructed (the “minimal glycolysis” [MG] strain). Remarkably, a combination of a quantitative systems approach and semiquantitative analysis in a wide array of growth environments revealed the absence of a phenotypic response to the cumulative deletion of 13 glycolytic paralogs. This observation indicates that duplication of glycolytic genes is not a prerequisite for achieving the high glycolytic fluxes and fermentative capacities that are characteristic of S. cerevisiae and essential for many of its industrial applications and argues against gene dosage effects as a means of fixing minor glycolytic paralogs in the yeast genome. The MG strain was carefully designed and constructed to provide a robust prototrophic platform for quantitative studies and has been made available to the scientific community.     

Evidence for short-time divergence and long-time conservation of tissue-specific expression after gene duplication

Gene duplication is one of the main mechanisms by which genomes can acquire novel functions. It has been proposed that the retention of gene duplicates can be associated to processes of tissue expression divergence. These models predict that acquisition of divergent expression patterns should be acquired shortly after the duplication, and that larger divergence in tissue expression would be expected for paralogs, as compared to orthologs of a similar age. Many studies have shown that gene duplicates tend to have divergent expression patterns and that gene family expansions are associated with high levels of tissue specificity. However, the timeframe in which these processes occur have rarely been investigated in detail, particularly in vertebrates, and most analyses do not include direct comparisons of orthologs as a baseline for the expected levels of tissue specificity in absence of duplications. To assess the specific contribution of duplications to expression divergence, we combine here phylogenetic analyses and expression data from human and mouse. In particular, we study differences in spatial expression among human-mouse paralogs, specifically duplicated after the radiation of mammals, and compare them to pairs of orthologs in the same species. Our results show that gene duplication leads to increased levels of tissue specificity and that this tends to occur promptly after the duplication event.     

Evolution by Gene Duplication and Compensatory Advantageous Mutations

Relaxation of selective constraint is thought to play an important role for evolution by gene duplication, in connection with compensatory advantageous mutant substitutions. Models were investigated by incorporating gene duplication by unequal crossing over, selection, mutation and random genetic drift into Monte Carlo simulations. Compensatory advantageous mutations were introduced, and simulations were carried out with and without relaxation, when genes are redundant on chromosomes. Relaxation was introduced by assuming that deleterious mutants have no effect on fitness, so long as one or more genes free of such mutations remain in the array. Compensatory mutations are characterized by the intermediate deleterious step of their substitutions, and therefore relaxation by gene redundancy is important. Through extensive Monte Carlo simulations, it was found that compensatory mutant substitutions require relaxation in addition to gene duplication, when mutant effects are large. However when mutant effects are small, such that the product of selection coefficient and population size is around unity, evolution by compensatory mutation is enhanced by gene duplication even without relaxation.     

Analysis of sequence, map position, and gene expression reveals conserved essential genes for iron uptake in Arabidopsis and tomato

Arabidopsis (Arabidopsis thaliana) and tomato (Lycopersicon esculentum) show similar physiological responses to iron deficiency, suggesting that homologous genes are involved. Essential gene functions are generally considered to be carried out by orthologs that have remained conserved in sequence and map position in evolutionarily related species. This assumption has not yet been proven for plant genomes that underwent large genome rearrangements. We addressed this question in an attempt to deduce functional gene pairs for iron reduction, iron transport, and iron regulation between Arabidopsis and tomato. Iron uptake processes are essential for plant growth. We investigated iron uptake gene pairs from tomato and Arabidopsis, namely sequence, conserved gene content of the regions containing iron uptake homologs based on conserved orthologous set marker analysis, gene expression patterns, and, in two cases, genetic data. Compared to tomato, the Arabidopsis genome revealed more and larger gene families coding for the iron uptake functions. The number of possible homologous pairs was reduced if functional expression data were taken into account in addition to sequence and map position. We predict novel homologous as well as partially redundant functions of ferric reductase-like and iron-regulated transporter-like genes in Arabidopsis and tomato. Arabidopsis nicotianamine synthase genes encode a partially redundant family. In this study, Arabidopsis gene redundancy generally reflected the presumed genome duplication structure. In some cases, statistical analysis of conserved gene regions between tomato and Arabidopsis suggested a common evolutionary origin. Although involvement of conserved genes in iron uptake was found, these essential genes seem to be of paralogous rather than orthologous origin in tomato and Arabidopsis.     

Redundant and non-redundant roles of the trehalose-6-phosphate phosphatases in leaf growth,root hair specification and energy-responses in Arabidopsis

The Arabidopsis trehalose-6-phosphate phosphatase (TPP) gene family arose mainly from whole genome duplication events and consists of 10 genes (TPPA-J). All the members encode active TPP enzymes, possibly regulating the levels of trehalose-6-phosphate, an established signaling metabolite in plants. GUS activity studies revealed tissue-, cell- and stage-specific expression patterns for the different members of the TPP gene family. Here we list additional examples of the remarkable features of the TPP gene family. TPPA-J expression levels seem, in most of the cases, differently regulated in response to light, darkness and externally supplied sucrose. Disruption of the TPPB gene leads to Arabidopsis plants with larger leaves, which is the result of an increased cell number in the leaves. Arabidopsis TPPA and TPPG are preferentially expressed in atrichoblast cells. TPPA and TPPG might fulfill redundant roles during the differentiation process of root epidermal cells, since the tppa tppg double mutant displays a hairy root phenotype, while the respective single knockouts have a distribution of trichoblast and atrichoblast cells similar to the wild type. These new data portray redundant and non-redundant functions of the TPP proteins in regulatory pathways of Arabidopsis.     

Retention of protein complex membership by ancient duplicated gene products in budding yeast

To investigate functional divergence of gene duplicates, we examined the protein-protein interactions and coexistence in complexes of paralogs resulting from an ancient whole-genome duplication in yeast. Strikingly, half the surveyed paralog pairs were found to be co-clustered in protein complexes, and were more conserved and highly expressed than non-co-clustered paralogs; however, their discordant expression patterns and conservation rates indicate differential regulation of subfunctionalized paralogs. These results highlight the value of protein complex membership in studying functional divergence among gene duplicates.     

Widespread conservation of genetic redundancy during a billion years of eukaryotic evolution

Genetic redundancy means that two genes can perform the same function. Using a comprehensive phylogenetic analysis, we show here in both Saccharomyces cerevisiae and Caenorhabditis elegans that genetic redundancy is not just a transient consequence of gene duplication, but is often an evolutionary stable state. In multiple examples, genes have retained redundant functions since the divergence of the animal, plant and fungi kingdoms over a billion years ago. The stable conservation of genetic redundancy contrasts with the more rapid evolution of genetic interactions between unrelated genes and can be explained by theoretical models including a 'piggyback' mechanism in which overlapping redundant functions are co-selected with nonredundant ones.     

Multiple mechanisms promote the retained expression of gene duplicates in the tetraploid frog Xenopus laevis

Gene duplication provides a window of opportunity for biological variants to persist under the protection of a co-expressed copy with similar or redundant function. Duplication catalyzes innovation (neofunctionalization), subfunction degeneration (subfunctionalization), and genetic buffering (redundancy), and the genetic survival of each paralog is triggered by mechanisms that add, compromise, or do not alter protein function. We tested the applicability of three types of mechanisms for promoting the retained expression of duplicated genes in 290 expressed paralogs of the tetraploid clawed frog, Xenopus laevis. Tests were based on explicit expectations concerning the ka/ks ratio, and the number and location of nonsynonymous substitutions after duplication. Functional constraints on the majority of paralogs are not significantly different from a singleton ortholog. However, we recover strong support that some of them have an asymmetric rate of nonsynonymous substitution: 6% match predictions of the neofunctionalization hypothesis in that (1) each paralog accumulated nonsynonymous substitutions at a significantly different rate and (2) the one that evolves faster has a higher ka/ks ratio than the other paralog and than a singleton ortholog. Fewer paralogs (3%) exhibit a complementary pattern of substitution at the protein level that is predicted by enhancement or degradation of different functional domains, and the remaining 13% have a higher average ka/ks ratio in both paralogs that is consistent with altered functional constraints, diversifying selection, or activity-reducing mutations after duplication. We estimate that these paralogs have been retained since they originated by genome duplication between 21 and 41 million years ago. Multiple mechanisms operate to promote the retained expression of duplicates in the same genome, in genes in the same functional class, over the same period of time following duplication, and sometimes in the same pair of paralogs. None of these paralogs are superfluous; degradation or enhancement of different protein subfunctions and neofunctionalization are plausible hypotheses for the retained expression of some of them. Evolution of most X. laevis paralogs, however, is consistent with retained expression via mechanisms that do not radically alter functional constraints, such as selection to preserve post-duplication stoichiometry or temporal, quantitative, or spatial subfunctionalization.     

Commensurate distances and similar motifs in genetic congruence and protein interaction networks in yeast

Background  

In a genetic interaction, the phenotype of a double mutant differs from the combined phenotypes of the underlying single mutants. When the single mutants have no growth defect, but the double mutant is lethal or exhibits slow growth, the interaction is termed synthetic lethality or synthetic fitness. These genetic interactions reveal gene redundancy and compensating pathways. Recently available large-scale data sets of genetic interactions and protein interactions in Saccharomyces cerevisiae provide a unique opportunity to elucidate the topological structure of biological pathways and how genes function in these pathways.     

Functional analysis of gene duplications in Saccharomyces cerevisiae

Gene duplication can occur on two scales: whole-genome duplications (WGD) and smaller-scale duplications (SSD) involving individual genes or genomic segments. Duplication may result in functionally redundant genes or diverge in function through neofunctionalization or subfunctionalization. The effect of duplication scale on functional evolution has not yet been explored, probably due to the lack of global knowledge of protein function and different times of duplication events. To address this question, we used integrated Bayesian analysis of diverse functional genomic data to accurately evaluate the extent of functional similarity and divergence between paralogs on a global scale. We found that paralogs resulting from the whole-genome duplication are more likely to share interaction partners and biological functions than smaller-scale duplicates, independent of sequence similarity. In addition, WGD paralogs show lower frequency of essential genes and higher synthetic lethality rate, but instead diverge more in expression pattern and upstream regulatory region. Thus, our analysis demonstrates that WGD paralogs generally have similar compensatory functions but diverging expression patterns, suggesting a potential of distinct evolutionary scenarios for paralogs that arose through different duplication mechanisms. Furthermore, by identifying these functional disparities between the two types of duplicates, we reconcile previous disputes on the relationship between sequence divergence and expression divergence or essentiality.     

Robustness of Helicobacter pylori Infection Conferred by Context-Variable Redundancy among Cysteine-Rich Paralogs

Deletion of single genes from expanded gene families in bacterial genomes often does not elicit a phenotype thus implying redundancy or functional non-essentiality of paralogous genes. The molecular mechanisms that facilitate evolutionary maintenance of such paralogs despite selective pressures against redundancy remain mostly unexplored. Here, we investigate the evolutionary, genetic, and functional interaction between the Helicobacter pylori cysteine-rich paralogs hcpG and hcpC in the context of H. pylori infection of cultured mammalian cells. We find that in natural H. pylori populations both hcpG and hcpC are maintained by positive selection in a dual genetic relationship that switches from complete redundancy during early infection, whereby ΔhcpC or ΔhcpG mutants themselves show no growth defect but a significant growth defect is seen in the ΔhcpC,ΔhcpG double mutant, to quantitative redundancy during late infection wherein the growth defect of the ΔhcpC mutant is exacerbated in the ΔhcpC,ΔhcpG double mutant although the ΔhcpG mutant itself shows no defect. Moreover, during early infection both hcpG and hcpC are essential for optimal translocation of the H. pylori HspB/GroEL chaperone, but during middle-to-late infection hcpC alone is necessary and sufficient for HspB/GroEL translocation thereby revealing the lack of functional compensation among paralogs. We propose that evolution of context-dependent differences in the nature of genetic redundancy, and function, between hcpG and hcpC may facilitate their maintenance in H. pylori genomes, and confer robustness to H. pylori growth during infection of cultured mammalian cells.