Arboviral screening of invasive Aedes species in northeastern Turkey: West Nile virus circulation and detection of insect-only viruses

Background The recent reports of Aedes aegypti and Ae. albopictus populations in Turkey, in parallel with the territorial expansion identified in several surrounding countries, have raised concerns about the establishment and re-establishment of these invasive Aedes mosquitoes in Turkey. This cross-sectional study was performed to detect Aedes aegypti and Ae. albopictus in regions of recent incursions, and screen for viral pathogens known to be transmitted elsewhere by these species. Methodology Mosquitoes were collected at several locations in Artvin, Rize and Trabzon provinces of the Black Sea region during 2016–2017, identified morphologically, pooled and analyzed via generic or specific nucleic acid amplification assays. Viruses in positive pools were identified by product sequencing, cell culture inoculation and next generation sequencing (NGS) in selected specimens. Principal findings The study group comprised 791 specimens. Aedes albopictus was the most abundant species in all locations (89.6%), followed by Ae. aegypti (7.8%) and Culex pipiens (2.5%). Mosquitoes were screened for viruses in 65 pools where fifteen (23.1%) were reactive. The infecting strains was identified as West Nile virus (WNV) in 5 pools (7.7%) with Ae. albopictus or Cx. pipiens mosquitoes. The obtained WNV sequences phylogenetically grouped with local and global lineage 1 clade 1a viruses. In 4 (6.2%) and 6 (9.2%) pools, respectively, cell fusing agent virus (CFAV) and Aedes flavivirus (AEFV) sequences were characterized. NGS provided a near-complete AEFV genome in a pool of Ae. albopictus. The strain is provisionally called “AEFV-Turkey”, and functional analysis of the genome revealed several conserved motifs and regions associated with virus replication. Merida-like virus Turkey (MERDLVT), a recently-described novel rhabdovirus, was also co-detected in a Cx. pipiens pool also positive for WNV. Conclusions/Significance Invasive Aedes mosquitoes are established in certain locations of northeastern Turkey. Herein we conclusively show the role of these species in WNV circulation in the region. Biosurveillance is imperative to monitor the spread of these species further into Asia Minor and to detect possible introduction of pathogens.


Introduction
Infections due to mosquito-borne viruses have become a global health problem during the past two decades, due to their wide geographic spread and high human disease burden. This is directly associated with their widespread distribution and ecological changes related to vector mosquitoes as well as increases in international trade and travel [1,2]. Any list of the mosquito-borne viruses with significant impact that have emerged or re-emerged would include dengue virus (DENV), yellow fever virus (YFV), chikungunya virus (CHIKV) and recently, Zika virus (ZIKV) [2,3]. These viruses are transmitted to susceptible vertebrates with varying degrees of vector competence via the globally invasive Aedes (Stegomyia) aegypti (L.) and Ae. (Stg.) albopictus (Skuse) [4]. In Europe, Ae. albopictus is the most prolific invasive mosquito species, having greatly expanded its range across many countries since its original introduction to Albania in 1979 [4,5]. Despite a more constrained geographical distribution, Ae. aegypti has significantly recolonized parts of southern and southeastern Europe with populations found in Portugal and the Black Sea coast of Russia, Republic of Georgia and, most recently, in northeastern Turkey [4,5]. Risk mapping efforts have identified 215 countries or territories to be potentially suitable for the survival and establishment of invasive Aedes species, with disease outbreaks from more than half of the target regions [6]. Therefore, concerns of virus transmission that could initiate and sustain epidemics in countries infested by these mosquitoes have been raised, requiring vector and pathogen surveillance [3,4].
In Turkey, Ae. albopictus was initially recorded in eastern Thrace (Edirne province, bordering Greece) in 2011 [7], along with detections from several Balkan countries including Bulgaria and Romania [8]. Following DDT-based eradication efforts in the 1950s, Aedes aegypti was rarely identified around the Mediterranean basin, with only sporadic later reports from Turkey, Italy and Israel [8]. However, both these invasive species were recently identified at several locations from the coastal Black Sea region of Turkey, with anthropophagic adults and immature aquatic stages in used tyres stored outdoors [9]. Established populations of both species are present in neighboring Republic of Georgia suggesting local encroachment as the source of the Turkish populations. Aedes aegypti populations from various regions around the Black Sea demonstrate a high genetic differentiation and are hypothesized to represent expansions from remnant populations within the area [10]. Thus, a resurgence or reintroduction of viruses vectored by these species must be considered in NE Turkey, along with appropriate strategies for routine vector biosurveillance and control. Among mosquito-borne viral pathogens, West Nile virus (WNV) seems to be ubiquitous in Anatolia and Thrace regions, with several reported cases of human and equine infections [11][12][13][14]. Moreover, vector screening efforts have identified WNV in both Culex and Aedes mosquitoes in various regions [14][15][16]. So far, no human infections due to indigenous transmission of DENV, YFV, CHIKV or ZIKV have been documented in Turkey. However, serological evidence of sporadic exposure to DENV or an antigenically-related flavivirus [17], as well as imported cases of DENV, CHIKV and ZIKV have been reported [18][19][20]. In this study, we aimed to detect newly established populations of invasive mosquito species including Ae. aegypti and Ae. albopictus in regions of the Black Sea coast and screen them for viral pathogens known to be transmitted by these species.

Ethics statement
The study involved testing of field-collected mosquitoes for which no local or institutional ethics committee approval is required. Peridomestic collections were undertaken with informed consent and cooperation of the property owners, householders or local authorities.

Study area, specimen collection and identification
Mosquito sampling was undertaken at 32 locations in Artvin, Rize and Trabzon provinces of the Black Sea region from June through October, in 2016 and 2017 (Fig 1, S1 Table). Members of the trained entomology team performed the samplings indoors and outdoors at urban sites using human landing catch method, as described previously [21]. Hepa Filter Mouth Aspirators and Prokopack Aspirator (John W. Hock Company, Gainesville, FL, USA) were employed specimen collection [22]. All collected mosquitoes were transferred on ice, killed by freezing and identified to species level using morphological keys [23,24]. Subsequently, the specimens were pooled according to the collection site, species and sex and stored at -80˚C.

Specimen processing and barcoding
Pooled mosquitoes were disrupted by vortexing following the addition of stainless steel beads (Qiagen, Hilden, Germany) and 500 microliters of Eagle's minimal essential medium, supplemented with 5% fetal bovine serum, 1% L-glutamine and 1% penicillin-streptomycin. Subsequently, the pools were cleared by centrifugation at 4000 rpm for 4 minutes and the supernatants were collected. They were subjected to nucleic acid purificaton via High Pure Viral Nucleic Acid Kit (Roche Diagnostics, Mannheim, Germany), followed by reverse transcription for complementary DNA synthesis, using the RevertAid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Hennigsdorf, Germany) by random hexamer priming, as directed by the manufacturer.
Mosquito pools with detectable virus sequences were further subjected to DNA barcoding for the confirmation of the species identification. For this purpose, a portion of the cytochrome c oxidase I (COI) gene was amplified and sequenced using LCO1490 and HCO2198 primers [25]. (Red: countries/territories with Ae. aegypti and Ae. albopictus; Orange: countries/ territories with Ae. albopictus). Blue dots represent sampling locations. The baseline map has been prepared using Natural Earth raster + vector map data in the public domain (URL: www.naturalearthdata.com. Accessed: April 2019), which is freely available for personal, educational, and commercial use. Current information on Aedes species were obtained from the European Centre for Disease Prevention and Control websites (https://ecdc.europa.eu/en/publicationsdata/aedes-aegypti-current-known-distribution-june-2018; https://ecdc.europa.eu/en/publications-data/aedes-albopictus-current-known-distribution-june-2018; Accessed: December 2018).

Flavivirus and Zika Virus screening
A nested PCR assay with degenerated primers targeting the NS5 region of the flavivirus genome was used for screening in the mosquito pools [26]. The primers were reported to provide sensitive amplification of all insect-specific and major mosquito-borne flaviviruses including WNV, DENV, YFV, ZIKV, Saint Louis encephalitis virus and Usutu virus, strains with a detection limit of 40 TCID 50 per reaction [26]. Assay optimization was undertaken using purified and serially-diluted nucleic acids from WNV NY99-4132 (standard) and T2 (local) isolates, grown on African green monkey (Vero) cells (ATCC-CCL81).
All pools were further tested individually for Zika virus, using the previously described realtime PCR assay targeting the envelope glycoprotein coding region [27]. The primers ZIKV1086-1162c and probe ZIKV1107 were used for the single step amplification via QuantiNova Patho-gen+IC Kit (QIAgen, Hilden, Germany) in a Rotor-Gene 6000 instrument (Corbett Research, Australia). The assay was reported to have a detection threshold of 25 copies per reaction [23]. Zika virus strain UVE/ZIKV/1947/UG/MR766 RNA, obtained from the European Virus Archive, was used for assay optimization and as the positive template during screening.

Alphavirus screening
Mosquito pools were screened for alphaviruses, using the nested PCR employing degenerate primers targeting the nsP4 region, as described previously [28]. The assay could amplify several medically-important alphaviruses including CHIKV, Sindbis, O'nyong nyong, western equine encephalitis, eastern equine encephalitis, Venezuelan equine encephalitis, Semliki forest, Ross river and Barmah forest viruses. The detection threshold of the assay was reported as 1-10 copies (cloned fragment) or 25 pfu (cell-grown virus) [24]. CHIKV strain LR2006-OPY1 cDNA, obtained from the European Virus Archive was used for optimization and as a positive template during screening.

Rhabdovirus screening
We further screened the mosquito pools for the recently-characterized rhabdovirus, tentatively named Merida-like virus Turkey (MERDLVT) [29]. Two PCR assays, designed to amplify 481 and 160 base pairs from the L-and N-regions of the MERDLVT genome was used for screening the mosquito pools [29]. The assays were optimized using previously-collected positive mosquito pools.
Products of the screening assays were visualized under ultraviolet light following electrophoresis in 1.5-2.0% agarose gels, depending on the amplicon size. Precautions to prevent carry-over contamination were strictly followed and pre and post-PCR steps were performed in spatially-separated areas, with several non-template controls during each run.

Virus isolation
Aliquots of the mosquito pools positive in the screening assays were inoculated onto semi-confluent monolayers of Vero (ATCC-CCL81) and Ae. albopictus (C6/36, ATCC-CRL1660) cells, following filtration through 0.22 micrometer sterile membrane filters (Merck Millipore, Darmstadt, Germany). The cells were incubated at 37˚C and 28˚C respectively, and monitored daily for cytopathic effects. Blind passages to fresh monolayers and testing for viruses via the screening assays were carried out weekly.

Sanger and next generation sequencing (NGS)
Products of the mosquito COI barcoding and virus screening assays were cleaned up using PureLink PCR Purification Kit (Thermo Fisher Scientific, Hennigsdorf, Germany) and sequenced using forward-reverse primers of the particular assay and the BigDye Terminator v3.1 Cycle Sequencing Kit (Thermo Fisher Scientific) in an ABI PRISM 3500xL Dx genetic analyzer (Thermo Fisher Scientific).
Mosquito pools positive in the screening assays and with available aliquots were subjected to direct NGS. Following purification, the nucleic acids were reverse transcribed with random hexamer primers to double-stranded cDNA using SuperScript IV Reverse Transcriptase (Thermo Fisher Scientific, Hennigsdorf, Germany) and NEBNext mRNA Second Strand Synthesis Module (New England BioLabs, Frankfurt am Main, Germany). The Agilent 2100 Bioanalyzer (Agilent Technologies, Waldbronn, Germany) and Agencourt AMPure XP Reagent (Beckman Coulter Biosciences, Krefeld, Germany) were used for cleanup and estimation of yield and size distribution. Fragmentation, adaptor ligation and amplification were performed according to the manufacturer protocols using the NexteraXT DNA Library Preparation Kit (Illumina Inc., San Diego, CA, USA). Sequencing runs were performed on the Illumina MiSeq (Illumina Inc.) instrument in the paired end mode.

Sequence data analysis
Raw sequences obtained from virus screening were handled using Geneious software v11.1.5 (Biomatters Ltd, Auckland, New Zealand). Trimmomatic [30] was used for adaptor removal, trimming for quality and length with a phred score of 33 and a minimum length of 30 base pairs (bp). Acquired reads were aligned to the RefSeq viral nucleotide and protein genome database using MALT (MEGAN alignment tool, v0.3.8) and MEGAN (Metagenome Analyzer, v. 6.12.3) [31,32]. Aligned reads were extracted and assembled into contigs using Velvet (v.1.2.10) [33]. The contigs were mapped to closely related virus strains, checked for heterogeneity via visual inspection and pairwise identity values using Geneious software v11.1.5 (Biomatters Ltd, Auckland, New Zealand).
BLASTn, BLASTn optimized for highly similar sequences (MEGABLAST) and BLASTp algorithms were employed for nucleotide and deduced amino acid similarity searches in the public databases, implemented in the NCBI website (www.ncbi.nlm.nih.gov/blast/) [34]. Nucleotide and putative amino acid alignments and pairwise sequence comparisons were generated via the CLUSTAL W program, implemented within Geneious software [35]. Nucleotide identity plots were generated by SimPlot version 3.5.1 [36]. Conserved protein domain and motif searches were performed using the web search tool (http://www.ncbi.nlm.nih.gov/ structure/bwrpsb/bwrpsb.cgi) and MOTIF Search (http://www.genome.jp/tools/motif/) in the PFAM database [37,38]. Evolutionary history was inferred via the maximum-likelihood method based on the model estimated as the optimal substitution model individually for each alignment according to the Bayesian information criterion and conducted using MEGA6 [39].
Mosquitoes were screened for viruses in 65 pools. Fifteen of these (23.1%) were reactive in at least one screening assay. Generic flavivirus assay was positive in all reactive pools, but the

Flavivirus findings
The detected flavivirus was characterized via amplicon sequencing in PCR positive pools. WNV sequences were identified in 5 pools (5/65, 7.7%) comprising Ae. albopictus (4/5) and Cx. pipiens s.l. (1/5) mosquitoes ( Table 2). Sequences of 777-1,000 base pairs (bp) were characterized. Pairwise comparisons showed diversity rates of 0.2-0.7% and 0.4-1.6% in the nucleotide and deduced amino acids, respectively. In the maximum likelihood tree, all sequences grouped within WNV lineage 1 clade 1a sequences, and formed a distinct cluster with viruses from the American Continent, Israel, Tunisia, Hungary and Aegean coast of Anatolia (Turkey) (Fig 2). Interestingly, WNV sequences of mosquito, equine and avian origins from Turkey grouped within different clusters in the lineage 1 clade 1a viruses. In 4 of the tested pools (6.2%) comprising 2 Ae. albopictus and 2 Ae. aegypti specimens, sequences identified as cell fusing agent virus (CFAV) were detected. The sequences were partially overlapping 283-695 bp segments, covering the flavivirus NS5 region amplicon. They displayed 86-97.8% identity to the CFAV strain Galveston (GenBank accession NC001564) in pairwise comparisons and formed a well-supported phylogenetic group with CFAV strains of various origins in the maximum likelihood tree (Fig 3).
Another flavivirus identified following amplicon sequencing was Aedes flavivirus (AEFV). AEFV was detected in a total of 6 pools (9.2%), 5 comprising Ae. albopictus (5/6) and 1 with Ae. aegypti ( Table 2). The obtained sequences comprised 283-695 bp of partially overlapping segments of the target amplicon, with 0.9-8.8% intramural diversity in pairwise comparisons. Maximum likelihood analysis using the sequences from the pools KRD1, KRD32 and P11 revealed a separate clustering of these sequences with the closely-related AEFV strains, among other distinct insect flaviviruses, including CFAV (Fig 3).
Despite lack of virus isolation in cell culture inoculation of flavivirus PCR positive pools, we could obtain near-complete polyprotein coding region of the AEFV via direct NGS in available aliquots in the mosquito pool, KRD32 ( Table 2). The assembled sequence comprised 9,931 bp, with the deduced polyprotein of 3310 amino acids, that covered over 99% of the complete AEFV polyprotein. The sequence demonstrated pairwise diversity rates of 1-8.5% and 0.6-2.2% on the nucleotide and deduced amino acid levels, respectively; when compared to AEFV strains Narita-21 (GenBank accession AB488408), Bangkok (KJ741266), SPFLD-MO-2011-MP6 (KC181923) and the recently described La Tina virus isolate 49 (KY320649). Region-specific identities were further determined on the viral polyprotein, which revealed similarity rates of >90% on individual mature proteins ( Table 3). Maximum likelihood analysis of the near-complete genome revealed a tree topology, comparable to the partial NS5 tree, with well-supported grouping of AEFV isolates (Fig 4). A genome-wide region specific nucleotide-based comparison is further provided in the SimPlot graph (Fig 5). The virus is provisionally named as AEFV-Turkey, as it represents the first near-complete AEFV genome reported from Asia minor. The insect-specific flavivirus ribosomal frameshifting site, that results in a longer overlapping ORF in the NS2A-NS2B regions [40], was observed as GGATTTT heptanucleotide motif, encompassing the nucleotides 3277-3283 in the AEFV-Turkey genome. Motifs of flavivirus envelope glycoprotein with the central/dimerisation domains (residues 306-478), flavivirus non-structural protein 1 (NS1) (879-1,052), NS3 serine protease (1,471-1,607), DEAD-like helicase domain (1,627-1,768), methyl transferase (2,475-2,663) and flavivirus RNA-dependent RNA polymerase (2,963-3,309) were identified in the deduced viral polyprotein. The envelope fusion peptide motif, involved in viral endosomal fusion and cellular entry [41], was located in 361-374. residues of the deduced viral polyprotein and characterized as NRGWGTGCFEWGLG.

Rhabdovirus findings
MERDLVT specific PCRs were positive in a single pool of molecularly confirmed Cx. pipiens s.s. ( Table 2)

Discussion
Targeted pathogen monitoring in arthropods is an important part of surveillance programs for the circulation of vector-borne agents and likely to predict probable disease emergence in susceptible human or animal populations [1]. The recent detection of invasive Aedes mosquitoes in Turkey, in parallel with the territorial expansion identified in several countries, have raised concerns about the reestablishment of these species as well as emergence of associated viruses [6,9]. This study was carried out to investigate the presence of the invasive Aedes species around the sites of previous detection and fill the current information gap on the circulation of viruses potentially spread by these mosquitoes. Specimen collection encompassed mosquito active seasons for two consecutive years, utilizing human landing catches, considered as a the most highly sensitive and effective approach for sampling anthropophilic Aedes species [21]. We detected Ae. albopictus and Cx. pipiens s.l. in all provinces targeted for surveillance and throughout the screening period (Table 1). However, Ae. aegypti was lacking in the Trabzon province, despite being identified regularly in the remaining provinces. Aedes albopictus was detected in high frequencies in all sampling sites, and comprised 89.6% of the study cohort overall. These findings indicate that these invasive Aedes are firmly established in the coastal Black Sea region of northeastern Turkey, at least in certain locations. It has been previously suggested that further spread of these species to the ports around the Black Sea via ships and ferries is probable, as well as dispersion via ground transportation into major cities of Turkey [9]. Therefore, continuous and integrated surveillance of invasive mosquitoes are imperative to monitor the spread of these species further into Asia Minor and to implement effective control strategies as become necessary.
We screened the field-collected mosquitoes for a diverse spectrum of viruses, including pathogenic flavi and alphaviruses mainly transmitted by Aedes mosquitoes, via generic or realtime PCR assays. WNV was noted as the prominent pathogenic agent, detected in 7.7% of the mosquito pools ( Table 2). WNV circulates among various avian species and vector mosquitoes in nature. Humans and horses are exposed to the virus via infected mosquitoes and considered as dead-end hosts, due to the lack of prolonged and high-level viremia required to contribute  Table 3. Functional organization of the AEFV-Turkey genome and comparison with related viruses. Similarity rates are given in percent.  [14][15][16]43], in organ specimens from migratory birds [44], as well as in symptomatic humans and equine infections, occasionally presenting as outbreaks [12,13], However, most data on WNV circulation originated from Aegean, Mediterranean, Thrace and Central Anatolian locations with scarce information from the Black Sea region. We have identified WNV sequences in Ae. albopictus and Cx. pipiens pools collected from Artvin and Rize provinces in this study ( Table 2), clearly indicating virus circulation in the region. We have previously evaluated mosquito specimens from Artvin province in 2013, comprising mostly Culex spp., without WNV detection [14]. However, WNV neutralizing antibodies were detected in 9.9% of the duck sera originating from the neighboring Kars province, suggesting prior virus exposure in the region [14]. It is known that different mosquito species possess highly variable potential to acquire and transmit WNV and Culex species are accepted as the primary global transmission vector [45,46]. WNV has been detected several other genera of mosquitoes including Aedes, Anopheles, Coquillettidia, Culiseta, Mansonia, Mimomyia, Invasive Aedes-associated viruses in Turkey

AEFV-Turkey AEFV-Bangkok AEFV-Narita21 AEFV-SPFLD LaTina virus-LT96
Ochlerotatus, Psorophora, and Uranoteania, that can serve as bridge vectors for transmission from birds to humans and equines [45,46]. Aedes albopictus mosquitoes are competent vectors for WNV, however their contribution to virus circulation in the field is considered as limited, due to their feeding preference for humans and variations in WNV transmission rates [47]. This may explain the lack of documented clinical cases from the screened provinces so far. Of note is the relatively hign WNV incidence in A. albopictus pools, which suggests considerable intensity of virus circulation in the sampling locations. Therefore, a detailed surveillance is required to better understand WNV epidemiology in the Black Sea region. Five major WNV lineages have been described according to the genomic phylogenies, where lineage 1 is widely distributed throughout Africa, Asia and America [42]. However, other lineages also circulate in Europe where lineage 2 may cause human infections [48]. The partial WNV sequences characterized in mosquito pools in this study grouped phylogenetically with lineage 1 clade 1a sequences, which include the majority of the global lineage 1 strains as well as previously characterized sequences in Anatolia [11]. However, a significant WNV sequence diversity has also been documented in Turkey [16,42], which is represented in this study as differential clustering of mosquito, avian and equine sequences (Fig 2). This pronounced diversity, as well as occasional detection of lineage 2 strains [49], probably results from independent virus introductions and dispersion via migrating birds throughout Anatolia and Thrace [43].
Besides pathogenic flaviviruses, our screening provided information on insect-specific flaviviruses (ISFs) in mosquitoes of the Black Sea region of Anatolia. Forming a major phylogenetic group of Flaviviridae, the ISFs do not seem to infect vertebrates and replicate exclusively in mosquito derived cell lines [40]. ISFs demonstrate a widespread geographic distribution and have been detected in a wide range of mosquito species, including those that act as virus vectors. Therefore, they share identical ecologic niches with vector-borne pathogens and frequently co-circulate in given geographical areas [40]. We have identified two distinct ISFs, Invasive Aedes-associated viruses in Turkey CFAV and AEFV, in this study. They could be detected in several sampling sites during consecutive years in Ae. albopictus and Ae. aegypti mosquitoes ( Table 2), suggesting dissemineated circulation in the region. CFAV was the first ISF to be discovered, initially isolated from an infected Ae. aegypti cell line [50]. It has subsequently been detected in several Aedes and Culex spp. mosquitoes from locations in Asia, Africa and the American continent [40], and recently reported from Brasil [51]. AEFV is another ubiquituous ISF, previously also detected in mosquitoes from Italy and United States, following its initial characterization in Japan [52][53][54][55][56]. Mostly identified in Ae. albopictus and Aedes flavopictus mosquitoes, AEFV was also present in Cx. pipiens mosquitoes from Italy, suggesting transspecies infections similar to CFAV [55]. Both viruses were also isolated from laboratory colonies established from mosquitoes collected in Thailand and the United States [57]. In our study, the local AEFV strain provisionally named as "AEFV-Turkey", was detected in Ae. aegypti as well as in Ae. albopictus mosquitoes ( Table 2). The near complete genome of AEFV-Turkey revealed significant sequence similarities, identical organization with specific regions and functional markers with AEFVs ( Table 3, Fig 5). We have previously isolated and characterized several ISFs in Anatolia and Thrace, including Culex theileri flavivirus Turkey [58], Ochlerotatus caspius flavivirus Turkey [16], Anopheles flavivirus and evidence for novel ISFs [43]. Therefore, AEFV-Turkey becomes the latest addition to the list of ISFs known to circulate in the Anatolian mosquito fauna. Of particular interest is the probable interaction of ISFs with pathogenic flaviviruses, resulting in alterations in vector infection efficiency or transmission dynamics [40]. However, currently available information is insufficient for well-supported conclusions and detailed screening is likely to provide epidemiological data required for a better understanding of the ISF-pathogen interactions occuring in nature.
We have further detected MERDLVT in a pool of Cx. pipiens s.s. mosquitoes, also positive for WNV ( Table 2). MERDLVT and closely related MERDV are putative members of the family Rhabdoviridae, identified using NGS without isolation in field-collected mosquitoes from Mexico and Turkey, without succesful isolation of a viable strain [29,59]. MERDLVT has previously been detected in Cx. pipiens mosquitoes from several locations in Mediterranean and Aegean Anatolia as well as from Thrace [29,43]. The detection of MERDLVT within a relatively small Culex cohort indirectly suggests prominent virus circulation in the Black Sea region. MERDLVT has so far been observed exclusively in Culex mosquitoes [29,43]. However, mosquitoes with MERDV infection are not limited to Culex genera, and partial viral sequences were detected in Aedes taeniorhynchus and Aedes trivittatus mosquitoes among seven mosquitoes species in Mexico [59]. The MERDLVT sequences are generally well-conserved, as observed in the amplified sections of the viral L and N genes in this study, and phylogenetically-related to the previously identified sequences (Fig 4). We have previously documented the cocirculation of flaviviruses and MERDLVT, which is also observed in this study [43]. Similar to ISFs, the impact of MERDLVT or MERDV on vector survival and pathogen transmission remains currently unexplored and requires further investigation.
In conclusion, we have identified ongoing activity of invasive Aedes mosquitoes in Black Sea region of Anatolia. WNV circulation is documented for the first time in potential mosquito vectors in the region. No evidence of recently-emergent Zika, Chikungunya or other pathogenic flavi/alphavirus was observed. Insect-associated flavi and rhabdoviruses were detected, with near-complete genome of AEFV, reported initially from Anatolia.