String kernels for protein sequence comparisons: improved fold recognition.

BACKGROUND: The amino acid sequence of a protein is the blueprint from which its structure and ultimately function can be derived. Therefore, sequence comparison methods remain essential for the determination of similarity between proteins. Traditional approaches for comparing two protein sequences begin with strings of letters (amino acids) that represent the sequences, before generating textual alignments between these strings and providing scores for each alignment. When the similitude between the two protein sequences to be compared is low however, the quality of the corresponding sequence alignment is usually poor, leading to poor performance for the recognition of similarity. RESULTS: In this study, we develop an alignment free alternative to these methods that is based on the concept of string kernels. Starting from recently proposed kernels on the discrete space of protein sequences (Shen et al, Found. Comput. Math., 2013,14:951-984), we introduce our own version, SeqKernel. Its implementation depends on two parameters, a coefficient that tunes the substitution matrix and the maximum length of k-mers that it includes. We provide an exhaustive analysis of the impacts of these two parameters on the performance of SeqKernel for fold recognition. We show that with the right choice of parameters, use of the SeqKernel similarity measure improves fold recognition compared to the use of traditional alignment-based methods. We illustrate the application of SeqKernel to inferring phylogeny on RNA polymerases and show that it performs as well as methods based on multiple sequence alignments. CONCLUSION: We have presented and characterized a new alignment free method based on a mathematical kernel for scoring the similarity of protein sequences. We discuss possible improvements of this method, as well as an extension of its applications to other modeling methods that rely on sequence comparison.

Lineage-specific variations in the trigger loop modulate RNA proofreading by bacterial RNA polymerases.

RNA cleavage by bacterial RNA polymerase (RNAP) has been implicated in transcriptional proofreading and reactivation of arrested transcription elongation complexes but its molecular mechanism is less understood than the mechanism of nucleotide addition, despite both reactions taking place in the same active site. RNAP from the radioresistant bacterium Deinococcus radiodurans is characterized by highly efficient intrinsic RNA cleavage in comparison with Escherichia coli RNAP. We find that the enhanced RNA cleavage activity largely derives from amino acid substitutions in the trigger loop (TL), a mobile element of the active site involved in various RNAP activities. The differences in RNA cleavage between these RNAPs disappear when the TL is deleted, or in the presence of GreA cleavage factors, which replace the TL in the active site. We propose that the TL substitutions modulate the RNA cleavage activity by altering the TL folding and its contacts with substrate RNA and that the resulting differences in transcriptional proofreading may play a role in bacterial stress adaptation.

Rapid phylogenetic and functional classification of short genomic fragments with signature peptides.

BACKGROUND: Classification is difficult for shotgun metagenomics data from environments such as soils, where the diversity of sequences is high and where reference sequences from close relatives may not exist. Approaches based on sequence-similarity scores must deal with the confounding effects that inheritance and functional pressures exert on the relation between scores and phylogenetic distance, while approaches based on sequence alignment and tree-building are typically limited to a small fraction of gene families. We describe an approach based on finding one or more exact matches between a read and a precomputed set of peptide 10-mers. RESULTS: At even the largest phylogenetic distances, thousands of 10-mer peptide exact matches can be found between pairs of bacterial genomes. Genes that share one or more peptide 10-mers typically have high reciprocal BLAST scores. Among a set of 403 representative bacterial genomes, some 20 million 10-mer peptides were found to be shared. We assign each of these peptides as a signature of a particular node in a phylogenetic reference tree based on the RNA polymerase genes. We classify the phylogeny of a genomic fragment (e.g., read) at the most specific node on the reference tree that is consistent with the phylogeny of observed signature peptides it contains. Using both synthetic data from four newly-sequenced soil-bacterium genomes and ten real soil metagenomics data sets, we demonstrate a sensitivity and specificity comparable to that of the MEGAN metagenomics analysis package using BLASTX against the NR database. Phylogenetic and functional similarity metrics applied to real metagenomics data indicates a signal-to-noise ratio of approximately 400 for distinguishing among environments. Our method assigns ~6.6 Gbp/hr on a single CPU, compared with 25 kbp/hr for methods based on BLASTX against the NR database. CONCLUSIONS: Classification by exact matching against a precomputed list of signature peptides provides comparable results to existing techniques for reads longer than about 300 bp and does not degrade severely with shorter reads. Orders of magnitude faster than existing methods, the approach is suitable now for inclusion in analysis pipelines and appears to be extensible in several different directions.

[Refinement of taxonomic position of Lactobacillus genus probiotic strains by 16S rDNA and rpoA gene sequencing].

AIM: Revision of the species identification of collection lactobacilli strains based on 16S rDNA and rpoA gene sequencing. MATERIALS AND METHODS: 52 lactobacilli cultures that present mostly Gamaleya Research Institute of Epidemiology and Microbiology (GIMC) collection were studied. 16S rDNA gene fragments were amplified by using Lb16a, Lb16b, 16S-midford, 16S-midrev primers. 2 different reverse primers were used for the analysis of rpoA gene depending on lactobacilli species. DNA fragments sequencing was performed with 3130 Genetic Analyzer (Applied Biosystems/Hitachi) with primers used for amplification. RESULTS: The effectiveness of sequencing of 2 targets for differentiation of species within lactobacilli phylogenetic groups was shown. Species diversity was demonstrated for GIMC lactobacilli strain collection that includes members of 9 species. All the strains marked previously as L. acidophilus were determined to belong to L. helveticus. Strains belonging to recently discovered L. farraginis species that has promising application in agriculture were detected. CONCLUSION: Genetic passports of original strains of 9 species of lactobacilli that are promising for further research.

tRNA-nucleotidyltransferases: highly unusual RNA polymerases with vital functions.

tRNA-nucleotidyltransferases are fascinating and unusual RNA polymerases responsible for the synthesis of the nucleotide triplet CCA at the 3'-terminus of tRNAs. As this CCA end represents an essential functional element for aminoacylation and translation, these polymerases (CCA-adding enzymes) are of vital importance in all organisms. With a possible origin of ancient telomerase-like activity, the CCA-adding enzymes obviously emerged twice during evolution, leading to structurally different, but functionally identical enzymes. The evolution as well as the unique polymerization features of these interesting proteins will be discussed in this review.

Molecular evidence for hybrid origins of the invasive polyploids Carthamus creticus and C. turkestanicus (Cardueae, Asteraceae).

The hexaploids Carthamus creticus and C. turkestanicus are noxious weeds with wide but non-overlapping Mediterranean distributions, and C. creticus, together with another polyploid, C. lanatus, have also invaded similar climatic regions in North and South America, South Africa and Australia. Here we infer their origins using sequences of the plastid intergenic spacer trnH-psbA and the intron trnK and three introns of nuclear low-copy genes of the RNA Polymerase family (RPD2 and the duplicated RPC2), as well as RAPD markers (random amplified polymorphic DNA). Phylogenetic analyses of the nuclear introns and additivity analysis of the RAPD markers support the hypotheses that the two hexaploids are allopolyploids sharing a tetraploid progenitor lineage represented by the broadly Mediterranean C. lanatus, combined with different diploid progenitor lineages consistent with the different geographic distributions of the hexaploids. Whereas C. leucocaulos from the south-eastern Greek Islands represents the diploid progenitor lineage of the western C. creticus, the Irano-Turanian C. glaucus represents the diploid progenitor lineage of the eastern C. turkestanicus. The plastid data suggest that the diploid lineages served as the maternal progenitors of the hexaploids.

Information theory based T7-like promoter models: classification of bacteriophages and differential evolution of promoters and their polymerases.

Molecular information theory was used to create sequence logos and promoter models for eight phages of the T7 group: T7, phiA1122, T3, phiYeO3-12, SP6, K1-5, gh-1 and K11. When these models were used to scan the corresponding genomes, a significant gap in the individual information distribution was observed between functional promoter sites and other sequences, suggesting that the models can be used to identify new T7-like promoters. When a combined 76-site model was used to scan the eight phages, 108 of the total 109 promoters were found, while none were found for other T7-like phages, phiKMV, P60, VpV262, SIO1, PaP3, Xp10, P-SSP7 and Ppu40, indicating that these phages do not belong to the T7 group. We propose that the T7-like transcription system, which consists of a phage-specific RNA polymerase and a set of conserved T7-like promoters, is a hallmark feature of the T7 group and can be used to classify T7-like phages. Phylogenetic trees of the T7-like promoter models and their corresponding RNA polymerases are similar, suggesting that the eight phages of the T7 group can be classified into five subgroups. However the SP6-like polymerases have apparently diverged from other polymerases more than their promoters have diverged from other promoters.

[Plastid RNA polymerases].

Plastids have a very interesting transcription apparatus that gives us an opportunity to investigate mono- and multisubunut RNA polymerase interaction under conditions of complex biogenesis of the organelles and the necessity to coordinate the expression of genes located in different cell compartments. The last decade has been a breakthrough in the study of chloroplast RNA polymerases. The most important advances are the discovery of nucleus-encoded RNA polymerase and nuclear gene family for sigma subunits of plastid-encoded RNA polymerase, the obtaining of knockout-mutants on the subunits of plastid-encoded RNA polymerase, and the revelation of this RNA polymerase reorganization during plastid biogenesis. The hypothesis concerning labor division between two plastid RNA polymerases has been proposed. The review focuses on the analysis of modern information about organization, function and evolution of plastid RNA polymerases.

Analysis of developing maize plastids reveals two mRNA stability classes correlating with RNA polymerase type.

The plastid genome is transcribed by two distinct RNA polymerases, the PEP encoded by the plastid genome and the NEP encoded in the nucleus. Initial models of plastid transcription held that the NEP is responsible for the transcription of housekeeping genes needed early in development, and that the PEP transcribes genes required for photosynthesis. Recently, this model was challenged by the discovery that all plastid genes are transcribed by NEP in PEP-deficient tobacco plastids, suggesting that mRNA turnover may have a strong role in previously observed transcription patterns. In this study, we provide evidence that the NEP enzyme level decreases as plastids mature. In contrast, production of mRNAs by NEP increases as plastids mature, yet their accumulations remain constant. These results suggest that as plastids mature NEP may become more active, and that mRNA turnover varies between transcripts synthesized by NEP and PEP.

[A role for Cajal bodies in assembly of the nuclear transcription machinery]. / Rol' telets Kakhala v sborke iadernogo apparata transkriptsii.

Cajal bodies (CBs) are small nuclear organelles that contain the three eukaryotic RNA polymerases and variety of factors involved in transcription and processing of all types of RNA. A number of these factors, as well as subunits of polymerase (pol) II itself, are rapidly and specifically targeted to CBs when injected into the cell. It is suggested that pol I, pol II and pol III transcription and processing complexes are preassembled in the CBs before transport to the sites of transcription on the chromosomes and in the nucleoli.

RpoA: a useful gene for phylogenetic analysis in diatoms.

The aim of this study was to compare the usefulness of two chloroplast-encoded genes (rpoA and rbcL) and the nuclear-encoded small subunit (SSU) ribosomal RNA for reconstructing phylogenetic relationships among diatoms at lower taxonomic levels. To this end, the rpoA and rbcL genes for selected centric and pennate diatoms were sequenced. The new rpoA and rbcL sequences, and an existing nuclear-encoded SSU rRNA data set, were subjected to weighted/unweighted parsimony, maximum likelihood, minimum evolution, and Bayesian analyses. All of the tree-building methods employed showed, based on the support values, that the rpoA gene was the most useful, relative to the rbcL and SSU rRNA genes, in determining phylogenetic relationships among the sampled diatoms. The support values for the relationships among the pennate lineages were, in many instances, greater in the rpoA trees than in the SSU rRNA trees. These results suggest that rpoA might be of value in determining phylogenetic relationships among pennate lineages.

N4 RNA polymerase II, a heterodimeric RNA polymerase with homology to the single-subunit family of RNA polymerases.

Bacteriophage N4 middle genes are transcribed by a phage-coded, heterodimeric, rifampin-resistant RNA polymerase, N4 RNA polymerase II (N4 RNAPII). Sequencing and transcriptional analysis revealed that the genes encoding the two subunits comprising N4 RNAPII are translated from a common transcript initiating at the N4 early promoter Pe3. These genes code for proteins of 269 and 404 amino acid residues with sequence similarity to the single-subunit, phage-like RNA polymerases. The genes encoding the N4 RNAPII subunits, as well as a synthetic construct encoding a fusion polypeptide, have been cloned and expressed. Both the individually expressed subunits and the fusion polypeptide reconstitute functional enzymes in vivo and in vitro.

Sequence analysis of an isolate from a fatal human infection of Australian bat lyssavirus.

Australian bat lyssavirus (ABLV), which occurs in pteropid and insectivorous bat populations, causes a rabies-like encephalitis in infected humans. We report the first complete sequence of an ABLV isolate obtained from a human who developed symptoms 27 months after being bitten by an infected flying fox. This isolate is the smallest lyssavirus to be sequenced, with a size of 11,918 nucleotides. Analyses of previously unsequenced regions and the complete genome confirm its close relationship with classical rabies viruses. In addition, a leucine zipper-like motif, not present in the other lyssaviruses, was found in the conserved domain I of the polymerase protein. This is the first report of a lyssavirus to vary in an 11-nucleotide, strictly conserved, complementary terminal sequence. This region is thought to encode important cis-acting regulatory signals; ABLV variation indicates a greater degree of flexibility than was thought for lyssaviruses in this region. A comparison of the pteropid and insectivorous isolates of ABLV indicates considerable differences between the two viruses. If the divergence of the two occurred on the Australian mainland, ABLV may have been endemic to Australia well before European colonisation.

Inhibition of transcription of class II, but not class III and I, genes in Xenopus postblastular embryos overexpressed with the TBP-binding protein, Dr1 (NC2beta).

Dr1 (NC2beta) is known to effectively repress transcription of class II genes, and much less effectively class III genes, but not class I genes through its binding to the TATA-binding protein (TBP), which is the major component of the basal transcription factor for polymerases II, III, and I. Previously, we isolated Xenopus Dr1 cDNA, and demonstrated that its mRNA is transcribed in oocytes and is inherited into early embryos, but its level decreases in later stages. Here, we overexpressed Xenopus Dr1 in Xenopus embryos and, found that the overexpression significantly reduces the levels of poly(A), cytoskeletal actin and histone H4 mRNAs, and the labeling of heterogeneous mRNA-like RNA in postblastular embryos, without affecting tRNA and rRNA syntheses. These results indicate that the overexpressed Dr1 specifically down-regulates the transcription of class II, but not class III and I, genes, and suggest that Dr1 plays an important role in the control of transcription in Xenopus embryogenesis.

Characterization of a gene encoding a single-subunit bacteriophage-type RNA polymerase from maize which is alternatively spliced.

Single-subunit RNA polymerases belonging to the T3/T7 bacteriophage family are thought to be common throughout eukaryotes. We report the isolation and characterization of a nucleus-encoded single-subunit RNA polymerase gene from maize. This gene is highly homologous to other single-subunit RNA polymerase genes from Arabidopsis, Chenopodium. yeast and Neurospora crassa involved in organellar transcription. Genomic Southern analysis reveals 10 to 15 hybridising fragments, suggesting that maize contains a small gene family. The isolated gene contains 19 exons and its genomic structure is highly conserved when compared to the three Arabidopsis homologues. Unlike the case in Arabidopsis, intron-12 of the maize bacteriophage-type RNA polymerase gene is alternatively spliced. Quantitative RT-PCR revealed that the resultant alternatively spliced transcript represents approximately 21 to 26% of the total polymerase mRNA in maize coleoptiles. The orthologous wheat bacteriophage-type RNA polymerase is also alternatively spliced and the intron exhibits 78% identity to maize intron-12. The conservation in alternative splicing between wheat and maize and its absence from Arabidopsis suggest a functional requirement for the alternatively spliced product.

The rpoN (sigma54) gene from Listeria monocytogenes is involved in resistance to mesentericin Y105, an antibacterial peptide from Leuconostoc mesenteroides.

To gain insight into the mode of action of mesentericin Y105, a bacteriocin bactericidal agent against Listeria monocytogenes, we undertook to identify the listerial factors mediating this susceptibility by using a genetic approach. Transposon mutants resistant to the bacteriocin were obtained. One of them corresponded to a transposon insertion in a gene (rpoN) encoding a putative protein (447 amino acids) with strong homologies to alternative transcriptional sigma54 factors, including that of Bacillus subtilis (38% identity). Complementation experiments with the wild-type rpoN gene demonstrated that the insertion in rpoN was responsible for the resistance phenotype in L. monocytogenes. Moreover, expression of the L. monocytogenes rpoN gene in an rpoN mutant strain of B. subtilis promoted transcription of a sigma54-dependent operon in the presence of the associated regulator. These results demonstrate that the L. monocytogenes rpoN gene encodes a new sigma54 factor.

Molecular systematic studies of eubacteria, using sigma70-type sigma factors of group 1 and group 2.

Sigma factors of the sigma70 family were used as a phylogenetic tool to compare evolutionary relationships among eubacteria. Several new sigma factor genes were cloned and sequenced to increase the variety of available sequences. Forty-two group 1 sigma factor sequences of various species were analyzed with the help of a distance matrix method to establish a phylogenetic tree. The tree derived by using sigma factors yielded subdivisions, including low-G+C and high-G+C gram-positive bacteria, cyanobacteria, and the alpha, beta, gamma, and delta subdivisions of proteobacteria, consistent with major bacterial groups found in trees derived from analyses with other molecules. However, some groupings (e.g., the chlamydiae, mycoplasmas, and green sulfur bacteria) are found in different positions than for trees obtained by using other molecular markers. A direct comparison to the most extensively used molecule in systematic studies, small-subunit rRNA, was made by deriving trees from essentially the same species set and using similar phylogenetic methods. Differences and similarities based on the two markers are discussed. Additionally, 31 group 2 sigma factors were analyzed in combination with the group 1 proteins in order to detect functional groupings of these alternative sigma factors. The data suggest that promoters recognized by the major vegetative sigma factors of eubacteria will contain sequence motifs and spacing very similar to those for the sigma70 sigma factors of Escherichia coli.

Structural and mechanistic relationships between nucleic acid polymerases.

A superfamily of nucleic acid polymerases that includes the pol I and pol alpha classes of DNA-directed DNA polymerases, mitochondrial and phage DNA-directed RNA polymerases, and most RNA-directed polymerases may be defined on the basis of the occurrence of conserved sequence motifs and tertiary structure similarities between HIV-1 reverse transcriptase, DNA polymerase I and T7 RNA polymerase. Although sequence or structural similarities do not yet justify inclusion of the multi-subunit DNA-directed RNA polymerases in this superfamily, mechanistic similarities suggest a deep relationship between these and the simpler T7-like RNA polymerases.