Rapid divergency of rodent CD99 orthologs: implications for the evolution of the pseudoautosomal region

The human pseudoautosomal region 1 (PAR1) is essential for the obligatory X-Y crossover in male meiosis. Despite its critical role, comparative studies of human and mouse pseudoautosomal genes have been limited owing to the scarcity of genes conserved between the two species. Human CD99 is a 32-kDa cell surface protein that is encoded by the MIC2 gene localized to the PAR1. Although several sequences such as CD99L2, PBDX, and CD99L1 are related to CD99, its murine ortholog, Cd99, has not yet been identified. Here we report a novel mouse Cd99, designated D4, which shows overall sequence homology to CD99, with the highest conservation between the two genes being found in the transmembrane regions. In addition, the D4 protein displays biochemical characteristics, functional homology, and expression patterns similar to those of CD99. The D4 gene is localized on an autosome, chromosome 4, reflecting a common mapping feature with other mouse orthologs of human PAR1 genes. Furthermore, a phylogenetic analysis of CD99-related genes confirmed that the D4 gene is indeed an ortholog of CD99 and exhibits the accelerated evolution pattern of CD99 orthologs, as compared to the CD99L2 orthologs. On the basis of these findings, we suggest that CD99 belongs to the ancient PAR genes, and that the rapid interspecies divergence of its present sequence and map position is due to a high recombination frequency and the occurrence of chromosomal translocation, supporting the addition-attrition hypothesis for PAR evolution.     

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.     

ZIPK: a unique case of murine-specific divergence of a conserved vertebrate gene

Zipper interacting protein kinase (ZIPK, also known as death-associated protein kinase 3 [DAPK3]) is a Ser/Thr kinase that functions in programmed cell death. Since its identification eight years ago, contradictory findings regarding its intracellular localization and molecular mode of action have been reported, which may be attributed to unpredicted differences among the human and rodent orthologs. By aligning the sequences of all available ZIPK orthologs, from fish to human, we discovered that rat and mouse sequences are more diverged from the human ortholog relative to other, more distant, vertebrates. To test experimentally the outcome of this sequence divergence, we compared rat ZIPK to human ZIPK in the same cellular settings. We found that while ectopically expressed human ZIPK localized to the cytoplasm and induced membrane blebbing, rat ZIPK localized exclusively within nuclei, mainly to promyelocytic leukemia oncogenic bodies, and induced significantly lower levels of membrane blebbing. Among the unique murine (rat and mouse) sequence features, we found that a highly conserved phosphorylation site, previously shown to have an effect on the cellular localization of human ZIPK, is absent in murines but not in earlier diverging organisms. Recreating this phosphorylation site in rat ZIPK led to a significant reduction in its promyelocytic leukemia oncogenic body localization, yet did not confer full cytoplasmic localization. Additionally, we found that while rat ZIPK interacts with PAR-4 (also known as PAWR) very efficiently, human ZIPK fails to do so. This interaction has clear functional implications, as coexpression of PAR-4 with rat ZIPK caused nuclear to cytoplasm translocation and induced strong membrane blebbing, thus providing the murine protein a possible adaptive mechanism to compensate for its sequence divergence. We have also cloned zebrafish ZIPK and found that, like the human and unlike the murine orthologs, it localizes to the cytoplasm, and fails to bind the highly conserved PAR-4 protein. This further supports the hypothesis that murine ZIPK underwent specific divergence from a conserved consensus. In conclusion, we present a case of species-specific divergence occurring in a specific branch of the evolutionary tree, accompanied by the acquisition of a unique protein–protein interaction that enables conservation of cellular function.     

Reverse polarization in amino acid and nucleotide substitution patterns between human-mouse orthologs of two compositional extrema.

Genome-wide analysis of sequence divergence patterns in 12,024 human-mouse orthologous pairs reveals, for the first time, that the trends in nucleotide and amino acid substitutions in orthologs of high and low GC composition are highly asymmetric and polarized to opposite directions. The entire dataset has been divided into three groups on the basis of the GC content at third codon sites of human genes: high, medium, and low. High-GC orthologs exhibit significant bias in favor of the replacements, Thr --> Ala, Ser --> Ala, Val --> Ala, Lys --> Arg, Asn --> Ser, Ile --> Val etc., from mouse to human, whereas in low-GC orthologs, the reverse trends prevail. In general, in the high-GC group, residues encoded by A/U-rich codons of mouse proteins tend to be replaced by the residues encoded by relatively G/C-rich codons in their human orthologs, whereas the opposite trend is observed among the low-GC orthologous pairs. The medium-GC group shares some trends with high-GC group and some with low-GC group. The only significant trend common in all groups of orthologs, irrespective of their GC bias, is (Asp)(Mouse) --> (Glu)(Human) replacement. At the nucleotide level, high-GC orthologs have undergone a large excess of (A/T)(Mouse) --> (G/C)(Human) substitutions over (G/C)(Mouse) --> (A/T)(Human) at each codon position, whereas for low-GC orthologs, the reverse is true.     

Human SNPs reveal no evidence of frequent positive selection

We compared the single-nucleotide polymorphisms (SNPs) in humans in 182 housekeeping and 148 tissue-specific genes. SNPs were divided into rare and common polymorphisms based on their frequencies. We found that housekeeping genes tend to be less polymorphic than tissue-specific genes for both rare and common SNPs. Using mouse as a second species for computing sequence divergences, we found no evidence of positive selection: for both housekeeping and tissue-specific genes, the ratio of nonsynonymous to synonymous common SNPs per site showed no significant difference from that of divergence. Similarly, we observed no evidence of positive selection for the 289 and 149 genes that have orthologs available for divergence calculation between humans and chimpanzees and between humans and Old World monkeys, respectively. A comparison with previous SNP studies suggests that approximately 20% of the nonsynonymous SNPs in the human population are nearly neutral and that positive selection in the human genome might not be as frequent as previously thought.     

Testing the ortholog conjecture with comparative functional genomic data from mammals

A common assumption in comparative genomics is that orthologous genes share greater functional similarity than do paralogous genes (the "ortholog conjecture"). Many methods used to computationally predict protein function are based on this assumption, even though it is largely untested. Here we present the first large-scale test of the ortholog conjecture using comparative functional genomic data from human and mouse. We use the experimentally derived functions of more than 8,900 genes, as well as an independent microarray dataset, to directly assess our ability to predict function using both orthologs and paralogs. Both datasets show that paralogs are often a much better predictor of function than are orthologs, even at lower sequence identities. Among paralogs, those found within the same species are consistently more functionally similar than those found in a different species. We also find that paralogous pairs residing on the same chromosome are more functionally similar than those on different chromosomes, perhaps due to higher levels of interlocus gene conversion between these pairs. In addition to offering implications for the computational prediction of protein function, our results shed light on the relationship between sequence divergence and functional divergence. We conclude that the most important factor in the evolution of function is not amino acid sequence, but rather the cellular context in which proteins act.     

In search of the vertebrate phylotypic stage: a molecular examination of the developmental hourglass model and von Baer's third law

In 1828, Karl von Baer proposed a set of four evolutionary "laws" pertaining to embryological development. According to von Baer's third law, young embryos from different species are relatively undifferentiated and resemble one another but as development proceeds, distinguishing features of the species begin to appear and embryos of different species progressively diverge from one another. An expansion of this law, called "the hourglass model," has been proposed independently by Denis Duboule and Rudolf Raff in the 1990s. According to the hourglass model, ontogeny is characterized by a starting point at which different taxa differ markedly from one another, followed by a stage of reduced intertaxonomic variability (the phylotypic stage), and ending in a von-Baer-like progressive divergence among the taxa. A possible "translation" of the hourglass model into molecular terminology would suggest that orthologs expressed in stages described by the tapered part of the hourglass should resemble one another more than orthologs expressed in the expansive parts that precede or succeed the phylotypic stage. We tested this hypothesis using 1,585 mouse genes expressed during 26 embryonic stages, and their human orthologs. Evolutionary divergence was estimated at different embryonic stages by calculating pairwise distances between corresponding orthologous proteins from mouse and human. Two independent datasets were used. One dataset contained genes that are expressed solely in a single developmental stage; the second was made of genes expressed at different developmental stages. In the second dataset the genes were classified according to their earliest stage of expression. We fitted second order polynomials to the two datasets. The two polynomials displayed minima as expected from the hourglass model. The molecular results suggest, albeit weakly, that a phylotypic stage (or period) indeed exists. Its temporal location, sometimes between the first-somites stage and the formation of the posterior neuropore, was in approximate agreement with the morphologically defined phylotypic stage. The molecular evidence for the later parts of the hourglass model, i.e., for von Baer's third law, was stronger than that for the earlier parts.     

Extracellular Ionic Locks Determine Variation in Constitutive Activity and Ligand Potency between Species Orthologs of the Free Fatty Acid Receptors FFA2 and FFA3

Free fatty acid receptors 2 and 3 (FFA2 and FFA3) are G protein-coupled receptors for short chain free fatty acids (SCFAs). They respond to the same set of endogenous ligands but with distinct rank-order of potency such that acetate (C2) has been described as FFA2-selective, whereas propionate (C3) is non-selective. Although C2 was confirmed to be selective for human FFA2 over FFA3, this ligand was not selective between the mouse orthologs. Moreover, although C3 was indeed not selective between the human orthologs, it displayed clear selectivity for mouse FFA3 over mouse FFA2. This altered selectivity to C2 and C3 resulted from broad differences in SCFAs potency at the mouse orthologs. In studies to define the molecular basis for these observations, marked variation in ligand-independent constitutive activity was identified using a [³⁵S]GTPγS assay. The orthologs with higher potency for the SCFAs, human FFA2 and mouse FFA3, displayed high constitutive activity in this assay, whereas the orthologs with lower potency for the agonist ligands, mouse FFA2 and human FFA3, did not. Sequence alignments of the second extracellular loop identified single negatively charged residues in FFA2 and FFA3 not conserved between species and predicted to form ionic lock interactions with arginine residues within the FFA2 or FFA3 agonist binding pocket to regulate constitutive activity and SCFA potency. Reciprocal mutation of these residues between species orthologs resulted in the induction (or repression) of constitutive activity and in most cases also yielded corresponding changes in SCFA potency.     

Comparative mapping of human Chromosome 14q11.2-q13 genes with mouse homologous gene regions

An examination of the synteny blocks between mouse and human chromosomes aids in understanding the evolution of chromosome divergence between these two species. We comparatively mapped the human (HSA) Chromosome (Chr) 14q11.2-q13 cytogenetic region with the intervals of orthologous genes on mouse (MMU) chromosomes. A lack of conserved gene order was identified between the human cytogenetic region and the interval of orthologs on MMU 12. The evolutionary breakpoint junction was defined within 2.5 Mb, where the conserved synteny of genes on HSA 14 changes from MMU 12 to MMU 14. At the evolutionary breakpoint junction, a human EST (GI: 1114654) with identity to the human and mouse BCL2 interacting gene, BNIP3, was mapped to mouse Chr 3. New gene homologs of LAMB1, MEOX2, NRCAM, and NZTF1 were identified on HSA 7 and on the proximal cytogenetic region of HSA 14 by mapping mouse genes recently reported to be genetically linked within the relevant MMU 12 interval. This study contributes to the identification of homology relationships between the genes of HSA 14q11.2-q13 and mouse Chr 3, 12, and 14. Received: 16 March 2000 / Accepted: 16 June 2000     

Cloning, Pharmacology, and Tissue Distribution of G-Protein-Coupled Receptor GPR105 (KIAA0001) Rodent Orthologs

It has recently been shown that UDP-glucose is a potent agonist of the orphan G-protein-coupled receptor (GPCR) KIAA0001. Here we report cloning and analysis of the rat and mouse orthologs of this receptor. In accordance with GPCR nomenclature, we have renamed the cDNA clone, KIAA0001, and its orthologs GPR105 to reflect their functionality as G-protein-coupled receptors. The rat and mouse orthologs show 80% and 83% amino acid identity, respectively, to the human GPR105 protein. We demonstrate by genomic Southern blot analysis that there are no genes in the mouse or rat genomes with higher sequence similarity. Chromosomal mapping shows that the mouse and human genes are located on syntenic regions of chromosome 3. Further analyses of the rat and mouse GPR105 proteins show that they are activated by the same agonists as the human receptor, responding to UDP-glucose and closely related molecules with similar affinities. The mouse and rat receptors are widely expressed, as is the human receptor. Thus we conclude that we have identified the rat and mouse orthologs of the human gene GPR105.     

Defining the role of essential genes in human disease

A greater understanding of the causes of human disease can come from identifying characteristics that are specific to disease genes. However, a full understanding of the contribution of essential genes to human disease is lacking, due to the premise that these genes tend to cause developmental abnormalities rather than adult disease. We tested the hypothesis that human orthologs of mouse essential genes are associated with a variety of human diseases, rather than only those related to miscarriage and birth defects. We segregated human disease genes according to whether the knockout phenotype of their mouse ortholog was lethal or viable, defining those with orthologs producing lethal knockouts as essential disease genes. We show that the human orthologs of mouse essential genes are associated with a wide spectrum of diseases affecting diverse physiological systems. Notably, human disease genes with essential mouse orthologs are over-represented among disease genes associated with cancer, suggesting links between adult cellular abnormalities and developmental functions. The proteins encoded by essential genes are highly connected in protein-protein interaction networks, which we find correlates with an over-representation of nuclear proteins amongst essential disease genes. Disease genes associated with essential orthologs also are more likely than those with non-essential orthologs to contribute to disease through an autosomal dominant inheritance pattern, suggesting that these diseases may actually result from semi-dominant mutant alleles. Overall, we have described attributes found in disease genes according to the essentiality status of their mouse orthologs. These findings demonstrate that disease genes do occupy highly connected positions in protein-protein interaction networks, and that due to the complexity of disease-associated alleles, essential genes cannot be ignored as candidates for causing diverse human diseases.     

Resolving the ortholog conjecture: orthologs tend to be weakly, but significantly, more similar in function than paralogs

The function of most proteins is not determined experimentally, but is extrapolated from homologs. According to the "ortholog conjecture", or standard model of phylogenomics, protein function changes rapidly after duplication, leading to paralogs with different functions, while orthologs retain the ancestral function. We report here that a comparison of experimentally supported functional annotations among homologs from 13 genomes mostly supports this model. We show that to analyze GO annotation effectively, several confounding factors need to be controlled: authorship bias, variation of GO term frequency among species, variation of background similarity among species pairs, and propagated annotation bias. After controlling for these biases, we observe that orthologs have generally more similar functional annotations than paralogs. This is especially strong for sub-cellular localization. We observe only a weak decrease in functional similarity with increasing sequence divergence. These findings hold over a large diversity of species; notably orthologs from model organisms such as E. coli, yeast or mouse have conserved function with human proteins.     

Discovery of gene families and alternatively spliced variants by RecA-mediated cloning

Probing the functional complexity of the human genome will require new gene cloning techniques, not only to discover intraspecies gene homologs and interspecies gene orthologs, but also to identify alternatively spliced gene variants. We report homologous cDNA cloning methods that allow cloning of gene family members, genes from different species, and alternatively spliced gene variants. We cloned human 14-3-3 gene family members using DNA probes with as much as 35% sequence divergence, cloned alternatively spliced gene forms of Rad51D, and cloned a novel splice form of the human 14-3-3 theta gene with a unique expression pattern. Interspecies gene cloning was demonstrated for the mouse Rad51C and mouse beta-actin genes using human gene probes. The gene family cloning method is fast, efficient, and free from PCR errors; moreover, it exploits the abilities of RecA protein to pair homologous or partially homologous DNA sequences stably in kinetically trapped, multistranded DNA hybrids that can be used for subsequent gene clone enrichment.     

Global analysis of small molecule binding to related protein targets

We report on the integration of pharmacological data and homology information for a large scale analysis of small molecule binding to related targets. Differences in small molecule binding have been assessed for curated pairs of human to rat orthologs and also for recently diverged human paralogs. Our analysis shows that in general, small molecule binding is conserved for pairs of human to rat orthologs. Using statistical tests, we identified a small number of cases where small molecule binding is different between human and rat, some of which had previously been reported in the literature. Knowledge of species specific pharmacology can be advantageous for drug discovery, where rats are frequently used as a model system. For human paralogs, we demonstrate a global correlation between sequence identity and the binding of small molecules with equivalent affinity. Our findings provide an initial general model relating small molecule binding and sequence divergence, containing the foundations for a general model to anticipate and predict within-target-family selectivity.     

Density peaks of paralog pairs in human and mouse genomes

Paralog gene trees, which reflect the increase of genomic complexity in the evolution, can be complicated and ambiguous. A simpler complementary approach is analysis of density distribution of paralog pairs. It can reveal general features of genome evolution, which may be hidden in the forest of gene trees. It is known that distribution of human paralog pairs along the axis of protein divergence between pair members forms two main peaks. Here I show that there are three main peaks in the mouse genome. Thus, the multimodality of paralog pair distribution seems to be a fundamental feature of mammalian genomes. Despite the great diversity of domains presented in small amounts or in multidomain architectures with a few predominant domains, both in human and mouse the first peak consists mostly of gene pairs with zinc finger domains or olfactory receptor domain. In the mouse the olfactory receptor predominates, which stipulates the three-peak distribution (since in the olfactory receptors the second peak is closer to the first peak than in other genes). The mammalian-wide zinc finger orthologs are biased towards the second peak. Thus, the marsupial orthologs are nearly absent in the first peak of human and mouse. The gene pairs in the first peak show a lower ratio of nonsynonymous to synonymous substitutions, which suggests that their evolution is more constrained. The plausible explanation is that they are in subfunctionalization state (partition of initial function of ancestral gene), whereas the second peak contains gene pairs that are already in neofunctionalization state (acquiring of novel functions). These data suggest that the adaptive radiation of mammals was accompanied by a burst of duplication of zinc finger genes, which are located in the first (most recent) peak of paralog pairs.