Analysis of among-site variation in substitution patterns

Substitution patterns among nucleotides are often assumed to be constant in phylogenetic analyses. Although variation in the average rate of substitution among sites is commonly accounted for, variation in the relative rates of specific types of substitution is not. Here, we review details of methodologies used for detecting and analyzing differences in substitution processes among predefined groups of sites. We describe how such analyses can be performed using existing phylogenetic tools, and discuss how new phylogenetic analysis tools we have recently developed can be used to provide more detailed and sensitive analyses, including study of the evolution of mutation and substitution processes. As an example we consider the mitochondrial genome, for which two types of transition deaminations (C⇒T and A⇒G) are strongly affected by single-strandedness during replication, resulting in a strand asymmetric mutation process. Since time spent single-stranded varies along the mitochondrial genome, their differential mutational response results in very different substitution patterns in different regions of the genome

Phylogenetic Codivergence Supports Coevolution of Mimetic Heliconius Butterflies

The unpalatable and warning-patterned butterflies _Heliconius erato_ and _Heliconius melpomene_ provide the best studied example of mutualistic Müllerian mimicry, thought – but rarely demonstrated – to promote coevolution. Some of the strongest available evidence for coevolution comes from phylogenetic codivergence, the parallel divergence of ecologically associated lineages. Early evolutionary reconstructions suggested codivergence between mimetic populations of _H. erato_ and _H. melpomene_, and this was initially hailed as the most striking known case of coevolution. However, subsequent molecular phylogenetic analyses found discrepancies in phylogenetic branching patterns and timing (topological and temporal incongruence) that argued against codivergence. We present the first explicit cophylogenetic test of codivergence between mimetic populations of _H. erato_ and _H. melpomene_, and re-examine the timing of these radiations. We find statistically significant topological congruence between multilocus coalescent population phylogenies of _H. erato_ and _H. melpomene_, supporting repeated codivergence of mimetic populations. Divergence time estimates, based on a Bayesian coalescent model, suggest that the evolutionary radiations of _H. erato_ and _H. melpomene_ occurred over the same time period, and are compatible with a series of temporally congruent codivergence events. This evidence supports a history of reciprocal coevolution between Müllerian co-mimics characterised by phylogenetic codivergence and parallel phenotypic change

Loss of Wolbachia infection during colonisation in the invasive Argentine ant Linepithema humile

WOLBACHIA are maternally inherited bacteria, which are very common in arthropods and nematodes. Wolbachia infection may affect host reproduction through feminisation, parthenogenesis, male-killing, cytoplasmic incompatibility and increased fecundity. Previous studies showing discrepancies between the phylogenies of Wolbachia and its arthropod hosts indicate that infection is frequently lost, but the causes of symbiont extinction have so far remained elusive. Here, we report data showing that colonisation of new habitats is a possible mechanism leading to the loss of infection. The presence and prevalence of Wolbachia were studied in three native and eight introduced populations of the Argentine ant Linepithema humile. The screening shows that the symbiont is common in the three native L. humile populations analysed. In contrast, Wolbachia was detected in only one of the introduced populations. The loss of infection associated with colonisation of new habitats may result from drift (founder effect) or altered selection pressures in the new habitat. Furthermore, a molecular phylogeny based on sequences of the Wolbachia wsp gene indicates that L. humile has been infected by a single strain. Horizontal transmission of the symbiont may be important in ants as suggested by the sequence similarity of strains in the three genera Linepithema, Acromyrmex, and Solenopsis native from South and Central America

GenomeBlast: a web tool for small genome comparison

BACKGROUND: Comparative genomics has become an essential approach for identifying homologous gene candidates and their functions, and for studying genome evolution. There are many tools available for genome comparisons. Unfortunately, most of them are not applicable for the identification of unique genes and the inference of phylogenetic relationships in a given set of genomes. RESULTS: GenomeBlast is a Web tool developed for comparative analysis of multiple small genomes. A new parameter called "coverage" was introduced and used along with sequence identity to evaluate global similarity between genes. With GenomeBlast, the following results can be obtained: (1) unique genes in each genome; (2) homologous gene candidates among compared genomes; (3) 2D plots of homologous gene candidates along the all pairwise genome comparisons; and (4) a table of gene presence/absence information and a genome phylogeny. We demonstrated the functions in GenomeBlast with an example of multiple herpesviral genome analysis and illustrated how GenomeBlast is useful for small genome comparison. CONCLUSION: We developed a Web tool for comparative analysis of small genomes, which allows the user not only to identify unique genes and homologous gene candidates among multiple genomes, but also to view their graphical distributions on genomes, and to reconstruct genome phylogeny. GenomeBlast runs on a Linux server with 4 CPUs and 4 GB memory. The online version of GenomeBlast is available to public by using a Web browser with the URL

9-Genes Reinforce the Phylogeny of Holometabola and Yield Alternate Views on the Phylogenetic Placement of Strepsiptera

Background: The extraordinary morphology, reproductive and developmental biology, and behavioral ecology of twisted wing parasites (order Strepsiptera) have puzzled biologists for centuries. Even today, the phylogenetic position of these enigmatic “insects from outer space” [1] remains uncertain and contentious. Recent authors have argued for the placement of Strepsiptera within or as a close relative of beetles (order Coleoptera), as sister group of flies (order Diptera), or even outside of Holometabola.Methodology/Principal Findings Here, we combine data from several recent studies with new data (for a total of 9 nuclear genes and ∼13 kb of aligned data for 34 taxa), to help clarify the phylogenetic placement of Strepsiptera. Our results unequivocally support the monophyly of Neuropteroidea ( = Neuropterida + Coleoptera) + Strepsiptera, but recover Strepsiptera either derived from within polyphagan beetles (order Coleoptera), or in a position sister to Neuropterida. All other supra-ordinal- and ordinal-level relationships recovered with strong nodal support were consistent with most other recent studies. Conclusions/Significance: These results, coupled with the recent proposed placement of Strepsiptera sister to Coleoptera, suggest that while the phylogenetic neighborhood of Strepsiptera has been identified, unequivocal placement to a specific branch within Neuropteroidea will require additional study.Organismic and Evolutionary Biolog

Comparative genomics of the class 4 histone deacetylase family indicates a complex evolutionary history

BACKGROUND: Histone deacetylases are enzymes that modify core histones and play key roles in transcriptional regulation, chromatin assembly, DNA repair, and recombination in eukaryotes. Three types of related histone deacetylases (classes 1, 2, and 4) are widely found in eukaryotes, and structurally related proteins have also been found in some prokaryotes. Here we focus on the evolutionary history of the class 4 histone deacetylase family. RESULTS: Through sequence similarity searches against sequenced genomes and expressed sequence tag data, we identified members of the class 4 histone deacetylase family in 45 eukaryotic and 37 eubacterial species representative of very distant evolutionary lineages. Multiple phylogenetic analyses indicate that the phylogeny of these proteins is, in many respects, at odds with the phylogeny of the species in which they are found. In addition, the eukaryotic members of the class 4 histone deacetylase family clearly display an anomalous phyletic distribution. CONCLUSION: The unexpected phylogenetic relationships within the class 4 histone deacetylase family and the anomalous phyletic distribution of these proteins within eukaryotes might be explained by two mechanisms: ancient gene duplication followed by differential gene losses and/or horizontal gene transfer. We discuss both possibilities in this report, and suggest that the evolutionary history of the class 4 histone deacetylase family may have been shaped by horizontal gene transfers

Complete Sequences of the Mitochondrial DNA of the Wild Gracilariopsis lemaneiformis and Two Mutagenic Cultivated Breeds (Gracilariaceae, Rhodophyta)

The complete mitochondrial DNA (mtDNA) of Gracilariopsis lemaneiformis was sequenced (25883 bp) and mapped to a circular model. The A+T composition was 72.5%. Forty six genes and two potentially functional open reading frames were identified. They include 24 protein-coding genes, 2 rRNA genes, 20 tRNA genes and 2 ORFs (orf60, orf142). There is considerable sequence synteny across the five red algal mtDNAs falling into Florideophyceae including Gr. lemaneiformis in this study and previously sequenced species. A long stem-loop and a hairpin structure were identified in intergenic regions of mt genome of Gr. lemaneiformis, which are believed to be involved with transcription and replication. In addition, the mtDNAs of two mutagenic cultivated breeds (“981” and “07-2”) were also sequenced. Compared with the mtDNA of wild Gr. lemaneiformis, the genome size and gene length and order of three strains were completely identical except nine base mutations including eight in the protein-coding genes and one in the tRNA gene. None of the base mutations caused frameshift or a premature stop codon in the mtDNA genes. Phylogenetic analyses based on mitochondrial protein-coding genes and rRNA genes demonstrated Gracilariopsis andersonii had closer phylogenetic relationship with its parasite Gracilariophila oryzoides than Gracilariopsis lemaneiformis which was from the same genus of Gracilariopsis

Daily Sampling of an HIV-1 Patient with Slowly Progressing Disease Displays Persistence of Multiple env Subpopulations Consistent with Neutrality

The molecular evolution of HIV-1 is characterized by frequent substitutions, indels and recombination events. In addition, a HIV-1 population may adapt through frequency changes of its variants. To reveal such population dynamics we analyzed HIV-1 subpopulation frequencies in an untreated patient with stable, low plasma HIV-1 RNA levels and close to normal CD4+ T-cell levels. The patient was intensively sampled during a 32-day period as well as approximately 1.5 years before and after this period (days −664, 1, 2, 3, 11, 18, 25, 32 and 522). 77 sequences of HIV-1 env (approximately 3100 nucleotides) were obtained from plasma by limiting dilution with 7–11 sequences per time point, except day −664. Phylogenetic analysis using maximum likelihood methods showed that the sequences clustered in six distinct subpopulations. We devised a method that took into account the relatively coarse sampling of the population. Data from days 1 through 32 were consistent with constant within-patient subpopulation frequencies. However, over longer time periods, i.e. between days 1…32 and 522, there were significant changes in subpopulation frequencies, which were consistent with evolutionarily neutral fluctuations. We found no clear signal of natural selection within the subpopulations over the study period, but positive selection was evident on the long branches that connected the subpopulations, which corresponds to >3 years as the subpopulations already were established when we started the study. Thus, selective forces may have been involved when the subpopulations were established. Genetic drift within subpopulations caused by de novo substitutions could be resolved after approximately one month. Overall, we conclude that subpopulation frequencies within this patient changed significantly over a time period of 1.5 years, but that this does not imply directional or balancing selection. We show that the short-term evolution we study here is likely representative for many patients of slow and normal disease progression

Molecular Evolution of Human Immunodeficiency Virus Type 1 upon Transmission between Human Leukocyte Antigen Disparate Donor-Recipient Pairs

BACKGROUND: To address evolution of HIV-1 after transmission, we studied sequence dynamics in and outside predicted epitopes of cytotoxic T lymphocytes (CTL) in subtype B HIV-1 variants that were isolated from 5 therapy-naive horizontal HLA-disparate donor-recipient pairs from the Amsterdam Cohort Studies on HIV-1 infection and AIDS. METHODOLOGY/PRINCIPAL FINDINGS: In the first weeks after transmission, the majority of donor-derived mutations in and outside donor-HLA-restricted epitopes in Gag, Env, and Nef, were preserved in the recipient. Reversion to the HIV-1 subtype B consensus sequence of mutations in- and outside donor-HLA-restricted CTL epitopes, and new mutations away from the consensus B sequence mostly within recipient-HLA-restricted epitopes, contributed equally to the early sequence changes. In the subsequent period (1-2 years) after transmission, still only a low number of both reverting and forward mutations had occurred. During subsequent long-term follow-up, sequence dynamics were dominated by forward mutations, mostly (50-85%) in recipient-HLA-restricted CTL epitopes. At the end of long-term follow-up, on average 43% of the transmitted CTL escape mutations in donor-HLA-restricted epitopes had reverted to the subtype B consensus sequence. CONCLUSIONS/SIGNIFICANCE: The relatively high proportion of long-term preserved mutations after transmission points to a lack of back selection even in the absence of CTL pressure, which may lead to an accumulating loss of critical CTL epitopes. Our data are supportive for a continuous adaptation of HIV-1 to host immune pressures which may have implications for vaccine design

Craniodental Morphology and Systematics of a New Family of Hystricognathous Rodents (Gaudeamuridae) from the Late Eocene and Early Oligocene of Egypt

BACKGROUND: Gaudeamus is an enigmatic hystricognathous rodent that was, until recently, known solely from fragmentary material from early Oligocene sites in Egypt, Oman, and Libya. Gaudeamus' molars are similar to those of the extant cane rat Thryonomys, and multiple authorities have aligned Gaudeamus with Thryonomys to the exclusion of other living and extinct African hystricognaths; recent phylogenetic analyses have, however, also suggested affinities with South American caviomorphs or Old World porcupines (Hystricidae). METHODOLOGY/PRINCIPAL FINDINGS: Here we describe the oldest known remains of Gaudeamus, including largely complete but crushed crania and complete upper and lower dentitions. Unlike younger Gaudeamus species, the primitive species described here have relatively complex occlusal patterns, and retain a number of plesiomorphic features. Unconstrained parsimony analysis nests Gaudeamus and Hystrix within the South American caviomorph radiation, implying what we consider to be an implausible back-dispersal across the Atlantic Ocean to account for Gaudeamus' presence in the late Eocene of Africa. An analysis that was constrained to recover the biogeographically more plausible hypothesis of caviomorph monophyly does not place Gaudeamus as a stem caviomorph, but rather as a sister taxon of hystricids. CONCLUSIONS/SIGNIFICANCE: We place Gaudeamus species in a new family, Gaudeamuridae, and consider it likely that the group originated, diversified, and then went extinct over a geologically brief period of time during the latest Eocene and early Oligocene in Afro-Arabia. Gaudeamurids are the only known crown hystricognaths from Afro-Arabia that are likely to be aligned with non-phiomorph members of that clade, and as such provide additional support for an Afro-Arabian origin of advanced stem and basal crown members of Hystricognathi