Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “functional divergence”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2Linked to original sources

Natural history and functional divergence of protein tyrosine kinases.

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

Animals↗

Assessing functional divergence in EF-1alpha and its paralogs in eukaryotes and archaebacteria.

A number of methods have recently been published that use phylogenetic information extracted from large multiple sequence alignments to detect sites that have changed properties in related protein families. In this study we use such methods to assess functional divergence between eukaryotic EF-1alpha (eEF-1alpha), archaebacterial EF-1alpha (aEF-1alpha) and two eukaryote-specific EF-1alpha paralogs-eukaryotic release factor 3 (eRF3) and Hsp70 subfamily B suppressor 1 (HBS1). Overall, the evolutionary modes of aEF-1alpha, HBS1 and eRF3 appear to significantly differ from that of eEF-1alpha. However, functionally divergent (FD) sites detected between aEF-1alpha and eEF-1alpha only weakly overlap with sites implicated as putative EF-1beta or aminoacyl-tRNA (aa-tRNA) binding residues in EF-1alpha, as expected based on the shared ancestral primary translational functions of these two orthologs. In contrast, FD sites detected between eEF-1alpha and its paralogs significantly overlap with the putative EF-1beta and/or aa-tRNA binding sites in EF-1alpha. In eRF3 and HBS1, these sites appear to be released from functional constraints, indicating that they bind neither eEF-1beta nor aa-tRNA. These results are consistent with experimental observations that eRF3 does not bind to aa-tRNA, but do not support the 'EF-1alpha-like' function recently proposed for HBS1. We re-assess the available genetic data for HBS1 in light of our analyses, and propose that this protein may function in stop codon-independent peptide release.

Amino Acid Sequence↗

Functional divergence in protein (family) sequence evolution.

As widely used today to infer 'function', the homology search is based on the neutral theory that sites of greatest functional significance are under the strongest selective constraints as well as lowest evolutionary rates, and vice versa. Therefore, site-specific rate changes (or altered selective constraints) are related to functional divergence during protein (family) evolution. In this paper, we review our recent work about this issue. We show a great deal of functional information can be obtained from the evolutionary perspective, which can in turn be used to facilitate high throughput functional assays. The emergence of evolutionary functional genomics is also indicated. The related software DIVERGE can be obtained from http://xgu1.zool.iastate.edu.

Amino Acid Sequence↗

Adaptive evolution and functional divergence of pepsin gene family.

In vertebrates, a large proportion of genes is organized in gene families. Paralogous gene groups generated by gene duplication are related by homology, high degree of sequence identity and similar structural architecture of their products. Aspartic proteinases form a widely distributed protein superfamily including cathepsins, pepsins, renin and napsin. In the present study, the nucleotide sequences coding for various pepsins in 30 vertebrate species have been used to derive a gene phylogeny. Gene duplication and losses have been inferred from a reconciled tree, reconstructed by combining information from gene tree and species tree. Our findings based on the results of the relative rate ratio test and maximum likelihood analysis suggest that each round of gene duplication is characterized by adaptive evolution, although instances of evolution under positive selection have been found also long after divergence of gene families. The results of functional divergence analysis provided statistical evidence for shifted evolutionary rate after gene duplication.

Adaptation, Physiological↗

Molecular phylogeny, evolution, and functional divergence of the LSD1-like gene family: inference from the rice genome.

The identification of LSD1-like genes in parasite, green algae, moss, pine, and monocot and dicot species allowed us to trace the phylogenetic history of this gene family. Computational analysis showed that the diversification of members of this family could be dated back to the early stage of plant evolution. The evolution of plant LSD1-like genes was possibly shaped by two duplication events. These proteins, which contain three copies of the LSD1 zinc finger (zf-LSD1) domain within their entire polypeptides and play crucial roles in modulating disease defense and cell death, resulted from the second duplication. A gain of zf-LSD1 domain model was reasonable for explaining the origination of three-zf-LSD1 domain-containing proteins. The zf-LSD1 domain phylogeny showed that the middle (M) and C-terminal (C) domains originated from a common ancestor; the N-terminal (N) domain might be more ancient than the former two. The divergence of the N, M, and C domains was well before the monocot-dicot split. Coevolution analysis revealed that four intramolecular domain pairs, including the N domain and the interregion between the M and the C domains (INTER2), the M and C domain, the N- and C-terminus, and the M domain and C-terminus, possibly coevolved during the evolution of three-zf-LSD1 domain-containing proteins. The three zf-LSD1 domains are evolutionary conserved. Thus, the differences at the N- and C-terminus would be crucial for functional specificity of LSD1 genes. Strong functional constraints should work on the zf-LSD1 domains, whereas reduced functional constraint was found in the INTER2 region. Functional divergence analysis showed that three-zf-LSD1 domain-containing proteins were significantly functionally divergent from those proteins containing only one zf-LSD1 domain, a result demonstrating that shifted evolutionary rates between the two clusters were significantly different from each other.

Amino Acid Sequence↗

The hallmark of AGC kinase functional divergence is its C-terminal tail, a cis-acting regulatory module.

The catalytic activities of eukaryotic protein kinases (EPKs) are regulated by movement of the C-helix, movement of the N and C lobes upon ATP binding, and movement of the activation loop upon phosphorylation. Statistical analysis of the selective constraints associated with AGC kinase functional divergence reveals conserved interactions between these regulatory regions and three regions of the C-terminal tail (C-tail): the N-lobe tether (NLT), the active-site tether (AST), and the C-lobe tether (CLT). The NLT serves as a docking site for an upstream kinase PDK1 and, upon activation, positions the C-helix within the ATP binding pocket. The AST directly interacts with the ATP binding pocket, and the CLT interacts with the interlobe linker and the alphaC-beta4 loop, which appears to serve as a hinge for C-helix movement. The C-tail is a hallmark of AGC functional divergence inasmuch as most of the conserved core residues that distinguish AGC kinases from other EPKs are associated with the NLT, AST, or CLT. Moreover, several AGC catalytic core conserved residues that interact with the C-tail strikingly diverge from the canonical residues observed at corresponding positions in nearly all other EPKs, suggesting that the catalytic core may have coevolved with the C-tail in AGC kinases. These observations, along with the fact that the C-tail is needed for catalytic activity suggests that the C-tail is a cis-acting regulatory module that can also serve as a regulatory "handle," to which trans-acting cellular components can bind to modulate activity.

Amino Acid Sequence↗

Functional divergence of six isoforms of antifungal peptide Drosomycin in Drosophila melanogaster.

Drosomycin (Drs) gene encodes a 44-residue inducible antifungal peptide, Drosomycin, in Drosophila melanogaster. Six genes, Drs-lC, Drs-lD, Drs-lE, Drs-lF, Drs-lG and Drs-lI, show homology to the Drs form in a multigene family on the 3rd chromosome of D. melanogaster. It is the first experimental demonstration that the six members in the Drs family act as functional genes. To further delineate the functional divergence of these six members, their cDNA sequences were cloned respectively into the pET-3C vector and expressed in the E. coli. The antifungal activity of the expression products was assayed using the Cerletti's method. The results showed a difference among the six isoforms in antifungal activity against the tested fungal strains: in which Drs was most effective and showed antifungal activity to all seven fungal strains, whereas isoform Drs-lC was effective to six strains, Drs-lD was effective to five strains, Drs-lG was effective to four strains, and Drs-lE and Drs-lF were effective to only three strains. Drs-lI had no activity against any tested fungal strains. By comparing the variable residue sites of these six isoforms to that of Drosomycin in the three-dimensional structure, we suggested that the reduction in the antifungal activity was due to the variable residues that were not in the alpha-helix. In addition, two inserted residues (RV) in Drs-lI may affect the dimensional structure and resulted in a functional change. These results may explain the evolution of the Drosomycin multigene family and its functional divergence.

Amino Acid Sequence↗

Mathematical modeling for functional divergence after gene duplication.

In this paper, I present a statistical framework for modeling the functional divergence after gene duplication. A rate-component model to describe the rate covariation among homologous genes of a gene family is implemented when a phylogenetic tree is known. The Markov chain model is rigorous but may require a huge amount of computational time when the number of sequences is large. On the other hand, the Poisson-based model is mathematically analytical so that computation is very fast even for a large dataset. Moreover, under the posterior framework, we have developed a site-specific profile for predicting important amino acid residues responsible for these functional differences between member genes of a gene family. Our study may have great potential for functional genomics because it is cost-effective, and these predictions can be further tested by biological experimentation.

Biological Evolution↗

Asymmetric Functional Divergence of alx4a and alx4b in Iridophore Differentiation and Cranial Development in Nile Tilapia.

Neural crest cells give rise to the craniofacial skeleton and multiple pigment cell lineages, yet how duplicated developmental regulators partition their ancestral functions after teleost-specific whole-genome duplication remains unclear. Here, we employed CRISPR/Cas9 to generate alx4a and alx4b single and double mutants in Nile tilapia (Oreochromis niloticus). By integrating phenotype, skeleton, transcriptome, quantitative PCR, and AlphaFold-based structural modeling analyses, we revealed their functional divergence. Loss of alx4a caused a regionally restricted reduction in iridophore-derived reflectance and abnormal cranial morphology, whereas alx4b single mutants showed no obvious phenotype under the conditions examined. By contrast, double mutants exhibited an almost complete loss of iridophore-derived structural coloration and substantially more severe cranial defects, accompanied by reduced calcein labeling in the opercular region, consistent with altered cranial mineralization. Skin transcriptomic and quantitative PCR analyses revealed marked downregulation of pnp4a and tfec, which are associated with iridophore differentiation and coloration, whereas no significant expression differences were detected for the iridophore survival-related genes ltk and mpv17. AlphaFold2-assisted HDOCK protein-DNA modeling yielded more favorable docking metrics for Alx4a than for Alx4b with the pnp4a promoter, supporting a potential Alx4a-pnp4a promoter interaction that requires experimental validation. In contrast, no significant genotype-dependent differences were detected in the measured abundance of melanophores, xanthophores, or erythrophores, and no obvious difference in gross dorsal-fin spine formation was observed under the conditions examined. Together, these findings reveal unequal functional contributions of alx4a and alx4b, with alx4a acting as the dominant paralog in iridophore-associated structural coloration and both paralogs contributing unequally to cranial development, and support pnp4a as a candidate downstream gene associated with Alx4a activity.

Animals↗

ago1 and dcr1, two core components of the RNA interference pathway, functionally diverge from rdp1 in regulating cell cycle events in Schizosaccharomyces pombe.

In the fission yeast Schizosaccharomyces pombe, three genes that function in the RNA interference (RNAi) pathway, ago1+, dcr1+, and rdp1+, have recently been shown to be important for timely formation of heterochromatin and accurate chromosome segregation. In the present study, we present evidence that null mutants for ago1+ and dcr1+ but not rdp1+, exhibit abnormal cytokinesis, cell cycle arrest deficiencies, and mating defects. Subsequent analyses showed that ago1+ and dcr1+ are required for regulated hyperphosphorylation of Cdc2 when encountering genotoxic insults. Because rdp1+ is dispensable for this process, the functions of ago1+ and dcr1+ in this pathway are presumably independent of their roles in RNAi-mediated heterochromatin formation and chromosome segregation. This was further supported by the finding that ago1+ is a multicopy suppressor of the S-M checkpoint deficiency and cytokinesis defects associated with loss of Dcr1 function, but not for the chromosome segregation defects of this mutant. Accordingly, we conclude that Dcr1-dependent production of small interfering RNAs is not required for enactment and/or maintenance of certain cell cycle checkpoints and that Ago1 and Dcr1 functionally diverge from Rdp1 to control cell cycle events in fission yeast. Finally, exogenous expression of hGERp95/EIF2C2/hAgo2, a human Ago1 homolog implicated in posttranscriptional gene silencing, compensated for the loss of ago1+ function in S. pombe. This suggests that PPD proteins may also be important for regulation of cell cycle events in higher eukaryotes.

Adenosine Triphosphatases↗

Two functionally divergent p53-responsive elements in the rat bradykinin B2 receptor promoter.

Although p53 is known to have dual functions as a transcriptional activator and repressor, there has not been an example where both p53-activating and -repressing elements reside within one target promoter. Previous work from this laboratory defined two different p53 response elements, termed P1 and P2, located at nucleotide positions -70 and -707, respectively, in the rat bradykinin B2 receptor promoter. In this study, through manipulation of the DNA sequence and context, we demonstrate opposing roles for P1 and P2 as transcriptional activator and repressor, respectively. Deletion of P1 abrogates p53-mediated activation. P1 maintains its role as an activator upon relocation to the P2 site and activates transcription from a heterologous promoter construct. Thus, P1 is a bona fide positive p53-response element. In contrast, deletion of P2 enhances P1-induced activation. P2 represses transcription when substituted for P1 or when relocated midway between P1 and P2. P2-mediated repression is sequence-dependent, because it is reversed to activation when P2 is substituted by the P1 or p53 consensus sequences. Moreover, site-directed mutagenesis that converts P2 to a higher affinity p53-binding site results in transcriptional activation rather than repression. Surprisingly, P2 strongly activates a heterologous promoter. Thus, P2-mediated transcriptional repression is both sequence- and context-dependent. Investigations into the mechanisms of P2-mediated repression indicate that it is trichostatin-insensitive and unaffected by CBP or mutation of the minimal repression C-terminal domain of p53. However, gel shift assays suggest that p53 competes with other transcriptional activators for binding to overlapping binding sequences within the P2 element. In conclusion, this study provides a rare example of a transcription factor having two divergent functional effects that are sequence- and context-dependent. The interplay of P1 and P2 may be important in the regulation of bradykinin B2 receptor gene expression in response to inflammatory stress and during development.

Animals↗

Annotation transfer for genomics: measuring functional divergence in multi-domain proteins.

Annotation transfer is a principal process in genome annotation. It involves "transferring" structural and functional annotation to uncharacterized open reading frames (ORFs) in a newly completed genome from experimentally characterized proteins similar in sequence. To prevent errors in genome annotation, it is important that this process be robust and statistically well-characterized, especially with regard to how it depends on the degree of sequence similarity. Previously, we and others have analyzed annotation transfer in single-domain proteins. Multi-domain proteins, which make up the bulk of the ORFs in eukaryotic genomes, present more complex issues in functional conservation. Here we present a large-scale survey of annotation transfer in these proteins, using scop superfamilies to define domain folds and a thesaurus based on SWISS-PROT keywords to define functional categories. Our survey reveals that multi-domain proteins have significantly less functional conservation than single-domain ones, except when they share the exact same combination of domain folds. In particular, we find that for multi-domain proteins, approximate function can be accurately transferred with only 35% certainty for pairs of proteins sharing one structural superfamily. In contrast, this value is 67% for pairs of single-domain proteins sharing the same structural superfamily. On the other hand, if two multi-domain proteins contain the same combination of two structural superfamilies the probability of their sharing the same function increases to 80% in the case of complete coverage along the full length of both proteins, this value increases further to > 90%. Moreover, we found that only 70 of the current total of 455 structural superfamilies are found in both single and multi-domain proteins and only 14 of these were associated with the same function in both categories of proteins. We also investigated the degree to which function could be transferred between pairs of multi-domain proteins with respect to the degree of sequence similarity between them, finding that functional divergence at a given amount of sequence similarity is always about two-fold greater for pairs of multi-domain proteins (sharing similarity over a single domain) in comparison to pairs of single-domain ones, though the overall shape of the relationship is quite similar. Further information is available at http://partslist.org/func or http://bioinfo.mbb.yale.edu/partslist/func.

Computational Biology↗

Convergent loss of an anciently duplicated, functionally divergent RH2 opsin gene in the fugu and Tetraodon pufferfish lineages.

We describe the complete opsin gene families from the sequenced fugu and Tetraodon pufferfish genomes. We report the convergent loss of function of an anciently duplicated, functionally divergent RH2 or "green-sensitive" opsin gene in both pufferfish lineages, designated RH2-2. In fugu, RH2-2 apparently ceased to function very recently following a transposon-induced deletion that truncated the N-terminal 115 amino acids from the translated protein. Although a lack of frameshift or nonsense mutations in the fugu RH2-2 pseudogene suggests that the gene was lost very recently in this lineage, we were unable to detect any evidence of a selective sweep associated with the fixation of the truncated allele from population data. Interspecific comparison of the remaining fugu RH2-2 coding sequence paradoxically indicates that the gene was under strong purifying selection until the truncation occurred.

Amino Acid Sequence↗

Natural selection during functional divergence to LMP7 and proteasome subunit X (PSMB5) following gene duplication.

The LMP7 and PSMB5 genes were created through an ancient gene duplication event of their ancestral locus. These proteins contain an active site of proteolysis, and LMP7 replaces PSMB5 as a component of the 20S proteasome after stimulation of cells by interferon-gamma. Replacement of PSMB5 by LMP7 changes the profile of the products of 20S proteasome processing, predisposing digested peptides for transport to and display by the immune system. The purpose of this study is to investigate evolutionary forces influencing functional divergence between LMP7 and PSMB5 following duplication. Levels of synonymous and nonsynonymous substitution rates are estimated to infer differences in levels of natural selection. Estimates of substitution rates indicate that natural selection elevated rates of nonsynonymous substitution in LMP7 following gene duplication, whereas PSMB5 experienced an increase in substitution rate that was not likely due to diversifying natural selection following duplication. Following initial divergence, nearly neutral mutations have dominated gene evolution in both lineages. The LMP7 gene locus provides a rare example of a protein with specialized function arising from duplication and divergence of a housekeeping protein by way of natural selection.

Animals↗

Functional divergence prediction from evolutionary analysis: a case study of vertebrate hemoglobin.

It is a central assumption of evolution that gene duplications provide the genetic raw material from which to create proteins with new functions. The increasing availability in multigene family sequences that has resulted from genome projects has inspired the creation of novel in silico approaches to predict details of protein function. The underlying principle of all such approaches is to compare the evolutionary properties of homologous sequence positions in paralogous proteins. It has been proposed that the positions that show switches in substitution rate over time-i.e., "heterotachous sites," are good indicators of functional divergence. Here, we analyzed the alpha and beta paralogous subunits of hemoglobin in search for such signatures. We found as many heterotachous sites in comparisons between groups of paralogous subunits (alpha/beta) as between orthologous ones (alpha/alpha, beta/beta). Thus, the importance of substitution rate shifts as predictors of specialization between protein subfamilies might be reconsidered. Instead, such shifts may reflect a more general process of protein evolution, consistent with the fact that they can be compatible with function conservation. As an alternative, we focused on those residues showing highly constrained states in two sequence groups, but different in each group, and we named them CBD (for "constant but different"). As opposed to heterotachous positions, CBD sites were markedly overrepresented in paralogous (alpha/beta) comparisons, as opposed to orthologous ones (alpha/alpha, beta/beta), identifying them as likely signatures of functional specialization between the two subunits. When superimposed onto the three-dimensional structure of hemoglobin, CBD positions consistently appeared to cluster preferentially on inter-subunit surfaces, two contact areas crucial to function in vertebrate tetrameric hemoglobin. The identification and analysis of CBD sites by complementing structural information with evolutionary data may represent a promising direction for future studies dealing with the functional characterization of a growing number of multigene families identified by complete genome analyses.

Animals↗

Gene duplication, the evolution of novel gene functions, and detecting functional divergence of duplicates in silico.

Duplication of genes increases the amount of genetic material on which evolution can work and has been considered of major importance for the development of biological novelties or to explain important transitions that have occurred during biological evolution. Recently, much research has been devoted to the study of the evolutionary and functional divergence of duplicated genes. Since the majority of genes are part of gene families, there is considerable interest in predicting differences in function between duplicates and assessing the functional redundancy of genes within gene families. In this review, we discuss the strengths and limitations of both older and novel approaches to investigate the evolution of duplicated genes in silico.

Algorithms↗

Functional divergence between histone deacetylases in fission yeast by distinct cellular localization and in vivo specificity.

Histone deacetylases (HDACs) are important for gene regulation and the maintenance of heterochromatin in eukaryotes. Schizosaccharomyces pombe was used as a model system to investigate the functional divergence within this conserved enzyme family. S. pombe has three HDACs encoded by the hda1(+), clr3(+), and clr6(+) genes. Strains mutated in these genes have previously been shown to display strikingly different phenotypes when assayed for viability, chromosome loss, and silencing. Here, conserved differences in the substrate binding pocket identify Clr6 and Hda1 as class I HDACs, while Clr3 belongs in the class II family. Furthermore, these HDACs were shown to have strikingly different subcellular localization patterns. Hda1 was localized to the cytoplasm, while most of Clr3 resided throughout the nucleus. Finally, Clr6 was localized exclusively on the chromosomes in a spotted pattern. Interestingly, Clr3, the only HDAC present in the nucleolus, was required for ribosomal DNA (rDNA) silencing. Clr3 presumably acts directly on heterochromatin, since it colocalized with the centromere, mating-type region, and rDNA as visualized by in situ hybridization. In addition, Clr3 could be cross-linked to mat3 in chromatin immunoprecipitation experiments. Western analysis of bulk histone preparations indicated that Hda1 (class I) had a generally low level of activity in vivo and Clr6 (class I) had a high level of activity and broad in vivo substrate specificity, whereas Clr3 (class II) displayed its main activity on acetylated lysine 14 of histone H3. Thus, the distinct functions of the S. pombe HDACs are likely explained by their distinct cellular localization and their different in vivo specificities.

Acetylation↗