脊椎动物球蛋白家族新成员--胞红蛋白

胞红蛋白(Cygb)是广泛存在于脊椎动物各类组织细胞中的球蛋白家族成员。Cygb为六配位体血红素单体蛋白质,能可逆结合O2、CO等气体配体,参与细胞内某些需氧酶反应。     

Evolution of the Vertebrate Resistin Gene Family

Resistin (encoded by Retn) was previously identified in rodents as a hormone associated with diabetes; however human resistin is instead linked to inflammation. Resistin is a member of a small gene family that includes the resistin-like peptides (encoded by Retnl genes) in mammals. Genomic searches of available genome sequences of diverse vertebrates and phylogenetic analyses were conducted to determine the size and origin of the resistin-like gene family. Genes encoding peptides similar to resistin were found in Mammalia, Sauria, Amphibia, and Actinistia (coelacanth, a lobe-finned fish), but not in Aves or fish from Actinopterygii, Chondrichthyes, or Agnatha. Retnl originated by duplication and transposition from Retn on the early mammalian lineage after divergence of the platypus, but before the placental and marsupial mammal divergence. The resistin-like gene family illustrates an instance where the locus of origin of duplicated genes can be identified, with Retn continuing to reside at this location. Mammalian species typically have a single copy Retn gene, but are much more variable in their numbers of Retnl genes, ranging from 0 to 9. Since Retn is located at the locus of origin, thus likely retained the ancestral expression pattern, largely maintained its copy number, and did not display accelerated evolution, we suggest that it is more likely to have maintained an ancestral function, while Retnl, which transposed to a new location, displays accelerated evolution, and shows greater variability in gene number, including gene loss, likely evolved new, but potentially lineage-specific, functions.     

Phylogenetic and Functional Analysis of the Vertebrate Cytochrome P450 2 Family

Cytochrome P450 (CYP) proteins compose a highly diverse superfamily found in all domains of life. These proteins are enzymes involved in metabolism of endogenous and exogenous compounds. In vertebrates, the CYP2 family is one of the largest, most diverse and plays an important role in mammalian drug metabolism. However, there are more than 20 vertebrate CYP2 subfamilies with uncertain evolution and fairly discrete subfamily composition within vertebrate classes, hindering extrapolation of knowledge across subfamilies. To better understand CYP2 diversity, a phylogenetic analysis of 196 CYP2 protein sequences from 16 species was performed using a maximum likelihood approach and Bayesian inference. The analyses included the CYP2 compliment from human, fugu, zebrafish, stickleback, medaka, cow, and dog genomes. Additional sequences were included from rabbit, marsupial, platypus, chicken, frog, and salmonid species. Three CYP2 sequences from the tunicate Ciona intestinalis were utilized as the outgroup. Results indicate a single ancestral vertebrate CYP2 gene and monophyly of all CYP2 subfamilies. Two subfamilies (CYP2R and CYP2U) pre-date vertebrate diversification, allowing direct comparison across vertebrate classes, while all other subfamilies originated during vertebrate diversification, often within specific vertebrate lineages. Analysis of site-specific evolution indicates that some substrate recognition sites (SRS) previously proposed for CYP genes do not have elevated rates of evolution, suggesting that these regions of the protein are not necessarily important in recognition of CYP2 substrates. Type II functional divergence analysis identified multiple residues in the active site of CYP2F, CYP2A, and CYP2B proteins that have undergone radical biochemical changes and may be functionally important.     

Domain Expansion and Functional Diversification in Vertebrate Reproductive Proteins

The rapid evolution of fertilization proteins has generated remarkable diversity in molecular structure and function. Glycoproteins of vertebrate egg coats contain multiple zona pellucida (ZP)-N domains (1–6 copies) that facilitate multiple reproductive functions, including species-specific sperm recognition. In this report, we integrate phylogenetics and machine learning to investigate how ZP-N domains diversify in structure and function. The most C-terminal ZP-N domain of each paralog is associated with another domain type (ZP-C), which together form a “ZP module.” All modular ZP-N domains are phylogenetically distinct from nonmodular or free ZP-N domains. Machine learning–based classification identifies eight residues that form a stabilizing network in modular ZP-N domains that is absent in free domains. Positive selection is identified in some free ZP-N domains. Our findings support that strong purifying selection has conserved an essential structural core in modular ZP-N domains, with the relaxation of this structural constraint allowing free N-terminal domains to functionally diversify.     

Ascidian and Amphioxus Adh Genes Correlate Functional and Molecular Features of the ADH Family Expansion During Vertebrate Evolution

The alcohol dehydrogenase (ADH) family has evolved into at least eight ADH classes during vertebrate evolution. We have characterized three prevertebrate forms of the parent enzyme of this family, including one from an urochordate (Ciona intestinalis) and two from cephalochordates (Branchiostoma floridae and Branchiostoma lanceolatum). An evolutionary analysis of the family was performed gathering data from protein and gene structures, exon–intron distribution, and functional features through chordate lines. Our data strongly support that the ADH family expansion occurred 500 million years ago, after the cephalochordate/vertebrate split, probably in the gnathostome subphylum line of the vertebrates. Evolutionary rates differ between the ancestral, ADH3 (glutathione-dependent formaldehyde dehydrogenase), and the emerging forms, including the classical alcohol dehydrogenase, ADH1, which has an evolutionary rate 3.6-fold that of the ADH3 form. Phylogenetic analysis and chromosomal mapping of the vertebrate Adh gene cluster suggest that family expansion took place by tandem duplications, probably concurrent with the extensive isoform burst observed before the fish/tetrapode split, rather than through the large-scale genome duplications also postulated in early vertebrate evolution. The absence of multifunctionality in lower chordate ADHs and the structures compared argue in favor of the acquisition of new functions in vertebrate ADH classes. Finally, comparison between B. floridae and B. lanceolatum Adhs provides the first estimate for a cephalochordate speciation, 190 million years ago, probably concomitant with the beginning of the drifting of major land masses from the Pangea. Received: 10 April 2001 / Accepted: 23 May 2001     

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.     

Evolution of prokaryotic homologues of the eukaryotic SEFIR protein domain

SEF/IL17 receptor (SEFIR) domains are mainly found in IL17 receptors (IL17Rs) and their adaptor proteins CIKS (connection to IKK and SAPK/JNK), which exert a host defense role in numbers of infectious diseases and promote inflammatory pathology in autoimmunity. Exploring the evolutionary pathway of SEFIR domains will provide further insight into their functions. Here, we have identified 84 SEFIR domain-containing proteins from more than 1400 prokaryotic genomes. As most SEFIR domain-containing bacterial genomes possess a single SEFIR encoding gene and the SEFIR protein domain forms homodimeric complexes like the Toll/IL1 receptor (TIR) domain, the single bacterial SEFIR proteins may receive binding partners from other organisms. Through comparative and phylogenetic sequence analyses, we show that bacterial SEFIR domain is more similar to that of vertebrate CIKS than IL17R, and it possibly emerges via a lateral gene transfer (LGT) from animals. In addition, our secondary and three-dimensional structural predictions of SEFIR domains reveal that human and pathogenic bacterial SEFIR domains share similar structural and electrostatic features. Our findings provide important clues for further experimental researches on determining the functions of SEFIR proteins in pathogenic prokaryotes.     

八肋游仆虫Rab家族新成员Eo-rab-1N基因的克隆与序列分析

Rab蛋白家族属于小分子GTP结合蛋白家族Ras超家族中最大的亚家族,主要在囊泡运输中起作用。本实验运用PCR、RT-PCR等技术,从八肋游仆虫中克隆到一种新的rab基因。序列分析结果表明：在大核中,该基因全长884bp,除去两端的端粒与非编码区,该基因在大核中由723bp组成。从小核中克隆相应的基因片段,此基因片段序列与大核中序列一致,表明该基因在小核中无内部删除序列的存在。通过RT-PCR,从mRNA获得的该基因的开放读框为663bp,表明该基因在转录过程中有内含子的删除。大核基因序列和cDNA序列比较,发现60bp的内含子序列位于大核基因的153~212bp之间,并符合一类内含子GU-AG剪切规则。在遗传密码使用上,该基因内部含有2个TGA,在游仆虫中编码半胱氨酸。同时首次发现,八肋游仆虫基因使用TAG作为终止密码子。NCBI上序列比对表明该基因翻译的蛋白与其它物种Rab1蛋白的同源性达49%~52%,因此我们将它命名为Eo-rab-1N,GenBank登录号为DQ105562。Eo-rab-1N与其他物种的Rab1蛋白构建进化树,发现该蛋白的进化与物种的进化保持一致,表明该基因在细胞中具有重要功能。     

脊椎动物Sox基因家族的系统发生分析

脊椎动物Sox基因是一类高度保守的基因家族,它们编码一类转录因子,参与多种发育过程的调控,这个家族的共有特征是Sox蛋白具有一个约79个氨基酸,可与DNA特异结合的HMG盒区,该基因家族成员众多复杂,对它们的基因结构,功能及进化关系的研究有助于对该基因家族的全面认识,利用已克隆的全部脊椎动物全长核苷酸／蛋白质数据,对这些数据进行序列比对及进化树的构建,分析Sox基因家族成员的分类及分子进化模式。     

Evolution of the Vertebrate Cytosolic Malate Dehydrogenase Gene Family: Duplication and Divergence in Actinopterygian Fish

Abstract A general correlation between neural expression and negative charge in isozymes suggests charge represents an adaptation to the neural environment. Interestingly, a notable exception exists in teleost fish. Two cytosolic malate dehydrogenase (MDH) isozymes have different spatial expression patterns in certain fishes: one is expressed in all tissues and the second is expressed primarily in the eye and skeletal muscle. While the neural MDH isozyme is negatively charged, the difference in charge between the two isozymes is not as pronounced as that observed in other gene families (e.g., triosephosphate isomerase and lactate dehydrogenase). Most tetrapods express a single cytosolic MDH isozyme, and it has been demonstrated recently that the pair of isozymes found in teleosts results from a gene duplication sometime after the separation of teleosts and tetrapods, although the exact timing of this duplication has not been inferred. Phylogenetic analyses suggest that the duplication of teleost isozymes occurred during the radiation of actinopterygian fish, consistent with the timing of duplication at other loci. Using inferred amino acid sequences, we examine the pattern of change following the duplication and across the rest of the MDH gene tree. Comparison between the MDH gene family and another gene family that shows a larger charge differential among members (triosephosphate isomerase) indicates that the smaller charge difference between MDH isozymes is best explained by greater constraint on amino acid change directly following the duplication, not greater constraint across the entire gene tree. This difference in constraint might result from the wider pattern of expression of the “neural” MDH isozyme.     

The Rho GTPase Family Genes in Bivalvia Genomes: Sequence,Evolution and Expression Analysis

Background

Rho GTPases are important members of the Ras superfamily, which represents the largest signaling protein family in eukaryotes, and function as key molecular switches in converting and amplifying external signals into cellular responses. Although numerous analyses of Rho family genes have been reported, including their functions and evolution, a systematic analysis of this family has not been performed in Mollusca or in Bivalvia, one of the most important classes of Mollusca.

Results

In this study, we systematically identified and characterized a total set (Rho, Rac, Mig, Cdc42, Tc10, Rnd, RhoU, RhoBTB and Miro) of thirty Rho GTPase genes in three bivalve species, including nine in the Yesso scallop Patinopecten yessoensis, nine in the Zhikong scallop Chlamys farreri, and twelve in the Pacific oyster Crassostrea gigas. Phylogenetic analysis and interspecies comparison indicated that bivalves might possess the most complete types of Rho genes in invertebrates. A multiple RNA-seq dataset was used to investigate the expression profiles of bivalve Rho genes, revealing that the examined scallops share more similar Rho expression patterns than the oyster, whereas more Rho mRNAs are expressed in C. farreri and C. gigas than in P. yessoensis. Additionally, Rho, Rac and Cdc42 were found to be duplicated in the oyster but not in the scallops. Among the expanded Rho genes of C. gigas, duplication pairs with high synonymous substitution rates (Ks) displayed greater differences in expression.

Conclusion

A comprehensive analysis of bivalve Rho GTPase family genes was performed in scallop and oyster species, and Rho genes in bivalves exhibit greater conservation than those in any other invertebrate. This is the first study focusing on a genome-wide characterization of Rho GTPase genes in bivalves, and the findings will provide a valuable resource for a better understanding of Rho evolution and Rho GTPase function in Bivalvia.     

Characterization and Functional Analysis of the Calmodulin-Binding Domain of Rac1 GTPase

Rac1, a member of the Rho family of small GTPases, has been shown to promote formation of lamellipodia at the leading edge of motile cells and affect cell migration. We previously demonstrated that calmodulin can bind to a region in the C-terminal of Rac1 and that this interaction is important in the activation of platelet Rac1. Now, we have analyzed amino acid residue(s) in the Rac1-calmodulin binding domain that are essential for the interaction and assessed their functional contribution in Rac1 activation. The results demonstrated that region 151-164 in Rac1 is essential for calmodulin binding. Within the 151-164 region, positively-charged amino acids K153 and R163 were mutated to alanine to study impact on calmodulin binding. Mutant form of Rac1 (K153A) demonstrated significantly reduced binding to calmodulin while the double mutant K153A/R163A demonstrated complete lack of binding to calmodulin. Thrombin or EGF resulted in activation of Rac1 in CHRF-288-11 or HeLa cells respectively and W7 inhibited this activation. Immunoprecipitation studies demonstrated that higher amount of CaM was associated with Rac1 during EGF dependent activation. In cells expressing mutant forms of Rac1 (K153A or K153A/R163A), activation induced by EGF was significantly decreased in comparison to wild type or the R163A forms of Rac1. The lack of Rac1 activation in mutant forms was not due to an inability of GDP-GTP exchange or a change in subcelllular distribution. Moreover, Rac1 activation was decreased in cells where endogenous level of calmodulin was reduced using shRNA knockdown and increased in cells where calmodulin was overexpressed. Docking analysis and modeling demonstrated that K153 in Rac1 interacts with Q41 in calmodulin. These results suggest an important role for calmodulin in the activation of Rac1 and thus, in cytoskeleton reorganization and cell migration.     

Role of Hsp70 ATPase Domain Intrinsic Dynamics and Sequence Evolution in Enabling its Functional Interactions with NEFs

Catalysis of ADP-ATP exchange by nucleotide exchange factors (NEFs) is central to the activity of Hsp70 molecular chaperones. Yet, the mechanism of interaction of this family of chaperones with NEFs is not well understood in the context of the sequence evolution and structural dynamics of Hsp70 ATPase domains. We studied the interactions of Hsp70 ATPase domains with four different NEFs on the basis of the evolutionary trace and co-evolution of the ATPase domain sequence, combined with elastic network modeling of the collective dynamics of the complexes. Our study reveals a subtle balance between the intrinsic (to the ATPase domain) and specific (to interactions with NEFs) mechanisms shared by the four complexes. Two classes of key residues are distinguished in the Hsp70 ATPase domain: (i) highly conserved residues, involved in nucleotide binding, which mediate, via a global hinge-bending, the ATPase domain opening irrespective of NEF binding, and (ii) not-conserved but co-evolved and highly mobile residues, engaged in specific interactions with NEFs (e.g., N57, R258, R262, E283, D285). The observed interplay between these respective intrinsic (pre-existing, structure-encoded) and specific (co-evolved, sequence-dependent) interactions provides us with insights into the allosteric dynamics and functional evolution of the modular Hsp70 ATPase domain.     

The Thalidomide-Binding Domain of Cereblon Defines the CULT Domain Family and Is a New Member of the β-Tent Fold

Despite having caused one of the greatest medical catastrophies of the last century through its teratogenic side-effects, thalidomide continues to be an important agent in the treatment of leprosy and cancer. The protein cereblon, which forms an E3 ubiquitin ligase compex together with damaged DNA-binding protein 1 (DDB1) and cullin 4A, has been recently indentified as a primary target of thalidomide and its C-terminal part as responsible for binding thalidomide within a domain carrying several invariant cysteine and tryptophan residues. This domain, which we name CULT (cereblon domain of unknown activity, binding cellular ligands and thalidomide), is also found in a family of secreted proteins from animals and in a family of bacterial proteins occurring primarily in δ-proteobacteria. Its nearest relatives are yippee, a highly conserved eukaryotic protein of unknown function, and Mis18, a protein involved in the priming of centromeres for recruitment of CENP-A. Searches for distant homologs point to an evolutionary relationship of CULT, yippee, and Mis18 to proteins sharing a common fold, which consists of two four-stranded β-meanders packing at a roughly right angle and coordinating a zinc ion at their apex. A β-hairpin inserted into the first β-meander extends across the bottom of the structure towards the C-terminal edge of the second β-meander, with which it forms a cradle-shaped binding site that is topologically conserved in all members of this fold. We name this the β-tent fold for the striking arrangement of its constituent β-sheets. The fold has internal pseudosymmetry, raising the possibility that it arose by duplication of a subdomain-sized fragment.     

Coevolutionary Landscape of Kinase Family Proteins: Sequence Probabilities and Functional Motifs

The protein kinase catalytic domain is one of the most abundant domains across all branches of life. Although kinases share a common core function of phosphoryl-transfer, they also have wide functional diversity and play varied roles in cell signaling networks, and for this reason are implicated in a number of human diseases. This functional diversity is primarily achieved through sequence variation, and uncovering the sequence-function relationships for the kinase family is a major challenge. In this study we use a statistical inference technique inspired by statistical physics, which builds a coevolutionary “Potts” Hamiltonian model of sequence variation in a protein family. We show how this model has sufficient power to predict the probability of specific subsequences in the highly diverged kinase family, which we verify by comparing the model’s predictions with experimental observations in the Uniprot database. We show that the pairwise (residue-residue) interaction terms of the statistical model are necessary and sufficient to capture higher-than-pairwise mutation patterns of natural kinase sequences. We observe that previously identified functional sets of residues have much stronger correlated interaction scores than are typical.