Using the COG Database to Improve Gene Recognition in Complete Genomes

A complete understanding of the biology of an organism necessarily starts with knowledge of its genetic makeup. Proteins encoded in a genome must be identified and characterized, and the presence or absence of specific sets of proteins must be noted in order to determine the possible biochemical pathways or functional systems utilized by that organism. The COG database presents a set of tools suited to these purposes, including the ability to select protein families (COGs) that contain proteins from a specified set of species. The selection is based upon a phylogenetic pattern, which is a shorthand representation of the presence or absence of a particular species in a COG. Here we present the use of phylogenetic patterns as a means to perform targeted searches for undetected protein-coding genes in complete genomes.     

Putative DHHC-Cysteine-Rich Domain S-Acyltransferase in Plants

Protein S-acyltransferases (PATs) containing Asp-His-His-Cys within a Cys-rich domain (DHHC-CRD) are polytopic transmembrane proteins that are found in eukaryotic cells and mediate the S-acylation of target proteins. S-acylation is an important secondary and reversible modification that regulates the membrane association, trafficking and function of target proteins. However, little is known about the characteristics of PATs in plants. Here, we identified 804 PATs from 31 species with complete genomes. The analysis of the phylogenetic relationships suggested that all of the PATs fell into 8 groups. In addition, we analysed the phylogeny, genomic organization, chromosome localisation and expression pattern of PATs in Arabidopsis, Oryza sative, Zea mays and Glycine max. The microarray data revealed that PATs genes were expressed in different tissues and during different life stages. The preferential expression of the ZmPATs in specific tissues and the response of Zea mays to treatments with phytohormones and abiotic stress demonstrated that the PATs play roles in plant growth and development as well as in stress responses. Our data provide a useful reference for the identification and functional analysis of the members of this protein family.     

An insight into domain combinations

Domains are the building blocks of all globular proteins, and are units of compact three-dimensional structure as well as evolutionary units. There is a limited repertoire of domain families, so that these domain families are duplicated and combined in different ways to form the set of proteins in a genome. Proteins are gene products. The processes that produce new genes are duplication and recombination as well as gene fusion and fission. We attempt to gain an overview of these processes by studying the structural domains in the proteins of seven genomes from the three kingdoms of life: Eubacteria, Archaea and Eukaryota. We use here the domain and superfamily definitions in Structural Classification of Proteins Database (SCOP) in order to map pairs of adjacent domains in genome sequences in terms of their superfamily combinations. We find 624 out of the 764 superfamilies in SCOP in these genomes, and the 624 families occur in 585 pairwise combinations. Most families are observed in combination with one or two other families, while a few families are very versatile in their combinatorial behaviour. This type of pattern can be described by a scale-free network. Finally, we study domain repeats and we compare the set of the domain combinations in the genomes to those in PDB, and discuss the implications for structural genomics.     

Comprehensive sequence analysis of horseshoe crab cuticular proteins and their involvement in transglutaminase-dependent cross-linking

Arthropod cuticles play an important role as the first barrier against invading pathogens. We extensively determined the sequences of horseshoe crab cuticular proteins. Proteins extracted from a part of the ventral side of the cuticle were purified by chitin-affinity chromatography, and separated by two-dimensional SDS/PAGE. Proteins appearing on the gel were designated high molecular mass chitin-binding proteins, and these proteins were then grouped into classes based on their approximate isoelectric points and predominant amino acid compositions. Members of groups designated basic G, basic Y, and acidic S groups contained a so-called Rebers and Riddiford consensus found in arthropod cuticular proteins. Proteins designated acidic DE25 and DE29 each contained a Cys-rich domain with sequences similar to those of insect peritrophic matrix proteins and chitinases. In contrast, basic QH4 and QH10 contained no consensus sequences found in known chitin-binding proteins. Alternatively, a low molecular mass chitin-binding fraction was prepared by size exclusion chromatography, and 15 low molecular mass chitin-binding proteins, named P1 through P15, were isolated. With the exception of P9 and P15, all were found to be identical to known antimicrobial peptides. P9 consisted of a Kunitz-type chymotrypsin inhibitor sequence, and P15 contained a Cys-rich motif found in insulin-like growth factor-binding proteins. Interestingly, we observed transglutaminase-dependent polymerization of nearly all high molecular mass chitin-binding proteins, a finding suggests that transglutaminase-dependent cross-linking plays an important role in host defense in the arthropod cuticle, analogous to that observed in the epidermal cornified cell envelope in mammals.     

Evolution of patchily distributed proteins shared between eukaryotes and prokaryotes: Dictyostelium as a case study

Protein families are often patchily distributed in the tree of life; they are present in distantly related organisms, but absent in more closely related lineages. This could either be the result of lateral gene transfer between ancestors of organisms that encode them, or losses in the lineages that lack them. Here a novel approach is developed to study the evolution of patchily distributed proteins shared between prokaryotes and eukaryotes. Proteins encoded in the genome of cellular slime mold Dictyostelium discoideum and a restricted number of other lineages, including at least one prokaryote, were identified. Analyses of the phylogenetic distribution of 49 such patchily distributed protein families showed conflicts with organismal phylogenies; 25 are shared with the distantly related amoeboflagellate Naegleria (Excavata), whereas only two are present in the more closely related Entamoeba. Most protein families show unexpected topologies in phylogenetic analyses; eukaryotes are polyphyletic in 85% of the trees. These observations suggest that gene transfers have been an important mechanism for the distribution of patchily distributed proteins across all domains of life. Further studies of this exchangeable gene fraction are needed for a better understanding of the origin and evolution of eukaryotic genes and the diversification process of eukaryotes.     

The soybean vegetative storage proteins VSP alpha and VSP beta are acid phosphatases active on polyphosphates.

The soybean vegetative storage protein genes (vspA, and vspB) are regulated in a complex manner developmentally and in response to external stimuli such as wounding and water deficit. The proteins accumulate to almost one-half the amount of soluble leaf protein when soybean plants are continually depodded and have been identified as storage proteins because of their abundance and pattern of expression in plant tissues. We have shown that purified VSP homodimers (VSP alpha and VSP beta) and heterodimers (VSP alpha/beta) possess acid phosphatase activity (alpha = 0.3-0.4 units/mg; beta = 2-4 units/mg; alpha/beta = 7-10 units/mg). Specific activities were determined by monitoring o-carboxyphenyl phosphate (0.7 mM) cleavage at pH 5.5 (VSP alpha) or pH 5.0 (VSP alpha/beta and VSP beta) in 0.15 M sodium acetate buffer at 25 degrees C. These enzymes are active over a broad pH range, maintaining greater than 40% of maximal activity from pH 4.0 to 6.5 and having maximal activity at pH 5.0-5.5. They are inactivated by sodium fluoride, sodium molybdate, and heating at 70 degrees C for 10 min. These phosphatases can liberate Pi from several different substrates, including napthyl acid phosphate, carboxyphenyl phosphate, sugar-phosphates, glyceraldehyde 3-phosphate, dihydroxyacetone phosphate, phosphoenolpyruvate, ATP, ADP, PPi, and short chain polyphosphates. VSP alpha/beta cleaved phosphoenolpyruvate, ATP, ADP, PPi, and polyphosphates most efficiently. Apparent Km and Vmax values at 25 degrees C and pH 5.0 were 42 microM and 2.0 mumol/min/mg, 150 microM and 4.2 mumol/min/mg, and 420 microM and 4.1 mumol/min/mg, for tetrapolyphosphate, pyrophosphate, and phosphoenolpyruvate, respectively.     

Reconstitution of the Very Short Patch Repair Pathway from Escherichia coli

The Escherichia coli very short patch (VSP) repair pathway corrects thymidine-guanine mismatches that result from spontaneous hydrolytic deamination damage of 5-methyl cytosine. The VSP repair pathway requires the Vsr endonuclease, DNA polymerase I, a DNA ligase, MutS, and MutL to function at peak efficiency. The biochemical roles of most of these proteins in the VSP repair pathway have been studied extensively. However, these proteins have not been studied together in the context of VSP repair in an in vitro system. Using purified components of the VSP repair system in a reconstitution reaction, we have begun to develop an understanding of the role played by each of these proteins in the VSP repair pathway and have gained insights into their interactions. In this report we demonstrate an in vitro reconstitution of the VSP repair pathway using a plasmid DNA substrate. Surprisingly, the repair track length can be modulated by the concentration of DNA ligase. We propose roles for MutL and MutS in coordination of this repair pathway.     

The Crystal Structure of Arabidopsis VSP1 Reveals the Plant Class C-Like Phosphatase Structure of the DDDD Superfamily of Phosphohydrolases

Arabidopsis thaliana vegetative storage proteins, VSP1 and VSP2, are acid phosphatases and belong to the haloacid dehalogenase (HAD) superfamily. In addition to their potential nutrient storage function, they were thought to be involved in plant defense and flower development. To gain insights into the architecture of the protein and obtain clues about its function, we have tested their substrate specificity and solved the structure of VSP1. The acid phosphatase activities of these two enzymes require divalent metal such as magnesium ion. Conversely, the activity of these two enzymes is inhibited by vanadate and molybdate, but is resistant to inorganic phosphate. Both VSP1 and VSP2 did not exhibit remarkable activities to any physiological substrates tested. In the current study, we presented the crystal structure of recombinant VSP1 at 1.8 Å resolution via the selenomethionine single-wavelength anomalous diffraction (SAD). Specifically, an α-helical cap domain on the top of the α/β core domain is found to be involved in dimerization. In addition, despite of the low sequence similarity between VSP1 and other HAD enzymes, the core domain of VSP1 containing conserved active site and catalytic machinery displays a classic haloacid dehalogenase fold. Furthermore, we found that VSP1 is distinguished from bacterial class C acid phosphatase P4 by several structural features. To our knowledge, this is the first study to reveal the crystal structure of plant vegetative storage proteins.     

Divergent polo box domains underpin the unique kinetoplastid kinetochore

Kinetochores are macromolecular machines that drive eukaryotic chromosome segregation by interacting with centromeric DNA and spindle microtubules. While most eukaryotes possess conventional kinetochore proteins, evolutionarily distant kinetoplastid species have unconventional kinetochore proteins, composed of at least 19 proteins (KKT1–19). Polo-like kinase (PLK) is not a structural kinetochore component in either system. Here, we report the identification of an additional kinetochore protein, KKT20, in Trypanosoma brucei. KKT20 has sequence similarity with KKT2 and KKT3 in the Cys-rich region, and all three proteins have weak but significant similarity to the polo box domain (PBD) of PLK. These divergent PBDs of KKT2 and KKT20 are sufficient for kinetochore localization in vivo. We propose that the ancestral PLK acquired a Cys-rich region and then underwent gene duplication events to give rise to three structural kinetochore proteins in kinetoplastids.     

Structure and function of cytokinin oxidase/dehydrogenase genes of maize,rice,<Emphasis Type="Italic"> Arabidopsis</Emphasis> and other species

Cytokinin oxidases/dehydrogenases (CKX) catalyze the irreversible degradation of the cytokinins isopentenyladenine, zeatin, and their ribosides in a single enzymatic step by oxidative side chain cleavage. To date the sequences of 17 fully annotated CKX genes are known, including two prokaryotic genes. The CKX gene families of Arabidopsis thaliana and rice comprise seven and at least ten members, respectively. The main features of CKX genes and proteins are summarized in this review. Individual proteins differ in their catalytic properties, their subcellular localization and their expression domains. The evolutionary development of cytokinin-catabolizing gene families and the individual properties of their members indicate an important role for the fine-tuned control of catabolism to assure proper regulation of cytokinin functions. The use of CKX genes as a tool in studies of cytokinin biology and biotechnological applications is discussed.     

Epigenetic silencing of Plasmodium falciparum genes linked to erythrocyte invasion

The process of erythrocyte invasion by merozoites of Plasmodium falciparum involves multiple steps, including the formation of a moving junction between parasite and host cell, and it is characterised by the redundancy of many of the receptor-ligand interactions involved. Several parasite proteins that interact with erythrocyte receptors or participate in other steps of invasion are encoded by small subtelomerically located gene families of four to seven members. We report here that members of the eba, rhoph1/clag, acbp, and pfRh multigene families exist in either an active or a silenced state. In the case of two members of the rhoph1/clag family, clag3.1 and clag3.2, expression was mutually exclusive. Silencing was clonally transmitted and occurred in the absence of detectable DNA alterations, suggesting that it is epigenetic. This was demonstrated for eba-140. Our data demonstrate that variant or mutually exclusive expression and epigenetic silencing in Plasmodium are not unique to genes such as var, which encode proteins that are exported to the surface of the erythrocyte, but also occur for genes involved in host cell invasion. Clonal variant expression of invasion-related ligands increases the flexibility of the parasite to adapt to its human host.     

Accurate quantification of functional analogy among close homologs

Correctly evaluating functional similarities among homologous proteins is necessary for accurate transfer of experimental knowledge from one organism to another, and is of particular importance for the development of animal models of human disease. While the fact that sequence similarity implies functional similarity is a fundamental paradigm of molecular biology, sequence comparison does not directly assess the extent to which two proteins participate in the same biological processes, and has limited utility for analyzing families with several parologous members. Nevertheless, we show that it is possible to provide a cross-organism functional similarity measure in an unbiased way through the exclusive use of high-throughput gene-expression data. Our methodology is based on probabilistic cross-species mapping of functionally analogous proteins based on Bayesian integrative analysis of gene expression compendia. We demonstrate that even among closely related genes, our method is able to predict functionally analogous homolog pairs better than relying on sequence comparison alone. We also demonstrate that the landscape of functional similarity is often complex and that definitive "functional orthologs" do not always exist. Even in these cases, our method and the online interface we provide are designed to allow detailed exploration of sources of inferred functional similarity that can be evaluated by the user.