Each mosquito was homogenized in 300µl of saline buffer. The homogenate was centrifuged at 800g for 2min, and the supernatant was retained. Total DNA was extracted from 100µl of supernatant using the RiboPrep kit (Amplisense, Moscow, Russia) according to the manufacturer's instructions. PCR was carried out in a 25μl volume with an Encyclo Plus PCR kit (Evrogene, Moscow, Russia). The COI gene fragments were amplified using primers 1454F and 2160R and 2027F and 2886R according to the original protocol (25 (link)). The amplification products were sequenced with forward and reverse primers for each primers pair. Sequencing of amplicons was performed using an ABI PRISM 310 sequencer and a BigDye Termination kit as recommended by Applied Biosystems (United States). The sequences were edited manually with Chromas and then aligned using MAFFT v7 (https://mafft.cbrc.jp). The total length of the COI gene sequences was 1317–1433bp. The nucleotide sequences were deposited in GenBank under accession numbers MZ 501500-MZ501561. Sequences MH817490–MH817558 from our preliminary study (38 ) were also included in the analysis. A total 131 sequences from individual mosquitoes were used for the genetic analysis.