Enrichment of phenotype among biological forms of Anopheles stephensi Liston through establishment of isofemale lines
Parasites & Vectors volume 16, Article number: 79 (2023)
Vector management programs rely on knowledge of the biology and genetic make-up of mosquitoes. Anopheles stephensi is a major invasive urban malaria vector, distributed throughout the Indian subcontinent and Middle East, and has recently been expanding its range in Africa. With the existence of three biological forms, distinctly identifiable based on the number of ridges on eggs and varying vectorial competence, An. stephensi is a perfect species for developing isofemale lines, which can be tested for insecticide susceptibility and vectorial competence of various biological forms.
We describe key steps involved in establishment and validation of isofemale lines. Isofemale colonies were further used for the characterization of insecticide susceptibility and differential vector competence. The results were statistically evaluated through descriptive and inferential statistics using Vassar Stat and Prism GraphPad software packages.
Through a meticulous selection process, we overcame an initial inbreeding depression and found no significant morphometric differences in wings and egg size between the parental and respective isofemale lines in later generations. IndCh and IndInt strains showed variations in resistance to different insecticides belonging to all four major classes. We observed a significant change in vectorial competence between the respective isofemale and parental lines.
Isofemale lines can be a valuable resource for characterizing and enhancing several genotypic and phenotypic traits. This is the first detailed report of the establishment of two isofemale lines of type and intermediate biological forms in Anopheles stephensi. The work encompasses characterization of fitness traits among two lines through a transgenerational study. Furthermore, isofemale colonies were established and used to characterize insecticide susceptibility and vector competence. The study provides valuable insights into differential susceptibility status of the parental and isofemale lines to different insecticides belonging to the same class. Corroborating an earlier hypothesis, we demonstrate the high vector competence of the type form relative to the intermediate form using homozygous lines. Using these lines, it is now possible to study host-parasite interactions and identify factors that might be responsible for altered susceptibility and increased vector competence in An. stephensi biological forms that would also pave the way for developing better vector management strategies.
The adaptive evolution of mosquitoes (considered as the deadliest animal) dates back to the time of dinosaurs from the Mesozoic to species in the Holocene era . Owing to the diversity of mosquitoes and their rapid rate of evolution [2,3,4], many research groups across the world work on various aspects of the mosquito vector including its biology, molecular genetics, physiology, population genetics, developmental biology, evolutionary biology, behavior and so on. During these studies, colonization of mosquitoes under laboratory conditions for the purpose of experimentation becomes essential. In recent years, several researchers have developed techniques to rear mosquito species that were previously considered difficult to maintain under laboratory conditions [5, 6]. One such effort led to the colonization of Anopheles stephensi under controlled environmental conditions. The rearing protocols with optimal temperature, light and humidity are established for standard insectary operations [7,8,9]. An. stephensi exhibits a greater magnitude of variation across populations, demanding the establishment of isofemale lines with greater genetic homogeneity, enabling the study of the quantitative and qualitative traits under laboratory conditions . The molecular variations at the genomic, transcriptomic and proteomic levels can be better established in isofemale lines . Isofemale lines have been previously used for studying behavior , morphometry , insecticide resistance [14, 15], pathogen susceptibility and speciation among mosquitoes and other insect species. By enrichment of phenotypes, isofemale lines enable the identification of associated genotypes, which are otherwise diluted in their heterozygous parent lines . Establishing an isofemale line has been shown to delineate the non-additive genetic variance due to local adaptations, bottleneck events and other epistatic effects .
Furthermore, when isofemale lines are maintained as isofemale colonies for many generations, they can be used for any live material experimental studies . During the initial stages of colonization, insect populations pass through a genetic bottleneck . Quantitative traits tend to vary greatly because of adaptation to the rearing conditions. However, since insectary operations are well standardized, maintaining the isofemale lines under these regimes would yield reliable the quantitative traits determined from the mean values. Though there are publications related to the establishment of isofemale lines, the techniques and success vary depending on the species (Drosophila , Trichogramma , Chilo sp.  and mosquitoes [13, 21]) depending on their biology. There are also many instances of unsuccessful establishment of isofemale lines due to improper mating, poor reproductive success, inbreeding depression, etc. . Among many hurdles, inbreeding depression is a daunting challenge that can be mitigated with careful line establishment practices and selection of the founder females in the filial generation with key fitness traits [13, 20, 23, 24]. In the present article, we outline a detailed procedure for the establishment of two isofemale lines from insectary colonized populations—a type form “IndCh” derived from TIGS-2 (T2) and an intermediate form “IndInt” derived from TIGS-6 (T6) . We report the characterization of the following parameters: (i) fitness, (ii) homozygosity across filial generations, (iii) insecticide susceptibility and (iv) vectorial competence for both these lines.
Materials and methods
Collection and colonization of mosquitoes
Mosquitoes were collected from their natural habitats through larval or adult sampling as per the previously established protocols . The T6 strain was originally collected from the semi-urban area of Sriramanahalli, Bangalore, Karnataka (12.972442°N, 77.580643°E), at the larval stage in 2016. The T2 strain was collected from the urban area of Anna Nagar, Chennai, Tamil Nadu (13.018410°N, 80.223068°E), at the adult stage in 2016. The rearing was standardized as per previously established protocols . The adults were identified using the identification keys of Nagpal and Sharma . On day 6 or 7, the adults were blood-fed through a standard membrane feeding system (approval ref. no. TIGS 2nd IBSC Oct 2018). For T2 populations, the resting blood-fed adults were collected inside concrete houses. The blood-fed adults were spotted in the corners of the room, resting on dark curtains. Around 60 adults were collected with the help of aspirators between 6 and 8 a.m. during summer. T6 populations were collected in fat larval stages from fresh water in new open cement tanks, primarily used to provide water for cows. Underground bore well water is used to replenish these tanks on a weekly basis. Around 250 larvae were collected from two nearby tanks.
Establishment of isofemale lines
The isofemale progenies were developed based on a method previously described by Ghosh and Shetty . About 20–25 gravid females (G0) were separated from their parental colonies (T2 and T6) maintained in the TIGS insectary. On day 3, each gravid female was transferred carefully to a single ovicup for egg laying. The individual females were allowed to lay eggs undisturbed for ~ 24–48 h under laboratory conditions. The single G0 female that laid the most eggs and had the highest percent hatchability was selected for the establishment of isofemale lines. Eggs collected from the single female were checked for egg ridge number and allowed to hatch. Larvae were provided with larval food and reared to adults as described earlier. Emerged adult siblings (G1) were kept in a rearing cage and allowed to inbreed (sibling mating). Around 5–10 mated females were randomly separated for continuing the filial generations. The same protocol was followed in subsequent filial generations (20 generations for IndCh and 23 generations for IndInt)(Fig. 1). Key life table parameters like fecundity (number of eggs produced per female), egg morphology, percent egg hatchability, number of larvae and pupae, sex ratio (male:female), longevity of adult females (lifespan) and frequency of blood meals taken by a female were recorded in every generation to check the fitness of the isofemale lines.
Establishment of isofemale colony
Both IndCh and IndInt isofemale lines were maintained as two separate isofemale colonies after 20 generations and 23 generations, respectively, instead of single female progeny . These two lines are currently being maintained in TIGS insectary. In both the colonies, males and females were allowed to inbreed within the filial generation, blood-fed and allowed to lay eggs in masses. All the key fitness parameters were monitored and recorded. These two isofemale colonies were utilized for insecticide-susceptibility assays and vectorial competence studies with Plasmodium berghei and Plasmodium falciparum.
Different developmental stages of the mosquitoes were carefully measured for their morphometric characteristics. In the present study, we have measured the eggs and wings in parental and isofemale lines.
In total, five parameters—egg ridge number, shape, length and width of egg, and egg float—from 15 randomly selected individuals were analyzed. Eggs were placed on a wet filter paper and measured under the microscope with an ocular micrometer (Unilab GE-34, Binocular Research Microscope, India) and recorded. The egg shape, ridge numbers and floats were studied following earlier protocols [30, 31]. In isofemale lines, egg ridges were counted in every alternate generation.
Measurement of wing length
A total of 15 wing samples of both male and female mosquitoes were selected for wing measurement. The wings of male and female mosquitoes were carefully dissected 10 days after eclosion. The mosquitoes were anesthetized on CO2 pads and the right wing dissected under a microscope (Olympus, SZX2-ILLK, Germany). Additional care was taken to ensure the wings were not damaged or folded during mounting. The wings were measured from the distal end of the alula to the tip between R1 and R2 veins, excluding the fringe scale as they are normally curved or wither off during mounting . Wings were photographed and recorded under a stereomicroscope (Leica MZ10F, Germany). The photographed wings were measured using Leica Application Suite X (LAS X) software. The measurements were carried out independently by two observers to minimize error.
Estimation of homozygosity
The Fastq files for the three reference genomes were downloaded from NCBI SRA (sequence read archive). The three genomes (IndCh = SRR15146350, IndInt = SRR15603373, STE2 = SRR1168951) were mapped against the UCI 2.0 reference genome using Bowtie2 to call variants and generate a VCF file. These files were filtered for DP = 3, QUAL = 10 and Chr 2, 3, X and MAF (minimum allele frequency) = 0.01 and merged. The merged VCF file was converted to a matrix file containing numeric genotype information (0, 1, 2) representing the three possible genotypes for each SNP. The information was used to quantify heterozygosity-homozygosity in the genomes.
Five-day-old adult female mosquitoes of the parental and isofemale colonies were used for insecticide-susceptibility assays as per WHO guidelines. Insecticide-impregnated papers of discriminating doses and insecticide resistance monitoring test kits were ordered from Vector Control Research, Universiti Sains Malaysia (a WHO referral center for these kits and impregnated papers) .
Each replicate consisted of 25 females, and four such replicates were used for the experiment. One corresponding control was tested with the control papers supplied with the kit for respective classes of insecticides. They were exposed to eight insecticides representing all four major insecticide classes (Table 1). Mosquitoes used for DDT, permethrin and deltamethrin were pre-exposed to PBO (synergist) to understand the role of the metabolic mechanism of resistance. Assays were performed with 1 h exposure followed by 24 h recovery time, where the mosquitoes were provided with 10% glucose solution soaked on cotton pads. The mortality was recorded 24 h post-exposure, and all the mosquitoes (dead and alive tested) were kept in individual microfuge tubes and preserved at – 20 °C for future analysis. All the susceptibility assays were carried out at the TIGS insectary at 27 ± 1 °C and 75 ± 5% RH.
Vectorial competence studies
The vector competence studies were carried out in parental and isofemale colonies. The competence was evaluated in both in vivo and in vitro assays using the P. berghei and P. falciparum parasites.
Albino mice (BALB/c) were provided by the Animal Care Resource Center (ACRC) of the National Centre for Biological Sciences (NCBS), Bangalore, India. All animal procedures were approved and conducted according to the Instem IAEC under project INS-IAE-2019/01(N). The rodent malarial parasite P. berghei (ANKA strain-MRA-311, BEI Resources, USA) was used to infect mice by injection of 100–150 µl intraperitoneally per mouse . Infected mice, upon reaching 5–6% parasitemia and 1:1 or 1:2 male:female gametocyte ratio, were anesthetized and exposed to 5–6-day-old female mosquitoes as detailed below. Prior to keeping mice for feeding mosquitoes, the final counts of parasitemia and exflagellation were recorded .
After eclosion, adult female mosquitoes were provided with 10% glucose solution for 4–6 days. Forty mosquitoes in two replicate cups were kept fasting for 6 h prior to blood feeding. The mosquito cups were covered with a black cloth and left undisturbed for 30–40 min during blood feeding. The fully engorged females were separated from unfed and half-fed mosquitoes. The cups were maintained in a bioenvironmental chamber (Percival insect chamber I-30VL, Perry, IA, USA) at 19 °C and 75% RH (12:12 h day and night cycles) for 13–14 days. A cotton wool pad was soaked with 10% d-glucose; 0.05% para-aminobenzoic acid (PABA) solution was changed every alternate day until dissection. The mosquito midguts were dissected on the 14th day using PBS and stained with 1% mercurochrome. Midguts were examined for the presence of oocysts and recorded under a microscope (Nikon eclipse Ni-U upright microscope, Japan).
Plasmodium falciparum [PfNF54 Line E, procured from BEI resources (MRA-1000)] was revived for standard membrane feeding assay (SMFA). The culture was initiated at 0.15–0.2% asexual parasitemia and 5% hematocrit in 10 ml complete medium (RPMI-1640 with 6 g/l HEPES, 50 mg/l hypoxanthine, 2.5 g/l sodium bicarbonate and 10% human serum) (ethics approval ref. no. inStem/IEC-12/001, dated 20/03/2019). The cultures were maintained in an atmosphere of 5% O2, 5% CO2 and 90% N2 for 16–18 days with daily media change. For the mosquito feeding experiment, staggered cultures (14 and 17 days old) were selected based on their stage-V gametocytemia and exflagellation activities (above 10 ex-flagellating centers per 40× field) and pooled .
During the day of SMFA, 4–6-day-old female mosquitoes (40 nos., in duplicate) were placed in paper cups and maintained with 5% d-glucose overnight. Sugar cotton was removed 4 h prior to blood feeding, and mosquitoes were kept in the dark till SMFA . They were fed with mature stage-V P. falciparum NF54 gametocytes as mentioned earlier. The temperatures were maintained at 37 °C during SMFA. After 30 min of blood-feeding, fully engorged mosquitoes were maintained in an environmental chamber set at 26 °C and 75% RH (Percival I-30VL insect chamber, Perry, IA, USA). Mosquitoes were provided fresh wet cotton balls soaked in 10% d-glucose and 0.05% PABA solution daily. The midguts were dissected on the 9th day and stained with 1% mercurochrome, and the number of oocysts per midgut was recorded using a light microscope (Nikon Eclipse Ni-U upright microscope, Japan).
Descriptive and inferential statistical analyses were used in the present study. Various life table parameters such as fecundity, hatchability, pupation percentage, eclosion percentage and male:female ratio were recorded and are presented in graphical forms. For eggs and wing measurements, statistical t-test analysis for independent or correlated samples was performed using Vassar Stat software (http://vassarstats.net/). The mean and SEM were compared between populations. A p-value > 0.05 was considered non-significant for each parameter. Insecticide susceptibility was evaluated using the Prism GraphPad software package. The vectorial competence between the mosquito populations was analyzed using a non-parametric Mann-Whitney test using Prism GraphPad.
Transgenerational fitness during sibling mating
In each generation, several key life table parameters such as fecundity, percent egg hatchability, number of larvae and pupae, sex ratio (male:female), longevity of adult females (lifespan) and frequency of blood meal were recorded. A greater number of females was selected to start the isofemale line as we expected higher mortality during the initial filial generations due to inbreeding [37, 38]. The average number of eggs was around 80–100 per individual female, and fecundity of both the lines normalized around the 15th generation (Fig. 2A). A similar trend was also observed in hatching percentage (Fig. 2B). Post normalization, the hatchability of both isofemale lines improved to 85%. The development of larvae to pupae (pupation %) (Fig. 2C), followed by the transformation of the pupae to adults (eclosion %) (Fig. 2D), was affected during sibling mating. There were instances where sharp dips were observed owing to the mortality of the larvae and pupae. Most of the pupal mortality was observed during eclosion as also observed in other studies . The male:female ratio was calculated in each generation. The findings suggest (Fig. 3A, B) that the ratios tend to get skewed during the first few filial generations and balance out after successive generations (1:1.5) . There was no significant change in the longevity of adult females (lifespan) or in their blood-feeding behavior.
Morphometrics of eggs and wings
To ascertain the morphometric changes that arise during and after sibling mating, several parameters such as egg shape, size, float numbers, float ridge and wing length were observed and recorded (Table 2). The egg was slightly boat-shaped, black in color and blunt at the anterior and posterior ends (Fig. 4). The floats, which are an important morphological marker for differentiating the biological forms of different Anopheles species, showed minor variations in size. The maximum float lengths of IndCh and IndInt isofemale forms were 229.41 ± 4.81 and 209.33 ± 10.62 μm and the float widths 71.17 ± 6.96 and 69.33 ± 7.98 μm, respectively. The float length-to-width ratios of IndCh and IndInt were 3.24 ± 0.08 and 3.15 ± 0.15, respectively. The maximum float lengths of IndCh and IndInt parental strains were 218.09 ± 3.49 and 210.58 ± 6.21 μm and the float widths 71.17 ± 6.96 and 69.33 ± 7.98 μm, respectively. The float length-to-width ratios of IndCh and IndInt parental strains were 3.40 ± 0.09 and 3.14 ± 0.12, respectively (Table 2). The float size also influenced the number of ridges. For type and intermediate isofemale eggs, the number of float ridges was 20.81 ± 0.08 (range 20–21) and 16.66 ± 0.13 (range 16–17), respectively. For type and intermediate parental strains, it was 20.8 ± 0.11 (range 19–21) and 17.35 ± 0.31 (range 16–18). Many other parameters related to morphometry of eggs are presented in Table 2. Overall, there were no significant differences in all the parameters of eggs when compared between the isofemale lines of IndCh (F = 1.03, p > 0.05) and IndInt (F = 1.01, p > 0.05) and also when compared with their respective parental lines.
There was no significant difference between the male and female mosquito wing lengths (Table 3) and between the parental line and their isofemale forms. The wing measurement values were non-significant when compared between the isofemale lines of IndCh (F = 1.32, p > 0.05 for males; F = 2.43, p > 0.05 for females) and IndInt (F = 1.31, p > 0.05 for males; F = 1.12, p > 0.05 for females) and with their respective parental lines.
Homozygosity achieved in isofemale lines
To decipher the levels of homozygosity achieved through sibling mating across filial generations, we obtained 100× sequence coverage using Illumina reads from DNA pooled from 50 individuals from the isofemale IndCh and IndInt strains homogenized for 5 and 14 generations, respectively. As a negative control, we used the dataset generated from a pool of 50 + individuals directly from a laboratory colony  of the STE2 strain without genetic homogenization. The low ratios of heterozygous over homozygous SNPs seen for the two isofemale lines showed the decrease in heterozygosity achieved by isofemale lines around the fifth generation over pools of individuals from laboratory samples (Fig. 5).
It should be mentioned that while IndCh is homozygous for standard forms of all known inversions, STE2 is heterozygous for the 2Rb genotype and IndInt is heterozygous for the 3Li genotype even after 14 generations of isofemale selection . Interestingly, the isofemale line for IndInt shows increased heterozygosity for 3Li relative to its parent line (unpublished data), perhaps amplifying the reduction in vectorial competence in the isofemale over its parental line (T6, Figs. 7, 8). Even in the case of IndCh, the percentage of insects with the standard 2Rb increased over generations of isofemale lines , perhaps explaining the increased vectorial competence of the isofemale line over that of the parental line (T2, Figs. 7, 8).
Susceptibility assay across four classes of insecticides
Insecticide susceptibility was determined for both isofemales and their parental colonies across four classes of insecticides  (Fig. 6). The number of surviving individuals at 24 h post-exposure was calculated to understand the susceptibility/resistance status. The results suggest both the isofemale lines showed higher susceptibility to DDT than their parental lines. In case of dieldrin, IndCh showed a decreased mortality percentage of 71.6% compared to T2, which exhibits possible resistance. In contrast, the IndInt showed greater sensitivity (93% mortality) than T6 (68.2% mortality) as per WHO standards . The change in the susceptibility might be due to different target sites of these two insecticides. In organophosphates, all the populations showed high mortality (> 90%), except T6. For fenitrothion, both isofemale lines exhibited significant resistance (< 30% mortality) compared to their respective parental lines (non-parametric, Mann-Whitney test; p < 0.0001). For carbamate and synthetic pyrethroids, all four mosquito colonies showed confirmed resistance (mortality < 90%, range 60–87%) (Fig. 6).
Vectorial competence studies
Four colonies (T2, IndCh, T6, IndInt) were fed with P. berghei-infected blood of BALB/c mice and evaluated for their infection rate. The infection rate was calculated for each colony as the percentage of mosquitoes with oocysts in their midguts. There was not much variation in mean infection rates of the type form T2 (82.73%) and IndCh (85.07%) (Fig. 7). However, it was interesting to note there were considerable differences in the mean infection rates between the intermediate form of T6 (92.86%) and IndInt (70.83%).
Among the four colonies, the highest oocyst range was observed in T2 (1-1050), followed by IndCh (1-1003), T6 (1-870) and IndInt (1-150) (Fig. 7B). The mean number of oocysts per positive midgut in T2, IndCh, T6 and IndInt was 52.47 (SD ± SEM = 120.3 ± 10.20), 235.6 (266.4 ± 32.55), 112.4 (168.6 ± 17.03) and 14.54 (24.39 ± 2.87), respectively, whereas the median values were 10.00, 190.00, 50.50, 5.00 for T2, IndCh, T6 and IndInt, respectively. Furthermore, a comparison between parental colonies and isofemale colonies revealed a significant difference in their vector competence ( p < 0.0001, non-parametric, Mann-Whitney test).
In the case of P. falciparum, there was little variation in mean infection rates of the type forms, T2 (80.40%) and IndCh (78.57%) (Fig. 8B). The mean infection rates differed between the intermediate form of T6 (79.62%) and IndInt (69.57%). Among the four colonies, the highest oocyst range was observed in T6 (1-231), followed by IndInt (1-125), IndCh (1-106) and T2 (1-85). The mean number of oocysts per positive midgut in T2, IndCh, T6 and IndInt was 16.34 (SD ± SEM = 17.94 ± 1.13), 26.39 (34.70 ± 2.77), 18.39 (19.56 ± 1.65) and 13.30 (19.63 ± 1.67), respectively, whereas the median values were 12, 14, 16, 6 for T2, IndCh, T6 and IndInt, respectively.
Furthermore, a non-parametric Mann-Whitney comparative test between T2-T6 colonies ( p < 0.0563) and T2-IndCh colonies ( p < 0.3968) did not show any significant difference. However, the homozygous forms of the isofemale colonies, IndCh and IndInt showed a significant difference ( p < 0.0045). Significant difference was also observed between T6 and IndInt ( p < 0.0003). These results corroborate earlier findings that the type form is a better vector than the intermediate form [34, 41] (Fig. 8B).
The present study provides a detailed overview of the key considerations followed during the establishment of the isofemale line in the two biological forms of An. stephensi. We have meticulously characterized several life cycle parameters in every generation during the creation of both isofemale lines to maintain their authenticity. During establishment, key emphasis was given to restore some characteristics of the parent population after the initial bottleneck to achieve normalization of fecundity, hatching, pupation, eclosion and male:female ratio. Accordingly, we found that in both isofemale lines, homogeneity in these parameters was achieved within ~ 13 generations (Fig. 2).
Even though homogeneity in the biological parameters was achieved in ~ 13 generations, the impact of this on other traits, such as insecticide resistance and vectorial competence, was dramatic. This is the first time to our knowledge that a discrete study has been conducted by establishing two homozygous isofemale lines from the type and intermediate forms of An. stephensi and assessing their vectorial competence against two different malaria parasites, P. berghei and P. falciparum. The findings suggest that there is significant enrichment in vectorial competence between the parental and corresponding isofemale lines. When mosquitoes were exposed to P. berghei-infected mice, there was increased susceptibility to the type (IndCh) compared with intermediate form (IndInt). Furthermore, a significant difference in vectorial competence was seen between the parental and corresponding isofemale lines. This might be due to the high levels of homozygosity in certain alleles enhancing the vectorial competence.
For the P. falciparum-infected blood meal, there was no significant difference between the parental and corresponding isofemale lines. However, significant differences were found between the two isofemale lines belonging to the type and intermediate forms, supporting earlier reports of the type form being a better vector than the intermediate form . Establishing an isofemale line has been shown to delineate the non-additive genetic variance due to local adaptations, bottleneck events and other epistatic effects . We conclude that isofemale lines derived from type and intermediate forms achieve enrichment of characteristic phenotypes regarding vectorial competence expected for each form.
We also sequenced DNA from both IndCh and IndInt lines after 5 and 14 generations to assess the homogeneity achieved at the genome level. As a negative control, we used a public dataset with 100× coverage of STE2 strain. The ratio of the heterozygous versus homozygous variants in these reads relative to a gold-standard reference genome  clearly shows that genome level homozygosity is achieved within as few as five generations in IndCh. Furthermore, the increased homozygosity in the organism provides an opportunity to decipher the genome architecture at very high resolution [43, 44]. While the high-resolution genome assemblies of the UCI and IndCh strains were achieved after 5 and 6 generations of isofemale lines, the IndInt line was homogenized over 14 generations before attempting an assembly because, as shown in Fig. 2A, increased generations were needed to normalize fecundity. This could be because, unlike the UCI and IndCh strains that have no heterozygous paracentric inversions, IndInt retains a large paracentric inversion even after 14 generations. Furthermore, our data show (unpublished data) that the percentage of this inversion increased considerably with more generations relative to that of the parent population. This may be a notable example of the complexity of maintaining a given trait while retaining fitness.
As a comparative account between P. berghei and P. falciparum, we also observed higher infection rates in both the parental and isofemale lines for P. berghei model. However, we did not find a significant difference between the two parental lines T2 and T6 and also the parental line T2 and its isofemale form for P. falciparum infections. We speculate that there might be differential immune response that could have happened during homogenization, leading to the possibility of some anti-Plasmodial genes being differentially expressed. In our earlier work, we have described few structural variations (2Rb and 3Li inversions) and genes associated with vector competence during the genomic characterization of these isofemale lines, IndCh and IndInt [44, 59]. Since functional genomics was not in the scope of this study, we did not explore it further. However, probing this in greater detail would shed some light on delineating the specific host-parasite interactions.
In the present study, we observed significant variations in insecticide resistance between the founding population and its corresponding isofemale progeny. This might be a result of the segregation of the genotypes resulting in the better demarcation of the susceptible/resistant phenotype compared with the parental lines. This can be further studied for a better understanding of underlying mechanisms of resistance, which is relatively easier to identify in isofemale lines because of their homozygosity.
In all four colonies including the two parent and two isofemale lines, the levels of resistance varied across the two insecticides belonging to the same class, such as DDT and Dieldrin, belonging to the organochlorine class. For Dieldrin, the parental line T2 exhibited greater sensitivity than its isofemale IndCh, whereas the parental line T6 showed greater resistance than its isofemale form IndInt. The differences in the modes of action of these two insecticides might be the reason for this trend. For example, DDT resistance is associated with mutation(s) in the target site of the voltage-gated sodium channel (vgsc) gene and metabolic resistance [45, 46], whereas dieldrin resistance is due to a single point mutation in GABA receptors . Also, among the organophosphates, resistance was observed to fenitrothion, but susceptibility was seen for malathion among all four colonies. For fenitrothion, there was a significant difference in the susceptibility status between the parental and their respective isofemale lines. Earlier studies had suggested that many strains of An. stephensi collected from different parts of India were susceptible to fenitrothion . However, for carbamates and synthetic pyrethroids, all four populations showed resistance to both classes of the insecticides. It is interesting to note that although the insectary-colonized mosquitoes and their isofemale lines were not under any selection pressure from insecticides and were grown without any insecticide exposure, they continued to maintain this resistance.
The observed variations in insecticide susceptibility might be due to large structural variants such as inversions. Recent advances have revealed the possible role of structural variations like duplications and inversions in insecticide resistance, which would be interesting to dissect among isofemale colonies in the future . The role of structural variation has also been implicated in different mechanisms of copy number variation, transposable element insertions and tandem duplication of different metabolic genes [49, 50]. Characterization of the isofemale lines during their establishment as reference genomes revealed that there is a homozygous standard 2Rb form in the IndCh isofemale line  and a heterozygous 3Li inversion in the IndInt isofemale line . The functional characterization of the 3Li region revealed 36 isoforms of cuticle-forming genes in the IndInt isofemale line. Earlier studies suggest that cuticle thickening is associated with pyrethroid and insecticide resistance in other Anopheles species . Similarly, the 2Rb region has 1353 predicted genes, which include several genes associated with insecticide resistance such as ACE1 as well as tandem clusters of GST and Cyp450 paralogs in the IndCh isofemale line .
Variations in vector competence among the biological forms
Although the various biological forms can be clearly delineated using the number of ridges on the egg float, there have been efforts to develop molecular markers to identify the various forms. These include mitochondrial genes cytochrome oxidase 1 (COI) and cytochrome oxidase 2 (COII) , the rDNA-ITS2 and domain-3 (D3) marker , odorant-binding protein AnsteObp1 marker  and microsatellite markers  to characterize the different biological forms. A single SNP has been previously reported  to distinguish between the various forms of An. stephensi found in Iran. When we tested the validity of this SNP for Indian An. stephensi populations (insectary-colonized and wild populations) that were collected from five different locations, we found that most of the individuals from the wild had the SNP (G) in AsteObp1 (intron I, position 91 bp) reported to be diagnostic of the type form, 18/105 individuals had an A at this same position, 18/105 had an A/G and 1/105 had an A/C in that position, likely representing population-level variations that have not been reported earlier (unpublished data). Also, many of the other markers are not universal and may vary from region to region and cannot resolve the biological forms accurately . Here, we clearly demonstrate that the number of egg ridges remains the gold standard and applicable marker to develop isofemale lines until new molecular markers are identified.
There are also reports of a third form with the lowest egg-float ridge number known as ‘var. mysorensis’ [25, 31, 34, 57], but this has not been established or characterized in the current study. The mysorensis form is mostly found in rural areas and breeds in fresh water sources like wells, ponds, cement tanks, streams, etc. Earlier studies have shown var. mysorensis to exhibit variations in their insecticide susceptibility status  and vectorial competence . The mysorensis form is highly zoophilic and has very poor vectorial competence. A thorough characterization of this third biological form might reveal cues that can be exploited in vector management strategies.
Overall, the findings of this article show the benefit of developing isofemale lines from the insectary-colonized mosquito populations. Furthermore, we showed that as few as five filial generations of sibling mating may be sufficient in An. stephensi to achieve greater levels of homozygosity. However, more sibling matings may be required (> 13 generations) to establish isofemale lines with optimal fitness. The isofemale lines thus developed could be used to elucidate mechanisms underlying insecticide resistance, which is more pronounced than in the parental lines. The isofemale lines should also prove useful in understanding the basis for the greater vectorial competence of the type variant relative to the intermediate form for both P. berghei and P. falciparum. By assembling high-quality reference genomes for IndInt and IndCh, we have been able to hypothesize the genetic factors/structural variants that might account for the differential vectorial competence .
Availability of data and materials
All data generated or analyzed during this study are included in this published article.
Jacobson RL, Poinar G Jr, Poinar R. What bugged the dinosaurs? in insects, disease and death in the cretaceous. Parasites Vectors. 2008;1:6.
Escobar D, et al. Distribution and phylogenetic diversity of Anopheles species in malaria endemic areas of Honduras in an elimination setting. Parasit Vectors. 2020;13:333.
Tainchum K, et al. Anopheles species diversity and distribution of the malaria vectors of Thailand. Trends Parasitol. 2015;31:109–19.
Singh V, et al. Why is it important to study malaria epidemiology in India? Trends Parasitol. 2009;25:452–7.
Phasomkusolsil S, et al. Colonization and maintenance of anopheles belenrae and anopheles pullus from the Republic Of Korea. J Am Mosq Control Assoc. 2018;34:260–4.
Zerpa N, et al. Colonization and laboratory maintenance of anopheles albimanus wiedemann in venezuela. Rev Inst Med Trop Sao Paulo. 1998;40:173–6.
Khan I, Farid A, Zeb A. Development of inexpensive and globally available larval diet for rearing anopheles stephensi (Diptera: Culicidae) mosquitoes. Parasit Vectors. 2013;6:90.
Gerberg EJ, et al. Manual for mosquito rearing and experimental techniques. Lake Charles: American Mosquito Control Association Inc; 1994.
Nasirian H, Ladonni H. Artificial bloodfeeding of Anopheles stephensi on a membrane apparatus with human whole blood. J Am Mosq Control Assoc. 2006;22:54–6.
Charlesworth D, Willis JH. The genetics of inbreeding depression. Nat Rev Genet. 2009;10:783–96.
Madi-Ravazzi L, et al. Fluctuating asymmetry of meristic traits: an isofemale line analysis in an invasive drosophilid. Zaprionus indianus Genetica. 2017;145:307–17.
Chinchilla-Ramirez M, et al. Genetic variation in the feeding behavior of isofemale lines of nesidiocoris tenuis. Insects. 2020;11:513.
Jirakanjanakit N, Leemingsawat S, Dujardin JP. The geometry of the wing of Aedes (Stegomyia) aegypti in isofemale lines through successive generations. Infect Genet Evol. 2008;8:414–21.
Epelboin Y, et al. CYP450 core involvement in multiple resistance strains of Aedes aegypti from French Guiana highlighted by proteomics, molecular and biochemical studies. PLoS ONE. 2021;16:e0243992.
Witzig C, et al. Identifying permethrin resistance loci in malaria vectors by genetic mapping. Parasitology. 2013;140:1468–77.
Churchill FB. William Johannsen and the genotype concept. J Hist Biol. 1974;7:5–30.
David JR, et al. Isofemale lines in Drosophila: an empirical approach to quantitative trait analysis in natural populations. Heredity (Edinb). 2005;94:3–12.
Zhang C, et al. A new method to obtain isofemale lines for monitoring Chilo suppressalis resistance to Bt toxin. J Invertebr Pathol. 2019;163:8–10.
Lainhart W, et al. Changes in genetic diversity from field to laboratory during colonization of anopheles darlingi root (Diptera: Culicidae). Am J Trop Med Hyg. 2015;93:998–1001.
Coelho A Jr, et al. Laboratory performance predicts the success of field releases in inbred lines of the egg parasitoid trichogramma pretiosum (Hymenoptera: Trichogrammatidae). PLoS ONE. 2016;11:e0146153.
Tardieux I, et al. Analysis of inheritance of oral susceptibility of Aedes aegypti (Diptera: Culicidae) to dengue-2 virus using isofemale lines. J Med Entomol. 1991;28:518–21.
Bubliy OA, Loeschcke V, Imasheva AG. Genetic variation of morphological traits in Drosophila melanogaster under poor nutrition: isofemale lines and offspring–parent regression. Heredity (Edinb). 2001;86:363–9.
Dembo E, et al. A user-friendly method to assess Anopheles stephensi (Diptera: Culicidae) vector fitness: fecundity. J Med Entomol. 2014;51:831–6.
Yeap HL, et al. Body size and wing shape measurements as quality indicators of Aedes aegypti mosquitoes destined for field release. Am J Trop Med Hyg. 2013. https://doi.org/10.4269/ajtmh.12-0719.
Subbarao SK, et al. Egg-float ridge number in Anopheles stephensi: ecological variation and genetic analysis. Med Vet Entomol. 1987;1:265–71.
Tilak R, Dutta Gupta KK, Verma AK. Vector databank in the Indian armed forces. Med J Armed Forces India. 2008;64:36–9.
Ghosh SK, Ghosh C. New Ways to Tackle Malaria, in Vector-Borne Diseases - Recent Developments in Epidemiology and Control. London: IntechOpen; 2020.
Nagpal, B.N. and V.P. Sharma, Indian Anophelines. 1995: Science Publishers.
Ghosh C, Shetty NJ. Mode of inheritance of fenitrothion resistance in Anopheles stephensi Liston. J Cytol Genet. 1999;34:141–6.
Shetty N, et al. Egg float ridge number in Anopheles stephensi: ecological variation. J Parasit Dis. 1999;23:45–8.
Tyagi V, et al. Morphometric and morphological appraisal of the eggs of Anopheles stephensi (Diptera: Culicidae) from India. J Vector Borne Dis. 2017;54:151–6.
Wilke AB, et al. Morphometric wing characters as a tool for mosquito identification. PLoS ONE. 2016;11:e0161643.
Richards SL, et al. Assessing insecticide resistance in adult mosquitoes: perspectives on current methods. Environ Health Insights. 2020;14:1178630220952790.
Basseri HR, et al. Comparative susceptibility of different biological forms of Anopheles stephensi to Plasmodium berghei ANKA strain. PLoS ONE. 2013;8:e75413.
Patra KP, et al. Alga-produced malaria transmission-blocking vaccine candidate Pfs25 formulated with a human use-compatible potent adjuvant induces high-affinity antibodies that block Plasmodium falciparum infection of mosquitoes. Infect Immun. 2015;83:1799–808.
Habtewold T, et al. Streamlined SMFA and mosquito dark-feeding regime significantly improve malaria transmission-blocking assay robustness and sensitivity. Malar J. 2019;18:24.
Duffy E, et al. Inbreeding alters intersexual fitness correlations in Drosophila simulans. Ecol Evol. 2014;4:3330–8.
Ross PA, Endersby-Harshman NM, Hoffmann AA. A comprehensive assessment of inbreeding and laboratory adaptation in Aedes aegypti mosquitoes. Evol Appl. 2019;12:572–86.
Alto BW, Juliano SA. Precipitation and temperature effects on populations of Aedes albopictus (Diptera: Culicidae): implications for range expansion. J Med Entomol. 2001;38:646–56.
Jiang X, et al. Genome analysis of a major urban malaria vector mosquito Anopheles stephensi. Genome biol. 2014;15:459. https://doi.org/10.1186/s13059-014-0459-2.
Shinzawa N, et al. Phenotypic dissection of a plasmodium-refractory strain of malaria vector anopheles stephensi: the reduced susceptibility to P. berghei and P. yoelii. PLOS ONE. 2013;8:e63753.
Ngowo HS, et al. Fitness characteristics of the malaria vector Anopheles funestus during an attempted laboratory colonization. Malar J. 2021;20:148.
Chakraborty M, et al. Author Correction: hidden genomic features of an invasive malaria vector, Anopheles stephensi, revealed by a chromosome-level genome assembly. BMC Biol. 2022;20:96.
Thakare A, et al. The genome trilogy of Anopheles stephensi, an urban malaria vector, reveals structure of a locus associated with adaptation to environmental heterogeneity. Sci Rep. 2022;12:3610.
Bloomquist JR, Miller TA. Sodium channel neurotoxins as probes of the knockdown resistance mechanism. Neurotoxicology. 1986;7:217–23.
Enayati AA, et al. Molecular evidence for a kdr-like pyrethroid resistance mechanism in the malaria vector mosquito Anopheles stephensi. Med Vet Entomol. 2003;17:138–44.
Thompson M, Steichen JC, Ffrench-Constant RH. Conservation of cyclodiene insecticide resistance-associated mutations in insects. Insect Mol Biol. 1993;2:149–54.
Weetman D, Djogbenou LS, Lucas E. Copy number variation (CNV) and insecticide resistance in mosquitoes: evolving knowledge or an evolving problem? Curr Opin Insect Sci. 2018;27:82–8.
Shetty NJ, et al. Chromosomal inversions among insecticide-resistant strains of Anopheles stephensi Liston, a malaria mosquito. Parasitol Res. 2013;112:3851–7.
Battlay P, et al. Structural variants and selective sweep foci contribute to insecticide resistance in the drosophila genetic reference panel. G3 (Bethesda). 2018;8:3489–97.
Wood O, et al. Cuticle thickening associated with pyrethroid resistance in the major malaria vector Anopheles funestus. Parasit Vectors. 2010;3:67.
Chavshin AR, et al. Molecular characterization, biological forms and sporozoite rate of Anopheles stephensi in southern Iran. Asian Pac J Trop Biomed. 2014;4:47–51.
Alam MT, et al. The type and mysorensis forms of the Anopheles stephensi (Diptera: Culicidae) in India exhibit identical ribosomal DNA ITS2 and domain-3 sequences. Parasitol Res. 2008;103:75–80.
Firooziyan S, Dinparast Djadid N, Gholizadeh S. Speculation on the possibility for introducing Anopheles stephensi as a species complex: preliminary evidence based on odorant binding protein 1 intron I sequence. Malar J. 2018;17:366.
Dube VM, Gakhar SKJIS. Genetic differentiation between three ecological variants (‘type ’, ‘ mysorensis ’ and ‘ intermediate ’) of malaria vector Anopheles stephensi (Diptera: Culicidae). Insect Sci. 2010;17:335–43.
Singh OP, Mishra S, Sharma G, Sindhania A, Kaur T, Sreehari U, et al. Evaluation of intron-1 of odorant-binding protein-1 of Anopheles stephensi as a marker for the identification of biological forms or putative sibling species. PLoS ONE. 2022;17:e0270760. https://doi.org/10.1371/journal.pone.0270760.
Sweet WC, Rao BA. Races of A. stephensi liston, 1901. Indian Med Gazette. 1937;72:665–74.
Anyanwu GI, et al. Variation in cuticular hydrocarbons among strains of Anopheles (Cellia) stephensi liston possibly related to prior insecticide exposure. Ann Trop Med Parasitol. 1997;91:649–59.
Srinivasan S, Ghosh C, Das S, Thakare A, Singh S, Ganesh A, et al. Identification of a TNF-TNFR-like system in malaria vectors (Anopheles stephensi) likely to influence Plasmodium resistance. Sci Reports. 2022;12:19079. https://doi.org/10.1038/s41598-022-23780-y
World Health Organization. Test Procedures for Insecticide Resistance Monitoring in Malaria Vector Mosquitoes, 2nd ed.; World Health Organization: Geneva, Switzerland, 2016b.
The authors thank Dr. Rakesh Mishra, Director of the Tata Institute for Genetics and Society (TIGS), India, for his support. We are thankful to Dr. Susanta Kumar Ghosh, Senior Scientist (retired) of NIMR-ICMR, Bangalore, for guiding the establishment of TIGS insectary. We are grateful to Dr. Kailash Patra, R&D Director at Intrivo Diagnostics, California, USA, for giving training on parasite cultures and Mr. Joydeep Roy for his technical support. Animal work in the NCBS/inStem Animal Care and Resource Center was partially supported by the National Mouse Research Resource (NaMoR) grant #BT/PR5981/MED/31/181/2012;2013-2016 and 102/IFD/SAN/5003/2017-2018 from the Department of Biotechnology, GoI. The authors thank Tata Institute for Genetics and Society (TIGS) India for funding researchers involved in this work and IBAB for bioinformatics work.
This research is funded by internal grants from Tata Institute for Genetics and Society in Bangalore, India.
Ethics approval and consent to participate
Ethics approval ref. no. inStem/IEC-12/001, dated 20/03/2019. Standard membrane feeding system (approval ref. no. TIGS 2nd IBSC Oct 2018).
Consent for publication
All authors have provided consent for publication.
The authors have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Ghosh, C., Kumar, N., Kushwah, R.B.S. et al. Enrichment of phenotype among biological forms of Anopheles stephensi Liston through establishment of isofemale lines. Parasites Vectors 16, 79 (2023). https://doi.org/10.1186/s13071-023-05696-2