Extensive Diversity of RNA Viruses in Australian Ticks

Virome of bat-infesting arthropods: highly divergent viruses in different vectors

Bats are reservoirs of important zoonotic viruses like Nipah and SARS viruses. However, whether the blood-sucking arthropods on the body surface of bats also carry these viruses, and the relationship between viruses carried by the blood-sucking arthropods and viruses carried by bats, have not been reported. This study collected 686 blood-sucking arthropods on the body surface of bats from Yunnan Province, China between 2012 and 2015, and they included wingless bat flies, bat flies, ticks, mites, and fleas. The viruses carried by these arthropods were analyzed using meta-transcriptomic approach, and 144 highly diverse positive-sense single-stranded RNA, negative-sense single-stranded RNA, and double-stranded RNA viruses were found, of which 138 were potentially new viruses. These viruses were classified into 14 different virus families or orders, including Bunyavirales, Mononegavirales, Reoviridae, and Picornavirales. Further analyses found that Bunyavirales were the most abundant virus group (84% of total virus RNA) in ticks, whereas narnaviruses were the most abundant (52-92%) in the bat flies and wingless bat flies libraries, followed by solemoviruses (1-29%) and reoviruses (0-43%). These viruses were highly structured based on the arthropod types. It is worth noting that no bat-borne zoonotic viruses were found in the virome of bat-infesting arthropod, seemly not supporting that bat surface arthropods are vectors of zoonotic viruses carried by bats. IMPORTANCE Bats are reservoir of many important viral pathogens. To evaluate whether bat-parasitic blood-sucking arthropods participate in the circulation of these important viruses, it is necessary to conduct unbiased virome studies on these arthropods. We evaluated five types of blood-sucking parasitic arthropods on the surface of bats in Yunnan, China and identified a variety of viruses, some of which had high prevalence and abundance level, although there is limited overlap in virome between distant arthropods. While most of the virome discovered here are potentially arthropod-specific viruses, we identified three possible arboviruses, including one orthobunyavirus and two vesiculoviruses (family Rhabdoviridae), suggesting bat-parasitic arthropods carry viruses with risk of spillage, which warrants further study.

Viromes and surveys of RNA viruses in camel-derived ticks revealing transmission patterns of novel tick-borne viral pathogens in Kenya

Tick-borne viruses (TBVs) capable of transmitting between ticks and hosts have been increasingly recognized as a global public health concern. In this study, Hyalomma ticks and serum samples from camels were collected using recorded sampling correlations in eastern Kenya. Viromes of pooled ticks were profiled by metagenomic sequencing, revealing a diverse community of viruses related to at least 11 families. Five highly abundant viruses, including three novel viruses (Iftin tick virus, Mbalambala tick virus [MATV], and Bangali torovirus [BanToV]) and new strains of previously identified viruses (Bole tick virus 4 [BLTV4] and Liman tick virus [LMTV]), were characterized in terms of genome sequences, organizations, and phylogeny, and their molecular prevalence was investigated in individual ticks. Moreover, viremia and antibody responses to these viruses have been investigated in camels. MATV, BLTV4, LMTV, and BanToV were identified as viral pathogens that can potentially cause zoonotic diseases. The transmission patterns of these viruses were summarized, suggesting three different types according to the sampling relationships between viral RNA-positive ticks and camels positive for viral RNA and/or antibodies. They also revealed the frequent transmission of BanToV and limited but effective transmission of other viruses between ticks and camels. Furthermore, follow-up surveys on TBVs from tick, animal, and human samples with definite sampling relationships are suggested. The findings revealed substantial threats from the emerging TBVs and may guide the prevention and control of TBV-related zoonotic diseases in Kenya and in other African countries.

Discovery and Surveillance of Tick-Borne Pathogens

Within the past 30 yr molecular assays have largely supplanted classical methods for detection of tick-borne agents. Enhancements provided by molecular assays, including speed, throughput, sensitivity, and specificity, have resulted in a rapid increase in the number of newly characterized tick-borne agents. The use of unbiased high throughput sequencing has enabled the prompt identification of new pathogens and the examination of tick microbiomes. These efforts have led to the identification of hundreds of new tick-borne agents in the last decade alone. However, little is currently known about the majority of these agents beyond their phylogenetic classification. Our article outlines the primary methods involved in tick-borne agent discovery and the current status of our understanding of tick-borne agent diversity.

Exploration of the Ixodes ricinus virosphere unveils an extensive virus diversity including novel coltiviruses and other reoviruses

Recent metagenomics studies have revealed several tick species to host a variety of previously undiscovered RNA viruses. Ixodes ricinus, which is known to be a vector for many viral, bacterial, and protozoan pathogens, is the most prevalent tick species in Europe. For this study, we decided to investigate the virosphere of Belgian I. ricinus ticks. High-throughput sequencing of tick pools collected from six different sampling sites revealed the presence of viruses belonging to many different viral orders and families, including Mononegavirales, Bunyavirales, Partitiviridae, and Reoviridae. Of particular interest was the detection of several new reoviruses, two of which cluster together with members of the genus Coltivirus. This includes a new strain of Eyach virus, a known causative agent of tick-borne encephalitis. All genome segments of this new strain are highly similar to those of previously published Eyach virus genomes, except for the fourth segment, encoding VP4, which is markedly more dissimilar, potentially indicating the occurrence of a genetic reassortment. Further polymerase chain reaction-based screening of over 230 tick pools for 14 selected viruses showed that most viruses could be found in all six sampling sites, indicating the wide spread of these viruses throughout the Belgian tick population. Taken together, these results illustrate the role of ticks as important virus reservoirs, highlighting the need for adequate tick control measures.

Interactive Gene Expression Between Metarhizium anisopliae JEF-290 and Longhorned Tick Haemaphysalis longicornis at Early Stage of Infection

The longhorned tick, Haemaphysalis longicornis (Acari: Ixodidae), is a hard tick and a vector for severe fever with thrombocytopenia syndrome (SFTS) virus. The number of patients infected with SFTS is rapidly increasing. Recently, the invertebrate pathogen Metarhizium anisopliae JEF-290 was reported to be useful to control the tick as an alternative to chemical acaricides, which are not easily applicable in human living areas where the tick is widely spread. In this study, we analyzed how the tick and the fungal pathogen interact at the transcriptional level. Field-collected tick nymphs were treated with JEF-290 conidia at 1 × 10⁸ conidia/ml. In the early stage of infection with 2.5% mortality, the infected ticks were subjected to RNA sequencing, and non-infected ticks and fungal masses served as controls. Fungus and tick genes were mostly up-regulated at the early stage of infection. In the gene set enrichment analysis of the infecting fungus, catabolic processes that included lipids, phospholipids, and detoxification processes, the response to oxidative stress, and toxic substances were significantly up-regulated. In this fungal up-regulation, various lipase, antioxidant enzyme, and hydrolase genes were highly transcribed. The gene set enrichment analysis of the infected tick showed that many peptide synthesis processes including translation, peptide metabolism, ribonucleotide metabolism, and energy production processes that included ATP generation and ADP metabolism were significantly up-regulated. Structurally, mitochondria and ribosome subunit genes in ticks were highly transcribed to upregulate these processes. Together these results indicate that JEF-290 initiates process that infects the tick while the tick actively defends against the fungal attack. This work provides background to improve our understanding of the early stage of fungal infection in longhorned tick.

Applications of Blocker Nucleic Acids and Non-Metazoan PCR Improves the Discovery of the Eukaryotic Microbiome in Ticks

Ticks serve as important vectors of a variety of pathogens. Recently, the viral and prokaryotic microbiomes in ticks have been explored using next-generation sequencing to understand the physiology of ticks and their interactions with pathogens. However, analyses of eukaryotic communities in ticks are limited, owing to the lack of suitable methods. In this study, we developed new methods to selectively amplify microeukaryote genes in tick-derived DNA by blocking the amplification of the 18S rRNA gene of ticks using artificial nucleic acids: peptide nucleic acids (PNAs) and locked nucleic acids (LNAs). In addition, another PCR using non-metazoan primers, referred to as UNonMet-PCR, was performed for comparison. We performed each PCR using tick-derived DNA and sequenced the amplicons using the Illumina MiSeq platform. Almost all sequences obtained by conventional PCR were derived from ticks, whereas the proportion of microeukaryotic reads and alpha diversity increased upon using the newly developed method. Additionally, the PNA- or LNA-based methods were suitable for paneukaryotic analyses, whereas the UNonMet-PCR method was particularly sensitive to fungi. The newly described methods enable analyses of the eukaryotic microbiome in ticks. We expect the application of these methods to improve our understanding of the tick microbiome.

Unbiased Characterization of the Microbiome and Virome of Questing Ticks

Due to their vector capacity, ticks are ectoparasites of medical and veterinary significance. Modern sequencing tools have facilitated tick-associated microbiota studies, but these have largely focused on bacterial pathogens and symbionts. By combining 16S rRNA gene sequencing with total RNA-sequencing methods, we aimed to determine the complete microbiome and virome of questing, female Ixodes holocyclus recovered from coastal, north-eastern New South Wales (NSW), Australia. We present, for the first time, a robust and unbiased method for the identification of novel microbes in ticks that enabled us to identify bacteria, viruses, fungi and eukaryotic pathogens. The dominant bacterial endosymbionts were Candidatus Midichloria sp. Ixholo1 and Candidatus Midichloria sp. Ixholo2. Candidatus Neoehrlichia australis and Candidatus Neoehrlichia arcana were also recovered, confirming that these bacteria encompass I. holocyclus’ core microbiota. In addition, seven virus species were detected—four previously identified in I. holocyclus and three novel species. Notably, one of the four previously identified virus species has pathogenic potential based on its phylogenetic relationship to other tick-associated pathogens. No known pathogenic eukaryotes or fungi were identified. This study has revealed the microbiome and virome of female I. holocyclus from the environment in north-eastern NSW. We propose that future tick microbiome and virome studies utilize equivalent methods to provide an improved representation of the microbial diversity in ticks globally.

The Efficacy of Ultrasonic Pest Repellent Devices against the Australian Paralysis Tick, Ixodes holocyclus (Acari: Ixodidae)

Ultrasonic pest repellers are often promoted as a means of protecting people and pets from the bites of hematophagous arthropods, such as ticks. However, to date, there has been no published research on the effectiveness of these devices against the Australian paralysis tick, Ixodes holocyclus Neumann. The purpose of this study was to test the effectiveness of nine ultrasonic devices against female I. holocyclus. Two arenas were constructed, one for the test (with the ultrasonic device) and one for the control (no device). Each arena had a test and an escape chamber, connected by a corridor. Twenty ticks were placed in each test chamber. After the ultrasonic device was operated for 1 h, the number of ticks in both chambers was recorded. Ten replicates were conducted for each device. The average number of ticks that moved from the test to the escape chamber was greater in all the test arenas, with three devices being statistically different from the control. However, the highest percent of ticks that escaped was only 19.5%. This amount is insufficient to offer adequate protection against tick bites and this study adds further weight to previous investigations that ultrasonic devices should not be employed in pest management.

Sequence diversity and evolution of a group of iflaviruses associated with ticks

We studied a group of tick-associated viruses with characteristics of members of the family Iflaviridae, a family of viruses frequently found in arthropods. Our aim was to gain insight into the evolutionary dynamics of this group of viruses, which may be linked to the biology of ticks. We explored assembled RNA-Seq data sets for different species of ticks. We identified members of five different iflavirus species, four of them novel, and discovered nine new genome sequences, including variants. Five variants represented a virus species associated with Ixodes ricinus. Unexpectedly, a sequence found in the Ixodes scapularis cell line ISE6 was nearly identical to the sequences of I. ricinus variants, suggesting a contamination of this cell line by I. ricinus material. Analysing patterns of substitutions between these variants, we detected a strong excess of synonymous mutations, suggesting evolution under strong positive selection. The phylogenies of the viruses and of their tick hosts were not congruent, suggesting recurrent host changes across tick genera during their evolution. Overall, our work constitutes a step in the understanding of the interactions between this family of viruses and ticks.

The virome of Rhipicephalus, Dermacentor and Haemaphysalis ticks from Eastern Romania includes novel viruses with potential relevance for public health

Ticks are involved in the transmission of various pathogens and several tick-borne diseases cause significant problems for the health of humans and livestock. The composition of viral communities in ticks and their interactions with pathogens, is poorly understood, particularly in Eastern Europe, an area that represents a major hub for animal-arthropod vectors exchanges (e.g., via bird migrations). The aim of this study was to describe the virome of Dermacentor sp., Rhipicephalus sp. and Haemaphysalis sp. ticks collected from relatively little studied regions of Romania (Iasi and Tulcea counties) located at the intersection of various biotopes, countries and routes of migrations. We also focused the study on viruses that could potentially have relevance for human and animal health. In 2019, more than 500 ticks were collected from the vegetation and from small ruminants and analysed by high-throughput transcriptome sequencing. Among the viral communities infecting Romanian ticks, viruses belonging to the Flaviviridae, Phenuiviridae and Nairoviridae families were identified and full genomes were derived. Phylogenetic analyses placed them in clades where mammalian isolates are found, suggesting that these viruses could constitute novel arboviruses. The characterization of these communities increase the knowledge of the diversity of viruses in Eastern Europe and provides a basis for further studies about the interrelationship between ticks and tick-borne viruses.

Molecular detection of Rickettsia raoultii, Rickettsia tamurae, and associated pathogens from ticks parasitizing water deer (Hydropotes inermis argyropus) in South Korea

Most defined Rickettsiales, which can be transmitted by ticks, are known to be important zoonotic pathogens. Some of these pathogens can cause severe diseases in humans, including anaplasmosis, rickettsioses, and ehrlichiosis. Previous studies in South Korea have investigated tick-borne pathogens (TBPs) residing in ticks found on grassy vegetation and animals. However, there is limited phylogenetic information on TBPs in ticks parasitizing Korean water deer (KWD; Hydropotes inermis argyropus). This study evaluated the prevalence, risk factors (regions, tick stages, and tick species), and coinfections of TBPs in ticks parasitizing KWD. Were collected a total of 283 hard ticks, including Haemaphysalis longicornis, Haemaphysalis flava, and Ixodes nipponensis from KWD in South Korea from 2013 to 2017. In 173 tested tick pools, genes for seven TBPs, namely Rickettsia raoultii (20 %), Rickettsia tamurae (1 %), Candidatus Rickettsia longicornii (31 %), Ehrlichia canis (3 %), Anaplasma capra (3 %), Anaplasma bovis (2 %), and Anaplasma sp. (1 %), were detected. The unidentified Anaplasma sp. isolates revealed a 98.4 %-99.3 % sequence identity with Anaplasma sp. in GenBank sequences obtained from ticks. To the best of our knowledge, this is the first study to report the presence of the emerging human pathogen R. tamurae in South Korea. These results should increase awareness about the need for continued development of epidemiological control measures, and medical and veterinary communities must be informed of their high infection potential and clinical complexity in humans.

Tick-virus interactions: Current understanding and future perspectives

Ticks are the primary vector of arboviruses in temperate climates worldwide. They are both the vector of these pathogens to humans and an integral component of the viral sylvatic cycle. Understanding the tick-pathogen interaction provides information about the natural maintenance of these pathogens and informs the development of countermeasures against human infection. In this review, we discuss currently available information on tick-viral interactions within the broader scope of general tick immunology. While the tick immune response to several pathogens has been studied extensively, minimal work centres on responses to viral infection. This is largely due to the high pathogenicity of tick-borne viruses; this necessitates high-containment laboratories or low-pathogenicity substitute viruses. This has biased most research towards tick-borne flaviviruses. More work is required to fully understand the role of tick-virus interaction in sylvatic cycling and transmission of diverse tick-borne viruses.

Parasites of wombats (family Vombatidae), with a focus on ticks and tick-borne pathogens

Ticks (Arachnida: Acari) are vectors for pathogens and the biggest threat to animal health. Many Australian ticks are associated with pathogens that impact humans, domestic animals and livestock. However, little is known about the presence or impact of tick-borne pathogens in native Australian wildlife. Wombats are particularly susceptible to the effects of the ectoparasite Sarcoptes scabiei which causes sarcoptic mange, the reason for which is unknown. Factors such as other ectoparasites and their associated pathogens may play a role. A critical understanding of the species of ectoparasites that parasitise wombats and their pathogens, and particularly ticks, is therefore warranted. This review describes the ectoparasites of wombats, pathogens known to be associated with those ectoparasites, and related literature gaps. Pathogens have been isolated in most tick species that typically feed on wombats; however, there are minimal molecular studies to determine the presence of pathogens in any other wombat ectoparasites. The development of next-generation sequencing (NGS) technologies allows us to explore entire microbial communities in ectoparasite samples, allowing fast and accurate identification of potential pathogens in many samples at once. These new techniques have highlighted the diversity and uniqueness of native ticks and their microbiomes, including pathogens of potential medical and veterinary importance. An increased understanding of all ectoparasites that parasitise wombats, and their associated pathogens, requires further investigation.

Molecular detection and genomic characterization of diverse hepaciviruses in African rodents

Hepatitis C virus (HCV; genus Hepacivirus) represents a major public health problem, infecting about three per cent of the human population. Because no animal reservoir carrying closely related hepaciviruses has been identified, the zoonotic origins of HCV still remain unresolved. Motivated by recent findings of divergent hepaciviruses in rodents and a plausible African origin of HCV genotypes, we have screened a large collection of small mammals samples from seven sub-Saharan African countries. Out of 4,303 samples screened, eighty were found positive for the presence of hepaciviruses in twenty-nine different host species. We, here, report fifty-six novel genomes that considerably increase the diversity of three divergent rodent hepacivirus lineages. Furthermore, we provide strong evidence for hepacivirus co-infections in rodents, which were exclusively found in four sampled species of brush-furred mice. We also detect evidence of recombination within specific host lineages. Our study expands the available hepacivirus genomic data and contributes insights into the relatively deep evolutionary history of these pathogens in rodents. Overall, our results emphasize the importance of rodents as a potential hepacivirus reservoir and as models for investigating HCV infection dynamics.

Soybean Thrips (Thysanoptera: Thripidae) Harbor Highly Diverse Populations of Arthropod, Fungal and Plant Viruses

Soybean thrips (Neohydatothrips variabilis) are one of the most efficient vectors of soybean vein necrosis virus, which can cause severe necrotic symptoms in sensitive soybean plants. To determine which other viruses are associated with soybean thrips, the metatranscriptome of soybean thrips, collected by the Midwest Suction Trap Network during 2018, was analyzed. Contigs assembled from the data revealed a remarkable diversity of virus-like sequences. Of the 181 virus-like sequences identified, 155 were novel and associated primarily with taxa of arthropod-infecting viruses, but sequences similar to plant and fungus-infecting viruses were also identified. The novel viruses were predicted to have positive-sense RNA, negative-stranded RNA, double-stranded RNA, and single-stranded DNA genomes. The assembled sequences included 100 contigs that represented at least 95% coverage of a virus genome or genome segment. Sequences represented 12 previously described arthropod viruses including eight viruses reported from Hubei Province in China, and 12 plant virus sequences of which six have been previously described. The presence of diverse populations of plant viruses within soybean thrips suggests they feed on and acquire viruses from multiple host plant species that could be transmitted to soybean. Assessment of the virome of soybean thrips provides, for the first time, information on the diversity of viruses present in thrips.

Targeted Metagenomics Offers Insights into Potential Tick-Borne Pathogens

Tick-borne disease pathogen identification remains a diagnostic challenge due to the multiple tests necessary for diagnosis. Targeted metagenomic sequencing is an evolving method to detect multiple different pathogens with a single test. In this issue of the Journal of Clinical Microbiology, L. Kingry, S. Sheldon, S. Oatman, B. Pritt, et al. (J Clin Microbiol 58:e00147-20, 2020, https://doi.org/10.1128/JCM.00147-20) explore 16S rRNA gene-targeted metagenomics to identify bacterial pathogens in 881 of more than 13,000 specimens submitted for tick-borne pathogen testing, giving new insights into suspected tick-borne pathogens.

Unmapped RNA Virus Diversity in Termites and their Symbionts

Despite their ecological importance, nothing is known about the diversity and abundance of RNA viruses in termites (Termitoidae). We used a metatranscriptomics approach to determine the RNA virome structure of 50 diverse species of termite that differ in both phylogenetic position and colony composition. From these samples, we identified 67 novel RNA viruses, characterized their genomes, quantified their abundance and inferred their evolutionary history. These viruses were found within or similar to those from the Togaviridae, Iflaviridae, Polycipiviridae, Flaviviridae, Leviviridae, Narnaviridae, Mitoviridae, Lispivirdae, Phasmaviridae, Picobirnaviridae and Partitiviridae. However, all viruses identified were novel and divergent, exhibiting only 20% to 45% amino acid identity to previously identified viruses. Our analysis suggested that 17 of the viruses identified were termite-infecting, with the remainder likely associated with the termite microbiome or diet. Unclassified sobemo-like and bunya-like viruses dominated termite viromes, while most of the phylogenetic diversity was provided by the picobirna- and mitovirus-like viruses. Of note was the identification of a novel flavi-like virus most closely related to those found in marine vertebrates and invertebrates. Notably, the sampling procedure had the strongest association with virome composition, with greater RNA virome diversity in libraries prepared from whole termite bodies than those that only sampled heads.

Full-length genome sequence of segmented RNA virus from ticks was obtained using small RNA sequencing data

BACKGROUND

In 2014, a novel tick-borne virus of the Flaviviridae family was first reported in the Mogiana region of Brazil and named the Mogiana tick virus (MGTV). Thereafter, the Jingmen tick virus (JMTV), Kindia tick virus (KITV), and Guangxi tick virus (GXTV)-evolutionarily related to MGTV-were reported.

RESULTS

In the present study, we used small RNA sequencing (sRNA-seq) to detect viruses in ticks and discovered a new MGTV strain in Amblyomma testudinarium ticks collected in China's Yunnan Province in 2016. We obtained the full-length genome sequence of this MGTV strain Yunnan2016 (GenBank: MT080097, MT080098, MT080099 and MT080100) and recommended it for its inclusion in the NCBI RefSeq database for future studies on MGTV, JMTV, KITV and GXTV. Phylogenetic analysis showed that MGTV, JMTV, KITV and GXTV are monophyletic and belong to a MGTV group. Furthermore, this MGTV group of viruses may be phylogenetically related to geographical regions that were formerly part of the supercontinents Gondwana and Laurasia.

CONCLUSIONS

To the best of our knowledge, this is the first study in which 5' and 3' sRNAs were used to generate full-length genome sequences of, but not limited to, RNA viruses. We also demonstrated the feasibility of using the sRNA-seq based method for the detection of viruses in pooled two and even possible one small ticks. MGTV may preserve the characteristic of ancient RNA viruses, which can be used to study the origin and evolution of RNA viruses. In addition, MGTV can be used as novel species for studies in phylogeography.

Meta-transcriptomic identification of Trypanosoma spp. in native wildlife species from Australia

Background

Wildlife species carry a remarkable diversity of trypanosomes. The detection of trypanosome infection in native Australian fauna is central to understanding their diversity and host-parasite associations. The implementation of total RNA sequencing (meta-transcriptomics) in trypanosome surveillance and diagnosis provides a powerful methodological approach to better understand the host species distribution of this important group of parasites.

Methods

We implemented a meta-transcriptomic approach to detect trypanosomes in a variety of tissues (brain, liver, lung, skin, gonads) sampled from native Australian wildlife, comprising four marsupials (koala, Phascolarctos cinereus; southern brown bandicoot, Isoodon obesulus; swamp wallaby, Wallabia bicolor; bare-nosed wombat, Vombatus ursinus), one bird (regent honeyeater, Anthochaera phrygia) and one amphibian (eastern dwarf tree frog, Litoria fallax). Samples corresponded to both clinically healthy and diseased individuals. Sequencing reads were de novo assembled into contigs and annotated. The evolutionary relationships among the trypanosomatid sequences identified were determined through phylogenetic analysis of 18S rRNA sequences.

Results

We detected trypanosome sequences in all six species of vertebrates sampled, with positive samples in multiple organs and tissues confirmed by PCR. Phylogenetic analysis indicated that the trypanosomes infecting marsupials were related to those previously detected in placental and marsupial mammals, while the trypanosome in the regent honeyeater grouped with avian trypanosomes. In contrast, we provide the first evidence for a trypanosome in the eastern dwarf tree frog that was phylogenetically distinct from those described in other amphibians.

Conclusions

To our knowledge, this is the first meta-transcriptomic analysis of trypanosomes in native Australian wildlife, expanding the known genetic diversity of these important parasites. We demonstrated that RNA sequencing is sufficiently sensitive to detect low numbers of Trypanosoma transcripts and from diverse hosts and tissues types, thereby representing an effective means to detect trypanosomes that are divergent in genome sequence.

Depicting the RNA Virome of Hematophagous Arthropods from Belgrade, Serbia

Hematophagous arthropods are important vectors for zoonotic pathogens. To date, a huge number of viruses have been identified in these arthropods, with a considerable proportion of them being human pathogens. However, the viromes of hematophagous arthropods are still largely unresearched. In this study, a number of arthropods were collected from Belgrade, Serbia including mosquitoes, ticks and bedbugs. The viromes of these arthropods were identified and characterized using Illumina MiSeq sequencing. In total, 21 viruses belonging to 11 families were characterized, with 11 of them representing novel species. These results may contribute to our knowledge of RNA viruses in arthropods and the discovery of novel human pathogens.

Novel hepaci- and pegi-like viruses in native Australian wildlife and non-human primates

The Flaviviridae family of positive-sense RNA viruses contains important pathogens of humans and other animals, including Zika virus, dengue virus, and hepatitis C virus. The Flaviviridae are currently divided into four genera-Hepacivirus, Pegivirus, Pestivirus, and Flavivirus-each with a diverse host range. Members of the genus Hepacivirus are associated with an array of animal species, including humans, non-human primates, other mammalian species, as well as birds and fish, while the closely related pegiviruses have been identified in a variety of mammalian taxa, also including humans. Using a combination of total RNA and whole-genome sequencing we identified four novel hepaci-like viruses and one novel variant of a known hepacivirus in five species of Australian wildlife. The hosts infected comprised native Australian marsupials and birds, as well as a native gecko (Gehyra lauta). From these data we identified a distinct marsupial clade of hepaci-like viruses that also included an engorged Ixodes holocyclus tick collected while feeding on Australian long-nosed bandicoots (Perameles nasuta). Distinct lineages of hepaci-like viruses associated with geckos and birds were also identified. By mining the SRA database we similarly identified three new hepaci-like viruses from avian and primate hosts, as well as two novel pegi-like viruses associated with primates. The phylogenetic history of the hepaci- and pegi-like viruses as a whole, combined with co-phylogenetic analysis, provided support for virus-host co-divergence over the course of vertebrate evolution, although with frequent cross-species virus transmission. Overall, our work highlights the diversity of the Hepacivirus and Pegivirus genera as well as the uncertain phylogenetic distinction between.

Circumpolar diversification of the Ixodes uriae tick virome

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.

As host populations diverge, so may those microorganisms, including viruses, that are dependent on those hosts. To examine this key issue in host-microbe evolution we compared the co-phylogenies of the seabird tick, Ixodes uriae, and their RNA viruses sampled from the far northern and southern hemispheres. Despite the huge geographic distance between them, phylogeographic analysis reveals that the same and closely related viruses were found both within and between the northern and southern circumpolar regions, most likely reflecting transfer by virus-infected migratory birds. In contrast, genomic data suggested that the Ix. uriae populations were phylogenetically distinct between the northern and southern hemispheres. This work emphasises the importance of migratory birds and ticks as vectors and sources of virus dispersal and introduction at both the local and global scales.

Characterization of the complete mitochondrial genome of Amblyomma ovale, comparative analyses and phylogenetic considerations

In this study, we sequenced two complete mitochondrial genomes of Amblyomma ovale, a tick of public health importance. Sequencing two distinct individuals, the resulting mitochondrial genomes were 14,756 and 14,760 bp in length and maintained the same gene order previously reported in Amblyomma. These were combined with RNA-seq derived mitochondrial sequences from three additional species, Amblyomma aureolatum, Amblyomma maculatum, and Amblyomma moreliae, to carry out mitogenome comparative and evolutionary analyses against all previously published tick mitochondrial genomes. We described a derivative genome rearrangement that isolates Ixodes from the remaining Ixodidae and consists of both a reverse translocation as well as an event of Tandem Duplication Random Loss. Genetic distance analyses indicated that cox2, nd1, nd5, and 16S are good candidates for future population studies in A. ovale. The phylogenetic analyses corroborated the utility of complete mitochondrial genomes as phylogenetic markers within the group. This study further supplements the genome information available for Amblyomma and facilitates future evolutionary and population genetic studies within the genus.

Tolerance of Honey Bees to Varroa Mite in the Absence of Deformed Wing Virus

The global spread of the parasitic mite Varroa destructor has emphasized the significance of viruses as pathogens of honey bee (Apis mellifera) populations. In particular, the association of deformed wing virus (DWV) with V. destructor and its devastating effect on honey bee colonies has led to that virus now becoming one of the most well-studied insect viruses. However, there has been no opportunity to examine the effects of Varroa mites without the influence of DWV. In Papua New Guinea (PNG), the sister species, V. jacobsoni, has emerged through a host-shift to reproduce on the local A. mellifera population. After initial colony losses, beekeepers have maintained colonies without chemicals for more than a decade, suggesting that this bee population has an unknown mite tolerance mechanism. Using high throughput sequencing (HTS) and target PCR detection, we investigated whether the viral landscape of the PNG honey bee population is the underlying factor responsible for mite tolerance. We found A. mellifera and A. cerana from PNG and nearby Solomon Islands were predominantly infected by sacbrood virus (SBV), black queen cell virus (BQCV) and Lake Sinai viruses (LSV), with no evidence for any DWV strains. V. jacobsoni was infected by several viral homologs to recently discovered V. destructor viruses, but Varroa jacobsoni rhabdovirus-1 (ARV-1 homolog) was the only virus detected in both mites and honey bees. We conclude from these findings that A. mellifera in PNG may tolerate V. jacobsoni because the damage from parasitism is significantly reduced without DWV. This study also provides further evidence that DWV does not exist as a covert infection in all honey bee populations, and remaining free of this serious viral pathogen can have important implications for bee health outcomes in the face of Varroa.

Evaluation of Viral RNA Recovery Methods in Vectors by Metagenomic Sequencing

Identification and characterization of viral genomes in vectors including ticks and mosquitoes positive for pathogens of great public health concern using metagenomic next generation sequencing (mNGS) has challenges. One such challenge is the ability to efficiently recover viral RNA which is typically dependent on sample processing. We evaluated the quantitative effect of six different extraction methods in recovering viral RNA in vectors using negative tick homogenates spiked with serial dilutions of tick-borne encephalitis virus (TBEV) and surrogate Langat virus (LGTV). Evaluation was performed using qPCR and mNGS. Sensitivity and proof of concept of optimal method was tested using naturally positive TBEV tick homogenates and positive dengue, chikungunya, and Zika virus mosquito homogenates. The amount of observed viral genome copies, percentage of mapped reads, and genome coverage varied among different extractions methods. The developed Method 5 gave a 120.8-, 46-, 2.5-, 22.4-, and 9.9-fold increase in the number of viral reads mapping to the expected pathogen in comparison to Method 1, 2, 3, 4, and 6, respectively. Our developed Method 5 termed ROVIV (Recovery of Viruses in Vectors) greatly improved viral RNA recovery and identification in vectors using mNGS. Therefore, it may be a more sensitive method for use in arbovirus surveillance.

Comparative Analysis of RNA Virome Composition in Rabbits and Associated Ectoparasites

Ectoparasites play an important role in the transmission of many vertebrate-infecting viruses, including Zika and dengue viruses. Although it is becoming increasingly clear that invertebrate species harbor substantial virus diversity, it is unclear how many of the viruses carried by invertebrates have the potential to infect vertebrate species. We used the European rabbit (Oryctolagus cuniculus) as a model species to compare virome compositions in a vertebrate host and known associated ectoparasite mechanical vectors, in this case, fleas and blowflies. In particular, we aimed to infer the extent of viral transfer between these distinct types of host. Our analysis revealed that despite extensive viral diversity in both rabbits and associated ectoparasites, and the close interaction of these vertebrate and invertebrate species, biological viral transmission from ectoparasites to vertebrate species is rare. We did, however, find evidence to support the idea of a role of blowflies in transmitting viruses without active replication in the insect.

Ectoparasites play an important role in virus transmission among vertebrates. Little, however, is known about the nature of those viruses that pass between invertebrates and vertebrates. In Australia, flies and fleas support the mechanical transmission of two viral biological controls against wild rabbits—rabbit hemorrhagic disease virus (RHDV) and myxoma virus. We compared virome compositions in rabbits and these ectoparasites, sequencing total RNA from multiple tissues and gut contents of wild rabbits, fleas collected from these rabbits, and flies trapped sympatrically. Meta-transcriptomic analyses identified 50 novel viruses from multiple RNA virus families. Rabbits and their ectoparasites were characterized by markedly different viromes, with virus abundance greatest in flies. Although viral contigs from six virus families/groups were found in both rabbits and ectoparasites, they clustered in distinct host-dependent lineages. A novel calicivirus and a picornavirus detected in rabbit cecal content were vertebrate specific; the newly detected calicivirus was distinct from known rabbit caliciviruses, while the picornavirus clustered with sapeloviruses. Several picobirnaviruses were also identified that fell in diverse phylogenetic positions, compatible with the idea that they are associated with bacteria. Further comparative analysis revealed that the remaining viruses found in rabbits, and all those from ectoparasites, were likely associated with invertebrates, plants, and coinfecting endosymbionts. While no full genomes of vertebrate-associated viruses were detected in ectoparasites, small numbers of reads from rabbit astrovirus, RHDV, and other lagoviruses were present in flies. This supports a role for flies in the mechanical transmission of RHDV, while their involvement in astrovirus transmission merits additional exploration.

IMPORTANCE Ectoparasites play an important role in the transmission of many vertebrate-infecting viruses, including Zika and dengue viruses. Although it is becoming increasingly clear that invertebrate species harbor substantial virus diversity, it is unclear how many of the viruses carried by invertebrates have the potential to infect vertebrate species. We used the European rabbit (Oryctolagus cuniculus) as a model species to compare virome compositions in a vertebrate host and known associated ectoparasite mechanical vectors, in this case, fleas and blowflies. In particular, we aimed to infer the extent of viral transfer between these distinct types of host. Our analysis revealed that despite extensive viral diversity in both rabbits and associated ectoparasites, and the close interaction of these vertebrate and invertebrate species, biological viral transmission from ectoparasites to vertebrate species is rare. We did, however, find evidence to support the idea of a role of blowflies in transmitting viruses without active replication in the insect.

Sustained RNA virome diversity in Antarctic penguins and their ticks

Despite its isolation and extreme climate, Antarctica is home to diverse fauna and associated microorganisms. It has been proposed that the most iconic Antarctic animal, the penguin, experiences low pathogen pressure, accounting for their disease susceptibility in foreign environments. There is, however, a limited understanding of virome diversity in Antarctic species, the extent of in situ virus evolution, or how it relates to that in other geographic regions. To assess whether penguins have limited microbial diversity we determined the RNA viromes of three species of penguins and their ticks sampled on the Antarctic peninsula. Using total RNA sequencing we identified 107 viral species, comprising likely penguin associated viruses (n = 13), penguin diet and microbiome associated viruses (n = 82), and tick viruses (n = 8), two of which may have the potential to infect penguins. Notably, the level of virome diversity revealed in penguins is comparable to that seen in Australian waterbirds, including many of the same viral families. These data run counter to the idea that penguins are subject to lower pathogen pressure. The repeated detection of specific viruses in Antarctic penguins also suggests that rather than being simply spill-over hosts, these animals may act as key virus reservoirs.

Newly identified viral genomes in pangolins with fatal disease

Epizootic pathogens pose a major threat to many wildlife species, particularly in the context of rapidly changing environments. Pangolins (order Pholidota) are highly threatened mammals, in large part due to the trade in illegal wildlife. During July to August 2018 four sick wild pangolins (three Manis javanica and one Manis pentadactyla) exhibiting a variety of clinical symptoms were rescued by the Jinhua Wildlife Protection Station in Zhejiang province, China. Although three of these animals died, fortunately one recovered after 2 weeks of symptomatic treatment. Using meta-transcriptomics combined with reverse transcription polymerase chain reaction (RT-PCR), we identified two novel RNA viruses in two of the dead pangolins. Genomic analysis revealed that these viruses were most closely related to pestiviruses and coltiviruses, although still highly genetically distinct, with more than 48 and 25 per cent sequence divergence at the amino acid level, respectively. We named these Dongyang pangolin virus (DYPV) and Lishui pangolin virus (LSPV) based on the sampling site and hosts. Although coltiviruses (LSPV) are known to be transmitted by ticks, we found no evidence of LSPV in ticks sampled close to where the pangolins were collected. In addition, although DYPV was present in nymph ticks (Amblyomma javanense) collected from a diseased pangolin, they were not found in the local tick population. Epidemiological investigation revealed that both novel viruses might have been imported following the illegal international trade of pangolins. Hence, these data indicate that illegal wildlife trafficking not only threatens the status of pangolin populations, but may also spread epizootic pathogens.

Metatranscriptomic virome assessment of Rhipicephalus microplus from Colombia

Ticks (Ixodida) are hematophagous ectoparasites that harbor and transmit diverse species of viruses, some of which cause serious diseases with worldwide veterinary and human health consequences. Rhipicephalus microplus is an important cattle tick in Colombia, where it causes significant economic losses. Despite the importance of this tick, its viral profile is unknown. RNA sequencing was used in this study as a surveillance method for virus detection in R. microplus. Most of the viral origin contigs were assigned to two putative viruses: one chuvirus (Wuhan tick virus 2) and one phlebovirus-like (Lihan tick virus). In addition, viral contigs corresponding to two jingmenviruses previously reported in R. microplus from China and Brazil were detected, as well as a novel putative tymovirus, named here as Antioquia tymovirus-like 1 (ATV-like 1). The presence of some of these viruses across numerous regions in the world could have several explanations, including i) a long-term association between those viruses and R. microplus and ii) a consequence of livestock historical trade. Our results shed new light on the virus diversity of this tick species and provide a basis for further studies on the evolutionary history and pathogenic potential of these interesting viruses.

RNA Viruses of Amblyomma variegatum and Rhipicephalus microplus and Cattle Susceptibility in the French Antilles

Ticks transmit a wide variety of pathogens including bacteria, parasites and viruses. Over the last decade, numerous novel viruses have been described in arthropods, including ticks, and their characterization has provided new insights into RNA virus diversity and evolution. However, little is known about their ability to infect vertebrates. As very few studies have described the diversity of viruses present in ticks from the Caribbean, we implemented an RNA-sequencing approach on Amblyomma variegatum and Rhipicephalus microplus ticks collected from cattle in Guadeloupe and Martinique. Among the viral communities infecting Caribbean ticks, we selected four viruses belonging to the Chuviridae, Phenuiviridae and Flaviviridae families for further characterization and designing antibody screening tests. While viral prevalence in individual tick samples revealed high infection rates, suggesting a high level of exposure of Caribbean cattle to these viruses, no seropositive animals were detected. These results suggest that the Chuviridae- and Phenuiviridae-related viruses identified in the present study are more likely tick endosymbionts, raising the question of the epidemiological significance of their occurrence in ticks, especially regarding their possible impact on tick biology and vector capacity. The characterization of these viruses might open the door to new ways of preventing and controlling tick-borne diseases.

Discovery of Jogalong virus, a novel hepacivirus identified in a Culex annulirostris (Skuse) mosquito from the Kimberley region of Western Australia

The discovery of hepaciviruses in non-human hosts has accelerated following the advancement of high-throughput sequencing technology. Hepaciviruses have now been described in reptiles, fish, birds, and an extensive array of mammals. Using metagenomic sequencing on pooled samples of field-collected Culex annulirostris mosquitoes, we discovered a divergent hepacivirus-like sequence, named Jogalong virus, from the Kimberley region in northern Western Australia. Using PCR, we screened the same 300 individual mosquitoes and found just a single positive sample (1/300, 0.33%). Phylogenetic analysis of the hepacivirus NS5B protein places Jogalong virus within the genus Hepacivirus but on a distinct and deeply rooted monophyletic branch shared with duck hepacivirus, suggesting a notably different evolutionary history. Vertebrate barcoding PCR targeting two mitochondrial genes, cytochrome c oxidase subunit I and cytochrome b, indicated that the Jogalong virus-positive mosquito had recently fed on the tawny frogmouth (Podargus strigoides), although it is currently unknown whether this bird species contributes to the natural ecology of this virus.

Revisiting new tick-associated viruses: what comes next?

Tick-borne viral infections continue to cause diseases with considerable impact on humans, livestock, companion animals and wildlife. Many lack specific therapeutics and vaccines are available for only a few. Tick-borne viruses will continue to emerge, facilitated by anthroponotic factors related to the modern lifestyle. We persistently identify and are obliged to cope with new examples of emerging tick-borne viral diseases and novel viruses today. Many new strains have been detected in vertebrates and arthropods, some causing severe diseases likely to challenge public and veterinary health. This manuscript aims to provide a narrative overview of recently-described tick-associated viruses, with perspectives on changing paradigms in identification, screening and control.

Sensitivity and specificity of metatranscriptomics as an arbovirus surveillance tool

The ability to identify all the viruses within a sample makes metatranscriptomic sequencing an attractive tool to screen mosquitoes for arboviruses. Practical application of this technique, however, requires a clear understanding of its analytical sensitivity and specificity. To assess this, five dilutions (1:1, 1:20, 1:400, 1:8,000 and 1:160,000) of Ross River virus (RRV) and Umatilla virus (UMAV) isolates were spiked into subsamples of a pool of 100 Culex australicus mosquitoes. The 1:1 dilution represented the viral load of one RRV-infected mosquito in a pool of 100 mosquitoes. The subsamples underwent nucleic acid extraction, mosquito-specific ribosomal RNA depletion, and Illumina HiSeq sequencing. The viral load of the subsamples was also measured using reverse transcription droplet digital PCR (RT-ddPCR) and quantitative PCR (RT-qPCR). Metatranscriptomic sequencing detected both RRV and UMAV in the 1:1, 1:20 and 1:400 subsamples. A high specificity was achieved, with 100% of RRV and 99.6% of UMAV assembled contigs correctly identified. Metatranscriptomic sequencing was not as sensitive as RT-qPCR or RT-ddPCR; however, it recovered whole genome information and detected 19 other viruses, including four first detections for Australia. These findings will assist arbovirus surveillance programs in utilising metatranscriptomics in routine surveillance activities to enhance arbovirus detection.

Identification of new viral variants specific to the honey bee mite Varroa destructor

Large-scale colony losses among managed Western honey bees have become a serious threat to the beekeeping industry in the last decade. Multiple factors contribute to these losses, but the impact of Varroa destructor parasitism is by far the most important, along with the contribution of some pathogenic viruses vectored by the mite. So far, more than 20 viruses have been identified infecting the honey bee, most of them RNA viruses. They may be maintained either as covert infections or causing severe symptomatic infections, compromising the viability of the colony. In silico analysis of available transcriptomic data obtained from mites collected in the USA and Europe, as well as additional investigation with new samples collected locally, allowed the description of three RNA viruses, two of them variants of the previously described VDV-2 and VDV-3 and the other a new species reported here for the first time. Our results showed that these viruses were widespread among samples and that they were present in the mites as well as in the bees but with differences in the relative abundance and prevalence. However, we have obtained strong evidence showing that these three viruses were able to replicate in the mite, but not in the bee, suggesting that they are selectively infecting the mite. This opens the door to future applications that may help controlling the mite through biological control approaches.

Viral Diversity of Tick Species Parasitizing Cattle and Dogs in Trinidad and Tobago

Ticks are vectors of a wide variety of pathogens that are implicated in mild to severe disease in humans and other animals. Nonetheless, the full range of tick-borne pathogens is unknown. Viruses, in particular, have been neglected in discovery efforts targeting tick-borne agents. High throughput sequencing was used to characterize the virome of 638 ticks, including Rhipicephalus microplus (n = 320), Rhipicephalus sanguineus (n = 300), and Amblyomma ovale (n = 18) collected throughout Trinidad and Tobago in 2017 and 2018. Sequences representing nine viruses were identified, including five novel species within Tymovirales, Bunyavirales, Chuviridae, Rhabdoviridae, and Flaviviridae. Thereafter the frequency of detection of viral sequences in individual tick species was investigated.

Identification of diverse arthropod associated viruses in native Australian fleas

Fleas are important vectors of zoonotic disease. However, little is known about the natural diversity and abundance of flea viruses, particularly in the absence of disease associations, nor the evolutionary relationships among those viruses found in different parasitic vector species. Herein, we present the first virome scale study of fleas, based on the meta-transcriptomic analysis of 52 fleas collected along the eastern coast of Australia. Our analysis revealed 18 novel RNA viruses belonging to nine viral families with diverse genome organizations, although the majority (72%) possessed single-stranded positive-sense genomes. Notably, a number of the viruses identified belonged to the same phylogenetic groups as those observed in ticks sampled at the same locations, although none were likely associated with mammalian infection. Overall, we identified high levels of genomic diversity and abundance of viruses in the flea species studied, and established that fleas harbor viruses similar to those seen to other vectors.

The Ecology of New Constituents of the Tick Virome and Their Relevance to Public Health

Ticks are vectors of several pathogens that can be transmitted to humans and their geographic ranges are expanding. The exposure of ticks to new hosts in a rapidly changing environment is likely to further increase the prevalence and diversity of tick-borne diseases. Although ticks are known to transmit bacteria and viruses, most studies of tick-borne disease have focused upon Lyme disease, which is caused by infection with Borrelia burgdorferi. Until recently, ticks were considered as the vectors of a few viruses that can infect humans and animals, such as Powassan, Tick-Borne Encephalitis and Crimean–Congo hemorrhagic fever viruses. Interestingly, however, several new studies undertaken to reveal the etiology of unknown human febrile illnesses, or to describe the virome of ticks collected in different countries, have uncovered a plethora of novel viruses in ticks. Here, we compared the virome compositions of ticks from different countries and our analysis indicates that the global tick virome is dominated by RNA viruses. Comparative phylogenetic analyses of tick viruses from these different countries reveals distinct geographical clustering of the new tick viruses. Some of these new tick RNA viruses (notably severe fever with thrombocytopenia syndrome virus and Heartland virus) were found to be associated with serious human diseases. Their relevance to public health remains unknown. It is plausible that most of these newly identified tick viruses are of endogenous origin or are restricted in their transmission potential, but the efforts to identify new tick viruses should continue. Indeed, future research aimed at defining the origin, the ecology and the spillover potential of this novel viral biodiversity will be critical to understand the relevance to public health.

Extensive Diversity of RNA Viruses in Australian Ticks

Each year a growing number of individuals along the east coast of Australia experience debilitating disease following tick bites. As there is no evidence for the presence of the causative agent of Lyme disease, Borrelia burgdorferisensu lato, in Australian ticks, the etiological basis of this disease syndrome remains controversial. To characterize the viruses associated with Australian ticks, particularly those that might be associated with mammalian infection, we performed unbiased RNA sequencing on 146 ticks collected across two locations along the coast of New South Wales, Australia. This revealed 19 novel RNA viruses from a diverse set of families. Notably, three of these viruses clustered with known mammalian viruses, including a novel coltivirus that was related to the human pathogen Colorado tick fever virus.

Understanding the microbiome of ticks in Australia is of considerable interest given the ongoing debate over whether Lyme disease and its causative agent, the bacterium Borrelia burgdorferisensu lato, are present in Australia. The diversity of bacteria infecting Australian ticks has been studied using both culture- and metagenomics-based techniques. However, little is known about the virome of Australian ticks, including whether this includes viruses with the potential to infect mammals. We used a meta-transcriptomics approach to reveal the diversity and evolution of viruses from Australian ticks collected from two locations on the central east coast of Australia, including metropolitan Sydney. From this we identified 19 novel RNA viruses belonging to 12 families, as well as 1 previously described RNA virus. The majority of these viruses were related to arthropod-associated viruses, suggesting that they do not utilize mammalian hosts. However, two novel viruses discovered in ticks feeding on bandicoot marsupials clustered closely within the mammal-associated hepacivirus and pestivirus groups (family Flaviviridae). Another bandicoot tick yielded a novel coltivirus (family Reoviridae), a group of largely tick-associated viruses containing the known human pathogen Colorado tick fever virus and its relative, Eyach virus. Importantly, our transcriptomic data provided no evidence for the presence of B. burgdorferisensu lato in any tick sample, providing further evidence against the presence of Lyme disease in Australia. In sum, this study reveals that Australian ticks harbor a diverse virome, including some viruses that merit additional screening in the context of emerging infectious disease.

IMPORTANCE Each year a growing number of individuals along the east coast of Australia experience debilitating disease following tick bites. As there is no evidence for the presence of the causative agent of Lyme disease, Borrelia burgdorferisensu lato, in Australian ticks, the etiological basis of this disease syndrome remains controversial. To characterize the viruses associated with Australian ticks, particularly those that might be associated with mammalian infection, we performed unbiased RNA sequencing on 146 ticks collected across two locations along the coast of New South Wales, Australia. This revealed 19 novel RNA viruses from a diverse set of families. Notably, three of these viruses clustered with known mammalian viruses, including a novel coltivirus that was related to the human pathogen Colorado tick fever virus.

Mapping Arbovirus-Vector Interactions Using Systems Biology Techniques

Studying how arthropod-borne viruses interact with their arthropod vectors is critical to understanding how these viruses replicate and are transmitted. Until recently, these types of studies were limited in scale because of the lack of classical tools available to study virus-host interaction for non-model viruses and non-model organisms. Advances in systems biology "-omics"-based techniques such as next-generation sequencing (NGS) and mass spectrometry can rapidly provide an unbiased view of arbovirus-vector interaction landscapes. In this mini-review, we discuss how arbovirus-vector interaction studies have been advanced by systems biology. We review studies of arbovirus-vector interactions that occur at multiple time and length scales, including intracellular interactions, interactions at the level of the organism, viral and vector populations, and how new techniques can integrate systems-level data across these different scales.