Search

Ticks harbour a high diversity of viruses, bacteria and protozoa. The soft tick Carios vespertilionis (Argasidae) is a common ectoparasite of bats in the Palearctic region and is suspected to be vector and reservoir of viruses and other microbial species in bat populations, some of which may act as zoonotic agents for human disease. The Soprano pipistrelle (Pipistrellus pygmaeus, Vespertilionidae) is widely distributed in Europe, where it can be found inside or close to human habitation. We used meta-transcriptomic sequencing to determine the RNA virome and common microbiota in blood-fed C. vespertilionis ticks collected from a Soprano pipistrelle bat roosting site in south-central Sweden. Our analyses identified 16 viruses from 11 virus families, of which 15 viruses were novel. For the first time in Sweden we identified Issuk-Kul virus, a zoonotic arthropod-borne virus previously associated with outbreaks of acute febrile illness in humans. Probable bat-associated and tick-borne viruses were classified within the families Nairoviridae, Caliciviridae and Hepeviridae, while other invertebrate-associated viruses included members of the Dicistroviridae, Iflaviridae, Nodaviridae, Partitiviridae, Permutotetraviridae, Polycipiviridae and Solemoviridae. Similarly, we found abundant bacteria in C. vespertilionis, including genera with known tick-borne bacteria, such as Coxiella spp. and Rickettsia spp. These findings demonstrate the remarkable diversity of RNA viruses and bacteria present in C. vespertilionis and highlight the importance of bat-associated ectoparasite surveillance as an effective and non-invasive means to track viruses and bacteria circulating in bats and ticks.

Encephalitis is most often caused by a variety of infectious agents identified through diagnostic tests utilizing cerebrospinal fluid. We investigated the clinical characteristics and potential aetiological agents of unexplained encephalitis through metagenomic sequencing of residual clinical samples from multiple tissue types and independent clinical review. Forty-three specimens were collected from 18 encephalitis cases with no cause identified by the Australian Childhood Encephalitis study. Samples were subjected to total RNA sequencing (‘metatranscriptomics’) to determine the presence and abundance of potential pathogens, and to describe the possible aetiologies of unexplained encephalitis. Using this protocol, we identified five RNA and two DNA viruses associated with human infection from both non-sterile and sterile sites, which were confirmed by PCR. These comprised two human rhinoviruses, two human seasonal coronaviruses, two polyomaviruses and one picobirnavirus. Human rhinovirus and seasonal coronaviruses may be responsible for five of the encephalitis cases. Immune-mediated encephalitis was considered likely in six cases and metatranscriptomics did not identify a possible pathogen in these cases. The aetiology remained unknown in nine cases. Our study emphasizes the importance of respiratory viruses in the aetiology of unexplained child encephalitis and suggests that non-central-nervous-system sampling in encephalitis clinical guidelines and protocols could improve the diagnostic yield.

The endosymbiont bacteria of the genus Wolbachia are associated with multiple mutualistic effects on insect biology, including nutritional and antiviral properties. Members of the genus Wolbachia naturally occur in fly species of the genus Drosophila, providing an operational model host for studying how virome composition may be affected by its presence. Drosophila simulans populations can carry a variety of strains of members of the genus Wolbachia , with the wAu strain associated with strong antiviral protection under experimental conditions. We used D. simulans sampled from the Perth Hills, Western Australia, to investigate the potential virus protective effect of the wAu strain of Wolbachia on individual wild-caught flies. Our data revealed no appreciable variation in virus composition and abundance between individuals infected or uninfected with Wolbachia associated with the presence or absence of wAu. However, it remains unclear whether wAu might affect viral infection and host survival by increasing tolerance rather than inducing complete resistance. These data also provide new insights into the natural virome diversity of D. simulans. Despite the small number of individuals sampled, we identified a repertoire of RNA viruses, including nora virus, galbut virus, thika virus and La Jolla virus, that have been identified in other species of the genus Drosophila. Chaq virus-like sequences associated with galbut virus were also detected. In addition, we identified five novel viruses from the families Reoviridae, Tombusviridae, Mitoviridae and Bunyaviridae. Overall, this study highlights the complex interaction between Wolbachia and RNA virus infections and provides a baseline description of the natural virome of D. simulans.

The identification of SARS-CoV-2-like viruses in Malayan pangolins (Manis javanica) has focused attention on these endangered animals and the viruses they carry. We successfully isolated a novel respirovirus from the lungs of a dead Malayan pangolin. Similar to murine respirovirus, the full-length genome of this novel virus was 15 384 nucleotides comprising six genes in the order 3′–(leader)–NP-P-M-F-HN-l-(trailer)−5’. Phylogenetic analysis revealed that this virus belongs to the genus Respirovirus and is most closely related to murine respirovirus. Notably, animal infection experiments indicated that the pangolin virus is highly pathogenic and transmissible in mice, with inoculated mice having variable clinical symptoms and a fatality rate of 70.37 %. The virus was found to replicate in most tissues with the exception of muscle and heart. Contact transmission of the virus was 100 % efficient, although the mice in the contact group displayed milder symptoms, with the virus mainly being detected in the trachea and lungs. The isolation of a novel respirovirus from the Malayan pangolin provides new insight into the evolution and distribution of this important group of viruses and again demonstrates the potential infectious disease threats faced by endangered pangolins.

Although RNA viruses are likely to exist in every species of cellular life, our knowledge of their biodiversity and evolution has been focused on those that can cause disease in humans and in economically important animals and plants. Recently published studies, however, have shown that every healthy organism can carry viruses. In this study, we present the unbiased discovery of highly divergent RNA viruses in European rodents and rabbits using Next Generation Sequencing (NGS). Tissue samples from different organs were collected from UK rodents and French rabbits and were initially screened for viruses with PCR. Following up on preliminary data, the positive samples were sent for NGS to acquire full genomes and perform unbiased virus discovery. Our findings encompass a number of novel viruses including astrovirus, rotavirus A, hantavirus, picornavirus, coronavirus, paramyxovirus etc. Among the novel viruses, the picornavirus interestingly showed 50 % and 34 % similarity to the closest relative for the non-structural and the structural protein respectively at amino acid level, suggesting a potentially novel genus within the family. This study shows the presence of highly divergent RNA viruses in European rodents and rabbits enriching thus our current knowledge of their origin and evolution. Finally, this study shows that these animal species can be the reservoirs of RNA viruses notorious for cross-species transmission, increasing the risk of a spill-over to humans or livestock.

Arboretum virus (ABTV) and Puerto Almendras virus (PTAMV) are two mosquito-associated rhabdoviruses isolated from pools of Psorophora albigenu and Ochlerotattus fulvus mosquitoes, respectively, collected in the Department of Loreto, Peru, in 2009. Initial tests suggested that both viruses were novel rhabdoviruses and this was confirmed by complete genome sequencing. Analysis of their 11 482 nt (ABTV) and 11 876 (PTAMV) genomes indicates that they encode the five canonical rhabdovirus structural proteins (N, P, M, G and L) with an additional gene (U1) encoding a small hydrophobic protein. Evolutionary analysis of the L protein indicates that ABTV and PTAMV are novel and phylogenetically distinct rhabdoviruses that cannot be classified as members of any of the eight currently recognized genera within the family Rhabdoviridae, highlighting the vast diversity of this virus family.

Kolente virus (KOLEV) is a rhabdovirus originally isolated from ticks and a bat in Guinea, West Africa, in 1985. Although tests at the time of isolation suggested that KOLEV is a novel rhabdovirus, it has remained largely uncharacterized. We assembled the complete genome sequence of the prototype strain DakAr K7292, which was found to encode the five canonical rhabdovirus structural proteins (N, P, M, G and L) with alternative ORFs (>180 nt) in the P and L genes. Serologically, KOLEV exhibited a weak antigenic relationship with Barur and Fukuoka viruses in the Kern Canyon group. Phylogenetic analysis revealed that KOLEV represents a distinct and divergent lineage that shows no clear relationship to any rhabdovirus except Oita virus, although with limited phylogenetic resolution. In summary, KOLEV represents a novel species in the family Rhabdoviridae.

Polyomaviruses (PyVs) have been identified in a wide range of avian and mammalian species. However, little is known about their occurrence, genetic diversity and evolutionary history in bats, even though bats are important reservoirs for many emerging viral pathogens. This study screened 380 specimens from 35 bat species from Kenya and Guatemala for the presence of PyVs by semi-nested pan-PyV PCR assays. PyV DNA was detected in 24 of the 380 bat specimens. Phylogenetic analysis revealed that the bat PyV sequences formed 12 distinct lineages. Full-genome sequences were obtained for seven representative lineages and possessed similar genomic features to known PyVs. Strikingly, this evolutionary analysis revealed that the bat PyVs were paraphyletic, suggestive of multiple species jumps between bats and other mammalian species, such that the theory of virus–host co-divergence for mammalian PyVs as a whole could be rejected. In addition, evidence was found for strong heterogeneity in evolutionary rate and potential recombination in a number of PyV complete genomes, which complicates both phylogenetic analysis and virus classification. In summary, this study revealed that bats are important reservoirs of PyVs and that these viruses have a complex evolutionary history.

Revealing the frequency and determinants of reassortment among RNA genome segments is fundamental to understanding basic aspects of the biology and evolution of the influenza virus. To estimate the extent of genomic reassortment in influenza viruses circulating in North American swine, we performed a phylogenetic analysis of 139 whole-genome viral sequences sampled during 1998–2011 and representing seven antigenically distinct viral lineages. The highest amounts of reassortment were detected between the H3 and the internal gene segments (PB2, PB1, PA, NP, M and NS), while the lowest reassortment frequencies were observed among the H1γ, H1pdm and neuraminidase segments, particularly N1. Less reassortment was observed among specific haemagglutinin–neuraminidase combinations that were more prevalent in swine, suggesting that some genome constellations may be evolutionarily more stable.

The genus Flavivirus, which contains approximately 70 single-stranded, positive-sense RNA viruses, represents a unique model for studying the evolution of vector-borne disease, as it includes viruses that are mosquito-borne, tick-borne or have no known vector. Both theoretical work and field studies suggest the existence of a large number of undiscovered flaviviruses. Recently, the first isolation of cell fusing agent virus (CFAV) was reported from a natural mosquito population in Puerto Rico, and sequences related to CFAV have been discovered in mosquitoes from Thailand. CFAV had previously been isolated from a mosquito cell line in 1975 and represented the only known ‘insect-only’ flavivirus, appearing to replicate in insect cells alone. A second member of the ‘insect-only’ group, Kamiti River virus (KRV), was isolated from Kenyan mosquitoes in 2003. A third tentative member of the ‘insect-only’ group, Culex flavivirus (CxFV), was first isolated in 2007 from Japan and further strains have subsequently been reported from the Americas. We report the discovery, isolation and characterization of two novel ‘insect-only’ flaviviruses from Entebbe, Uganda: a novel lineage tentatively designated Nakiwogo virus (NAKV) and a new strain of CxFV. The individual mosquitoes from which these strains were isolated, identified retrospectively by using a reference molecular phylogeny generated using voucher specimens from the region, were Mansonia africana nigerrima and Culex quinquefasciatus, respectively. This represents the first isolation, to our knowledge, of a novel insect-only flavivirus from a Mansonia species and the first isolation of a strain of CxFV from Africa.

The burden of rabies in Africa is estimated at 24 000 human deaths year−1, almost all of which result from infection with dog rabies viruses (RABV). To investigate the evolutionary dynamics of RABV in western and central Africa, 92 isolates sampled from 27 African countries over 29 years were collected and sequenced. This revealed that RABV currently circulating in dogs in this region fell into a single lineage designated ‘Africa 2’. A detailed analysis of the phylogeographical structure of this Africa 2 lineage revealed strong population subdivision at the country level, with only limited movement of virus among localities, including a possible east-to-west spread across Africa. In addition, Bayesian coalescent analysis suggested that the Africa 2 lineage was introduced into this region of Africa only recently (probably

Human (HMPV) and avian (AMPV) metapneumoviruses are closely related viruses that cause respiratory tract illnesses in humans and birds, respectively. Although HMPV was first discovered in 2001, retrospective studies have shown that HMPV has been circulating in humans for at least 50 years. AMPV was first isolated in the 1970s, and can be classified into four subgroups, A–D. AMPV subgroup C is more closely related to HMPV than to any other AMPV subgroup, suggesting that HMPV has emerged from AMPV-C upon zoonosis. Presently, at least four genetic lineages of HMPV circulate in human populations – A1, A2, B1 and B2 – of which lineages A and B are antigenically distinct. We used a Bayesian Markov Chain Monte Carlo (MCMC) framework to determine the evolutionary and epidemiological dynamics of HMPV and AMPV-C. The rates of nucleotide substitution, relative genetic diversity and time to the most recent common ancestor (TMRCA) were estimated using large sets of sequences of the nucleoprotein, the fusion protein and attachment protein genes. The sampled genetic diversity of HMPV was found to have arisen within the past 119–133 years, with consistent results across all three genes, while the TMRCA for HMPV and AMPV-C was estimated to have existed around 200 years ago. The relative genetic diversity observed in the four HMPV lineages was low, most likely reflecting continual population bottlenecks, with only limited evidence for positive selection.

A comparative analysis was performed of the 3′-untranslated region (UTR) of Dengue virus (DENV) sampled from Bangkok, Thailand, over a 30 year period and representing all four serotypes. Considerable genetic variation was observed both within and among serotypes. Notably, a full-length version of the critical 3′-long stable hairpin structure was absent from some isolates, suggesting the occurrence of complex structural interactions within the 3′-UTR, including the influence of upstream mutations. The Thai sequences were then combined with 61 globally sampled isolates of DENV taken from patients with either dengue fever or severe dengue disease. No consistent association was found between 3′-UTR secondary structure and the clinical outcome of DENV infection, although some evidence for a trend in this direction was observed in DENV-2. It was concluded that the 3′-UTR is not the sole determinant of DENV virulence in nature, although variation in secondary structure may greatly influence viral fitness.

A maximum-likelihood approach was used to analyse selection pressures acting on genes from all four serotypes of dengue virus (DEN). A number of amino acid positions were identified within the envelope (E) glycoprotein that have been subject to relatively weak positive selection in both DEN-3 and DEN-4, as well as in two of the five genotypes of DEN-2. No positive selection was detected in DEN-1. In accordance with the function of the E protein as the major antigenic determinant of DEN, the majority of these sites were located in, or near to, potential T- or B-cell epitopes. A smaller number of selected sites was located in other well-defined functional domains of the E protein, suggesting that cell tropism and virus-mediated membrane fusion may also confer fitness advantages to DEN in nature. Several positively selected amino acid substitutions were also identified in the NS2B and NS5 genes of DEN-2, although the cause of this selection is unclear, whereas the capsid, membrane and non-structural genes NS1, NS2A, NS3 and NS4 were all subject to strong functional constraints. Hence, evidence was found for localized adaptive evolution in natural isolates of DEN, revealing that selection pressures differ among serotypes, genotypes and viral proteins.

We compared the extent of positive selection acting on acute and persistent strains of measles virus (MV). Far stronger positive selection was found in the fusion (F) and haemagglutinin (H) genes from subacute sclerosing panencephalitis (SSPE) compared to acute MV cases. Most of the positively selected sites identified in these surface glycoprotein genes from SSPE cases correspond to structural, functional or antigenic areas, and could not be explained by the effects of cell passaging. The correlations between selected sites and functional studies of MV are discussed in detail with reference to the maintenance of persistent infection. No positive selection was found in the matrix (M) gene from acute cases of MV and the effects of including hypermutated SSPE M gene sequences in phylogenetic inference were also explored. Finally, using H gene data, we estimated the rate of molecular evolution for SSPE strains as 3·4×10−4substitutions/site/year, which is similar to previous estimates obtained for acute strains.

Variants of hepatitis C virus (HCV) have been classified by nucleotide sequence comparisons in different regions of the genome. Many investigators have defined the ranges of sequence similarity values or evolutionary distances corresponding to divisions of HCV into types, subtypes and isolates. Using these criteria, novel variants of HCV from Vietnam, Thailand and Indonesia have been classified as types 7, 8, 9, 10 and 11, many of which can be further subdivided into between two to four subtypes. In this study, this distance-based method of virus classification was compared with phylogenetic analysis and statistical measures to establish the confidence of the groupings. Using bootstrap resampling of phylogenetic trees in several subgenomic regions (core, E1, NS5) and with complete genomic sequences, we found that one set of novel HCV variants ('types 7, 8, 9 and 11') consistently grouped together into a single clade that also contained type 6a, while 'type 10a' grouped with type 3. In contrast, no robust higher-order groupings were observed between any of the other five previously described HCV genotypes (types 1-5). In each subgenomic region, the distribution of pairwise distances between members of the type 6 clade were consistently bi-modal and therefore provided no justification for classification of these variants into the three proposed categories (type, subtype, isolate). Based on these results, we propose that a more useful classification would regard all these variants as subtypes of type 6 or type 3, even though the level of sequence diversity within the clade was greater than observed for other genotypes. Classification by phylogenetic relatedness rules out simple sequence similarity measurements as a method for assigning HCV genotypes, but provides a more appropriate description of the evolutionary and epidemiological history of a virus.

Genotypes of hepatitis C virus (HCV) present within 104 samples from HCV-infected individuals from Africa, the Middle East, the Indian subcontinent and South-East Asia were identified by sequence comparisons in the core and NS-5 regions. Relatively short sequences (such as the 222 bp fragment of NS-5) provided effective discrimination of types, subtypes and isolates, and produced equivalent relationships between genotypes as were found upon comparison of longer sequences of NS-5, of the core region, and by comparison of the limited number of complete genomic sequences currently available. Measurement of evolutionary distances in the core and NS-5 regions allowed 79 of the 104 samples to be identified as examples of known genotypes, while 17 of the remainder could be provisionally classified as new subtypes of types 1 (Nigeria), 2 (Gambia), 3 (India, Pakistan and Bangladesh) and 4 (Middle East) on the basis of sequence comparison in core and NS-5 (n = 9) or provisionally using core alone (n = 8). The remaining sequences from Thailand made up two groups showing no close similarity to any of the six major genotypes classified to date, although one corresponded to an as yet unclassified variant of HCV also found in Thailand. However, phylogenetic analysis of the core and NS-5 regions indicated a distant relationship between these sequences with variants found in Vietnam and with type 6a, and collectively they formed a diverse single phylogenetic group. The existence of great diversity within a single genotype was also found amongst type 3 sequences in the Indian subcontinent, amongst type 4 variants in Central Africa and the Middle East, and amongst type variants in Nigeria. These findings may provide clues for understanding the origins and mechanisms of transmission of HCV.

A previous phylogenetic study suggested that mammalian gammaretroviruses may have originated in bats. Here we report the discovery of RNA transcripts from two putative endogenous gammaretroviruses in frugivorous (Rousettus leschenaultii retrovirus, RlRV) and insectivorous (Megaderma lyra retrovirus, MlRV) bat species. Both genomes possess a large deletion in pol, indicating that they are defective retroviruses. Phylogenetic analysis places RlRV and MlRV within the diversity of mammalian gammaretroviruses, with the former falling closer to porcine endogenous retroviruses and the latter to Mus dunni endogenous virus, koala retrovirus and gibbon ape leukemia virus. Additional genomic mining suggests that both microbat (Myotis lucifugus) and megabat (Pteropus vampyrus) genomes harbour many copies of endogenous retroviral forms related to RlRV and MlRV. Furthermore, phylogenetic analysis reveals the presence of three genetically diverse groups of endogenous gammaretroviruses in bat genomes, with M. lucifugus possessing members of all three groups. Taken together, this study indicates that bats harbour distinct gammaretroviruses and may have played an important role as reservoir hosts during the diversification of mammalian gammaretroviruses.

The genetic diversity present in populations of RNA viruses is likely to be strongly modulated by aspects of their life history, including mode of transmission. However, how transmission mode shapes patterns of intra- and inter-host genetic diversity, particularly when acting in combination with de novo mutation, population bottlenecks and the selection of advantageous mutations, is poorly understood. To address these issues, this study performed ultradeep sequencing of zucchini yellow mosaic virus in a wild gourd, Cucurbita pepo ssp. texana, under two infection conditions: aphid vectored and mechanically inoculated, achieving a mean coverage of approximately 10 000×. It was shown that mutations persisted during inter-host transmission events in both the aphid vectored and mechanically inoculated populations, suggesting that the vector-imposed transmission bottleneck is not as extreme as previously supposed. Similarly, mutations were found to persist within individual hosts, arguing against strong systemic bottlenecks. Strikingly, mutations were seen to go to fixation in the aphid-vectored plants, suggestive of a major fitness advantage, but remained at low frequency in the mechanically inoculated plants. Overall, this study highlights the utility of ultradeep sequencing in providing high-resolution data capable of revealing the nature of virus evolution, particularly as the full spectrum of genetic diversity within a population may not be uncovered without sequence coverage of at least 2500-fold.

Papillomaviruses (PVs) infect a wide range of vertebrates and have diversified into multiple genetic types, some of which have serious consequences for human health. Although PVs have to date only been characterized as exogenous viral forms, here we report the observation of an endogenous viral element (EPVLoa) in the genome of the platypus (Ornithorhynchus anatinus) that is related to PVs. Further data mining for endogenous PV-like elements is therefore warranted.

GII.4 noroviruses are a major cause of acute gastroenteritis in humans. A new variant of GII.4, the 2008 variant, has recently increased its prevalence on a global scale. A previous study of this variant in Japan suggested that it might be of recombinant origin, with a breakpoint at the ORF1–ORF2 junction. Here, examination of the evolutionary origin of the 2008 variant based on a larger sample of worldwide GII.4 norovirus sequences revealed a more complex pattern of recombination between the 2006a- and 2006b-like variants of genotype GII.4, involving the P2 antigenic domain. Double (termed ‘2008i’) and triple (termed ‘2008ii’) recombinant forms of 2008 variants were identified. This study highlights the possible importance of intra-genotypic recombination over antigenic regions in driving norovirus evolution, and is suggestive of a process analogous to the antigenic shift of influenza A virus by reassortment.

Isolates of hepatitis C virus (HCV) show considerable nucleotide sequence variability throughout the genome. Comparisons of complete genome sequences have been used as the basis of classification of HCV into a number of genotypes that show 67 to 77% sequence similarity. In order to investigate whether sequence relationships between genotypes are equivalent in different regions of the genome, we have carried out formal sequence analysis of variants in the 5' non-coding region (5'NCR) and in the genes encoding the core protein, an envelope protein (E1) and a non-structural protein (NS-5). In the E1 region, variants grouped into a series of six major genotypes and a series of subtypes that could be matched to the phylogenetic groupings previously observed for the NS-5 region. Furthermore, core and E1 sequences showed three non-overlapping ranges of sequence similarity corresponding to those between different genotypes, subtypes and isolates previously described in NS-5. Each major genotype could also be reliably identified by sequence comparisons in the well conserved 5'NCR, although many subtypes, such as 1a/1b, 2a/2c and some of those of type 4, could not be reliably distinguished from each other in this region. These data indicate that subgenomic regions such as E1 and NS-5 contain sufficient phylogenetic information for the identification of each of the 11 or 12 known types and subtypes of HCV. No evidence was found for variants of HCV that had sequences of one genotype in the 5'NCR but of a different one in the E1 or NS-5 region. This suggests that recombination between different HCV types is rare or non-existent and does not currently pose a problem in the use of subgenomic regions in classification.