Nanopore-based metagenomics reveal a new Rickettsia in Europe.

Accurate identification of tick-borne bacteria, including those associated with rickettsioses, pose significant challenges due to the polymicrobial and polyvectoral nature of the infections. We aimed to carry out a comparative evaluation of a non-targeted metagenomic approach by nanopore sequencing (NS) and commonly used PCR assays amplifying Rickettsia genes in field-collected ticks. The study included a total of 310 ticks, originating from Poland (44.2 %) and Bulgaria (55.8 %). Samples comprised 7 species, the majority of which were Ixodes ricinus (62.9 %), followed by Dermacentor reticulatus (21.2 %). Screening was carried out in 55 pools, using total nucleic acid extractions from individual ticks. NS and ompA/gltA PCRs identified Rickettsia species in 47.3 % and 54.5 % of the pools, respectively. The most frequently detected species were Rickettsia asiatica (27.2 %) and Rickettsia raoultii (21.8 %), followed by Rickettsia monacensis (3.6 %), Rickettsia helvetica (1.8 %), Rickettsia massiliae (1.8 %) and Rickettsia tillamookensis (1.8 %). Phylogeny construction on mutS, uvrD, argS and virB4 sequences and a follow-up deep sequencing further supported R. asiatica identification, documented in Europe for the first time. NS further enabled detection of Anaplasma phagocytophilum (9.1 %), Coxiella burnetii (5.4 %) and Neoehrlichia mikurensis (1.8 %), as well as various endosymbionts of Rickettsia and Coxiella. Co-detection of multiple rickettsial and non-rickettsial bacteria were observed in 16.4 % of the pools with chromosome and plasmid-based contigs. In conclusion, non-targeted metagenomic sequencing was documented as a robust strategy capable of providing a broader view of the tick-borne bacterial pathogen spectrum.

Novel clades of tick-borne pathogenic nairoviruses in Europe.

Members of the Orthonairovirus genus (family Nairoviridae) include many tick-borne viruses of significant human and animal health impact, with several recently-documented pathogenic viruses lacking sufficient epidemiological information. We screened 215 adult ticks of seven species collected in Bulgaria, Georgia, Latvia and Poland for orthonairoviruses, followed by nanopore sequencing (NS) for genome characterization. Initial generic amplification revealed Sulina virus (SULV, Orthonairovirus sulinaense), for which an updated amplification assay was used, revealing an overall prevalence of 2.7% in Ixodes ricinus ticks from Latvia. Three complete and additional partial SULV genomes were generated, that consistently formed a separate, distinct clade with further intragroup divergence in the maximum likelihood analyses. Comparisons with previously described viruses from Romania exhibited similar genome topologies, albeit with divergent motifs and cleavage sites on the glycoprotein precursor. Preliminary evidence of recombination involving the S segment was documented, in addition to variations in predicted viral glycoproteins. Generic screening further identified Tacheng tick virus 1 (TCTV1, Orthonairovirus tachengense), with documented human infections, in Dermacentor reticulatus ticks from Poland, with a prevalence of 0.9%. Subsequent NS and assembly provided the first complete TCTV1 genome outside of China, where it was originally described. Phylogenetic analysis of virus genome segments revealed TCTV1-Poland as a discrete taxon within the TCTV1 cluster in the Orthonairovirus genus, representing a geographically segregated clade. Comparable genome topology with TCTV1 from China was observed, aside from minor variations in the M segment. Similar to SULV, TCTV1 exhibited several mismatches on previously described screening primer binding sites, likely to prevent amplification. These findings indicate presence of novel TCTV1 and SULV clades in Eastern Europe, confirming the expansion of orthonairoviruses with pathogenic potential.

Metagenomics analysis reveals presence of the Merida-like virus in Georgia.

Arbovirus surveillance is fundamental for the discovery of novel viruses and prevention of febrile vector-borne illnesses. Vector-borne pathogens can rapidly expand and adapt in new geographic and environmental conditions. In this study, metagenomic surveillance was conducted to identify novel viruses in the Country of Georgia. A total of 521 mosquitoes were captured near a military training facility and pooled from species Culex pipiens (Linnaeus) (87%) and Aedes albopictus (Skuse) (13%). We decided to further analyze the Culex pipiens mosquitoes, due to the more extensive number of samples collected. Our approach was to utilize an unbiased total RNA-seq for pathogen discovery in order to explore the mosquito virome. The viral reads from this analysis were mostly aligned to Insect-specific viruses from two main families, the Iflaviridae; a positive-stranded RNA virus and the Rhabdoviridae; a negative- and single-stranded RNA virus. Our pathogen discovery analysis revealed viral reads aligning to the Merida-like virus Turkey (MERDLVT) strain among the Rhabdoviridae. To further validate this result, we conducted a BLAST sequence comparison analysis of our samples with the MERDLVT strain. Our positive samples aligned to the MERDLVT strain with 96-100% sequence identity and 99.7-100% sequence coverage. A bootstrapped maximum-likelihood phylogenetic tree was used to evaluate the evolutionary relationships among these positive pooled specimens with the (MERDLVT) strain. The Georgia samples clustered most closely with two strains from Turkey, the Merida-like virus KE-2017a isolate 139-1-21 and the Merida-like virus Turkey isolate P431. Collectively, these results show the presence of the MERDLVT strain in Georgia.

Impact of nanopore-based metagenome sequencing on tick-borne virus detection.

Introduction: We evaluated metagenomic nanopore sequencing (NS) in field-collected ticks and compared findings from amplification-based assays. Methods: Forty tick pools collected in Anatolia, Turkey and screened by broad-range or nested polymerase chain reaction (PCR) for Crimean-Congo Hemorrhagic Fever Virus (CCHFV) and Jingmen tick virus (JMTV) were subjected to NS using a standard, cDNA-based metagenome approach. Results: Eleven viruses from seven genera/species were identified. Miviruses Bole tick virus 3 and Xinjiang mivirus 1 were detected in 82.5 and 2.5% of the pools, respectively. Tick phleboviruses were present in 60% of the pools, with four distinct viral variants. JMTV was identified in 60% of the pools, where only 22.5% were PCR-positive. CCHFV sequences characterized as Aigai virus were detected in 50%, where only 15% were detected by PCR. NS produced a statistically significant increase in detection of these viruses. No correlation of total virus, specific virus, or targeted segment read counts was observed between PCR-positive and PCR-negative samples. NS further enabled the initial description of Quaranjavirus sequences in ticks, where human and avian pathogenicity of particular isolates had been previously documented. Discussion: NS was observed to surpass broad-range and nested amplification in detection and to generate sufficient genome-wide data for investigating virus diversity. It can be employed for monitoring pathogens in tick vectors or human/animal clinical samples in hot-spot regions for examining zoonotic spillover.

The expanding range of emerging tick-borne viruses in Eastern Europe and the Black Sea Region.

We analysed both pooled and individual tick samples collected from four countries in Eastern Europe and the Black Sea region, using metagenome-based nanopore sequencing (NS) and targeted amplification. Initially, 1337 ticks, belonging to 11 species, were screened in 217 pools. Viruses (21 taxa) and human pathogens were detected in 46.5% and 7.3%, respectively. Tick-borne viral pathogens comprised Tacheng Tick Virus 2 (TTV2, 5.9%), Jingmen Tick Virus (JMTV, 0.9%) and Tacheng Tick Virus 1 (TTV1, 0.4%). An association of tick species with individual virus taxa was observed, with the exception of TTV2, which was observed in both Dermacentor and Haemaphysalis species. Individual ticks from pools with pathogen detection were then further screened by targeted amplification and then NS, which provided extensive genome data and revealed probable pathogen Haseki Tick Virus (HTV, 10.2%). Two distinct TTV2 clades were observed in phylogenetic analysis, one of which included closely related Dermacentor reticulatus Uukuviruses. JMTV detection indicated integrated virus sequences. Overall, we observed an expansion of newly documented pathogenic tick-borne viruses into Europe, with TTV1 being identified on the continent for the first time. These viruses should be included in the diagnostic assessment of symptomatic cases associated with tick bites and vector surveillance efforts. NS is shown as a useful tool for monitoring tick-associated pathogens in pooled or individual samples.

Discovery of Rickettsia spp. in mosquitoes collected in Georgia by metagenomics analysis and molecular characterization.

Arthropods have a broad and expanding worldwide presence and can transmit a variety of viral, bacterial, and parasite pathogens. A number of Rickettsia and Orientia species associated with ticks, fleas, lice, and mites have been detected in, or isolated from, patients with febrile illness and/or animal reservoirs throughout the world. Mosquitoes are not currently considered vectors for Rickettsia spp. pathogens to humans or to animals. In this study, we conducted a random metagenome next-generation sequencing (NGS) of 475 pools of Aedes, Culex, and Culiseta species of mosquitoes collected in Georgia from 2018 to 2019, identifying rickettsial gene sequences in 33 pools of mosquitoes. We further confirmed the findings of the Rickettsia by genus-specific quantitative PCR (qPCR) and multi-locus sequence typing (MLST). The NGS and MLST results indicate that Rickettsia spp. are closely related to Rickettsia bellii, which is not known to be pathogenic in humans. The results, together with other reports of Rickettsia spp. in mosquitoes and the susceptibility and transmissibility experiments, suggest that mosquitoes may play a role in the transmission cycle of Rickettsia spp.

Development of conventional and real-time reverse transcription polymerase chain reaction assays to detect Tembusu virus in Culex tarsalis mosquitoes.

Tembusu virus (TMUV) is an important emerging arthropod-borne virus that may cause encephalitis in humans and has been isolated in regions of southeast Asia, including Malaysia, Thailand, and China. Currently, detection and identification of TMUV are limited to research laboratories, because quantitative rapid diagnostic assays for the virus do not exist. We describe the development of sensitive and specific conventional and real-time quantitative reverse transcription polymerase chain reaction assays for detecting TMUV RNA in infected cell culture supernatant and Culex tarsalis mosquitoes. We used this assay to document the replication of TMUV in Cx. tarsalis, where titers increased 1,000-fold 5 days after inoculation. These assays resulted in the detection of virus-specific RNA in the presence of copurified mosquito nucleic acids. The use of these rapid diagnostic assays may have future applications for field pathogen surveillance and may assist in early detection, diagnosis, and control of the associated arthropod-borne pathogens.

Determinants of Anopheles seasonal distribution patterns across a forest to periurban gradient near Iquitos, Peru.

As part of a field ecology study of arbovirus and malaria activity in the Amazon Basin, Loreto Department, Peru, we collected mosquitoes landing on humans at a forest site and inside and outside of residences and military barracks at periurban, rural, and village sites. We collected 11 Anopheles spp. from these four sites. An. darlingi, the principal malaria vector in the region, accounted for 98.7% of all Anopheles spp. collected at Puerto Almendra. Peaks in landing activity occurred during the December and April collection periods. However, the percent of sporozoite-positive Anopheles spp. was highest 1-2 months later, when landing activity decreased to approximately 10% of the peak activity periods. At all sites, peak landing activity occurred about 2 hours after sunset. These data provide a better understanding of the taxonomy, population density, and seasonal and habitat distribution of potential malaria vectors within the Amazon Basin region.