Metavirome of 31 tick species provides a compendium of 1,801 RNA virus genomes

Differential susceptibility of geographically distinct Ixodes ricinus populations to tick-borne encephalitis virus and louping ill virus

Tick-borne encephalitis virus (TBEV) is an emerging pathogen in the Netherlands. Multiple divergent viral strains are circulating and the focal distribution of TBEV remains poorly understood. This may, however, be explained by differences in the susceptibility of tick populations for specific viruses and viral strains, and by viral strains having higher infection success in their local tick population. We investigated this hypothesis by exposing Dutch Ixodes ricinus ticks to two different TBEV strains: TBEV-NL from the Netherlands and TBEV-Neudoerfl from Austria. In addition, we exposed ticks to louping Ill virus (LIV), which is endemic to large parts of the United Kingdom and Ireland, but has not been reported in the Netherlands. Ticks were collected from two locations in the Netherlands: one location without evidence of TBEV circulation and one location endemic for the TBEV-NL strain. Ticks were infected in a biosafety level 3 laboratory using an artificial membrane feeding system. Ticks collected from the region without evidence of TBEV circulation had lower infection rates for TBEV-NL as compared to TBEV-Neudoerfl. Vice versa, ticks collected from the TBEV-NL endemic region had higher infection rates for TBEV-NL compared to TBEV-Neudoerfl. In addition, LIV infection rates were much lower in Dutch ticks compared to TBEV, which may explain why LIV is not present in the Netherlands. Our findings show that ticks from two distinct geographical populations differ in their susceptibility to TBEV strains, which could be the result of differences in the genetic background of the tick populations.

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.

Molecular characterization and induced changes of histone acetyltransferases in the tick Haemaphysalis longicornis in response to cold stress

BACKGROUND

Epigenetic modifications of histones play important roles in the response of eukaryotic organisms to environmental stress. However, many histone acetyltransferases (HATs), which are responsible for histone acetylation, and their roles in mediating the tick response to cold stress have yet to be identified. In the present study, HATs were molecularly characterized and their associations with the cold response of the tick Haemaphysalis longicornis explored.

METHODS

HATs were characterized by using polymerase chain reaction (PCR) based on published genome sequences, followed by multiple bioinformatic analyses. The differential expression of genes in H. longicornis under different cold treatment conditions was evaluated using reverse transcription quantitative PCR (RT-qPCR). RNA interference was used to explore the association of HATs with the cold response of H. longicornis.

RESULTS

Two HAT genes were identified in H. longicornis (Hl), a GCN5-related N-acetyltransferase (henceforth HlGNAT) and a type B histone acetyltransferase (henceforth HlHAT-B), which are respectively 960 base pairs (bp) and 1239 bp in length. Bioinformatics analysis revealed that HlGNAT and HlHAT-B are unstable hydrophilic proteins characterized by the presence of the acetyltransferase 16 domain and Hat1_N domain, respectively. RT-qPCR revealed that the expression of HlGNAT and HlHAT-B decreased after 3 days of cold treatment, but gradually increased with a longer period of cold treatment. The mortality rate following knockdown of HlGNAT or HlHAT-B by RNA interference, which was confirmed by RT-qPCR, significantly increased (P < 0.05) when H. longicornis was treated at the lowest lethal temperature (- 14 °C) for 2 h.

CONCLUSIONS

The findings demonstrate that HATs may play a crucial role in the cold response of H. longicornis. Thus further research is warranted to explore the mechanisms underlying the epigenetic regulation of the cold response in ticks.

Metagenomic analysis of individual mosquito viromes reveals the geographical patterns and drivers of viral diversity

Mosquito transmitted viruses are responsible for an increasing burden of human disease. Despite this, little is known about the diversity and ecology of viruses within individual mosquito hosts. Here, using a meta-transcriptomic approach, we determined the viromes of 2,438 individual mosquitoes (81 species), spanning ~4,000 km along latitudes and longitudes in China. From these data we identified 393 viral species associated with mosquitoes, including 7 (putative) species of arthropod-borne viruses (that is, arboviruses). We identified potential mosquito species and geographic hotspots of viral diversity and arbovirus occurrence, and demonstrated that the composition of individual mosquito viromes was strongly associated with host phylogeny. Our data revealed a large number of viruses shared among mosquito species or genera, enhancing our understanding of the host specificity of insect-associated viruses. We also detected multiple virus species that were widespread throughout the country, perhaps reflecting long-distance mosquito dispersal. Together, these results greatly expand the known mosquito virome, linked viral diversity at the scale of individual insects to that at a country-wide scale, and offered unique insights into the biogeography and diversity of viruses in insect vectors.

Metagenomic sequencing reveals viral diversity of mosquitoes from Shandong Province, China

Mosquitoes carry a large number of known and unknown viruses, some of which could cause serious diseases in humans or animals. Metagenomic sequencing for mosquito viromes is crucial for understanding the evolutionary history of viruses and preventing emerging mosquito-borne diseases. We collected 1,598 mosquitoes belonging to four species from five counties in Shandong Province, China in 2021. They were grouped by species and sampling locations and subjected to metagenomic next-generation sequencing for the analysis of the viromes. A total of 233,317,352 sequencing reads were classified into 30 viral families and an unclassified group. Comparative analysis showed that mosquitoes in Shandong Province generally possessed host-specific virome. We detected mosquito-borne viruses including Japanese encephalitis virus, Getah virus, and Kadipiro virus in Culex tritaeniorhynchus and Anopheles sinensis samples. Phylogenetic analysis showed that these pathogenic viruses may have existed in mosquitoes in Shandong Province for a long time. Meanwhile, we identified 22 novel viruses belonging to seven families and the genus Negevirus. Our study comprehensively described the viromes of several common mosquito species in Shandong Province, China, and demonstrated the major role of host species in shaping mosquito viromes. Furthermore, the metagenomic data provided valuable epidemiological information on multiple mosquito-borne viruses, highlighting the potential risk of infection transmission.

IMPORTANCE

Mosquitoes are known as the source of various pathogens for humans and animals. Culex tritaeniorhynchus, Armigeres subalbatus, and Anopheles sinensis have been found to transmit the Getah virus, which has recently caused increasing infections in China. Cx. tritaeniorhynchus and Culex pipiens are the main vectors of Japanese encephalitis virus and have caused epidemics of Japanese encephalitis in China in past decades. These mosquitoes are widely present in Shandong Province, China, leading to a great threat to public health and the breeding industry. This study provided a comprehensive insight into the viromes of several common mosquito species in Shandong Province, China. The metagenomic sequencing data revealed the risks of multiple pathogenic mosquito-borne viruses, including Japanese encephalitis virus, Getah virus, and Kadipiro virus, which are of great importance for preventing emerging viral epidemics.

Viromes of Haemaphysalis longicornis reveal different viral abundance and diversity in free and engorged ticks

Haemaphysalis longicornis ticks, commonly found in East Asia, can transmit various pathogenic viruses, including the severe fever with thrombocytopenia syndrome virus (SFTSV) that has caused febrile diseases among humans in Hubei Province. However, understanding of the viromes of H. longicornis was limited, and the prevalence of viruses among H. longicornis ticks in Hubei was not well clarified. This study investigates the viromes of both engorged (fed) and free (unfed) H. longicornis ticks across three mountainous regions in Hubei Province from 2019 to 2020. RNA-sequencing analysis identified viral sequences that were related to 39 reference viruses belonging to unclassified viruses and seven RNA viral families, namely Chuviridae, Nairoviridae, Orthomyxoviridae, Parvoviridae, Phenuiviridae, Rhabdoviridae, and Totiviridae. Viral abundance and diversity in these ticks were analysed, and phylogenetic characteristics of the Henan tick virus (HNTV), Dabieshan tick virus (DBSTV), Okutama tick virus (OKTV), and Jingmen tick virus (JMTV) were elucidated based on their full genomic sequences. Prevalence analysis demonstrated that DBSTV was the most common virus found in individual H. longicornis ticks (12.59%), followed by HNTV (0.35%), whereas JMTV and OKTV were not detected. These results improve our understanding of H. longicornis tick viromes in central China and highlight the role of tick feeding status and geography in shaping the viral community. The findings of new viral strains and their potential impact on public health raise the need to strengthen surveillance efforts for comprehensively assessing their spillover potentials.

Virome diversity shaped by genetic evolution and ecological landscape of Haemaphysalis longicornis

BACKGROUND

Haemaphysalis longicornis is drawing attentions for its geographic invasion, extending population, and emerging disease threat. However, there are still substantial gaps in our knowledge of viral composition in relation to genetic diversity of H. longicornis and ecological factors, which are important for us to understand interactions between virus and vector, as well as between vector and ecological elements.

RESULTS

We conducted the meta-transcriptomic sequencing of 136 pools of H. longicornis and identified 508 RNA viruses of 48 viral species, 22 of which have never been reported. Phylogenetic analysis of mitochondrion sequences divided the ticks into two genetic clades, each of which was geographically clustered and significantly associated with ecological factors, including altitude, precipitation, and normalized difference vegetation index. The two clades showed significant difference in virome diversity and shared about one fifth number of viral species that might have evolved to "generalists." Notably, Bandavirus dabieense, the pathogen of severe fever with thrombocytopenia syndrome was only detected in ticks of clade 1, and half number of clade 2-specific viruses were aquatic-animal-associated.

CONCLUSIONS

These findings highlight that the virome diversity is shaped by internal genetic evolution and external ecological landscape of H. longicornis and provide the new foundation for promoting the studies on virus-vector-ecology interaction and eventually for evaluating the risk of H. longicornis for transmitting the viruses to humans and animals. Video Abstract.

Ecoclimate drivers shape virome diversity in a globally invasive tick species

Spillovers of viruses from animals to humans occur more frequently under warmer conditions, particularly arboviruses. The invasive tick species Haemaphysalis longicornis, the Asian longhorned tick, poses a significant public health threat due to its global expansion and its potential to carry a wide range of pathogens. We analyzed meta-transcriptomic data from 3595 adult H. longicornis ticks collected between 2016 and 2019 in 22 provinces across China encompassing diverse ecological conditions. Generalized additive modeling revealed that climate factors exerted a stronger influence on the virome of H. longicornis than other ecological factors, such as ecotypes, distance to coastline, animal host, tick gender, and antiviral immunity. To understand how climate changes drive the tick virome, we performed a mechanistic investigation using causality inference with emphasis on the significance of this process for public health. Our findings demonstrated that higher temperatures and lower relative humidity/precipitation contribute to variations in animal host diversity, leading to increased diversity of the tick virome, particularly the evenness of vertebrate-associated viruses. These findings may explain the evolution of tick-borne viruses into generalists across multiple hosts, thereby increasing the probability of spillover events involving tick-borne pathogens. Deep learning projections have indicated that the diversity of the H. longicornis virome is expected to increase in 81.9% of regions under the SSP8.5 scenario from 2019 to 2030. Extension of surveillance should be implemented to avert the spread of tick-borne diseases.

Genetic diversity of RNA viruses infecting invertebrate pests of rice

Invertebrate species are a natural reservoir of viral genetic diversity, and invertebrate pests are widely distributed in crop fields. However, information on viruses infecting invertebrate pests of crops is limited. In this report, we describe the deep metatranscriptomic sequencing of 88 invertebrate samples covering all major invertebrate pests in rice fields. We identified 296 new RNA viruses and 13 known RNA viruses. These viruses clustered within 31 families, with many highly divergent viruses constituting potentially new families and genera. Of the identified viruses, 13 RNA viruses clustered within the Fiersviridae family of bacteriophages, and 48 RNA viruses clustered within families and genera of mycoviruses. We detected known rice viruses in novel invertebrate hosts at high abundances. Furthermore, some novel RNA viruses have genome structures closely matching to known plant viruses and clustered within genera of several plant virus species. Forty-five potential insect pathogenic RNA viruses were detected in invertebrate species. Our analysis revealed that host taxonomy plays a major role and geographical location plays an important role in structuring viral diversity. Cross-species transmission of RNA viruses was detected between invertebrate hosts. Newly identified viral genomes showed extensive variation for invertebrate viral families or genera. Together, the large-scale metatranscriptomic analysis greatly expands our understanding of RNA viruses in rice invertebrate species, the results provide valuable information for developing efficient strategies to manage insect pests and virus-mediated crop diseases.

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.

Metatranscriptomics Reveals the RNA Virome of Ixodes Persulcatus in the China-North Korea Border, 2017

In recent years, numerous viruses have been identified from ticks, and some have been linked to clinical cases of emerging tick-borne diseases. Chinese northeast frontier is tick infested. However, there is a notable lack of systematic monitoring efforts to assess the viral composition in the area, leaving the ecological landscape of viruses carried by ticks not clear enough. Between April and June 2017, 7101 ticks were collected to perform virus surveillance on the China-North Korea border, specifically in Tonghua, Baishan, and Yanbian. A total of 2127 Ixodes persulcatus were identified. Further investigation revealed the diversity of tick-borne viruses by transcriptome sequencing of Ixodes persulcatus. All ticks tested negative for tick-borne encephalitis virus. Transcriptome sequencing expanded 121 genomic sequence data of 12 different virus species from Ixodes persulcatus. Notably, a new segmented flavivirus, named Baishan Forest Tick Virus, were identified, closely related to Alongshan virus and Harz mountain virus. Therefore, this new virus may pose a potential threat to humans. Furthermore, the study revealed the existence of seven emerging tick-borne viruses dating back to 2017. These previously identified viruses included Mudanjiang phlebovirus, Onega tick phlebovirus, Sara tick phlebovirus, Yichun mivirus, and three unnamed viruses (one belonging to the Peribunyaviridae family and the other two belonging to the Phenuiviridae family). The existence of these emerging tick-borne viruses in tick samples collected in 2017 suggests that their history may extend further than previously recognized. This study provides invaluable insights into the virome of Ixodes persulcatus in the China-North Korea border region, enhancing our ongoing efforts to manage the risks associated with tick-borne viruses.

RNA virus diversity in three parasitoid wasps of tephritid flies: insights from novel and known species

IMPORTANCE

Parasitoid wasp populations have developed persistent beneficial symbiotic relationships with several viruses through repeated evolution. However, there have been limited reports on RNA viruses in parasitoid wasps of tephritid flies, a significant pest group affecting fruits and vegetables. This study explores the diversity of RNA viruses in three parasitoid wasps of tephritid flies and highlights the potential biological significance of specific viruses in Diachasmimorpha longicaudata. These findings have important implications for the development of sustainable pest management strategies and the enhancement of artificial rearing techniques for parasitoid wasps.

Hazard potential of Swiss Ixodes ricinus ticks: Virome composition and presence of selected bacterial and protozoan pathogens

Ticks play an important role in transmitting many different emerging zoonotic pathogens that pose a significant threat to human and animal health. In Switzerland and abroad, the number of tick-borne diseases, in particular tick-borne encephalitis (TBE), has been increasing over the last few years. Thus, it remains essential to investigate the pathogen spectrum of ticks to rapidly detect emerging pathogens and initiate the necessary measures. To assess the risk of tick-borne diseases in different regions of Switzerland, we collected a total of 10'286 ticks from rural and urban areas in ten cantons in 2021 and 2022. Ticks were pooled according to species, developmental stage, gender, and collection site, and analyzed using next generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR). The metagenomic analysis revealed for the first time the presence of Alongshan virus (ALSV) in Swiss ticks. Interestingly, the pool-prevalence of ALSV was higher than that of tick-borne encephalitis virus (TBEV). Furthermore, several TBEV foci have been identified and pool prevalence of selected non-viral pathogens determined.

Metagenomic analysis of individual mosquitos reveals the ecology of insect viruses

Mosquito transmitted viruses are responsible for an increasing burden of human disease. Despite this, little is known about the diversity and ecology of viruses within individual mosquito hosts. Using a meta-transcriptomic approach, we analysed the virome of 2,438 individual mosquitos (79 species), spanning ~4000 km along latitudes and longitudes in China. From these data we identified 393 core viral species associated with mosquitos, including seven (putative) arbovirus species. We identified potential species and geographic hotspots of viral richness and arbovirus occurrence, and demonstrated that host phylogeny had a strong impact on the composition of individual mosquito viromes. Our data revealed a large number of viruses shared among mosquito species or genera, expanding our knowledge of host specificity of insect-associated viruses. We also detected multiple virus species that were widespread throughout the country, possibly facilitated by long-distance mosquito migrations. Together, our results greatly expand the known mosquito virome, linked the viral diversity at the scale of individual insects to that at a country-wide scale, and offered unique insights into the ecology of viruses of insect vectors.

Virome analysis of an ectomycorrhizal fungus Suillus luteus revealing potential evolutionary implications

Suillus luteus is a widespread edible ectomycorrhizal fungus that holds significant importance in both ecological and economic value. Mycoviruses are ubiquitous infectious agents hosted in different fungi, with some known to exert beneficial or detrimental effects on their hosts. However, mycoviruses hosted in ectomycorrhizal fungi remain poorly studied. To address this gap in knowledge, we employed next-generation sequencing (NGS) to investigate the virome of S. luteus. Using BLASTp analysis and phylogenetic tree construction, we identified 33 mycovirus species, with over half of them belonging to the phylum Lenarviricota, and 29 of these viruses were novel. These mycoviruses were further grouped into 11 lineages, with the discovery of a new negative-sense single-stranded RNA viral family in the order Bunyavirales. In addition, our findings suggest the occurrence of cross-species transmission (CST) between the fungus and ticks, shedding light on potential evolutionary events that have shaped the viral community in different hosts. This study is not only the first study to characterize mycoviruses in S. luteus but highlights the enormous diversity of mycoviruses and their implications for virus evolution.