Search

Malayan pangolin SARS-CoV-2-related coronavirus (SARSr-CoV-2) is closely related to SARS-CoV-2. However, little is known about its pathogenicity in pangolins. Using CT scans we show that SARSr-CoV-2 positive Malayan pangolins are characterized by bilateral ground-glass opacities in lungs in a similar manner to COVID-19 patients. Histological examination and blood gas tests are indicative of dyspnea. SARSr-CoV-2 infected multiple organs in pangolins, with the lungs the major target, and histological expression data revealed that ACE2 and TMPRSS2 were co-expressed with viral RNA. Transcriptome analysis indicated that virus-positive pangolins were likely to have inadequate interferon responses, with relative greater cytokine and chemokine activity in the lung and spleen. Notably, both viral RNA and viral proteins were detected in three pangolin fetuses, providing initial evidence for vertical virus transmission. In sum, our study outlines the biological framework of SARSr-CoV-2 in pangolins, revealing striking similarities to COVID-19 in humans.

Most of our understanding of the ecology and evolution of avian influenza A virus (AIV) in wild birds is derived from studies conducted in the northern hemisphere on waterfowl, with a substantial bias towards dabbling ducks. However, relevant environmental conditions and patterns of avian migration and reproduction are substantially different in the southern hemisphere. Through the sequencing and analysis of 333 unique AIV genomes collected from wild birds collected over 15 years we show that Australia is a global sink for AIV diversity and not integrally linked with the Eurasian gene pool. Rather, AIV are infrequently introduced to Australia, followed by decades of isolated circulation and eventual extinction. The number of co-circulating viral lineages varies per subtype. AIV haemagglutinin (HA) subtypes that are rarely identified at duck-centric study sites (H8-12) had more detected introductions and contemporary co-circulating lineages in Australia. Combined with a lack of duck migration beyond the Australian-Papuan region, these findings suggest introductions by long-distance migratory shorebirds. In addition, on the available data we found no evidence of directional or consistent patterns in virus movement across the Australian continent. This feature corresponds to patterns of bird movement, whereby waterfowl have nomadic and erratic rainfall-dependant distributions rather than consistent intra-continental migratory routes. Finally, we detected high levels of virus gene segment reassortment, with a high diversity of AIV genome constellations across years and locations. These data, in addition to those from other studies in Africa and South America, clearly show that patterns of AIV dynamics in the Southern Hemisphere are distinct from those in the temperate north.

At the end of 2019 Wuhan witnessed an outbreak of "atypical pneumonia" that later developed into a global pandemic. Metagenomic sequencing rapidly revealed the causative agent of this outbreak to be a novel coronavirus denoted SARS-CoV-2. To provide a snapshot of the pathogens in pneumonia-associated respiratory samples from Wuhan prior to the emergence of SARS-CoV-2, we collected bronchoalveolar lavage fluid samples from 408 patients presenting with pneumonia and acute respiratory infections at the Central Hospital of Wuhan between 2016 and 2017. Unbiased total RNA sequencing was performed to reveal their "total infectome", including viruses, bacteria and fungi. We identified 35 pathogen species, comprising 13 RNA viruses, 3 DNA viruses, 16 bacteria and 3 fungi, often at high abundance and including multiple co-infections (13.5%). SARS-CoV-2 was not present. These data depict a stable core infectome comprising common respiratory pathogens such as rhinoviruses and influenza viruses, an atypical respiratory virus (EV-D68), and a single case of a sporadic zoonotic pathogen-Chlamydia psittaci. Samples from patients experiencing respiratory disease on average had higher pathogen abundance than healthy controls. Phylogenetic analyses of individual pathogens revealed multiple origins and global transmission histories, highlighting the connectedness of the Wuhan population. This study provides a comprehensive overview of the pathogens associated with acute respiratory infections and pneumonia, which were more diverse and complex than obtained using targeted PCR or qPCR approaches. These data also suggest that SARS-CoV-2 or closely related viruses were absent from Wuhan in 2016-2017.

Identifying the animal reservoirs from which zoonotic viruses will likely emerge is central to understanding the determinants of disease emergence. Accordingly, there has been an increase in studies attempting zoonotic "risk assessment." Herein, we demonstrate that the virological data on which these analyses are conducted are incomplete, biased, and rapidly changing with ongoing virus discovery. Together, these shortcomings suggest that attempts to assess zoonotic risk using available virological data are likely to be inaccurate and largely only identify those host taxa that have been studied most extensively. We suggest that virus surveillance at the human-animal interface may be more productive.

Ticks (order: Ixodida) are a highly diverse and ecologically important group of ectoparasitic blood-feeding organisms. One such species, the seabird tick (Ixodes uriae), is widely distributed around the circumpolar regions of the northern and southern hemispheres. It has been suggested that Ix. uriae spread from the southern to the northern circumpolar region millions of years ago and has remained isolated in these regions ever since. Such a profound biographic subdivision provides a unique opportunity to determine whether viruses associated with ticks exhibit the same evolutionary patterns as their hosts. To test this, we collected Ix. uriae specimens near a Gentoo penguin (Pygoscelis papua) colony at Neko harbour, Antarctica, and from migratory birds-the Razorbill (Alca torda) and the Common murre (Uria aalge)-on Bonden island, northern Sweden. Through meta-transcriptomic next-generation sequencing we identified 16 RNA viruses, seven of which were novel. Notably, we detected the same species, Ronne virus, and two closely related species, Bonden virus and Piguzov virus, in both hemispheres indicating that there have been at least two cross-circumpolar dispersal events. Similarly, we identified viruses discovered previously in other locations several decades ago, including Gadgets Gully virus, Taggert virus and Okhotskiy virus. By identifying the same or closely related viruses in geographically disjunct sampling locations we provide evidence for virus dispersal within and between the circumpolar regions. In marked contrast, our phylogenetic analysis revealed no movement of the Ix. uriae tick hosts between the same locations. Combined, these data suggest that migratory birds are responsible for the movement of viruses at both local and global scales.

BackgroundIn this large-scale cluster-randomized controlled trial (cRCT) we sought to assess the effectiveness of facemasks against viral respiratory infections.Methods and resultsOver three consecutive Hajj seasons (2013, 2014, 2015) pilgrims' tents in Makkah were allocated to 'facemask' or 'no facemask' group. Fifty facemasks were offered to participants in intervention tents, to be worn over four days, and none were offered to participants in control tents. All participants recorded facemask use and respiratory symptoms in health diaries. Nasal swabs were collected from the symptomatic for virus detection by reverse transcription polymerase chain reaction. Clinical symptoms and laboratory results were analyzed by 'intention- to-treat' and 'per-protocol'. A total of 7687 adult participants from 318 tents were randomized: 3864 from 149 tents to the intervention group, and 3823 from 169 tents to the control group. Participants were aged 18 to 95 (median 34, mean 37) years, with a male to female ratio of 1:1.2. Overall, respiratory viruses were detected in 277 of 650 (43%) nasal/pharyngeal swabs collected from symptomatic pilgrims. Common viruses were rhinovirus (35.1%), influenza (4.5%) and parainfluenza (1.7%). In the intervention arm, respectively 954 (24.7%) and 1842 (47.7%) participants used facemasks daily and intermittently, while in the control arm, respectively 546 (14.3%) and 1334 (34.9%) used facemasks daily and intermittently. By intention-to-treat analysis, facemask use did not seem to be effective against laboratory-confirmed viral respiratory infections (odds ratio [OR], 1.4; 95% confidence interval [CI], 0.9 to 2.1, p = 0.18) nor against clinical respiratory infection (OR, 1.1; 95% CI, 0.9 to 1.4, p = 0.40). Similarly, in a per-protocol analysis, facemask use did not seem to be effective against laboratory-confirmed viral respiratory infections (OR 1.2, 95% CI 0.9-1.7, p = 0.26) nor against clinical respiratory infection (OR 1.3, 95% CI 1.0-1.8, p = 0.06).ConclusionThis trial was unable to provide conclusive evidence on facemask efficacy against viral respiratory infections most likely due to poor adherence to protocol.

Dengue virus (DENV) transmission by mosquitoes is a time-dependent process that begins with the consumption of an infectious blood-meal. DENV infection then proceeds stepwise through the mosquito from the midgut to the carcass, and ultimately to the salivary glands, where it is secreted into saliva and then transmitted anew on a subsequent bite. We examined viral kinetics in tissues of the Aedes aegypti mosquito over a finely graded time course, and as per previous studies, found that initial viral dose and serotype strain diversity control infectivity. We also found that a threshold level of virus is required to establish body-wide infections and that replication kinetics in the early and intermediate tissues do not predict those of the salivary glands. Our findings have implications for mosquito GMO design, modeling the contribution of transmission to vector competence and the role of mosquito kinetics in the overall DENV epidemiological landscape.

Eukaryotes of the genus Plasmodium cause malaria, a parasitic disease responsible for substantial morbidity and mortality in humans. Yet, the nature and abundance of any viruses carried by these divergent eukaryotic parasites is unknown. We investigated the Plasmodium virome by performing a meta-transcriptomic analysis of blood samples taken from patients suffering from malaria and infected with P. vivax, P. falciparum or P. knowlesi. This resulted in the identification of a narnavirus-like sequence, encoding an RNA polymerase and restricted to P. vivax samples, as well as an associated viral segment of unknown function. These data, confirmed by PCR, are indicative of a novel RNA virus that we term Matryoshka RNA virus 1 (MaRNAV-1) to reflect its analogy to a "Russian doll": a virus, infecting a parasite, infecting an animal. Additional screening revealed that MaRNAV-1 was abundant in geographically diverse P. vivax derived from humans and mosquitoes, strongly supporting its association with this parasite, and not in any of the other Plasmodium samples analyzed here nor Anopheles mosquitoes in the absence of Plasmodium. Notably, related bi-segmented narnavirus-like sequences (MaRNAV-2) were retrieved from Australian birds infected with a Leucocytozoon-a genus of eukaryotic parasites that group with Plasmodium in the Apicomplexa subclass hematozoa. Together, these data support the establishment of two new phylogenetically divergent and genomically distinct viral species associated with protists, including the first virus likely infecting Plasmodium parasites. As well as broadening our understanding of the diversity and evolutionary history of the eukaryotic virosphere, the restriction to P. vivax may be of importance in understanding P. vivax-specific biology in humans and mosquitoes, and how viral co-infection might alter host responses at each stage of the P. vivax life-cycle.

The development of high-throughput genome sequencing enables accurate measurements of levels of sub-consensus intra-host virus genetic diversity and analysis of the role played by natural selection during cross-species transmission. We analysed the natural and experimental evolution of rabies virus (RABV), an important example of a virus that is able to make multiple host jumps. In particular, we (i) analyzed RABV evolution during experimental host switching with the goal of identifying possible genetic markers of host adaptation, (ii) compared the mutational changes observed during passage with those observed in natura, and (iii) determined whether the colonization of new hosts or tissues requires adaptive evolution in the virus. To address these aims, animal infection models (dog and fox) and primary cell culture models (embryo brain cells of dog and fox) were developed and viral variation was studied in detail through deep genome sequencing. Our analysis revealed a strong unidirectional host evolutionary effect, as dog-adapted rabies virus was able to replicate in fox and fox cells relatively easily, while dogs or neuronal dog cells were not easily susceptible to fox adapted-RABV. This suggests that dog RABV may be able to adapt to some hosts more easily than other host variants, or that when RABV switched from dogs to red foxes it lost its ability to adapt easily to other species. Although no difference in patterns of mutation variation between different host organs was observed, mutations were common following both in vitro and in vivo passage. However, only a small number of these mutations also appeared in natura, suggesting that adaptation during successful cross-species virus transmission is a complex, multifactorial evolutionary process.

[This corrects the article DOI: 10.1371/journal.ppat.1006750.].

[This corrects the article DOI: 10.1371/journal.ppat.1006780.].

Hepatitis B virus (HBV) is a ubiquitous viral pathogen associated with large-scale morbidity and mortality in humans. However, there is considerable uncertainty over the time-scale of its origin and evolution. Initial shotgun data from a mid-16th century Italian child mummy, that was previously paleopathologically identified as having been infected with Variola virus (VARV, the agent of smallpox), showed no DNA reads for VARV yet did for hepatitis B virus (HBV). Previously, electron microscopy provided evidence for the presence of VARV in this sample, although similar analyses conducted here did not reveal any VARV particles. We attempted to enrich and sequence for both VARV and HBV DNA. Although we did not recover any reads identified as VARV, we were successful in reconstructing an HBV genome at 163.8X coverage. Strikingly, both the HBV sequence and that of the associated host mitochondrial DNA displayed a nearly identical cytosine deamination pattern near the termini of DNA fragments, characteristic of an ancient origin. In contrast, phylogenetic analyses revealed a close relationship between the putative ancient virus and contemporary HBV strains (of genotype D), at first suggesting contamination. In addressing this paradox we demonstrate that HBV evolution is characterized by a marked lack of temporal structure. This confounds attempts to use molecular clock-based methods to date the origin of this virus over the time-frame sampled so far, and means that phylogenetic measures alone cannot yet be used to determine HBV sequence authenticity. If genuine, this phylogenetic pattern indicates that the genotypes of HBV diversified long before the 16th century, and enables comparison of potential pathogenic similarities between modern and ancient HBV. These results have important implications for our understanding of the emergence and evolution of this common viral pathogen.

The factors that determine the pattern and rate of spread of influenza virus at a continental-scale are uncertain. Although recent work suggests that influenza epidemics in the United States exhibit a strong geographical correlation, the spatiotemporal dynamics of influenza in Australia, a country and continent of approximately similar size and climate complexity but with a far smaller population, are not known. Using a unique combination of large-scale laboratory-confirmed influenza surveillance comprising >450,000 entries and genomic sequence data we determined the local-level spatial diffusion of this important human pathogen nationwide in Australia. We used laboratory-confirmed influenza data to characterize the spread of influenza virus across Australia during 2007-2016. The onset of established epidemics varied across seasons, with highly synchronized epidemics coinciding with the emergence of antigenically distinct viruses, particularly during the 2009 A/H1N1 pandemic. The onset of epidemics was largely synchronized between the most populous cities, even those separated by distances of >3000 km and those that experience vastly diverse climates. In addition, by analyzing global phylogeographic patterns we show that the synchronized dissemination of influenza across Australian cities involved multiple introductions from the global influenza population, coupled with strong domestic connectivity, rather than through the distinct radial patterns of geographic dispersal that are driven by work-flow transmission as observed in the United States. In addition, by comparing the spatial structure of influenza A and B, we found that these viruses tended to occupy different geographic regions, and peak in different seasons, perhaps indicative of moderate cross-protective immunity or viral interference effects. The highly synchronized outbreaks of influenza virus at a continental-scale revealed here highlight the importance of coordinated public health responses in the event of the emergence of a novel, human-to-human transmissible, virus.

The co-evolution of myxoma virus (MYXV) and the European rabbit occurred independently in Australia and Europe from different progenitor viruses. Although this is the canonical study of the evolution of virulence, whether the genomic and phenotypic outcomes of MYXV evolution in Europe mirror those observed in Australia is unknown. We addressed this question using viruses isolated in the United Kingdom early in the MYXV epizootic (1954-1955) and between 2008-2013. The later UK viruses fell into three distinct lineages indicative of a long period of separation and independent evolution. Although rates of evolutionary change were almost identical to those previously described for MYXV in Australia and strongly clock-like, genome evolution in the UK and Australia showed little convergence. The phenotypes of eight UK viruses from three lineages were characterized in laboratory rabbits and compared to the progenitor (release) Lausanne strain. Inferred virulence ranged from highly virulent (grade 1) to highly attenuated (grade 5). Two broad disease types were seen: cutaneous nodular myxomatosis characterized by multiple raised secondary cutaneous lesions, or an amyxomatous phenotype with few or no secondary lesions. A novel clinical outcome was acute death with pulmonary oedema and haemorrhage, often associated with bacteria in many tissues but an absence of inflammatory cells. Notably, reading frame disruptions in genes defined as essential for virulence in the progenitor Lausanne strain were compatible with the acquisition of high virulence. Combined, these data support a model of ongoing host-pathogen co-evolution in which multiple genetic pathways can produce successful outcomes in the field that involve both different virulence grades and disease phenotypes, with alterations in tissue tropism and disease mechanisms.

The cross-species transmission of viruses from one host species to another is responsible for the majority of emerging infections. However, it is unclear whether some virus families have a greater propensity to jump host species than others. If related viruses have an evolutionary history of co-divergence with their hosts there should be evidence of topological similarities between the virus and host phylogenetic trees, whereas host jumping generates incongruent tree topologies. By analyzing co-phylogenetic processes in 19 virus families and their eukaryotic hosts we provide a quantitative and comparative estimate of the relative frequency of virus-host co-divergence versus cross-species transmission among virus families. Notably, our analysis reveals that cross-species transmission is a near universal feature of the viruses analyzed here, with virus-host co-divergence occurring less frequently and always on a subset of viruses. Despite the overall high topological incongruence among virus and host phylogenies, the Hepadnaviridae, Polyomaviridae, Poxviridae, Papillomaviridae and Adenoviridae, all of which possess double-stranded DNA genomes, exhibited more frequent co-divergence than the other virus families studied here. At the other extreme, the virus and host trees for all the RNA viruses studied here, particularly the Rhabdoviridae and the Picornaviridae, displayed high levels of topological incongruence, indicative of frequent host switching. Overall, we show that cross-species transmission plays a major role in virus evolution, with all the virus families studied here having the potential to jump host species, and that increased sampling will likely reveal more instances of host jumping.

Pregnancy complications are poorly represented in the archeological record, despite their importance in contemporary and ancient societies. While excavating a Byzantine cemetery in Troy, we discovered calcified abscesses among a woman’s remains. Scanning electron microscopy of the tissue revealed ‘ghost cells’, resulting from dystrophic calcification, which preserved ancient maternal, fetal and bacterial DNA of a severe infection, likely chorioamnionitis. Gardnerella vaginalis and Staphylococcus saprophyticus dominated the abscesses. Phylogenomic analyses of ancient, historical, and contemporary data showed that G. vaginalis Troy fell within contemporary genetic diversity, whereas S. saprophyticus Troy belongs to a lineage that does not appear to be commonly associated with human disease today. We speculate that the ecology of S. saprophyticus infection may have differed in the ancient world as a result of close contacts between humans and domesticated animals. These results highlight the complex and dynamic interactions with our microbial milieu that underlie severe maternal infections.

The natural evolution of rabies virus (RABV) provides a potent example of multiple host shifts and an important opportunity to determine the mechanisms that underpin viral emergence. Using 321 genome sequences spanning an unprecedented diversity of RABV, we compared evolutionary rates and selection pressures in viruses sampled from multiple primary host shifts that occurred on various continents. Two major phylogenetic groups, bat-related RABV and dog-related RABV, experiencing markedly different evolutionary dynamics were identified. While no correlation between time and genetic divergence was found in bat-related RABV, the evolution of dog-related RABV followed a generally clock-like structure, although with a relatively low evolutionary rate. Subsequent molecular clock dating indicated that dog-related RABV likely underwent a rapid global spread following the intensification of intercontinental trade starting in the 15th century. Strikingly, although dog RABV has jumped to various wildlife species from the order Carnivora, we found no clear evidence that these host-jumping events involved adaptive evolution, with RABV instead characterized by strong purifying selection, suggesting that ecological processes also play an important role in shaping patterns of emergence. However, specific amino acid changes were associated with the parallel emergence of RABV in ferret-badgers in Asia, and some host shifts were associated with increases in evolutionary rate, particularly in the ferret-badger and mongoose, implying that changes in host species can have important impacts on evolutionary dynamics.

The emergence of severe acute respiratory syndrome coronavirus (SARS-CoV) in 2002 and Middle East respiratory syndrome coronavirus (MERS-CoV) in 2012 has generated enormous interest in the biodiversity, genomics and cross-species transmission potential of coronaviruses, especially those from bats, the second most speciose order of mammals. Herein, we identified a novel coronavirus, provisionally designated Rousettus bat coronavirus GCCDC1 (Ro-BatCoV GCCDC1), in the rectal swab samples of Rousettus leschenaulti bats by using pan-coronavirus RT-PCR and next-generation sequencing. Although the virus is similar to Rousettus bat coronavirus HKU9 (Ro-BatCoV HKU9) in genome characteristics, it is sufficiently distinct to be classified as a new species according to the criteria defined by the International Committee of Taxonomy of Viruses (ICTV). More striking was that Ro-BatCoV GCCDC1 contained a unique gene integrated into the 3'-end of the genome that has no homologs in any known coronavirus, but which sequence and phylogeny analyses indicated most likely originated from the p10 gene of a bat orthoreovirus. Subgenomic mRNA and cellular-level observations demonstrated that the p10 gene is functional and induces the formation of cell syncytia. Therefore, here we report a putative heterologous inter-family recombination event between a single-stranded, positive-sense RNA virus and a double-stranded segmented RNA virus, providing insights into the fundamental mechanisms of viral evolution.

Infectious haematopoietic necrosis virus (IHNV) and viral haemorrhagic septicaemia virus (VHSV) are the most important viral pathogens impacting rainbow trout farming. These viruses are persistent in Italy, where they are responsible for severe disease outbreaks (epizootics) that affect the profitability of the trout industry. Despite the importance of IHNV and VHSV, little is known about their evolution at a local scale, although this is likely to be important for virus eradication and control. To address this issue we performed a detailed molecular evolutionary and epidemiological analysis of IHNV and VHSV in trout farms from northern Italy. Full-length glycoprotein gene sequences of a selection of VHSV (n=108) and IHNV (n=89) strains were obtained. This revealed that Italian VHSV strains belong to sublineages Ia1 and Ia2 of genotype Ia and are distributed into 7 genetic clusters. In contrast, all Italian IHNV isolates fell within genogroup E, for which only a single genetic cluster was identified. More striking was that IHNV has evolved more rapidly than VHSV (mean rates of 11 and 7.3 × 10-4 nucleotide substitutions per site, per year, respectively), indicating that these viruses exhibit fundamentally different evolutionary dynamics. The time to the most recent common ancestor of both IHNV and VHSV was consistent with the first reports of these pathogens in Italy. By combining sequence data with epidemiological information it was possible to identify different patterns of virus spread among trout farms, in which adjacent facilities can be infected by either genetically similar or different viruses, and farms located in different water catchments can be infected by identical strains. Overall, these findings highlight the importance of combining molecular and epidemiological information to identify the determinants of IHN and VHS spread, and to provide data that is central to future surveillance strategies and possibly control.

The development of novel approaches that combine epidemiological and genomic data provides new opportunities to reveal the spatiotemporal dynamics of infectious diseases and determine the processes responsible for their spread and maintenance. Taking advantage of detailed epidemiological time series and viral sequence data from more than 20 years reported by the National Reference Centre for Rabies of Bangui, the capital city of Central African Republic, we used a combination of mathematical modeling and phylogenetic analysis to determine the spatiotemporal dynamics of rabies in domestic dogs as well as the frequency of extinction and introduction events in an African city. We show that although dog rabies virus (RABV) appears to be endemic in Bangui, its epidemiology is in fact shaped by the regular extinction of local chains of transmission coupled with the introduction of new lineages, generating successive waves of spread. Notably, the effective reproduction number during each wave was rarely above the critical value of 1, such that rabies is not self-sustaining in Bangui. In turn, this suggests that rabies at local geographic scales is driven by human-mediated dispersal of RABV among sparsely connected peri-urban and rural areas as opposed to dispersion in a relatively large homogenous urban dog population. This combined epidemiological and genomic approach enables development of a comprehensive framework for understanding disease persistence and informing control measures, indicating that control measures are probably best targeted towards areas neighbouring the city that appear as the source of frequent incursions seeding outbreaks in Bangui.

The 14th–18th century pandemic of Yersinia pestis caused devastating disease outbreaks in Europe for almost 400 years. The reasons for plague’s persistence and abrupt disappearance in Europe are poorly understood, but could have been due to either the presence of now-extinct plague foci in Europe itself, or successive disease introductions from other locations. Here we present five Y. pestis genomes from one of the last European outbreaks of plague, from 1722 in Marseille, France. The lineage identified has not been found in any extant Y. pestis foci sampled to date, and has its ancestry in strains obtained from victims of the 14th century Black Death. These data suggest the existence of a previously uncharacterized historical plague focus that persisted for at least three centuries. We propose that this disease source may have been responsible for the many resurgences of plague in Europe following the Black Death.

The winter seasonality of influenza A virus in temperate climates is one of the most widely recognized, yet least understood, epidemiological patterns in infectious disease. Central to understanding what drives the seasonal emergence of this important human pathogen is determining what becomes of the virus during the non-epidemic summer months. Herein, we take a step towards elucidating the seasonal emergence of influenza virus by determining the evolutionary relationship between populations of influenza A virus sampled from opposite hemispheres. We conducted a phylogenetic analysis of 487 complete genomes of human influenza A/H3N2 viruses collected between 1999 and 2005 from Australia and New Zealand in the southern hemisphere, and a representative sub-sample of viral genome sequences from 413 isolates collected in New York state, United States, representing the northern hemisphere. We show that even in areas as relatively geographically isolated as New Zealand's South Island and Western Australia, global viral migration contributes significantly to the seasonal emergence of influenza A epidemics, and that this migration has no clear directional pattern. These observations run counter to suggestions that local epidemics are triggered by the climate-driven reactivation of influenza viruses that remain latent within hosts between seasons or transmit at low efficiency between seasons. However, a complete understanding of the seasonal movements of influenza A virus will require greatly expanded global surveillance, particularly of tropical regions where the virus circulates year-round, and during non-epidemic periods in temperate climate areas.

The factors that determine the characteristic seasonality of influenza remain enigmatic. Current models predict that occurrences of influenza outside the normal surveillance season within a temperate region largely reflect the importation of viruses from the alternate hemisphere or from equatorial regions in Asia. To help reveal the drivers of seasonality we investigated the origins and evolution of influenza viruses sampled during inter-seasonal periods in Australia. To this end we conducted an expansive phylogenetic analysis of 9912, 3804, and 3941 hemagglutinnin (HA) sequences from influenza A/H1N1pdm, A/H3N2, and B, respectively, collected globally during the period 2009-2014. Of the 1475 viruses sampled from Australia, 396 (26.8% of Australian, or 2.2% of global set) were sampled outside the monitored temperate influenza surveillance season (1 May - 31 October). Notably, rather than simply reflecting short-lived importations of virus from global localities with higher influenza prevalence, we documented a variety of more complex inter-seasonal transmission patterns including "stragglers" from the preceding season and "heralds" of the forthcoming season, and which included viruses sampled from clearly temperate regions within Australia. We also provide evidence for the persistence of influenza B virus between epidemic seasons, in which transmission of a viral lineage begins in one season and continues throughout the inter-seasonal period into the following season. Strikingly, a disproportionately high number of inter-seasonal influenza transmission events occurred in tropical and subtropical regions of Australia, providing further evidence that climate plays an important role in shaping patterns of influenza seasonality.

RNA viruses exhibit substantial structural, ecological and genomic diversity. However, genome size in RNA viruses is likely limited by a high mutation rate, resulting in the evolution of various mechanisms to increase complexity while minimising genome expansion. Here we conduct a large-scale analysis of the genome sequences of 99 animal rhabdoviruses, including 45 genomes which we determined de novo, to identify patterns of genome expansion and the evolution of genome complexity. All but seven of the rhabdoviruses clustered into 17 well-supported monophyletic groups, of which eight corresponded to established genera, seven were assigned as new genera, and two were taxonomically ambiguous. We show that the acquisition and loss of new genes appears to have been a central theme of rhabdovirus evolution, and has been associated with the appearance of alternative, overlapping and consecutive ORFs within the major structural protein genes, and the insertion and loss of additional ORFs in each gene junction in a clade-specific manner. Changes in the lengths of gene junctions accounted for as much as 48.5% of the variation in genome size from the smallest to the largest genome, and the frequency with which new ORFs were observed increased in the 3' to 5' direction along the genome. We also identify several new families of accessory genes encoded in these regions, and show that non-canonical expression strategies involving TURBS-like termination-reinitiation, ribosomal frame-shifts and leaky ribosomal scanning appear to be common. We conclude that rhabdoviruses have an unusual capacity for genomic plasticity that may be linked to their discontinuous transcription strategy from the negative-sense single-stranded RNA genome, and propose a model that accounts for the regular occurrence of genome expansion and contraction throughout the evolution of the Rhabdoviridae.

Lymphoproliferative disease virus (LPDV) is a poorly understood, oncogenic avian retrovirus of domestic turkeys that has historically been restricted to Europe and Israel. However, a recent study reported LPDV in multiple wild turkey diagnostic cases from throughout the eastern United States of America (USA). To better understand the distribution of LPDV in the eastern USA, we surveyed 1,164 reportedly asymptomatic hunter-harvested wild turkeys from 17 states for the presence of LPDV proviral DNA by PCR. In total, 564/1,164 (47%) turkeys were positive for LPDV. Wild turkeys from each state had a relatively high prevalence of LPDV, although statewide prevalence varied from 26 to 83%. Phylogenetic analysis revealed two major clades of LPDV in the USA, although one was at a low frequency suggesting restricted transmission, as well as significant clustering by state of isolation. To determine the best tissue to target for diagnostic purposes, liver, spleen, and bone marrow were tested from a subset of 15 hunter-harvested wild turkeys and 20 wild turkey diagnostic cases. Overall, bone marrow provided the highest level of detection for both hunter-harvested turkeys and diagnostic cases. The sensitivity of LPDV detection between tissues was not significantly different for diagnostic cases, but was for hunter-harvested birds. These results indicate that LPDV infection is common and widespread in wild turkey populations throughout the eastern USA, even without overt signs of disease.

Dengue virus (DENV) transmission is ubiquitous throughout the tropics. More than 70% of the current global dengue disease burden is borne by people who live in the Asia-Pacific region. We sequenced the E gene of DENV isolated from travellers entering Western Australia between 2010-2012, most of whom visited Indonesia, and identified a diverse array of DENV1-4, including multiple co-circulating viral lineages. Most viruses were closely related to lineages known to have circulated in Indonesia for some time, indicating that this geographic region serves as a major hub for dengue genetic diversity. Most notably, we identified a new lineage of DENV-2 (Cosmopolitan genotype) that emerged in Bali in 2011-2012. The spread of this lineage should clearly be monitored. Surveillance of symptomatic returned travellers provides important and timely information on circulating DENV serotypes and genotypes, and can reveal the herald wave of dengue and other emerging infectious diseases.

A complex interplay of viral, host, and ecological factors shapes the spatio-temporal incidence and evolution of human influenza viruses. Although considerable attention has been paid to influenza A viruses, a lack of equivalent data means that an integrated evolutionary and epidemiological framework has until now not been available for influenza B viruses, despite their significant disease burden. Through the analysis of over 900 full genomes from an epidemiological collection of more than 26,000 strains from Australia and New Zealand, we reveal fundamental differences in the phylodynamics of the two co-circulating lineages of influenza B virus (Victoria and Yamagata), showing that their individual dynamics are determined by a complex relationship between virus transmission, age of infection, and receptor binding preference. In sum, this work identifies new factors that are important determinants of influenza B evolution and epidemiology.

Although arthropods are important viral vectors, the biodiversity of arthropod viruses, as well as the role that arthropods have played in viral origins and evolution, is unclear. Through RNA sequencing of 70 arthropod species we discovered 112 novel viruses that appear to be ancestral to much of the documented genetic diversity of negative-sense RNA viruses, a number of which are also present as endogenous genomic copies. With this greatly enriched diversity we revealed that arthropods contain viruses that fall basal to major virus groups, including the vertebrate-specific arenaviruses, filoviruses, hantaviruses, influenza viruses, lyssaviruses, and paramyxoviruses. We similarly documented a remarkable diversity of genome structures in arthropod viruses, including a putative circular form, that sheds new light on the evolution of genome organization. Hence, arthropods are a major reservoir of viral genetic diversity and have likely been central to viral evolution.

Canine parvovirus (CPV) emerged as a new pandemic pathogen of dogs in the 1970s and is closely related to feline panleukopenia virus (FPV), a parvovirus of cats and related carnivores. Although both viruses have wide host ranges, analysis of viral sequences recovered from different wild carnivore species, as shown here, demonstrated that>95% were derived from CPV-like viruses, suggesting that CPV is dominant in sylvatic cycles. Many viral sequences showed host-specific mutations in their capsid proteins, which were often close to sites known to control binding to the transferrin receptor (TfR), the host receptor for these carnivore parvoviruses, and which exhibited frequent parallel evolution. To further examine the process of host adaptation, we passaged parvoviruses with alternative backgrounds in cells from different carnivore hosts. Specific mutations were selected in several viruses and these differed depending on both the background of the virus and the host cells in which they were passaged. Strikingly, these in vitro mutations recapitulated many specific changes seen in viruses from natural populations, strongly suggesting they are host adaptive, and which were shown to result in fitness advantages over their parental virus. Comparison of the sequences of the transferrin receptors of the different carnivore species demonstrated that many mutations occurred in and around the apical domain where the virus binds, indicating that viral variants were likely selected through their fit to receptor structures. Some of the viruses accumulated high levels of variation upon passage in alternative hosts, while others could infect multiple different hosts with no or only a few additional mutations. Overall, these studies demonstrate that the evolutionary history of a virus, including how long it has been circulating and in which hosts, as well as its phylogenetic background, has a profound effect on determining viral host range.

Arthropod-borne viruses are a major cause of emerging disease with significant public health and economic impacts. However, the factors that determine their activity and seasonality are not well understood. In Australia, a network of sentinel cattle herds is used to monitor the distribution of several such viruses and to define virus-free regions. Herein, we utilize these serological data to describe the seasonality, and its drivers, of three economically important animal arboviruses: bluetongue virus, Akabane virus and bovine ephemeral fever virus. Through epidemiological time-series analyses of sero-surveillance data of 180 sentinel herds between 2004-2012, we compared seasonal parameters across latitudes, ranging from the tropical north (-10°S) to the more temperate south (-40°S). This analysis revealed marked differences in seasonality between distinct geographic regions and climates: seasonality was most pronounced in southern regions and gradually decreased as latitude decreased toward the Equator. Further, we show that both the timing of epidemics and the average number of seroconversions have a strong geographical component, which likely reflect patterns of vector abundance through co-varying climatic factors, especially temperature and rainfall. Notably, despite their differences in biology, including insect vector species, all three viruses exhibited very similar seasonality. By revealing the factors that shape spatial and temporal distributions, our study provides a more complete understanding of arbovirus seasonality that will enable better risk predictions.

Host-range shifts in influenza virus are a major risk factor for pandemics. A key question in the study of emerging zoonoses is how the evolution of transmission efficiency interacts with heterogeneity in contact patterns in the new host species, as this interplay influences disease dynamics and prospects for control. Here we use a synergistic mixture of models and data to tease apart the evolutionary and demographic processes controlling a host-range shift in equine H3N8-derived canine influenza virus (CIV). CIV has experienced 15 years of continuous transfer among dogs in the United States, but maintains a patchy distribution, characterized by sporadic short-lived outbreaks coupled with endemic hotspots in large animal shelters. We show that CIV has a high reproductive potential in these facilities (mean R(0) = 3.9) and that these hotspots act as refugia from the sparsely connected majority of the dog population. Intriguingly, CIV has evolved a transmission efficiency that closely matches the minimum required to persist in these refugia, leaving it poised on the extinction/invasion threshold of the host contact network. Corresponding phylogenetic analyses show strong geographic clustering in three US regions, and that the effective reproductive number of the virus (R(e)) in the general dog population is close to 1.0. Our results highlight the critical role of host contact structure in CIV dynamics, and show how host contact networks could shape the evolution of pathogen transmission efficiency. Importantly, efficient control measures could eradicate the virus, in turn minimizing the risk of future sustained transmission among companion dogs that could represent a potential new axis to the human-animal interface for influenza.

Hemorrhagic fevers (HF) caused by viruses and bacteria are a major public health problem in China and characterized by variable clinical manifestations, such that it is often difficult to achieve accurate diagnosis and treatment. The causes of HF in 85 patients admitted to Dandong hospital, China, between 2011-2012 were determined by serological and PCR tests. Of these, 34 patients were diagnosed with Huaiyangshan hemorrhagic fever (HYSHF), 34 with Hemorrhagic Fever with Renal Syndrome (HFRS), one with murine typhus, and one with scrub typhus. Etiologic agents could not be determined in the 15 remaining patients. Phylogenetic analyses of recovered bacterial and viral sequences revealed that the causative infectious agents were closely related to those described in other geographical regions. As these diseases have no distinctive clinical features in their early stage, only 13 patients were initially accurately diagnosed. The distinctive clinical features of HFRS and HYSHF developed during disease progression. Enlarged lymph nodes, cough, sputum, and diarrhea were more common in HYSHF patients, while more HFRS cases presented with headache, sore throat, oliguria, percussion pain kidney area, and petechiae. Additionally, HYSHF patients displayed significantly lower levels of white blood cells (WBC), higher levels of creations kinase (CK) and alanine aminotransferase (ALT), while HFRS patients presented with an elevation of blood urea nitrogen (BUN) and creatinine (CREA). These clinical features will assist in the accurate diagnosis of both HYSHF and HFRS. Overall, our data reveal the complexity of pathogens causing HFs in a single Chinese hospital, and highlight the need for accurate early diagnosis and a better understanding of their distinctive clinical features.

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.

Hantaviruses are among the most important zoonotic pathogens of humans and the subject of heightened global attention. Despite the importance of hantaviruses for public health, there is no consensus on their evolutionary history and especially the frequency of virus-host co-divergence versus cross-species virus transmission. Documenting the extent of hantavirus biodiversity, and particularly their range of mammalian hosts, is critical to resolving this issue. Here, we describe four novel hantaviruses (Huangpi virus, Lianghe virus, Longquan virus, and Yakeshi virus) sampled from bats and shrews in China, and which are distinct from other known hantaviruses. Huangpi virus was found in Pipistrellus abramus, Lianghe virus in Anourosorex squamipes, Longquan virus in Rhinolophus affinis, Rhinolophus sinicus, and Rhinolophus monoceros, and Yakeshi virus in Sorex isodon, respectively. A phylogenetic analysis of the available diversity of hantaviruses reveals the existence of four phylogroups that infect a range of mammalian hosts, as well as the occurrence of ancient reassortment events between the phylogroups. Notably, the phylogenetic histories of the viruses are not always congruent with those of their hosts, suggesting that cross-species transmission has played a major role during hantavirus evolution and at all taxonomic levels, although we also noted some evidence for virus-host co-divergence. Our phylogenetic analysis also suggests that hantaviruses might have first appeared in Chiroptera (bats) or Soricomorpha (moles and shrews), before emerging in rodent species. Overall, these data indicate that bats are likely to be important natural reservoir hosts of hantaviruses.

Large-scale codon re-encoding represents a powerful method of attenuating viruses to generate safe and cost-effective vaccines. In contrast to specific approaches of codon re-encoding which modify genome-scale properties, we evaluated the effects of random codon re-encoding on the re-emerging human pathogen Chikungunya virus (CHIKV), and assessed the stability of the resultant viruses during serial in cellulo passage. Using different combinations of three 1.4 kb randomly re-encoded regions located throughout the CHIKV genome six codon re-encoded viruses were obtained. Introducing a large number of slightly deleterious synonymous mutations reduced the replicative fitness of CHIKV in both primate and arthropod cells, demonstrating the impact of synonymous mutations on fitness. Decrease of replicative fitness correlated with the extent of re-encoding, an observation that may assist in the modulation of viral attenuation. The wild-type and two re-encoded viruses were passaged 50 times either in primate or insect cells, or in each cell line alternately. These viruses were analyzed using detailed fitness assays, complete genome sequences and the analysis of intra-population genetic diversity. The response to codon re-encoding and adaptation to culture conditions occurred simultaneously, resulting in significant replicative fitness increases for both re-encoded and wild type viruses. Importantly, however, the most re-encoded virus failed to recover its replicative fitness. Evolution of these viruses in response to codon re-encoding was largely characterized by the emergence of both synonymous and non-synonymous mutations, sometimes located in genomic regions other than those involving re-encoding, and multiple convergent and compensatory mutations. However, there was a striking absence of codon reversion (

Intra-host sequence data from RNA viruses have revealed the ubiquity of defective viruses in natural viral populations, sometimes at surprisingly high frequency. Although defective viruses have long been known to laboratory virologists, their relevance in clinical and epidemiological settings has not been established. The discovery of long-term transmission of a defective lineage of dengue virus type 1 (DENV-1) in Myanmar, first seen in 2001, raised important questions about the emergence of transmissible defective viruses and their role in viral epidemiology. By combining phylogenetic analyses and dynamical modeling, we investigate how evolutionary and ecological processes at the intra-host and inter-host scales shaped the emergence and spread of the defective DENV-1 lineage. We show that this lineage of defective viruses emerged between June 1998 and February 2001, and that the defective virus was transmitted primarily through co-transmission with the functional virus to uninfected individuals. We provide evidence that, surprisingly, this co-transmission route has a higher transmission potential than transmission of functional dengue viruses alone. Consequently, we predict that the defective lineage should increase overall incidence of dengue infection, which could account for the historically high dengue incidence reported in Myanmar in 2001-2002. Our results show the unappreciated potential for defective viruses to impact the epidemiology of human pathogens, possibly by modifying the virulence-transmissibility trade-off, or to emerge as circulating infections in their own right. They also demonstrate that interactions between viral variants, such as complementation, can open new pathways to viral emergence.

Aquatic birds harbor diverse influenza A viruses and are a major viral reservoir in nature. The recent discovery of influenza viruses of a new H17N10 subtype in Central American fruit bats suggests that other New World species may similarly carry divergent influenza viruses. Using consensus degenerate RT-PCR, we identified a novel influenza A virus, designated as H18N11, in a flat-faced fruit bat (Artibeus planirostris) from Peru. Serologic studies with the recombinant H18 protein indicated that several Peruvian bat species were infected by this virus. Phylogenetic analyses demonstrate that, in some gene segments, New World bats harbor more influenza virus genetic diversity than all other mammalian and avian species combined, indicative of a long-standing host-virus association. Structural and functional analyses of the hemagglutinin and neuraminidase indicate that sialic acid is not a ligand for virus attachment nor a substrate for release, suggesting a unique mode of influenza A virus attachment and activation of membrane fusion for entry into host cells. Taken together, these findings indicate that bats constitute a potentially important and likely ancient reservoir for a diverse pool of influenza viruses.

BACKGROUND: Longquan City, Zhejiang province, China, has been seriously affected by hemorrhagic fever with renal syndrome (HFRS) since the first cases were registered in 1974. To understand the epidemiology and emergence of HFRS in Longquan, which may be indicative of large parts of rural China, we studied long-term incidence patterns and performed a molecular epidemiological investigation of the causative hantaviruses in human and rodent populations. METHOD/PRINCIPAL FINDINGS: During 1974-2011, 1866 cases of HFRS were recorded in Longquan, including 20 deaths. In 2011, the incidence of HFRS remained high, with 19.61 cases/100,000 population, despite the onset of vaccination in 1997. During 1974-1998, HFRS cases in Longquan occurred mainly in winter, while in the past decade the peak of HFRS has shifted to the spring. Notably, the concurrent prevalence of rodent-borne hantaviruses in the region was also high. Phylogenetic analyses of viral sequences recovered from rodents in Longquan revealed the presence of novel genetic variants of Gou virus (GOUV) in Rattus sp. rats and Hantaan virus (HTNV) in the stripe field mice, respectively. Strikingly, viral sequences sampled from infected humans were very closely related to those from rodents. CONCLUSIONS/SIGNIFICANCE: HFRS represents an important public health problem in Longquan even after years of preventive measures. Our data suggest that continual spillover of the novel genetic variant of GOUV and the new genetic lineage of HTNV are responsible for the high prevalence of HFRS in humans. In addition, this is the first report of GOUV associated with human HFRS cases, and our data suggest that GOUV is now the major cause of HFRS in this region.

Revealing the patterns and determinants of the spread of dengue virus (DENV) at local scales is central to understanding the epidemiology and evolution of this major human pathogen. We performed a phylogenetic analysis of the envelope (E) genes of DENV-1, -2, -3, and -4 isolates (involving 97, 23, 5, and 74 newly collected sequences, respectively) sampled from school-based cohort and village-based cluster studies in Kamphaeng Phet, Thailand, between 2004 and 2007. With these data, we sought to describe the spatial and temporal patterns of DENV spread within a rural population where a future vaccine efficacy trial is planned. Our analysis revealed considerable genetic diversity within the study population, with multiple lineages within each serotype circulating for various lengths of time during the study period. These results suggest that DENV is frequently introduced into both semi-urban and rural areas in Kamphaeng Phet from other populations. In contrast, the persistence of viral lineages across sampling years was observed less frequently. Analysis of phylogenetic clustering indicated that DENV transmission was highly spatially and temporally focal, and that it occurred in homes rather than at school. Overall, the strength of temporal clustering suggests that seasonal bottlenecks in local DENV populations facilitate the invasion and establishment of viruses from outside of the study area, in turn reducing the extent of lineage persistence.

Dengue virus transmission occurs in both epidemic and endemic cycles across tropical and sub-tropical regions of the world. Incidence is particularly high in much of Southeast Asia, where hyperendemic transmission plagues both urban and rural populations. However, endemicity has not been established in some areas with climates that may not support year-round viral transmission. An understanding of how dengue viruses (DENV) enter these environments and whether the viruses persist in inapparent local transmission cycles is central to understanding how dengue emerges in areas at the margins of endemic transmission. Dengue is highly endemic in tropical southern Vietnam, while increasingly large seasonal epidemics have occurred in northern Viet Nam over the last decade. We have investigated the spread of DENV-1 throughout Vietnam to determine the routes by which the virus enters northern and central regions of the country. Phylogeographic analysis of 1,765 envelope (E) gene sequences from Southeast Asia revealed frequent movement of DENV between neighboring human populations and strong local clustering of viral lineages. Long-distance migration of DENV between human population centers also occurred regularly and on short time-scales, indicating human-mediated viral invasion into northern Vietnam. Human populations in southern Vietnam were found to be the primary source of DENV circulating throughout the country, while central and northern Vietnam acted as sink populations, likely due to reduced connectedness to other populations in the case of the central regions and to the influence of temperature variability on DENV replication and vector survival and competence in the north. Finally, phylogeographic analyses suggested that viral movement follows a gravity model and indicates that population immunity and physical and economic connections between populations may play important roles in shaping patterns of DENV transmission.

The ability of influenza A viruses (IAVs) to cross species barriers and evade host immunity is a major public health concern. Studies on the phylodynamics of IAVs across different scales - from the individual to the population - are essential for devising effective measures to predict, prevent or contain influenza emergence. Understanding how IAVs spread and evolve during outbreaks is critical for the management of epidemics. Reconstructing the transmission network during a single outbreak by sampling viral genetic data in time and space can generate insights about these processes. Here, we obtained intra-host viral sequence data from horses infected with equine influenza virus (EIV) to reconstruct the spread of EIV during a large outbreak. To this end, we analyzed within-host viral populations from sequences covering 90% of the infected yards. By combining gene sequence analyses with epidemiological data, we inferred a plausible transmission network, in turn enabling the comparison of transmission patterns during the course of the outbreak and revealing important epidemiological features that were not apparent using either approach alone. The EIV populations displayed high levels of genetic diversity, and in many cases we observed distinct viral populations containing a dominant variant and a number of related minor variants that were transmitted between infectious horses. In addition, we found evidence of frequent mixed infections and loose transmission bottlenecks in these naturally occurring populations. These frequent mixed infections likely influence the size of epidemics.

Both multiple sequence alignment and phylogenetic analysis are problematic in the "twilight zone" of sequence similarity (≤ 25% amino acid identity). Herein we explore the accuracy of phylogenetic inference at extreme sequence divergence using a variety of simulated data sets. We evaluate four leading multiple sequence alignment (MSA) methods (MAFFT, T-COFFEE, CLUSTAL, and MUSCLE) and six commonly used programs of tree estimation (Distance-based: Neighbor-Joining; Character-based: PhyML, RAxML, GARLI, Maximum Parsimony, and Bayesian) against a novel MSA-independent method (PHYRN) described here. Strikingly, at "midnight zone" genetic distances (~7% pairwise identity and 4.0 gaps per position), PHYRN returns high-resolution phylogenies that outperform traditional approaches. We reason this is due to PHRYN's capability to amplify informative positions, even at the most extreme levels of sequence divergence. We also assess the applicability of the PHYRN algorithm for inferring deep evolutionary relationships in the divergent DANGER protein superfamily, for which PHYRN infers a more robust tree compared to MSA-based approaches. Taken together, these results demonstrate that PHYRN represents a powerful mechanism for mapping uncharted frontiers in highly divergent protein sequence data sets.

The attenuation of myxoma virus (MYXV) following its introduction as a biological control into the European rabbit populations of Australia and Europe is the canonical study of the evolution of virulence. However, the evolutionary genetics of this profound change in host-pathogen relationship is unknown. We describe the genome-scale evolution of MYXV covering a range of virulence grades sampled over 49 years from the parallel Australian and European epidemics, including the high-virulence progenitor strains released in the early 1950s. MYXV evolved rapidly over the sampling period, exhibiting one of the highest nucleotide substitution rates ever reported for a double-stranded DNA virus, and indicative of a relatively high mutation rate and/or a continually changing selective environment. Our comparative sequence data reveal that changes in virulence involved multiple genes, likely losses of gene function due to insertion-deletion events, and no mutations common to specific virulence grades. Hence, despite the similarity in selection pressures there are multiple genetic routes to attain either highly virulent or attenuated phenotypes in MYXV, resulting in convergence for phenotype but not genotype.

Influenza viruses are characterized by an ability to cross species boundaries and evade host immunity, sometimes with devastating consequences. The 2009 pandemic of H1N1 influenza A virus highlights the importance of pigs in influenza emergence, particularly as intermediate hosts by which avian viruses adapt to mammals before emerging in humans. Although segment reassortment has commonly been associated with influenza emergence, an expanded host-range is also likely to be associated with the accumulation of specific beneficial point mutations. To better understand the mechanisms that shape the genetic diversity of avian-like viruses in pigs, we studied the evolutionary dynamics of an Eurasian Avian-like swine influenza virus (EA-SIV) in naïve and vaccinated pigs linked by natural transmission. We analyzed multiple clones of the hemagglutinin 1 (HA1) gene derived from consecutive daily viral populations. Strikingly, we observed both transient and fixed changes in the consensus sequence along the transmission chain. Hence, the mutational spectrum of intra-host EA-SIV populations is highly dynamic and allele fixation can occur with extreme rapidity. In addition, mutations that could potentially alter host-range and antigenicity were transmitted between animals and mixed infections were commonplace, even in vaccinated pigs. Finally, we repeatedly detected distinct stop codons in virus samples from co-housed pigs, suggesting that they persisted within hosts and were transmitted among them. This implies that mutations that reduce viral fitness in one host, but which could lead to fitness benefits in a novel host, can circulate at low frequencies.

Dengue is one of the most important infectious diseases of humans and has spread throughout much of the tropical and subtropical world. Despite this widespread dispersal, the determinants of dengue transmission in endemic populations are not well understood, although essential for virus control. To address this issue we performed a phylogeographic analysis of 751 complete genome sequences of dengue 1 virus (DENV-1) sampled from both rural (Dong Thap) and urban (Ho Chi Minh City) populations in southern Viet Nam during the period 2003-2008. We show that DENV-1 in Viet Nam exhibits strong spatial clustering, with likely importation from Cambodia on multiple occasions. Notably, multiple lineages of DENV-1 co-circulated in Ho Chi Minh City. That these lineages emerged at approximately the same time and dispersed over similar spatial regions suggests that they are of broadly equivalent fitness. We also observed an important relationship between the density of the human host population and the dispersion rate of dengue, such that DENV-1 tends to move from urban to rural populations, and that densely populated regions within Ho Chi Minh City act as major transmission foci. Despite these fluid dynamics, the dispersion rates of DENV-1 are relatively low, particularly in Ho Chi Minh City where the virus moves less than an average of 20 km/year. These low rates suggest a major role for mosquito-mediated dispersal, such that DENV-1 does not need to move great distances to infect a new host when there are abundant susceptibles, and imply that control measures should be directed toward the most densely populated urban environments.

The emergence and rapid global spread of the swine-origin H1N1/09 pandemic influenza A virus in humans underscores the importance of swine populations as reservoirs for genetically diverse influenza viruses with the potential to infect humans. However, despite their significance for animal and human health, relatively little is known about the phylogeography of swine influenza viruses in the United States. This study utilizes an expansive data set of hemagglutinin (HA1) sequences (n = 1516) from swine influenza viruses collected in North America during the period 2003-2010. With these data we investigate the spatial dissemination of a novel influenza virus of the H1 subtype that was introduced into the North American swine population via two separate human-to-swine transmission events around 2003. Bayesian phylogeographic analysis reveals that the spatial dissemination of this influenza virus in the US swine population follows long-distance swine movements from the Southern US to the Midwest, a corn-rich commercial center that imports millions of swine annually. Hence, multiple genetically diverse influenza viruses are introduced and co-circulate in the Midwest, providing the opportunity for genomic reassortment. Overall, the Midwest serves primarily as an ecological sink for swine influenza in the US, with sources of virus genetic diversity instead located in the Southeast (mainly North Carolina) and South-central (mainly Oklahoma) regions. Understanding the importance of long-distance pig transportation in the evolution and spatial dissemination of the influenza virus in swine may inform future strategies for the surveillance and control of influenza, and perhaps other swine pathogens.

Although overall pollinator populations have declined over the last couple of decades, the honey bee (Apis mellifera) malady, colony collapse disorder (CCD), has caused major concern in the agricultural community. Among honey bee pathogens, RNA viruses are emerging as a serious threat and are suspected as major contributors to CCD. Recent detection of these viral species in bumble bees suggests a possible wider environmental spread of these viruses with potential broader impact. It is therefore vital to study the ecology and epidemiology of these viruses in the hymenopteran pollinator community as a whole. We studied the viral distribution in honey bees, in their pollen loads, and in other non-Apis hymenopteran pollinators collected from flowering plants in Pennsylvania, New York, and Illinois in the United States. Viruses in the samples were detected using reverse transcriptase-PCR and confirmed by sequencing. For the first time, we report the molecular detection of picorna-like RNA viruses (deformed wing virus, sacbrood virus and black queen cell virus) in pollen pellets collected directly from forager bees. Pollen pellets from several uninfected forager bees were detected with virus, indicating that pollen itself may harbor viruses. The viruses in the pollen and honey stored in the hive were demonstrated to be infective, with the queen becoming infected and laying infected eggs after these virus-contaminated foods were given to virus-free colonies. These viruses were detected in eleven other non-Apis hymenopteran species, ranging from many solitary bees to bumble bees and wasps. This finding further expands the viral host range and implies a possible deeper impact on the health of our ecosystem. Phylogenetic analyses support that these viruses are disseminating freely among the pollinators via the flower pollen itself. Notably, in cases where honey bee apiaries affected by CCD harbored honey bees with Israeli Acute Paralysis virus (IAPV), nearby non-Apis hymenopteran pollinators also had IAPV, while those near apiaries without IAPV did not. In containment greenhouse experiments, IAPV moved from infected honey bees to bumble bees and from infected bumble bees to honey bees within a week, demonstrating that the viruses could be transmitted from one species to another. This study adds to our present understanding of virus epidemiology and may help explain bee disease patterns and pollinator population decline in general.

Understanding the role of humans in the dispersal of predominantly animal pathogens is essential for their control. We used newly developed Bayesian phylogeographic methods to unravel the dynamics and determinants of the spread of dog rabies virus (RABV) in North Africa. Each of the countries studied exhibited largely disconnected spatial dynamics with major geopolitical boundaries acting as barriers to gene flow. Road distances proved to be better predictors of the movement of dog RABV than accessibility or raw geographical distance, with occasional long distance and rapid spread within each of these countries. Using simulations that bridge phylodynamics and spatial epidemiology, we demonstrate that the contemporary viral distribution extends beyond that expected for RABV transmission in African dog populations. These results are strongly supportive of human-mediated dispersal, and demonstrate how an integrated phylogeographic approach will turn viral genetic data into a powerful asset for characterizing, predicting, and potentially controlling the spatial spread of pathogens.

A better description of the extent and structure of genetic diversity in dengue virus (DENV) in endemic settings is central to its eventual control. To this end we determined the complete coding region sequence of 187 DENV-2 genomes and 68 E genes from viruses sampled from Vietnamese patients between 1995 and 2009. Strikingly, an episode of genotype replacement was observed, with Asian 1 lineage viruses entirely displacing the previously dominant Asian/American lineage viruses. This genotype replacement event also seems to have occurred within DENV-2 in Thailand and Cambodia, suggestive of a major difference in viral fitness. To determine the cause of this major evolutionary event we compared both the infectivity of the Asian 1 and Asian/American genotypes in mosquitoes and their viraemia levels in humans. Although there was little difference in infectivity in mosquitoes, we observed significantly higher plasma viraemia levels in paediatric patients infected with Asian 1 lineage viruses relative to Asian/American viruses, a phenotype that is predicted to result in a higher probability of human-to-mosquito transmission. These results provide a mechanistic basis to a marked change in the genetic structure of DENV-2 and more broadly underscore that an understanding of DENV evolutionary dynamics can inform the development of vaccines and anti-viral drugs.