Evolutionary distances in the twilight zone -- a rational kernel approach

Phylogenetic tree reconstruction is traditionally based on multiple sequence alignments (MSAs) and heavily depends on the validity of this information bottleneck. With increasing sequence divergence, the quality of MSAs decays quickly. Alignment-free methods, on the other hand, are based on abstract string comparisons and avoid potential alignment problems. However, in general they are not biologically motivated and ignore our knowledge about the evolution of sequences. Thus, it is still a major open question how to define an evolutionary distance metric between divergent sequences that makes use of indel information and known substitution models without the need for a multiple alignment. Here we propose a new evolutionary distance metric to close this gap. It uses finite-state transducers to create a biologically motivated similarity score which models substitutions and indels, and does not depend on a multiple sequence alignment. The sequence similarity score is defined in analogy to pairwise alignments and additionally has the positive semi-definite property. We describe its derivation and show in simulation studies and real-world examples that it is more accurate in reconstructing phylogenies than competing methods. The result is a new and accurate way of determining evolutionary distances in and beyond the twilight zone of sequence alignments that is suitable for large datasets.Comment: to appear in PLoS ON

Potential conservation of circadian clock proteins in the phylum Nematoda as revealed by bioinformatic searches

Although several circadian rhythms have been described in C. elegans, its molecular clock remains elusive. In this work we employed a novel bioinformatic approach, applying probabilistic methodologies, to search for circadian clock proteins of several of the best studied circadian model organisms of different taxa (Mus musculus, Drosophila melanogaster, Neurospora crassa, Arabidopsis thaliana and Synechoccocus elongatus) in the proteomes of C. elegans and other members of the phylum Nematoda. With this approach we found that the Nematoda contain proteins most related to the core and accessory proteins of the insect and mammalian clocks, which provide new insights into the nematode clock and the evolution of the circadian system.Fil: Romanowski, Andrés. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Parque Centenario. Instituto de Investigaciones Bioquímicas de Buenos Aires. Fundación Instituto Leloir. Instituto de Investigaciones Bioquímicas de Buenos Aires; Argentina. Universidad Nacional de Quilmes. Departamento de Ciencia y Tecnología. Laboratorio de Cronobiología; ArgentinaFil: Garavaglia, Matías Javier. Universidad Nacional de Quilmes. Departamento de Ciencia y Tecnología. Laboratorio de Ing.genética y Biolog.molecular y Celular. Area Virus de Insectos; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Goya, María Eugenia. Universidad Nacional de Quilmes. Departamento de Ciencia y Tecnología. Laboratorio de Cronobiología; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Ghiringhelli, Pablo Daniel. Universidad Nacional de Quilmes. Departamento de Ciencia y Tecnología. Laboratorio de Ing.genética y Biolog.molecular y Celular. Area Virus de Insectos; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Golombek, Diego Andres. Universidad Nacional de Quilmes. Departamento de Ciencia y Tecnología. Laboratorio de Cronobiología; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentin

Accelerated Profile HMM Searches

Profile hidden Markov models (profile HMMs) and probabilistic inference methods have made important contributions to the theory of sequence database homology search. However, practical use of profile HMM methods has been hindered by the computational expense of existing software implementations. Here I describe an acceleration heuristic for profile HMMs, the “multiple segment Viterbi” (MSV) algorithm. The MSV algorithm computes an optimal sum of multiple ungapped local alignment segments using a striped vector-parallel approach previously described for fast Smith/Waterman alignment. MSV scores follow the same statistical distribution as gapped optimal local alignment scores, allowing rapid evaluation of significance of an MSV score and thus facilitating its use as a heuristic filter. I also describe a 20-fold acceleration of the standard profile HMM Forward/Backward algorithms using a method I call “sparse rescaling”. These methods are assembled in a pipeline in which high-scoring MSV hits are passed on for reanalysis with the full HMM Forward/Backward algorithm. This accelerated pipeline is implemented in the freely available HMMER3 software package. Performance benchmarks show that the use of the heuristic MSV filter sacrifices negligible sensitivity compared to unaccelerated profile HMM searches. HMMER3 is substantially more sensitive and 100- to 1000-fold faster than HMMER2. HMMER3 is now about as fast as BLAST for protein searches

Spag16, an Axonemal Central Apparatus Gene, Encodes a Male Germ Cell Nuclear Speckle Protein that Regulates SPAG16 mRNA Expression

Spag16 is the murine orthologue of Chlamydomonas reinhardtii PF20, a protein known to be essential to the structure and function of the “9+2” axoneme. In Chlamydomonas, the PF20 gene encodes a single protein present in the central pair of the axoneme. Loss of PF20 prevents central pair assembly/integrity and results in flagellar paralysis. Here we demonstrate that the murine Spag16 gene encodes two proteins: 71 kDa SPAG16L, which is found in all murine cells with motile cilia or flagella, and 35 kDa SPAG16S, representing the C terminus of SPAG16L, which is expressed only in male germ cells, and is predominantly found in specific regions within the nucleus that also contain SC35, a known marker of nuclear speckles enriched in pre-mRNA splicing factors. SPAG16S expression precedes expression of SPAG16L. Mice homozygous for a knockout of SPAG16L alone are infertile, but show no abnormalities in spermatogenesis. Mice chimeric for a mutation deleting the transcripts for both SPAG16L and SPAG16S have a profound defect in spermatogenesis. We show here that transduction of SPAG16S into cultured dispersed mouse male germ cells and BEAS-2B human bronchial epithelial cells increases SPAG16L expression, but has no effect on the expression of several other axoneme components. We also demonstrate that the Spag16L promoter shows increased activity in the presence of SPAG16S. The distinct nuclear localization of SPAG16S and its ability to modulate Spag16L mRNA expression suggest that SPAG16S plays an important role in the gene expression machinery of male germ cells. This is a unique example of a highly conserved axonemal protein gene that encodes two protein products with different functions

The YARHG Domain: An Extracellular Domain in Search of a Function

We have identified a new bacterial protein domain that we hypothesise binds to peptidoglycan. This domain is called the YARHG domain after the most highly conserved sequence-segment. The domain is found in the extracellular space and is likely to be composed of four alpha-helices. The domain is found associated with protein kinase domains, suggesting it is associated with signalling in some bacteria. The domain is also found associated with three different families of peptidases. The large number of different domains that are found associated with YARHG suggests that it is a useful functional module that nature has recombined multiple times

Diffractive Dijet Production at sqrt(s)=630 and 1800 GeV at the Fermilab Tevatron

We report a measurement of the diffractive structure function $F_{jj}^D$ of the antiproton obtained from a study of dijet events produced in association with a leading antiproton in $\bar pp$ collisions at $\sqrt s=630$ GeV at the Fermilab Tevatron. The ratio of $F_{jj}^D$ at $\sqrt s=630$ GeV to $F_{jj}^D$ obtained from a similar measurement at $\sqrt s=1800$ GeV is compared with expectations from QCD factorization and with theoretical predictions. We also report a measurement of the $\xi$ ($x$-Pomeron) and $\beta$ ($x$ of parton in Pomeron) dependence of $F_{jj}^D$ at $\sqrt s=1800$ GeV. In the region $0.035<\xi<0.095$, $|t|<1$ GeV$^2$ and $\beta<0.5$, $F_{jj}^D(\beta,\xi)$ is found to be of the form $\beta^{-1.0\pm 0.1} \xi^{-0.9\pm 0.1}$, which obeys $\beta$-$\xi$ factorization.Comment: LaTeX, 9 pages, Submitted to Phys. Rev. Letter

A Study of B0 -> J/psi K(*)0 pi+ pi- Decays with the Collider Detector at Fermilab

We report a study of the decays B0 -> J/psi K(*)0 pi+ pi-, which involve the creation of a u u-bar or d d-bar quark pair in addition to a b-bar -> c-bar(c s-bar) decay. The data sample consists of 110 1/pb of p p-bar collisions at sqrt{s} = 1.8 TeV collected by the CDF detector at the Fermilab Tevatron collider during 1992-1995. We measure the branching ratios to be BR(B0 -> J/psi K*0 pi+ pi-) = (8.0 +- 2.2 +- 1.5) * 10^{-4} and BR(B0 -> J/psi K0 pi+ pi-) = (1.1 +- 0.4 +- 0.2) * 10^{-3}. Contributions to these decays are seen from psi(2S) K(*)0, J/psi K0 rho0, J/psi K*+ pi-, and J/psi K1(1270)

Phamerator: a bioinformatic tool for comparative bacteriophage genomics

Background: Bacteriophage genomes have mosaic architectures and are replete with small open reading frames of unknown function, presenting challenges in their annotation, comparative analysis, and representation.Results: We describe here a bioinformatic tool, Phamerator, that assorts protein-coding genes into phamilies of related sequences using pairwise comparisons to generate a database of gene relationships. This database is used to generate genome maps of multiple phages that incorporate nucleotide and amino acid sequence relationships, as well as genes containing conserved domains. Phamerator also generates phamily circle representations of gene phamilies, facilitating analysis of the different evolutionary histories of individual genes that migrate through phage populations by horizontal genetic exchange.Conclusions: Phamerator represents a useful tool for comparative genomic analysis and comparative representations of bacteriophage genomes. © 2011 Cresawn et al; licensee BioMed Central Ltd

The Feasibility, Appropriateness, Meaningfulness, and Effectiveness of Parenting and Family Support Programs Delivered in the Criminal Justice System: A Systematic Review

Children whose parents are involved in the criminal justice system (CJS) are at increased risk of developing social, emotional, and behavioural difficulties and are more likely than their peers to become involved in the CJS themselves. Parenting behaviour and parent-child relationships have the potential to affect children’s outcomes with positive parenting practices having the potential to moderate some of the negative outcomes associated with parental involvement in the CJS. However, many parents in the CJS may lack appropriate role models to support the development of positive parenting beliefs and practices. Parenting programs offer an opportunity for parents to enhance their parenting knowledge and behaviours and improve relationships with children. Quantitative and qualitative evidence pertaining to the implementation and effectiveness of parenting programs delivered in the CJS was included. Five databases were searched and a total of 1145 articles were identified of which 29 met the review inclusion criteria. Overall, programs were found to significantly improve parenting attitudes; however, evidence of wider effects is limited. Additionally, the findings indicate that parenting programs can be meaningful for parents. Despite this, a number of challenges for implementation were found including the transient nature of the prison population and a lack of parent-child contact. Based on these findings, recommendations for the future development and delivery of programs are discussed

Species by Environment Interactions Affect Pyrrolizidine Alkaloid Expression in Senecio jacobaea, Senecio aquaticus, and Their Hybrids

We examined the effects of water and nutrient availability on the expression of the defense pyrrolizidine alkaloids (PAs) in Senecio jacobaea and S. aquaticus. Senecio jacobaea, and S. aquaticus are adapted to different natural habitats, characterized by differing abiotic conditions and different selection pressures from natural enemies. We tested if PA concentration and diversity are plastic over a range of water and nutrient treatments, and also whether such plasticity is dependent on plant species. We also tested the hypothesis that hybridization may contribute to PA diversity within plants, by comparing PA expression in parental species to that in artificially generated F1 hybrids, and also in later generation natural hybrids between S. jacobaea and S. aquaticus. We showed that total PA concentration in roots and shoots is not dependent on species, but that species determines the pattern of PA diversification. Pyrrolizidine alkaloid diversity and concentration are both dependent on environmental factors. Hybrids produce a putatively novel PA, and this PA is conserved in natural hybrids, that are backcrossed to S. jacobaea. Natural hybrids that are backcrossed several times to S. jacobaea are with regard to PA diversity significantly different from S. jacobaea but not from S. aquaticus, while F1 hybrids are in all cases more similar to S. jacobaea. These results collectively suggest that PA diversity is under the influence of natural selection