Chikungunya virus (CHIKV), an arthritogenic alphavirus, is transmitted to humans by infected Aedes (Ae.) aegypti and Ae.albopictus mosquitoes. In the study, reverse-transcription PCR (RT PCR) and virus isolation detected CHIKV in patient samples and also in adult Ae.albopictus mosquitoes that was derived from larvae collected during a chikungunya (CHIK) outbreak in Kerala in 2009. The CHIKV strains involved in the outbreak were the East, Central and South African (ECSA) genotype that had the E1 A226V mutation. The viral strains from the mosquitoes and CHIK patients from the same area showed a close relationship based on phylogenetic analysis. Genetic characterization by partial sequencing of non-structural protein 2 (nsP2; 378 bp), envelope E1 (505 bp) and E2 (428 bp) identified one critical mutation in the E2 protein coding region of these CHIKV strains. This novel, non-conservative mutation, L210Q, consistently present in both human and mosquito-derived samples studied, was within the region of the E2 protein (amino acids E2 200-220) that determines mosquito cell infectivity in many alpha viruses. Our results show the involvement of Ae. albopictus in this outbreak in Kerala and appearance of CHIKV with novel genetic changes. Detection of virus in adult mosquitoes, emerged in the laboratory from larvae, also points to the possibility of transovarial transmission (TOT) of mutant CHIKV strains in mosquitoes.
Chikungunya virus (CHIKV) is an alphavirus of the Togaviridae family and is an important re-emerging pathogen. It has been responsible for major fever epidemics in many parts of the world [1, 2]. The disease, chikungunya (CHIK), is characterized by high fever, headache, myalgia, severe and prolonged arthralgia, and erythematous skin rashes . In general, it is considered as a self-limiting illness. However, recent outbreaks of CHIK exhibited unusual severity, neurological complications and suspected mortality [3–6]. The disease is transmitted by the bite of Aedes ( Ae.) aegypti and Ae. albopictus mosquitoes. Studies have shown that Ae. albopictus facilitates rapid transmission of the new strains of CHIKV that had adaptive mutations in the viral genome [7, 8].
CHIK epidemic has caused considerable morbidity in recent years in India [9, 10]. Kerala, in South India, was one among the worst affected states [11–14]. Abundance of Ae.albopictus in many parts of the state was implicated for the rapid spread of the infection . Recent studies carried out in CHIKV from Kerala [11, 12, 14] have revealed novel genetic changes in the virus isolates from 2006-2008 outbreaks. Reports on virus isolation from mosquito vectors from the region are currently not available. The aim of the present work was to look for novel genetic changes in the isolates from 2009 by sequence analysis of selected genomic regions, and also to look for CHIKV in Ae. albopictus mosquitoes
The study was done during a fever outbreak in May-September 2009 in Kozhikkode district of northern Kerala (Figure 1). All the patients included in the study had classical symptoms of CHIK . Samples were obtained from the outpatient department of three Primary Health Centres (Olavanna, Beypore and Chaliyum) in the district. 2-5ml of whole blood was collected from patients who were clinically diagnosed with CHIK and had a history of fever of 1-5 days duration. Samples were transported to the laboratory in wet-ice; serum was separated and stored in aliquots at -80°C. Standard ethical and bio-safety guidelines were followed, and informed consent was obtained from all the patients prior to blood withdrawal.
For virus detection in mosquitoes, households of CHIK patients, whose serum samples were confirmed in the laboratory by RT-PCR, were subsequently visited and larval sampling was done. Stagnant water collected in discarded articles such as coconut shells, broken earthern-wares, plastic bottles and damaged drains were searched for Ae. albopictus larvae. Third and fourth instar larvae and pupae were phenotypically identified in situ using standard keys and these were collected and transferred to containers with fresh water. Four households each in Olavanna and Chaliyum, and three households in Beypore were surveyed. Larvae and pupae collected from each location were made into a single pool. In the laboratory, these three pools were independently reared in bowls with water, kept in mosquito cages at an ambient temperature of 25-30°C and a relative humidity 60-70%. The newly emerged adult mosquitoes were collected and frozen at -20°C for 30 minutes. Whole-mosquito tissue extracts were prepared by homogenizing pools of adult mosquitoes [each pool with 30 individual mosquitoes (both males and females) representing a single location]. Frozen mosquitoes were homogenized in 700 μl of Dulbecos Modified Eagle's Medium (DMEM) using a micropestle. These were then clarified by centrifugation at 800 × g at 4°C and sterilized by filtering through 0.2 μM membrane filter (Millex GV, Millipore) and used for RNA isolation.
RNA isolation from the 70 patient serum samples and the three extracts from mosquito samples were carried out using QIAamp Viral RNA Mini kit (Qiagen, GmBH, Hilden) exactly as per the kit protocol. Single-step RT PCR was done using 10 μl of the isolated RNA from all the samples using Fidelitaq RT-PCR kit (USB, Cleveland, Ohio), as previously described . PCR primers (Table 1; Figure 2) for CHIKV detection PCR were designed based on earlier reports  and on the conserved genomic regions of local strains of CHIKV . The conditions for RT PCR were: a reverse transcription step at 50°C for 45 min; followed by 35 cycles of thermal cycling, which included denaturation at 95°C for 1 min, annealing at 55°C for 1 min, and an extension at 68°C for 2 min. Extreme care was taken to avoid PCR-contamination, by carrying out the pre-and post amplification steps in laboratories located in separate buildings and also by including a non-template control in all amplifications.
Details of the primers used for PCR amplification in the study.
Sequence (5'→3'); location with respect to S27 sequence (GenBank Accession AF369024)
RT PCR for CHIKV detection in patient and adult mosquitoes derived from larvae
For nucleotide sequencing and phylogenetic analysis, 15 clinical samples and all the three mosquito-derived samples were used. Only clinical samples that gave a high intensity amplicon in the primary detection PCR were selected to ensure that sufficient DNA would be available for sequencing reactions. Five clinical samples each from Olavanna, Beypore and Chaliyum were used, making a total of 15 samples. Selected regions of the CHIKV genome (nucleotide position, with respect to S27 reference sequence AF369024: nsP2 3134-3636; E2 8832-9332; E1 10246-10539; Table 1; Figure 2) were amplified by RT PCR as described above using new sets of primers (Table1; Figure 2). These specific regions were chosen as they showed nucleotide variability and novel mutations in our previous study with the local strains of CHIKV , making them suitable for phylogenetic analysis. Purified PCR products were directly subjected to automated DNA sequencing as per manufacturer's directions in an ABI-Prism 3730 Genetic analyzer (PE Applied Biosystems, Foster City, CA). The sequences were aligned with corresponding CHIKV sequences obtained from NCBI GenBank using Clustal W program of MEGA3.1  software, with Kimura-2 distance correction. To get representation from different gene segments in the evolution of the CHIKV strains, the partial sequences of nsP2, E2 and E1 genes were arranged in tandem to obtain a 1311 bp sequence (Figure 3a), which was then used for phylogenetic analysis. The phylogeny was reconstructed by Neighbor-Joining method with 10,000 bootstrap replications using the MEGA 3.1 program. 100 μL of the mosquito extracts or patient serum samples were used for CHIKV isolation in confluent monolayer of Vero cells cultured in 75cm2 flasks, as per standardized protocols . The titration of CHIKV in the infected cultures was done by plaque assay using a carboxymethyl-cellulose overlay method  on Vero cells.
CHIKV RNA was detected in 49 out of the 70 patient samples (70%) and in adult mosquitoes derived from larvae from Chaliyum and Olavanna by RT PCR (Figure 4). All the three mosquito derived samples were positive for CHIKV, as indicated by cytopathic effects and RT PCR (Figure 4), in the 3rd passage of virus isolation in Vero cell monolayer cultures. In plaque assays, the culture supernatants from these infected cells had a virus titre of 2.0 × 1011, 3.3 × 1010, 1.4 × 1010 plaque forming units (pfu) ml-1 for samples from Olavanna, Chaliyum and Beypore, respectively.
Analysis of the partial nucleotide sequences of nsP2 (378 bp; position 3246-3623), E1 (position 10427-10931) and E2 (position 8893- 9320) revealed a few random nucleotide changes in the CHKV isolates studied (Additional File 1) with respect to the corresponding sequences of the previous isolates from Kerala . The nucleotide change T3297C observed in the 2007 & 2008 Kerala isolates, causing an L539 S mutation in the nsP2 protein, was absent in CHIKV strains of the present outbreak. A novel substitution (T3296C) was consistently observed in a few strains from patients (RGCB711, RGCB730, and RGCB755) and in all the three isolates from mosquito samples. However, this was a synonymous substitution. The E1 sequence of all the strains had the C10670T substitution resulting in the A226V mutation identified in the recent isolates of CHIKV [3, 14]. Another new substitution (E1 G10864A) detected consistently in all the mosquito-derived strains and two of the clinical isolates (RGCB711 & RGCB755) can result in an amino acid change of V291I. Two nucleotide substitutions (A9114G and T9170A) were observed in the E2 coding region of all the strains studied from the outbreak. The latter substitution resulted in an amino acid change L210Q in the predicted sequence of amino acids of the E2 protein. Phylogenetic analysis revealed that the strains involved in the outbreak were closely related to the East-Central South African genotype of the CHIKV (Figure 3b). The gene sequences of CHIKV obtained from mosquito and patient samples formed a close cluster, distinct from the strains isolated previously from Kerala , rest of India and other parts of the world. This show a common genetic origin of the virus strains from patients and mosquitoes in this outbreak.
Apart from these genetic changes, an interesting observation in the study was the detection of CHIKV from adult mosquitoes derived from larval samples. Considering that these mosquitoes were freshly emerged in the laboratory from the larvae collected from areas encountering a CHIK outbreak and did not have a blood-meal, the possibility of acquiring the virus through transovarial transmission (TOT) can be thought of. Even though TOT has been proven in flaviviruses [19–23], the occurrence of this phenomenon in alphaviruses is still inconclusive [24–27]. Studies using a Réunion Island isolate of the CHIKV from 2006 outbreak [Strain 06.21; GenBank: AM258992] could not demonstrate vertical transmission in the mosquito vector . The mosquito infectivity of alphaviruses is modulated by mutations in specific viral proteins [28–31]. Amino acid residues 200-220 of the E2 protein determine the cellular receptor tropism and mid-gut infectivity in Ae. aegypti mosquitoes [28, 30]. An E2 I211T mutation was found to strongly enhance Ae.albopictus infectivity of CHIKV strains with the E1 A226V change . Both the mutations were present in the isolates in this study and also in the recent Indian isolates [10, 12, 14] (Figure 5). Interestingly, the novel mutation in E2 (L210Q) that was detected exclusively in these 2009 CHIKV strains was adjacent to the E2-211 position. This substitution of the aliphatic amino acid leucine with glutamine, an amino acid with polar side chains, can have critical effects on local protein structure. One of the predicted effects of such amino acid changes is the exposure of buried protein surfaces. Possibly, this may alter the interaction of E2 with other proteins, particularly with cellular receptors, and may change the tissue tropism. However, more studies are required to understand the effects of the L210Q mutation.
The results from this study, along with the previous observations [11, 12, 14], indicate a constant genomic evolution of the CHIKV strains circulating in Kerala. The availability of large numbers of Ae.albopictus vector mosquitoes  and an immunologically naïve human population unexposed to CHIK in different parts of the state might facilitate recurrent infections and viral evolution. Emergence of newer strains with altered virulence and transmission potential is a possible out come of the long term viral persistence in the community. Further entomological and virological studies with these new CHIKV strains would help to understand the changing epidemiology of this re-emerging virus.
The authors are thankful to the medical staff of the primary health centres (PHC) in Olavanna, Beypore and Chaliyum for the help extended for the patient sample collection. The financial assistance by Department of Biotechnology, Government of India as intramural funding and the encouragement and support by the Director, RGCB, are gratefully acknowledged.
Molecular Virology Laboratory, Rajiv Gandhi Centre for Biotechnology (RGCB)
State Disease Control and Monitoring Cell (SDCMC), National Rural Health Mission (NRHM), Government of Kerala
Department of Community Medicine, Medical College
Powers AM, Logue CH: Changing patterns of chikungunya virus: re-emergence of a zoonotic arbovirus.J Gen Virol 2007, 88:2363–2377.PubMedView Article
Staples JE, Breiman RF, Powers AM: Chikungunya fever: an epidemiological review of a re-emerging infectious disease.Clin Infect Dis 2009, 49:942–948.PubMedView Article
Schuffenecker I, Iteman I, Michault A, Murri S, Frangeul L, Vaney MC, Lavenir R, Pardigon N, Reynes JM, Pettinelli F, et al.: Genome microevolution of chikungunya viruses causing the Indian Ocean outbreak.PLoS Med 2006, 3:e263.PubMedView Article
Rezza G, Nicoletti L, Angelini R, Romi R, Finarelli AC, Panning M, Cordioli P, Fortuna C, Boros S, Magurano F, et al.: Infection with chikungunya virus in Italy: an outbreak in a temperate region.Lancet 2007, 370:1840–1846.PubMedView Article
Chandak NH, Kashyap RS, Kabra D, Karandikar P, Saha SS, Morey SH, Purohit HJ, Taori GM, Daginawala HF: Neurological complications of Chikungunya virus infection.Neurol India 2009, 57:177–180.PubMedView Article
Das T, Jaffar-Bandjee MC, Hoarau JJ, Krejbich Trotot P, Denizot M, Lee-Pat-Yuen G, Sahoo R, Guiraud P, Ramful D, Robin S, et al.: Chikungunya fever: CNS infection and pathologies of a re-emerging arbovirus.Prog Neurobiol 2009, 91:121–9.PubMedView Article
Tsetsarkin KA, Vanlandingham DL, McGee CE, Higgs S: A single mutation in chikungunya virus affects vector specificity and epidemic potential.PLoS Pathog 2007, 3:e201.PubMedView Article
Vazeille M, Moutailler S, Coudrier D, Rousseaux C, Khun H, Huerre M, Thiria J, Dehecq JS, Fontenille D, Schuffenecker I, et al.: Two Chikungunya isolates from the outbreak of La Reunion (Indian Ocean) exhibit different patterns of infection in the mosquito,Aedes albopictus.PLoS One 2007, 2:e1168.PubMedView Article
Arankalle VA, Shrivastava S, Cherian S, Gunjikar RS, Walimbe AM, Jadhav SM, Sudeep AB, Mishra AC: Genetic divergence of Chikungunya viruses in India (1963–2006) with special reference to the 2005–2006 explosive epidemic.J Gen Virol 2007, 88:1967–76.PubMedView Article
Kumar NP, Joseph R, Kamaraj T, Jambulingam P: A226V mutation in virus during the 2007 chikungunya outbreak in Kerala, India.J Gen Virol 2008, 89:1945–1948.PubMedView Article
Santhosh SR, Dash PK, Parida MM, Khan M, Tiwari M, Lakshmana Rao PV: Comparative full genome analysis revealed E1: A226V shift in 2007 Indian Chikungunya virus isolates.Virus Res 2008, 135:36–41.PubMedView Article
Kannan M, Rajendran R, Sunish IP, Balasubramaniam R, Arunachalam N, Paramsivan R, Tewari SC, Samuel PP, Tyagi BK: A study on chikungunya outbreak during 2007 in Kerala, south India.Indian J Med Res 2009, 129:311–315.PubMed
Sreekumar E, Issac A, Nair S, Hariharan R, Janki MB, Arathy DS, Regu R, Mathew T, Anoop M, Niyas KP, Pillai MR: Genetic characterization of 2006–2008 isolates of Chikungunya virus from Kerala, South India, by whole genome sequence analysis.Virus Genes 2010, 40:14–27.PubMedView Article
Hasebe F, Parquet MC, Pandey BD, Mathenge EG, Morita K, Balasubramaniam V, Saat Z, Yusop A, Sinniah M, Natkunam S, Igarashi A: Combined detection and genotyping of Chikungunya virus by a specific reverse transcription-polymerase chain reaction.J Med Virol 2002, 67:370–374.PubMedView Article
Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment.Brief Bioinform 2004, 5:150–163.PubMedView Article
Escobar-Herrera J, Medina-Ramirez FJ, Gutierrez-Escolano AL: A carboxymethyl-cellulose plaque assay for feline calicivirus.J Virol Methods 2007, 146:393–396.PubMedView Article
Rosen L, Shroyer DA, Tesh RB, Freier JE, Lien JC: Transovarial transmission of dengue viruses by mosquitoesAedes albopictusandAedes aegypti.Am J Trop Med Hyg 1983, 32:1108–1119.PubMed
Gunther J, Martinez-Munoz JP, Perez-Ishiwara DG, Salas-Benito J: Evidence of vertical transmission of dengue virus in two endemic localities in the state of Oaxaca, Mexico.Intervirology 2007, 50:347–352.PubMedView Article
Mishra AC, Mourya DT: Transovarial transmission of West Nile virus inCulex vishnuimosquito.Indian J Med Res 2001, 114:212–214.PubMed
Rosen L, Shroyer DA, Lien JC: Transovarial transmission of Japanese encephalitis virus byCulex tritaeniorhynchusmosquitoes.Am J Trop Med Hyg 1980, 29:711–712.PubMed
Aitken TH, Tesh RB, Beaty BJ, Rosen L: Transovarial transmission of yellow fever virus by mosquitoes (Aedes aegypti).Am J Trop Med Hyg 1979, 28:119–121.PubMed
Mourya DT: Absence of transovarial transmission of Chikungunya virus inAedes aegypti&Ae. albopictusmosquitoes.Indian J Med Res 1987, 85:593–595.PubMed
Vazeille M, Mousson L, Failloux AB: Failure to demonstrate experimental vertical transmission of the epidemic strain of Chikungunya virus inAedes albopictusfrom La Reunion Island, Indian Ocean.Mem Inst Oswaldo Cruz 2009, 104:632–635.PubMedView Article
Thavara U, Tawatsin A, Pengsakul T, Bhakdeenuan P, Chanama S, Anantapreecha S, Molito C, Chompoosri J, Thammapalo S, Sawanpanyalert P, Siriyasatien P: Outbreak of chikungunya fever in Thailand and virus detection in field population of vector mosquitoes,Aedes aegypti(L.) andAedes albopictusSkuse (Diptera: Culicidae).Southeast Asian J Trop Med Public Health 2009, 40:951–962.PubMed
Myles KM, Pierro DJ, Olson KE: Deletions in the putative cell receptor-binding domain of Sindbis virus strain MRE16 E2 glycoprotein reduce midgut infectivity inAedes aegypti.J Virol 2003, 77:8872–8881.PubMedView Article
Pierro DJ, Powers EL, Olson KE: Genetic determinants of Sindbis virus strain TR339 affecting midgut infection in the mosquitoAedes aegypti.J Gen Virol 2007, 88:1545–1554.PubMedView Article
Pierro DJ, Powers EL, Olson KE: Genetic determinants of Sindbis virus mosquito infection are associated with a highly conserved alphavirus and flavivirus envelope sequence.J Virol 2008, 82:2966–2974.PubMedView Article
Tsetsarkin KA, McGee CE, Volk SM, Vanlandingham DL, Weaver SC, Higgs S: Epistatic roles of E2 glycoprotein mutations in adaption of chikungunya virus toAedes albopictusandAe. aegyptimosquitoes.PLoS One 2009, 4:e6835.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.