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. The type I interferon (IFN) response is triggered upon sensing of an incoming pathogen in an infected cell and results in the expression of hundreds of IFN-stimulated genes (ISGs, collectively referred to as ‘the interferome’). Studies on the interferome have been carried out mainly in human cells and therefore often lack the power to understand comparative evolutionary aspects of this critical pathway. In this study, we characterized the interferome in several animal species (including humans) using a single experimental framework. This approach allowed us to identify fundamental properties of the innate immune system. In particular, we revealed 62 ‘core’ ISGs, up-regulated in response to IFN in all vertebrates, highlighting the ancestral functions of the IFN system. In addition, we show that many genes repressed by the IFN response normally function as regulators of cell transcription. ISGs shared by multiple species have a higher propensity than other genes to exist as multiple copies in the genome. Importantly, we observed that genes have arisen as ISGs throughout evolution. Hence, every animal species possesses a unique repertoire of ISGs that includes core and lineage-specific genes. Collectively, our data provide a framework on which it will be possible to test the role of the IFN response in pathogen emergence and cross-species transmission.

Advances in DNA sequencing have dramatically increased the rate at which new viruses are being identified. However, the host species associations of most virus sequences identified in metagenomic samples are difficult to determine. Our analysis indicates that viruses proposed to infect vertebrates (in some cases being linked to human disease) may in fact be restricted to arthropod hosts. The detection of these sequences in vertebrate samples may reflect their widespread presence in the environment as viruses of parasitic arthropods.

Parvoviruses (family Parvoviridae) are small, single-stranded DNA viruses. Many parvoviral pathogens of medical, veterinary and ecological importance have been identified. In this study, we used high-throughput sequencing (HTS) to investigate the diversity of parvoviruses infecting wild and domestic animals in Brazil. We identified 21 parvovirus sequences (including twelve nearly complete genomes and nine partial genomes) in samples derived from rodents, bats, opossums, birds and cattle in Pernambuco, São Paulo, Paraná and Rio Grande do Sul states. These sequences were investigated using phylogenetic and distance-based approaches and were thereby classified into eight parvovirus species (six of which have not been described previously), representing six distinct genera in the subfamily Parvovirinae. Our findings extend the known biogeographic range of previously characterized parvovirus species and the known host range of three parvovirus genera (Dependovirus, Aveparvovirus and Tetraparvovirus). Moreover, our investigation provides a window into the ecological dynamics of parvovirus infections in vertebrates, revealing that many parvovirus genera contain well-defined sub-lineages that circulate widely throughout the world within particular taxonomic groups of hosts.

Circoviruses (family Circoviridae) are small, non-enveloped viruses that have short, single-stranded DNA genomes. Circovirus sequences are frequently recovered in metagenomic investigations, indicating that these viruses are widespread, yet they remain relatively poorly understood. Endogenous circoviral elements (CVe) are DNA sequences derived from circoviruses that occur in vertebrate genomes. CVe are a useful source of information about the biology and evolution of circoviruses. In this study, we screened 362 vertebrate genome assemblies in silico to generate a catalog of CVe loci. We identified a total of 179 CVe sequences, most of which have not been reported previously. We show that these CVe loci reflect at least 19 distinct germline integration events. We determine the structure of CVe loci, identifying some that show evidence of potential functionalization. We also identify orthologous copies of CVe in snakes, fish, birds, and mammals, allowing us to add new calibrations to the timeline of circovirus evolution. Finally, we observed that some ancient CVe group robustly with contemporary circoviruses in phylogenies, with all sequences within these groups being derived from the same host class or order, implying a hitherto underappreciated stability in circovirus-host relationships. The openly available dataset constructed in this investigation provides new insights into circovirus evolution, and can be used to facilitate further studies of circoviruses and CVe.

Genomic data contribute invaluable information to the epidemiological investigation of pathogens of public health importance. However, whole-genome sequencing (WGS) of bacteria typically relies on culture, which represents a major hurdle for generating such data for a wide range of species for which culture is challenging. In this study, we assessed the use of culture-free target-enrichment sequencing as a method for generating genomic data for two bacterial species: (1)Bacillus anthracis,which causes anthrax in both people and animals and whose culture requires high-level containment facilities; and (2)Mycoplasma amphoriforme, a fastidious emerging human respiratory pathogen. We obtained high-quality genomic data for both species directly from clinical samples, with sufficient coverage (>15×) for confident variant calling over at least 80% of the baited genomes for over two thirds of the samples tested. Higher qPCR cycle threshold (Ct) values (indicative of lower pathogen concentrations in the samples), pooling libraries prior to capture, and lower captured library concentration were all statistically associated with lower capture efficiency. TheCtvalue had the highest predictive value, explaining 52 % of the variation in capture efficiency. Samples withCtvalues ≤30 were over six times more likely to achieve the threshold coverage than those with aCt> 30. We conclude that target-enrichment sequencing provides a valuable alternative to standard WGS following bacterial culture and creates opportunities for an improved understanding of the epidemiology and evolution of many clinically important pathogens for which culture is challenging.

Genomic sequencing has revolutionized our understanding of bacterial disease epidemiology, but remains underutilized for zoonotic pathogens in remote endemic settings. Anthrax, caused by the spore-forming bacterium Bacillus anthracis , remains a threat to human and animal health and rural livelihoods in low- and middle-income countries. While the global genomic diversity of B. anthracis has been well-characterized, there is limited information on how its populations are genetically structured at the scale at which transmission occurs, critical for understanding the pathogen’s evolution and transmission dynamics. Using a uniquely rich dataset, we quantified genome-wide SNPs among 73 B. anthracis isolates derived from 33 livestock carcasses sampled over 1 year throughout the Ngorongoro Conservation Area, Tanzania, a region hyperendemic for anthrax. Genome-wide SNPs distinguished 22 unique B. anthracis genotypes (i.e. SNP profiles) within the study area. However, phylogeographical structure was lacking, as identical SNP profiles were found throughout the study area, likely the result of the long and variable periods of spore dormancy and long-distance livestock movements. Significantly, divergent genotypes were obtained from spatio-temporally linked cases and even individual carcasses. The high number of SNPs distinguishing isolates from the same host is unlikely to have arisen during infection, as supported by our simulation models. This points to an unexpectedly wide transmission bottleneck for B. anthracis , with an inoculum comprising multiple variants being the norm. Our work highlights that inferring transmission patterns of B. anthracis from genomic data will require analytical approaches that account for extended and variable environmental persistence, as well as co-infection.