Protein Tyrosine Phosphatases from Amphioxus, Hagfish, and Ray: Divergence of Tissue-Specific Isoform Genes in the Early Evolution of Vertebrates

Since separation from fungi and plants, multicellular animals evolved a variety of gene families involved in cell-cell communication from a limited number of ancestral precursors by gene duplications in two separate periods of animal evolution. In the very early evolution of animals before the separation of parazoans and eumetazoans, animals underwent extensive gene duplications by which different subtypes (subfamilies) with distinct functions diverged. The multiplicity of members (isoforms) in the same subtype increased by further gene duplications (isoform duplications) in the first half of chordate evolution before the fish-tetrapod split; different isoforms are virtually identical in structure and function but differ in tissue distribution. From cloning and phylogenetic analyses of four subfamilies of the protein tyrosine kinase (PTK) family, we recently showed extensive isoform duplications in a limited period around or just before the cyclostome-gnathostome split. To obtain a reliable estimate for the divergence time of vertebrate isoforms, we have conducted isolation of cDNAs encoding the protein tyrosine phosphatases (PTPs) from Branchiostoma belcheri, an amphioxus, Eptatretus burgeri, a hagfish, and Potamotrygon motoro, a ray. We obtained 33 different cDNAs in total, most of which belong to known PTP subfamilies. The phylogenetic analyses of five subfamilies based on the maximum likelihood method revealed frequent isoform duplications in a period around or just before the gnathostome-cyclostome split. An evolutionary implication was discussed in relation to the Cambrian explosion.     

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     

The monosaccharide transporter gene family in Arabidopsis and rice: a history of duplications, adaptive evolution, and functional divergence

Current hypotheses of gene duplicate divergence propose that surviving members of a gene duplicate pair may evolve, under conditions of purifying or nearly neutral selection, in one of two ways: with new function arising in one duplicate while the other retains original function (neofunctionalization [NF]) or partitioning of the original function between the 2 paralogs (subfunctionalization [SF]). More recent studies propose that SF followed by NF (subneofunctionalization [SNF]) explains the divergence of many duplicate genes. In this analysis, we evaluate these hypotheses in the context of the large monosaccharide transporter (MST) gene families in Arabidopsis and rice. MSTs have an ancient origin, predating plants, and have evolved in the seed plant lineage to comprise 7 subfamilies. In Arabidopsis, 53 putative MST genes have been identified, with one subfamily greatly expanded by tandem gene duplications. We searched the rice genome for members of the MST gene family and compared them with the MST gene family in Arabidopsis to determine subfamily expansion patterns and estimate gene duplicate divergence times. We tested hypotheses of gene duplicate divergence in 24 paralog pairs by comparing protein sequence divergence rates, estimating positive selection on codon sites, and analyzing tissue expression patterns. Results reveal the MST gene family to be significantly larger (65) in rice with 2 subfamilies greatly expanded by tandem duplications. Gene duplicate divergence time estimates indicate that early diversification of most subfamilies occurred in the Proterozoic (2500-540 Myr) and that expansion of large subfamilies continued through the Cenozoic (65-0 Myr). Two-thirds of paralog pairs show statistically symmetric rates of sequence evolution, most consistent with the SF model, with half of those showing evidence for positive selection in one or both genes. Among 8 paralog pairs showing asymmetric divergence rates, most consistent with the NF model, nearly half show evidence of positive selection. Positive selection does not appear in any duplicate pairs younger than approximately 34 Myr. Our data suggest that the NF, SF, and SNF models describe different outcomes along a continuum of divergence resulting from initial conditions of relaxed constraint after duplication.     

Multiple receptor-like kinase cDNAs from liverwort Marchantia polymorpha and two charophycean green algae, Closterium ehrenbergii and Nitella axillaris: Extensive gene duplications and gene shufflings in the early evolution of streptophytes

Plant receptor-like kinases (RLKs) comprise a large family with more than several hundred members in vascular plants. The RLK family is thought to have diverged specifically in the plant kingdom, and no family member has been identified in other lineages except for animals and Plasmodium, both of which have RLK related families of small size. To know the time of divergence of RLK family members by gene duplications and domain shufflings, comprehensive isolations of RLK cDNAs were performed from a nonvascular plant, liverwort Marchantia polymorpha and two charophycean green algae, Closterium ehrenbergii, and Nitella axillaris, thought to be the closest relatives to land plants. We obtained twenty-nine, fourteen, and thirteen RLK related cDNAs from M. polymorpha, C. ehrenbergii, and N. axillaris, respectively. The amino acid sequences of these RLKs were compared with those of vascular plants, and phylogenetic trees were inferred by GAMT, a genetic algorithm-based maximum likelihood (ML) method that outputs multiple trees, together with best one. The inferred ML trees revealed ancient gene duplications generating subfamilies with different domain organizations, which occurred extensively at least before the divergence of vascular and nonvascular plants. Rather it remains possible that the extensive gene duplications occurred during the early evolution of streptophytes. Multicellular-specific somatic embryogenesis receptor kinase (SERK) involved in somatic embryogenesis was found in a unicellular alga C. ehrenbergii, suggesting the evolution of SERK by gene recruitment of a unicellular gene.     

The Evolution of the Calpain Family as Reflected in Paralogous Chromosome Regions

Calpains, the Ca²⁺-dependent intracellular proteinases, are involved in the regulation of distinct cellular pathways including signal transduction and processing, cytoskeleton dynamics, and muscle homeostasis. To investigate the evolutionary origin of diverse calpain subfamilies, a phylogenetic study was carried out. The topology of the calpain phylogenetic tree has shown that some of the gene duplications occurred before the divergence of the protostome and deuterostome lineages. Other gene doublings, leading to vertebrate-specific calpain forms, took place during early chordate evolution and coincided with genome duplications as disclosed by the localization of calpain genes to paralogous chromosome regions in the human genome. On the basis of the phylogenetic tree, the time of gene duplications, and the localization of calpain genes, we propose a model of tandem and chromosome duplications for the evolution of vertebrate-specific calpain forms. The data presented here are consistent with scenarios proposed for the evolution of other multigene families. Received: 17 November 1998 / Accepted: 30 April 1999     

Kinesin-related genes from diplomonad,sponge, amphioxus,and cyclostomes: divergence pattern of kinesin family and evolution of giardial membrane-bounded organella

To understand the question of whether divergence of eukaryotic genes by gene duplications and domain shufflings proceeded gradually or intermittently during evolution, we have cloned and sequenced Giardia lamblia cDNAs encoding kinesins and kinesin-related proteins and have obtained 13 kinesin-related cDNAs, some of which are likely homologs of vertebrate kinesins involved in vesicle transfer to ER, Golgi, and plasma membrane. A phylogenetic tree of the kinesin family revealed that most gene duplications that gave rise to different kinesin subfamilies with distinct functions have been completed before the earliest divergence of extant eukaryotes. This suggests that the complex endomembrane system has arisen very early in eukaryotic evolution, and the diminutive ER and Golgi apparatus recognized in the giardial cells, together with the absence of mitochondria, might be characters acquired secondarily during the evolution of parasitism. To understand the divergence pattern of the kinesin family in the lineage leading to vertebrates, seven more Unc104-related cDNAs have been cloned from sponge, amphioxus, hagfish, and lamprey. The divergence pattern of the animal Unc104/KIF1 subfamily is characterized by two active periods in gene duplication interrupted by a considerably long period of silence, instead of proceeding gradually: animals underwent extensive gene duplications before the parazoan-eumetazoan split. In the early evolution of vertebrates around the cyclostome-gnathostome split, further gene duplications occurred, by which a variety of genes with similar structures over the entire regions were generated. This pattern of divergence is similar to those of animal genes involved in cell-cell communication and developmental control.     

The temporal distribution of gene duplication events in a set of highly conserved human gene families

Using a data set of protein translations associated with map positions in the human genome, we identified 1520 mapped highly conserved gene families. By comparing sharing of families between genomic windows, we identified 92 potentially duplicated blocks in the human genome containing 422 duplicated members of these families. Using branching order in the phylogenetic trees, we timed gene duplication events in these families relative to the primate-rodent divergence, the amniote-amphibian divergence, and the deuterostome-protostome divergence. The results showed similar patterns of gene duplication times within duplicated blocks and outside duplicated blocks. Both within and outside duplicated blocks, numerous duplications were timed prior to the deuterostome-protostome divergence, whereas others occurred after the amniote-amphibian divergence. Thus, neither gene duplication in general nor duplication of genomic blocks could be attributed entirely to polyploidization early in vertebrate history. The strongest signal in the data was a tendency for intrachromosomal duplications to be more recent than interchromosomal duplications, consistent with a model whereby tandem duplication-whether of single genes or of genomic blocks-may be followed by eventual separation of duplicates due to chromosomal rearrangements. The rate of separation of tandemly duplicated gene pairs onto separated chromosomes in the human lineage was estimated at 1.7 x 10(-9) per gene-pair per year.     

Phylogenetic analysis and evolution of aromatic amino acid hydroxylase

A study was performed to investigate the phylogenetic relationship among AAAH members and to statistically evaluate sequence conservation and functional divergence. In total, 161 genes were identified from 103 species. Phylogenetic analysis showed that well-conserved subfamilies exist. Exon-intron structure analysis showed that the gene structures of AAAH were highly conserved across some different lineage species, while some species-specific introns were also found. The dynamic distribution of ACT domain suggested one gene fusion event has occurred in eukaryota. Significant functional divergence was found between some subgroups. Analysis of the site-specific profiles revealed critical amino acid residues for functional divergence. This study highlights the molecular evolution of this family and may provide a starting point for further experimental verifications.     

Extensive duplications of phototransduction genes in early vertebrate evolution correlate with block (chromosome) duplications

Many gene families in mammals have members that are expressed more or less uniquely in the retina or differentially in specific retinal cell types. We describe here analyses of nine such gene families with regard to phylogenetic relationships and chromosomal location. The families are opsins, G proteins (alpha, beta, and gamma subunits), phosphodiesterases type 6, cyclic nucleotide-gated channels, G-protein-coupled receptor kinases, arrestins, and recoverins. The results suggest that multiple new gene copies arose in all of these families very early in vertebrate evolution during a period with extensive gene duplications. Many of the new genes arose through duplications of large chromosome regions (blocks of genes) or even entire chromosomes, as shown by linkage with other gene families. Some of the phototransduction families belong to the same duplicated regions and were thus duplicated simultaneously. We conclude that gene duplications in early vertebrate evolution probably helped facilitate the specialization of the retina and the subspecialization of different retinal cell types.     

Species-specific size expansion and molecular evolution of the oleosins in angiosperms

Oleosins are hydrophobic plant proteins thought to be important for the formation of oil bodies, which supply energy for seed germination and subsequent seedling growth. To better understand the evolutionary history and diversity of the oleosin gene family in plants, especially angiosperms, we systematically investigated the molecular evolution of this family using eight representative angiosperm species. A total of 73 oleosin members were identified, with six members in each of four monocot species and a greater but variable number in the four eudicots. A phylogenetic analysis revealed that the angiosperm oleosin genes belonged to three monophyletic lineages. Species-specific gene duplications, caused mainly by segmental duplication, led to the great expansion of oleosin genes and occurred frequently in eudicots after the monocot–eudicot divergence. Functional divergence analyses indicate that significant amino acid site-specific selective constraints acted on the different clades of oleosins. Adaptive evolution analyses demonstrate that oleosin genes were subject to strong purifying selection after their species-specific duplications and that rapid evolution occurred with a high degree of evolutionary dynamics in the pollen-specific oleosin genes. In conclusion, this study serves as a foundation for genome-wide analyses of the oleosins. These findings provide insight into the function and evolution of this gene family in angiosperms and pave the way for studies in other plants.     

Duplication and adaptive evolution of the COR15 genes within the highly cold-tolerant Draba lineage (Brassicaceae)

Plants have evolved diverse adaptive mechanisms that enable them to tolerate abiotic stresses, to varying degrees, and such stresses may have strongly influenced evolutionary changes at levels ranging from molecular to morphological. Previous studies on these phenomena have focused on the adaptive evolution of stress-related orthologous genes in specific lineages. However, heterogenetic evolution of the paralogous genes following duplication has only been examined in a very limited number of stress-response gene families. The COR15 gene encodes a low molecular weight protein that plays an important role in protecting plants from cold stresses. Although two different copies of this gene have been found in the model species, Arabidopsis thaliana, evolutionary patterns of this small gene family in plants have not been previously explored. In this study, we cloned COR15-like sequences and performed evolutionary analyses of these sequences (including those previously reported) in the highly cold-tolerant Draba lineage and related lineages of Brassicaceae. Our phylogenetic analyses indicate that all COR15-like sequences clustered into four clades that corresponded well to the morphological lineages. Gene conversions were found to have probably occurred before/during the divergence of Brassica and Draba lineage. However, repeated, independent duplications of this gene have occurred in different lineages of Brassicaceae. Further comparisons of all sequences suggest that there have been significant inter-lineage differences in evolutionary rates between the duplicated and original genes. We assessed the likelihood that the differences between two well-supported gene subfamilies that appear to have originated from a single duplication, COR15a and COR15b, within the Draba lineage have been driven by adaptive evolution. Comparisons of their non-synonymous/synonymous substitution ratios and rates of predicted amino acid changes indicate that these two gene groups are evolving under different selective pressures and may be functionally divergent. This functional divergence was confirmed by comparing site-specific shifts in evolution indexes of the two groups of predicted proteins. The evidence of differential selection and possible functional divergence suggests that the duplication may be of adaptive significance, with possible implications for the explosive diversification of the Draba lineage during the cooling Quaternary stages and the following worldwide colonization of arid alpine and artic regions.     

Analysis of lamprey and hagfish genes reveals a complex history of gene duplications during early vertebrate evolution

It has been proposed that two events of duplication of the entire genome occurred early in vertebrate history (2R hypothesis). Several phylogenetic studies with a few gene families (mostly Hox genes and proteins from the MHC) have tried to confirm these polyploidization events. However, data from a single locus cannot explain the evolutionary history of a complete genome. To study this 2R hypothesis, we have taken advantage of the phylogenetic position of the lamprey to study the history of gene duplications in vertebrates. We selected most gene families that contain several paralogous genes in vertebrates and for which lamprey genes and an out-group are known in databases. In addition, we isolated members of the nuclear receptor superfamily in lamprey. Hagfish genes were also analyzed and found to confirm the lamprey gene analysis. Consistent with the 2R hypothesis, the phylogenetic analysis of 33 selected gene families, dispersed through the whole genome, revealed that one period of gene duplication arose before the lamprey-gnathostome split and this was followed by a second period of gene duplication after the lamprey-gnathostome split. Nevertheless, our analysis suggests that numerous gene losses and other gene-genome duplications occurred during the evolution of the vertebrate genomes. Thus, the complexity of all the paralogy groups present in vertebrates should be explained by the contribution of genome duplications (2R hypothesis), extra gene duplications, and gene losses.     

Origin and evolution of eukaryotic chaperonins: phylogenetic evidence for ancient duplications in CCT genes

Chaperonins are oligomeric protein-folding complexes which are divided into two distantly related structural classes. Group I chaperonins (called GroEL/cpn60/hsp60) are found in bacteria and eukaryotic organelles, while group II chaperonins are present in archaea and the cytoplasm of eukaryotes (called CCT/TriC). While archaea possess one to three chaperonin subunit-encoding genes, eight distinct CCT gene families (paralogs) have been characterized in eukaryotes. We are interested in determining when during eukaryotic evolution the multiple gene duplications producing the CCT subunits occurred. We describe the sequence and phylogenetic analysis of five CCT genes from TRICHOMONAS: vaginalis and seven from GIARDIA: lamblia, representatives of amitochondriate protist lineages thought to have diverged early from other eukaryotes. Our data show that the gene duplications producing the eight CCT paralogs took place prior to the organismal divergence of TRICHOMONAS: and GIARDIA: from other eukaryotes. Thus, these divergent protists likely possess completely hetero-oligomeric CCT complexes like those in yeast and mammalian cells. No close phylogenetic relationship between the archaeal chaperonins and specific CCT subunits was observed, suggesting that none of the CCT gene duplications predate the divergence of archaea and eukaryotes. The duplications producing the CCTdelta and CCTepsilon subunits, as well as CCTalpha, CCTbeta, and CCTeta, are the most recent in the CCT gene family. Our analyses show significant differences in the rates of evolution of archaeal chaperonins compared with the eukaryotic CCTs, as well as among the different CCT subunits themselves. We discuss these results in light of current views on the origin, evolution, and function of CCT complexes.     

Sponge homologs of vertebrate protein tyrosine kinases and frequent domain shufflings in the early evolution of animals before the parazoan–eumetazoan split

The protein tyrosine kinases (PTKs) diverged specifically in animal lineages by gene duplications and domain shufflings to form a large protein family comprising diverse subfamilies with distinct domain organizations and functions. On the basis of a phylogenetic tree inferred from a comparison of the shared kinase domains, we previously showed that gene duplications that gave rise to diverse subfamilies predate the divergence of parazoans and eumetazoans. There is, however, still a possibility that, although the kinase domain duplications are ancient events, the domain shufflings that gave rise to different subfamilies with distinct domain organization are more recent event than the kinase domain duplications. To clarify this problem, we have determined the complete sequences of 15 sponge PTKs and have compared the domain organizations of these sponge PTKs and those of eumetazoans. For each of ten sponge PTKs out of 15 analyzed here, a possible eumetazoan (human and Drosophila) ortholog has been identified. The sponge and eumetazoan orthologs are virtually identical in domain organization and belong to the same subfamily in the PTK family tree for each of ten orthologous pairs, except for one subfamily in which a considerable deletions and/or insertions of domains are observed. This result suggests that most, if not all, of the domain shufflings, together with gene duplications, are very old, going back to dates before the parazoan–eumetazoan split, the earliest divergence among extant animal phyla.     

A genetic model for the evolution of the glycosphingolipids

The glycosphingolipids have been found in many animal tissues, but the complexity of their molecular structure varies considerably among the different phyla. Relatively simple structures have been found in invertebrate species, while the most complex have been demonstrated in brain tissue of modern fishes and amphibians. The data on the phylogenetic distribution of the glycosphingolipids has been interpreted to indicate that a significant number of gene duplications, involving many different structural genes, may have occurred during a few specific periods of vertebrate evolution. The transition from invertebrate to jawless vertebrate, the divergence of rays and skates from true sharks, the advent of modern bony fishes and the transition from aquatic to terrestrial vertebrates, each could have veen accompained by duplications of genes involved in the synthesis and degradation of glycosphingolipids. The evolutionary study of such a multi-enzyme system may be one means to detect alterations in the genome as a whole. The apparent correspondence in time of these gene duplications involved in glycosphingolipid metabolism and periods of rapid vertebrate evolution which may have been accompanied by significant increases in the amount of cellular DNA suggests that such changes may have occurred via the mechanism of tetraploidization.     

Complex evolutionary history and diverse domain organization of SET proteins suggest divergent regulatory interactions

? Plants and animals possess very different developmental processes, yet share conserved epigenetic regulatory mechanisms, such as histone modifications. One of the most important forms of histone modification is methylation on lysine residues of the tails, carried out by members of the SET protein family, which are widespread in eukaryotes. ? We analyzed molecular evolution by comparative genomics and phylogenetics of the SET genes from plant and animal genomes, grouping SET genes into several subfamilies and uncovering numerous gene duplications, particularly in the Suv, Ash, Trx and E(z) subfamilies. ? Domain organizations differ between different subfamilies and between plant and animal SET proteins in some subfamilies, and support the grouping of SET genes into seven main subfamilies, suggesting that SET proteins have acquired distinctive regulatory interactions during evolution. We detected evidence for independent evolution of domain organization in different lineages, including recruitment of new domains following some duplications. ? More recent duplications in both vertebrates and land plants are probably the result of whole-genome or segmental duplications. The evolution of the SET gene family shows that gene duplications caused by segmental duplications and other mechanisms have probably contributed to the complexity of epigenetic regulation, providing insights into the evolution of the regulation of chromatin structure.     

Using gene-history and expression analyses to assess the involvement of LGI genes in human disorders

Mutations in the leucine-rich, glioma-inactivated 1 gene, LGI1, cause autosomal-dominant lateral temporal lobe epilepsy via unknown mechanisms. LGI1 belongs to a subfamily of leucine-rich repeat genes comprising four members (LGI1-LGI4) in mammals. In this study, both comparative developmental as well as molecular evolutionary methods were applied to investigate the evolution of the LGI gene family and, subsequently, of the functional importance of its different gene members. Our phylogenetic studies suggest that LGI genes evolved early in the vertebrate lineage. Genetic and expression analyses of all five zebrafish lgi genes revealed duplications of lgi1 and lgi2, each resulting in two paralogous gene copies with mostly nonoverlapping expression patterns. Furthermore, all vertebrate LGI1 orthologs experience high levels of purifying selection that argue for an essential role of this gene in neural development or function. The approach of combining expression and selection data used here exemplarily demonstrates that in poorly characterized gene families a framework of evolutionary and expression analyses can identify those genes that are functionally most important and are therefore prime candidates for human disorders.     

Comparable Rates of Gene Loss and Functional Divergence after Genome Duplications Early in Vertebrate Evolution

Duplicated genes are an important source of new protein functions and novel developmental and physiological pathways. Whereas most models for fate of duplicated genes show that they tend to be rapidly lost, models for pathway evolution suggest that many duplicated genes rapidly acquire novel functions. Little empirical evidence is available, however, for the relative rates of gene loss vs. divergence to help resolve these contradictory expectations. Gene families resulting from genome duplications provide an opportunity to address this apparent contradiction. With genome duplication, the number of duplicated genes in a gene family is at most 2(n), where n is the number of duplications. The size of each gene family, e.g., 1, 2, 3, . . . , 2(n), reflects the patterns of gene loss vs. functional divergence after duplication. We focused on gene families in humans and mice that arose from genome duplications in early vertebrate evolution and we analyzed the frequency distribution of gene family size, i.e., the number of families with two, three or four members. All the models that we evaluated showed that duplicated genes are almost as likely to acquire a new and essential function as to be lost through acquisition of mutations that compromise protein function. An explanation for the unexpectedly high rate of functional divergence is that duplication allows genes to accumulate more neutral than disadvantageous mutations, thereby providing more opportunities to acquire diversified functions and pathways.