Viral pathogen detection in U.S. game-farm mallard (Anas platyrhynchos) flags spillover risk to wild birds.

The threat posed by emerging infectious diseases is a major concern for global public health, animal health and food security, and the role of birds in transmission is increasingly under scrutiny. Each year, millions of mass-reared game-farm birds are released into the wild, presenting a unique and a poorly understood risk to wild and susceptible bird populations, and to human health. In particular, the shedding of enteric pathogens through excrement into bodies of water at shared migratory stop-over sites, and breeding and wintering grounds, could facilitate multi-species long-distance pathogen dispersal and infection of high numbers of naive endemic birds annually. The Mallard (Anas platyrhynchos) is the most abundant of all duck species, migratory across much of its range, and an important game species for pen-rearing and release. Major recent population declines along the US Atlantic coast has been attributed to game-farm and wild mallard interbreeding and the introduction maladaptive traits into wild populations. However, pathogen transmission and zoonosis among game-farms Mallard may also impact these populations, as well as wildlife and human health. Here, we screened 16 game-farm Mallard from Wisconsin, United States, for enteric viral pathogens using metatranscriptomic data. Four families of viral pathogens were identified - Picobirnaviridae (Genogroup I), Caliciviridae (Duck Nacovirus), Picornaviridae (Duck Aalivirus) and Sedoreoviridae (Duck Rotavirus G). To our knowledge, this is the first report of Aalivirus in the Americas, and the first report of Calicivirus outside domestic chicken and turkey flocks in the United States. Our findings highlight the risk of viral pathogen spillover from peri-domestically reared game birds to naive wild bird populations.

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.

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.

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.

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.

High Levels of Diversity in Anopheles Subgenus Kerteszia Revealed by Species Delimitation Analyses.

The Anopheles subgenus Kerteszia is a poorly understood group of mosquitoes that includes several species of medical importance. Although there are currently twelve recognized species in the subgenus, previous studies have shown that this is likely to be an underestimate of species diversity. Here, we undertake a baseline study of species delimitation using the barcode region of the mtDNA COI gene to explore species diversity among a geographically and taxonomically diverse range of Kerteszia specimens. Beginning with 10 of 12 morphologically identified Kerteszia species spanning eight countries, species delimitation analyses indicated a high degree of cryptic diversity. Overall, our analyses found support for at least 28 species clusters within the subgenus Kerteszia. The most diverse taxon was Anopheles neivai, a known malaria vector, with eight species clusters. Five other species taxa showed strong signatures of species complex structure, among them Anopheles bellator, which is also considered a malaria vector. There was some evidence for species structure within An. homunculus, although the results were equivocal across delimitation analyses. The current study, therefore, suggests that species diversity within the subgenus Kerteszia has been grossly underestimated. Further work will be required to build on this molecular characterization of species diversity and will rely on genomic level approaches and additional morphological data to test these species hypotheses.

Global Distribution of Aedes aegypti and Aedes albopictus in a Climate Change Scenario of Regional Rivalry.

Arboviral mosquito vectors are key targets for the surveillance and control of vector-borne diseases worldwide. In recent years, changes to the global distributions of these species have been a major research focus, aimed at predicting outbreaks of arboviral diseases. In this study, we analyzed a global scenario of climate change under regional rivalry to predict changes to these species' distributions over the next century. Using occurrence data from VectorMap and environmental variables (temperature and precipitation) from WorldClim v. 2.1, we first built fundamental niche models for both species with the boosted regression tree modelling approach. A scenario of climate change on their fundamental niche was then analyzed. The shared socioeconomic pathway scenario 3 (regional rivalry) and the global climate model Geophysical Fluid Dynamics Laboratory Earth System Model v. 4.1 (GFDL-ESM4.1; gfdl.noaa.gov) were utilized for all analyses, in the following time periods: 2021-2040, 2041-2060, 2061-2080, and 2081-2100. Outcomes from these analyses showed that future climate change will affect Ae. aegypti and Ae. albopictus distributions in different ways across the globe. The Northern Hemisphere will have extended Ae. aegypti and Ae. albopictus distributions in future climate change scenarios, whereas the Southern Hemisphere will have the opposite outcomes. Europe will become more suitable for both species and their related vector-borne diseases. Loss of suitability in the Brazilian Amazon region further indicated that this tropical rainforest biome will have lower levels of precipitation to support these species in the future. Our models provide possible future scenarios to help identify locations for resource allocation and surveillance efforts before a significant threat to human health emerges.

A Novel Coronavirus and a Broad Range of Viruses in Kenyan Cave Bats.

BACKGROUND AND METHODS: To investigate virus diversity in hot zones of probable pathogen spillover, 54 oral-fecal swabs were processed from five bat species collected from three cave systems in Kenya, using metagenome sequencing. RESULTS: Viruses belonging to the Astroviridae, Circoviridae, Coronaviridae, Dicistroviridae, Herpesviridae and Retroviridae were detected, with unclassified viruses. Retroviral sequences were prevalent; 74.1% of all samples were positive, with distinct correlations between virus, site and host bat species. Detected retroviruses comprised Myotis myotis, Myotis ricketti, Myotis daubentonii and Galidia endogenous retroviruses, murine leukemia virus-related virus and Rhinolophus ferrumequinum retrovirus (RFRV). A near-complete genome of a local RFRV strain with identical genome organization and 2.8% nucleotide divergence from the prototype isolate was characterized. Bat coronavirus sequences were detected with a prevalence of 24.1%, where analyses on the ORF1ab region revealed a novel alphacoronavirus lineage. Astrovirus sequences were detected in 25.9%of all samples, with considerable diversity. In 9.2% of the samples, other viruses including Actinidia yellowing virus 2, bat betaherpesvirus, Bole tick virus 4, Cyclovirus and Rhopalosiphum padi virus were identified. CONCLUSIONS: Further monitoring of bats across Kenya is essential to facilitate early recognition of possibly emergent zoonotic viruses.

Ticks (Acari: Ixodidae) and Associated Pathoge Collected From Domestic Animals and Vegetation in Stann Creek District, Southeastern Belize, Central America.

Data on the prevalence and distribution of ticks and tick-borne diseases in Belize are lacking. Ticks (n = 564) collected from dogs, horses, and vegetation in two villages in Stann Creek District in southeastern Belize in 2018, were molecularly identified and screened for tick-borne nonviral human pathogens. The identity of 417 ticks was molecularly confirmed by DNA barcoding as Rhipicephalus sanguineus (Latreille) (66.43%), Amblyomma ovale Koch (15.59%), Dermacentor nitens Neumann (11.51%), Amblyomma sp. ADB0528 (3.6%), and the remainder being small records (2.87%) of Amblyomma coelebs Neumann, Amblyomma imitator Kohls, Amblyomma tapirellum Dunn, Amblyomma auricularium Conil, and Amblyomma maculatum Koch. Individual tick extracts were screened for the presence of Rickettsia spp., Babesia spp., Babesia microti, Borrelia spp., Ehrlichia spp., and Anaplasma spp. using available conventional polymerase chain reaction (PCR) assays. Rickettsia parkeri strain Atlantic Rainforest was identified in five specimens of A. ovale, and one other unidentified tick, all collected from dogs. Another unidentified tick-also collected from a dog-tested positive for an undefined but previously detected Ehrlichia sp. With the exception of D. nitens, all eight other tick species identified in this study were collected on dogs, suggesting that dogs could be usefully employed as sentinel animals for tick surveillance in Belize.

Phylogenetic analysis of the Neotropical Albitarsis Complex based on mitogenome data.

BACKGROUND: Some of the most important malaria vectors in South America belong to the Albitarsis Complex (Culicidae; Anophelinae; Anopheles). Understanding the origin, nature, and geographical distribution of species diversity in this important complex has important implications for vector incrimination, control, and management, and for modelling future responses to climate change, deforestation, and human population expansion. This study attempts to further explore species diversity and evolutionary history in the Albitarsis Complex by undertaking a characterization and phylogenetic analysis of the mitogenome of all 10 putative taxa in the Albitarsis Complex. METHODS: Mitogenome assembly and annotation allowed for feature comparison among Albitarsis Complex and Anopheles species. Selection analysis was conducted across all 13 protein-coding genes. Maximum likelihood and Bayesian inference methods were used to construct gene and species trees, respectively. Bayesian methods were also used to jointly estimate species delimitation and species trees. RESULTS: Gene composition and order were conserved across species within the complex. Unique signatures of positive selection were detected in two species-Anopheles janconnae and An. albitarsis G-which may have played a role in the recent and rapid diversification of the complex. The COI gene phylogeny does not fully recover the mitogenome phylogeny, and a multispecies coalescent-based phylogeny shows that considerable uncertainty exists through much of the mitogenome species tree. The origin of divergence in the complex dates to the Pliocene/Pleistocene boundary, and divergence within the distinct northern South American clade is estimated at approximately 1 million years ago. Neither the phylogenetic trees nor the delimitation approach rejected the 10-species hypothesis, although the analyses could not exclude the possibility that four putative species with scant a priori support (An. albitarsis G, An. albitarsis H, An. albitarsis I, and An. albitarsis J), represent population-level, rather than species-level, splits. CONCLUSION: The lack of resolution in much of the species tree and the limitations of the delimitation analysis warrant future studies on the complex using genome-wide data and the inclusion of additional specimens, particularly from two putative species, An. albitarsis I and An. albitarsis J.

Molecular species delimitation reveals high diversity in the mosquito Anopheles tessellatus Theobald, 1901 (Diptera, Culicidae) across its range.

Anopheles tessellatus is a potentially important vector found across South, East and Southeast Asia. While it was formerly considered a formidable vector of human Plasmodium and filarial parasites in the Maldives, and of lesser importance as a vector of human Plasmodium in Sri Lanka and parts of Indonesia, it is currently of little or unknown health importance in many other parts of its range. This study describes the genetic diversity and evolutionary relationships among An. tessellatus populations in nine Asian countries at the COI gene using maximum-likelihood and Bayesian phylogenetic inference tree and cluster-based species delimitation approaches. These analyses reveal exceptional levels of genetic diversity in An. tessellatus populations across its known range, and identify up to six putative species in the newly determined Tessellatus Complex. The existence of such cryptic diversity has potentially important consequences for vector management and disease control. Differences in the ecologies and life histories among these species may have considerable impact on vectorial capacity and may go some way towards explaining why An. tessellatus s.l. has such varying degrees of public health importance across its range.

Exploring malaria vector diversity on the Amazon Frontier.

BACKGROUND: Deforestation in the Amazon and the social vulnerability of its settler communities has been associated with increased malaria incidence. The feeding biology of the most important malaria vectors in the region, notably Nyssorhynchus darlingi, compounds efforts to control vectors and reduce transmission of what has become known as "Frontier Malaria". Exploring Anophelinae mosquito diversity is fundamental to understanding the species responsible for transmission and developing appropriate management and intervention strategies for malaria control in the Amazon River basin. METHODS: This study describes Anophelinae mosquito diversity from settler communities affected by Frontier Malaria in the states of Acre, Amazonas and Rondônia by analysing COI gene data using cluster and tree-based species delimitation approaches. RESULTS: In total, 270 specimens from collection sites were sequenced and these were combined with 151 reference (GenBank) sequences in the analysis to assist in species identification. Conservative estimates found that the number of species collected at these sites was between 23 (mPTP partition) and 27 (strict ABGD partition) species, up to 13 of which appeared to be new. Nyssorhynchus triannulatus and Nyssorhynchus braziliensis displayed exceptional levels of intraspecific genetic diversity but there was little to no support for putative species complex status. CONCLUSIONS: This study demonstrates that Anophelinae mosquito diversity continues to be underestimated in poorly sampled areas where frontier malaria is a major public health concern. The findings will help shape future studies of vector incrimination and transmission dynamics in these areas and support efforts to develop more effective vector control and transmission reduction strategies in settler communities in the Amazon River basin.

Spatio-temporal genetic structure of Anopheles gambiae in the Northwestern Lake Victoria Basin, Uganda: implications for genetic control trials in malaria endemic regions.

BACKGROUND: Understanding population genetic structure in the malaria vector Anopheles gambiae (s.s.) is crucial to inform genetic control and manage insecticide resistance. Unfortunately, species characteristics such as high nucleotide diversity, large effective population size, recent range expansion, and high dispersal ability complicate the inference of genetic structure across its range in sub-Saharan Africa. The ocean, along with the Great Rift Valley, is one of the few recognized barriers to gene flow in this species, but the effect of inland lakes, which could be useful sites for initial testing of genetic control strategies, is relatively understudied. Here we examine Lake Victoria as a barrier between the Ugandan mainland and the Ssese Islands, which lie up to 60 km offshore. We use mitochondrial DNA (mtDNA) from populations sampled in 2002, 2012 and 2015, and perform Bayesian cluster analysis on mtDNA combined with microsatellite data previously generated from the same 2002 mosquito DNA samples. RESULTS: Hierarchical analysis of molecular variance and Bayesian clustering support significant differentiation between the mainland and lacustrine islands. In an mtDNA haplotype network constructed from this and previous data, haplotypes are shared even between localities separated by the Rift Valley, a result that more likely reflects retention of shared ancestral polymorphism than contemporary gene flow. CONCLUSIONS: The relative genetic isolation of An. gambiae on the Ssese Islands, their small size, level terrain and ease of access from the mainland, the relative simplicity of the vectorial system, and the prevalence of malaria, are all attributes that recommend these islands as possible sites for the testing of genetic control strategies.

Phylogenetic analysis and DNA-based species confirmation in Anopheles (Nyssorhynchus).

Specimens of neotropical Anopheles (Nyssorhynchus) were collected and identified morphologically. We amplified three genes for phylogenetic analysis-the single copy nuclear white and CAD genes, and the COI barcode region. Since we had multiple specimens for most species we were able to test how well the single or combined genes were able to corroborate morphologically defined species by placing the species into exclusive groups. We found that single genes, including the COI barcode region, were poor at confirming species, but that the three genes combined were able to do so much better. This has implications for species identification, species delimitation, and species discovery, and we caution that single genes are not enough. Higher level groupings were partially resolved with some well-supported groupings, whereas others were found to be either polyphyletic or paraphyletic. There were examples of known groups, such as the Myzorhynchella Section, which were poorly supported with single genes but were well supported with combined genes. From this we can infer that more sequence data will be needed in order to show more higher-level groupings with good support. We got unambiguously good support (0.94-1.0 Bayesian posterior probability) from all DNA-based analyses for a grouping of An. dunhami with An. nuneztovari and An. goeldii, and because of this and because of morphological similarities we propose that An. dunhami be included in the Nuneztovari Complex. We obtained phylogenetic corroboration for new species which had been recognised by morphological differences; these will need to be formally described and named.

Mosquito (Diptera: Culicidae) assemblages associated with Nidularium and Vriesea bromeliads in Serra do Mar, Atlantic Forest, Brazil.

BACKGROUND: The most substantial and best preserved area of Atlantic Forest is within the biogeographical sub-region of Serra do Mar. The topographic complexity of the region creates a diverse array of microclimates, which can affect species distribution and diversity inside the forest. Given that Atlantic Forest includes highly heterogeneous environments, a diverse and medically important Culicidae assemblage, and possible species co-occurrence, we evaluated mosquito assemblages from bromeliad phytotelmata in Serra do Mar (southeastern Brazil). METHODS: Larvae and pupae were collected monthly from Nidularium and Vriesea bromeliads between July 2008 and June 2009. Collection sites were divided into landscape categories (lowland, hillslope and hilltop) based on elevation and slope. Correlations between bromeliad mosquito assemblage and environmental variables were assessed using multivariate redundancy analysis. Differences in species diversity between bromeliads within each category of elevation were explored using the Renyi diversity index. Univariate binary logistic regression analyses were used to assess species co-occurrence. RESULTS: A total of 2,024 mosquitoes belonging to 22 species were collected. Landscape categories (pseudo-F value = 1.89, p = 0.04), bromeliad water volume (pseudo-F = 2.99, p = 0.03) and bromeliad fullness (Pseudo-F = 4.47, p < 0.01) influenced mosquito assemblage structure. Renyi diversity index show that lowland possesses the highest diversity indices. The presence of An. homunculus was associated with Cx. ocellatus and the presence of An. cruzii was associated with Cx. neglectus, Cx. inimitabilis fuscatus and Cx. worontzowi. Anopheles cruzii and An. homunculus were taken from the same bromeliad, however, the co-occurrence between those two species was not statistically significant. CONCLUSIONS: One of the main findings of our study was that differences in species among mosquito assemblages were influenced by landscape characteristics. The bromeliad factor that influenced mosquito abundance and assemblage structure was fullness. The findings of the current study raise important questions about the role of An. homunculus in the transmission of Plasmodium in Serra do Mar, southeastern Atlantic Forest.

Phylogenetic relationships among species of Anopheles (Nyssorhynchus) (Diptera, Culicidae) based on nuclear and mitochondrial gene sequences.

Phylogenetic relationships among 21 species of mosquitoes in subgenus Nyssorhynchus were inferred from the nuclear white and mitochondrial NADH dehydrogenase subunit 6 (ND6) genes. Bayesian phylogenetic methods found that none of the three Sections within Nyssorhynchus (Albimanus, Argyritarsis, Myzorhynchella) were supported in all analyses, although Myzorhynchella was found to be monophyletic at the combined genes. Within the Albimanus Section the monophyly of the Strodei Subgroup was strongly supported and within the Myzorhynchella Section Anopheles antunesi and An. lutzii formed a strongly supported monophyletic group. The epidemiologically significant Albitarsis Complex showed evidence of paraphyly (relative to An. lanei-Myzorhynchella) and discordance across gene trees, and the previously synonomized species of An. dunhami and An. goeldii were recovered as sister species. Finally, there was evidence of complexes in several species, including An. antunesi, An. deaneorum, and An. strodei.