Abstract

Understanding the population dynamics of mosquito populations through life table analysis and insecticide susceptibility is important to assess the likely impact of vector control strategies as well as to aid the design of novel interventions. Variation in the life tables and other biological data was compared for two populations of Culex quinquefasciatus Say 1823 from geographically isolated regions, Gorakhpur and Pune from India. Under a standardized rearing regime and constant laboratory conditions, mosquitoes were reared and biological attributes of these populations were compared. Development and survival of immature and adult stages of Culex quinquefasciatus were found significantly different in Gorakhpur and Pune populations. Principal component analysis of morphological data revealed that the two populations form significantly different clusters which can be differentiated from each other based on siphon, saddle, anal gills, and pecten teeth related variables. Insecticide susceptibility results suggest that the larvae from both areas were more susceptible to deltamethrin as compared to DDT and malathion. The current study provides baseline information on survivorship, morphological variation and insecticide susceptibility of Culex quinquefasciatus. The results obtained in this study suggest that different geographical areas with contrasting habitats have significant influence on survival and reproductive strategies of Culex quinquefasciatus.

1. Introduction

Feasible mosquito intervention strategies require in-depth knowledge of the population dynamics and insecticide susceptibility status. To better understand the population dynamics of the vector species, understanding the parameters such as development time, rate of survival, and fecundity of particular importance. A life table describes the development, survival, and fecundity of a cohort and provides basic data on population growth parameters. Life table study also provide, concise information on inherent differences in the survivorship and reproductive strategies of populations under different ecological regimes [1] and can help in explaining why certain species proliferate in particular environments.

Culex quinquefasciatus Say is a cosmopolitan mosquito species found in tropical, subtropical, and warm temperate regions. Cx. quinquefasciatus acts as a vector of filarial worms [2, 3], protozoan parasites [4], and various arboviruses [59]. In India, Cx. quinquefasciatus is the principal vector of bancroftial filariasis [10]. Rapid urbanization and industrialization without proper drainage facilities are responsible for the establishment of the Cx. quinquefasciatus in various cities of India [11]. Of the estimated global 128 million lymphatic filarial cases, 91 per cent are caused by Wuchereria bancrofti Cobbold [12]. India alone contributes 40 per cent to the world’s lymphatic filariasis disease burden [13]. The estimates of disease burden due to filariasis suggest that 2.06 million disability adjusted life years (DALYs) are lost in India [14].

Numerous studies have been carried out on Cx. quinquefasciatus biology and insecticide susceptibility [1528]. Recently, Suman et al. [25, 28] reported significant variations in morphology among eggs and life table attributes from the four strains of Cx. quinquefasciatus. Comparison of biological attributes such as insecticide susceptibility, life table analysis, and morphological variation of Cx. quinquefasciatus of endemic and nonendemic areas of lymphatic filariasis is not well documented.

Mosquitoes inhabiting in distinct location normally adapt to local conditions. Gorakhpur and Pune are geographically well separated (~1500 km) (Figure 1) and having different climatic conditions. Uttar Pradesh (Gorakhpur) is one of the lymphatic filariasis endemic states of the India whereas Maharashtra (Pune) is nonendemic to lymphatic filariasis. The environment and climatic features of these two ecotypes are also different. It seems possible that environmental stress of these two ecoclimatic regions influence their life history traits. We hypothesized that these mosquito strains will differ in their life table and biological attributes because they appear to be adapted to their respective local environments. To test this hypothesis, we investigated the life table analysis, morphological variation, and insecticide susceptibility between these two populations of Cx. quinquefasciatus from geographically well-separated areas in India, namely, Gorakhpur and Pune.

2. Materials and Methods

2.1. Ethics Statement

All animals were handled in strict accordance with good animal practice as defined by Institutional Animal Ethics Committee (IAEC) affiliated with National Institute of Virology (NIV), Pune, India. All animal work was approved by the IAEC. Animal experiments were carried out in strict compliance with Committee for the Purpose of Control and Supervision of Experiment on Animals (CPCSEA) guidelines, India.

2.2. Mosquito Collection and Rearing

Adult and larvae of Cx. quinquefasciatus were sampled from two populations to examine the possibility of population differences in the endemic and nonendemic areas of lymphatic filariasis infection. Uttar Pradesh is an endemic area of JE virus and lymphatic filariasis infection and Maharashtra is nonendemic to JE virus and lymphatic filariasis infections. Gorakhpur (26°45′ N, 83°24′ E) is situated in the Tarai belt near Himalaya ranges in Uttar Pradesh state bordering Nepal whereas Pune (18°31′ N, 73°55′ E) is situated near Western Ghats of India, in Maharashtra state. Gorakhpur is situated in the Indo-Gangetic Plain region that receives 1200–1300 mm of rainfall per year. Annual temperature ranges between 8°C in winter and 42°C in summer. The vegetation of this region is tropical, moist, and deciduous. Pune is situated in a semiarid region that receives 700–800 mm of rainfall per year. Annual temperature ranges between 12°C in winter and 38°C in winter. The vegetation of this region is dry deciduous. Larvae were collected from temporary and semipermanent ground water pools and fields. Morphometric analysis was carried out on these field collected mosquitoes. In the present study, to better understand baseline parameters of the life cycle of Culex quinquefasciatus, colonies of this species from two different ecological areas were monitored under laboratory conditions for life table analysis and insecticide susceptibility. Individuals of Cx. quinquefasciatus were identified using keys of Barraud [29]. Voucher specimens were deposited in National Institute of Virology museum, Pune, India.

2.3. Life Table Analysis

Rearing of mosquito in laboratory environment imposes a certain degree of selection pressure on mosquito biology. Mosquito colonies were maintained in standardized environmental conditions which were favorable to mosquito survival. Hence, we assume that life-table parameters derived from data collected from the colonized wild strains represent maximum estimates of their life-table parameters and are likely to reflect true differences between geographically isolated strains. Similar assumptions were made earlier by Reisen et al. [30] and Suman et al. [28] to compare the life-table attributes of geographically distinct strains of Culex mosquitoes.

Colonies of each population were maintained at National Institute of Virology Pune in standard laboratory conditions at 28 ± 1°C, 70 ± 5% relative humidity (RH) and light: dark (LD) 12 : 12 h. These colonies were maintained for up to seven generations (approximately 4 months). The adults were offered 10% glucose solution soaked in cotton pads daily, as a source of nutrition. A fowl was provided as a source of blood twice a week. Fowls were maintained in accordance with approved guidelines of the Committee for the Purpose of Control and Supervision of Experiments on Animals (CPCSEA) India. Egg rafts were collected from the individual mosquitoes in plastic cups (60 ml) containing water and lined with filter papers, in the cage. The single egg raft was transferred from the filter papers with the aid of mounting needles into plastic bowls (10 cm height and 25 cm diameter) where they were held for 24 to 72 h for hatching. Twenty-four hours after hatching, the larvae were maintained at a density of 100 larvae per tray (45 × 30 × 10 cm) containing 2 L of dechlorinated water. The larvae were fed with mixture of yeast extract and dog biscuit powder, sprinkled on the water surface. On every alternate day, the water from the culture tray was changed carefully until pupation. The pupae were separated from the larvae daily and placed in plastic bowls (10 cm height and 30 cm diameter) half-filled with water. These plastic bowls with pupae were placed in adult-holding cages for emergence. Adults after emergence were offered 10% glucose solution soaked in cotton pads daily. Fecundity was estimated based on the average numbers of eggs laid per female. Egg-hatching rate was estimated from the number of eggs hatched per female. Duration of the preadult development periods was determined by observing each pan daily, and all larval skins were removed, scored to instar, and counted. There were two replicates for each experiment described above and the whole study was repeated three times, resulting in the monitoring of about 600 larvae from first instar stage to adulthood.

2.4. Survival Analysis

One hundred newly emerged larvae were placed in basins. All basins were screened with white, insect-proof, nylon netting to prevent colonization by other mosquitoes and predators. Each day the number of larvae alive and their stage of development were counted and recorded. Adult mosquitoes that emerged were counted daily, and the sexes were recorded. Standard life table analysis was performed on the data (Table 1). Computation of stage-specific survivorship was according to the formula: where is the number of individuals entering stage and is the number of individuals that entered the preceding stage.

2.5. Morphological Analysis

Field collected larvae were used for morphometric analysis. Eight morphological characters and three ratios were scored for morphological analysis of fourth instar larvae (Table 3). Siphon index, anal gill index, and siphon saddle ratio are used in traditional taxonomy to distinguish among species and subspecies. However, Fakoorziba and Vijayan [31] and Kanojia et al. [32] have used these ratios even to distinguish between different populations. The characters of fourth-instar larvae were measured using micrometric oculars with the least count 0.01 mm.

2.6. Adult Mosquito Bioassays

Insecticide susceptibility tests were carried out using the standard WHO protocol, insecticide susceptibility test kits, and impregnated papers. Two- to three-day-old nonblood-fed adult female Cx. quinquefasciatus were tested. Three batches of 20 mosquitoes each were exposed to test papers impregnated with DDT (4%), Deltamethrin (0.05%), and Malathion (5%). Controls included batches of mosquitoes from each site exposed to untreated papers.

2.7. Larval Bioassays

Batches of fourth instars were exposed to insecticide solution (DDT, Deltamethrin and Malathion) diluted in distilled water. Twenty larvae per concentration and five concentrations per test, providing mortalities ranging from 0 to 100%, were run. Temperature was maintained at 28°C throughout the test, and larval mortality was recorded after 1 h exposure and then corrected for control mortality [33]. Three replicates with insect larvae from different rearing batches were made at different times, and the results were pooled. The data were subject to a log-probit analysis [34] to determine the LC50 and LC95 values as well as their 95% confidence intervals (CI).

2.8. Analysis of Data

Survival was evaluated by scoring the number of dead mosquitoes/larvae every day, to estimate mosquito stage specific lifespan. Larvae those who are entered in pupal stage were considered as censored data. The survival rate for each mosquito population was described using life table analysis. These survival curves were then compared using the logrank test. Pairwise comparison of each character was carried out using -test, with the Bonferroni adjustment to the probabilities (as 11 comparisons were made therefore we used as our cutoff value). Data were analyzed using multivariate statistics. Multivariate analysis offers the advantage of taking into account all the variables in a single analysis, thus making it possible to assess variation in the morphological characters of larvae from different collection sites. Principal component analysis (PCA) was done to create uncorrelated principal components from the original variables. Factor scores estimate the actual values of individual observations for the factors, whereas correlation between variables and factors is called as factor loading. Principal components were extracted from the covariance matrix. Mean values (centroids) and 95% asymptotic confidence limits of the scores of individual larvae on the first two principal components were computed for each sample. First two principal factors, which explained maximum variation in the data, were analyzed to understand the variation in the morphology of different individuals.

2.9. Analysis of Bioassay Data

Bioassay data were pooled and analyzed (LC50 and 95% LC values) by using SPSS ver. 16 program (SPSS Inc., Chicago, IL).

2.9.1. Abbott’s Formula

Bioassay data were not considered if the control mortalities were >20%. Actual mortalities were calculated against control mortalities by using Abbott’s formula [33]. Corrected mortality = (% survival in control experiment − % survival in treated experiment)/(% survival in control experiment) × 100.

2.9.2. Log-Probit Regression Lines

Log-probit regression lines were obtained by plotting percentage mortalities ( axis—probability) against concentration ( axis—log). LC50 (dosage which kills 50% of the population) and LC95 (dosage which kills 95% of the population) values were obtained. Fitness of the regression lines was statistically tested using values [35].

3. Results

3.1. Preadult Development and Survival

The mean duration and survival rates of different life stages of Cx. quinquefasciatus of Gorakhpur and Pune are depicted in Table 1. The percentage of egg hatching was lower in Gorakhpur (86.5 ± 7.3) as compared to Pune (91.5 ± 3.8). Larval-to-adult survival rate of Cx. quinquefasciatus in Pune (58.063.61) was significantly higher than Gorakhpur (28% ± 0.76) (-test ). Pune population of Cx. quinquefasciatus had slightly longer larval development time (14.57 days) when compared with that of Gorakhpur Cx. quinquefasciatus population (13.33 days). Larval survival rate of Pune population (68%) was higher than Gorakhpur population (42%) of Cx. quinquefasciatus (-test ). Pupal survival rates of both populations were not significantly different (-test ). No significant difference was observed in male female emergence ratio (-test ). Survival rate of Cx. quinquefasciatus adults of Pune was higher than those of Gorakhpur (logrank test; on 1 degrees of freedom, ) (Figure 2).

3.2. Fecundity

Comparison of egg production of Gorakhpur population and Pune population of Cx. quinquefasciatus is reported in Table 2. Gorakhpur population produced bigger egg raft (151–186 eggs) as compared to Pune population (121–160 eggs). The number of eggs per raft were higher in Gorakhpur population (169 ± 14.5) as compared to Pune population (139.75 ± 15.96) (Table 2). The number of eggs/female life span was higher in Gorakhpur population (676 ± 32) (Table 2). Eggs/female/day was ranged from 32.19 (1.68) in Gorakhpur population and 26.59 (2.45) in Pune population (Table 2).

3.3. Phenotypic Variation among the Populations

The pairwise comparison of morphological characters using unpaired -test revealed that characters like length of anal gills, length of siphon, denticles on the apical pecten teeth, siphon saddle ratio, and anal gill index were significantly different in both populations studied (Table 3). PCA of the morphological analysis extracted four factors with eigenvalues more than one. Cumulatively, these factors explained 79.61% of the total variability in the data. The first factor (F1) explained 31.97% of the total variability while the second factor (F2) explained 19.73% of total variability and together the first two factors explained 51.70% of the total variability. Characters which showed high factor loading on F1 included length of siphon, width of anal gill, siphon index, and siphon/saddle ratio, while characters which showed high factor loading on F2 included length of anal gill, length of saddle, and anal gill index. PCA could distinguish two significant clusters for Cx. quinquefasciatus of Pune, Cx. quinquefasciatus of Gorakhpur (Figure 3). Our null hypothesis that there is no significant difference in morphometry of Cx. quinquefasciatus populations was rejected.

3.4. Insecticide Susceptibility

Based on WHO criteria [36], Cx. quinquefasciatus from Gorakhpur and Pune showed complete susceptibility to deltamethrin (100% mortality), potential resistance to malathion (~80% mortality), and complete resistance to DDT (Less than 80% mortality). No significant difference was noted in the toxicity of diagnostic concentrations of DDT, deltamethrin, and malathion between the Gorakhpur and Pune populations of Cx. quinquefasciatus (-test ). The larvicidal potency of DDT, deltramethrin, and malathion against Cx. quinquefasciatus is depicted in Table 4. Gorakhpur population showed higher LC50 value for DDT and malathion as compared to Pune population. No significant difference was observed for LC50 values of deltamethrin for both the populations.

4. Discussion

The effects of geographical and environmental events on life table attributes and morphological structure of local populations of Cx. quinquefasciatus were analyzed in this study. Our null hypothesis that there is no significant difference between the life cycle parameters, morphology, and profiles of Gorakhpur and Pune populations was rejected. Gorakhpur and Pune populations of Cx. quinquefasciatus are significantly different with respect to life table attributes and morphological characters.

Gorakhpur and Pune are ~1500 km apart, Gorakhpur area having the higher range of annual temperatures, lower relative humidity for most of the year, and higher annual rain fall as compared to Pune. Accordingly, we might expect that environmental conditions in Pune would favour survival and reproductive potential of mosquitoes. The population from Gorakhpur area exposed to harsh conditions for survival in nature than Pune strain had a shorter life expectancy (20–23 days in Gorakhpur vs 24–27 days in Pune), laid eggs on most days (18–20 days in Gorakhpur vs 20–22 days in Pune), and constructed larger egg rafts (169 ± 14.5 eggs/raft in Gorakhpur vs 139.75 ± 15.96 eggs/raft in Pune), indicating that Gorakhpur strain invested the most energy in rapid reproduction. The higher survivorship of adult Cx. quinquefasciatus in the Pune compared to the Gorakhpur is likely caused because of favorable habitat and environmental conditions in Pune.

The percentage of eggs hatching was not significantly different in Gorakhpur (86.5 ± 7.3) population as compared to Pune (91.5 ± 3.8) (-test ). The percentage of eggs hatching has been reported to be 82.1% in Culex pipiens fatigans (synonymous with Cx. quinquefasciatus; Gomez et al. [37]), 79.6% in Cx. quinquefasciatus (Suleman and Reisen [38]), and 80.5%–95.6% in Cx. quinquefasciatus (Suman et al. [28]; five populations from India) which are close to the values reported in the present study (Table 1). The mean development time from egg to adult was found to vary little, with showing a range of ≤1 day (15.4 days in Gorakhpur and 16.5 days in Pune). These values were higher than the developmental period reported for Cx. quinquefasciatus collected from Bikaner, Jamnagar, and Bhatinda areas of India [28]. It has been documented that variations in immature mosquito developmental rates depend on habitat conditions [28, 3942]. The results obtained in this study also suggest that habitat and environmental conditions play important role in biological attributes of mosquitoes.

Life expectancy is an important attribute of vector capacity as it plays important role in transmission pathogens. Studies on lymphatic filariasis transmission indicate that microfilariae require 16–17 days in the mosquito vector to reach the infective stage [43]. This period shortens with higher temperature and humidity [44] but increases up to 42 days at lower temperatures [45]. As life table analysis suggests that Gorakhpur population invest more energy in reproduction it is likely to take blood meal more frequently. The repeated biting may cause spread of diseases rapidly. Higher ambient temperatures and frequent biting nature of Cx. quinquefasciatus of Gorakhpur might be responsible for frequent cases of lymphatic filariasis in this area.

Multivariate analysis of morphological characters revealed that two populations of Cx. quinquefasciatus differ significantly in morphological characters (Table 2; Figure 3(a)). Principal component analysis of the data suggested that siphon, saddle, anal gills and pecten teeth related variables were the most important distinguishing characters (Figure 3(b)). Anal gills are longer in Gorakhpur population as compared to Pune population (Table 3). As salt concentration of water increases the length of anal gills starts reducing [46]. The salt concentration of Pune water might be a probable reason behind the reduced anal gill length. The causes of morphological differences between populations are often quite difficult to explain. Phenotype is under the double control of environmental conditions and genotype, but morphological changes can be rapid when different environmental conditions occur [47]. We hypothesize that transition of saddle and anal gill length in Cx. quinquefasciatus populations is a consequence of environmental conditions selection.

The history of insecticide resistance in Culex species suggests that members of Culex genus have a notorious reputation for developing resistance to insecticides, including organophosphates, carbamates, and pyrethroids [17, 21, 22, 48, 49]. In India, first evidence of development of DDT resistance in Cx. quinquefasciatus was observed in 1952 from a village near Delhi [50]. Later, several areas like Nagpur, Pune, Patna, and Rajahmundry reported resistance in Culex mosquitoes to various insecticides such as BHC, fenitrothion, DDT, malathion, and temephos [15, 16, 18, 20, 23, 51]. Investigation by Mukhopadhyay et al. [23] showed no mortality in Cx. quinquefasciatus adults against the 4% DDT and 5% malathion. About 8% and 14% mortality were found in Cx. quinquefasciatus larvae in 0.125 and 3.125 mg/l concentration of DDT and malathion, respectively. The present study concludes that the larvae of both the populations of Cx. quinquefasciatus are highly susceptible to deltamethrin. However, Cx. quinquefasciatus populations in these areas showed resistance to DDT and malathion. Therefore, as a baseline study it may be concluded that in future prudent use of vector control strategies should be applied to reduce the Cx. quinquefasciatus nuisance in Gorakhpur and Pune.

In summary, this study provides baseline information on survivorship, morphological variation, and insecticide susceptibility of Culex quinquefasciatus populations. The climate conditions of Pune seems more permissive to the development and survival of Culex quinquefasciatus. The insecticide susceptibility results revealed that in the future different control measures will be required to reduce the Cx. quinquefasciatus populations. This knowledge, coupled with enhanced understanding of the ecology of Cx. quinquefasciatus, will enable locally specific, efficient deployment of public health resources in reducing the burden of Cx. quinquefasciatus.

Conflict of Interests

The authors have declared that no competing interests exist.

Funding

The study was supported by the Indian Council of Medical Research, Government of India. The funding agency had no role in study design, data collection and analysis, decision to publish, or preparation of the paper.

Acknowledgments

The authors thank Dr. A. C. Mishra, Director, National Institute of Virology, Pune, for the facilities and encouragement. They appreciate the technical assistance provided by the field staff of NIV. They are thankful to Dr. Neelesh Dahanukar, IISER, Pune, for the valuable suggestions and help in statistical analysis.