Population Genetics of Culex tritaeniorhynchus (Diptera: Culicidae) in Türkiye

Purpose Mosquitoes are important vectors of pathogens that can affect humans and animals. Culex tritaeniorhynchus is an important vector of arboviruses such as Japanese encephalitis virus, West Nile virus among various human and animal communities. These diseases are of major public health concern and can have huge economic and health burdens in prevalent countries. Although populations of this important mosquito species have been detected in the Mediterranean and Aegean regions of Türkiye; little is known about its population structure. Our study is to examine the population genetics and genetic composition of Cx. tritaeniorhynchus mosquitoes collected from several localities using cytochrome oxidase subunit I (COI) and the NADH dehydrogenase subunit 5 genes (ND5). This is the first extensive study of Cx. tritaeniorhynchus in the mainland Türkiye with sampling spanning many of provinces. Methods In this study, DNA extraction, amplification of mitochondrial COI and ND5 genes and population genetic analyses were performed on ten geographic populations of Culex tritaeniorhynchus in the Aegean and Mediterranean region of Türkiye. Results Between 2019 and 2020, 96 samples were collected from 10 geographic populations in the Aegean and Mediterranean regions; they were molecularly analyzed and 139 sequences (50 sequence for COI and 89 sequence for ND5) were used to determine the population structure and genetic diversity. For ND5 gene region, the samples produced 24 haplotypes derived from 15 variable sites and for COI gene region, 43 haplotypes were derived from 17 variable sites. The haplotype for both gene regions was higher than nucleotide diversity. Haplotype phylogeny revealed two groups present in all populations. AMOVA test results show that the geographical populations were the same for all gene regions. Results suggest that Cx. tritaeniorhynchus is a native population in Türkiye, the species is progressing towards speciation and there is no genetic differentiation between provinces and regions. Conclusion This study provides useful information on the molecular identifcation and genetic diversity of Cx. tritaeniorhynchus; these results are important to improve mosquito control programs.


Introduction
Culex is an important and diverse genus that encompasses more than 770 species grouped into twenty-six subgenera, several of which are of major public health concern [1].Culex tritaeniorhynchus Giles is a cosmopolitan species widely distributed in Africa, Southeast Asia, Middle East fever, West Nile virus (WNV) and Tembusu virus among various human and animal communities [12].These Culextransmitted disease can leave huge economic and health burdens in prevalent countries.WNV and Japanese encephalitis (JE) are one of the serious transmitted infections.They are caused by RNA viruses in the Flaviviridae family and similarly affect the central nervous system of host.WNV is prevalent in Africa, Europe, the Middle East, North America and West Asia [13,14].There have been reports of WNV cases and outbreaks in several European countries with the largest occurring in Israel, Greece, Romania, and Russia; between 2010 and 2022.There has been over 5,800 reported human cases and 378 deaths reported [15,16].JE, in contrast, is endemic in rural agricultural areas in the East Asia, with an estimated 70,000 cases and a mortality rate of 25-30% reported annually.Children aged ≤ 12 years are the most affected group [17,18].Although no JEV has been reported in Türkiye, there have been outbreaks of WNV; 47 WNV cases with 10 fatalities in 2010 [19].In the following period, 7, 23 and 10 WNV cases were reported in Türkiye in 2017, 2018 and 2019, respectively [20].Therefore, effective vector control, early detection of new vector species and ongoing vector surveillance are very important in combating these arboviral diseases.The success of vector control efforts to accurately identify the targeted mosquito species, depends on a correct understanding of their biology and ecology.Additionally, many taxonomists evaluate behavioral and population biology data together in identifying and classifying a species [21].
In recent years, numerous molecular phylogeny and population genetics studies have estimated the patterns of relatedness among and between important mosquito species [22][23][24][25][26]. Molecular markers have been developed to distinguish and identify mosquito species, replacing the limitations of morphological identification, which can be challenging for closely related species and requires welltrained technicians [27].DNA barcoding markers are frequently used for mosquito species differentiation.These studies are important as they provide valuable information on past biogeographic events and potential associations with life history traits of insect vectors.In Türkiye several studies have looked at the population structure of mosquito species such as Aedes zammitii [28], Ae. phoeniciae [29], Ae. albopictus [30], Ae. aegypti [31], Cx. pipiens [32].Most recently, Gunay et al. [7] investigated the taxonomic status of the Culex fauna in Türkiye using both morphological and DNA barcoding techniques.This work successfully determined or confirmed the presence of Cx. tritaeniorhynchus and 15 other Culex species, including four newly discovered species; however, no reports exist on the genetic structure of Cx. tritaeniorhynchus.Population genetic studies are important for the development of vector control measures, especially genetic control, to prevent or reduce the effects of epidemic diseases.mtCOI DNA barcode regions are commonly used to distinguish mosquito species [33,34].Ergunay et al. [35] explored the prevalence and distribution of WNV in potential vectors (mosquitoes) and diverse animal species across 15 Turkish provinces between 2011 and 2013.Notably, WNV RNA was detected in a limited number of mosquito pools, including Aedes caspius, Cx pipiens, Cx. quinquefasciatus, and Cx.perexiguus.Interestingly, all identified WNV strains belonged to lineage 1 clade 1a.Other field studies in Anatolia and Thrace regions also found WNV lineage 1 in similar mosquito species and migratory birds [36][37][38].These findings suggest widespread WNV activity across the investigated regions of Türkiye.
Mitochondrial DNA (mtDNA) is commonly used in insect ecology research.Its benefits include ease of use and specific features.Several studies worldwide have explored the population genetics of Cx. tritaeniorhynchus [10,12,39,40].These studies have reported high haplotype diversity and low nucleotide diversity in different Cx.tritaeniorhynchus populations (mainland China), which are attributed to factors such as rapid population expansion, genetic variation, and bottleneck effects [41,42].Additionally, it is considered that Cx. tritaeniorhynchus can cover distances (~ 2-7.5 km), so there might be more gene flow among populations [42].Elevation is considered an important factor affecting the distribution of this species Cx. tritaeniorhynchus [12].A recent study in Cambodia employed COI sequences to investigate the phylogeny and spatial distribution of Cx. tritaeniorhynchus [43].Similarly, a study in Thailand used geometric morphometrics and DNA barcoding on Culex vishnui Subgroup, which includes Cx. tritaeniorhynchus, to inform effective vector control of JE in the region [44].While COI remains the gold standard for barcoding and species identification [45,46], incorporating additional markers like ND5 can provide valuable complementary information about genetic diversity and population structure within a species and can offer insights into evolutionary rates or neutrality tests [47,48].This proves particularly beneficial when studying under-investigated regions or species with limited COI data.However, we acknowledge a limitation: the absence of ND5 data from other geographic regions like Cambodia and Thailand [43,44], hinders direct comparisons with those populations.To address this, we propose that future studies in this region consider including ND5 analysis alongside COI to facilitate broader comparisons and a more comprehensive understanding of the species' distribution and evolutionary history.To our knowledge, no detailed study has been conducted on the population genetics of the Cx.tritaeniorhynchus populations in Türkiye.Additionally, no studies have investigated the population genetics of Cx. tritaeniorhynchus using the mtND5 gene region.
This study aims to be the first to investigate the population structure and genetic variation of Cx. tritaeniorhynchus populations in the Aegean and Mediterranean regions of Türkiye using both mtCOI and mtND5, providing valuable information on the population dynamics of this species.

Specimen Collection
Cx. tritaeniorhynchus samples were collected from different districts between the months of May and September, over the 2019-2020 period (Table 1; Fig. 1).Ten geographic populations were selected from all known areas of this species has been detected in.CDC-light traps placed outside barns overnight and mouth aspirators were used to collect adults flying in and around outbuildings.
The adult specimens were transported to the laboratory in paper cups covered with gauze and were used in systematic evaluations under stereomicroscopes and identification key [49].Identified samples were labelled and stored within 95% ethanol at 20 o C prior to DNA extraction.

DNA Extraction, Amplification, and Sequencing
Total DNA extraction of individuals belonging to each population was performed using Invitrogen PureLink genomic tritaeniorhynchus populations were assessed using DnaSP 5.0 software [54] for each gene region.Nucleotide diversity (π) and haplotype diversity (Hd) are two major indicators to calculate the diversity of species populations among different geographical strains.Pairwise FSTs values of the genetic differentiation among the populations, selective neutrality tests including Fu'sFS statistics and Tajima's D tests, and analysis of molecular variance analysis (AMOVA) were performed in Arlequin 3.5 software [55].Median-joining Network analysis to show the relationship between haplotypes were performed using Network 10.2 software [56].Also, all p values were corrected using Holm's correction method.A phylogenetic tree was constructed using the neighbor-joining algorithm in MEGA7 software.Culex pipiens mtCOI and mtND5 sequences from GenBank database was used as an outgroup.The analysis was run on 1000 replicates for inferred bootstrap consensus.The best fit model was selected using MODEL TEST 3.0 software [57,58].

Results
Culex tritaeniorhynchus samples were obtained from 10 geographic populations after sampling mosquitoes in the Aegean and Mediterranean regions.DNA from the mtCOI and mtND5 gene regions from a total of 96 Cx. tritaeniorhynchus samples were amplified and after sequence analysis, only 50 mtCOI sequence data and 89 mtND5 sequence data were suitable and used for phylogenetic studies (Haplotypes' GenBank accession numbers for mtCOI: PP447202-PP447244; for mtND5: PP459986-460,009).
According to the data from 89 mtND5 (410 bp) nucleotide sequences, the G + C ratio of these sequences was calculated as 0.240 and found lower than the A + T ratio.All sequences generated 24 haplotypes generating from 15 polymorphic sites, and Hd and pi values were determined as 0.841 and 0.01783, respectively.The highest number of haplotypes were determined in the Antalya population (h = 9, Hd = 0.886 and pi = 0.01681), the lowest haplotype number was determined in the Mersin population (h = 6, Hd = 0.889 and pi = 0.01857).Tajima's D and Fu's FS tests were conducted for neutrality analysis to detect the history of geographical population size; Fu's FS tests were positive in the majority of populations (except Antalya population), the values were found to be statistically insignificant.Tajima's D test results were positive for all populations and this value was statistically significant for Osmaniye, Hatay and Adana populations.The genetic diversity values of the mtND5 gene region are summarized in Table 2.
For mtCOI gene region, 50 (final edited product size 450 bp) nucleotide sequences were used for haplotype analysis.The G + C ratio of these sequences was ranged DNA isolation kit as instructed by manufacturer.The total number of samples collected, and the number used for molecular studies is given in Table 2. Extracted DNA was stored at − 20°C until PCR amplification.Two mitochondrial gene regions: COI and ND5 were used for molecular studies as described by Folmer et al. [50] and Birungi and Munstermann [51].The universal LCO1490 and HC02198 barcoding primers were used to amplify the barcode region of the mtCOI gene (658-bp after primer removal), and PCR amplification conditions comprised initial denaturing at 95°C for 5 min, then 34 cycles of 95°C for 30 s, 48°C for 30 min and 72°C for 45 s, followed by a 5-min extension at 72°C and a 4°C hold.The mtND5 gene (approximately 450 bp) was amplified using primers ND5-F and ND5-R.Reactions comprised initial denaturing at 98°C, 2 min, then 30 cycles of 72°C, 45°C, 30 s, 72°C, 45 s, followed by a 5-min extension at 72°C and a 4°C hold.PCR master mix was prepared using Ampliqon Taq 2x master mix according to the manufacturer's instructions with a final volume of 25 µL each.The PCR reactions contained 12.5 µl PCR mix, 0.25 µl 20 µM each of primers, 1 µl template DNA and 11 µl ddH 2 O.The PCR reactions were performed using a T100 Thermal Cycler (Bio-Rad, Hercules, CA).Subsequently, amplified fragments were loaded on 1% agarose gel, stained with Safeview™ and visualized under UV light.Visualized products were purified by PCR purification kit and the cleaned PCR products sent out to Macrogen, Amsterdam, Netherlands for sequencing.

Population Genetics, Differentiation and Data Analyses
The ND5 and COI sequences obtained using forward and reverse PCR primers were visualized and edited using the BioEdit ver 7.0.9software [52] and MEGA7 software [53], respectively.Sequences were purified based on high-quality read peaks for each gene region using MEGA7and aligned using the Clustal W algorithm.To minimize positional dissimilarities between sequences of varying lengths, the alignments were manually adjusted after trimming the start and end points.All missing data and gaps within the sequences were removed.The sequence arrangement was then manually reassessed in an attempt to minimize the positional dissimilarity.All missing data and gaps within the sequences were removed.The sequences were compared within themselves and with the previous DNA sequences obtained from the Genbank, and sequences with more than 97% similarity were used in the study.Multiple sequence alignment of sequences was performed using the Clustal W algorithm in MEGA7 software [53].

Genetic Relationships among Geographical Populations in Türkiye
Genetic relationships among 6 geographical populations were determined using pairwise Fst distance.Fst distance values ranged between − 0,06268 and 0,06835 for ND5 gene region.The Fst values for all population pairs were detected as insignificant for mtND5 gene region (Table 3).
Fst distance values ranged between − 01980 and 0.29660 for mt mtCOI gene region.While between Osmaniye-Mugla, Hatay-Adana, Adana-Antalya population pairs, significant differences were determined between the populations, the differences between other populations were found to be insignificant (p > 0.05) (Table 4).between 0.240 and 0.315.All sequences generated 43 haplotypes originating from 17 polymorphic sites.Calculated Hd and Pi values were 0,992 and 0,01433, respectively.The highest haplotype number was determined in Antalya (h = 14) similar as in the ND5 gene region, while the lowest haplotype number was determined in Adana (h = 5) population.Fu's FS tests for neutrality tests were negative in most populations (except Adana and Mugla population), these values were found to be statistically significant (p < 0,05).Tajima's D test results were positive for all populations and this value was statistically significant like mtND5 gene region.Between populations, although the values were positive, they were not significant (p > 0,05) and the Tajima's D value could not be calculated for Adana and Mugla populations.The genetic diversity values of the mtCOI gene region are summarized in Table 2. Türkiye were also represented in these two groups.All the samples outside of Türkiye were in Cgroup 1. Also, a phylogenetic tree was constructed using concatenated haplotype from each population.The sequences concatenated in Sequence Matrix 1.7 software.The compatibility of the sequences was assessed with using a partition homogeneity test in PAUP 4.0 software (Fig. 4).

Haplotype Network Analyses
The relationship between haplotypes were analyzed using the median-joining method in Network 10.2 software.In the network, each circle represents one haplotype, and the circle size is related to the frequency of occurrence of the haplotype.The haplotype network was made with 43 haplotypes for the mtCOI gene region which were clustered into 2 main groups, supporting the NJ tree of mtND5 gene region.Osmaniye, Hatay, Antalya and Mugla populations, included in CGroup 1, was represented by all populations except Adana population.CGroup2 comprised of all populations.Hap 17, Hap 20, Hap 27 was represented higher than other haplotype in CGroup 1 for CGroup 3, the main haplotype was determined as Hap5, and Hap27, Hap17, Hap29 were the other highest frequency haplotypes (Fig. 5).

Phylogenetic Analysis of Culex tritaeniorhynchus Haplotypes
For the mtND5 gene region, Cx. tritaeniorhynchus samples from this study, two samples from China (GenBank number: KT851544, NC_028616) and a Culex pipiens (Gen-Bank number: MN389459) as an outgroup were obtained from the GenBank and used in obtaining a phylogenetic Neighbor-joining (NJ) tree using the T92 + G + I method.NJ tree revealed two main groups.The first group (Expressed as NGroup1 on the Fig. 2 For the mtCOI gene region, some of the haplotypes obtained showed high similarity with the GenBank samples.However, none of GenBank samples idented 100% with our samples.Phylogenetic tree was constructed using the samples with high blasted GenBank samples and haplotypes, 2 main groups emerged similar to the mtND5 gene region (Fig. 3).Moreover, GenBank samples recorded from

Analysis of Molecular Variance (AMOVA) Analysis
We performed the AMOVA test to determine whether there is a difference between the Aegean and Mediterranean populations of Türkiye.For mtCOI gene region, variance component proportion of within populations is high (92.41%),while among geographical strains (Aegean and Mediterranean) and within groups proportion is low (4.05% and 3.54% respectively).F-statistics, among groups, among populations within groups, and within populations were Fsc = 0.03688, Fst = 0.07593, and Fct = 0.04055, respectively, and were found statistically insignificant (Table 5).
AMOVA analysis for the mtND5 gene region gave similar results to the mtCOI gene region; variance component proportion of within populations is high (101.38%),while among geographical strains (Aegean and Mediterranean) and within group proportion is low (-0.20% and − 1.18% respectively).F-statistics, among groups, among populations within groups, and within populations were Fsc = 0.01174, Fst = 0.01378, and Fst = 0.00202, respectively, and were found statistically insignificant (Table 5).

Discussion
Molecular phylogeny is used to study various aspects of evolution, including gene duplication rates, diversification patterns, genetic variation, recombination events, population dynamics, and organismal relationships.This approach often involves combining molecular data with other sources of information [21].There have been relatively few comparative studies on Cx. tritaeniorhynchus population diversity, or variation in vectorial capacity, with the majority characterizing Cx. tritaeniorhynchus populations within Asia [40,41,[59][60][61][62]. Historically this is where this species has had the greatest occurrence, abundance, and caused the greatest impact on human health as the major vector of JEV [63].Studies on this species are limited in Türkiye.The mtCOI gene has been the most frequently used for previous studies of Cx. tritaeniorhynchus [40,41,[59][60][61]64] given it contains areas of highly conserved sequence, in combination with sufficiently diverse areas, allowing species discrimination and investigation of maternal inheritance patterns [22,64,65].As a first, our study examined the population genetics and genetic composition of Turkish Cx. tritaeniorhynchus mosquitoes collected from all known areas this species has been detected in using mtCOI and mtND5 sequences.
Results on nucleotide diversity and haplotype diversity, two major indicators to calculate the diversity of species populations among different geographical strains, showed that there was a high genetic variation detected in Cx. tritaeniorhynchus populations in Türkiye.Similarly, Rajavel et According to the mtND5 gene region, haplotype network was supported NJ tree results and the haplotypes was clustered on two main groups (NGroup 1 and NGroup 2).NGroup 1 and NGroup 2 were represented by all populations.In NGgroup 1, Hap 6 was determined as the main haplotype, and Hap 3 and Hap 1 was detected other high frequency haplotypes.For NGroup 2, the main haplotype was determined as Hap 2 (Fig. 6).and these values were statistically significant for Osmaniye, Hatay and Adana populations [68,69].Fu's FS statistical values were similarly positive except for the Antalya population.The Antalya population may indicate a re-expansion, possibly due to frequent control studies against mosquito species, because both the mtCOI and mtND5 gene regions gave negative and significant Fu's FS statistical values.On the other hand, Tajima's D values for the mtCOI gene region were positive, as in the mtND5 gene region, but all of the values were statistically insignificant.Fu's FS statistical values was negative in all populations (except Mugla and Adana populations) in other populations, and these values were statistically significant.This is evidence of an excess of low-frequency haplotypes as expected from population growth or may result from secondary contact between prior allopatric populations [70].The different results observed in Tajima's D and Fu's FS statistics may possibly be due to the previous observation of variation in the mtCOI gene and, accordingly, the insufficient number of samples (Table 2).Both NJ trees and Haplotype network trees formed two different groups and all groups were represented by all populations.This situation may have been caused by allopatric al. [61] and Airi et al. [66] in India, Jeffries et al. [67] in Asian, African and Middle Eastern populations determined high haplotype diversity (~ 0,9−1); while relatively low haplotype diversity (between ~ 0.5-0.8)was determined in the European continent where the species is newly entered.In another study done by Jeffries et al. [67], our results are not compatible with them, the reason for this situation may be due to the fact that not enough samples were used in the study (only 13 samples).
Previous studies have identified 28 [64], 303 [40] and 444 haplotypes [67].Analysis of regional population groups identified haplotypes in Asia  was determined in the regions where it is native [40,41,67].While high gene flux was determined in the populations for the mtND5 gene region among the populations, significant differences were determined between the Osmaniye, Mugla and Adana Antalya pairs for the mtCOI gene region.In fact, this is quite logical given the geographical distances.On the other hand, it is not possible to explain the significant speciation and could be shown to cause negative Fu' FS values for the mtCOI gene region.However, neutrality tests obtained from the mtND5 gene region did not support this situation.
In general, high variation was detected in all populations of the species.On the other hand, high variation is frequently seen in native species [32].In other studies, high variation is likely to provide more opportunities for range expansion and to increase gene exchange among different populations [40].
Despite the presence of Cx. tritaeniorhynchus first being reported in Europe-specifically Albania-in 1960 [77], further published European occurrence reports were scarce until the 2000s, with the species recorded in Greece; including from coastal marsh in Marathon near Athens [72], rice fields in Messolonghi, western Greece [11] and an urban area in Epirus, northwestern Greece [78].Recent extensive entomological surveys carried out in Albania have identified the presence of Cx. tritaeniorhynchus within multiple areas across the country.Concurrently, there has been an increasing trend of incursion, outbreaks and circulation of mosquito-borne arboviruses such as WNV in Europe, with many becoming established and endemic in multiple countries [15,79].
Our study is the first to assess the distribution pattern and perform molecular screening of Cx. tritaeniorhynchus populations in the Aegean and Mediterranean regions of Türkiye.In this study, distribution patterns of Cx. tritaeniorhynchus were identified, and genetic characterization of the mosquitoes is reported for the first time in Türkiye.The study targeted the Aegean and Mediterranean regions due to their contrasting geographies: the Aegean's gentle landscape (plains, hills, and valleys) versus the Mediterranean's rugged terrain dominated by the Taurus and Amanos Mountains.While both regions share a Mediterranean climate, the Aegean is cooler and wetter.This geographical and environmental variations can influence mosquito breeding and distribution.Also the Taurus Mountains might act as a barrier, promoting genetic differentiation between Cx. tritaeniorhynchus mosquito populations in these regions.These findings will enhance our understanding of vector-pathogen dynamics in Türkiye and inform researchers about the possibility of the occurrence of a possible outbreak.It may also facilitate control efforts in many aspects.Türkiye is the bridge between Asia and Europe or Middle East to Eastern part of Europe and many refugees originating from arboviral disease endemic areas use this route to reach continental Europe.The presence and spread of Cx. tritaeniorhynchus in Türkiye is an important public health difference between Hatay-Adana couples with geographical distance.The reason for this situation may be the geographical barrier created by Taurus and Amanos mountains.In the studies, it has been determined that the mtCOI gene region is more variant and mores ensitive than the mtND5 gene region [30].Phylogenetic tree and haplotype networks formed two main groups for mtND5 and a third intermediate group was determined between these two main groups.The fact that these groups are represented by all populations is probably due to gene flux between populations of the species, both actively and passively, which causes gene flux between populations.Almost the same situation has been identified for the mtCOI gene region by Xie et al. [40] in Chinese populations.This grouping may be geared towards the speciation of the species.These results were similar in the mtCOI gene region line [7] studied with samples collected from nearby geographical regions (Doğankent Adana), and it was stated in their studies that this situation was related to speciation.This grouping may be important in terms of contributing to the vector competence of the species [7].Seeing these variations within the species necessitates research on the basis of subspecies.The effect of geographical regions that may cause this situation was investigated with the AMOVA test.However, no significant differences could be determined between the groups for both the mtCOI and mtND5 gene regions.
There have been numerous recent reports of Cx. tritaeniorhynchus in countries where it had previously not been reported [71][72][73][74] highlighting a trend towards expansion of its known geographical range.Determining the phylogenetic origins of maternal lineages of Cx. tritaeniorhynchus can provide some insight into possible movement patterns when compared across countries and regions.During normal daily activity, Cx. tritaeniorhynchus are estimated to have an average flight distance of just under 70 m, however, some studies have found that during long-distance wind-assisted dispersal, they are estimated to migrate between 200 and 500 km [75].The adults over winter and it is thought this species may use a combination of long-distance migration and hibernation in situ, as strategies to survive unfavorable conditions in temperate regions [76].The ability to disperse over such long distances and adapt to variable conditions threat not only for the introduction of emerging arboviruses but also for autochthonous West Nile virus circulation.In these aspects, the introduction of national entomological surveillance programs and applying efficient vector control measures are critical factors in reducing the risk of future autochthonous arboviral disease.

Conclusions
To our knowledge this is the first published study about genetic diversity of Cx. tritaeniorhynchus populations in Türkiye.Availability of a broader range of genetic data, with wide coverage of informative genes will be valuable in further understanding the phylogeography, divergence, range expansions and evolution of this species.Also, Cx. tritaeniorhynchus populations will greatly expand the utility for comparison and potential for understanding this mosquito species and its contribution to vector-borne disease transmission.
This study also adds valuable information about this mosquito species acting as a significant vector for Japanese encephalitis in various Asian continents.The generated mtCOI sequences could be used as reference nucleotide sequences of the respective haplotypes in future mosquito identification studies and will facilitate the conspecific comparison to reveal the appropriate reason for high intra-specific divergence.Furthermore, different climatic conditions and geographical barriers can be responsible for the genetic variations in as pecies.This information is further used for the effective implementation of region-specific vector control strategies.

Fig. 2
Fig. 2 Neighbor-joining tree of Culex tritaeniorhynchus haplotypes, Genbank samples and outgroup using T92 + G + I distance model for mtND5 gene region.Bootstrap values were shown under the respective nodes

Fig. 5
Fig. 5 Network of the mtCOI haplotypes of Culex tritaeniorhynchus

Fig. 6
Fig. 6 Network of the mtND5 haplotypes of Culex tritaeniorhynchus

Table 1
Mosquito sampling localities and sample numbers

Table 2
Genetic diversity of Culex tritaeniorhynchus, using mtCOI and mtND5 gene sequences Number of samples; S = Number of variable sites; H = Number of haplotypes; Hd = Haplotypes diversity; p = Nucleotide diversity

Table 3
Pairwise Fst's of 6 populations using partial sequences of mtND5.Statistically significant values were marked in bold

Table 4
Pairwise Fst's of 6 populations using partial sequences of mtCOI.Statistically significant values were marked in bold

Table 5
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.Alten B, Simsek F, Aldemir A, Linton YM (2015)Barcoding Turkish Culex mosquitoes to facilitate arbovirus vector incrimination studies reveals hidden diversity and new potential vectors.Acta Trop 143:112-120.https://doi.org/10.1016/j.actatropica.2014.10.013 8. Touray M, Bakirci S, Ulug D, Gulsen SH, Cimen H, Yavasoglu SI, Hazir S (2023) Arthropod vectors of disease agents: their role in public and veterinary health in Turkiye and their control measures.Acta Trop 106893.https://doi.org/10.1016/j.actatropica.2023.1068939. Jupp PG, Kemp A, Grobbelaar A, Leman P, Burt FJ, Alahmed M (2002) The epidemic of Rift Valley fever in Saudi Arabia: Mosquito vector studies.Med Vet Entomol 16(3):245-252.https:// doi.org/10.1046/j.1365-2915.2002.00371.x10.Sallam MF, Al AAM, Abdel-Dayem MS, Abdullah MAR, Edward PR (2013) Ecological niche modeling and land cover risk areas for rift valley fever vector, Culex Tritaeniorhynchus Giles in Jazan, Saudi Arabia.PLoS ONE 8(6):e65786.https://doi.org/10.1371/journal.pone.006578611.Lytra I, Emmanouel N (2014) Study of Culex Tritaeniorhynchus and species composition of mosquitoes in a rice field in Greece.Acta Trop 134(1):66-71.https://doi.org/10.1016/j.actatropica.2014.02.018 12. Liu B, Gao X, Ma J, Jiao Z, Xiao J, Wang H (2018) Influence of host and environmental factors on the distribution of the Japanese encephalitis vector Culex Tritaeniorhynchus in China.IJERPH 15(9):184.https://doi.org/10.3390/ijerph15091848