Phylogenetic analysis and positive-selection site detecting of vascular endothelial growth factor family in vertebrates

Vascular endothelial growth factor (VEGF), known to play an important role in vascular homeostasis, vascular integrity and angiogenesis, is little known about the evolutionary relationship of its five members especially the role of gene duplication and natural selection in the evolution of the VEGF family. In this study, seventy-five full-length cDNA sequences from 33 vertebrate species were extracted from the NCBI's GenBank, UniProt protein database and the Ensembl database. By phylogenetic analyses, we investigated the origin, conservation, and evolution of the VEGFs. Five VEGF family members in vertebrates might be formed by gene duplication. The inferred evolutionary transitions that separate members which belong to different gene clusters correlated with changes in functional properties. Selection analysis and protein structure analysis were combined to explain the relationship of the site-specific evolution in the vertebrate VEGF family. Eleven positive selection sites, one transmembrane region and the active sites were detected in this process.     

Evolutionary history of the vertebrate mitogen activated protein kinases family

Background

The mitogen activated protein kinases (MAPK) family pathway is implicated in diverse cellular processes and pathways essential to most organisms. Its evolution is conserved throughout the eukaryotic kingdoms. However, the detailed evolutionary history of the vertebrate MAPK family is largely unclear.

Methodology/Principal Findings

The MAPK family members were collected from literatures or by searching the genomes of several vertebrates and invertebrates with the known MAPK sequences as queries. We found that vertebrates had significantly more MAPK family members than invertebrates, and the vertebrate MAPK family originated from 3 progenitors, suggesting that a burst of gene duplication events had occurred after the divergence of vertebrates from invertebrates. Conservation of evolutionary synteny was observed in the vertebrate MAPK subfamilies 4, 6, 7, and 11 to 14. Based on synteny and phylogenetic relationships, MAPK12 appeared to have arisen from a tandem duplication of MAPK11 and the MAPK13-MAPK14 gene unit was from a segmental duplication of the MAPK11-MAPK12 gene unit. Adaptive evolution analyses reveal that purifying selection drove the evolution of MAPK family, implying strong functional constraints of MAPK genes. Intriguingly, however, intron losses were specifically observed in the MAPK4 and MAPK7 genes, but not in their flanking genes, during the evolution from teleosts to amphibians and mammals. The specific occurrence of intron losses in the MAPK4 and MAPK7 subfamilies might be associated with adaptive evolution of the vertebrates by enhancing the gene expression level of both MAPK genes.

Conclusions/Significance

These results provide valuable insight into the evolutionary history of the vertebrate MAPK family.     

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.     

Molecular Evolution and Structural Features of IRAK Family Members

The interleukin-1 receptor-associated kinase (IRAK) family comprises critical signaling mediators of the TLR/IL-1R signaling pathways. IRAKs are Ser/Thr kinases. There are 4 members in the vertebrate genome (IRAK1, IRAK2, IRAKM, and IRAK4) and an IRAK homolog, Pelle, in insects. IRAK family members are highly conserved in vertebrates, but the evolutionary relationship between IRAKs in vertebrates and insects is not clear. To investigate the evolutionary history and functional divergence of IRAK members, we performed extensive bioinformatics analysis. The phylogenetic relationship between IRAK sequences suggests that gene duplication events occurred in the evolutionary lineage, leading to early vertebrates. A comparative phylogenetic analysis with insect homologs of IRAKs suggests that the Tube protein is a homolog of IRAK4, unlike the anticipated protein, Pelle. Furthermore, the analysis supports that an IRAK4-like kinase is an ancestral protein in the metazoan lineage of the IRAK family. Through functional analysis, several potentially diverged sites were identified in the common death domain and kinase domain. These sites have been constrained during evolution by strong purifying selection, suggesting their functional importance within IRAKs. In summary, our study highlighted the molecular evolution of the IRAK family, predicted the amino acids that contributed to functional divergence, and identified structural variations among the IRAK paralogs that may provide a starting point for further experimental investigations.     

Molecular evolution of proglucagon

The vertebrate proglucagon gene encodes glucagon, and the two glucagon-like peptides GLP-1 and GLP-2. To better understand the origin and diversification of the distinct hormonal roles of the three glucagon-like sequences encoded by the proglucagon gene, we have examined the evolution of this gene. The structure of proglucagon has been largely maintained within vertebrates. Duplication of the proglucagon gene or duplications of sequences within the proglucagon gene are rare. All proglucagon gene duplications are likely to be the result of genome duplication events. Examination of the rates of amino acid sequence evolution of each hormone reveals that they have not evolved in a uniform manner. Each hormone has evolved in an episodic fashion, suggesting that the selective constraints acting upon the sequence vary between, and within, vertebrate classes. Changes in selection on a sequence often reflect changes in the function of the sequence, such as the change in function of GLP-1 from a glucagon-like hormone in fish to an incretin in mammals. We found that the GLP-2 sequence underwent rapid sequence evolution in the early mammal lineage, therefore we have concluded that mammalian GLP-2 has acquired a new biological function that is not found in other vertebrates. Comparisons of the hormone sequences show that many amino acid residues that are functionally important in mammalian hormones are not conserved through vertebrate evolution. This observation suggests that the sequences involved in hormone action change through evolution.     

Comparative genomic analysis reveals evolutionary characteristics and patterns of microRNA clusters in vertebrates

MicroRNAs (miRNAs) are a class of small non-coding RNAs that can play important regulatory roles in many important biological processes. Although clustering patterns of miRNA clusters have been uncovered in animals, the origin and evolution of miRNA clusters in vertebrates are still poorly understood. Here, we performed comparative genomic analyses to construct 51 sets of orthologous miRNA clusters (SOMCs) across seven test vertebrate species, a collection of miRNA clusters from two or more species that are likely to have evolved from a common ancestral miRNA cluster, and used these to systematically examine the evolutionary characteristics and patterns of miRNA clusters in vertebrates. We found that miRNA clusters are continuously generated, and most of them tend to be conserved and maintained in vertebrate genomes, although some adaptive gains and losses of miRNA cluster have occurred during evolution. Furthermore, miRNA clusters appeared relatively early in the evolutionary history might suffer from more complicated adaptive gain-and-loss than those young miRNA clusters. Detailed analysis showed that genomic duplication events of ancestral miRNAs or miRNA clusters are likely to be major driving force and apparently contribute to origin and evolution of miRNA clusters. Comparison of conserved with lineage-specific miRNA clusters revealed that the contribution of duplication events for the formation of miRNA cluster appears to be more important for conserved miRNA clusters than lineage-specific. Our study provides novel insights for further exploring the origins and evolution of miRNA clusters in vertebrates at a genome scale.     

Evolutionary growth process of highly conserved sequences in vertebrate genomes

Genome sequence comparison between evolutionarily distant species revealed ultraconserved elements (UCEs) among mammals under strong purifying selection. Most of them were also conserved among vertebrates. Because they tend to be located in the flanking regions of developmental genes, they would have fundamental roles in creating vertebrate body plans. However, the evolutionary origin and selection mechanism of these UCEs remain unclear. Here we report that UCEs arose in primitive vertebrates, and gradually grew in vertebrate evolution. We searched for UCEs in two teleost fishes, Tetraodon nigroviridis and Oryzias latipes, and found 554 UCEs with 100% identity over 100 bps. Comparison of teleost and mammalian UCEs revealed 43 pairs of common, jawed-vertebrate UCEs (jUCE) with high sequence identities, ranging from 83.1% to 99.2%. Ten of them retain lower similarities to the Petromyzon marinus genome, and the substitution rates of four non-exonic jUCEs were reduced after the teleost-mammal divergence, suggesting that robust conservation had been acquired in the jawed vertebrate lineage. Our results indicate that prototypical UCEs originated before the divergence of jawed and jawless vertebrates and have been frozen as perfect conserved sequences in the jawed vertebrate lineage. In addition, our comparative sequence analyses of UCEs and neighboring regions resulted in a discovery of lineage-specific conserved sequences. They were added progressively to prototypical UCEs, suggesting step-wise acquisition of novel regulatory roles. Our results indicate that conserved non-coding elements (CNEs) consist of blocks with distinct evolutionary history, each having been frozen since different evolutionary era along the vertebrate lineage.     

Evolutionary Dynamics of the Interferon-Induced Transmembrane Gene Family in Vertebrates

Vertebrate interferon-induced transmembrane (IFITM) genes have been demonstrated to have extensive and diverse functions, playing important roles in the evolution of vertebrates. Despite observance of their functionality, the evolutionary dynamics of this gene family are complex and currently unknown. Here, we performed detailed evolutionary analyses to unravel the evolutionary history of the vertebrate IFITM family. A total of 174 IFITM orthologous genes and 112 pseudogenes were identified from 27 vertebrate genome sequences. The vertebrate IFITM family can be divided into immunity-related IFITM (IR-IFITM), IFITM5 and IFITM10 sub-families in phylogeny, implying origins from three different progenitors. In general, vertebrate IFITM genes are located in two loci, one containing the IFITM10 gene, and the other locus containing IFITM5 and various numbers of IR-IFITM genes. Conservation of evolutionary synteny was observed in these IFITM genes. Significant functional divergence was detected among the three IFITM sub-families. No gene duplication or positive selection was found in IFITM5 sub-family, implying the functional conservation of IFITM5 in vertebrate evolution, which is involved in bone formation. No IFITM5 locus was identified in the marmoset genome, suggesting a potential association with the tiny size of this monkey. The IFITM10 sub-family was divided into two groups: aquatic and terrestrial types. Functional divergence was detected between the two groups, and five IFITM10-like genes from frog were dispersed into the two groups. Both gene duplication and positive selection were observed in aquatic vertebrate IFITM10-like genes, indicating that IFITM10 might be associated with the adaptation to aquatic environments. A large number of lineage- and species-specific gene duplications were observed in IR-IFITM sub-family and positive selection was detected in IR-IFITM of primates and rodents. Because primates have experienced a long history of viral infection, such rapid expansion and positive selection suggests that the evolution of primate IR-IFITM genes is associated with broad-spectrum antiviral activity.     

Evolutionary Genomics and Adaptive Evolution of the Hedgehog Gene Family (Shh,Ihh and Dhh) in Vertebrates

The Hedgehog (Hh) gene family codes for a class of secreted proteins composed of two active domains that act as signalling molecules during embryo development, namely for the development of the nervous and skeletal systems and the formation of the testis cord. While only one Hh gene is found typically in invertebrate genomes, most vertebrates species have three (Sonic hedgehog – Shh; Indian hedgehog – Ihh; and Desert hedgehog – Dhh), each with different expression patterns and functions, which likely helped promote the increasing complexity of vertebrates and their successful diversification. In this study, we used comparative genomic and adaptive evolutionary analyses to characterize the evolution of the Hh genes in vertebrates following the two major whole genome duplication (WGD) events. To overcome the lack of Hh-coding sequences on avian publicly available databases, we used an extensive dataset of 45 avian and three non-avian reptilian genomes to show that birds have all three Hh paralogs. We find suggestions that following the WGD events, vertebrate Hh paralogous genes evolved independently within similar linkage groups and under different evolutionary rates, especially within the catalytic domain. The structural regions around the ion-binding site were identified to be under positive selection in the signaling domain. These findings contrast with those observed in invertebrates, where different lineages that experienced gene duplication retained similar selective constraints in the Hh orthologs. Our results provide new insights on the evolutionary history of the Hh gene family, the functional roles of these paralogs in vertebrate species, and on the location of mutational hotspots.     

Molecular evolution of the ankyrin gene family

Ankyrins are membrane adaptor molecules that play important roles in coupling integral membrane proteins to the spectrin-based cytoskeleton network. Human mutations of ankyrin genes lead to severe genetic diseases such as fatal cardiac arrhythmias and hereditary spherocytosis. To elucidate the evolutionary history of ankyrins, we have identified novel ankyrin sequences in insect, fish, frog, chicken, dog, and chimpanzee genomes and explored the phylogenetic relationships of the ankyrin gene family. Our data demonstrate that duplication of ankyrin genes occurred at two different stages. The first duplication resulted from an independent evolution event specific in Arthropoda after its divergence from Chordata. Following the separation from Urochordata, expansion of ankyrins in vertebrates involved ancestral genome duplications. We did not find evidence of coordinated arrangements of gene families of ankyrin-associated membrane proteins on paralogous chromosomes. In addition, evolution of the 24 ANK-repeats strikingly correlated with the exon boundary sites of ankyrin genes, which might have occurred before its duplication in vertebrates. Such correlation is speculated to bring functional diversity and complexity. Moreover, based on the phylogenetic analysis of the ANK-repeat domain, we put forward a novel model for the putative primordial ankyrin that contains the fourth six-ANK-repeat subdomain and the spectrin-binding domain. These findings will provide guides for future studies concerning structure, function, evolutionary origins of ankyrins, and possibly other cytoskeletal proteins.     

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.     

猪GATA-3基因的电子克隆、表达及其分子进化分析

基于猪的表达标签数据库电子克隆了猪的GATA-3基因,并通过RT-PCR验证了猪的GATA-3基因的核苷酸序列。测序结果显示猪的GATA-3基因的核苷酸长度由1,760bp个碱基组成,包括1,335bp的开放阅读框,编码产物为由444个氨基酸残基组成的多肽。通过半定量RT-PCR检测了GATA-3 mRNA在大白猪的各个组织的表达情况。GATA转录因子家族通常有2个Ⅳ型锌指蛋白结构,根据Ⅳ型锌指蛋白结构序列,使用Mega3．1软件构建了分子进化关系树。系统发育分析表明所有的脊椎动物的GATA转录因子都起源于共同的祖先,拓扑结构也表明进化过程中有多种事件发生,包括基因复制和Ⅳ型锌指蛋白结构域重组,根据所得数据有利于进一步了解GATA家族基因趋同和趋异的进化途径。     

Lamprey as an evo-devo model: lessons from comparative embryology and molecular phylogenetics

Lamprey, the living jawless vertebrate, has been regarded as one of the most primitive groups of vertebrates. The evolutionary phylogenetic position of the lamprey promises to provide hints about the origin of the vertebrate genome as well as the origin of the body plan, a part of which may be written in the genome. Since the lamprey split from the gnathostome lineage early in the history of vertebrates, the shared developmental mechanisms in lampreys and gnathostomes can be regarded as possessed by the hypothetical common ancestor of these animals, whereas the gnathostome-specific developmental mechanisms that are absent from lampreys indicate that they are relatively new, added to the developmental program only after the split of gnathostomes. Thus, the sequential establishment of the gnathostome body plan is inherently related to the history of genomic duplication events. In this review, recent molecular developmental and evolutionary molecular research on the living lampreys are summarized and discussed, taking vertebrate comparative morphology and embryology into consideration.     

Selection on coding regions determined Hox7 genes evolution

The important role of Hox genes in determining the regionalization of the body plan of the vertebrates makes them invaluable candidates for evolutionary analyses regarding functional and morphological innovation. Gene duplication and gene loss led to a variable number of Hox genes in different vertebrate lineages. The evolutionary forces determining the conservation or loss of Hox genes are poorly understood. In this study, we show that variable selective pressures acted on Hox7 genes in different evolutionary lineages, with episodes of positive selection occurring after gene duplications. Tests for functional divergence in paralogs detected significant differentiation in a region known to modulate HOX7 protein activity. Our results show that both positive and negative selection on coding regions are influencing Hox7 genes evolution.     

Lactate dehydrogenase (LDH) gene duplication during chordate evolution: the cDNA sequence of the LDH of the tunicate Styela plicata

L-Lactate dehydrogenase (L-LDH, E.C. 1.1.1.27) is encoded by two or three loci in all vertebrates examined, with the exception of lampreys, which have a single LDH locus. Biochemical characterizations of LDH proteins have suggested that a gene duplication early in vertebrate evolution gave rise to Ldh-A and Ldh-B and that an additional locus, Ldh-C arose in a number of lineages more recently. Although some phylogenetic studies of LDH protein sequences have supported this pattern of gene duplication, others have contradicted it. In particular, a number of studies have suggested that Ldh-C represents the earliest divergence among vertebrate LDHs and that it may have diverged from the other loci well before the origin of vertebrates. Such hypotheses make explicit statements about the relationship of vertebrate and invertebrate LDHs, but to date, no closely related invertebrate LDH sequences have been available for comparison. We have attempted to provide further data on the timing of gene duplications leading to multiple vertebrate LDHs by determining the cDNA sequence of the LDH of the tunicate Styela plicata. Phylogenetic analyses of this and other LDH sequences provide strong support for the duplications giving rise to multiple vertebrate LDHs having occurred after vertebrates diverged from tunicates. The timing of these LDH duplications is consistent with data from a number of other gene families suggesting widespread gene duplication near the origin of vertebrates. With respect to the relationships among vertebrate LDHs, our data are not consistent with previous claims that Ldh-C represented the earliest divergence. However, the precise relationships among some of the main lineages of vertebrate LDHs were not resolved in our analyses.