The host innate immune response mediated by type I interferon (IFN) and the resulting up-regulation of hundreds of interferon-stimulated genes (ISGs) provide an immediate barrier to virus infection. Studies of the type I ‘interferome’ have mainly been carried out at a single species level, often lacking the power necessary to understand key evolutionary features of this pathway. Here, using a single experimental platform, we determined the properties of the interferomes of multiple vertebrate species and developed a webserver to mine the dataset. This approach revealed a conserved ‘core’ of 62 ISGs, including genes not previously associated with IFN, underscoring the ancestral functions associated with this antiviral host response. We show that gene expansion contributes to the evolution of the IFN system and that interferomes are shaped by lineage-specific pressures. Consequently, each mammal possesses a unique repertoire of ISGs, including genes common to all mammals and others unique to their specific species or phylogenetic lineages. An analysis of genes commonly down-regulated by IFN suggests that epigenetic regulation of transcription is a fundamental aspect of the IFN response. Our study provides a resource for the scientific community highlighting key paradigms of the type I IFN response.

Microsatellites are widely used in population genetics to uncover recent evolutionary events. They are typically genotyped using capillary sequencer, which capacity is usually limited to 9, at most 12 loci for each run, and which analysis is a tedious task that is performed by hand. With the rise of next-generation sequencing (NGS), a much larger number of loci and individuals are available from sequencing: for example, on a single run of a GS Junior, 28 loci from 96 individuals are sequenced with a 30X cover. We have developed an algorithm to automatically and efficiently genotype microsatellites from a collection of reads sorted by individual (e.g. specific PCR amplifications of a locus or a collection of reads that encompass a locus of interest). As the sequencing and the PCR amplification introduce artefactual insertions or deletions, the set of reads from a single microsatellite allele shows several length variants. The algorithm infers, without alignment, the true unknown allele(s) of each individual from the observed distributions of microsatellites length of all individuals. MicNeSs, a python implementation of the algorithm, can be used to genotype any microsatellite locus from any organism and has been tested on 454 pyrosequencing data of several loci from fruit flies (a model species) and red deers (a nonmodel species). Without any parallelization, it automatically genotypes 22 loci from 441 individuals in 11 hours on a standard computer. The comparison of MicNeSs inferences to the standard method shows an excellent agreement, with some differences illustrating the pros and cons of both methods.

Understanding multi-host pathogen maintenance and transmission dynamics is critical for disease control. However, transmission dynamics remain enigmatic largely because they are difficult to observe directly, particularly in wildlife. Here, we investigate the transmission dynamics of canine parvovirus (CPV) using state–space modelling of 20 years of CPV serology data from domestic dogs and African lions in the Serengeti ecosystem. We show that, although vaccination reduces the probability of infection in dogs, and despite indirect enhancement of population seropositivity as a result of vaccine shedding, the vaccination coverage achieved has been insufficient to prevent CPV from becoming widespread. CPV is maintained by the dog population and has become endemic with approximately 3.5-year cycles and prevalence reaching approximately 80%. While the estimated prevalence in lions is lower, peaks of infection consistently follow those in dogs. Dogs exposed to CPV are also more likely to become infected with a second multi-host pathogen, canine distemper virus. However, vaccination can weaken this coupling, raising questions about the value of monovalent versus polyvalent vaccines against these two pathogens. Our findings highlight the need to consider both pathogen- and host-level community interactions when seeking to understand the dynamics of multi-host pathogens and their implications for conservation, disease surveillance and control programmes.

Evolutionary events co-occurring along phylogenetic trees usually point to complex adaptive phenomena, sometimes implicating epistasis. While a number of methods have been developed to account for co-occurrence of events on the same internal or external branch of an evolutionary tree, there is a need to account for the larger diversity of possible relative positions of events in a tree. Here we propose a method to quantify to what extent two or more evolutionary events are associated on a phylogenetic tree. The method is applicable to any discrete character, like substitutions within a coding sequence or gains/losses of a biological function. Our method uses a general approach to statistically test for significant associations between events along the tree, which encompasses both events inseparable on the same branch, and events genealogically ordered on different branches. It assumes that the phylogeny and themapping of branches is known without errors. We address this problem from the statistical viewpoint by a linear algebra representation of the localization of the evolutionary events on the tree.We compute the full probability distribution of the number of paired events occurring in the same branch or in different branches of the tree, under a null model of independence where each type of event occurs at a constant rate uniformly inthephylogenetic tree. The strengths andweaknesses of themethodare assessed via simulations;we then apply the method to explore the loss of cell motility in intracellular pathogens.

An evolutionary process is reflected in the sequence of changes of any trait (e.g. morphological or molecular) through time. Yet, a better understanding of evolution would be procured by characterizing correlated evolution, or when two or more evolutionary processes interact. Previously developed parametric methods often require significant computing time as they rely on the estimation of many parameters. Here we propose a minimal likelihood framework modelling the joint evolution of two traits on a known phylogenetic tree. The type and strength of correlated evolution is characterized by a few parameters tuning mutation rates of each trait and interdependencies between these rates. The framework can be applied to study any discrete trait or character ranging from nucleotide substitution to gain or loss of a biological function. More specifically, it can be used to 1) test for independence between two evolutionary processes, 2) identify the type of interaction between them and 3) estimate parameter values of the most likely model of interaction. In the current implementation, the method takes as input a phylogenetic tree with discrete evolutionary events mapped on its branches. The method then maximizes the likelihood for one or several chosen scenarios. The strengths and limits of the method, as well as its relative power compared to a few other methods, are assessed using both simulations and data from 16S rRNA sequences in a sample of 54 γ-enterobacteria. We show that, even with datasets of fewer than 100 species, the method performs well in parameter estimation and in evolutionary model selection.

No abstract

This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

11Serology is a core component of the surveillance and management of viral zoonoses. Virus 12 neutralization tests are a gold standard serological diagnostic, but requirements for large volumes of 13 serum and high biosafety containment can limit widespread use. Here, focusing on Rabies lyssavirus, 14 a globally important zoonosis, we developed a pseudotype micro-neutralization rapid fluorescent 15focus inhibition test (pmRFFIT) that overcomes these limitations. Specifically, we adapted an existing 16 micro-neutralization test to use a green fluorescent protein-tagged murine leukemia virus 17 pseudotype in lieu of pathogenic rabies virus, reducing the need for specialized reagents for antigen 18 detection and enabling use in low-containment laboratories. We further used statistical analysis to 19 generate rapid, quantitative predictions of the probability and titer of rabies virus neutralizing 20antibodies from microscopic imaging of neutralization outcomes. Using 47 serum samples from 21 domestic dogs with neutralizing antibody titers estimated using the fluorescent antibody virus 22 neutralization test (FAVN), pmRFFIT showed moderate sensitivity (78.79%) and high specificity 23(84.62%). Despite small conflicts, titer predictions were correlated across tests repeated on different 24 dates both for dog samples (r = 0.93), and for a second dataset of sera from wild common vampire 25 bats (r = 0.72, N = 41), indicating repeatability. Our test uses a starting volume of 3.5 µL of serum, 26 estimates titers from a single dilution of serum rather than requiring multiple dilutions and end point 27 titration, and may be adapted to target neutralizing antibodies against alternative lyssavirus species. 28The pmRFFIT enables high-throughput detection of rabies virus neutralizing antibodies in low-29 biocontainment settings and is suited to studies in wild or captive animals where large serum volumes 30 cannot be obtained. 31