Phylogenetic analysis of the vertebrate galectin family

Galectins form a family of structurally related carbohydrate binding proteins (lectins) that have been identified in a large variety of metazoan phyla. They are involved in many biological processes such as morphogenesis, control of cell death, immunological response, and cancer. To elucidate the evolutionary history of galectins and galectin-like proteins in chordates, we have exploited three independent lines of evidence: (i) location of galectin encoding genes (LGALS) in the human genome; (ii) exon-intron organization of galectin encoding genes; and (iii) sequence comparison of carbohydrate recognition domains (CRDs) of chordate galectins. Our results suggest that a duplication of a mono-CRD galectin gene gave rise to an original bi-CRD galectin gene, before or early in chordate evolution. The N-terminal and C-terminal CRDs of this original galectin subsequently diverged into two different subtypes, defined by exon-intron structure (F4-CRD and F3-CRD). We show that all vertebrate mono-CRD galectins known to date belong to either the F3- or F4- subtype. A sequence of duplication and divergence events of the different galectins in chordates is proposed.     

Concerted evolution of exons and introns in the MHC-linked tenascin-X gene of mammals.

Two modes of evolution of repeated domains in proteins have been described: (1) a conservative mode, whereby individual domains are conserved across gene duplication and speciation events, and (2) a concerted mode, whereby repeat domains become homogenized within a gene, presumably by intragenic partial duplication and/or gene conversion. The evolution of repeated EGF-like and fibronection-type-III-like (Fn-III) domains in the vertebrate extracellular matrix proteins tenascin-X (TNX) and tenascin-C (TNC) was studied by comparisons between human and mouse orthologs and between the paralogous TNC and TNX genes. The EGF-like repeats have largely been homogenized within each gene by concerted evolution since the duplication of the two genes but have been conserved since the divergence of rodents and primates. The Fn-III domains of TNC have likewise mainly evolved in a conservative fashion since the divergence of rodents and primates. In contrast, the Fn-III repeats of TNX fall into three distinct categories with regard to mode of evolution: (1) The three C-terminal repeats have been conserved since before duplication of the TNX and TNC genes. (2) Certain other repeats have been homogenized within each gene since gene duplication but have been conserved since the divergence of rodents and primates. (3) Still other repeats have evolved in a concerted fashion in rodent and primate lineages since their divergence. Remarkably, certain introns adjacent to the exons encoding these concertedly evolving Fn-III repeats have themselves evolved in a concerted fashion. This is the first known example of concerted evolution of repeated introns within a protein-coding gene.     

Fgfrl1, a fibroblast growth factor receptor-like gene, is found in the cephalochordate Branchiostoma floridae but not in the urochordate Ciona intestinalis

FGFRL1 is a novel member of the fibroblast growth factor receptor family that controls the formation of musculoskeletal tissues. Some vertebrates, including man, cow, dog, mouse, rat and chicken, possess a single copy the FGFRL1 gene. Teleostean fish have two copies, fgfrl1a and fgfrl1b, because they have undergone a whole genome duplication. Vertebrates belong to the chordates, a phylum that also includes the subphyla of the cephalochordates (e.g. Branchiostoma floridae) and urochordates (tunicates, e.g. Ciona intestinalis). We therefore investigated whether other chordates might also possess an FGFRL1 related gene. In fact, a homologous gene was found in B. floridae (amphioxus). The corresponding protein showed 60% sequence identity with the human protein and all sequence motifs identified in the vertebrate proteins were also conserved in amphioxus Fgfrl1. In contrast, the genome of the urochordate C. intestinalis and those from more distantly related invertebrates including the insect Drosophila melanogaster and the nematode Caenorhabditis elegans did not appear to contain any related sequences. Thus, the FGFRL1 gene might have evolved just before branching of the vertebrate lineage from the other chordates.     

Gene duplication and the origins of morphological complexity in pancrustacean eyes,a genomic approach

Background  

Duplication and divergence of genes and genetic networks is hypothesized to be a major driver of the evolution of complexity and novel features. Here, we examine the history of genes and genetic networks in the context of eye evolution by using new approaches to understand patterns of gene duplication during the evolution of metazoan genomes. We hypothesize that 1) genes involved in eye development and phototransduction have duplicated and are retained at higher rates in animal clades that possess more distinct types of optical design; and 2) genes with functional relationships were duplicated and lost together, thereby preserving genetic networks. To test these hypotheses, we examine the rates and patterns of gene duplication and loss evident in 19 metazoan genomes, including that of Daphnia pulex - the first completely sequenced crustacean genome. This is of particular interest because the pancrustaceans (hexapods+crustaceans) have more optical designs than any other major clade of animals, allowing us to test specifically whether the high amount of disparity in pancrustacean eyes is correlated with a higher rate of duplication and retention of vision genes.     

Phylogenetic and evolutionary analysis of the septin protein family in metazoan

Septins, a conserved family of cytoskeletal GTP-binding proteins, were presented in diverse eukaryotes. Here, a comprehensive phylogenetic and evolutionary analysis for septin proteins in metazoan was carried out. First, we demonstrated that all septin proteins in metazoan could be clustered into four subgroups, and the representative homologue of every subgroup was presented in the non-vertebrate chordate Ciona intestinalis, indicating that the emergence of the four septin subgroups should have occurred prior to divergence of vertebrates and invertebrates, and the expansion of the septin gene number in vertebrates was mainly by the duplication of pre-existing genes rather than by the appearance of new septin subgroup. Second, the direct orthologues of most human septins existed in zebrafish, which suggested that human septin gene repertoire was mainly formed by as far as before the split between fishes and land vertebrates. Third, we found that the evolutionary rate within septin family in mammalian lineage varies significantly, human SEPT1, SEPT 10, SEPT 12, and SEPT 14 displayed a relative elevated evolutionary rate compared with other septin members. Our data will provide new insights for the further function study of this protein family.     

Evolution after Whole-Genome Duplication: Teleost MicroRNAs

MicroRNAs (miRNAs) are important gene expression regulators implicated in many biological processes, but we lack a global understanding of how miRNA genes evolve and contribute to developmental canalization and phenotypic diversification. Whole-genome duplication events likely provide a substrate for species divergence and phenotypic change by increasing gene numbers and relaxing evolutionary pressures. To understand the consequences of genome duplication on miRNA evolution, we studied miRNA genes following the teleost genome duplication (TGD). Analysis of miRNA genes in four teleosts and in spotted gar, whose lineage diverged before the TGD, revealed that miRNA genes were retained in ohnologous pairs more frequently than protein-coding genes, and that gene losses occurred rapidly after the TGD. Genomic context influenced retention rates, with clustered miRNA genes retained more often than nonclustered miRNA genes and intergenic miRNA genes retained more frequently than intragenic miRNA genes, which often shared the evolutionary fate of their protein-coding host. Expression analyses revealed both conserved and divergent expression patterns across species in line with miRNA functions in phenotypic canalization and diversification, respectively. Finally, major strands of miRNA genes experienced stronger purifying selection, especially in their seeds and 3′-complementary regions, compared with minor strands, which nonetheless also displayed evolutionary features compatible with constrained function. This study provides the first genome-wide, multispecies analysis of the mechanisms influencing metazoan miRNA evolution after whole-genome duplication.     

Gene expansion and retention leads to a diverse tyrosine kinase superfamily in amphioxus

Tyrosine kinase (TK) proteins play a central role in cellular behavior and development of animals. The expansion of this superfamily is regarded as a key event in the evolution of the complex signaling pathways and gene networks of metazoans and is a prominent example of how shuffling of protein modules may generate molecular novelties. Using the intron/exon structure within the TK domain (TK intron code) as a complementary tool for the assignment of orthology and paralogy, we identified and studied the 118 TK proteins of the amphioxus Branchiostoma floridae genome to elucidate TK gene family evolution in metazoans and chordates in particular. Unlike all characterized metazoans to date, amphioxus has members of all known widespread TK families, with not a single loss. Putting amphioxus TKs in an evolutionary context, including new data from the cnidarian Nematostella vectensis, the echinoderm Strongylocentrotus purpuratus, and the ascidian Ciona intestinalis, we suggest new evolutionary histories for different TK families and draw a new global picture of gene loss/gain in the different phyla. Surprisingly, our survey also detected an unprecedented expansion of a group of closely related TK families, including TIE, FGFR, PDGFR, and RET, due most probably to massive gene duplication and exon shuffling. Based on their highly similar intron/exon structure at the TK domain, we suggest that this group of TK families constitute a superfamily of TK proteins, which we termed EXpanding TK, after their seemingly unique propensity to gene duplication and exon shuffling, not only in amphioxus but also across all metazoan groups. Due to this extreme tendency to both retention and expansion of TK genes, amphioxus harbors the richest and most diverse TK repertoire among all metazoans studied so far, retaining most of the gene complement of its ancestors, but having evolved its own repertoire of genetic novelties.     

Evaluation of whether accelerated protein evolution in chordates has occurred before, after, or simultaneously with gene duplication

Gene duplication and loss are predicted to be at least of the order of the substitution rate and are key contributors to the development of novel gene function and overall genome evolution. Although it has been established that proteins evolve more rapidly after gene duplication, we were interested in testing to what extent this reflects causation or association. Therefore, we investigated the rate of evolution prior to gene duplication in chordates. Two patterns emerged; firstly, branches, which are both preceded by a duplication and followed by a duplication, display an elevated rate of amino acid replacement. This is reflected in the ratio of nonsynonymous to synonymous substitution (mean nonsynonymous to synonymous nucleotide substitution rate ratio [Ka:Ks]) of 0.44 compared with branches preceded by and followed by a speciation (mean Ka:Ks of 0.23). The observed patterns suggest that there can be simultaneous alteration in the selection pressures on both gene duplication and amino acid replacement, which may be consistent with co-occurring increases in positive selection, or alternatively with concurrent relaxation of purifying selection. The pattern is largely, but perhaps not completely, explained by the existence of certain families that have elevated rates of both gene duplication and amino acid replacement. Secondly, we observed accelerated amino acid replacement prior to duplication (mean Ka:Ks for postspeciation preduplication branches was 0.27). In some cases, this could reflect adaptive changes in protein function precipitating a gene duplication event. In conclusion, the circumstances surrounding the birth of new proteins may frequently involve a simultaneous change in selection pressures on both gene-copy number and amino acid replacement. More precise modeling of the relative importance of preduplication, postduplication, and simultaneous amino acid replacement will require larger and denser genomic data sets from multiple species, allowing simultaneous estimation of lineage-specific fluctuations in mutation rates and adaptive constraints.     

Perforin evolved from a gene duplication of MPEG1, followed by a complex pattern of gene gain and loss within Euteleostomi

ABSTRACT: BACKGROUND: The pore-forming protein perforin is central to the granule-exocytosis pathway used by cytotoxic lymphocytes to kill abnormal cells. Although this mechanism of killing is conserved in bony vertebrates, cytotoxic cells are present in other chordates and invertebrates, and their cytotoxic mechanism has not been elucidated. In order to understand the evolution of this pathway, here we characterize the origins and evolution of perforin. RESULTS: We identified orthologs and homologs of human perforin in all but one species analysed from Euteleostomi, and present evidence for an earlier ortholog in Gnathostomata but not in more primitive chordates. In placental mammals perforin is a single copy gene, but there are multiple perforin genes in all lineages predating marsupials, except birds. Our comparisons of these many-to-one homologs of human perforin show that they mainly arose from lineagespecific gene duplications in multiple taxa, suggesting acquisition of new roles or different modes of regulation. We also present evidence that perforin arose from duplication of the ancient MPEG1 gene, and that it shares a common ancestor with the functionally related complement proteins. CONCLUSIONS: The evolution of perforin in vertebrates involved a complex pattern of gene, as well as intron, gain and loss. The primordial perforin gene arose at least 500 million years ago, at around the time that the major histocompatibility complex-T cell receptor antigen recognition system was established. As it is absent from primitive chordates and invertebrates, cytotoxic cells from these lineages must possess a different effector molecule or cytotoxic mechanism.     

Origins, lineage-specific expansions, and multiple losses of tyrosine kinases in eukaryotes

Tyrosine kinases are important components of metazoan signaling pathways, and their mutant forms are implicated in various malignancies. Searching the sequences from the genomes of 28 eukaryotes and the GenBank, we found tyrosine kinases not only in metazoans but also in the green algae Chlamydomonas reinhardtii, the potato late blight pathogen Phytophthora infestans, and the protozoan pathogen Entamoeba histolytica, contrary to the current view that tyrosine kinases are animal-specific. Based on a phylogenetic analysis, we divided this gene family into 43 subfamilies and found that at least 19 tyrosine kinases were likely present in the common ancestor of chordates, arthropods, and nematodes. Interestingly, most of the subfamilies have conserved domain organizations among subfamily members but have undergone different degrees of expansion during the evolution of metazoans. In particular, a large number of duplications occurred in the lineage leading to the common ancestor of Tagifugu and mammals after its split from the Ciona lineage about 450 to 550 MYA. The timing of expansion coincides with proposed large-scale duplication event in the chordate lineage. Furthermore, gene losses have occurred in most subfamilies. Interestingly, different subfamilies have similar net gain rates in the chordates studied. However, the tyrosine kinases in mouse and human or in fruit fly and mosquito mostly have a one-to-one relationship between species, indicating that static periods of 90 Myr or longer in tyrosine kinase evolution have followed large expansion events.     

Mean time to resolution of gene duplication

The mean time to resolution of gene duplication (T_r) is studied in this paper under the double null recessive (DNR) and haplo-insufficient (HI) models within the same analytical and simulation framework. We show that when population size is not too small (more precisely Nμ > 0.1), T_r for unlinked duplication is usually larger than that for linked and T_r for unlinked duplication under the HI model might be greatly prolonged, which were consistent with previous observations. Furthermore, by analytical approach we here indicate the primary underlying mechanism is that the frequency of the original (or wild-type) chromosomal haplotype of the linked duplication decreases nearly exponential to zero with time while that of the unlinked decreases quickly to an quasi-equilibrium; and this phenomenon is particularly profound under the HI model, because the quasi-equilibrium frequency of the original chromosomal haplotype (x₀) under the HI model is higher than that under the DNR model. These results suggest that recombination and HI model might jointly contribute to the marked prolongation of T_r even in a modest population. The prolonged T_r and higher quasi-equilibrium frequency of the original allele at both duplicated loci might have offered more opportunities for the emergence of novel genes.     

Phylogenetic Analysis of the NEEP21/Calcyon/P19 Family of Endocytic Proteins: Evidence for Functional Evolution in the Vertebrate CNS

Endocytosis and vesicle trafficking are required for optimal neural transmission. Yet, little is currently known about the evolution of neuronal proteins regulating these processes. Here, we report the first phylogenetic study of NEEP21, calcyon, and P19, a family of neuronal proteins implicated in synaptic receptor endocytosis and recycling, as well as in membrane protein trafficking in the somatodendritic and axonal compartments of differentiated neurons. Database searches identified orthologs for P19 and NEEP21 in bony fish, but not urochordate or invertebrate phyla. Calcyon orthologs were only retrieved from mammalian databases and distant relatives from teleost fish. In situ localization of the P19 zebrafish ortholog, and extant progenitor of the gene family, revealed a CNS specific expression pattern. Based on non-synonymous nucleotide substitution rates, the calcyon genes appear to be under less intense negative selective pressure. Indeed, a functional group II WW domain binding motif was detected in primate and human calcyon, but not in non-primate orthologs. Sequencing of the calcyon gene from 80 human subjects revealed a non-synonymous single nucleotide polymorphism that abrogated group II WW domain protein binding. Altogether, our data indicate the NEEP21/calcyon/P19 gene family emerged, and underwent two rounds of gene duplication relatively late in metazoan evolution (but early in vertebrate evolution at the latest). As functional studies suggest NEEP21 and calcyon play related, but distinct roles in regulating vesicle trafficking at synapses, and in neurons in general, we propose the family arose in chordates to support a more diverse range of synaptic and behavioral responses.     

Phylogenetic analysis of three complete gap junction gene families reveals lineage-specific duplications and highly supported gene classes

Gap junctions, composed of connexin proteins in chordates, are the most ubiquitous form of intercellular communication. Complete connexin gene families have been identified from human (20) and mouse (19), revealing significant diversity in gap junction channels. We searched current databases and identified 37 putative zebrafish connexin genes, almost twice the number found in mammals. Phylogenetic comparison of entire connexin gene families from human, mouse, and zebrafish revealed 23 zebrafish relatives of 16 mammalian connexins, and 14 connexins apparently unique to zebrafish. We found evidence for duplication events in all genomes, as well as evidence for recent tandem duplication events in the zebrafish, indicating that the complexity of the connexin family is growing. The identification of a third complete connexin gene family provides novel insight into the evolution of connexins, and sheds light into the phenotypic evolution of intercellular communication via gap junctions.     

Identification and expression analysis of BMP signaling inhibitors genes of the DAN family in amphioxus

Bone morphogenetic proteins (BMPs) are members of the Transforming Growth Factor-β (TGF-β) family implicated in many developmental processes in metazoans such as embryo axes specification. Their wide variety of actions is in part controlled by inhibitors that impede the interaction of BMPs with their specific receptors. Here, we focused our attention on the Differential screening-selected gene Aberrative in Neuroblastoma (DAN) family of inhibitors. Although they are well-characterized in vertebrates, few data are available for this family in other metazoan species. In order to understand the evolution of potential developmental roles of these inhibitors in chordates, we identified the members of this family in the cephalochordate amphioxus, and characterized their expression patterns during embryonic development. Our data suggest that the function of Cerberus/Dand5 subfamily genes is conserved among chordates, whereas Gremlin1/2 and NBL1 subfamily genes seem to have acquired divergent expression patterns in each chordate lineage. On the other hand, the expression of Gremlin in the amphioxus neural plate border during early neurulation strengthens the hypothesis of a conserved neural plate border gene network in chordates.     

Tunicates Illuminate the Enigmatic Evolution of Chordate Metallothioneins by Gene Gains and Losses,Independent Modular Expansions,and Functional Convergences

To investigate novel patterns and processes of protein evolution, we have focused in the metallothioneins (MTs), a singular group of metal-binding, cysteine-rich proteins that, due to their high degree of sequence diversity, still represents a “black hole” in Evolutionary Biology. We have identified and analyzed more than 160 new MTs in nonvertebrate chordates (especially in 37 species of ascidians, 4 thaliaceans, and 3 appendicularians) showing that prototypic tunicate MTs are mono-modular proteins with a pervasive preference for cadmium ions, whereas vertebrate and cephalochordate MTs are bimodular proteins with diverse metal preferences. These structural and functional differences imply a complex evolutionary history of chordate MTs—including de novo emergence of genes and domains, processes of convergent evolution, events of gene gains and losses, and recurrent amplifications of functional domains—that would stand for an unprecedented case in the field of protein evolution.     

Evolution and diversity of axon guidance Robo receptor family genes

The diversity of axon guidance (AG) receptors reflects gains in complexity of the animal nervous system during evolution. Members of the Roundabout (Robo) family of receptors interact with Slit proteins and play important roles in many developmental processes, including AG and neural crest cell migration. There are four members of the Robo gene family. However, the evolutionary history of Robo family genes remain obscure. We analyzed the distribution of Robo family members in metazoan species ranging in complexity from hydras to humans. We undertook a phylogenetic analysis in metazoans, synteny analysis, and ancestral chromosome mapping in vertebrates, and detected selection pressure and functional divergence among four mammalian Robo paralogs. Based on our analysis, we proposed that the ancestral Robo gene could have undergone a tandem duplication in the vertebrate ancestor; then one round of whole genome duplication events occurred before the divergence of ancestral lamprey and gnathostome, generating four paralogs in early vertebrates. Robo4 paralog underwent segmental loss in the following evolutionary process. Our results showed that Robo3 paralog is under more powerful purifying selection pressure compared with other three paralogs, which could correlate with its unique expression pattern and function. Furthermore, we found four sites under positive selection pressure on the Ig1‐2 domains of Robo4 that might interfere with its binding to Slits ligand. Diverge analysis at the amino acid level showed that Robo4 paralog have relatively greater functional diversifications than other Robo paralogs. This coincides with the fact that Robo4 predominantly functions in vascular endothelial cells but not the nervous system.     

Comparative genomics in hemiascomycete yeasts: evolution of sex, silencing, and subtelomeres

The recent release of sequences of several unexplored yeast species that cover an evolutionary range comparable to the entire phylum of chordates offers us a unique opportunity to investigate how genes involved in adaptation have been shaped by evolution. We have examined how three different sets of genes, all related to adaptative processes at the genomic level, have evolved in hemiascomycetes: (1) the mating-type genes that govern sexuality, (2) the silencing genes that are connected to regulation of mating-type cassettes and to telomere position effect, and (3) the gene families found repeated in subtelomeric regions.We report new combinations of mating-type genes and cassettes in hemiascomycetous species; we show that silencing proteins diverge rapidly. We have also found that in all species studied, subtelomeric gene families exist and are specific to each species.     

Expansion of protein domain repeats

Many proteins, especially in eukaryotes, contain tandem repeats of several domains from the same family. These repeats have a variety of binding properties and are involved in protein–protein interactions as well as binding to other ligands such as DNA and RNA. The rapid expansion of protein domain repeats is assumed to have evolved through internal tandem duplications. However, the exact mechanisms behind these tandem duplications are not well-understood. Here, we have studied the evolution, function, protein structure, gene structure, and phylogenetic distribution of domain repeats. For this purpose we have assigned Pfam-A domain families to 24 proteomes with more sensitive domain assignments in the repeat regions. These assignments confirmed previous findings that eukaryotes, and in particular vertebrates, contain a much higher fraction of proteins with repeats compared with prokaryotes. The internal sequence similarity in each protein revealed that the domain repeats are often expanded through duplications of several domains at a time, while the duplication of one domain is less common. Many of the repeats appear to have been duplicated in the middle of the repeat region. This is in strong contrast to the evolution of other proteins that mainly works through additions of single domains at either terminus. Further, we found that some domain families show distinct duplication patterns, e.g., nebulin domains have mainly been expanded with a unit of seven domains at a time, while duplications of other domain families involve varying numbers of domains. Finally, no common mechanism for the expansion of all repeats could be detected. We found that the duplication patterns show no dependence on the size of the domains. Further, repeat expansion in some families can possibly be explained by shuffling of exons. However, exon shuffling could not have created all repeats.     

Expression of engrailed proteins in arthropods, annelids, and chordates

engrailed is a homeobox gene that has an important role in Drosophila segmentation. Genes homologous to engrailed have been identified in several other organisms. Here we describe a monoclonal antibody that recognizes a conserved epitope in the homeodomain of engrailed proteins of a number of different arthropods, annelids, and chordates; we use this antibody to isolate the grasshopper engrailed gene. In Drosophila embryos, the antibody reveals engrailed protein in the posterior portion of each segment during segmentation, and in a segmentally reiterated subset of neuronal cells during neurogenesis. Other arthropods, including grasshopper and two crustaceans, have similar patterns of engrailed expression. However, these patterns of expression are not shared by the annelids or chordates we examined. Our results provide the most comprehensive view that has been obtained of how expression patterns of a regulatory gene vary during evolution. On the basis of these patterns, we suggest that engrailed is a gene whose ancestral function was in neurogenesis and whose function was co-opted during the evolution of segmentation in the arthropods, but not in the annelids and chordates.