The Ghost of Selection Past: Rates of Evolution and Functional Divergence of Anciently Duplicated Genes

The duplication of genes and even complete genomes may be a prerequisite for major evolutionary transitions and the origin of evolutionary novelties. However, the evolutionary mechanisms of gene evolution and the origin of novel gene functions after gene duplication have been a subject of many debates. Recently, we compiled 26 groups of orthologous genes, which included one gene from human, mouse, and chicken, one or two genes from the tetraploid Xenopus and two genes from zebrafish. Comparative analysis and mapping data showed that these pairs of zebrafish genes were probably produced during a fish-specific genome duplication that occurred between 300 and 450 Mya, before the teleost radiation (Taylor et al. 2001). As discussed here, many of these retained duplicated genes code for DNA binding proteins. Different models have been developed to explain the retention of duplicated genes and in particular the subfunctionalization model of Force et al. (1999) could explain why so many developmental control genes have been retained. Other models are harder to reconcile with this particular set of duplicated genes. Most genes seem to have been subjected to strong purifying selection, keeping properties such as charge and polarity the same in both duplicates, although some evidence was found for positive Darwinian selection, in particular for Hox genes. However, since only the cumulative pattern of nucleotide substitutions can be studied, clear indications of positive Darwinian selection or neutrality may be hard to find for such anciently duplicated genes. Nevertheless, an increase in evolutionary rate in about half of the duplicated genes seems to suggest that either positive Darwinian selection has occurred or that functional constraints have been relaxed at one point in time during functional divergence. Received: 4 January 2001 / Accepted: 29 March 2001     

The pseudokinase domain is required for suppression of basal activity of Jak2 and Jak3 tyrosine kinases and for cytokine-inducible activation of signal transduction

Janus (Jak) tyrosine kinases contain a tyrosine kinase (JH1) domain adjacent to a catalytically inactive pseudokinase domain (JH2). The JH2 domain has been implicated in regulation of Jak activity, but its function remains poorly understood. Here, we found that the JH2 domain negatively regulates the activity of Jak2 and Jak3. Deletion of JH2 resulted in increased tyrosine phosphorylation of the Jak2- and Jak3-JH2 deletion mutants as well as of coexpressed STAT5. In cytokine receptor signaling, the deletion of the Jak2- and Jak3-JH2 domains resulted in interferon-gamma and interleukin-2-independent STAT activation, respectively. However, cytokine stimulations did not further induce the JH2 deletion mutant-mediated STAT activation. The deletion of the Jak2 JH2 domain also abolished interferon-gamma-inducible kinase activation, although it did not affect the reciprocal Jak1-Jak2 interaction in 293T cells. Chimeric constructs, where the JH2 domains were swapped between Jak2 and Jak3, retained low basal activity and cytokine inducible signaling, indicating functional conservation between the two JH2 domains. However, the basal activity of Jak2 was significantly lower than that of Jak3, suggesting differences in the regulation of Jak2 and Jak3 activity. In conclusion, we found that the JH2 domain has a conserved function in Jak2 and Jak3. The JH2 domain is required for two distinct functions in cytokine signaling: (i) inhibition of the basal activity of Jak2 and Jak3, and (ii) cytokine-inducible activation of signaling. The Jak-JH2 deletion mutants are catalytically active, activate STAT5, and interact with another Jak kinase, but the JH2 domain is required to connect these signaling events to receptor activation. Thus, we propose that the JH2 domain contributes to both the uninduced and ligand-induced Jak-receptor complex, where it acts as a cytokine-inducible switch to regulate signal transduction.     

Evolutionary history and complementary selective relaxation of the duplicated PI genes in grasses

Gene duplication plays an important role in the evolution of organisms by allowing functional innovation and the divergence of duplicate genes. Previous studies found two PI-like genes in grass species, suggesting functional divergence between the paralogous copies. Here, we reconstructed the evolutionary history of two PI genes from major lineages of grasses and other monocot species, and demonstrated that two PI genes (PI1 and PI2) arose from a whole genome duplication that occurred in a common ancestor of extant grasses. Molecular evolutionary analyses at the family and tribal levels found strong purifying selection acting on two genes in grasses, consistent with the conserved class B function of the PI genes. Importantly, we detected different patterns of selective relaxation between the duplicated PI genes although no signature of positive selection was found. Likelihood ratio tests revealed that the ω ratio for M domain is significantly higher in PI1 than in PI2 but that for K domain is significantly higher in PI2 than in PI1. These findings imply that complementary selective relaxation occurs in two PI genes after duplication, and provide additional molecular evidence for the subfunctionalization of the duplicated PI genes in grasses.     

Phylogenies of Developmentally Important Proteins Do Not Support the Hypothesis of Two Rounds of Genome Duplication Early in Vertebrate History

It has been proposed that two rounds of duplication of the entire genome (polyploidization) occurred early in vertebrate history (the 2R hypothesis); and the observation that certain gene families important in regulating development have four members in vertebrates, as opposed to one in Drosophila, has been adduced as evidence in support of this hypothesis. However, such a pattern of relationship can be taken as support of the 2R hypothesis only if (1) the four vertebrate genes can be shown to have diverged after the origin of vertebrates, and (2) the phylogeny of the four vertebrate genes (A–D) exhibits a topology of the form (AB) (CD), rather than (A) (BCD). In order to test the 2R hypothesis, I constructed phylogenies for nine protein families important in development. Only one showed a topology of the form (AB) (CD), and that received weak statistical support. In contrast, four phylogenies showed topologies of the form (A) (BCD) with statistically significant support. Furthermore, in two cases there was significant support for duplication of the vertebrate genes prior to the divergence of deuterostomes and protostomes: in one case there was significant support for duplication of the vertebrate genes at least prior to the divergence of vertebrates and urochordates, and in one case there was weak support for duplication of the vertebrate genes prior to the divergence of deuterostomes and protostomes. Taken together with other recently published phylogenies of developmentally important genes, these results provide strong evidence against the 2R hypothesis. Received: 22 December 1997 / Accepted: 5 October 1998     

Nucleotide Substitution Rate of Mammalian Mitochondrial Genomes

We present here for the first time a comprehensive study based on the analysis of closely related organisms to provide an accurate determination of the nucleotide substitution rate in mammalian mitochondrial genomes. This study examines the evolutionary pattern of the different functional mtDNA regions as accurately as possible on the grounds of available data, revealing some important ``genomic laws.' The main conclusions can be summarized as follows. (1) High intragenomic variability in the evolutionary dynamic of mtDNA was found. The substitution rate is strongly dependent on the region considered, and slow- and fast-evolving regions can be identified. Nonsynonymous sites, the D-loop central domain, and tRNA and rRNA genes evolve much more slowly than synonymous sites and the two peripheral D-loop region domains. The synonymous rate is fairly uniform over the genome, whereas the rate of nonsynonymous sites depends on functional constraints and therefore differs considerably between genes. (2) The commonly accepted statement that mtDNA evolves more rapidly than nuclear DNA is valid only for some regions, thus it should be referred to specific mitochondrial components. In particular, nonsynonymous sites show comparable rates in mitochondrial and nuclear genes; synonymous sites and small rRNA evolve about 20 times more rapidly and tRNAs about 100 times more rapidly in mitochondria than in their nuclear counterpart. (3) A species-specific evolution is particularly evident in the D-loop region. As the divergence times of the organism pairs under consideration are known with sufficient accuracy, absolute nucleotide substitution rates are also provided. Received: 11 May 1998 / Accepted: 2 September 1998     

Extensive Gene Duplication in the Early Evolution of Animals Before the Parazoan–Eumetazoan Split Demonstrated by G Proteins and Protein Tyrosine Kinases from Sponge and Hydra

To know whether genes involved in cell–cell communication typical of multicellular animals dramatically increased in concert with the Cambrian explosion, the rapid evolutionary burst in the major groups of animals, and whether these genes exist in the sponge lacking cell cohesiveness and coordination typical of eumetazoans, we have carried out cloning of the G-protein α subunit (Gα) and the protein tyrosine kinase (PTK) cDNAs from Ephydatia fluviatilis (freshwater sponge) and Hydra magnipapillata strain 105 (hydra). We obtained 13 Gα and 20 PTK cDNAs. Generally animal gene families diverged first by gene duplication (subtype duplication) that gave rise to diverse subtypes with different primary functions, followed by further gene duplication in the same subtype (isoform duplication) that gave rise to isoform genes with virtually identical function. Phylogenetic trees of Gα and PTK families including cDNAs from sponge and hydra revealed that most of the present-day subtypes had been established in the very early evolution of animals before the parazoan–eumetazoan split, the earliest branching among the extant animal phyla, by extensive subtype duplication: for PTK and Gα families, 23 and 9 subtype duplications were observed in the early stage before the parazoan–eumetazoan split, respectively, and after that split, only 2 and 1 subtype duplications were found, respectively. After the separation from arthropods, vertebrates underwent frequent isoform duplications before the fish–tetrapod split. Furthermore, rapid amino acid changes appear to have occurred in concert with the extensive subtype duplication and isoform duplication. Thus the pattern of gene diversification during animal evolution might be characterized by bursts of gene duplication interrupted by considerably long periods of silence, instead of proceeding gradually, and there might be no direct link between the Cambrian explosion and the extensive gene duplication that generated diverse functions (subtypes) of these families. Received: 4 November 1998 / Accepted: 17 November 1998     

The Evolutionary History of Prosaposin: Two Successive Tandem-Duplication Events Gave Rise to the Four Saposin Domains in Vertebrates

Prosaposin is a multifunctional protein encoded by a single-copy gene. It contains four saposin domains (A, B, C, and D) occurring as tandem repeats connected by linker sequences. Because the saposin domains are similar to one another, it is deduced that they were created by sequential duplications of an ancestral domain. There are two types of evolutionary scenarios that may explain the creation of the four-domain gene: (1) two rounds of tandem internal gene duplication and (2) three rounds of duplications. An evolutionary and phylogenetic analysis of saposin DNA and amino acid sequences from human, mouse, rat, chicken, and zebrafish indicates that the first evolutionary scenario is the most likely. Accordingly, an ancestral saposin-unit duplication produced a two-domain gene, which, subsequently, underwent a second complete tandem duplication to give rise to the present four-domain structure of the prosaposin gene. Received: 8 February 2001 / Accepted: 29 June 2001     

Phylogenetic Analyses Reveal Ancient Duplication of Estrogen Receptor Isoforms

To determine the origin and evolutionary significance of a recently discovered isoform of the estrogen receptor (ERβ), we examined the phylogenetic relationship of ERβ to the well-known α isoform (ERα) and other steroid receptors. Our phylogenetic analyses traced the origin of ERβ to a single duplication event at least 450 million years ago. Since this duplication, the evolution of both ER isoforms has apparently been constrained such that 80% of the amino acid positions in the DNA binding domain (DBD) and 53% of the ligand binding domain (LBD) have remained unchanged. Using the phylogenetic tree, we determined the amount of evolutionary change that had occurred in two ER isoforms. The DBD and the LBD had lower rates of evolutionary change compared to the NH₂ terminal domain. However, even with strong selective constraints on the DBD and LBD, our phylogenetic analyses demonstrate two clearly separate phylogenetic histories for ERα and ERβ dating back several hundred million years. The ancient duplication of ER and the parallel evolution of the two ER isoforms suggest that, although ERα and ERβ share a substantial degree of sequence identity, they play unique roles in vertebrate physiology and reproduction. Received: 19 January 1999 / Accepted: 26 May 1999     

Protein Tyrosine Kinase cDNAs from Amphioxus, Hagfish, and Lamprey: Isoform Duplications Around the Divergence of Cyclostomes and Gnathostomes

Animals evolved a variety of gene families involved in cell–cell communication and developmental control by gene duplication and domain shuffling. Each family is made up of several subtypes or subfamilies with distinct structures and functions, which diverged by gene duplications and domain shufflings before the divergence of parazoans and eumetazoans. Since the separation from protostomes, vertebrates expanded the multiplicity of members (isoforms) in the same subfamily by further gene duplications in their early evolution before the fish–tetrapod split. To know the dates of isoform duplications more closely, we have conducted isolation and sequencing cDNAs encoding the fibroblast growth factor receptor, Eph, src, and platelet-derived growth factor receptor subtypes belonging to the protein tyrosine kinase family from Branchiostoma belcheri, an amphioxus, Eptatretus burgeri, a hagfish, and Lampetra reissneri, a lamprey. From a phylogenetic tree of each subfamily inferred from a maximum likelihood (ML) method, together with a bootstrap analysis based on the ML method, we have shown that the isoform duplications frequently occurred in the early evolution of vertebrates around or just before the divergence of cyclostomes and gnathostomes by gene duplications and possibly chromosomal duplications. Received: 28 April 1998 / Accepted: 30 June 1999     

Determinants for the interaction between Janus kinase 2 and protein phosphatase 2A

We have shown that the serine/threonine protein phosphatase 2A (PP2A) associates with the Jak2 tyrosine kinase in a myeloid progenitor line. In this study, we characterized the regions of Jak2 and PP2A responsible for association and evaluated the functional consequences of association. We demonstrate that PP2A interacts with truncated forms of Jak2 containing the JH1 catalytic domain. Using GST fusion proteins, we show that the isolated JH1 and JH3 domains of Jak2 bind directly to PP2A. Jak2 contains putative PP2A binding sequences (LXXLL) in the JH1 domain (residues 1078-1082) and in the JH3 domain (residues 474-478). Mutation of the LXXLL sequence in the JH1 domain decreased PP2A binding in vitro, while mutation of the similar JH3 sequence did not affect PP2A binding. We analyzed full-length Jak2 bearing the LXXLL mutation in Cos-7 cells for association with PP2A. The JH1 mutation impaired Jak2 activity and had a modest effect on PP2A binding. Finally, we show that a mutant form of the PP2A catalytic subunit lacking a site for phosphorylation (Y307F) binds more tightly to Jak2 than wild-type PP2A, consistent with a model where phosphorylation disrupts the Jak2-PP2A interaction.     

Autoinhibition of Jak2 tyrosine kinase is dependent on specific regions in its pseudokinase domain

Jak tyrosine kinases have a unique domain structure containing a kinase domain (JH1) adjacent to a catalytically inactive pseudokinase domain (JH2). JH2 is crucial for inhibition of basal Jak activity, but the mechanism of this regulation has remained elusive. We show that JH2 negatively regulated Jak2 in bacterial cells, indicating that regulation is an intrinsic property of Jak2. JH2 suppressed basal Jak2 activity by lowering the V(max) of Jak2, whereas JH2 did not affect the K(m) of Jak2 for a peptide substrate. Three inhibitory regions (IR1-3) within JH2 were identified. IR3 (residues 758-807), at the C terminus of JH2, directly inhibited JH1, suggesting an inhibitory interaction between IR3 and JH1. Molecular modeling of JH2 showed that IR3 could form a stable alpha-helical fold, supporting that IR3 could independently inhibit JH1. IR2 (725-757) in the C-terminal lobe of JH2, and IR1 (619-670), extending from the N-terminal to the C-terminal lobe, enhanced IR3-mediated inhibition of JH1. Disruption of IR3 either by mutations or a small deletion increased basal Jak2 activity, but abolished interferon-gamma-inducible signaling. Together, the results provide evidence for autoinhibition of a Jak family kinase and identify JH2 regions important for autoregulation of Jak2.     

Tyrosine kinase 2 interacts with and phosphorylates receptor for activated C kinase-1, a WD motif-containing protein

Receptor for activated C kinase (Rack)-1 is a protein kinase C-interacting protein, and contains a WD repeat but has no enzymatic activity. In addition to protein kinase C, Rack-1 also binds to Src, phospholipase Cgamma, and ras-GTPase-activating proteins. Thus, Rack-1 is thought to function as a scaffold protein that recruits specific signaling elements. In a cytokine signaling cascade, Rack-1 has been reported to interact with the IFN-alphabeta receptor and Stat1. In addition, we show here that Rack-1 associates with a member of Jak, tyrosine kinase 2 (Tyk2). Rack-1 interacts weakly with the kinase domain and interacts strongly with the pseudokinase domain of Tyk2. Rack-1 associates with Tyk2 via two regions, one in the N terminus and one in the middle portion (aa 138-203) of Rack-1. Jak activation causes the phosphorylation of tyrosine 194 on Rack-1. After phosphorylation, Rack-1 is translocated toward the perinuclear region. In addition to functioning as a scaffolding protein, these results raise the possibility that Rack-1 functions as a signaling molecule in cytokine signaling cascades.     

Origin and Evolution of Paralogous rRNA Gene Clusters Within the Flatworm Family Dugesiidae (Platyhelminthes, Tricladida)

Analysis of the 18S rDNA sequences of five species of the family Dugesiidae (phylum Platyhelminthes, suborder Tricladida, infraorder Paludicola) and eight species belonging to families Dendrocoelidae and Planaridae and to the infraorder Maricola showed that members of the family Dugesiidae have two types of 18S rDNA genes, while the rest of the species have only one. The duplication event also affected the ITS-1, 5.8S, ITS-2 region and probably the 28S gene. The mean divergence value between the type I and the type II sequences is 9% and type II 18S rDNA genes are evolving 2.3 times more rapidly than type I. The evolutionary rates of type I and type II genes were calibrated from biogeographical data, and an approximate date for the duplication event of 80–120 million years ago was calculated. The type II gene was shown, by RT-PCR, to be transcribed in adult individuals of Schmidtea polychroa, though at very low levels. This result, together with the fact that most of the functionally important positions for small-subunit rRNA in prokaryotes have been conserved, indicates that the type II gene is probably functional. Received: 24 March 1998 / Accepted: 17 March 1999     

Molecular Evolution of Nitrogen Fixation: The Evolutionary History of the nifD, nifK, nifE, and nifN Genes

The pairs of nitrogen fixation genes nifDK and nifEN encode for the α and β subunits of nitrogenase and for the two subunits of the NifNE protein complex, involved in the biosynthesis of the FeMo cofactor, respectively. Comparative analysis of the amino acid sequences of the four NifD, NifK, NifE, and NifN in several archaeal and bacterial diazotrophs showed extensive sequence similarity between them, suggesting that their encoding genes constitute a novel paralogous gene family. We propose a two-step model to reconstruct the possible evolutionary history of the four genes. Accordingly, an ancestor gene gave rise, by an in-tandem paralogous duplication event followed by divergence, to an ancestral bicistronic operon; the latter, in turn, underwent a paralogous operon duplication event followed by evolutionary divergence leading to the ancestors of the present-day nifDK and nifEN operons. Both these paralogous duplication events very likely predated the appearance of the last universal common ancestor. The possible role of the ancestral gene and operon in nitrogen fixation is also discussed. Received: 21 June 1999 / Accepted: 1 March 2000     

Evolutionary Patterns of Gene Families Generated in the Early Stage of Vertebrates

In this paper we have analyzed 49 vertebrate gene families that were generated in the early stage of vertebrates and/or shortly before the origin of vertebrates, each of which consists of three or four member genes. We have dated the first (T₁) and second (T₂) gene duplications of 26 gene families with 3 member genes. The means of T₁ (594 mya) and T₂ (488 mya) are largely consistent to a well-cited version of two-round (2R) genome duplication theory. Moreover, in most cases, the time interval between two successive gene duplications is large enough that the fate of duplicate genes generated by the first gene duplication was likely to be determined before the second one took place. However, the phylogenetic pattern of 23 gene families with 4 members is complicated; only 5 of them are predicted by 2R model, but 11 families require an additional gene (or genome) duplication. For the rest (7 families), at least one gene duplication event had occurred before the divergence between vertebrate and Drosophila, indicating a possible misleading of the 4:1 rule (member gene ratio between vertebrates and invertebrates). Our results show that Ohno's 2R conjecture is valid as a working hypothesis for providing a most parsimonious explanation. Although for some gene families, additional gene duplication is needed, the credibility of the third genome duplication (3R) remains to be investigated. Received: 13 December 1999 / Accepted: 7 April 2000     

Natural history and functional divergence of protein tyrosine kinases

Cellular signaling is important for many biological processes including growth, differentiation, adhesion, motility and apoptosis. The protein tyrosine kinase (PTK) supergene family is the key mediator in cellular signaling in metazoans, directly associated with a variety of human diseases. All PTKs contain a highly conserved catalytic kinase domain, in spite of variable multi-domain structures. Within each PTK gene family, members exhibit functional divergence in substrate-specificity or temporal/tissue-specific expression, although their primary function is conserved. After conducting phylogenetic analysis on major PTK gene families, we found that the expanding of each PTK family was likely caused by gene or genome duplication event(s) that occurred before the emergence of teleosts but after the vertebrate-amphioxus split. We further investigated the evolutionary pattern of functional divergence after gene duplication in those gene families. Our results show that site-specific shifted evolutionary rate (altered functional constraint) is a common pattern in PTK gene family evolution.     

Genomic Organization Around the Centromeric End of the HLA Class I Region: Large-Scale Sequence Analysis

We previously sequenced two regions around the centromeric end of HLA class I and the boundary between class I and class III. In this paper we analyze the two regions of about 385 kb and confirm, giving a new line of evidence, that the following two pairs of the genomic segments were duplicated in evolution: (i) a 43-kb genomic segment including the HLA-B gene showing the highest polymorphism among the classical HLA class I loci (class Ia) and a 40-kb segment including the HLA-C locus showing the lowest polymorphism and (ii) a 52-kb segment including the MIC (MHC class I chain related gene) B and a 35-kb segment including MICA. We also found that repetitive elements such as SINEs, LINEs, and LTRs occupy as much as 47% of nucleotides in this 385-kb region. This unusually high content of repetitive elements indicates that repeat-mediated rearrangements have frequently occurred in the evolutionary history of the HLA class Ia region. Analysis of LINE compositions within the two pairs of duplicated segments revealed that (i) LINEs in these regions had been dispersed prior to both the duplication of the HLA-B and -C loci and the duplication of the MICB and MICA loci, and (ii) the divergence of the HLA-B and -C loci occurred prior to the duplication of the MICA and MICB loci. To find novel genes responsible for HLA class I-associated or other diseases, we performed computer analysis applying GenScan and GRAIL to GenBank's dbEST. As a result, at least five as yet uncharacterized genes were newly mapped on the HLA class I centromeric region studied. These novel genes should be analyzed further to determine their relationships to diseases associated with this region. Received: 16 June 1998 / Accepted: 18 August 1998     

Rapid evolution through gene duplication and subfunctionalization of the testes-specific alpha4 proteasome subunits in Drosophila

Gene duplication is an important mechanism for acquiring new genes and creating genetic novelty in organisms. Evidence suggests that duplicated genes are retained at a much higher rate than originally thought and that functional divergence of gene copies is a major factor promoting their retention in the genome. We find that two Drosophila testes-specific alpha4 proteasome subunit genes (alpha4-t1 and alpha4-t2) have a higher polymorphism within species and are significantly more diverged between species than the somatic alpha4 gene. Our data suggest that following gene duplication, the alpha4-t1 gene experienced relaxed selective constraints, whereas the alpha4-t2 gene experienced positive selection acting on several codons. We report significant heterogeneity in evolutionary rates among all three paralogs at homologous codons, indicating that functional divergence has coincided with genic divergence. Reproductive subfunctionalization may allow for a more rapid evolution of reproductive traits and a greater specialization of testes function. Our data add to the increasing evidence that duplicated genes experience lower selective constraints and in some cases positive selection following duplication. Newly duplicated genes that are freer from selective constraints may provide a mechanism for developing new interactions and a pathway for the evolution of new genes.