Elucidation of subfamily segregation and intramolecular coevolution of the olfactomedin-like proteins by comprehensive phylogenetic analysis and gene expression pattern assessment

The categorization of genes by structural distinctions relevant to biological characteristics is very important for understanding of gene functions and predicting functional implications of uncharacterized genes. It was absolutely necessary to deploy an effective and efficient strategy to deal with the complexity of the large olfactomedin-like (OLF) gene family sharing sequence similarity but playing diversified roles in many important biological processes, as the simple highest-hit homology analysis gave incomprehensive results and led to inappropriate annotation for some uncharacterized OLF members. In light of evolutionary information that may facilitate the classification of the OLF family and proper association of novel OLF genes with characterized homologs, we performed phylogenetic analysis on all 116 OLF proteins currently available, including two novel members cloned by our group. The OLF family segregated into seven subfamilies and members with similar domain compositions or functional properties all fell into relevant subfamilies. Furthermore, our Northern blot analysis and previous studies revealed that the typical human OLF members in each subfamily exhibited tissue-specific expression patterns, which in turn supported the segregation of the OLF subfamilies with functional divergence. Interestingly, the phylogenetic tree topology for the OLF domains alone was almost identical with that of the full-length tree representing the unique phylogenetic feature of full-length OLF proteins and their particular domain compositions. Moreover, each of the major functional domains of OLF proteins kept the same phylogenetic feature in defining similar topology of the tree. It indicates that the OLF domain and the various domains in flanking non-OLF regions have coevolved and are likely to be functionally interdependent. Expanded by a plausible gene duplication and domain couplings scenario, the OLF family comprises seven evolutionarily and functionally distinct subfamilies, in which each member shares similar structural and functional characteristics including the composition of coevolved and interdependent domains. The phylogenetically classified and preliminarily assessed subfamily framework may greatly facilitate the studying on the OLF proteins. Furthermore, it also demonstrated a feasible and reliable strategy to categorize novel genes and predict the functional implications of uncharacterized proteins based on the comprehensive phylogenetic classification of the subfamilies and their relevance to preliminary functional characteristics.     

Molecular evolution of glutamate receptors: a primitive signaling mechanism that existed before plants and animals diverged.

We performed a genealogical analysis of the ionotropic glutamate receptor (iGluR) gene family, which includes the animal iGluRs and the newly isolated glutamate receptor-like genes (GLR) of plants discovered in Arabidopsis. Distance measures firmly placed the plant GLR genes within the iGluR clade as opposed to other ion channel clades and indicated that iGluRs may be a primitive signaling mechanism that predated the divergence of animals and plants. Moreover, phylogenetic analyses using both parsimony and neighbor joining indicated that the divergence of animal iGluRs and plant GLR genes predated the divergence of iGluR subtypes (NMDA vs. AMPA/KA) in animals. By estimating the congruence of the various glutamate receptor gene regions, we showed that the different functional domains, including the two ligand-binding domains and the transmembrane regions, have coevolved, suggesting that they assembled together before plants and animals diverged. Based on residue conservation and divergence as well as positions of residues with respect to functional domains of iGluR proteins, we attempted to examine structure-function relationships. This analysis defined M3 as the most highly conserved transmembrane domain and identified potential functionally important conserved residues whose function can be examined in future studies.     

Topology and dynamics of the 10 kDa C-terminal domain of DnaK in solution

Hsp70 molecular chaperones contain three distinct structural domains, a 44 kDa N-terminal ATPase domain, a 17 kDa peptide-binding domain, and a 10 kDa C-terminal domain. The ATPase and peptide binding domains are conserved in sequence and are functionally well characterized. The function of the 10 kDa variable C-terminal domain is less well understood. We have characterized the secondary structure and dynamics of the C-terminal domain from the Escherichia coli Hsp70, DnaK, in solution by high-resolution NMR. The domain was shown to be comprised of a rigid structure consisting of four helices and a flexible C-terminal subdomain of approximately 33 amino acids. The mobility of the flexible region is maintained in the context of the full-length protein and does not appear to be modulated by the nucleotide state. The flexibility of this region appears to be a conserved feature of Hsp70 architecture and may have important functional implications. We also developed a method to analyze 15N nuclear spin relaxation data, which allows us to extract amide bond vector directions relative to a unique diffusion axis. The extracted angles and rotational correlation times indicate that the helices form an elongated, bundle-like structure in solution.     

Molecular evolution of the HSP70 multigene family

Eukaryotic genomes encode multiple 70-kDa heat-shock proteins (HSP70s). The Saccharomyces cerevisiae HSP70 family is comprised of eight members. Here we present the nucleotide sequence of the SSA3 and SSB2 genes, completing the nucleotide sequence data for the yeast HSP70 family. We have analyzed these yeast sequences as well as 29 HSP70s from 24 additional eukaryotic and prokaryotic species. Comparison of the sequences demonstrates the extreme conservation of HSP70s; proteins from the most distantly related species share at least 45% identity and more than one-sixth of the amino acids are identical in the aligned region (567 amino acids) among all proteins analyzed. Phylogenetic trees constructed by two independent methods indicate that ancient molecular and cellular events have given rise to at least four monophyletic groups of eukaryotic HSP70 proteins. Each group of evolutionarily similar HSP70s shares a common intracellular localization and is presumed to be comprised of functional homologues; these include heat-shock proteins of the cytoplasm, endoplasmic reticulum, mitochondria, and chloroplasts. HSP70s localized in mitochondria and plastids are most similar to the DnaK HSP70 homologues in purple bacteria and cyanobacteria, respectively, which is consistent with the proposed prokaryotic origin of these organelles. The analyses indicate that the major eukaryotic HSP70 groups arose prior to the divergence of the earliest eukaryotes, roughly 2 billion years ago. In some cases, as exemplified by the SSA genes encoding the cytoplasmic HSP70s of S. cerevisiae, more recent duplication events have given rise to subfamilies within the major groups. The S. cerevisiae SSB proteins comprise a unique subfamily not identified in other species to date. This subfamily appears to have resulted from an ancient gene duplication that occurred at approximately the same time as the origin of the major eukaryotic HSP70 groups. Correspondence to: E.A. Craig     

Phylogenetic analysis and evolution of aromatic amino acid hydroxylase

A study was performed to investigate the phylogenetic relationship among AAAH members and to statistically evaluate sequence conservation and functional divergence. In total, 161 genes were identified from 103 species. Phylogenetic analysis showed that well-conserved subfamilies exist. Exon-intron structure analysis showed that the gene structures of AAAH were highly conserved across some different lineage species, while some species-specific introns were also found. The dynamic distribution of ACT domain suggested one gene fusion event has occurred in eukaryota. Significant functional divergence was found between some subgroups. Analysis of the site-specific profiles revealed critical amino acid residues for functional divergence. This study highlights the molecular evolution of this family and may provide a starting point for further experimental verifications.     

Evolution of thesrc-related protein tyrosine kinases

A phylogenetic analysis ofsrc-related protein tyrosine kinases (PTKs) showed that one group of these genes is quite ancient in the animals, its divergence predating the divergence of the diploblast and triploblast phyla. Three other major groupings of genes were found to predate the divergence of protostome and deuterostome phyla. Most knownsrc-related PTKs of mammals were found to belong to five well-differentiated families: srcA, srcB, abl, csk, and tec. One srcA gene (fyn) has an alternatively spliced seventh exon which shows a different pattern of relationship from the remainder of the gene; this suggests that this exon may have been derived by a recombinational event with another gene, perhaps one related tofgr. The recently published claim that mammalian members of this family expressed in the nervous system evolve more slowly at nonsynonymous nucleotide sites than do those expressed in the immune system was not supported by an analysis of 13 pairs of human and mouse orthologues. Rather, T-cell-specificsrc-related PTKs were found to have higher rates of nonsynonymous substitution than were those having broader expression. This effect was particularly marked in the peptide binding site of the SH2 domain. While the SH2 binding site was highly conserved among paralogous mammalian members of the srcA and srcB subfamilies, no such effect was seen in the comparison of paralogous members of the csk and tec subfamilies. This suggests that, while the peptide binding function of SH2 is conserved within both srcA and srcB subfamilies, paralogous members of the csk and tec subfamilies have diverged functionally with respect to peptide recognition by SH2.     

Feeding truncated heat shock protein 70s protect Artemia franciscana against virulent Vibrio campbellii challenge

The 70 kDa heat shock proteins (Hsp70s) are highly conserved in evolution, leading to striking similarities in structure and composition between eukaryotic Hsp70s and their homologs in prokaryotes. The eukaryotic Hsp70 like the DnaK (Escherichia coli equivalent Hsp70) protein, consist of three functionally distinct domains: an N-terminal 44-kDa ATPase portion, an 18-kDa peptide-binding domain and a C-terminal 10-kDa fragment. Previously, the amino acid sequence of eukaryotic (the brine shrimp Artemia franciscana) Hsp70 and DnaK proteins were shown to share a high degree of homology, particularly in the peptide-binding domain (59.6%, the putative innate immunity-activating portion) compared to the N-terminal ATPase (48.8%) and the C-terminal lid domains (19.4%). Next to this remarkable conservation, these proteins have been shown to generate protective immunity in Artemia against pathogenic Vibrio campbellii. This study, aimed to unravel the Vibrio-protective domain of Hsp70s in vivo, demonstrated that gnotobiotically cultured Artemia fed with recombinant C-terminal fragment (containing the conserved peptide binding domain) of Artemia Hsp70 or DnaK protein were well protected against subsequent Vibrio challenge. In addition, the prophenoloxidase (proPO) system, at both mRNA and protein activity levels, was also markedly induced by these truncated proteins, suggesting epitope(s) responsible for priming the proPO system and presumably other immune-related genes, consequently boosting Artemia survival upon challenge with V. campbellii, might be located within this conserved region of the peptide binding domain.     

美国白蛾热激蛋白70基因的鉴定及功能分析

【目的】探究热激蛋白70(heat shock protein 70, HSP70)基因在美国白蛾Hyphantria cunea抵御高温胁迫过程中的作用,为揭示美国白蛾的扩散机制以及预测潜在分布区提供理论支撑。【方法】利用PCR技术克隆美国白蛾HSP70基因,并进行生物信息学分析;利用qPCR技术检测新蜕皮的美国白蛾4龄第2天幼虫在25, 30, 35和40℃处理下HSP70基因的表达特性;构建美国白蛾HSP70的原核表达载体,并在大肠杆菌Escherichia coli BL21感受态细胞中诱导表达该蛋白,利用Ni²⁺-His柱纯化蛋白,并通过Western blot验证该蛋白;利用体外实验测定原核表达获得的重组蛋白在高温胁迫下的ATP酶活性。【结果】克隆获得美国白蛾2条HSP70基因HcHSP70(GenBank登录号: MT995848)和HcHSC70(GenBank登录号: MT261583),其ORF长度分别为1 917和2 061 bp,分别编码637和687个氨基酸,分子质量分别约为69.66和74.96 kD,等电点分别为5.90和5.96;氨基酸序列结构均含有3个高度保守的区域GIDLGTTYS, IFDLGGGTFDVSIL和VGGSTRIPKVQ,符合热激蛋白70家族的特征;蛋白三维结构均由N端ATPase功能域和C端底物结合功能域所组成。系统进化树显示HcHSP70与鳞翅目HSP70家族其他成员聚为一支,而HcHSC70与鳞翅目HSC70家族的其他成员聚为另一支。qPCR结果表明,新蜕皮的美国白蛾4龄第2天幼虫在热胁迫下HcHSP70的相对表达量显著上调,且在35℃处理2 h下达到峰值,而热胁迫下HcHSC70的表达响应较微弱。成功构建了HcHSP70的原核表达载体并在体外诱导表达。纯化后的重组蛋白HcHSP70具备ATPase活性,且在高温胁迫下酶活性稳定。【结论】本研究克隆获得美国白蛾2条HSP70基因HcHSP70和HcHSC70,明确了两条基因在高温处理下的表达特性;成功对美国白蛾HcHSP70进行原核表达及纯化,并证明重组蛋白HcHSP70在高温下具有稳定的ATPase活性,推测其在美国白蛾应对高温胁迫中发挥着重要作用。     

黄粉虫热休克蛋白70基因的克隆、序列分析与表达（英文）

为给黄粉虫Tenebrio molitor抗逆机理研究提供理论依据, 本研究采用PCR和RACE法从黄粉虫幼虫中克隆出一个热休克蛋白70基因Tmhsp70, 并运用半定量RT-PCR法检测其在黄粉虫不同发育阶段的mRNA表达水平。结果表明： 克隆出的Tmhsp70 序列全长2 282 bp, 具有一个富含A的115 bp 5′ 非翻译区和一个1 935 bp的开放阅读框及一个富含A、 T的232 bp 3′-非翻译区。5′-非翻译区含有7个热休克元件nGAAn, 3′-非翻译区末端有长22 bp的Poly(A)尾。Tmhsp70编码的黄粉虫热休克蛋白（TmHSP70）具有3个典型的HSP70特征基序（IDLGTTYS, IFDLGGGTFDVSIL和IVLVGGSTRIPKIQQ）和1个胞质HSP70末端特征基序（EEVD）, 无信号肽和跨膜区域, 包含2个主要的结构域, 即： N-端42 kDa的高度保守ATPase功能域和C-端18 kDa的保守多肽结合功能域。ATPase功能域的三级结构由2个大球形亚功能域组成, 具有1个核苷酸结合中心; 多肽结合功能域形成1个双层4股β-折叠片样的三明治结构和2个α-螺旋, 内含1个多肽结合通道。此外, 黄粉虫Tmhsp70 mRNA的表达具有热激诱导和发育调控的特征。半定量RT-PCR分析表明, 42℃热激1 h的黄粉虫各发育阶段Tmhsp70 mRNA的表达量上升了1.4~26.9倍。25℃下1日龄黄粉虫蛹中的Tmhsp70 mRNA 表达量要高于其余各发育阶段的累积表达量; 42℃热激1 h 后90日龄幼虫中的Tmhsp70 mRNA 表达量最丰富, 既高于30日龄和60日龄幼虫中的累积表达量, 也高于15日龄和30日龄成虫中的累积表达量。这些结果为进一步研究黄粉虫热休克蛋白的结构、 功能和表达调控及其与抗逆性的关系奠定了基础。     

Comparative analysis of vertebrate <Emphasis Type="Italic">PEPT1</Emphasis> and <Emphasis Type="Italic">PEPT2</Emphasis> genes

The plasma membrane transport proteins belong to SoLute Carrier 15 (SLC15) family and two members of this family have been characterized extensively in higher vertebrates, namely PEPT1 and PEPT2. Despite many efforts have made to define a pharmacophore model for efficient binding and transporting of substrates, there is not a comprehensive study performed to elucidate the evolutionary mechanisms among the SLC15 family members and to statistically evaluate sequence conservation and functional divergence between members. In this study, we compared and contrasted the rates and patterns of molecular evolution of 2 PEPT genes. Phylogenetic tree assembly with all available vertebrate PEPTs suggests that the PEPTs originated by duplications and diverged from a common protein at the base of the eukaryotic tree. Topological structure demonstrates both members share the similar hydrophobic domains (TMDs), which have been constrained by purifying selection. Although both genes show qualitatively similar patterns, their rates of evolution differ significantly due to an increased rate of synonymous substitutions in the structural domains in one copy, suggesting substantial differences in functional constraint on each gene. Site-specific profiles were established by posterior probability analysis revealing significantly divergent regions mainly locate at the hydrophobic region between predicted transmembrane domains 9 and 10 of the proteins. Thus, these results provide the evidence that several amino acid residues with reduced selective constraints are largely responsible for functional divergence between the paralogous PEPTs. These findings may provide a starting point for further experimental verifications.     

Structural and functional differences of SWIRM domain subtypes

SWIRM is a conserved domain found in several chromatin-associated proteins. Based on their sequences, the SWIRM family members can be classified into three subfamilies, which are represented by Swi3, LSD1, and Ada2. Here we report the SWIRM structure of human MYb-like, Swirm and Mpn domain-containing protein-1 (MYSM1). The MYSM1 SWIRM structure forms a compact HTH-related fold comprising five alpha-helices, which best resembles the Swi3 SWIRM structure, among the known SWIRM structures. The MYSM1 and Swi3 SWIRM structures are more similar to the LSD1 structure than the Ada2alpha structure. The SWIRM domains of MYSM1 and LSD1 lacked DNA binding activity, while those of Ada2alpha and the human Swi3 counterpart, SMARCC2, bound DNA. The dissimilarity in the DNA-binding ability of the MYSM1 and SMARCC2 SWIRM domains might be due to a couple of amino acid differences in the last helix. These results indicate that the SWIRM family has indeed diverged into three structural subfamilies (Swi3/MYSM1, LSD1, and Ada2 types), and that the Swi3/MYSM1-type subfamily has further diverged into two functionally distinct groups. We also solved the structure of the SANT domain of MYSM1, and demonstrated that it bound DNA with a similar mode to that of the c-Myb DNA-binding domain.     

Functional characterization of recombinant human HSP70 domains and interdomain interactions

ATPase and peptide-binding activity of recombinant human heat shock proteins HSP70_A1B and HSC70 and two hybrid proteins derived from them was investigated. UV-spectral recorded data were used to characterize conformational rearrangements induced by domain replacement or HSP70-peptide interaction. It was shown that the N-terminal domain dramatically affects the substrate specificity of the C-terminal peptide-binding domain, which puts forward a new hypothesis for HSP70 chaperone machinery. On the other hand, the peptide-binding domain affected the ATPase activity of the recombinant proteins. There was a linear relationship between the ATPase activity and the peptide complex percentage. This connection can be used for quantification of HSP70 complexes with unlabeled peptides.     

Stch encodes the 'ATPase core' of a microsomal stress 70 protein.

The stress70 protein chaperone family plays a central role in the processing of cytosolic and secretory proteins. We have cloned a human cDNA, designated Stch, that is conserved in rat tissues and which encodes a novel microsome-associated member of the stress70 protein chaperone family. Stch mRNA is constitutively expressed in all human cell types and is induced by incubation with the calcium ionophore A23187, but not by exposure to heat shock. Inspection of the predicted amino acid sequence reveals that the STCH product contains a unique hydrophobic leader sequence and shares homology within the amino terminal domains of the stress70 gene family, but has a 50 residue insertion within the ATP-binding domains and truncates the carboxyl terminal peptide-binding region. Immunofluorescent and subcellular analyses show that STCH migrates predominantly as a 60 kDa species and is enriched in a membrane-bound microsome fraction. In contrast to purified BiP and dnaK, however, STCH demonstrates ATPase activity that is independent of peptide stimulation. Stch, therefore, encodes a calcium-inducible, microsome-associated ATPase activity with properties similar to a proteolytically cleaved N-terminal HSC70/BiP fragment. This truncated stress70 molecule may allow increased diversity in cellular responses to protein processing requirements.     

A computational study of elongation factor G (EFG) duplicated genes: diverged nature underlying the innovation on the same structural template

Background

Elongation factor G (EFG) is a core translational protein that catalyzes the elongation and recycling phases of translation. A more complex picture of EFG''s evolution and function than previously accepted is emerging from analyzes of heterogeneous EFG family members. Whereas the gene duplication is postulated to be a prominent factor creating functional novelty, the striking divergence between EFG paralogs can be interpreted in terms of innovation in gene function.

Methodology/Principal Findings

We present a computational study of the EFG protein family to cover the role of gene duplication in the evolution of protein function. Using phylogenetic methods, genome context conservation and insertion/deletion (indel) analysis we demonstrate that the EFG gene copies form four subfamilies: EFG I, spdEFG1, spdEFG2, and EFG II. These ancient gene families differ by their indispensability, degree of divergence and number of indels. We show the distribution of EFG subfamilies and describe evidences for lateral gene transfer and recent duplications. Extended studies of the EFG II subfamily concern its diverged nature. Remarkably, EFG II appears to be a widely distributed and a much-diversified subfamily whose subdivisions correlate with phylum or class borders. The EFG II subfamily specific characteristics are low conservation of the GTPase domain, domains II and III; absence of the trGTPase specific G2 consensus motif “RGITI”; and twelve conserved positions common to the whole subfamily. The EFG II specific functional changes could be related to changes in the properties of nucleotide binding and hydrolysis and strengthened ionic interactions between EFG II and the ribosome, particularly between parts of the decoding site and loop I of domain IV.

Conclusions/Significance

Our work, for the first time, comprehensively identifies and describes EFG subfamilies and improves our understanding of the function and evolution of EFG duplicated genes.     

Functional Evolution of Phosphatidylethanolamine Binding Proteins in Soybean and Arabidopsis

Gene duplication provides resources for novel gene functions. Identification of the amino acids responsible for functional conservation and divergence of duplicated genes will strengthen our understanding of their evolutionary course. Here, we conducted a systemic functional investigation of phosphatidylethanolamine binding proteins (PEBPs) in soybean (Glycine max) and Arabidopsis thaliana. Our results demonstrated that after the ancestral duplication, the lineage of the common ancestor of the FLOWERING LOCUS T (FT) and TERMINAL FLOWER1 (TFL1) subfamilies functionally diverged from the MOTHER OF FT AND TFL1 (MFT) subfamily to activate flowering and repress flowering, respectively. They also underwent further specialization after subsequent duplications. Although the functional divergence increased with duplication age, we observed rapid functional divergence for a few pairs of young duplicates in soybean. Association analysis between amino acids and functional variations identified critical amino acid residues that led to functional differences in PEBP members. Using transgenic analysis, we validated a subset of these differences. We report clear experimental evidence for the functional evolution of the PEBPs in the MFT, FT, and TFL1 subfamilies, which predate the origin of angiosperms. Our results highlight the role of amino acid divergence in driving evolutionary novelty after duplication.