Phospho-RNA sequencing with circAID-p-seq.

Most RNA footprinting approaches that require ribonuclease cleavage generate RNA fragments bearing a phosphate or cyclic phosphate group at their 3' end. Unfortunately, current library preparation protocols rely only on a 3' hydroxyl group for adaptor ligation or poly-A tailing. Here, we developed circAID-p-seq, a PCR-free library preparation for selective 3' phospho-RNA sequencing. As a proof of concept, we applied circAID-p-seq to ribosome profiling, which is based on sequencing of RNA fragments protected by ribosomes after endonuclease digestion. CircAID-p-seq, combined with the dedicated computational pipeline circAidMe, facilitates accurate, fast and highly efficient sequencing of phospho-RNA fragments from eukaryotic cells and tissues. We used circAID-p-seq to portray ribosome occupancy in transcripts, providing a versatile and PCR-free strategy to possibly unravel any endogenous 3'-phospho RNA molecules.

Quantitative frame analysis and the annotation of GC-rich (and other) prokaryotic genomes. An application to Anaeromyxobacter dehalogenans.

MOTIVATION: Graphical representations of contrasts in GC usage among codon frame positions (frame analysis) provide evidence of genes missing from the annotations of prokaryotic genomes of high GC content but the qualitative approach of visual frame analysis prevents its applicability on a genomic scale. RESULTS: We developed two quantitative methods for the identification and statistical characterization in sequence regions of three-base periodicity (hits) associated with open reading frame structures. The methods were implemented in the N-Profile Analysis Computational Tool (NPACT), which highlights in graphical representations inconsistencies between newly identified ORFs and pre-existing annotations of coding-regions. We applied the NPACT procedures to two recently annotated strains of the deltaproteobacterium Anaeromyxobacter dehalogenans, identifying in both genomes numerous conserved ORFs not included in the published annotation of coding regions. AVAILABILITY AND IMPLEMENTATION: NPACT is available as a web-based service and for download at http://genome.ufl.edu/npact. CONTACT: lucianob@ufl.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Archaeal ubiquitin-like SAMP3 is isopeptide-linked to proteins via a UbaA-dependent mechanism.

SAMP1 and SAMP2 are ubiquitin-like proteins that function as protein modifiers and are required for the production of sulfur-containing biomolecules in the archaeon Haloferax volcanii. Here we report a novel small archaeal modifier protein (named SAMP3) with a ß-grasp fold and C-terminal diglycine motif characteristic of ubiquitin that is functional in protein conjugation in Hfx. volcanii. SAMP3 conjugates were dependent on the ubiquitin-activating E1 enzyme homolog of archaea (UbaA) for synthesis and were cleaved by the JAMM/MPN+ domain metalloprotease HvJAMM1. Twenty-three proteins (28 lysine residues) were found to be isopeptide-linked to the C-terminal carboxylate of SAMP3, and 331 proteins were reproducibly found associated with SAMP3 in a UbaA-dependent manner based on tandem mass spectrometry (MS/MS) analysis. The molybdopterin (MPT) synthase large subunit homolog MoaE, found samp3ylated at conserved active site lysine residues in MS/MS analysis, was also shown to be covalently bound to SAMP3 by immunoprecipitation and tandem affinity purifications. HvJAMM1 was demonstrated to catalyze the cleavage of SAMP3 from MoaE, suggesting a mechanism of controlling MPT synthase activity. The levels of samp3ylated proteins and samp3 transcripts were found to be increased by the addition of dimethyl sulfoxide to aerobically growing cells. Thus, we propose a model in which samp3ylation is covalent and reversible and controls the activity of enzymes such as MPT synthase. Sampylation of MPT synthase may govern the levels of molybdenum cofactor available and thus facilitate the scavenging of oxygen prior to the transition to respiration with molybdenum-cofactor-containing terminal reductases that use alternative electron acceptors such as dimethyl sulfoxide. Overall, our study of SAMP3 provides new insight into the diversity of functional ubiquitin-like protein modifiers and the network of ubiquitin-like protein targets in Archaea.

Ancient Origin of Chaperonin Gene Paralogs Involved in Ciliopathies.

The Bardet-Biedl Syndrome (BBS) is a human developmental disorder that has been associated with fourteen BBS genes affecting the development of cilia. Three BBS genes are distant relatives of chaperonin proteins, a family of chaperones well known for the protein-folding role of their double-ringed complexes. Chaperonin-like BBS genes were originally thought to be vertebrate-specific, but related genes from different metazoan species have been identified as chaperonin-like BBS genes based on sequence similarity. Our phylogenetic analyses confirmed the classification of these genes in the chaperonin-like BBS gene family, and set the origin of the gene family earlier than the time of separation of Bilateria, Cnidaria, and Placozoa. By extensive searches of chaperonin-like genes in complete genomes representing several eukaryotic lineages, we discovered the presence of chaperonin-like BBS genes also in the genomes of Phytophthora and Pythium, belonging to the group of Oomycetes. This finding suggests that the chaperonin-like BBS gene family had already evolved before the origin of Metazoa, as early in eukaryote evolution as before separation of the lineages of Unikonts and Chromalveolates. The analysis of coding sequences indicated that chaperonin-like BBS proteins have evolved in all lineages under constraining selection. Furthermore, analysis of the predicted structural features suggested that, despite their high rate of divergence, chaperonin-like BBS proteins mostly conserve a typical chaperonin-like three-dimensional structure, but question their ability to assemble and function as chaperonin-like double-ringed complexes.

Chaperonin genes on the rise: new divergent classes and intense duplication in human and other vertebrate genomes.

BACKGROUND: Chaperonin proteins are well known for the critical role they play in protein folding and in disease. However, the recent identification of three diverged chaperonin paralogs associated with the human Bardet-Biedl and McKusick-Kaufman Syndromes (BBS and MKKS, respectively) indicates that the eukaryotic chaperonin-gene family is larger and more differentiated than previously thought. The availability of complete genome sequences makes possible a definitive characterization of the complete set of chaperonin sequences in human and other species. RESULTS: We identified fifty-four chaperonin-like sequences in the human genome and similar numbers in the genomes of the model organisms mouse and rat. In mammal genomes we identified, besides the well-known CCT chaperonin genes and the three genes associated with the MKKS and BBS pathological conditions, a newly-defined class of chaperonin genes named CCT8L, represented in human by the two sequences CCT8L1 and CCT8L2. Comparative analyses from several vertebrate genomes established the monophyletic origin of chaperonin-like MKKS and BBS genes from the CCT8 lineage. The CCT8L gene originated from a later duplication also in the CCT8 lineage at the onset of mammal evolution and duplicated in primate genomes. The functionality of CCT8L genes in different species was confirmed by evolutionary analyses and in human by expression data. Detailed sequence analysis and structural predictions of MKKS, BBS and CCT8L proteins strongly suggested that they conserve a typical chaperonin-like core structure but that they are unlikely to form a CCT-like oligomeric complex. The characterization of many newly-discovered chaperonin pseudogenes uncovered the intense duplication activity of eukaryotic chaperonin genes. CONCLUSIONS: In vertebrates, chaperonin genes, driven by intense duplication processes, have diversified into multiple classes and functionalities that extend beyond their well-known protein-folding role as part of the typical oligomeric chaperonin complex, emphasizing previous observations on the involvement of individual CCT monomers in microtubule elongation. The functional characterization of newly identified chaperonin genes will be a challenge for future experimental analyses.

A comprehensive classification and evolutionary analysis of plant homeobox genes.

The full complement of homeobox transcription factor sequences, including genes and pseudogenes, was determined from the analysis of 10 complete genomes from flowering plants, moss, Selaginella, unicellular green algae, and red algae. Our exhaustive genome-wide searches resulted in the discovery in each class of a greater number of homeobox genes than previously reported. All homeobox genes can be unambiguously classified by sequence evolutionary analysis into 14 distinct classes also characterized by conserved intron-exon structure and by unique codomain architectures. We identified many new genes belonging to previously defined classes (HD-ZIP I to IV, BEL, KNOX, PLINC, WOX). Other newly identified genes allowed us to characterize PHD, DDT, NDX, and LD genes as members of four new evolutionary classes and to define two additional classes, which we named SAWADEE and PINTOX. Our comprehensive analysis allowed us to identify several newly characterized conserved motifs, including novel zinc finger motifs in SAWADEE and DDT. Members of the BEL and KNOX classes were found in Chlorobionta (green plants) and in Rhodophyta. We found representatives of the DDT, WOX, and PINTOX classes only in green plants, including unicellular green algae, moss, and vascular plants. All 14 homeobox gene classes were represented in flowering plants, Selaginella, and moss, suggesting that they had already differentiated in the last common ancestor of moss and vascular plants.

hsp70 genes in the human genome: Conservation and differentiation patterns predict a wide array of overlapping and specialized functions.

BACKGROUND: Hsp70 chaperones are required for key cellular processes and response to environmental changes and survival but they have not been fully characterized yet. The human hsp70-gene family has an unknown number of members (eleven counted over ten years ago); some have been described but the information is incomplete and inconsistent. A coherent body of knowledge encompassing all family components that would facilitate their study individually and as a group is lacking. Nowadays, the study of chaperone genes benefits from the availability of genome sequences and a new protocol, chaperonomics, which we applied to elucidate the human hsp70 family. RESULTS: We identified 47 hsp70 sequences, 17 genes and 30 pseudogenes. The genes distributed into seven evolutionarily distinct groups with distinguishable subgroups according to phylogenetic and other data, such as exon-intron and protein features. The N-terminal ATP-binding domain (ABD) was conserved at least partially in the majority of the proteins but the C-terminal substrate-binding domain (SBD) was not. Nine proteins were typical Hsp70s (65-80 kDa) with ABD and SBD, two were lighter lacking partly or totally the SBD, and six were heavier (>80 kDa) with divergent C-terminal domains. We also analyzed exon-intron features, transcriptional variants and protein structure and isoforms, and modality and patterns of expression in various tissues and developmental stages. Evolutionary analyses, including human hsp70 genes and pseudogenes, and other eukaryotic hsp70 genes, showed that six human genes encoding cytosolic Hsp70s and 27 pseudogenes originated from retro-transposition of HSPA8, a gene highly expressed in most tissues and developmental stages. CONCLUSION: The human hsp70-gene family is characterized by a remarkable evolutionary diversity that mainly resulted from multiple duplications and retrotranspositions of a highly expressed gene, HSPA8. Human Hsp70 proteins are clustered into seven evolutionary Groups, with divergent C-terminal domains likely defining their distinctive functions. These functions may also be further defined by the observed differences in the N-terminal domain.

Chaperonomics, a new tool to study ageing and associated diseases.

The participation of molecular chaperones in the process of senescence and in the mechanisms of age-related diseases is currently under investigation in many laboratories. However, accurate, complete information about the number and diversity of chaperone genes in any given genome is scarce. Consequently, the results of efforts aimed at elucidating the role of chaperones in ageing and disease are often confusing and contradictory. To remedy this situation, we have developed chaperonomics, including means to identify and characterize chaperone genes and their families applicable to humans and model organisms. The problem is difficult because in eukaryotic organisms chaperones have evolved into complex multi-gene families. For instance, the occurrence of multiple paralogs in a single genome makes it difficult to interpret results if consideration is not given to the fact that similar but distinct chaperone genes can be differentially expressed in separate cellular compartments, tissues, and developmental stages. The availability of complete genome sequences allows implementation of chaperonomics with the purpose of understanding the composition of chaperone families in all cell compartments, their evolutionary and functional relations and, ultimately, their role in pathogenesis. Here, we present a series of concatenated, complementary procedures for identifying, characterizing, and classifying chaperone genes in genomes and for elucidating evolutionary relations and structural features useful in predicting functional properties. We illustrate the procedures with applications to the complex family of hsp70 genes and show that the kind of data obtained can provide a solid basis for future research.

Distinguishing features of delta-proteobacterial genomes.

We analyzed several features of five currently available delta-proteobacterial genomes, including two aerobic bacteria exhibiting predatory behavior and three anaerobic sulfate-reducing bacteria. The delta genomes are distinguished from other bacteria by several properties: (i) The delta genomes contain two "giant" S1 ribosomal protein genes in contrast to all other bacterial types, which encode a single or no S1; (ii) in most delta-proteobacterial genomes the major ribosomal protein (RP) gene cluster is near the replication terminus whereas most bacterial genomes place the major RP cluster near the origin of replication; (iii) the delta genomes possess the rare combination of discriminating asparaginyl and glutaminyl tRNA synthetase (AARS) together with the amido-transferase complex (Gat CAB) genes that modify Asp-tRNA(Asn) into Asn-tRNA(Asn) and Glu-tRNA(Gln) into Gln-tRNA(Gln); (iv) the TonB receptors and ferric siderophore receptors that facilitate uptake and removal of complex metals are common among delta genomes; (v) the anaerobic delta genomes encode multiple copies of the anaerobic detoxification protein rubrerythrin that can neutralize hydrogen peroxide; and (vi) sigma(54) activators play a more important role in the delta genomes than in other bacteria. delta genomes have a plethora of enhancer binding proteins that respond to environmental and intracellular cues, often as part of two-component systems; (vii) delta genomes encode multiple copies of metallo-beta-lactamase enzymes; (viii) a host of secretion proteins emphasizing SecA, SecB, and SecY may be especially useful in the predatory activities of Myxococcus xanthus; (ix) delta proteobacteria drive many multiprotein machines in their periplasms and outer membrane, including chaperone-feeding machines, jets for slime secretion, and type IV pili. Bdellovibrio replicates in the periplasm of prey cells. The sulfate-reducing delta proteobacteria metabolize hydrogen and generate a proton gradient by electron transport. The predicted highly expressed genes from delta genomes reflect their different ecologies, metabolic strategies, and adaptations.

Evolution of a protein-folding machine: genomic and evolutionary analyses reveal three lineages of the archaeal hsp70(dnaK) gene.

The stress chaperone protein Hsp70 (DnaK) (abbreviated DnaK) and its co-chaperones Hsp40(DnaJ) (or DnaJ) and GrpE are universal in bacteria and eukaryotes but occur only in some archaea clustered in the order 5'-grpE-dnaK-dnaJ-3' in a locus termed Locus I. Three structural varieties of Locus I, termed Types I, II, and III, were identified, respectively, in Methanosarcinales, in Thermoplasmatales and Methanothermobacter thermoautotrophicus, and in Halobacteriales. These Locus I types corresponded to three groups identified by phylogenetic trees of archaeal DnaK proteins including the same archaeal subdivisions. These archaeal DnaK groups were not significantly interrelated, clustering instead with DnaKs from three bacterial lineages, Methanosarcinales with Firmicutes, Thermoplasmatales and M. thermoautotrophicus with Thermotoga, and Halobacteriales with Actinobacteria, suggesting that the three archaeal types of Locus I were acquired by independent events of lateral gene transfer. These associations, however, lacked strong bootstrap support and were sensitive to dataset choice and tree-reconstruction method. Structural features of dnaK loci in bacteria revealed that Methanosarcinales and Firmicutes shared a similar structure, also common to most other bacterial groups. Structural differences were observed instead in Thermotoga compared to Thermoplasmatales and M. thermoautotrophicus, and in Actinobacteria compared to Halobacteriales. It was also found that the association between the DnaK sequences from Halobacteriales and Actinobacteria likely reflects common biases in their amino acid compositions. Although the loci structural features and the DnaK trees suggested the possibility of lateral gene transfer between Firmicutes and Methanosarcinales, the similarity between the archaeal and the ancestral bacterial loci favors the more parsimonious hypothesis that all archaeal sequences originated from a unique prokaryotic ancestor.

Protein length in eukaryotic and prokaryotic proteomes.

We analyzed length differences of eukaryotic, bacterial and archaeal proteins in relation to function, conservation and environmental factors. Comparing Eukaryotes and Prokaryotes, we found that the greater length of eukaryotic proteins is pervasive over all functional categories and involves the vast majority of protein families. The magnitude of these differences suggests that the evolution of eukaryotic proteins was influenced by processes of fusion of single-function proteins into extended multi-functional and multi-domain proteins. Comparing Bacteria and Archaea, we determined that the small but significant length difference observed between their proteins results from a combination of three factors: (i) bacterial proteomes include a greater proportion than archaeal proteomes of longer proteins involved in metabolism or cellular processes, (ii) within most functional classes, protein families unique to Bacteria are generally longer than protein families unique to Archaea and (iii) within the same protein family, homologs from Bacteria tend to be longer than the corresponding homologs from Archaea. These differences are interpreted with respect to evolutionary trends and prevailing environmental conditions within the two prokaryotic groups.

Genomic and proteomic comparisons between bacterial and archaeal genomes and related comparisons with the yeast and fly genomes.

Bacterial, archaeal, yeast, and fly genomes are compared with respect to predicted highly expressed (PHX) genes and several genomic properties. There is a striking difference in the status of PHX ribosomal protein (RP) genes where the archaeal genome generally encodes more RP genes and fewer PHX RPs compared with bacterial genomes. The increase in RPs in archaea and eukaryotes compared with that in bacteria may reflect a more complex set of interactions in archaea and eukaryotes in regulating translation, e.g., differences in structure requiring scaffolding of longer rRNA molecules, expanded interactions with the chaperone machinery, and, in eukaryotic interactions with endoplasmic reticulum components. The yeast genome is similar to fast-growing bacteria in PHX genes but also features several cytoskeletal genes, including actin and tropomyosin, and several signal transduction regulatory proteins from the 14.3.3 family. The most PHX genes of Drosophila encode cytoskeletal and exoskeletal proteins. We found that the preference of a microorganism for an anaerobic metabolism correlates with the number of PHX enzymes of the glycolysis pathway that well exceeds the number of PHX enzymes acting in the tricarboxylic acid cycle. Conversely, if the number of PHX enzymes of the tricarboxylic acid cycle well exceeds the PHX enzymes of glycolysis, an aerobic metabolism is preferred. Where the numbers are approximately commensurate, a facultative growth behavior prevails.

Predicted highly expressed genes in archaeal genomes.

Based primarily on 16S rRNA sequence comparisons, life has been broadly divided into the three domains of Bacteria, Archaea, and Eukarya. Archaea is further classified into Crenarchaea and Euryarchaea. Archaea generally thrive in extreme environments as assessed by temperature, pH, and salinity. For many prokaryotic organisms, ribosomal proteins (RP), transcription/translation factors, and chaperone genes tend to be highly expressed. A gene is predicted highly expressed (PHX) if its codon usage is rather similar to the average codon usage of at least one of the RP, transcription/translation factors, and chaperone gene classes and deviates strongly from the average gene of the genome. The thermosome (Ths) chaperonin family represents the most salient PHX genes among Archaea. The chaperones Trigger factor and HSP70 have overlapping functions in the folding process, but both of these proteins are lacking in most archaea where they may be substituted by the chaperone prefoldin. Other distinctive PHX proteins of Archaea, absent from Bacteria, include the proliferating cell nuclear antigen PCNA, a replication auxiliary factor responsible for tethering the catalytic unit of DNA polymerase to DNA during high-speed replication, and the acidic RP P0, which helps to initiate mRNA translation at the ribosome. Other PHX genes feature Cell division control protein 48 (Cdc48), whereas the bacterial septation proteins FtsZ and minD are lacking in Crenarchaea. RadA is a major DNA repair and recombination protein of Archaea. Archaeal genomes feature a strong Shine-Dalgarno ribosome-binding motif more pronounced in Euryarchaea compared with Crenarchaea.

Predicting coding potential from genome sequence: application to betaherpesviruses infecting rats and mice.

Prediction of protein-coding regions and other features of primary DNA sequence have greatly contributed to experimental biology. Significant challenges remain in genome annotation methods, including the identification of small or overlapping genes and the assessment of mRNA splicing or unconventional translation signals in expression. We have employed a combined analysis of compositional biases and conservation together with frame-specific G+C representation to reevaluate and annotate the genome sequences of mouse and rat cytomegaloviruses. Our analysis predicts that there are at least 34 protein-coding regions in these genomes that were not apparent in earlier annotation efforts. These include 17 single-exon genes, three new exons of previously identified genes, a newly identified four-exon gene for a lectin-like protein (in rat cytomegalovirus), and 10 probable frameshift extensions of previously annotated genes. This expanded set of candidate genes provides an additional basis for investigation in cytomegalovirus biology and pathogenesis.

The composition, structure and stability of a group II chaperonin are temperature regulated in a hyperthermophilic archaeon.

The hyperthermoacidophilic archaeon Sulfolobus shibatae contains group II chaperonins, known as rosettasomes, which are two nine-membered rings composed of three different 60 kDa subunits (TF55 alpha, beta and gamma). We sequenced the gene for the gamma subunit and studied the temperature-dependent changes in alpha, beta and gamma expression, their association into rosettasomes and their phylogenetic relationships. Alpha and beta gene expression was increased by heat shock (30 min, 86 degrees C) and decreased by cold shock (30 min, 60 degrees C). Gamma expression was undetectable at heat shock temperatures and low at normal temperatures (75-79 degrees C), but induced by cold shock. Polyacrylamide gel electrophoresis indicated that in vitro alpha and beta subunits form homo-oligomeric rosettasomes, and mixtures of alpha, beta and gamma form hetero-oligomeric rosettasomes. Transmission electron microscopy revealed that beta homo-oligomeric rosettasomes and all hetero-oligomeric rosettasomes associate into filaments. In vivo rosettasomes were hetero-oligomeric with an average subunit ratio of 1alpha:1beta:0.1gamma in cultures grown at 75 degrees C, a ratio of 1alpha:3beta:1gamma in cultures grown at 60 degrees C and a ratio of 2alpha:3beta:0gamma after 86 degrees C heat shock. Using differential scanning calorimetry, we determined denaturation temperatures (Tm) for alpha, beta and gamma subunits of 95.7 degrees C, 96.7 degrees C and 80.5 degrees C, respectively, and observed that rosettasomes containing gamma were relatively less stable than those with alpha and/or beta only. We propose that, in vivo, the rosettasome structure is determined by the relative abundance of subunits and not by a fixed geometry. Furthermore, phylogenetic analyses indicate that archaeal chaperonin subunits underwent multiple duplication events within species (paralogy). The independent evolution of these paralogues raises the possibility that chaperonins have functionally diversified between species.

Heterogeneity of genome and proteome content in bacteria, archaea, and eukaryotes.

Our analysis compares bacteria, archaea, and eukaryota with respect to a wide assortment of genome and proteome properties. These properties include ribosomal protein gene distributions, chaperone protein contrasts, major variation of transcription/translation factors, gene encoding pathways of energy metabolism, and predicted protein expression levels. Significant differences within and between the three domains of life include protein lengths, information processing procedures, many metabolic and lipid biosynthesis pathways, cellular controls, and regulatory proteins. Differences among genomes are influenced by lifestyle, habitat, physiology, energy sources, and other factors.

Amino acid runs in eukaryotic proteomes and disease associations.

We present a comparative proteome analysis of the five complete eukaryotic genomes (human, Drosophila melanogaster, Caenorhabditis elegans, Saccharomyces cerevisiae, Arabidopsis thaliana), focusing on individual and multiple amino acid runs, charge and hydrophobic runs. We found that human proteins with multiple long runs are often associated with diseases; these include long glutamine runs that induce neurological disorders, various cancers, categories of leukemias (mostly involving chromosomal translocations), and an abundance of Ca(2 +) and K(+) channel proteins. Many human proteins with multiple runs function in development and/or transcription regulation and are Drosophila homeotic homologs. A large number of these proteins are expressed in the nervous system. More than 80% of Drosophila proteins with multiple runs seem to function in transcription regulation. The most frequent amino acid runs in Drosophila sequences occur for glutamine, alanine, and serine, whereas human sequences highlight glutamate, proline, and leucine. The most frequent runs in yeast are of serine, glutamine, and acidic residues. Compared with the other eukaryotic proteomes, amino acid runs are significantly more abundant in the fly. This finding might be interpreted in terms of innate differences in DNA-replication processes, repair mechanisms, DNA-modification systems, and mutational biases. There are striking differences in amino acid runs for glutamine, asparagine, and leucine among the five proteomes.