SlidingBayes: exploring recombination using a sliding window approach based on Bayesian phylogenetic inference

Likelihood-Based Tests of Species Tree Hypotheses

Likelihood-based tests of phylogenetic trees are a foundation of modern systematics. Over the past decade, an enormous wealth and diversity of model-based approaches have been developed for phylogenetic inference of both gene trees and species trees. However, while many techniques exist for conducting formal likelihood-based tests of gene trees, such frameworks are comparatively underdeveloped and underutilized for testing species tree hypotheses. To date, widely used tests of tree topology are designed to assess the fit of classical models of molecular sequence data and individual gene trees and thus are not readily applicable to the problem of species tree inference. To address this issue, we derive several analogous likelihood-based approaches for testing topologies using modern species tree models and heuristic algorithms that use gene tree topologies as input for maximum likelihood estimation under the multispecies coalescent. For the purpose of comparing support for species trees, these tests leverage the statistical procedures of their original gene tree-based counterparts that have an extended history for testing phylogenetic hypotheses at a single locus. We discuss and demonstrate a number of applications, limitations, and important considerations of these tests using simulated and empirical phylogenomic data sets that include both bifurcating topologies and reticulate network models of species relationships. Finally, we introduce the open-source R package SpeciesTopoTestR (SpeciesTopology Tests in R) that includes a suite of functions for conducting formal likelihood-based tests of species topologies given a set of input gene tree topologies.

Phylogeny and character evolution in the Dacrymycetes, and systematics of Unilacrymaceae and Dacryonaemataceae fam. nov

We present a multilocus phylogeny of the class Dacrymycetes, based on data from the 18S, ITS, 28S, RPB1, RPB2, TEF-1α, 12S, and ATP6 DNA regions, with c. 90 species including the types of most currently accepted genera. A variety of methodological approaches was used to infer phylogenetic relationships among the Dacrymycetes, from a supermatrix strategy using maximum likelihood and Bayesian inference on a concatenated dataset, to coalescence-based calculations, such as quartet-based summary methods of independent single-locus trees, and Bayesian integration of single-locus trees into a species tree under the multispecies coalescent. We evaluate for the first time the taxonomic usefulness of some cytological phenotypic characters, i.e., vacuolar contents (vacuolar bodies and lipid bodies), number of nuclei of recently discharged basidiospores, and pigments, with especial emphasis on carotenoids. These characters, along with several others traditionally used for the taxonomy of this group (basidium shape, presence and morphology of clamp connections, morphology of the terminal cells of cortical/marginal hyphae, presence and degree of ramification of the hyphidia), are mapped on the resulting phylogenies and their evolution through the class Dacrymycetes discussed. Our analyses reveal five lineages that putatively represent five different families, four of which are accepted and named. Three out of these four lineages correspond to previously circumscribed and published families (Cerinomycetaceae, Dacrymycetaceae, and Unilacrymaceae), and one is proposed as the new family Dacryonaemataceae. Provisionally, only a single order, Dacrymycetales, is accepted within the class. Furthermore, the systematics of the two smallest families, Dacryonaemataceae and Unilacrymaceae, are investigated to the species level, using coalescence-based species delimitation on multilocus DNA data, and a detailed morphological study including morphometric analyses of the basidiospores. Three species are accepted in Dacryonaema, the type, Da. rufum, the newly combined Da. macnabbii (basionym Dacrymyces macnabbii), and a new species named Da. macrosporum. Two species are accepted in Unilacryma, the new U. bispora, and the type, U. unispora, the latter treated in a broad sense pending improved sampling across the Holarctic.

Nonstructural Protein L* Species Specificity Supports a Mouse Origin for Vilyuisk Human Encephalitis Virus

Vilyuisk human encephalitis virus (VHEV) is a picornavirus related to Theiler's murine encephalomyelitis virus (TMEV). VHEV was isolated from human material passaged in mice. Whether this VHEV is of human or mouse origin is therefore unclear. We took advantage of the species-specific activity of the nonstructural L* protein of theiloviruses to track the origin of TMEV isolates. TMEV L* inhibits RNase L, the effector enzyme of the interferon pathway. By using coimmunoprecipitation and functional RNase L assays, the species specificity of RNase L antagonism was tested for L* from mouse (DA) and rat (RTV-1) TMEV strains as well as for VHEV. Coimmunoprecipitation and functional assay data confirmed the species specificity of L* activity and showed that L* from rat strain RTV-1 inhibited rat but not mouse or human RNase L. Next, we showed that the VHEV L* protein was phylogenetically related to L* of mouse viruses and that it failed to inhibit human RNase L but readily antagonized mouse RNase L, unambiguously showing the mouse origin of VHEV.IMPORTANCE Defining the natural host of a virus can be a thorny issue, especially when the virus was isolated only once or when the isolation story is complex. The species Theilovirus includes Theiler's murine encephalomyelitis virus (TMEV), infecting mice and rats, and Saffold virus (SAFV), infecting humans. One TMEV strain, Vilyuisk human encephalitis virus (VHEV), however, was isolated from mice that were inoculated with cerebrospinal fluid of a patient presenting with chronic encephalitis. It is therefore unclear whether VHEV was derived from the human sample or from the inoculated mouse. The L* protein encoded by TMEV inhibits RNase L, a cellular enzyme involved in innate immunity, in a species-specific manner. Using binding and functional assays, we show that this species specificity even allows discrimination between TMEV strains of mouse and of rat origins. The VHEV L* protein clearly inhibited mouse but not human RNase L, indicating that this virus originates from mice.

Introgression browser: high-throughput whole-genome SNP visualization

Breeding by introgressive hybridization is a pivotal strategy to broaden the genetic basis of crops. Usually, the desired traits are monitored in consecutive crossing generations by marker-assisted selection, but their analyses fail in chromosome regions where crossover recombinants are rare or not viable. Here, we present the Introgression Browser (iBrowser), a bioinformatics tool aimed at visualizing introgressions at nucleotide or SNP (Single Nucleotide Polymorphisms) accuracy. The software selects homozygous SNPs from Variant Call Format (VCF) information and filters out heterozygous SNPs, multi-nucleotide polymorphisms (MNPs) and insertion-deletions (InDels). For data analysis iBrowser makes use of sliding windows, but if needed it can generate any desired fragmentation pattern through General Feature Format (GFF) information. In an example of tomato (Solanum lycopersicum) accessions we visualize SNP patterns and elucidate both position and boundaries of the introgressions. We also show that our tool is capable of identifying alien DNA in a panel of the closely related S. pimpinellifolium by examining phylogenetic relationships of the introgressed segments in tomato. In a third example, we demonstrate the power of the iBrowser in a panel of 597 Arabidopsis accessions, detecting the boundaries of a SNP-free region around a polymorphic 1.17 Mbp inverted segment on the short arm of chromosome 4. The architecture and functionality of iBrowser makes the software appropriate for a broad set of analyses including SNP mining, genome structure analysis, and pedigree analysis. Its functionality, together with the capability to process large data sets and efficient visualization of sequence variation, makes iBrowser a valuable breeding tool.

The role of selection in shaping diversity of natural M. tuberculosis populations

Mycobacterium tuberculosis (M.tb), the cause of tuberculosis (TB), is estimated to infect a new host every second. While analyses of genetic data from natural populations of M.tb have emphasized the role of genetic drift in shaping patterns of diversity, the influence of natural selection on this successful pathogen is less well understood. We investigated the effects of natural selection on patterns of diversity in 63 globally extant genomes of M.tb and related pathogenic mycobacteria. We found evidence of strong purifying selection, with an estimated genome-wide selection coefficient equal to -9.5 × 10(-4) (95% CI -1.1 × 10(-3) to -6.8 × 10(-4)); this is several orders of magnitude higher than recent estimates for eukaryotic and prokaryotic organisms. We also identified different patterns of variation across categories of gene function. Genes involved in transport and metabolism of inorganic ions exhibited very low levels of non-synonymous polymorphism, equivalent to categories under strong purifying selection (essential and translation-associated genes). The highest levels of non-synonymous variation were seen in a group of transporter genes, likely due to either diversifying selection or local selective sweeps. In addition to selection, we identified other important influences on M.tb genetic diversity, such as a 25-fold expansion of global M.tb populations coincident with explosive growth in human populations (estimated timing 1684 C.E., 95% CI 1620-1713 C.E.). These results emphasize the parallel demographic histories of this obligate pathogen and its human host, and suggest that the dominant effect of selection on M.tb is removal of novel variants, with exceptions in an interesting group of genes involved in transportation and defense. We speculate that the hostile environment within a host imposes strict demands on M.tb physiology, and thus a substantial fitness cost for most new mutations. In this respect, obligate bacterial pathogens may differ from other host-associated microbes such as symbionts.

Phylogenetic and epidemic modeling of rapidly evolving infectious diseases

Epidemic modeling of infectious diseases has a long history in both theoretical and empirical research. However the recent explosion of genetic data has revealed the rapid rate of evolution that many populations of infectious agents undergo and has underscored the need to consider both evolutionary and ecological processes on the same time scale. Mathematical epidemiology has applied dynamical models to study infectious epidemics, but these models have tended not to exploit--or take into account--evolutionary changes and their effect on the ecological processes and population dynamics of the infectious agent. On the other hand, statistical phylogenetics has increasingly been applied to the study of infectious agents. This approach is based on phylogenetics, molecular clocks, genealogy-based population genetics and phylogeography. Bayesian Markov chain Monte Carlo and related computational tools have been the primary source of advances in these statistical phylogenetic approaches. Recently the first tentative steps have been taken to reconcile these two theoretical approaches. We survey the Bayesian phylogenetic approach to epidemic modeling of infection diseases and describe the contrasts it provides to mathematical epidemiology as well as emphasize the significance of the future unification of these two fields.

Towards an accurate identification of mosaic genes and partial horizontal gene transfers

Many bacteria and viruses adapt to varying environmental conditions through the acquisition of mosaic genes. A mosaic gene is composed of alternating sequence polymorphisms either belonging to the host original allele or derived from the integrated donor DNA. Often, the integrated sequence contains a selectable genetic marker (e.g. marker allowing for antibiotic resistance). An effective identification of mosaic genes and detection of corresponding partial horizontal gene transfers (HGTs) are among the most important challenges posed by evolutionary biology. We developed a method for detecting partial HGT events and related intragenic recombination giving rise to the formation of mosaic genes. A bootstrap procedure incorporated in our method is used to assess the support of each predicted partial gene transfer. The proposed method can be also applied to confirm or discard complete (i.e. traditional) horizontal gene transfers detected by any HGT inferring method. While working on a full-genome scale, the new method can be used to assess the level of mosaicism in the considered genomes as well as the rates of complete and partial HGT underlying their evolution.

GiRaF: robust, computational identification of influenza reassortments via graph mining

Reassortments in the influenza virus—a process where strains exchange genetic segments—have been implicated in two out of three pandemics of the 20th century as well as the 2009 H1N1 outbreak. While advances in sequencing have led to an explosion in the number of whole-genome sequences that are available, an understanding of the rate and distribution of reassortments and their role in viral evolution is still lacking. An important factor in this is the paucity of automated tools for confident identification of reassortments from sequence data due to the challenges of analyzing large, uncertain viral phylogenies. We describe here a novel computational method, called GiRaF (Graph-incompatibility-based Reassortment Finder), that robustly identifies reassortments in a fully automated fashion while accounting for uncertainties in the inferred phylogenies. The algorithms behind GiRaF search large collections of Markov chain Monte Carlo (MCMC)-sampled trees for groups of incompatible splits using a fast biclique enumeration algorithm coupled with several statistical tests to identify sets of taxa with differential phylogenetic placement. GiRaF correctly finds known reassortments in human, avian, and swine influenza populations, including the evolutionary events that led to the recent ‘swine flu’ outbreak. GiRaF also identifies several previously unreported reassortments via whole-genome studies to catalog events in H5N1 and swine influenza isolates.

A complex hepatitis B virus (X/C) recombinant is common in Long An county, Guangxi and may have originated in southern China

Recently, a complex (X/C) hepatitis B virus (HBV) recombinant, first reported in 2000, was proposed as a new genotype; although this was refuted immediately because the strains differ by less than 8 % in nucleotide distance from genotype C. Over 13.5 % (38/281) of HBV isolates from the Long An cohort in China were not assigned to a specific genotype, using current genotyping tools to analyse surface ORF sequences, and these have about 98 % similarity to the X/C recombinants. To determine whether this close identity extends to the full-length sequences and to investigate the evolutionary history of the Long An X/C recombinants, 17 complete genome sequences were determined. They are highly similar (96–99 %) to the Vietnamese strains and, although some reach or exceed 8 % nucleotide sequence difference from all known genotypes, they cluster together in the same clade, separating in a phylogenetic tree from the genotype C branch. Analysis of recombination reveals that all but one of the Long An isolates resembles the Vietnamese isolates in that they result from apparent recombination between genotype C and a parent of unknown genotype (X), which shows similarity in part to genotype G. The exception, isolate QL523, has a greater proportion of genotype C parent. Phylogeographic analysis reveals that these recombinants probably arose in southern China and spread later to Vietnam and Laos.

The phylogenetic information profile of HIV-1 and the degradation effect of recombination

HIV-1, while known to recombine frequently and evolve rapidly, is one of the most sequenced organisms. The availability of many and long sequences (almost full-length) renders HIV-1 as a good model for studying theoretical predictions linked to evolution and phylogenetic inferences. Here we study the effects of rapid and through-recombination evolution on phylogenetic information in order to confirm theoretical predictions of the characteristics of phylogenetic information on a real dataset. Firstly we study the fluctuation of the phylogenetic information along the HIV-1 genome showing that genomic regions such as the first part and the last part of the pol gene contain less phylogenetic information, while the vpr, vpu and the first exon of the tat gene contain more phylogenetic information compared to the rest of the genome. Moreover, we provide evidence that phylogenetic information is correlated to the sequence similarity of the dataset used and is degraded by the effect of recombination.

The evolution of transposon repeat-induced point mutation in the genome of Colletotrichum cereale: reconciling sex, recombination and homoplasy in an ''asexual" pathogen

Mobile transposable elements are among the primary drivers of the evolution of eukaryotic genomes. For fungi, repeat-induced point mutation (RIP) silencing minimizes deleterious effects of transposons by mutating multicopy DNA during meiosis. In this study we identify five transposon species from the mitosporic fungus Colletotrichum cereale and report the signature pattern of RIP acting in a lineage-specific manner on 21 of 35 unique transposon copies, providing the first evidence for sexual recombination for this species. Sequence analysis of genomic populations of the retrotransposon Ccret2 showed repeated rounds of RIP mutation acting on different copies of the element. In the RIPped Ccret2 population, there were multiple inferences of incongruence primarily attributed to RIP-induced homoplasy. This study supports the view that the sequence variability of transposon populations in filamentous fungi reflects the activities of evolutionary processes that fall outside of typical phylogenetic or population genetic reconstructions.

Recombination confounds the early evolutionary history of human immunodeficiency virus type 1: subtype G is a circulating recombinant form

Human immunodeficiency virus type 1 (HIV-1) is classified in nine subtypes (A to D, F, G, H, J, and K), a number of subsubtypes, and several circulating recombinant forms (CRFs). Due to the high level of genetic diversity within HIV-1 and to its worldwide distribution, this classification system is widely used in fields as diverse as vaccine development, evolution, epidemiology, viral fitness, and drug resistance. Here, we demonstrate how the high recombination rates of HIV-1 may confound the study of its evolutionary history and classification. Our data show that subtype G, currently classified as a pure subtype, has in fact a recombinant history, having evolved following recombination between subtypes A and J and a putative subtype G parent. In addition, we find no evidence for recombination within one of the lineages currently classified as a CRF, CRF02_AG. Our analysis indicates that CRF02_AG was the parent of the recombinant subtype G, rather than the two having the opposite evolutionary relationship, as is currently proposed. Our results imply that the current classification of HIV-1 subtypes and CRFs is an artifact of sampling history, rather than reflecting the evolutionary history of the virus. We suggest a reanalysis of all pure subtypes and CRFs in order to better understand how high rates of recombination have influenced HIV-1 evolutionary history.

Spider monkey, Muriqui and Woolly monkey relationships revisited

The taxonomic relationships among the four genera of the Atelidae family, Alouatta (Howler), Ateles (Spider), Lagothrix (Woolly) and Brachyteles (Muriqui), have been the subject of great debate. In general, almost all authors agree with the assignment of Howler monkeys as the basal genus, either in its own tribe Alouattini or in the subfamily Alouattinae, but they disagree on the associations among the other members of the family. Muriquis have been grouped with Spider monkeys based on the fact that they share various behavioral and morphological characteristics. Cladistic analyses using morphological, biochemical, karyotype and behavioral characteristics depicted a phylogenetic tree that places Howler as the basal genus and the remaining genera in an unresolved politomy. More recent studies using molecular data have suggested that Muriqui and Woolly monkeys are sister groups. However, a recent study based on nuclear and mtDNA argued that politomy is what best represents the relationships among Spider, Woolly and Muriqui. To contribute to this debate we have added new data from two nuclear genes, Transferrin and von Willebrand Factor, and using an alignment of 17,997 bp we demonstrate that a total analysis strongly supports the Muriqui-Woolly clade. A gene-to-gene approach showed that four of the eight nuclear genes provide support for the Muriqui-Woolly clade, two strongly and two moderately, while none of the eight genes provide support for any alternative arrangement. The mitochondrial genes were not able to resolve the politomy. A possible reason for the difficulty in resolving atelid relationships may be the short period of time separating each cladogenetic event in the evolutionary process that shaped this family.

Cytochrome P450 1A4 and 1A5 in Common Cormorant (Phalacrocorax carbo): Evolutionary Relationships and Functional Implications Associated with Dioxin and Related Compounds

The present study characterized cytochrome P4501A (CYP1A) isoforms from common cormorant (Phalacrocorax carbo) with regard to their evolutionary relationships and their roles in disposition of dioxin and related compounds (DRCs). Two clones isolated from a cormorant liver cDNA library were named CYP1A4 and CYP1A5 on the basis of greatest overall amino acid identity shared with chicken (Gallus gallus) CYP1A4 (78%) and CYP1A5 (78%), respectively. Spatial heterogeneity in phylogenetic signal along the sequences strongly indicated that cormorant CYP1A4 and CYP1A5 have undergone partial interparalog gene conversion, similar to chicken and mammalian CYP1As. Phylogenetic analysis of a putatively unconverted region produced a tree topology consistent with the orthology of avian CYP1A5s with mammalian CYP1A2s and avian CYP1A4s with mammalian CYP1A1s. Hepatic CYP1A4 and CYP1A5 mRNA levels in wild cormorants from Lake Biwa, Japan, were quantified to examine the effects of DRCs on isoform-specific expression and to evaluate the toxicokinetics of DRCs in which CYP1A expression is involved. Both CYP1A4 and CYP1A5 mRNA levels were positively correlated with total tetrachlorodibenzo-p-dioxin toxic equivalents and concentrations of each congener in most cases in the liver, suggesting the induction of both enzymes through a shared transcriptional mechanism. The lack of correlation of 2,3,7,8-tetrachlorodibenzofuran and 3,3',4,4'-tetrachlorobiphenyl (PCB77) to CYP1A gene expression is likely due to the rapid metabolism of these two congeners. Liver-to-muscle concentration ratios for most DRC congeners except PCB77 and mono-ortho coplanar polychlorinated biphenyls significantly increased with an elevation of CYP1A4 and CYP1A5 mRNA levels. The present data suggest that hepatic sequestration of some DRCs occurs in cormorant via binding to either CYP1A5 or both CYP1A4 and CYP1A5.

A Revised Evolutionary History of the CYP1A Subfamily: Gene Duplication, Gene Conversion, and Positive Selection

Members of cytochrome P450 subfamily 1A (CYP1As) are involved in detoxification and bioactivation of common environmental pollutants. Understanding the functional evolution of these genes is essential to predicting and interpreting species differences in sensitivity to toxicity caused by such chemicals. The CYP1A gene subfamily comprises a single ancestral representative in most fish species and two paralogs in higher vertebrates, including birds and mammals. Phylogenetic analysis of complete coding sequences suggests that mammalian and bird paralog pairs (CYP1A1/2 and CYP1A4/5, respectively) are the result of independent gene duplication events. However, comparison of vertebrate genome sequences revealed that CYP1A genes lie within an extended region of conserved fine-scale synteny, suggesting that avian and mammalian CYP1A paralogs share a common genomic history. Algorithms designed to detect recombination between nucleotide sequences indicate that gene conversion has homogenized most of the length of the chicken CYP1A genes, as well as the 5' end of mammalian CYP1As. Together, these data indicate that avian and mammalian CYP1A paralog pairs resulted from a single gene duplication event and that extensive gene conversion is responsible for the exceptionally high degree of sequence similarity between CYP1A4 and CYP1A5. Elevated nonsynonymous/synonymous substitution ratios within a putatively unconverted stretch of approximately 250 bp suggests that positive selection may have reduced the effective rate of gene conversion in this region, which contains two substrate recognition sites. This work significantly alters our understanding of functional evolution in the CYP1A subfamily, suggesting that gene conversion and positive selection have been the dominant processes of sequence evolution.