FACT complex gene duplicates exhibit redundant and non-redundant functions in C. elegans.

FACT (facilitates chromatin transcription) is a histone chaperone complex important in genomic processes including transcription, DNA replication, and DNA repair. FACT is composed of two proteins, SSRP1 and SPT16, which are highly conserved across eukaryotes. While the mechanisms for FACT in nucleosome reorganization and its relationship to DNA processes is well established, how these roles impact coordination in multicellular animal development are less well understood. Here we characterize the genes encoding FACT complex proteins in the nematode C. elegans. We show that whereas C. elegans includes one SPT16 gene (spt-16), two genes (hmg-3 and hmg-4) encode SSRP1 proteins. Depletion of FACT complex genes interferes with embryonic cell division and cell cycle timing generally, with anterior pharynx development especially sensitive to these defects. hmg-3 and hmg-4 exhibit redundancy for these maternally-provided embryonic functions, but are each uniquely required zygotically for normal germline development. This work provides a framework to study FACT gene function in developmental processes, and identifies that distinct functional requirements for gene duplicates can be manifest within a single tissue.

Plasticity and innovation of regulatory mechanisms underlying seed oil content mediated by duplicated genes in the palaeopolyploid soybean.

Many plants have undergone whole genome duplication (WGD). However, how regulatory networks underlying a particular trait are reshaped in polyploids has not been experimentally investigated. Here we show that the regulatory pathways modulating seed oil content, which involve WRINKLED1 (WRI1), LEAFY COTYLEDON1 (LEC1), and LEC2 in Arabidopsis, have been modified in the palaeopolyploid soybean. Such modifications include functional reduction of GmWRI1b of the GmWRI1a/GmWRI1b homoeologous pair relevant to WRI1, complementary non-allelic dosage effects of the GmLEC1a/GmLEC1b homoeologous pair relevant to LEC1, pseudogenization of the singleton GmLEC2 relevant to LEC2, and the rise of the LEC2-like function of GmABI3b, contrasting to its homoeolog GmABI3a, which maintains the ABSCISIC ACID INSENSITIVE 3 (ABI3)-like function in modulating seed maturation and dormancy. The function of GmABI3b in modulating seed oil biosynthesis was fulfilled by direct binding to a RY (CATGCA) cis-regulatory element in the GmWRI1a promoter, which was absent in the GmWRI1b promoter, resulting in reduction of the GmWRI1b expression. Nevertheless, the three regulators each exhibited similar intensities of purifying selection to their respective duplicates since these pairs were formed by a WGD event that is proposed to have occurred approximately 13 million years ago (mya), suggesting that the differentiation in spatiotemporal expression between the duplicated genes is more likely to be the outcome of neutral variation in regulatory sequences. This study thus exemplifies the plasticity, dynamics, and novelty of regulatory networks mediated by WGD.

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

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

Developmental and environmental regulation of Aquaporin gene expression across Populus species: divergence or redundancy?

Aquaporins (AQPs) are membrane channels belonging to the major intrinsic proteins family and are known for their ability to facilitate water movement. While in Populus trichocarpa, AQP proteins form a large family encompassing fifty-five genes, most of the experimental work focused on a few genes or subfamilies. The current work was undertaken to develop a comprehensive picture of the whole AQP gene family in Populus species by delineating gene expression domain and distinguishing responsiveness to developmental and environmental cues. Since duplication events amplified the poplar AQP family, we addressed the question of expression redundancy between gene duplicates. On these purposes, we carried a meta-analysis of all publicly available Affymetrix experiments. Our in-silico strategy controlled for previously identified biases in cross-species transcriptomics, a necessary step for any comparative transcriptomics based on multispecies design chips. Three poplar AQPs were not supported by any expression data, even in a large collection of situations (abiotic and biotic constraints, temporal oscillations and mutants). The expression of 11 AQPs was never or poorly regulated whatever the wideness of their expression domain and their expression level. Our work highlighted that PtTIP1;4 was the most responsive gene of the AQP family. A high functional divergence between gene duplicates was detected across species and in response to tested cues, except for the root-expressed PtTIP2;3/PtTIP2;4 pair exhibiting 80% convergent responses. Our meta-analysis assessed key features of aquaporin expression which had remained hidden in single experiments, such as expression wideness, response specificity and genotype and environment interactions. By consolidating expression profiles using independent experimental series, we showed that the large expansion of AQP family in poplar was accompanied with a strong divergence of gene expression, even if some cases of functional redundancy could be suspected.

Why has butyrylcholinesterase been retained? Structural and functional diversification in a duplicated gene.

While acetylcholinesterase (EC 3.1.1.7) has a clearly defined role in neurotransmission, the functions of its sister enzyme butyrylcholinesterase (EC 3.1.1.8) are more obscure. Numerous mutations, many inactivating, are observed in the human butyrylcholinesterase gene, and the butyrylcholinesterase knockout mouse has an essentially normal phenotype, suggesting that the enzyme may be redundant. Yet the gene has survived for many millions of years since the duplication of an ancestral acetylcholinesterase early in vertebrate evolution. In this paper, we ask the questions: why has butyrylcholinesterase been retained, and why are inactivating mutations apparently tolerated? Butyrylcholinesterase has diverged both structurally and in terms of tissue and cellular expression patterns from acetylcholinesterase. Butyrylcholinesterase-like activity and enzymes have arisen a number of times in the animal kingdom, suggesting the usefulness of such enzymes. Analysis of the published literature suggests that butyrylcholinesterase has specific roles in detoxification as well as in neurotransmission, both in the brain, where it appears to control certain areas and functions, and in the neuromuscular junction, where its function appears to complement that of acetylcholinesterase. An analysis of the mutations in human butyrylcholinesterase and their relation to the enzyme's structure is shown. In conclusion, it appears that the structure of butyrylcholinesterase's catalytic apparatus is a compromise between the apparently conflicting selective demands of a more generalised detoxifier and the necessity for maintaining high catalytic efficiency. It is also possible that the tolerance of mutation in human butyrylcholinesterase is a consequence of the detoxification function. Butyrylcholinesterase appears to be a good example of a gene that has survived by subfunctionalisation.

Dual neofunctionalization of a rapidly evolving aquaporin-1 paralog resulted in constrained and relaxed traits controlling channel function during meiosis resumption in teleosts.

The preovulatory hydration of teleost oocytes is a unique process among vertebrates. The hydration mechanism is most pronounced in marine acanthomorph teleosts that spawn pelagic (floating) eggs; however, the molecular pathway for water influx remains poorly understood. Recently, we revealed that whole-genome duplication (WGD) resulted in teleosts harboring the largest repertoire of molecular water channels in the vertebrate lineage and that a duplicated aquaporin-1 paralog is implicated in the oocyte hydration process. However, the origin and function of the aquaporin-1 paralogs remain equivocal. By integrating the molecular phylogeny with synteny and structural analyses, we show here that the teleost aqp1aa and -1ab paralogs (previously annotated as aqp1a and -1b, respectively) arose by tandem duplication rather than WGD and that the Aqp1ab C-terminus is the most rapidly evolving subdomain within the vertebrate aquaporin superfamily. The functional role of Aqp1ab was investigated in Atlantic halibut, a marine acanthomorph teleost that spawns one of the largest pelagic eggs known. We demonstrate that Aqp1ab is required for full hydration of oocytes undergoing meiotic maturation. We further show that the rapid structural divergence of the C-terminal regulatory domain causes ex vivo loss of function of halibut Aqp1ab when expressed in amphibian oocytes but not in zebrafish or native oocytes. However, by using chimeric constructs of halibut Aqp1aa and -1ab and antisera specifically raised against the C-terminus of Aqp1ab, we found that this cytoplasmic domain regulates in vivo trafficking to the microvillar portion of the oocyte plasma membrane when intraoocytic osmotic pressure is at a maximum. Interestingly, by coinjecting polyA(+) mRNA from postvitellogenic halibut follicles, ex vivo intracellular trafficking of Aqp1ab is rescued in amphibian oocytes. These data reveal that the physiological role of Aqp1ab during meiosis resumption is conserved in teleosts, but the remarkable degeneracy of the cytoplasmic domain has resulted in alternative regulation of the trafficking mechanism.

Exposing the fitness contribution of duplicated genes.

Duplicate genes from the whole-genome duplication (WGD) in yeast are often dispensable--removing one copy has little or no phenotypic consequence. It is unknown, however, whether such dispensability reflects insignificance of the ancestral function or compensation from paralogs. Here, using precise competition-based measurements of the fitness cost of single and double deletions, we estimate the exposed fitness contribution of WGD duplicate genes in metabolism and bound the importance of their ancestral pre-duplication function. We find that the functional overlap between paralogs sufficiently explains the apparent dispensability of individual WGD genes. Furthermore, the lower bound on the fitness value of the ancestral function, which is estimated by the degree of synergistic epistasis, is at least as large as the average fitness cost of deleting single non-WGD genes. These results suggest that most metabolic functions encoded by WGD genes are important today and were also important at the time of duplication.

Preferential protection of protein interaction network hubs in yeast: evolved functionality of genetic redundancy.

The widely observed dispensability of duplicate genes is typically interpreted to suggest that a proportion of the duplicate pairs are at least partially redundant in their functions, thus allowing for compensatory affects. However, because redundancy is expected to be evolutionarily short lived, there is currently debate on both the proportion of redundant duplicates and their functional importance. Here, we examined these compensatory interactions by relying on a genome wide data analysis, followed by experiments and literature mining in yeast. Our data, thus, strongly suggest that compensated duplicates are not randomly distributed within the protein interaction network but are rather strategically allocated to the most highly connected proteins. This design is appealing because it suggests that many of the potentially vulnerable nodes that would otherwise be highly sensitive to mutations are often protected by redundancy. Furthermore, divergence analyses show that this association between redundancy and protein connectivity becomes even more significant among the ancient duplicates, suggesting that these functional overlaps have undergone purifying selection. Our results suggest an intriguing conclusion-although redundancy is typically transient on evolutionary time scales, it tends to be preserved among some of the central proteins in the cellular interaction network.

The degree of redundancy in metabolic genes is linked to mode of metabolism.

An understanding of the factors favoring the maintenance of duplicate genes in microbial genomes is essential for developing models of microbial evolution. A genome-scale flux-balance analysis of the metabolic network of Saccharomyces cerevisiae has suggested that gene duplications primarily provide increased enzyme dosage to enhance metabolic flux because the incidence of gene duplications in essential genes is no higher than that in nonessential genes. Here, we used genome-scale metabolic models to analyze the extent of genetic and biochemical redundancy in prokaryotes that are either specialists, with one major mode of energy generation, or generalists, which have multiple metabolic strategies for conservation of energy. Surprisingly, the results suggest that generalists, such as Escherichia coli and Bacillus subtilis, are similar to the eukaryotic generalist, S. cerevisiae, in having a low percentage (<10%) of essential genes and few duplications of these essential genes, whereas metabolic specialists, such as Geobacter sulfurreducens and Methanosarcina barkeri, have a high percentage (>30%) of essential genes and a high degree of genetic redundancy in these genes compared to nonessential genes. Furthermore, the specialist organisms appear to rely more on gene duplications rather than alternative-but-equivalent metabolic pathways to provide resilience to gene loss. Generalists rely more on alternative pathways. Thus, the concept that the role of gene duplications is to boost enzymatic flux rather than provide metabolic resilience may not be universal. Rather, the degree of gene duplication in microorganisms may be linked to mode of metabolism and environmental niche.

Differential expression of duplicated VAL-opsin genes in the developing zebrafish.

Non-visual opsins mediate various light-dependent physiological events. Our previous search for non-visual opsin genes in zebrafish led to the discovery of VAL-opsin (VAL-opsinA) in deep brain cells and retinal horizontal cells of the adult fish. In this study, we report the identification and characterization of its duplicated gene, VAL-opsinB, in zebrafish. A molecular phylogenetic analysis indicates that VAL-opsinB is orthologous to a previously reported salmon gene and that the duplication of the VAL-opsin gene occurred in the teleost lineage. The recombinant protein of zebrafish VAL-opsinB forms a green-sensitive photopigment when reconstituted with 11-cis-retinal. VAL-opsinB expression was detected in a limited number of cells of the brain and the eye, and the expression pattern is distinct from that of the VAL-opsinA gene. Such a differential expression pattern suggests that VAL-opsinA and VAL-opsinB are involved in different physiological events in zebrafish.

The evolutionary fate of recently duplicated retrogenes in mice.

Inferences about the evolutionary impact of gene duplications often rely on the analysis of their long-term outcome. The fate of the majority of them must, however, be decided shortly after duplication. Here we analysed the evolutionary pattern of 10 mouse genes very recently duplicated by retrotransposition, by sequencing the retroposed copy in five to 10 closely related mouse species. In all cases the retroposed copy experienced accelerated nonsynonymous evolution whereas the divergence pattern of the source copy appeared unaffected by the duplication, consistent with the neofunctionalization model. The analysis further revealed that most retrogenes, including pseudogenes, did not experience a period of relaxed neutral evolution, but have been submitted to purifying selection ever since their retroposition. We propose that these duplicates play a biochemical role but are not indispensable. Purifying selection prevents them from acquiring a negative role until they are lost or silenced. This period of unnecessary redundancy could in rare cases give the time for new functions to evolve.

Transcriptional reprogramming and backup between duplicate genes: is it a genomewide phenomenon?

Deleting a duplicate gene often results in a less severe phenotype than deleting a singleton gene, a phenomenon commonly attributed to functional compensation among duplicates. However, duplicate genes rapidly diverge in expression patterns after duplication, making functional compensation less probable for ancient duplicates. Case studies suggested that a gene may provide compensation by altering its expression upon removal of its duplicate copy. On the basis of this observation and a genomic analysis, it was recently proposed that transcriptional reprogramming and backup among duplicates is a genomewide phenomenon in the yeast Saccharomyces cerevisiae. Here we reanalyze the yeast data and show that the high dispensability of duplicate genes with low expression similarity is a consequence of expression similarity and gene dispensability, each being correlated with a third factor, the number of protein interactions per gene. There is little evidence supporting widespread functional compensation of divergently expressed duplicate genes by transcriptional reprogramming.

Duplicated Pax6 genes in Glomeris marginata (Myriapoda: Diplopoda), an arthropod with simple lateral eyes.

Composite (facetted) eyes comprised by several units, termed ommatidia, are an ancestral feature in the arthropods. Some arthropods, however, do not possess composite eyes, obviously by secondary reduction. Reductions on the level of conserved eye developmental genes are one possibility to reduce the visual system. The genes of the Pax6 family have been shown to be key regulators of visual system development in a wide variety of animals. Reduction of Pax6 expression may therefore be expected in a species with reduced eyes. Here I have investigated the myriapod Glomeris marginata that displays very simple eyes. Glomeris, however, possesses two Pax6 genes that, based on their sequence, are similar to Drosophila eyeless (ey) and twin of eyeless (toy), respectively. Both genes are highly expressed in the optic lobes and the ventral nerve cord of developing embryos. Furthermore, homologs of other high-ranking eye developmental genes like hedgehog, decapentaplegic, dachshund, and homothorax are expressed in the optic lobes. This indicates that eye reduction in Glomeris is not realized at the level of the Pax6 genes or other genes on the upper levels of the eye development network. I suggest instead that the simple eyes of Glomeris are the product of changes at a much lower level in the network, probably at the level of genes directly regulating ommatidia development or ommatidia number and arrangement.

Evolution in action: following function in duplicated floral homeotic genes.

Gene duplication plays a fundamental role in evolution by providing the genetic material from which novel functions can arise. Newly duplicated genes can be maintained by subfunctionalization (the duplicated genes perform different aspects of the original gene's function) and/or neofunctionalization (one of the genes acquires a novel function). PLENA in Antirrhinum and AGAMOUS in Arabidopsis are the canonical C-function genes that are essential for the specification of reproductive organs. These functionally equivalent genes encode closely related homeotic MADS-box transcription factors. Using genome synteny, we confirm phylogenetic analyses showing that PLENA and AGAMOUS are nonorthologous genes derived from a duplication in a common ancestor. Their respective orthologs, SHATTERPROOF in Arabidopsis and FARINELLI in Antirrhinum, have undergone independent subfunctionalization via changes in regulation and protein function. Surprisingly, the functional divergence between PLENA and FARINELLI, is morphologically manifest in both transgenic Antirrhinum and Arabidopsis. This provides a clear illustration of a random evolutionary trajectory for gene functions after a duplication event. Different members of a duplicated gene pair have retained the primary homeotic functions in different lineages, illustrating the role of chance in evolution. The differential ability of the Antirrhinum genes to promote male or female development provides a striking example of subfunctionalization at the protein level.

Duplicated genes evolve slower than singletons despite the initial rate increase.

BACKGROUND: Gene duplication is an important mechanism that can lead to the emergence of new functions during evolution. The impact of duplication on the mode of gene evolution has been the subject of several theoretical and empirical comparative-genomic studies. It has been shown that, shortly after the duplication, genes seem to experience a considerable relaxation of purifying selection. RESULTS: Here we demonstrate two opposite effects of gene duplication on evolutionary rates. Sequence comparisons between paralogs show that, in accord with previous observations, a substantial acceleration in the evolution of paralogs occurs after duplication, presumably due to relaxation of purifying selection. The effect of gene duplication on evolutionary rate was also assessed by sequence comparison between orthologs that have paralogs (duplicates) and those that do not (singletons). It is shown that, in eukaryotes, duplicates, on average, evolve significantly slower than singletons. Eukaryotic ortholog evolutionary rates for duplicates are also negatively correlated with the number of paralogs per gene and the strength of selection between paralogs. A tally of annotated gene functions shows that duplicates tend to be enriched for proteins with known functions, particularly those involved in signaling and related cellular processes; by contrast, singletons include an over-abundance of poorly characterized proteins. CONCLUSIONS: These results suggest that whether or not a gene duplicate is retained by selection depends critically on the pre-existing functional utility of the protein encoded by the ancestral singleton. Duplicates of genes of a higher biological import, which are subject to strong functional constraints on the sequence, are retained relatively more often. Thus, the evolutionary trajectory of duplicated genes appears to be determined by two opposing trends, namely, the post-duplication rate acceleration and the generally slow evolutionary rate owing to the high level of functional constraints.

Differential expression of duplicated genes for brain-type fatty acid-binding proteins (fabp7a and fabp7b) during early development of the CNS in zebrafish (Danio rerio).

A gene for the zebrafish brain-type fatty acid-binding protein (fabp7b) was identified and its structure defined. The zebrafish fabp7b gene spans 1479 bp and consists of four exons encoding 24, 58, 34 and 16 amino acids, respectively, which is identical to the structure of the fabp7a gene previously described. The complete fabp7b cDNA was isolated by 5' and 3' RACE and its nucleotide sequence determined. The deduced amino acid sequence of FABP7B encoded by the zebrafish fabp7b gene shares 82% identity with that of FABP7A encoded by the zebrafish fabp7a gene. A single transcription start site for the fabp7b gene was mapped by 5' RNA ligase-mediated RACE. Phylogenetic analysis indicated that the duplication of the fabp7 genes occurred in the fish lineage after their divergence from mammals. The zebrafish fabp7b gene was assigned to linkage group 20 by radiation hybrid mapping. Reverse transcription-polymerase chain reaction detected fabp7b transcripts in the same adult tissues as fabp7a transcripts. In the brain, levels of fabp7b transcripts were lower than fabp7a transcripts. Whole-mount in situ hybridization showed that the zebrafish fabp7a transcripts were distributed in the early developing central nervous system. In addition to being expressed in the developing brain and retina, zebrafish fabp7b mRNA was also detected in the swim bladder and pharynx during the embryonic to larval transitory phase.

Divergence in the spatial pattern of gene expression between human duplicate genes.

Microarray gene expression data provide a wealth of information for elucidating the mode and tempo of molecular evolution. In the present study,we analyze the spatial expression pattern of human duplicate gene pairs by using oligonucleotide microarray data,and study the relationship between coding sequence divergence and expression divergence. First,we find a strong positive correlation between the proportion of duplicate gene pairs with divergent expression (as presence or absence of expression in a tissue) and both synonymous (K(S)) and nonsynonymous divergence (K(A)). The divergence of gene expression between human duplicate genes is rapid, probably faster than that between yeast duplicates in terms of generations. Second,we compute the correlation coefficient (R) between the expression levels of duplicate genes in different tissues and find a significant negative correlation between R and K(S). There is also a negative correlation between R and K(A), when K(A)

Gene duplications: the gradual evolution of functional divergence.

Budding yeast provides a useful resource for studies of gene function. A new analysis of the fitness effects of deletion mutations in budding yeast reveals that genes that have duplicates create lower fitness losses when inactivated than do genes that are singletons.

GenomeHistory: a software tool and its application to fully sequenced genomes.

We present a publicly available software tool (http://www.unm.edu/~compbio/software/GenomeHistory) that identifies all pairs of duplicate genes in a genome and then determines the degree of synonymous and non-synonymous divergence between each duplicate pair. Using this tool, we analyze the relations between (i) gene function and the propensity of a gene to duplicate and (ii) the number of genes in a gene family and the family's rate of sequence evolution. We do so for the complete genomes of four eukaryotes (fission and budding yeast, fruit fly and nematode) and one prokaryote (Escherichia coli). For some classes of genes we observe a strong relationship between gene function and a gene's propensity to undergo duplication. Most notably, ribosomal genes and transcription factors appear less likely to undergo gene duplication than other genes. In both fission and budding yeast, we see a strong positive correlation between the selective constraint on a gene and the size of the gene family of which this gene is a member. In contrast, a weakly negative such correlation is seen in multicellular eukaryotes.