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) 相似文献   

The KA/KS Ratio Test for Assessing the Protein-Coding Potential of Genomic Regions: An Empirical and Simulation Study

Comparative genomics is a simple, powerful way to increase the accuracy of gene prediction. In this study, we show the utility of a simple test for the identification of protein-coding exons using human/mouse sequence comparisons. The test takes advantage of the fact that in the vast majority of coding regions, synonymous substitutions (K(S)) occur much more frequently than nonsynonymous ones (K(A)) and uses the K(A)/K(S) ratio as the criterion. We show the following: (1) most of the human and mouse exons are sufficiently long and have a suitable degree of sequence divergence for the test to perform reliably; (2) the test is suited for the identification of long exons and single exon genes, which are difficult to predict by current methods; (3) the test has a false-negative rate, lower than most of current gene prediction methods and a false-positive rate lower than all current methods; (4) the test has been automated and can be used in combination with other existing gene-prediction methods.     

Origin and evolution of new exons in rodents

Gene number difference among organisms demonstrates that new gene origination is a fundamental biological process in evolution. Exon shuffling has been universally observed in the formation of new genes. Yet to be learned are the ways new exons originate and evolve, and how often new exons appear. To address these questions, we identified 2695 newly evolved exons in the mouse and rat by comparing the expressed sequences of 12,419 orthologous genes between human and mouse, using 743,856 pig ESTs as the outgroup. The new exon origination rate is about 2.71 x 10(-3) per gene per million years. These new exons have markedly accelerated rates both of nonsynonymous substitutions and of insertions/deletions (indels). A much higher proportion of new exons have K(a)/K(s) ratios >1 (where K(a) is the nonsynonymous substitution rate and K(s) is the synonymous substitution rate) than do the old exons shared by human and mouse, implying a role of positive selection in the rapid evolution. The majority of these new exons have sequences unique in the genome, suggesting that most new exons might originate through "exonization" of intronic sequences. Most of the new exons appear to be alternative exons that are expressed at low levels.     

A burst of protein sequence evolution and a prolonged period of asymmetric evolution follow gene duplication in yeast

It is widely accepted that newly arisen duplicate gene pairs experience an altered selective regime that is often manifested as an increase in the rate of protein sequence evolution. Many details about the nature of the rate acceleration remain unknown, however, including its typical magnitude and duration, and whether it applies to both gene copies or just one. We provide initial answers to these questions by comparing the rate of protein sequence evolution among eight yeast species, between a large set of duplicate gene pairs that were created by a whole-genome duplication (WGD) and a set of genes that were returned to single-copy after this event. Importantly, we use a new method that takes into account the tendency for slowly evolving genes to be retained preferentially in duplicate. We show that, on average, proteins encoded by duplicate gene pairs evolved at least three times faster immediately after the WGD than single-copy genes to which they behave identically in non-WGD lineages. Although the high rate in duplicated genes subsequently declined rapidly, it has not yet returned to the typical rate for single-copy genes. In addition, we show that although duplicate gene pairs often have highly asymmetric rates of evolution, even the slower members of pairs show evidence of a burst of protein sequence evolution immediately after duplication. We discuss the contribution of neofunctionalization to duplicate gene preservation and propose that a form of subfunctionalization mediated by coding region activity-reducing mutations is likely to have played an important role.     

Molecular evolution of avian reovirus: evidence for genetic diversity and reassortment of the S-class genome segments and multiple cocirculating lineages

Nucleotide sequences of the S-class genome segments of 17 field-isolates and vaccine strains of avian reovirus (ARV) isolated over a 23-year period from different hosts, pathotypes, and geographic locations were examined and analyzed to define phylogenetic profiles and evolutionary mechanism. The S1 genome segment showed noticeably higher divergence than the other S-class genes. The sigma C-encoding gene has evolved into six distinct lineages. In contrast, the other S-class genes showed less divergence than that of the sigma C-encoding gene and have evolved into two to three major distinct lineages, respectively. Comparative sequence analysis provided evidence indicating extensive sequence divergence between ARV and other orthoreoviruses. The evolutionary trees of each gene were distinct, suggesting that these genes evolve in an independent manner. Furthermore, variable topologies were the result of frequent genetic reassortment among multiple cocirculating lineages. Results showed genetic diversity correlated more closely with date of isolation and geographic sites than with host species and pathotypes. This is the first evidence demonstrating genetic variability among circulating ARVs through a combination of evolutionary mechanisms involving multiple cocirculating lineages and genetic reassortment. The evolutionary rates and patterns of base substitutions were examined. The evolutionary rate for the sigma C-encoding gene and sigma C protein was higher than for the other S-class genes and other family of viruses. With the exception of the sigma C-encoding gene, which nonsynonymous substitutions predominate over synonymous, the evolutionary process of the other S-class genes can be explained by the neutral theory of molecular evolution. Results revealed that synonymous substitutions predominate over nonsynonymous in the S-class genes, even though genetic diversity and substitution rates vary among the viruses.     

Genomic regionality in rates of evolution is not explained by clustering of genes of comparable expression profile

In mammalian genomes, linked genes show similar rates of evolution, both at fourfold degenerate synonymous sites (K4) and at nonsynonymous sites (KA). Although it has been suggested that the local similarity in the synonymous substitution rate is an artifact caused by the inclusion of disparately evolving gene pairs, we demonstrate here that this is not the case: after removal of disparately evolving genes, both (1) linked genes and (2) introns from the same gene have more similar silent substitution rates than expected by chance. What causes the local similarity in both synonymous and nonsynonymous substitution rates? One class of hypotheses argues that both may be related to the observed clustering of genes of comparable expression profile. We investigate these hypotheses using substitution rates from both human-mouse and mouse-rat comparisons, and employing three different methods to assay expression parameters. Although we confirm a negative correlation of expression breadth with both K4 and KA, we find no evidence that clustering of similarly expressed genes explains the clustering of genes of comparable substitution rates. If gene expression is not responsible, what about other causes? At least in the human-mouse comparison, the local similarity in KA can be explained by the covariation of KA and K4. As regards K4, our results appear consistent with the notion that local similarity is due to processes associated with meiotic recombination.     

Evidence for natural selection in the HAVCR1 gene: high degree of amino-acid variability in the mucin domain of human HAVCR1 protein

The family of genes encoding T-cell immunoglobulin and mucin-domain containing proteins (Tim), which are cell-surface molecules expressed in CD4(+) T helper cells, has important roles in the immune system. Here, we report three unusual patterns of genetic variation in the human hepatitis A virus cellular receptor 1 gene (HAVCR1) that are similar to patterns observed in major histocompatibility complex loci. First, levels of polymorphism in exon 4 of HAVCR1 were exceptionally high in humans (nucleotide diversity (pi)=45.45 x 10(-4)). Second, nonsynonymous substitutions and insertion/deletion variants were more frequent than synonymous substitutions in that exon (10 out of 12 variants). The rate of the mean number of nucleotide substitutions at nonsynonymous sites to synonymous sites at HAVCR1-exon 4 is >1 (P(A)/P(S)=1.92 and pi(A)/pi(S)=2.23). Third, levels of divergence among human, chimp, and gorilla sequences were unusually high in HAVCR1-exon 4 sequences. These features suggest that patterns of variation in HAVCR1 have been shaped by both positive and balancing natural selection in the course of primate evolution. Evidence that the effects of natural selection are largely restricted to the mucin domain of HAVCR1 suggests that this region may be of particular evolutionary and epidemiological interest.     

Evolutionary strata on the mouse X chromosome correspond to strata on the human X chromosome

Lahn and Page previously observed that genes on the human X chromosome were physically arranged along the chromosome in "strata," roughly ordered by degree of divergence from related genes on the Y chromosome. They hypothesized that this ordering results from a historical series of suppressions of recombination along the mammalian Y chromosome, thereby allowing formerly recombining X and Y chromosomal genes to diverge independently. Here predictions of this hypothesis are confirmed in a nonprimate mammalian order, Rodentia, through an analysis of eight gene pairs from the X and Y chromosomes of the house mouse, Mus musculus. The mouse X chromosome has been rearranged relative to the human X, so strata were not found in the same physical order on the mouse X. However, based on synonymous evolutionary distances, X-linked genes in M. musculus fall into the same strata as orthologous genes in humans, as predicted. The boundary between strata 2 and 3 is statistically significant, but the boundary between strata 1 and 2 is not significant in mice. An analysis of smaller fragments of Smcy, Smcx, Zfy, and Zfx from seven species of Mus confirmed that the strata in Mus musculus were representative of the genus Mus.     

Prevalence and evolutionary origins of autoimmune susceptibility alleles in natural mouse populations

The evolutionary origin of genetic diversity in the SLAM/CD2 gene cluster, implicated in autoimmune lupus susceptibility in mice, was investigated by sequence analysis of exons from six members of the cluster in 48 wild mouse samples derived from the global mouse population. A total of 80 coding region SNPs were identified among the six genes analyzed, indicating that this gene cluster is highly polymorphic in natural mouse populations. Phylogenetic analyses of these allelic sequences revealed clustering of alleles derived from multiple Mus species and subspecies, indicating alleles at several SLAM/CD2 loci were present in ancestral Mus populations prior to speciation and have persisted as polymorphisms for more than 1 million years. Analyses of nonsynonymous/synonymous ratios using likelihood codon substitution models identified several segments in Cd229, Cd48 and Cd84 that were impacted by positive diversifying selective pressures. These findings support the interpretation that selection favoring the generation and retention of functional polymorphisms has played a role in the evolutionary origin of genetic polymorphisms that are predisposing to autoimmunity.     

Comparative evolutionary genomics of androgen-binding protein genes

Allelic variation within the mouse androgen-binding protein (ABP) alpha subunit gene (Abpa) has been suggested to promote assortative mating and thus prezygotic isolation. This is consistent with the elevated evolutionary rates observed for the Abpa gene, and the Abpb and Abpg genes whose products (ABPbeta and ABPgamma) form heterodimers with ABPalpha. We have investigated the mouse sequence that contains the three Abpa/b/g genes, and orthologous regions in rat, human, and chimpanzee genomes. Our studies reveal extensive "remodeling" of this region: Duplication rates of Abpa-like and Abpbg-like genes in mouse are >2 orders of magnitude higher than the average rate for all mouse genes; synonymous nucleotide substitution rates are twofold higher; and the Abpabg genomic region has expanded nearly threefold since divergence of the rodents. During this time, one in six amino acid sites in ABPbetagamma-like proteins appear to have been subject to positive selection; these may constitute a site of interaction with receptors or ligands. Greater adaptive variation among Abpbg-like sequences than among Abpa-like sequences suggests that assortative mating preferences are more influenced by variation in Abpbg-like genes. We propose a role for ABPalpha/beta/gamma proteins as pheromones, or in modulating odorant detection. This would account for the extraordinary adaptive evolution of these genes, and surrounding genomic regions, in murid rodents.     

Sequencing and analysis of 10,967 full-length cDNA clones from Xenopus laevis and Xenopus tropicalis reveals post-tetraploidization transcriptome remodeling

Sequencing of full-insert clones from full-length cDNA libraries from both Xenopus laevis and Xenopus tropicalis has been ongoing as part of the Xenopus Gene Collection Initiative. Here we present 10,967 full ORF verified cDNA clones (8049 from X. laevis and 2918 from X. tropicalis) as a community resource. Because the genome of X. laevis, but not X. tropicalis, has undergone allotetraploidization, comparison of coding sequences from these two clawed (pipid) frogs provides a unique angle for exploring the molecular evolution of duplicate genes. Within our clone set, we have identified 445 gene trios, each comprised of an allotetraploidization-derived X. laevis gene pair and their shared X. tropicalis ortholog. Pairwise dN/dS, comparisons within trios show strong evidence for purifying selection acting on all three members. However, dN/dS ratios between X. laevis gene pairs are elevated relative to their X. tropicalis ortholog. This difference is highly significant and indicates an overall relaxation of selective pressures on duplicated gene pairs. We have found that the paralogs that have been lost since the tetraploidization event are enriched for several molecular functions, but have found no such enrichment in the extant paralogs. Approximately 14% of the paralogous pairs analyzed here also show differential expression indicative of subfunctionalization.     

Genome duplication,a trait shared by 22000 species of ray-finned fish

Through phylogeny reconstruction we identified 49 genes with a single copy in man, mouse, and chicken, one or two copies in the tetraploid frog Xenopus laevis, and two copies in zebrafish (Danio rerio). For 22 of these genes, both zebrafish duplicates had orthologs in the pufferfish (Takifugu rubripes). For another 20 of these genes, we found only one pufferfish ortholog but in each case it was more closely related to one of the zebrafish duplicates than to the other. Forty-three pairs of duplicated genes map to 24 of the 25 zebrafish linkage groups but they are not randomly distributed; we identified 10 duplicated regions of the zebrafish genome that each contain between two and five sets of paralogous genes. These phylogeny and synteny data suggest that the common ancestor of zebrafish and pufferfish, a fish that gave rise to approximately 22000 species, experienced a large-scale gene or complete genome duplication event and that the pufferfish has lost many duplicates that the zebrafish has retained.     

A comparative map of the zebrafish genome

Zebrafish mutations define the functions of hundreds of essential genes in the vertebrate genome. To accelerate the molecular analysis of zebrafish mutations and to facilitate comparisons among the genomes of zebrafish and other vertebrates, we used a homozygous diploid meiotic mapping panel to localize polymorphisms in 691 previously unmapped genes and expressed sequence tags (ESTs). Together with earlier efforts, this work raises the total number of markers scored in the mapping panel to 2119, including 1503 genes and ESTs and 616 previously characterized simple-sequence length polymorphisms. Sequence analysis of zebrafish genes mapped in this study and in prior work identified putative human orthologs for 804 zebrafish genes and ESTs. Map comparisons revealed 139 new conserved syntenies, in which two or more genes are on the same chromosome in zebrafish and human. Although some conserved syntenies are quite large, there were changes in gene order within conserved groups, apparently reflecting the relatively frequent occurrence of inversions and other intrachromosomal rearrangements since the divergence of teleost and tetrapod ancestors. Comparative mapping also shows that there is not a one-to-one correspondence between zebrafish and human chromosomes. Mapping of duplicate gene pairs identified segments of 20 linkage groups that may have arisen during a genome duplication that occurred early in the evolution of teleosts after the divergence of teleost and mammalian ancestors. This comparative map will accelerate the molecular analysis of zebrafish mutations and enhance the understanding of the evolution of the vertebrate genome.     

A large-scale study of the random variability of a coding sequence: a study on the CFTR gene

Coding single nucleotide substitutions (cSNSs) have been studied on hundreds of genes using small samples (n(g) approximately 100-150 genes). In the present investigation, a large random European population sample (average n(g) approximately 1500) was studied for a single gene, the CFTR (Cystic Fibrosis Transmembrane conductance Regulator). The nonsynonymous (NS) substitutions exhibited, in accordance with previous reports, a mean probability of being polymorphic (q > 0.005), much lower than that of the synonymous (S) substitutions, but they showed a similar rate of subpolymorphic (q < 0.005) variability. This indicates that, in autosomal genes that may have harmful recessive alleles (nonduplicated genes with important functions), genetic drift overwhelms selection in the subpolymorphic range of variability, making disadvantageous alleles behave as neutral. These results imply that the majority of the subpolymorphic nonsynonymous alleles of these genes are selectively negative or even pathogenic.     

Effect of gene conversion on polymorphic patterns at major histocompatibility complex loci

Summary: The pattern of polymorphisms at MHC loci was studied by computer simulations and by DNA sequence analyses, with special reference to synonymous and non-synonymous divergence among alleles. A model of overdominance plus short-term selection was simulated for a gene family with allelic and non-allelic gene conversion. It was found that the ratio of non-synonymous-to-synonymous divergence among alleles decreases as the non-allelic conversion becomes more frequent. On the other hand, the ratio of synonymous divergence at the antigen recognition site to that at non-coding region increases by non-allelic: conversion. By comparing such a pattern with the observed values of mammalian genes, it is suggested that non-allelic conversion is less frequent in human genes than in genes of other mammals. A quantitative analysis on the patchwork pattern was attempted by using identity excess among nucleotide sites of a gene. Comparison of simulation results with DNA data analysis shows that the rate of allelic conversion is high.     

Evolution of alternative splicing after gene duplication

Alternative splicing and gene duplication are two major sources of proteomic function diversity. Here, we study the evolutionary trend of alternative splicing after gene duplication by analyzing the alternative splicing differences between duplicate genes. We observed that duplicate genes have fewer alternative splice (AS) forms than single-copy genes, and that a negative correlation exists between the mean number of AS forms and the gene family size. Interestingly, we found that the loss of alternative splicing in duplicate genes may occur shortly after the gene duplication. These results support the subfunctionization model of alternative splicing in the early stage after gene duplication. Further analysis of the alternative splicing distribution in human duplicate pairs showed the asymmetric evolution of alternative splicing after gene duplications; i.e., the AS forms between duplicates may differ dramatically. We therefore conclude that alternative splicing and gene duplication may not evolve independently. In the early stage after gene duplication, young duplicates may take over a certain amount of protein function diversity that previously was carried out by the alternative splicing mechanism. In the late stage, the gain and loss of alternative splicing seem to be independent between duplicates.