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

The extraordinary evolutionary history of the reticuloendotheliosis viruses

The reticuloendotheliosis viruses (REVs) comprise several closely related amphotropic retroviruses isolated from birds. These viruses exhibit several highly unusual characteristics that have not so far been adequately explained, including their extremely close relationship to mammalian retroviruses, and their presence as endogenous sequences within the genomes of certain large DNA viruses. We present evidence for an iatrogenic origin of REVs that accounts for these phenomena. Firstly, we identify endogenous retroviral fossils in mammalian genomes that share a unique recombinant structure with REVs—unequivocally demonstrating that REVs derive directly from mammalian retroviruses. Secondly, through sequencing of archived REV isolates, we confirm that contaminated Plasmodium lophurae stocks have been the source of multiple REV outbreaks in experimentally infected birds. Finally, we show that both phylogenetic and historical evidence support a scenario wherein REVs originated as mammalian retroviruses that were accidentally introduced into avian hosts in the late 1930s, during experimental studies of P. lophurae, and subsequently integrated into the fowlpox virus (FWPV) and gallid herpesvirus type 2 (GHV-2) genomes, generating recombinant DNA viruses that now circulate in wild birds and poultry. Our findings provide a novel perspective on the origin and evolution of REV, and indicate that horizontal gene transfer between virus families can expand the impact of iatrogenic transmission events

The All-Data-Based Evolutionary Hypothesis of Ciliated Protists with a Revised Classification of the Phylum Ciliophora (Eukaryota, Alveolata)

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/ The file attached is the published version of the article

A healthy school start - Parental support to promote healthy dietary habits and physical activity in children: Design and evaluation of a cluster-randomised intervention

<p>Abstract</p> <p>Background</p> <p>Childhood obesity is multi-factorial and determined to a large extent by dietary habits, physical activity and sedentary behaviours. Previous research has shown that school-based programmes are effective but that their effectiveness can be improved by including a parental component. At present, there is a lack of effective parental support programmes for improvement of diet and physical activity and prevention of obesity in children.</p> <p>Methods/Design</p> <p>This paper describes the rationale and design of a parental support programme to promote healthy dietary habits and physical activity in six-year-old children starting school. The study is performed in close collaboration with the school health care and is designed as a cluster-randomised controlled trial with a mixed methods approach. In total, 14 pre-school classes are included from a municipality in Stockholm county where there is large variation in socio-economic status between the families. The school classes are randomised to intervention (n = 7) and control (n = 7) groups including a total of 242 children. The intervention is based on social cognitive theory and consists of three main components: 1) a health information brochure; 2) two motivational interviewing sessions with the parents; and 3) teacher-led classroom activities with the children. The primary outcomes are physical activity in the children measured objectively by accelerometry, children's dietary and physical activity habits measured with a parent-proxy questionnaire and parents' self-efficacy measured by a questionnaire. Secondary outcomes are height, weight and waist circumference in the children. The duration of the intervention is six months and includes baseline, post intervention and six months follow-up measurements. Linear and logistic regression models will be used to analyse differences between intervention and control groups in the outcome variables. Mediator and moderator analysis will be performed. Participants will be interviewed.</p> <p>Discussion</p> <p>The results from this study will show if it is possible to promote a healthy lifestyle and a normal weight development among children from low-income districts with relatively limited efforts involving parents. Hopefully the study will provide new insights to the further development of effective programmes to prevent overweight and obesity in children.</p> <p>Trial registration</p> <p>ISRCTN: <a href="http://www.controlled-trials.com/ISRCTN32750699">ISRCTN32750699</a></p

Fatal Cases of Influenza A(H3N2) in Children: Insights from Whole Genome Sequence Analysis

During the Northern Hemisphere winter of 2003–2004 the emergence of a novel influenza antigenic variant, A/Fujian/411/2002-like(H3N2), was associated with an unusually high number of fatalities in children. Seventeen fatal cases in the UK were laboratory confirmed for Fujian/411-like viruses. To look for phylogenetic patterns and genetic markers that might be associated with increased virulence, sequencing and phylogenetic analysis of the whole genomes of 63 viruses isolated from fatal cases and non fatal “control” cases was undertaken. The analysis revealed the circulation of two main genetic groups, I and II, both of which contained viruses from fatal cases. No associated amino acid substitutions could be linked with an exclusive or higher occurrence in fatal cases. The Fujian/411-like viruses in genetic groups I and II completely displaced other A(H3N2) viruses, but they disappeared after 2004. This study shows that two A(H3N2) virus genotypes circulated exclusively during the winter of 2003–2004 in the UK and caused an unusually high number of deaths in children. Host factors related to immune state and differences in genetic background between patients may also play important roles in determining the outcome of an influenza infection

Long-Branch Attraction Bias and Inconsistency in Bayesian Phylogenetics

Bayesian inference (BI) of phylogenetic relationships uses the same probabilistic models of evolution as its precursor maximum likelihood (ML), so BI has generally been assumed to share ML's desirable statistical properties, such as largely unbiased inference of topology given an accurate model and increasingly reliable inferences as the amount of data increases. Here we show that BI, unlike ML, is biased in favor of topologies that group long branches together, even when the true model and prior distributions of evolutionary parameters over a group of phylogenies are known. Using experimental simulation studies and numerical and mathematical analyses, we show that this bias becomes more severe as more data are analyzed, causing BI to infer an incorrect tree as the maximum a posteriori phylogeny with asymptotically high support as sequence length approaches infinity. BI's long branch attraction bias is relatively weak when the true model is simple but becomes pronounced when sequence sites evolve heterogeneously, even when this complexity is incorporated in the model. This bias—which is apparent under both controlled simulation conditions and in analyses of empirical sequence data—also makes BI less efficient and less robust to the use of an incorrect evolutionary model than ML. Surprisingly, BI's bias is caused by one of the method's stated advantages—that it incorporates uncertainty about branch lengths by integrating over a distribution of possible values instead of estimating them from the data, as ML does. Our findings suggest that trees inferred using BI should be interpreted with caution and that ML may be a more reliable framework for modern phylogenetic analysis

Domain architecture evolution of pattern-recognition receptors

In animals, the innate immune system is the first line of defense against invading microorganisms, and the pattern-recognition receptors (PRRs) are the key components of this system, detecting microbial invasion and initiating innate immune defenses. Two families of PRRs, the intracellular NOD-like receptors (NLRs) and the transmembrane Toll-like receptors (TLRs), are of particular interest because of their roles in a number of diseases. Understanding the evolutionary history of these families and their pattern of evolutionary changes may lead to new insights into the functioning of this critical system. We found that the evolution of both NLR and TLR families included massive species-specific expansions and domain shuffling in various lineages, which resulted in the same domain architectures evolving independently within different lineages in a process that fits the definition of parallel evolution. This observation illustrates both the dynamics of the innate immune system and the effects of “combinatorially constrained” evolution, where existence of the limited numbers of functionally relevant domains constrains the choices of domain architectures for new members in the family, resulting in the emergence of independently evolved proteins with identical domain architectures, often mistaken for orthologs

Multi-Locus Phylogeographic and Population Genetic Analysis of Anolis carolinensis: Historical Demography of a Genomic Model Species

The green anole (Anolis carolinensis) has been widely used as an animal model in physiology and neurobiology but has recently emerged as an important genomic model. The recent sequencing of its genome has shed new light on the evolution of vertebrate genomes and on the process that govern species diversification. Surprisingly, the patterns of genetic diversity within natural populations of this widespread and abundant North American lizard remain relatively unknown. In the present study, we use 10 novel nuclear DNA sequence loci (N = 62 to 152) and one mitochondrial locus (N = 226) to delimit green anole populations and infer their historical demography. We uncovered four evolutionarily distinct and geographically restricted lineages of green anoles using phylogenetics, Bayesian clustering, and genetic distance methods. Molecular dating indicates that these lineages last shared a common ancestor ∼2 million years ago. Summary statistics and analysis of the frequency distributions of DNA polymorphisms strongly suggest range-wide expansions in population size. Using Bayesian Skyline Plots, we inferred the timing of population size expansions, which differ across lineages, and found evidence for a relatively recent and rapid westward expansion of green anoles across the Gulf Coastal Plain during the mid-Pleistocene. One surprising result is that the distribution of genetic diversity is not consistent with a latitudinal shift caused by climatic oscillations as is observed for many co-distributed taxa. This suggests that the most recent Pleistocene glacial cycles had a limited impact on the geographic distribution of the green anole at the northern limits of its range

Programmed DNA elimination of germline development genes in songbirds

In some eukaryotes, germline and somatic genomes differ dramatically in their composition. Here we characterise a major germline–soma dissimilarity caused by a germline-restricted chromosome (GRC) in songbirds. We show that the zebra finch GRC contains >115 genes paralogous to single-copy genes on 18 autosomes and the Z chromosome, and is enriched in genes involved in female gonad development. Many genes are likely functional, evidenced by expression in testes and ovaries at the RNA and protein level. Using comparative genomics, we show that genes have been added to the GRC over millions of years of evolution, with embryonic development genes bicc1 and trim71 dating to the ancestor of songbirds and dozens of other genes added very recently. The somatic elimination of this evolutionarily dynamic chromosome in songbirds implies a unique mechanism to minimise genetic conflict between germline and soma, relevant to antagonistic pleiotropy, an evolutionary process underlying ageing and sexual traits

Candidate chemoreceptor subfamilies differentially expressed in the chemosensory organs of the mollusc Aplysia

<p>Abstract</p> <p>Background</p> <p>Marine molluscs, as is the case with most aquatic animals, rely heavily on olfactory cues for survival. In the mollusc <it>Aplysia californica</it>, mate-attraction is mediated by a blend of water-borne protein pheromones that are detected by sensory structures called rhinophores. The expression of G protein and phospholipase C signaling molecules in this organ is consistent with chemosensory detection being via a G-protein-coupled signaling mechanism.</p> <p>Results</p> <p>Here we show that novel multi-transmembrane proteins with similarity to rhodopsin G-protein coupled receptors are expressed in sensory epithelia microdissected from the <it>Aplysia </it>rhinophore. Analysis of the <it>A. californica </it>genome reveals that these are part of larger multigene families that possess features found in metazoan chemosensory receptor families (that is, these families chiefly consist of single exon genes that are clustered in the genome). Phylogenetic analyses show that the novel <it>Aplysia </it>G-protein coupled receptor-like proteins represent three distinct monophyletic subfamilies. Representatives of each subfamily are restricted to or differentially expressed in the rhinophore and oral tentacles, suggesting that they encode functional chemoreceptors and that these olfactory organs sense different chemicals. Those expressed in rhinophores may sense water-borne pheromones. Secondary signaling component proteins Gα_q, Gα_i, and Gα_oare also expressed in the rhinophore sensory epithelium.</p> <p>Conclusion</p> <p>The novel rhodopsin G-protein coupled receptor-like gene subfamilies identified here do not have closely related identifiable orthologs in other metazoans, suggesting that they arose by a lineage-specific expansion as has been observed in chemosensory receptor families in other bilaterians. These candidate chemosensory receptors are expressed and often restricted to rhinophores and oral tentacles, lending support to the notion that water-borne chemical detection in <it>Aplysia </it>involves species- or lineage-specific families of chemosensory receptors.</p