- Open Access
Evidence of co-infection of chikungunya and densonucleosis viruses in C6/36 cell lines and laboratory infected Aedes aegypti (L.) mosquitoes
© Sivaram et al; licensee BioMed Central Ltd. 2010
- Received: 13 August 2010
- Accepted: 12 October 2010
- Published: 12 October 2010
Densonucleosis viruses are the etiological agents of insect's disease. We have reported the isolation of densovirus from India and its distribution among the natural populations of Aedes aegypti mosquitoes across the country. Since densonucleosis virus persistently infects mosquito populations, and is demonstrated to negatively affect multiplication of dengue virus in Aedes albopictus, it would be interesting to study if this virus has a role in determining the susceptibility of the vector mosquito Ae. aegypti to chikugunya virus.
Mosquito cell lines and adult Ae. aegypti mosquitoes infected with densovirus were superinfected with Chikungunya virus and both the viruses were quantitated by determining their genomic copy number by real time amplification. Comparison was made between the log of genomic copy numbers of the viruses in the presence and absence of each other.
The log of copy number of the viruses did not vary due to co-infection. Even though the RNA copy number of chikungunya virus increased over the period of time, no change was observed in the RNA copy number between the control and the co-infected group on any given day. Similarly, DNA copy number of densovirus also remained unchanged between the control and the co-infected groups.
Chikungunya virus neither stimulates the replication of densovirus nor is its own replication suppressed due to co-infection. Ae. aegypti mosquitoes with densovirus infection were as susceptible to infection by chikungunya virus as the uninfected mosquitoes.
- Post Infection
- Japanese Encephalitis
- Genomic Copy Number
- Chikungunya Virus
Densonucleosis viruses (DNV) are known to persistently infect the mosquito cell lines and mosquito populations in nature. They belong to the family Parvoviridae and genus Brevidensovirus. DNV have been isolated from several species of mosquitoes like Culex salinarius, Aedes aegypti[2, 3], Ae. albopictus, Anopheles minimus and mosquito cell lines . They are isometric, non-enveloped viruses of about 20 nm in diameter with single stranded linear DNA, packaged either as plus or a minus strand . Palindromic sequences that form hair pin structure are found in both the ends of the genome [8, 9]. The coding sequences include the genes for non structural proteins NS1 and NS2 occupying the 5' end of the strand and the genes for structural protein VP occupying the 3' end [3, 9, 10].
DNV is pathogenic to mosquitoes, resulting in loss of mobility, deformity and loss of pigmentation in the mosquito larvae. A very high rate of mortality (up to 90%) was reported when first instar Ae. aegypti larvae were exposed to Ae. albopictus DNV . The pathogenicity of the virus is dose dependent, where in at lower doses of infection with DNV, the larvae survive to form adults that carry the infection [4, 12]. Virus multiplication has also been observed in the ovaries of the infected mosquitoes, which enable them to transmit the virus vertically to the next generation [11, 13]. DNV has a role in significantly reducing the vectorial capacity of the mosquitoes by reducing mosquito lifespan . DNV has also been developed as a transducing vector used to introduce constructs that interfere with the arboviral infection in mosquitoes.
Vector competence refers to the intrinsic permissiveness of an arthropod vector for infection, replication and transmission of a vertebrate pathogen . Vector competence is governed by intrinsic (genetic) factors and extrinsic factors [15, 16]. The microflora present in the midgut of mosquito might also have a role in determining their susceptibility to the viruses . Since DNV has been seen to be persistently infecting mosquito populations, it would be interesting to study if this virus has a role in determining the susceptibility of the vector to other arboviruses. C6/36 mosquito cell cultures persistently infected with AalDNV showed markedly lower CPE than naive-cell cultures or acutely AalDNV infected cultures when super-challenged with DENV-2 . Paterson et al reported that persistent infection with the low virulence type of DNV did not prevent severe CPE due to super-challenge with a more virulent type. DNV infected Ae. albopictus mosquitoes have been reported to show reduced numbers of DENV-2 when super challenged as compared to those uninfected with DNV . Their study also showed that super-challenge with DENV-2 stimulated DNV replication. Recent studies show that mosquito cells can accommodate balanced, persistent co-infections with a DNV and DENV . When cells dually infected with DENV and DNV were super-challenged with Japanese encephalitis (JE), cultures were stable without signs of cytopathology, with 99% cells producing antigens of all the three viruses .
Chikungunya fever (CHIK) is a viral disease transmitted by Aedes mosquitoes. CHIKV belongs to the family Togaviridae, genus Alphavirus. The disease has a significant potential to spread globally given the wide distribution of its arthropod vector [23, 24]. In India CHIKV re-emerged after nearly 32 years in October 2005 in a very high magnitude. During this outbreak of CHIKV, our laboratory received mosquitoes from which we had isolated and characterized Ae. aegypti DNV and studied its distribution among different Ae. aegypti populations across India . Whether the susceptibility of Ae. aegypti mosquitoes to CHIKV is affected due to coinfection of CHIKV and DNV is unknown. The present communication focuses on the effect of co-infection of DNV and CHIKV in the mosquito cell line and in Ae. aegypti mosquitoes.
Approximately 10,000 fourth instar larvae from the DNV infected colony of Ae. aegypti maintained at National Institute of Virology (NIV), Pune were triturated in 150 ml of Phosphate Buffer Saline (PBS) using a tissue homogenizer (Fisher Scientific). The resulting slurry was centrifuged at 10,000 rpm for 30 min at 4°C to remove the cell debris. The supernatant was aliquoted and stored at -80°C until used. CHIKV isolate (Strain No: 601573) was obtained from virus repository of this institute and stock was prepared in Vero E6 cell lines. These viral stocks were used throughout the study. Both the virus stocks were quantitated by real time PCR.
C6/36 cells were grown in Mitsuhashi Maramorosch (MM) media , supplemented with 10% fetal bovine serum (FBS) and 1% penicillin - streptomycin. The cell line was maintained at 28°C.
The Ae. aegypti mosquitoes used for these studies were from an insectary maintained at NIV, Pune for the last 35 years. The mosquitoes were maintained at 28 ± 2°C and a relative humidity of 70-80%. The larvae were fed on sterilized yeast tablets and dog biscuits mixed in the proportion of 70 and 30%, respectively. The adults were reared in cages and were fed on 10% glucose. Female mosquitoes were fed with chicken blood every third day to obtain eggs.
Primer designing for SYBR GREEN real time PCR
Primers designed for real time PCR for DNV and CHIKV
Name of Primer
Denso SYBR Fwd
Denso SYBR Rev
CHIK Q Fwd
CHIK Q Rev
SYBR Green I quantitative real time PCR of DNV
DNA extracted from DNV stock using QiaAmp DNA Mini kit (Qiagen, Hilden, Germany) was amplified with DNV specific primers . The amplicon of 1.1 kb was cloned into E. coli cells using TA cloning kit (Invitrogen, Carlsbad, CA) following manufacturer's protocol. Blue white colony screening was used to select the transformed cells. Plasmids were extracted from positive clones using Plasmid Mini kit (Qiagen, Hilden, Germany) and the insert was sequenced on an ABI 3100 automated DNA sequencer using Big Dye terminator kit (Applied Biosystems). The concentration of the plasmid was determined using NanoDrop spectrophotometer (NanoDrop technologies, Inc., USA). The plasmid was serially diluted 10 fold and used as standard for real time PCR. The reaction mixture for real time amplification composed of 12.5 μl of Power SYBR Green PCR Master Mix (Applied Biosystems, Foster City, CA), 1 μl of each primer (10 pico moles) and 5 μl of template DNA. The mixture was made up to 25 μl by adding 5.5 μl of water. The real time PCR conditions for DNV were 10 min of initial hot start at 95°C, denaturation at 95°C for 15 sec and 1 min of annealing and extension at 60°C for 40 cycles. The reactions were carried out in ABI 7300 real time PCR system (Applied Bio-systems). The system software provided (SDS) was used to analyze the results by plotting standard curve.
SYBR Green I quantitative real time RT-PCR of CHIKV
RNA was extracted from CHIKV stock using QIAamp viral RNA mini kit (Qiagen, Hilden, Germany) and was reverse transcribed to cDNA using Reverse transcription system (Promega Corp., Madison, USA) according to the manufacturer's protocol. The reverse primer designed for real time amplification of CHIKV was used to synthesize cDNA. The forward primer used to amplify the cDNA contained a site of T7 polymerase which could be anchored to the amplicon to enable in vitro transcription. The amplicon of 250 bp was subjected to in vitro transcription using T7 RNA polymerase by Riboprobe in vitro transcription system (Promega Corp., Madison, USA) according to manufacturer's instructions. The transcript was quantitated using NanoDrop spectrophotometer. The RNA was reverse transcribed to cDNA and serially diluted 10 fold, which served as standard. The reaction mix for real time amplification was prepared similar to that of DNV. The real time amplification conditions were 10 min of initial hot start at 95°C, denaturation at 95°C for 15 sec and 30 secs of annealing at 50°C and 30 secs of extension at 62°C for 40 cycles. Results were analyzed by plotting standard curve using SDS software.
Determining the effect of co-infection of CHIKV and DNV in C6/36 cell line
C6/36 cells were counted using hemocytometer and were seeded in 6 well plates at 3000 cells per well. When the cells were 80% confluent, the media was removed and the cells were infected by adding DNV to the cell sheet. The quantity of DNV used to infect corresponded to 6.07 × 108 viral particles. The cells were incubated at 28°C for one hour with intermittent rocking of the plate at every 15 min to enable adsorption of the virus. After incubation, the inoculum was removed, and the cells were washed with sterile PBS, fed with 2 ml of fresh media supplemented with 2% FBS and returned to the incubator.
The cells infected with DNV were incubated for two days. On day 3 post DNV infection, the cells were superinfected with 5.73 × 105 RNA copy number of CHIKV. All infections and superinfections were done in triplicates. Appropriate controls were used which include cells infected with only CHIKV, only DNV and uninfected cells.
Virus was harvested from the infected cells every 24 hr post superinfection for 4 consecutive days by repeated freeze thaw cycles and centrifugation at 5000 rpm for 20 mins at 4°C. The supernatant was divided into two batches, one of which was used to quantitate DNV and the other to quantitate CHIKV. Nucleic acids extracted from these samples were amplified real time to quantify the viruses using specific primers. Standards were also amplified along with the samples and the viral genomic copy numbers were quantitated' by standard curve analysis. Each sample was amplified in triplicates by real time PCR. The copy number of genomes of DNV and CHIKV in the presence and absence of each other was compared.
Determining the effect of co-infection of CHIKV and DNV in Ae. aegypti mosquitoes
First instar Ae. aegypti larvae after one hour of hatching, (n = 3000) were counted and transferred to a beaker containing 180 ml of water and infected with DNV by adding 20 ml of virus stock to the water in which the larvae were reared. The quantity of DNV used for infection corresponded to 6.07 × 108 viral particles. After 48 hr, the larvae were transferred to a larger pan, fed on a mixture of yeast tablet and dog biscuits and were maintained under normal insectary conditions. Many larvae died during the course of development. Among the females that emerged, 5 individuals were randomly checked by PCR for DNV infection. The remaining females were counted and were infected orally by membrane feeding on heparinised blood containing 5.73 × 105 RNA copy number of CHIKV (5 IU/ml of blood-virus mixture) as described by Harada et al. The mosquitoes were allowed to feed for 45 minutes and the fully engorged mosquitoes were segregated and counted. They were maintained under normal insectary conditions with 10% glucose as a nutritional supplement. Appropriate controls were used which include mosquitoes infected with only CHIKV, only DNV and uninfected mosquitoes.
Twelve mosquitoes were collected and stored in 2 batches of 6 mosquitoes each on alternate day post infection starting from day 0 up to day 8. DNA was extracted from one batch to analyse for DNV and RNA from the other batch to analyse for CHIKV.
Nucleic acid extraction was done from individual mosquitoes by triturating the whole mosquito in Nuclease free water. The nucleic acids were then quantitated by real time amplification along with the standards using the corresponding primers and analyzed by plotting standard curve. The real time amplification of each sample was performed in triplicates.
SYBR Green I quantitative real time PCR of DNV
SYBR Green I quantitative real time PCR of CHIKV
Determining the effect of co-infection in mosquito C6/36 cells
Effect of co-infection of DNV and CHIKV in C6/36 cell lines. The values are approximated to three decimal points.
Days post infection
LOG of RNA copy number of CHIKV (± SD)
LOG of DNA copy number of DNV (± SD)
In presence of DNV
In absence of DNV
In presence of CHIKV
In absence of CHIKV
Determining the effect of co-infection in Ae. aegypti mosquitoes
Effect of co-infection of DNV and CHIKV in Ae. aegypt i mosquitoes. The values are approximated to three decimal points.
Days post infection
LOG of RNA copy number of CHIKV (± SD)
LOG of DNA copy number of DNV (± SD)
In presence of DNV
In absence of DNV
In presence of CHIKV
In absence of CHIKV
DNV are among the most important mosquito pathogenic viruses found in natural mosquito populations. Studies at the laboratory level show them to be highly detrimental to the mosquito populations when early instar larvae are infected. The virus has also been used as a transducing agent to deliver genes of interest into the mosquitoes. Infection with DNV reduces the life span of mosquitoes there by altering its vectorial capacity. There were also a few reports about the virus altering the vector competence of mosquitoes. When AalDNV infected cultures or Ae.albopictus mosquitoes were super-challenged with DENV-2, the cell lines showed lower CPE and the mosquitoes showed lesser number of DENV [18, 20]. Later studies showed that DENV and DNV can stably co-exist in cell lines. When these dually infected cells were super challenged with JEV, no signs of cytopathology were observed [21, 22]. Given the re-emergence of CHIKV and its potential to spread far and wide, it would be of significance to study whether infection with DNV alters the susceptibility of Ae. aegypti to CHIKV.
The mosquito cell lines and mosquitoes were first infected with a known quantity of DNV. This was followed by incubation for a certain interval of time so that DNV could establish itself in the system. Then the cells and individual mosquitoes were infected with CHIKV. Throughout this study we have quantitated the genomic copy number of the viruses by real time amplification, which in turn implies to the quantity of the virus in each sample. Direct methods of quantitation like plaque assays and TCID50 are not feasible for DNV  since they do not produce morphological changes in cells. For appropriate comparison, it is required that CHIKV is also quantitated through the same means.
This study shows that CHIKV neither triggers the replication of DNV nor is its own replication suppressed due to co-infection. Ae. aegypti mosquitoes with DNV infection were found to be as susceptible to CHIKV infection as uninfected or lowly infected natural population of mosquitoes . There was an increase in the copy number of CHIKV on each day post infection as expected. But the copy number of CHIKV remained more or less constant between the DNV infected and uninfected group on any given day i.e. the multiplication of CHIKV was not affected by presence of DNV in the same system. There was no increase in the levels of DNV over time. This could be because DNV attains maximum growth on day 3 post infection . In our study, we have done superinfection on day 3 post DNV infection. By this time DNV would have attained maximum growth and hence the level of DNA copy number remained more or less constant on the following days (table 2).
The difference of around 0.1-0.8 log in the genomic copy number of the viruses in the presence and absence of co-infection could be attributed to the individual mosquito rather than any real change in the virus particles. This change might not reflect the actual number of infectious particles. Statistical analysis of the data by t-test also suggests that the genomic copy number of neither of the viruses was altered significantly due to co-infection in vitro and in vivo (p > 0.01).
A phenomenon called viral accommodation has been reported in shrimp where in natural multiple infections are possible [29, 30]. We observed a similar phenomena here when mosquitoes were co-infected. The viruses could stably co-exist both in the cell lines and adult mosquitoes. This is further ascertained by the fact that DNV was found infecting natural populations of Ae. aegypti across India including states like Kerala, Andhra Pardesh, Karnataka and Gujarat , which faced severe CHIKV outbreak. The vector mosquitoes continued to transmit the arbovirus despite DNV infection. DNV being a DNA virus has its replication in nucleus whereas CHIKV being a RNA virus replicates in cytoplasm and probably because of this, their replication might not interfere with that of each other.
This study shows that neither CHIKV nor DNV influence the multiplication of each other in mosquito cell lines or mosquitoes.
The authors acknowledge Dr. A.C. Mishra, Director, National Institute of Virology, for his critical comments on the manuscript. The authors thank Mr. K. A. Pillay and Mr. Rajen Lakra for technical support during this study.
AS was supported by fellowship from Department of Biotechnology, Govt. of India for carrying out research work.
- Clark TB, Chapman HC: A polyhedrosis in Culex salinarius of Louisiana. J Invertebr Pathol. 1969, 13: 312-10.1016/0022-2011(69)90227-4.View ArticlePubMedGoogle Scholar
- Kittayapong P, Baisley KJ, O'Neill SL: A mosquito densovirus infecting Aedes aegypti and Aedes albopictus from Thailand. Am J Trop Med Hyg. 1999, 61: 612-617.PubMedGoogle Scholar
- Sivaram A, Barde PV, Kumar SR, Yadav P, Gokhale MD, Basu A, Mourya DT: Isolation and characterization of densonucleosis virus from Aedes aegypti mosquitoes and its distribution in India. Intervirology. 2009, 52 (1): 1-7. 10.1159/000210044.View ArticlePubMedGoogle Scholar
- Boublik Y, Jousset FX, Bergoin M: Complete nucleotide sequence and genomic organization of the Aedes albopictus parvovirus (AaPV) pathogenic for Aedes aegypti larvae. Virology. 1994, 200: 752-763. 10.1006/viro.1994.1239.View ArticlePubMedGoogle Scholar
- Rwegoshora RT, Baisley KJ, Kittayapong P: Seasonal and spatial variation in natural densovirus infection in Anopheles minimus S.L. in Thailand. Southeast Asian J Trop Med Publ Health. 2000, 31: 3-9.Google Scholar
- Jousset FX, Barreau C, Boublik Y, Cornet M: A parvo-like virus persistently infecting a C6/36 clone of Aedes albopictus mosquito cell line and pathogenic for Aedes aegypti larvae. Virus Res. 1993, 29: 99-114. 10.1016/0168-1702(93)90052-O.View ArticlePubMedGoogle Scholar
- Siegl G, Bates RC, Berns KI, Carter BJ, Kelly DC, Kurstak E, Tattersall P: Characteristics and taxonomy of Parvoviridae. Intervirology. 1985, 23 (2): 61-73. 10.1159/000149587.View ArticlePubMedGoogle Scholar
- Barreau C, Jousset FX, Bergoin M: Pathogenicity of the Aedes albopictus parvovirus (Aa PV), a denso-like virus, for Ae. aegypti mosquitoes. J Invertebr Pathol. 1996, 68: 299-309. 10.1006/jipa.1996.0100.View ArticlePubMedGoogle Scholar
- Chen S, Cheng L, Zhang Q, Lin W, Lu X, Brannan J, Zhou ZH, Zhang J: Genetic, biochemical, and structural characterization of a new densovirus isolated from a chronically infected Aedes albopictus C6/36 cell line. Virology. 2004, 318: 123-133. 10.1016/j.virol.2003.09.013.View ArticlePubMedGoogle Scholar
- Afanasiev BN, Galyov EE, Buchatsky LP, Kozlov YV: Nucleotide sequence and genomic organization of Aedes densonucleosis virus. Virology. 1991, 185: 323-336. 10.1016/0042-6822(91)90780-F.View ArticlePubMedGoogle Scholar
- Ledermann JP, Suchman EL, Black WC, Carlson JO: Infection and pathogenicity of mosquito densoviruses AeDNV, HeDNV and APeDNV in Aedes aegypti mosquitoes (Diptera: Culicidae). J Econ Entomol. 2004, 97: 1827-1835. 10.1603/0022-0493-97.6.1828.View ArticleGoogle Scholar
- Buchatsky LP, Bogdanova EN, Kuznetsova MA, Lebedinets NM, Kononko AG, Chabanenko AA, Podrezova LM: Field trial of the viral preparation Viroden on the preimaginal stages of blood-sucking mosquitoes. Meditsinskaya parazitologiya I Parazitare. Bolezni. 1987, 4: 69-71.Google Scholar
- Suchman EL, Kononko A, Plake E, Doehling M, Kleler B, Black WC, Buchatsky L, Carlson J: Effects of AeDNV infection on Aedes aegypti (L) lifespan and reproduction. Biol Control. 2006, 39: 465-473. 10.1016/j.biocontrol.2006.05.001.View ArticleGoogle Scholar
- Black WC, Severson DW: Genetics of vector competence. Biology of disease vectors. Edited by: Marquardt WC, et al. 2004, Elsivier Academic Press, California, 167-185.Google Scholar
- Reeves WC: Factors that influence the probability of epidemics of Western equine, St. Louis and California encephalitis in California. Vector Views. 1967, 14: 13-18.Google Scholar
- WHO: Vector Ecology. Tech. Rep. Ser. 1972, 501: 1-38.Google Scholar
- Mourya DT, Pidiyar V, Patole M, Gokhale MD, Shouche Y: Mid gut bacterial flora of Aedes aegypti affects the susceptibility of mosquitoes to dengue viruses. Dengue Bull. 2002, 26: 190-194.Google Scholar
- Burivong P, Pattanakitsakul SN, Thongrungkiat S, Malasit P, Flegel TW: Markedly reduced severity of Dengue virus infection in mosquito cell cultures persistently infected with Aedes albopictus densovirus (AalDNV). Virology. 2004, 329: 261-269.View ArticlePubMedGoogle Scholar
- Paterson A, Robinson E, Suchman E, Afanasiev B, Carlson J: Mosquito densonucleosis viruses cause dramatically different infection phenotypes in the C6/36 Aedes albopictus cell line. Virology. 2005, 337: 253-261. 10.1016/j.virol.2005.04.037.View ArticlePubMedGoogle Scholar
- Wei W, Shao D, Huang X, Li J, Chen H, Zhang Q, Zhang J: The pathogenicity of mosquito DNV (C6/36DNV) and its interaction with dengue virus type II in Aedes albopictus. Am J Trop Med Hyg. 2006, 5: 1118-1126.Google Scholar
- Kanthong N, Khemnu N, Sriurairatana S, Pattanakitsakul S, Malasit P, Flegel TW: Mosquito cells accommodate balanced, persistent co-infections with a densovirus and Dengue virus. Dev Comp Immunol. 2008, 32: 1063-1075. 10.1016/j.dci.2008.02.008.View ArticlePubMedGoogle Scholar
- Kanthong N, Khemnu N, Sriurairatana S, Pattanakitsakul S, Malasit P, Flegel TW: Persistent, triple-virus co-infections in mosquito cells. BMC microbiology. 2010, 10: 14-10.1186/1471-2180-10-14.PubMed CentralView ArticlePubMedGoogle Scholar
- Gratz NG: Emerging and resurging vector-borne diseases. Annu Rev Entomol. 1999, 44: 51-75. 10.1146/annurev.ento.44.1.51.View ArticlePubMedGoogle Scholar
- Charrel RN, de Lamballerie X, Raoult D: Chikungunya outbreaks-the globalization of vectorborne diseases. N Engl J Med. 2007, 356: 769-771. 10.1056/NEJMp078013.View ArticlePubMedGoogle Scholar
- Mitsuhashi J, Maramorosch K: Leafhopper tissue culture: embryonic, nymphal and imaginal tissues from aseptic insects. Contrib Boyce Thompson Inst. 1964, 22: 435-460.Google Scholar
- Hasebe F, Parquet MC, Pandey BD, Mathenge EGM, 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. 10.1002/jmv.10085.View ArticlePubMedGoogle Scholar
- Harada M, Matsuoka H, Suguri S: A convenient mosquito membrane feeding method. Med Entomol Zool. 1996, 47: 103-105.Google Scholar
- Azarkh A, Robinson E, Hirunkanokpun S, Afanasiev B, Kittayapong P, Carlson J, Corsini J: Mosquito Densonucleosis Virus non-structural protein NS2 is necessary for a productive infection. Virology. 2008, 374 (1): 128-137. 10.1016/j.virol.2007.11.035.PubMed CentralView ArticlePubMedGoogle Scholar
- Flegel TW: The shrimp response to viral pathogens. The new wave. Proceedings of the special session on sustainable shrimp aquaculture, World Aquaculture, Orlando. Edited by: Browdy CL, Jory DE. 2001, World Aquaculture Society, Boca Raton, 190-214.Google Scholar
- Flegel TW, Pasharawipas T: Active viral accommodation: a new concept for crustacean response to viral pathogens. Advances in Shrimp Biotechnology. Edited by: Flegel TW. 1998, National Center for Genetic Engineering and Biotechnology, Bangkok, 245-250.Google Scholar
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.