Unbiased Characterization of the Microbiome and Virome of Questing Ticks

Discovery of Novel Viruses in Culicoides Biting Midges in Chihuahua, Mexico

Biting midges (Culicoides) are vectors of many pathogens of medical and veterinary importance, but their viromes are poorly characterized compared to certain other hematophagous arthropods, e.g., mosquitoes and ticks. The goal of this study was to use metagenomics to identify viruses in Culicoides from Mexico. A total of 457 adult midges were collected in Chihuahua, northern Mexico, in 2020 and 2021, and all were identified as female Culicoides reevesi. The midges were sorted into five pools and homogenized. An aliquot of each homogenate was subjected to polyethylene glycol precipitation to enrich for virions, then total RNA was extracted and analyzed by unbiased high-throughput sequencing. We identified six novel viruses that are characteristic of viruses from five families (Nodaviridae, Partitiviridae, Solemoviridae, Tombusviridae, and Totiviridae) and one novel virus that is too divergent from all classified viruses to be assigned to an established family. The newly discovered viruses are phylogenetically distinct from their closest known relatives, and their minimal infection rates in female C. reevesi range from 0.22 to 1.09. No previously known viruses were detected, presumably because viral metagenomics had never before been used to study Culicoides from the Western Hemisphere. To conclude, we discovered multiple novel viruses in C. reevesi from Mexico, expanding our knowledge of arthropod viral diversity and evolution.

Identification of a novel member of the genus Laulavirus (family Phenuiviridae) from the entomopathogenic ascomycete fungus Cordyceps javanica

The virus family Phenuiviridae (order Hareavirales, comprising segmented negative-sense single stranded RNA viruses) has highly diverse members that are known to infect animals, plants, protozoans, and fungi. In this study, we identified a novel phenuivirus infecting a strain of the entomopathogenic fungus Cordyceps javanica isolated from a small brown plant hopper (Laodelphax striatellus), and this virus was tentatively named "Cordyceps javanica negative-strand RNA virus 1" (CjNRSV1). The CjNRSV1 genome consists of three negative-sense single stranded RNA segments (RNA1-3) with lengths of 7252, 2401, and 1117 nt, respectively. The 3'- and 5'-terminal regions of the RNA1, 2, and 3 segments have identical sequences, and the termini of the RNA segments are complementary to each other, reflecting a common characteristic of viruses in the order Hareavirales. RNA1 encodes a large protein (∼274 kDa) containing a conserved domain for the bunyavirus RNA-dependent RNA polymerase (RdRP) superfamily, with 57-80% identity to the RdRP encoded by phenuiviruses in the genus Laulavirus. RNA2 encodes a protein (∼79 kDa) showing sequence similarity (47-63% identity) to the movement protein (MP, a plant viral cell-to-cell movement protein)-like protein (MP-L) encoded by RNA2 of laulaviruses. RNA3 encodes a protein (∼28 kDa) with a conserved domain of the phenuivirid nucleocapsid protein superfamily. Phylogenetic analysis using the RdRPs of various phenuiviruses and other unclassified phenuiviruses showed CjNRSV1 to be grouped with established members of the genus Laulavirus. Our results suggest that CjNRSV1 is a novel fungus-infecting member of the genus Laulavirus in the family Phenuiviridae.

Metatranscriptomic investigation of single Ixodes pacificus ticks reveals diverse microbes, viruses, and novel mRNA-like endogenous viral elements

Ticks are increasingly important vectors of human and agricultural diseases. While many studies have focused on tick-borne bacteria, far less is known about tick-associated viruses and their roles in public health or tick physiology. To address this, we investigated patterns of bacterial and viral communities across two field populations of western black-legged ticks (Ixodes pacificus). Through metatranscriptomic analysis of 100 individual ticks, we quantified taxon prevalence, abundance, and co-occurrence with other members of the tick microbiome. In addition to commonly found tick-associated microbes, we assembled 11 novel RNA virus genomes from Rhabdoviridae, Chuviridae, Picornaviridae, Phenuiviridae, Reoviridae, Solemovidiae, Narnaviridae and two highly divergent RNA virus genomes lacking sequence similarity to any known viral families. We experimentally verified the presence of these in I. pacificus ticks across several life stages. We also unexpectedly identified numerous virus-like transcripts that are likely encoded by tick genomic DNA, and which are distinct from known endogenous viral element-mediated immunity pathways in invertebrates. Taken together, our work reveals that I. pacificus ticks carry a greater diversity of viruses than previously appreciated, in some cases resulting in evolutionarily acquired virus-like transcripts. Our findings highlight how pervasive and intimate tick-virus interactions are, with major implications for both the fundamental biology and vectorial capacity of I. pacificus ticks.

IMPORTANCE

Ticks are increasingly important vectors of disease, particularly in the United States where expanding tick ranges and intrusion into previously wild areas has resulted in increasing human exposure to ticks. Emerging human pathogens have been identified in ticks at an increasing rate, and yet little is known about the full community of microbes circulating in various tick species, a crucial first step to understanding how they interact with each and their tick host, as well as their ability to cause disease in humans. We investigated the bacterial and viral communities of the Western blacklegged tick in California and found 11 previously uncharacterized viruses circulating in this population.

Development and validation of a long-read metabarcoding platform for the detection of filarial worm pathogens of animals and humans

BACKGROUND

Filarial worms are important vector-borne pathogens of a large range of animal hosts, including humans, and are responsible for numerous debilitating neglected tropical diseases such as, lymphatic filariasis caused by Wuchereria bancrofti and Brugia spp., as well as loiasis caused by Loa loa. Moreover, some emerging or difficult-to-eliminate filarioid pathogens are zoonotic using animals like canines as reservoir hosts, for example Dirofilaria sp. 'hongkongensis'. Diagnosis of filariasis through commonly available methods, like microscopy, can be challenging as microfilaremia may wane below the limit of detection. In contrast, conventional PCR methods are more sensitive and specific but may show limited ability to detect coinfections as well as emerging and/or novel pathogens. Use of deep-sequencing technologies obviate these challenges, providing sensitive detection of entire parasite communities, whilst also being better suited for the characterisation of rare or novel pathogens. Therefore, we developed a novel long-read metabarcoding assay for deep-sequencing the filarial nematode cytochrome c oxidase subunit I gene on Oxford Nanopore Technologies' (ONT) MinION™ sequencer. We assessed the overall performance of our assay using kappa statistics to compare it to commonly used diagnostic methods for filarial worm detection, such as conventional PCR (cPCR) with Sanger sequencing and the microscopy-based modified Knott's test (MKT).

RESULTS

We confirmed our metabarcoding assay can characterise filarial parasites from a diverse range of genera, including, Breinlia, Brugia, Cercopithifilaria, Dipetalonema, Dirofilaria, Onchocerca, Setaria, Stephanofilaria and Wuchereria. We demonstrated proof-of-concept for this assay by using blood samples from Sri Lankan dogs, whereby we identified infections with the filarioids Acanthocheilonema reconditum, Brugia sp. Sri Lanka genotype and zoonotic Dirofilaria sp. 'hongkongensis'. When compared to traditionally used diagnostics, such as the MKT and cPCR with Sanger sequencing, we identified an additional filarioid species and over 15% more mono- and coinfections.

CONCLUSIONS

Our developed metabarcoding assay may show broad applicability for the metabarcoding and diagnosis of the full spectrum of filarioids from a wide range of animal hosts, including mammals and vectors, whilst the utilisation of ONT' small and portable MinION™ means that such methods could be deployed for field use.

Characterisation of putative novel tick viruses and zoonotic risk prediction

Tick-associated viruses remain a substantial zoonotic risk worldwide, so knowledge of the diversity of tick viruses has potential health consequences. Despite their importance, large amounts of sequences in public data sets from tick meta-genomic and -transcriptomic projects remain unannotated, sequence data that could contain undocumented viruses. Through data mining and bioinformatic analysis of more than 37,800 public meta-genomic and -transcriptomic data sets, we found 83 unannotated contigs exhibiting high identity with known tick viruses. These putative viral contigs were classified into three RNA viral families (Alphatetraviridae, Orthomyxoviridae and Chuviridae) and one DNA viral family (Asfarviridae). After manual checking of quality and dissimilarity towards other sequences in the data set, these 83 contigs were reduced to five contigs in the Alphatetraviridae from four putative viruses, four in the Orthomyxoviridae from two putative viruses and one in the Chuviridae which clustered with known tick-associated viruses, forming a separate clade within the viral families. We further attempted to assess which previously known tick viruses likely represent zoonotic risks and thus deserve further investigation. We ranked the human infection potential of 133 known tick-associated viruses using a genome composition-based machine learning model. We found five high-risk tick-associated viruses (Langat virus, Lonestar tick chuvirus 1, Grotenhout virus, Taggert virus and Johnston Atoll virus) that have not been known to infect human and two viral families (Nairoviridae and Phenuiviridae) that contain a large proportion of potential zoonotic tick-associated viruses. This adds to the knowledge of tick virus diversity and highlights the importance of surveillance of newly emerging tick-associated diseases.

Annual (2023) taxonomic update of RNA-directed RNA polymerase-encoding negative-sense RNA viruses (realm Riboviria: kingdom Orthornavirae: phylum Negarnaviricota)

In April 2023, following the annual International Committee on Taxonomy of Viruses (ICTV) ratification vote on newly proposed taxa, the phylum Negarnaviricota was amended and emended. The phylum was expanded by one new family, 14 new genera, and 140 new species. Two genera and 538 species were renamed. One species was moved, and four were abolished. This article presents the updated taxonomy of Negarnaviricota as now accepted by the ICTV.

The Troublesome Ticks Research Protocol: Developing a Comprehensive, Multidiscipline Research Plan for Investigating Human Tick-Associated Disease in Australia

In Australia, there is a paucity of data about the extent and impact of zoonotic tick-related illnesses. Even less is understood about a multifaceted illness referred to as Debilitating Symptom Complexes Attributed to Ticks (DSCATT). Here, we describe a research plan for investigating the aetiology, pathophysiology, and clinical outcomes of human tick-associated disease in Australia. Our approach focuses on the transmission of potential pathogens and the immunological responses of the patient after a tick bite. The protocol is strengthened by prospective data collection, the recruitment of two external matched control groups, and sophisticated integrative data analysis which, collectively, will allow the robust demonstration of associations between a tick bite and the development of clinical and pathological abnormalities. Various laboratory analyses are performed including metagenomics to investigate the potential transmission of bacteria, protozoa and/or viruses during tick bite. In addition, multi-omics technology is applied to investigate links between host immune responses and potential infectious and non-infectious disease causations. Psychometric profiling is also used to investigate whether psychological attributes influence symptom development. This research will fill important knowledge gaps about tick-borne diseases. Ultimately, we hope the results will promote improved diagnostic outcomes, and inform the safe management and treatment of patients bitten by ticks in Australia.

A global dataset of microbial community in ticks from metagenome study

Ticks are important vectors of various zoonotic pathogens that can infect animals and humans, and most documented tick-borne pathogens have a strong bias towards microorganisms with strong disease phenotypes. The recent development of next-generation sequencing (NGS) has enabled the study of microbial communities, referred to as microbiome. Herein, we undertake a systematic review of published literature to build a comprehensive global dataset of microbiome determined by NGS in field-collected ticks. The dataset comprised 4418 records from 76 literature involving geo-referenced occurrences for 46 species of ticks and 219 microorganism families, revealing a total of 83 emerging viruses identified from 24 tick species belonging to 6 tick genera since 1980. The viral, bacterial and eukaryotic composition was compared regarding the tick species, their live stage and types of the specimens, or the geographic location. The data can assist the further investigation of ecological, biogeographical and epidemiological features of the tick-borne disease.

Measurement(s)	microbial community
Technology Type(s)	Next Generation Sequencing
Factor Type(s)	tick
Sample Characteristic - Organism	tick
Sample Characteristic - Environment	microbial community
Sample Characteristic - Location	Whole world

Metatranscriptomic profiling reveals diverse tick-borne bacteria, protozoans and viruses in ticks and wildlife from Australia

Tick-borne zoonoses are emerging globally due to changes in climate and land use. While the zoonotic threats associated with ticks are well studied elsewhere, in Australia, the diversity of potentially zoonotic agents carried by ticks and their significance to human and animal health is not sufficiently understood. To this end, we used untargeted metatranscriptomics to audit the prokaryotic, eukaryotic and viral biomes of questing ticks and wildlife blood samples from two urban and rural sites in New South Wales, Australia. Ixodes holocyclus and Haemaphysalis bancrofti were the main tick species collected, and blood samples from Rattus rattus, Rattus fuscipes, Perameles nasuta and Trichosurus vulpecula were also collected and screened for tick-borne microorganisms using metatranscriptomics followed by conventional targeted PCR to identify important microbial taxa to the species level. Our analyses identified 32 unique tick-borne taxa, including 10 novel putative species. Overall, a wide range of tick-borne microorganisms were found in questing ticks including haemoprotozoa such as Babesia, Theileria, Hepatozoon and Trypanosoma spp., bacteria such as Borrelia, Rickettsia, Ehrlichia, Neoehrlichia and Anaplasma spp., and numerous viral taxa including Reoviridiae (including two coltiviruses) and a novel Flaviviridae-like jingmenvirus. Of note, a novel hard tick-borne relapsing fever Borrelia sp. was identified in questing H. bancrofti ticks which is closely related to, but distinct from, cervid-associated Borrelia spp. found throughout Asia. Notably, all tick-borne microorganisms were phylogenetically unique compared to their relatives found outside Australia, and no foreign tick-borne human pathogens such as Borrelia burgdorferi s.l. or Babesia microti were found. This work adds to the growing literature demonstrating that Australian ticks harbour a unique and endemic microbial fauna, including potentially zoonotic agents which should be further studied to determine their relative risk to human and animal health.

Extended characterisation of five archival tick-borne viruses provides insights for virus discovery in Australian ticks

Background

A subset of Australians who have been bitten by ticks experience a complex of chronic and debilitating symptoms which cannot be attributed to the known pathogenic species of bacteria present in Australia. As a result, there has been a renewed effort to identify and characterise viruses in Australian terrestrial ticks. Recent transcriptome sequencing of Ixodes and Amblyomma ticks has revealed the presence of multiple virus sequences. However, without virus isolates our ability to understand the host range and pathogenesis of newly identified viruses is limited. We have established a successful method for high-throughput virus discovery and isolation in mosquitoes using antibodies to double-stranded RNA. In this study we sought to characterise five archival tick-borne viruses to adapt our virus discovery protocol for Australian ticks.

Methods

We performed virus characterisation using a combination of bioinformatic sequence analysis and in vitro techniques including replication kinetics, antigenic profiling, virus purification and mass spectrometry.

Results

Our sequence analysis of Nugget virus, Catch-me-Cave virus and Finch Creek virus revealed marked genetic stability in isolates collected from the same location approximately 30 years apart. We demonstrate that the Ixodes scapularis-derived ISE6 cell line supports replication of Australian members of the Flaviviridae, Nairoviridae, Phenuiviridae and Reoviridae families, including Saumarez Reef virus (SREV), a flavivirus isolated from the soft tick Ornithodoros capensis. While antibodies against double-stranded RNA could be used to detect replication of a tick-borne reovirus and mosquito-borne flavivirus, the tick-borne flaviviruses Gadgets Gully virus and SREV could not be detected using this method. Finally, four novel virus-like sequences were identified in transcriptome sequencing of the Australian native tick Ixodes holocyclus.

Conclusions

Genetic and antigenic characterisations of archival viruses in this study confirm that three viruses described in 2002 represent contemporary isolates of virus species first identified 30 years prior. Our findings with antibodies to double-stranded RNA highlight an unusual characteristic shared by two Australian tick-borne flaviviruses. Finally, comparative growth kinetics analyses of Australian tick-borne members of the Flaviviridae, Nairoviridae, Phenuiviridae and Reoviridae families in ISE6 and BSR cells will provide a useful resource for isolation of Australian tick-borne viruses using existing cell lines.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05176-z.