The evolutionary demography of duplicate genes

Although gene duplication has generally been viewed as a necessary source of material for the origin of evolutionary novelties, the rates of origin, loss, and preservation of gene duplicates are not well understood. Applying steady-state demographic techniques to the age distributions of duplicate genes censused in seven completely sequenced genomes, we estimate the average rate of duplication of a eukaryotic gene to be on the order of 0.01/gene/million years, which is of the same order of magnitude as the mutation rate per nucleotide site. However, the average half-life of duplicate genes is relatively small, on the order of 4.0 million years. Significant interspecific variation in these rates appears to be responsible for differences in species-specific genome sizes that arise as a consequence of a quasi-equilibrium birth-death process. Most duplicated genes experience a brief period of relaxed selection early in their history and a minority exhibit the signature of directional selection, but those that survive more than a few million years eventually experience strong purifying selection. Thus, although most theoretical work on the gene-duplication process has focused on issues related to adaptive evolution, the origin of a new function appears to be a very rare fate for a duplicate gene. A more significant role of the duplication process may be the generation of microchromosomal rearrangements through reciprocal silencing of alternative copies, which can lead to the passive origin of post-zygotic reproductive barriers in descendant lineages of incipient species.     

Evolutionary dynamics of duplicated genes in plants

Gene duplication, arising from region-specific duplication or genome-wide polyploidization, is a prominent feature in plant genome evolution. Understanding the mechanisms generating duplicate gene copies and the subsequent dynamics among gene duplicates is vital because these investigations shed light on regional and genome-wide aspects of evolutionary forces shaping intra- and interspecific genome contents, evolutionary relationships, and interactions. This review discusses recent gene duplication analyses in plants, focusing on the molecular and evolutionary dynamics occurring at three different timescales following duplication: (1). initial establishment and persistence of cytotypes, (2). interactions among duplicate gene copies, and (3). longer term differentiation between duplicated genes. These relative time points are presented in terms of their potential adaptive significance and impact on plant evolutionary genomics research.     

Expression and evolutionary conservation of nanos-related genes in Hydra

The Drosophila gene nanos encodes two particular zinc finger motifs which are also found in germline-associated factors from nematodes to vertebrates. We cloned two nanos (nos)-related genes, Cnnos1 and Cnnos2 from Hydra magnipapillata. Using whole-mount in situ hybridization, the expression of Cnnos1 and Cnnos2 was examined. Cnnos1 was specifically expressed in multipotent stem cells and germline cells, but not in somatic cells. Cnnos2 was weakly expressed in germline cells and more specifically in the endoderm of the hypostome where it appears to be involved in head morphogenesis. In addition to structural conservation in the zinc finger domain of nanos-related genes, functional conservation of Cnnos1 was also demonstrated by the finding that a Cnnos1 transgene can partially rescue the nos ^RC phenotype that is defective in the egg production of Drosophila. Thus, the function of nanos-related genes in the germline appears to be well conserved from primitive to highly evolved metazoans. Received: 28 April 2000 / Accepted: 1 July 2000     

The evolutionary relationship of avian and mammalian myosin heavy-chain genes

Summary Sequence comparisons of avian and mammalian skeletal and cardiac myosin heavy-chain isoforms are used to examine the evolutionary relationships of sarcomeric myosin multigene families. Mammalian fast-myosin heavy-chain isoforms forms from different species, with comparable developmental expression, are more similar to each other than they are to other fast isoforms within the same genome. In contrast, the developmentally regulated chicken fast isoforms are more similar to each other than they are to myosin heavy-chain isoforms in other species. Extensive regions of nucleotide identity among the chicken fast myosin heavy chains and in the mouse and rat α- and β-cardiac myosin heavy-chain sequences suggest that geneconversion-like mechanisms have played a major role in the concerted evolution of these gene families. We also conclude that the chicken fast myosin heavy-chain multigene family has undergone recent expansion subsequent to the divergence of birds and mammals and that both the developmental regulation and the specialization of myosin isoforms have likely developed independently in birds and mammals.     

Preservation of duplicate genes by originalization

Neofunctionalization, subfunctionalization and increasing gene dosage were proposed to be the possible ways to explain duplicate-gene preservation in previous studies. However, in some natural populations, such as yeast Saccharomyces cerevisiae, a considerable proportion of the duplicate genes originated from ancient whole genomic duplication (WGD) is preserved till now, which cannot be sufficiently explained by these mechanisms. In this article, we present another possible way to explain this conundrum—originalization, by which duplicate genes are both preserved intact at a high frequency in the population under only purifying selection. With approximate equal rates of mutation at the two duplicated loci, analytical, numerical and simulation results consistently show that the mean time to nonfunctionalization for unlinked haploinsufficient gene duplication might become markedly prolonged, which results from originalization. These theoretical results imply that originalization might be an alternative effective and temporary way of preserving duplicate genes.     

Evolution of vertebrate central nervous system is accompanied by novel expression changes of duplicate genes

The evolution of the central nervous system (CNS) is one of the most striking changes during the transition from invertebrates to vertebrates. As a major source of genetic novelties, gene duplication might play an important role in the functional innovation of vertebrate CNS. In this study, we focused on a group of CNS-biased genes that duplicated during early vertebrate evolution. We investigated the tempo-spatial expression patterns of 33 duplicate gene families and their orthologs during the embryonic development of the vertebrate Xenopus laevis and the cephalochordate Brachiostoma belcheri. Almost all the identified duplicate genes are differentially expressed in the CNS in Xenopus embryos, and more than 50% and 30% duplicate genes are expressed in the telencephalon and mid-hindbrain boundary, respectively, which are mostly considered as two innovations in the vertebrate CNS. Interestingly, more than 50% of the amphioxus orthologs do not show apparent expression in the CNS in amphioxus embryos as detected by in situ hybridization, indicating that some of the vertebrate CNS-biased duplicate genes might arise from non-CNS genes in invertebrates. Our data accentuate the functional contribution of gene duplication in the CNS evolution of vertebrate and uncover an invertebrate non-CNS history for some vertebrate CNS-biased duplicate genes.     

Divergence pattern of duplicate genes in protein-protein interactions follows the power law

The impact of the biological network structures on the divergence between the two copies of one duplicate gene pair involved in the networks has not been documented on a genome scale. Having analyzed the most recently updated Database of Interacting Proteins (DIP) by incorporating the information for duplicate genes of the same age in yeast, we find that there was a highly significantly positive correlation between the level of connectivity of ancient genes and the number of shared partners of their duplicates in the protein-protein interaction networks. This suggests that duplicate genes with a low ancestral connectivity tend to provide raw materials for functional novelty, whereas those duplicate genes with a high ancestral connectivity tend to create functional redundancy for a genome during the same evolutionary period. Moreover, the difference in the number of partners between two copies of a duplicate pair was found to follow a power-law distribution. This suggests that loss and gain of interacting partners for most duplicate genes with a lower level of ancestral connectivity is largely symmetrical, whereas the "hub duplicate genes" with a higher level of ancient connectivity display an asymmetrical divergence pattern in protein-protein interactions. Thus, it is clear that the protein-protein interaction network structures affect the divergence pattern of duplicate genes. Our findings also provide insights into the origin and development of biological networks.     

Faster evolution of Z-linked duplicate genes in chicken

It has been shown that duplicate genes on the X chromosome evolve much faster than duplicate genes on autosomes in Drosophila melanogaster.However,whether this phenomenon is general and can be applied to other species is not known.Here we examined this issue in chicken that have heterogametic females(females have ZW sex chromosome).We compared sequence divergence of duplicate genes on the Z chromosome with those on autosomes.We found that duplications on the Z chromosome indeed evolved faster than those on autosomes and show distinct patterns of molecular evolution from autosomal duplications.Examination of the expression of duplicate genes revealed an enrichment of duplications on the Z chromosome having male-biased expression and an enrichment of duplications on the autosomes having female-biased expression.These results suggest an evolutionary trend of the recruitment of duplicate genes towards reproduction-specific function.The faster evolution of duplications on Z than on the autosomes is most likely contributed by the selective forces driving the fixation of adaptive mutations on Z.Therefore,the common phenomena observed in both flies and chicken suggest that duplicate genes on sex chromosomes have distinct dynamics and are more influenced by natural selection than antosomal duplications,regardless of the kind of sex determination systems.     

Polyploidy, evolutionary opportunity, and crop adaptation

The finding that even the smallest of plant genomes has incurred multiple genome-wide chromatin duplication events, some of which may predate the origins of the angiosperms and therefore shape all of flowering plant biology, adds new importance to the molecular analysis of polyploidization/diploidization cycles and their phenotypic consequences. Early clues as to the possible phenotypic consequences of polyploidy derive from recent QTL mapping efforts in a number of diverse crop plants of recent and well-defined polyploid origins. A small sampling examples of the role(s) of polyploidy in conferring crop adaptation from human needs include examples of (1) dosage effects of multiple alleles in autopolyploids, and (2) intergenomic heterosis conferring novel traits or transgressive levels of existing traits, associated with merging divergent genomes in a common allopolyploid nucleus. A particularly interesting manifestation of #2 is the evolution of complementary alleles at corresponding (homoeologous) loci in divergent polyploid taxa derived from a common ancestor. Burgeoning genomic data for both botanical models and major crops offer new avenues for investigation of the molecular and phenotypic consequences of polyploidy, promising new insights into the role of this important process in the evolution of botanical diversity.     

Genome-wide comparative analysis of NBS-encoding genes between Brassica species and Arabidopsis thaliana

Background

Plant disease resistance (R) genes with the nucleotide binding site (NBS) play an important role in offering resistance to pathogens. The availability of complete genome sequences of Brassica oleracea and Brassica rapa provides an important opportunity for researchers to identify and characterize NBS-encoding R genes in Brassica species and to compare with analogues in Arabidopsis thaliana based on a comparative genomics approach. However, little is known about the evolutionary fate of NBS-encoding genes in the Brassica lineage after split from A. thaliana.

Results

Here we present genome-wide analysis of NBS-encoding genes in B. oleracea, B. rapa and A. thaliana. Through the employment of HMM search and manual curation, we identified 157, 206 and 167 NBS-encoding genes in B. oleracea, B. rapa and A. thaliana genomes, respectively. Phylogenetic analysis among 3 species classified NBS-encoding genes into 6 subgroups. Tandem duplication and whole genome triplication (WGT) analyses revealed that after WGT of the Brassica ancestor, NBS-encoding homologous gene pairs on triplicated regions in Brassica ancestor were deleted or lost quickly, but NBS-encoding genes in Brassica species experienced species-specific gene amplification by tandem duplication after divergence of B. rapa and B. oleracea. Expression profiling of NBS-encoding orthologous gene pairs indicated the differential expression pattern of retained orthologous gene copies in B. oleracea and B. rapa. Furthermore, evolutionary analysis of CNL type NBS-encoding orthologous gene pairs among 3 species suggested that orthologous genes in B. rapa species have undergone stronger negative selection than those in B .oleracea species. But for TNL type, there are no significant differences in the orthologous gene pairs between the two species.

Conclusion

This study is first identification and characterization of NBS-encoding genes in B. rapa and B. oleracea based on whole genome sequences. Through tandem duplication and whole genome triplication analysis in B. oleracea, B. rapa and A. thaliana genomes, our study provides insight into the evolutionary history of NBS-encoding genes after divergence of A. thaliana and the Brassica lineage. These results together with expression pattern analysis of NBS-encoding orthologous genes provide useful resource for functional characterization of these genes and genetic improvement of relevant crops.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-3) contains supplementary material, which is available to authorized users.     

Anthocyanin biosynthetic genes in Brassica rapa

Background

Anthocyanins are a group of flavonoid compounds. As a group of important secondary metabolites, they perform several key biological functions in plants. Anthocyanins also play beneficial health roles as potentially protective factors against cancer and heart disease. To elucidate the anthocyanin biosynthetic pathway in Brassica rapa, we conducted comparative genomic analyses between Arabidopsis thaliana and B. rapa on a genome-wide level.

Results

In total, we identified 73 genes in B. rapa as orthologs of 41 anthocyanin biosynthetic genes in A. thaliana. In B. rapa, the anthocyanin biosynthetic genes (ABGs) have expanded and most genes exist in more than one copy. The anthocyanin biosynthetic structural genes have expanded through whole genome and tandem duplication in B. rapa. More structural genes located upstream of the anthocyanin biosynthetic pathway have been retained than downstream. More negative regulatory genes are retained in the anthocyanin biosynthesis regulatory system of B. rapa.

Conclusions

These results will promote an understanding of the genetic mechanism of anthocyanin biosynthesis, as well as help the improvement of the nutritional quality of B. rapa through the breeding of high anthocyanin content varieties.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-426) contains supplementary material, which is available to authorized users.     

Mass spectrometric analysis of expression of ATPase subunits encoded by duplicated genes in the 19S regulatory particle of rice 26S proteasome

The 26S proteasome consisting of a 20S proteasome and a pair of 19S regulatory particles (RP) plays important roles in degradation of the ubiquitinated protein in eukaryotic cells. The RP consists of six different ATPase subunits and, at least, 11 non-ATPase subunits. In rice, we previously identified duplicated genes encoding four ATPase subunits, OsRpt1, OsRpt2, OsRpt4, and OsRpt5. In this study, the genomic sequences of all rice ATPase subunits were identified from the rice genome database and the genomic structure of ATPase subunit genes was determined. The rice RP was purified, and the ATPase subunit isoforms encoded by three pairs of duplicated genes, OsRpt2a/OsRpt2b, OsRpt4a/OsRpt4b, and OsRpt5a/OsRpt5b, were identified in RP by using electrospray ionization quadrupole time-of-flight mass spectrometry. The relative amounts and the expression patterns of these ATPase subunit isoforms in the bran were found to be different from those of the callus, suggesting the presence of multiform 19S regulatory particles engaged in the tissue-specific protein metabolism.     

Characterization of genes with tissue-specific differential expression patterns in Populus

Like many plants, Populus has an evolutionary history in which several, both recent and more ancient, genome duplication events have occurred and, therefore, constitutes an excellent model system for studying the functional evolution of genes. In the present study, we have focused on the properties of genes with tissue-specific differential expression patterns in poplar. We identified the genes by analyzing digital expression profiles derived by mapping 90,000+ expressed sequence tags (ESTs) from 18 sources to the predicted genes of Populus. Our sequence analysis suggests that tissue-specific differentially expressed genes have less diverged paralogs than average, indicating that gene duplication events is an important event in the pathway leading to this type of expression pattern. The functional analysis showed that genes coding for proteins involved in processes of functional importance for the specific tissue(s) in which they are expressed and genes coding for regulatory or responsive proteins are most common among the differentially expressed genes, demonstrating that the expression differentiation process is under strong selective pressure. Thus, our data supports a model where gene duplication followed by gene specialization or expansion of the regulatory and responsive networks leads to tissue-specific differential expression patterns. We have also searched for clustering of genes with similar expression pattern into gene-expression neighborhoods within the Populus genome. However, we could not detect any major clustering among the analyzed genes with highly specific expression patterns. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Genome-wide identification of NBS genes in<Emphasis Type="Italic"> japonica</Emphasis> rice reveals significant expansion of divergent non-TIR NBS-LRR genes

A complete set of candidate disease resistance ( R) genes encoding nucleotide-binding sites (NBSs) was identified in the genome sequence of japonica rice ( Oryza sativa L. var. Nipponbare). These putative R genes were characterized with respect to structural diversity, phylogenetic relationships and chromosomal distribution, and compared with those in Arabidopsis thaliana. We found 535 NBS-coding sequences, including 480 non-TIR (Toll/IL-1 receptor) NBS-LRR (Leucine Rich Repeat) genes. TIR NBS-LRR genes, which are common in A. thaliana, have not been identified in the rice genome. The number of non-TIR NBS-LRR genes in rice is 8.7 times higher than that in A. thaliana, and they account for about 1% of all of predicted ORFs in the rice genome. Some 76% of the NBS genes were located in 44 gene clusters or in 57 tandem arrays, and 16 apparent gene duplications were detected in these regions. Phylogenetic analyses based both NBS and N-terminal regions classified the genes into about 200 groups, but no deep clades were detected, in contrast to the two distinct clusters found in A. thaliana. The structural and genetic diversity that exists among NBS-LRR proteins in rice is remarkable, and suggests that diversifying selection has played an important role in the evolution of R genes in this agronomically important species. (Supplemental material is available online at .)Communicated by R. HagemannThe first three authors contributed equally to this work     

Structure,evolution and expression of a second subfamily of protein phosphatase 2A catalytic subunit genes in the rice plant (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Protein phosphatase 2A (PP2A) is one of the major serine/threonine protein phosphatases in the cell and plays a variety of regulatory roles in metabolism and signal transduction. Previously, we described the structure and expression of two genes encoding PP2A catalytic subunits (PP2Ac)—OsPP2A-1 and OsPP2A-3—in the rice plant (Yu et al. 2003). Here, we report the isolation and characterisation of a second structurally distinguishable PP2Ac subfamily comprised of three additional isogenes, OsPP2A-2, OsPP2A-4 (each containing ten introns) and OsPP2A-5 (which contains nine introns). Northern blot analysis demonstrated that the three isogenes are ubiquitously expressed in all rice tissues during plant development, and differentially expressed in response to high salinity and the combined stresses of drought and heat. Phylogenetic analyses indicated that the two PP2Ac subfamilies are descended from two ancient lineages, which derived from gene duplications that occurred after the monocotyledon–dicotyledon split. In the second subfamily, it is proposed that two duplication events were involved; in which, the initial duplication of a ten-intron primordial gene yielded OsPP2A-2 and the progenitor of OsPP2A-4 and OsPP2A-5. The OsPP2A-4/OsPP2A-5 progenitor, in turn, underwent a second duplication event, resulting in the present day OsPP2A-4 and OsPP2A-5. It is proposed that loss of the 5′-most intron from OsPP2A-5 occurred after these two duplication events.     

Molecular evolution of duplicate copies of genes encoding cytosolic glutamine synthetase in Pisum sativum

Here, we describe two nearly identical expressed genes for cytosolic glutamine synthetase (GS3A and GS3B) in Pisum sativum L. RFLP mapping data indicates that the GS3A and GS3B genes are separate loci located on different chromosomes. DNA sequencing of the GS3A and GS3B genes revealed that the coding regions are 99% identical with only simple nucleotide substitutions resulting in three amino acid differences. Surprisingly, the non-coding regions (5 non-coding leader, the 11 introns, and 3 non-coding tail) all showed a high degree of identity (96%). In these non-coding regions, 25% of the observed differences between the GS3A and GS3B genes were deletions or duplications. The single difference in the 3 non-coding regions of the GS3A and GS3B genes was a 25 bp duplication of an AU-rich element in the GS3B gene. As the GS3B mRNA accumulates to lower levels than the GS3A gene, we tested whether this sequence which resembles an mRNA instability determinant functioned as such in the context of the GS mRNA. Using the GS3B 3 tail as part of a chimeric gene in transgenic plants, we showed that this AU-rich sequence has little effect on transgene mRNA levels. To determine whether the GS3A/GS3B genes represent a recent duplication, we examined GS3-like genes in genomic DNA of ancient relatives of P. sativum. We observed that several members of the Viceae each contain two genomic DNA fragments homologous to the GS3B gene, suggesting that this is an ancient duplication event. Gene conversion has been invoked as a possible mechanism for maintaining the high level of nucleotide similarity found between the GS3A and GS3B genes. Possible evolutionary reasons for the maintenance of these twin GS genes in pea, and the general duplication of genes for cytosolic GS in all plant species are discussed.     

Two patterns of genome organization in mammals: the chromosomal distribution of duplicate genes in human and mouse

Gene duplication occurs repeatedly in the evolution of genomes, and the rearrangement of genomic segments has also occurred repeatedly over the evolution of eukaryotes. We studied the interaction of these two factors in mammalian evolution by comparing the chromosomal distribution of multigene families in human and mouse. In both species, gene families tended to be confined to a single chromosome to a greater extent than expected by chance. The average number of families shared between chromosomes was nearly 60% higher in mouse than in human, and human chromosomes rarely shared large numbers of gene families with more than one or two other chromosomes, whereas mouse chromosomes frequently did so. A higher proportion of duplicate gene pairs on the same chromosome originated from recent duplications in human than in mouse, whereas a higher proportion of duplicate gene pairs on separate chromosomes arose from ancient duplications in human than in mouse. These observations are most easily explained by the hypotheses that (1) most gene duplications arise in tandem and are subsequently separated by segmental rearrangement events, and (2) that the process of segmental rearrangement has occurred at a higher rate in the lineage of mouse than in that of human.     

Gene loss and evolutionary rates following whole-genome duplication in teleost fishes

Teleost fishes provide the first unambiguous support for ancient whole-genome duplication in an animal lineage. Studies in yeast or plants have shown that the effects of such duplications can be mediated by a complex pattern of gene retention and changes in evolutionary pressure. To explore such patterns in fishes, we have determined by phylogenetic analysis the evolutionary origin of 675 Tetraodon duplicated genes assigned to chromosomes, using additional data from other species of actinopterygian fishes. The subset of genes, which was retained in double after the genome duplication, is enriched in development, signaling, behavior, and regulation functional categories. The evolutionary rate of duplicate fish genes appears to be determined by 3 forces: 1) fish proteins evolve faster than mammalian orthologs; 2) the genes kept in double after genome duplication represent the subset under strongest purifying selection; and 3) following duplication, there is an asymmetric acceleration of evolutionary rate in one of the paralogs. These results show that similar mechanisms are at work in fishes as in yeast or plants and provide a framework for future investigation of the consequences of duplication in fishes and other animals.