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.     

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.     

The characterisation of six ADAMTS proteases in the basal chordate Ciona intestinalis provides new insights into the vertebrate ADAMTS family

ADAMTS, constituting a recently discovered family of secreted zinc-dependent metalloproteases, have been shown to have critical physiological roles through identification of a number of natural animal and human gene mutations. The identification of six ADAMTS genes in the basal chordate Ciona intestinalis provides new insight into how, when and in what order the vertebrate orthologues have evolved. The phylogenetic assignments, based on sequences conserved across all genes, are supported by conserved domain structures within defined sub-families. The phylogeny and the frequent localisation of ADAMTS genes in paralogous regions of the genome are consistent with the vertebrate lineages having arisen by large scale or genome duplication. The high level of conservation in the protease active site of vertebrate orthologues within some sub-families suggests subfunctionalisation, whereas the greater divergence in others would favour the evolution of novel substrate specificities and these observations are borne-out where substrate-specificity is known. The expansion and sub-specialization of the ADAMTS family is a component of the increased complexity of extracellular matrix that is associated with the evolution of vertebrates.     

Evolutionary History of Cathepsin L (L-like) Family Genes in Vertebrates

Cathepsin L family, an important cysteine protease found in lysosomes, is categorized into cathepsins B, F, H, K, L, S, and W in vertebrates. This categorization is based on their sequence alignment and traditional functional classification, but the evolutionary relationship of family members is unclear. This study determined the evolutionary relationship of cathepsin L family genes in vertebrates through phylogenetic construction. Results showed that cathepsins F, H, S and K, and L and V were chronologically diverged. Tandem-repeat duplication was found to occur in the evolutionary history of cathepsin L family. Cathepsin L in zebrafish, cathepsins S and K in xenopus, and cathepsin L in mice and rats underwent evident tandem-repeat events. Positive selection was detected in cathepsin L-like members in mice and rats, and amino acid sites under positive selection pressure were calculated. Most of these sites appeared at the connection of secondary structures, suggesting that the sites may slightly change spatial structure. Severe positive selection was also observed in cathepsin V (L2) of primates, indicating that this enzyme had some special functions. Our work provided a brief evolutionary history of cathepsin L family and differentiated cathepsins S and K from cathepsin L based on vertebrate appearance. Positive selection was the specific cause of differentiation of cathepsin L family genes, confirming that gene function variation after expansion events was related to interactions with the environment and adaptability.     

Evolutionary analysis of vertebrate Notch genes

We have conducted an evolutionary analysis of Notch genes of the vertebrates Danio rerio and Mus musculus to examine the expansion and diversification of the Notch family during vertebrate evolution. The existence of multiple Notch genes in vertebrate genomes suggests that the increase in Notch signaling pathways may be necessary for the additional complexity observed in the vertebrate body plan. However, orthology relationships within the vertebrate Notch family indicate that biological functions are not fixed within orthologous groups. Phylogenetic reconstruction of the vertebrate Notch family suggests that the zebrafish notch1a and 1b genes resulted from a duplication occurring around the time of the teleost/mammalian divergence. There is also evidence that the mouse Notch4 gene is the result of a rapid divergence from a Notch3-like gene. Investigation of the ankyrin repeat region sequences showed there to be little evidence for gene conversion events between repeat units. However, relationships between repeats 2-5 suggest that these repeats are the result of a tandem duplication of a dual repeat unit. Selective pressure on maintenance of ankyrin repeat sequences indicated by relationships between the repeats suggests that specific repeats are responsible for particular biological activities, a finding consistent with mutational studies of the Caenorhabditis elegans gene glp-1. Sequence similarities between the ankyrin repeats and the region immediately C-terminal of the repeats further suggests that this region may be involved in the modulation of ankyrin repeat function.     

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.     

Maximum likelihood methods for detecting adaptive evolution after gene duplication

The rapid accumulation of genomic sequences in public databases will finally allow large scale studies of gene family evolution, including evaluation of the role of positive Darwinian selection following a duplication event. This will be possible because recent statistical methods of comparing synonymous and nonsynonymous substitution rates permit reliable detection of positive selection at individual amino acid sites and along evolutionary lineages. Here, we summarize maximum-likelihood based methods, and present a framework for their application to analysis of gene families. Using these methods, we investigated the role of positive Darwinian selection in the ECP-EDN gene family of primates and the Troponin C gene family of vertebrates. We also comment on the limitations of these methods and discuss directions for further improvements.     

Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages

The nonsynonymous (amino acid-altering) to synonymous (silent) substitution rate ratio (omega = d(N)/d(S)) provides a measure of natural selection at the protein level, with omega = 1, >1, and <1, indicating neutral evolution, purifying selection, and positive selection, respectively. Previous studies that used this measure to detect positive selection have often taken an approach of pairwise comparison, estimating substitution rates by averaging over all sites in the protein. As most amino acids in a functional protein are under structural and functional constraints and adaptive evolution probably affects only a few sites at a few time points, this approach of averaging rates over sites and over time has little power. Previously, we developed codon-based substitution models that allow the omega ratio to vary either among lineages or among sites. In this paper we extend previous models to allow the omega ratio to vary both among sites and among lineages and implement the new models in the likelihood framework. These models may be useful for identifying positive selection along prespecified lineages that affects only a few sites in the protein. We apply those branch-site models as well as previous branch- and site-specific models to three data sets: the lysozyme genes from primates, the tumor suppressor BRCA1 genes from primates, and the phytochrome (PHY) gene family in angiosperms. Positive selection is detected in the lysozyme and BRCA genes by both the new and the old models. However, only the new models detected positive selection acting on lineages after gene duplication in the PHY gene family. Additional tests on several data sets suggest that the new models may be useful in detecting positive selection after gene duplication in gene family evolution.     

Molecular evolution of growth hormone gene family in old world monkeys and hominoids

Growth hormone is a classic molecule in the study of the molecular clock hypothesis as it exhibits a relatively constant rate of evolution in most mammalian orders except primates and artiodactyls, where dramatically enhanced rate of evolution (25–50-fold) has been reported. The rapid evolution of primate growth hormone occurred after the divergence of tarsiers and simians, but before the separation of old world monkeys (OWM) from new world monkeys (NWM). Interestingly, this event of rapid sequence evolution coincided with multiple duplications of the growth hormone gene, suggesting gene duplication as a possible cause of the accelerated sequence evolution. Here we determined 21 different GH-like sequences from four species of OWM and hominoids. Combining with published sequences from OWM and hominoids, our analysis demonstrates that multiple gene duplications and several gene conversion events both occurred in the evolutionary history of this gene family in OWM/hominoids. The episode of recent duplications of CSH-like genes in gibbon is accompanied with rapid sequence evolution likely resulting from relaxation of purifying selection. GHN genes in both hominoids and OWM are under strong purifying selection. In contrast, CSH genes in both lineages are probably not. GHV genes in OWM and hominoids evolved at different evolutionary rates and underwent different selective constraints. Our results disclosed the complex history of the primate growth hormone gene family and raised intriguing questions on the consequences of these evolutionary events.     

Marsupials and monotremes possess a novel family of MHC class I genes that is lost from the eutherian lineage

Background

Major histocompatibility complex (MHC) class I genes are found in the genomes of all jawed vertebrates. The evolution of this gene family is closely tied to the evolution of the vertebrate genome. Family members are frequently found in four paralogous regions, which were formed in two rounds of genome duplication in the early vertebrates, but in some species class Is have been subject to additional duplication or translocation, creating additional clusters. The gene family is traditionally grouped into two subtypes: classical MHC class I genes that are usually MHC-linked, highly polymorphic, expressed in a broad range of tissues and present endogenously-derived peptides to cytotoxic T-cells; and non-classical MHC class I genes generally have lower polymorphism, may have tissue-specific expression and have evolved to perform immune-related or non-immune functions. As immune genes can evolve rapidly and are subject to different selection pressure, we hypothesised that there may be divergent, as yet unannotated or uncharacterised class I genes.

Results

Application of a novel method of sensitive genome searching of available vertebrate genome sequences revealed a new, extensive sub-family of divergent MHC class I genes, denoted as UT, which has not previously been characterized. These class I genes are found in both American and Australian marsupials, and in monotremes, at an evolutionary chromosomal breakpoint, but are not present in non-mammalian genomes and have been lost from the eutherian lineage. We show that UT family members are expressed in the thymus of the gray short-tailed opossum and in other immune tissues of several Australian marsupials. Structural homology modelling shows that the proteins encoded by this family are predicted to have an open, though short, antigen-binding groove.

Conclusions

We have identified a novel sub-family of putatively non-classical MHC class I genes that are specific to marsupials and monotremes. This family was present in the ancestral mammal and is found in extant marsupials and monotremes, but has been lost from the eutherian lineage. The function of this family is as yet unknown, however, their predicted structure may be consistent with presentation of antigens to T-cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1745-4) contains supplementary material, which is available to authorized users.     

Rapid expansion of killer cell immunoglobulin-like receptor genes in primates and their coevolution with MHC Class I genes

The gene family of killer cell immunoglobulin-like receptors (KIRs) in primates provides the first line of defense against virus infection and tumor transformation. Interacting with MHC class I molecules, KIRs can regulate the cytotoxic activity of natural killer (NK) cells and distinguish the tumor and virus infected cells from normal body cells. Phylogenetic analysis and comparison of domain structures identified three major groups of KIR genes (group I, II, and III genes). These groups of KIR genes, generated by a series of gene duplications, have acquired different MHC-binding specificity. Inference of ancestral KIR sequences suggested that the functional divergence of group I genes from group II genes occurred by positive selection at the MHC-binding sites after duplication. Our evolutionary study has shown that group I genes diverged from group II genes about 17 million years ago (Mya) apparently after separation of hominoids from Old World (OW) monkeys. Around the same time, gene duplication generating the class I MHC-C locus appears to have occurred. These findings suggest that KIR and MHC class I genes have coevolved as an interacting system. The KIR gene family has experienced a rapid expansion in primate species. The rate of expansion of this gene family seems to be one of the highest among all hominoid gene families. The KIR gene family is also subject to birth-and-death evolution.     

Cytochrome P450 1 genes in early deuterostomes (tunicates and sea urchins) and vertebrates (chicken and frog): origin and diversification of the CYP1 gene family

Cytochrome P450 family 1 (CYP1) proteins are important in a large number of toxicological processes. CYP1A and CYP1B genes are well known in mammals, but the evolutionary history of the CYP1 family as a whole is obscure; that history may provide insight into endogenous functions of CYP1 enzymes. Here, we identify CYP1-like genes in early deuterostomes (tunicates and echinoderms), and several new CYP1 genes in vertebrates (chicken, Gallus gallus and frog, Xenopus tropicalis). Profile hidden Markov models (HMMs) generated from vertebrate CYP1A and CYP1B protein sequences were used to identify 5 potential CYP1 homologs in the tunicate Ciona intestinalis genome. The C. intestinalis genes were cloned and sequenced, confirming the predicted sequences. Orthologs of 4 of these genes were found in the Ciona savignyi genome. Bayesian phylogenetic analyses group the tunicate genes in the CYP1 family, provisionally in 2 new subfamilies, CYP1E and CYP1F, which fall in the CYP1A and CYP1B/1C clades. Bayesian and maximum likelihood analyses predict functional divergence between the tunicate and vertebrate CYP1s, and regions within CYP substrate recognition sites were found to differ significantly in position-specific substitution rates between tunicates and vertebrates. Subsequently, 10 CYP1-like genes were found in the echinoderm Strongylocentrotus purpuratus (sea urchin) genome. Several of the tunicate and echinoderm CYP1-like genes are expressed during development. Canonical xenobiotic response elements are present in the upstream genomic sequences of most tunicate and sea urchin CYP1s, and both groups are predicted to possess an aryl hydrocarbon receptor (AHR), suggesting possible regulatory linkage of AHR and these CYPs. The CYP1 family has undergone multiple rounds of gene duplication followed by functional divergence, with at least one gene lost in mammals. This study provides new insight into the origin and evolution of CYP1 genes.     

Gene Conversion and Functional Divergence in the β-Globin Gene Family

Different models of gene family evolution have been proposed to explain the mechanism whereby gene copies created by gene duplications are maintained and diverge in function. Ohta proposed a model which predicts a burst of nonsynonymous substitutions following gene duplication and the preservation of duplicates through positive selection. An alternative model, the duplication–degeneration–complementation (DDC) model, does not explicitly require the action of positive Darwinian selection for the maintenance of duplicated gene copies, although purifying selection is assumed to continue to act on both copies. A potential outcome of the DDC model is heterogeneity in purifying selection among the gene copies, due to partitioning of subfunctions which complement each other. By using the d_N/d_S () rate ratio to measure selection pressure, we can distinguish between these two very different evolutionary scenarios. In this study we investigated these scenarios in the -globin family of genes, a textbook example of evolution by gene duplication. We assembled a comprehensive dataset of 72 vertebrate -globin sequences. The estimated phylogeny suggested multiple gene duplication and gene conversion events. By using different programs to detect recombination, we confirmed several cases of gene conversion and detected two new cases. We tested evolutionary scenarios derived from Ohtas model and the DDC model by examining selective pressures along lineages in a phylogeny of -globin genes in eutherian mammals. We did not find significant evidence for an increase in the ratio following major duplication events in this family. However, one exception to this pattern was the duplication of -globin in simian primates, after which a few sites were identified to be under positive selection. Overall, our results suggest that following gene duplications, paralogous copies of -globin genes evolved under a nonepisodic process of functional divergence.[Reviewing Editor: Martin Kreitman]     

Adaptive evolution of the uncoupling protein 1 gene contributed to the acquisition of novel nonshivering thermogenesis in ancestral eutherian mammals

Homeotherms possess various physiological mechanisms to maintain their body temperature, thus allowing them to adapt to various environments. Under cold conditions, most eutherian mammals upregulate heat production in brown adipose tissue (BAT), and uncoupling protein (UCP) 1 is an essential factor in BAT thermogenesis. The evolutionary origin of UCP1 was believed to have been a specific event occurring in eutherian lineages. Recently, however, the UCP1 ortholog was found in fishes, which uncovers a more ancient origin of this gene than previously believed. Here we investigate the evolutionary process of UCP1 by comparative genomic approach. We found that UCP1 evolved rapidly by positive Darwinian selection in the common ancestor of eutherians, although this gene arose in the ancestral vertebrate, since the orthologous genes were shared among most of the vertebrate species. Adaptive evolution occurred after the divergence between eutherians and marsupials, which is consistent with the fact that BAT has been found only in eutherians. Our findings indicate that positive Darwinian selection acted on UCP1 contributed to the acquisition of an efficient mechanism for body temperature regulation in primitive eutherians. Phylogenetic reconstruction of UCP1 with two paralogs (UCP2 and UCP3) among vertebrate species revealed that the gene duplication events which produced these three genes occurred in the common ancestor of vertebrates much earlier than the emergence of eutherians. Thus, our data demonstrate that novel gene function can evolve without de novo gene duplication event.     

Multidimensional scaling reveals the main evolutionary pathways of class A G-protein-coupled receptors

Class A G-protein-coupled receptors (GPCRs) constitute the largest family of transmembrane receptors in the human genome. Understanding the mechanisms which drove the evolution of such a large family would help understand the specificity of each GPCR sub-family with applications to drug design. To gain evolutionary information on class A GPCRs, we explored their sequence space by metric multidimensional scaling analysis (MDS). Three-dimensional mapping of human sequences shows a non-uniform distribution of GPCRs, organized in clusters that lay along four privileged directions. To interpret these directions, we projected supplementary sequences from different species onto the human space used as a reference. With this technique, we can easily monitor the evolutionary drift of several GPCR sub-families from cnidarians to humans. Results support a model of radiative evolution of class A GPCRs from a central node formed by peptide receptors. The privileged directions obtained from the MDS analysis are interpretable in terms of three main evolutionary pathways related to specific sequence determinants. The first pathway was initiated by a deletion in transmembrane helix 2 (TM2) and led to three sub-families by divergent evolution. The second pathway corresponds to the differentiation of the amine receptors. The third pathway corresponds to parallel evolution of several sub-families in relation with a covarion process involving proline residues in TM2 and TM5. As exemplified with GPCRs, the MDS projection technique is an important tool to compare orthologous sequence sets and to help decipher the mutational events that drove the evolution of protein families.     

Heterogeneous rate of protein evolution in serotonin genes

Serotonin (5-hydroxytryptamine) is a neurotransmitter crucial for cardiovascular, gastrointestinal, and brain function. It is also involved in several aspects of behavior and associated with a variety of personality disorders in humans. Its dual role as a crucial element in vital physiological functions (strictly evolutionary conserved) and in traits that differ substantially across species makes the evolution of serotonin function particularly interesting. We studied the evolution of serotonin function through the identification of the selective forces shaping the evolution of genes in its functional pathway in primates and rodents. Serotonin genes are highly conserved and show no signals of positive selection, suggesting functional constraint as the main force driving their evolution. They show, nevertheless, considerable differences in constraint between primates and rodents, with some genes showing dramatic differences between the 2 groups. These genes most likely represent cases of functional divergence between primates and rodents and point out to the relevance of using closely related species in gene-based evolutionary studies to avoid the effect of unrecognized functional differences between distant species. Within each group (rodents or primates), genes also show heterogeneity in evolution. Genes from the same gene family (with structure and function alike) tend to evolve at a similar rate, but this is not always the case. A few serotonin genes show substantial differences in constraint with the rest of members of their family, suggesting the presence of important and unrecognized functional differences among the genes, which may be involved in species-specific evolution.     

Natural mutations in IFITM3 modulate post‐translational regulation and toggle antiviral specificity

The interferon‐induced transmembrane (IFITM) proteins protect host cells from diverse virus infections. IFITM proteins also incorporate into HIV‐1 virions and inhibit virus fusion and cell‐to‐cell spread, with IFITM3 showing the greatest potency. Here, we report that amino‐terminal mutants of IFITM3 preventing ubiquitination and endocytosis are more abundantly incorporated into virions and exhibit enhanced inhibition of HIV‐1 fusion. An analysis of primate genomes revealed that IFITM3 is the most ancient antiviral family member of the IFITM locus and has undergone a repeated duplication in independent host lineages. Some IFITM3 genes in nonhuman primates, including those that arose following gene duplication, carry amino‐terminal mutations that modify protein localization and function. This suggests that “runaway” IFITM3 variants could be selected for altered antiviral activity. Furthermore, we show that adaptations in IFITM3 result in a trade‐off in antiviral specificity, as variants exhibiting enhanced activity against HIV‐1 poorly restrict influenza A virus. Overall, we provide the first experimental evidence that diversification of IFITM3 genes may boost the antiviral coverage of host cells and provide selective functional advantages.     

α/β Barrel evolution and the modular assembly of enzymes: Emerging trends in the flavin oxidase/dehydrogenase family

α/β barrels have an ill-defined origin. Evidence exists which favours their divergent evolution from a common ancestral barrel and convergent evolution to a stable fold. However, recent sequence and structural information for the flavin oxidase/dehydrogenase family of barrel enzymes indicate that sub-families of α/β barrels have evolved divergently. The modular fusion of barrel domains with core structures from other gene families has also contributed to the evolution of related but catalytically distinct enzyme molecules within each sub-family of the flavin oxidases/dehydrogenases. An analysis of the structures and sequences of the flavin oxidases/dehydrogenases has now enabled studies focusing on the evolutionary origins and modular assembly of this important family of proteins to be initiated.     

Functional characterization of the Arabidopsis ubiquitin-specific protease gene family reveals specific role and redundancy of individual members in development

Ubiquitin-specific proteases (UBPs) are a highly conserved family of proteins in eukaryotes, and play critical roles in protein de-ubiquitination. Here we report a systematic genetic and expression profiling analysis of the UBP gene family in the Arabidopsis thaliana genome. Mutation analysis of 25 of the 27 member genes representing 13 of the 14 sub-families of the UBP gene family revealed that single-gene mutants of three genes in two sub-families exhibit visible phenotypes. Two of these three genes belonging to the UBP15 sub-family were selected for further characterization. The ubp15 mutants display narrower, serrated and flat rosette leaves, partially due to a defect in cell proliferation, as well as other phenotypes such as early flowering, weak apical dominance and reduced fertility, while the line over-expressing UBP15 shows opposite phenotypes. We demonstrated that UPB15 has UBP activity in vitro , and that this biochemical activity is essential for its in vivo function. A genetic interaction analysis among members of this sub-family revealed that UBP15 and UBP16, but not UBP17, have functional redundancy. Our data thus suggest that distinct UBPs, even within a closely related sub-family, can function in different developmental pathways. Although there are clearly functional redundancies among related sub-family members, those redundancies cannot be inferred simply based on the amino acid identity of the family members.