4.1. Ethics statement
As the present study did not involve any human trial or higher vertebrates, the Institutional Animal Ethics Committee (IAEC) Department of Zoology, University of North Bengal (Regn. no. 840/GO/Re/S/04/CPCSEA) granted a waiver for ethics approval. The use of rat for blood feeding was also approved by the IAEC (approval no. IAEC/NBU/2019/19). All procedures were performed in accordance with relevant guidelines of the IAEC and ARRIVE (Animal Research: Reporting of In Vivo Experiments).
4.2. Study area and Sample collection
Dengue endemic Darjeeling district from northern part of West Bengal was surveyed and four different sites were selected for sampling. Larvae and pupae were collected from different breeding places and were transferred to plastic containers and brought to the laboratory. In the laboratory, mosquito larvae and pupae were reared up to F1 generation under controlled conditions (temperature 27°C±2°C; relative humidity 75%±10%). Standard identification key of larva and adult mosquitoes was used to identify the field population up to species level [33]. The sampling was done during June to November in 2020 and March to September in 2021. Since all the sampling was done from private land, prior permission was taken from the land owner for mosquito collection.
4.3. Insecticides used
Insecticide impregnated papers (4% DDT, 0.05% deltamethrin, 0.05% lambda-cyhalothrin, 0.75% permethrin) used for adult bioassay were purchased from Vector control unit, Universiti Sains Malaysia.
4.4. Insecticide susceptibility bioassay
The WHO (World Health Organization) adult bioassay protocol was followed for the detection of susceptibility status of the mosquito populations [34]. Seven replicates of 20 active 3-5 days non-blood fed female mosquitoes from each population were exposed to insecticide impregnated paper for an hour and were transferred to a retention tube containing cotton balls soaked in 10% sucrose solution and maintained at laboratory condition. For control, mosquitoes were placed in tubes containing paper impregnated with silicone oil. Mortality percentage was recorded after 24 hours post-exposure and was repeated thrice for every insecticide. In order to calculate the knockdown time (KDT), knocked down mosquitoes were counted after every 10 min during one hour insecticide exposure as per previous protocol [6]. The live and dead mosquitoes obtained from the adult bioassays were kept at -20°C and employed for DNA isolation.
4.5. Extraction of genomic DNA
Genomic DNA was extracted from both the resistant and susceptible mosquitoes following the High Salt protocol with minor modifications as described previously [35]. Individual mosquito was homogenized using digestion buffer and further incubated with proteinase K at 55-60°C for at least two hours. Afterwards, a 25:24:1 mixture of phenol, chloroform, and isoamyl alcohol was added to facilitate the partitioning of lipids and cellular debris into the organic phase, while isolated DNA remained in the aqueous phase. Following centrifugation, the aqueous phase containing the purified DNA was transferred to a clean tube for further analysis. Purity of the extracted DNA was checked by the SPECTROstar Nano fast scanning UV-Visible Microplate Reader (Make-BMG Labtech, Germany). DNA with an OD260/OD280 value between 1.8 - 2 was selected for kdr genotyping.
4.6. Allele-specific PCR (AS-PCR) assay for F1534C kdr mutation
DNA stock solutions were prepared at a concentration of 25 ng/μl and used for AS-PCR genotyping. The Polymerase chain reaction (PCR) involved one reverse primer: 5′-TCT GCT CGT TGA AGT TGT CGA T-3′ used for both alleles, and two forward allele-specific primers: 1534Phe: 5′-GCG GGC TCT ACT TTG TGT TCT TCA TCA TAT T-3′ and 1534Cys kdr allele: 5′-GCG GGC AGG GCG GCG GGG GCG GGG CCT CTA CTT TGT GTT CTT CAT CAT GTG-3′ with an annealing temperature of 60°C [13]. Each reaction was performed in a 25 μl volume consisting of 1.5 mM MgCl2, 1x PCR buffer (Promega, USA), 0.25 μM common reverse primer, 0.125 μM each mutation specific primer, 200 μM dNTP mixture (Promega, USA), 0.2 units Taq polymerase (Promega, USA) and 25 ng genomic DNA. The thermal cycling condition was set with an initial DNA denaturation step for two minutes at 94°C, followed by 35 cycles of denaturation for 30 sec at 94°C, annealing for 30 sec at mentioned temperature and extension at 30 sec at 72°C. PCR amplification products were loaded onto a 4% agarose gel and run for 1 hour at 100 V in TAE buffer and visualized by ethidium bromide staining under UV light. Since the primer used had GC tails of varying lengths, amplified products could be differentiated by base pair size.
4.7. Amplification and sequencing of IIIS6 vgsc gene of Ae. albopictus
PCR reaction was carried out for the partial amplification of IIIS6 vgsc gene. Primers used for these reactions were AaEx31P (5′-TCG CGG GAG GTA AGT TAT TG-3′,) and AaEx31Q (5′-GTT GAT GTG CGA TGG AAA TG-3′,) [13]. Reaction was carried out with 1X Go®Taq G2 Green Master Mix (Promega, USA) of 12.5μl, 1 μl of both forward and reverse primers, 2 μl of template DNA and 8.5 μl of nuclease free water in 25 μl reaction mixture. PCR condition was: 95 °C for 5 min, followed by 35 cycles of 95 °C for 30 s, 60 °C for 40 s and 72 °C for 1 min with a final extension step at 72 °C for 5 min. The amplified fragments of the expected size were purified using ExoSAP following manufacturer recommendations and were sequenced (Heredity Lifesciences Pvt. Ltd. Patia, Bhubaneswar-751024, Odisha, India). The sequences were analyzed with BioEdit software (v 7.0.9) and aligned with different homologous regions of vgsc gene sequences (KX371864, KX371865, and AB827824) of Ae. albopictus available in Gene Bank by using ClustalW software (v 2.0) [36,37].
4.8. Data Analysis
Mean mortality percentage against all the tested insecticides were calculated by using kyPlot 6.0. In WHO adult bioassays, control mortalities were below 10%, so no calculation of corrected mortality was needed. WHO 2016 criteria were followed to determine the resistance /susceptibility status [S=Susceptible (Mortality percentage=98-100%); R=Confirm Resistance (Mortality percentage < 90%); PR=Possible Resistance (Mortality percentage=90-98%)]. Knockdown times were determined by performing probit regression analysis in IBM SPSS (v21.0) at 95% confidence level. The web-based programme 'de FINETTI generator' version (v3.0.5) (2008) (https://finetti.meb.uni-bonn.de/) was used to compute genotype frequencies and their deviation from the Hardy-Weinberg equilibrium (HWE), which was shown within the de Finetti diagram. The diagram includes a triangular plot which represent the distribution of three genotypes in reference to one another. The curved line in the diagram represents the Hardy-Weinberg parabola that indicates the sites where alleles are in a state of HWE. The chi-square test is used to calculate the significance of the distance between the parabolic curve and the genotypes, which reflects the extent of divergence from the HWE.