A high recombination rate in eusocial Hymenoptera: evidence from the common wasp Vespula vulgaris
© Sirviö et al; licensee BioMed Central Ltd. 2011
Received: 23 May 2011
Accepted: 4 November 2011
Published: 4 November 2011
High recombination rates have previously been detected in two groups of eusocial insects; honeybees and ants. In this study we estimate recombination rate in a eusocial wasp Vespula vulgaris that represents a third phylogenetic lineage within eusocial hymenopterans.
A genetic linkage map of V. vulgaris based on 210 markers shows that the total map length is 2129 cM and the recombination rate is 9.7 cM/Mb (or 103 kb/cM). The present estimate in V. vulgaris is somewhat smaller than in the honeybee Apis mellifera and intermediate between the estimates from two ant species (Acromyrmex echinatior, Pogonomyrmex rugosus). Altogether, the estimates from these eusocial species are higher than in any other insect reported so far.
The four species (V. vulgaris, A. mellifera, A. echinatior, P. rugosus) are characterized by advanced eusociality with large colonies, clear queen-worker dimorphism and well developed task specialization. They also have colonies with a single, normally multiply inseminated (polyandrous) queen. Benefits of genotypic diversity within colonies (e.g. through improved task specialization or pathogen and parasite resistance) may have selected for both polyandry and high recombination rate in such advanced eusocial insects.
Genetic relatedness, the probability of sharing genes that are identical by descent (IBD), between interacting individuals is a central variable in the genetic models of social evolution . The coefficient of relatedness can be calculated on the basis of a pedigree by following the rules of Mendelian transmission. However, the segregation of alleles among the gametes produced by a diploid individual is a stochastic process and therefore the relatedness at a specific locus is a variable with a predicted mean and variance . The variance at a single locus depends on the segregation, and the variance among loci depends on the amount of recombination. Based on these principles, Sherman suggested that the evolution of insect sociality should benefit from a high recombination rate . His argument was that recombination reduces the variance of relatedness and makes it more difficult for nepotistic recognition alleles to invade the population. This argument applies mainly to a case where the society consists of a simple family with one single-mated mother queen, and the variation of relatedness among the offspring depends only on the segregation of alleles from this single mother. Templeton  pointed out that recombination can also favor social evolution under a quantitative genetic model (through effects on the variance of inclusive fitness within broods and on the nonadditive genetic components in quantitative traits) without the restrictions imposed by single mating and hypothetical recognition alleles.
At the time when Sherman and Templeton presented their hypotheses, the most conceivable way to estimate variation in the recombination rate among species was to compare chromosome numbers. The comparison showed that eusocial species have, on average, significantly higher haploid chromosome numbers than their non-social relatives, both in the Hymenoptera and in the Isoptera . Although the comparisons could be criticized for not correcting for phylogenetic non-independence, the differences were clear. Among ants, the chromosome numbers vary widely and there is no general trend of increasing chromosome numbers in different ant lineages . However, Schmid-Hempel  detected a positive correlation between the chromosome number and colony size in ants, indicating that advanced social life might be associated with recombination.
Sherman  also stated that increasing the crossing-over rate could have the same effect as increasing the chromosome number and predicted that "eusocial species should have higher crossing-over rates than nonsocial species". Since then, estimates of the crossing-over rate have been made in a number of insects, including eusocial species. The results from two species of honeybees, Apis mellifera [7–9] and A. florea , and two species of ants, Acromyrmex echinatior and Pogonomyrmex rugosus [11, 12], showed crossing-over rates that were higher than in any other animal looked at so far, thereby lending support to Sherman's prediction. Sherman  suggested that genetic diversity underlying caste and task specialization of workers could be important for the evolution of the recombination rate. Recent studies have emphasized this possibility [11, 13, 14] even though empirical evidence is lacking. The hypothesis is based on the general effect of recombination in increasing the diversity of multilocus genotypes among the progeny. As the division of labor among eusocial insect workers can have a genetic component [e.g. 15], the hypothesis suggests that recombination helps to produce a larger number of multilocus combinations and the colony thus avoids a too narrowly specialized workforce.
Even though the available studies show a high recombination rate in the honeybees and ants, the comparison still suffers from the small number of data points which are not phylogenetically independent. Hence, the aim of the present study was to estimate the recombination rate in a representative of yet another group of highly eusocial insects, the vespid wasps, where eusociality has evolved independently from ants and bees. Our results demonstrate that the wasp Vespula vulgaris is also characterized by a high recombination rate, thereby lending further support for the hypothesis that this feature is likely shared by most advanced eusocial species.
With the error detection option of Mapmaker, our map (1760 cM) decreased by 4.9%; thus our map spans 1674 cM. The map was unsaturated because the haploid chromosome number of Vespula vulgaris (n = 25 ) was exceeded by 13 linkage groups. Since the number of linkage groups exceeds the haploid chromosome number, we added 35 cM (our threshold distance for joining the linkage groups) for each gap to cover the distance to the thirteen additional linkage groups. This resulted in a final recombinational map size of 2129 cM. However, it is notable that this underestimates the real distance and the recombination rate because the gaps in the linkage map may exceed the used threshold value of 35 cM.
The physical genome sizes estimated from the neural tissues were 219.8 ± 1.3 Mb for V. vulgaris (N = 4) and 219.6 ± 1.4 for the closely related V. germanica (N = 7). Males and females gave identical 1C estimates. As the total map length of V. vulgaris was estimated to be 2129 cM, 1 cM corresponds to approximately 103 kb, resulting in an estimated recombination rate of 9.7 cM/Mb.
Recombination by crossing-over is a genetic trait and the recombination frequency commonly varies both within and between populations . As a result, the recombination rate can evolve and respond to natural selection. Laboratory experiments have demonstrated that the recombination rate can be selected readily, and directional selection has had a tendency to elevate recombination rates in domesticated animals and plants . The findings that the recombination rate in the honeybees and ants is higher than in other insects [7–11] have led to suggestions that it has been selected for in eusocial insects. Our present results from V. vulgaris are in agreement with this suggestion even though the role of selection still remains to be demonstrated.
Beye et al.  estimated that the average recombination rate in the honeybee (Apis mellifera) genome is 19 cM/Mb with little variation among its 16 chromosomes (but with considerable variation of the local recombination rate along the genome). A somewhat lower estimate (16.0 cM/Mb) was calculated by Wilfert et al. . Comparative genetic maps indicate that the recombination rate in another honeybee species, A. florae, is similar to that in A. mellifera , These estimates are four to five times higher than the estimate from the primitively eusocial bumblebee Bombus terrestris, 4.4 cM/Mb . We earlier constructed linkage maps and estimated the recombination rates as 14.0 cM/Mb in the harvester ant Pogonomyrmex rugosus  and 6.2 cM/Mb in the leaf-cutting ant Acromyrmex echinatior . The present estimate from the wasp V. vulgaris (9.7 cM/Mb) thus falls in between the two estimates from ants.
These estimates from the advanced eusocial bees, ants and wasps are higher than in other insects . The estimates in non-social hymenopterans are within the range 2.5 - 5.4 cM/Mb (4 species of parasitoid wasps) and in other insects 0.1 - 6.1 cM/Mb (15 species) . The recombination rates in the four advanced eusocial species (honey bee, leaf-cutter ant, seed harvester ant, yellow jacket wasp) are significantly larger than in the other hymenopterans (including the bumblebee, probability of no overlap in the values is P = 0.016) or in other insects in general (P = 0.0005). Bees, ants and wasps belong to aculeate Hymenoptera and are thus not phylogenetically independent. However, the lineages have diverged a long time ago and eusociality has evolved separately in them. According to Brady et al. , ants, bees and wasps had a common ancestor about 160 Mya, and the lineages leading to ants and eusocial wasps diverged about 140 Mya. The two ants in which the recombination rate has been estimated had their common ancestor about 80 Mya , and Apis mellifera and A. florea at least 8-10 Mya  and probably about 20 Mya . Even though we cannot exclude the possibility that a high recombination rate is an ancestral state, there has been ample time for selection to modify the rates if they had any adaptive significance.
The use of AFLP-marker data deserves discussion because the methodology has gained criticism e.g. due to frequent occurrence of non-homologous fragments with the same amplicon length [23, 24]. Study on Nasonia jewel wasp  revealed 41.5% shorter map size when SNP-markers were used instead of RAPD/AFLP markers . Similarly in Bombyx mor i silk moth recombination map estimates vary from 3432 cM (simple sequence repeat; ) to 1413.4 cM (SNP; ) depending on the method and the number of the markers. On the other hand the original estimate of the honeybee map size was based on RAPD markers (3500 cM ) and the subsequent estimates based on microsatellites (4000 cM ) or genome sequencing (4553 cM ) have not decreased it. There is thus no universal trend that the RAPD/AFLP markers would overestimate recombination, and the data from most insects used in our comparisons were obtained with these methods, making the results comparable.
Sherman  suggested that a high recombination rate could be adaptive in eusocial insects either because recombination equalizes the fractions of genomes shared by colony members or because it generates a larger number of different multilocus genotypes. Sherman particularly considered the advantage of genotypic diversity in caste and task specialization, and Schmid-Hempel  suggested that the same can also apply to defense against parasites. As noted by Schmid-Hempel, both recombination and multiple mating by females increase the genotypic diversity among the offspring and can be beneficial to eusocial insect colonies under selection by parasites. The difference between the two factors is that unlike recombination, mating with many males also increases allelic diversity. Multiple mating by queens is known in many eusocial insects but the average mating frequency is generally rather low . It is noteworthy that the species in which high recombination rates have been estimated, have all monogynous societies, i.e. societies with a single queen. They also have large colonies with clear queen-worker dimorphism and elaborate division of tasks among workers, and the queens are typically highly polyandrous. The estimated number of effective matings is up to 17.6 in the honeybee A. mellifera, 1.9 in Vespula vulgaris, 4.7 in Pogonomyrmex rugosus and 5.3 in Acromyrmex echinatior . These estimates are clearly higher than the mean estimates for eusocial insects in general. One could thus suggest that these species benefit from genotypic diversity within colonies and that this has selected both for polyandry and for a high recombination rate
Sherman  initially hypothesized an association between recombination rate and sociality because of effects on genomic multilocus relatedness. The point is that the expected relatedness among full sisters (r = 0.75) in a single locus is in fact a mean between complete identity (r = 1) when the sisters received an identical allele from the mother and 'half identity' (r = 0.5) when the sisters received different maternal alleles and share only the paternal allele. Lack of recombination could result in genetic cliques within which sisters are unusually highly related over many loci. If there is any kin discrimination within the societies and nepotistic behaviour based on this discrimination, such genetic cliques could lead to nepotistic conflicts and harm the function of the colony. Nepotistic behavior based on kin discrimination has been doubted but some evidence for it has been presented recently [32, 33]. Whether the effect of recombination on the distribution of pair-wise relatednesses among colony members could affect kin recognition and discrimination remains to be studied.
The data from three major groups of eusocial Hymenoptera (ants, bees and wasps) show high recombination rates, supporting the theoretical predictions on the positive association between sociality and recombination. These predictions are derived from two separate lines of reasoning. First, recombination equalizes multilocus relatedness among brood members and prevents kin discrimination. However, kin discrimination is controversial and there is little evidence for it (e.g. ) and it is doubtful whether it could be a factor selecting for recombination. We therefore incline to the second hypothesis, according to which recombination is an important factor generating multilocus genotypic diversity within a society [3, 30] and therefore highly relevant for the function of advanced eusocial colonies [11, 13, 14]. It can be hypothesized that benefits of intracolonial genotypic diversity has selected for both polyandry and high recombination rate. Hughes et al  found a negative correlation between polyandry and polygyny in data from 241 eusocial insect species. This suggests that polyandry and polygyny are alternative ways to increase the genetic diversity within societies. Consequently, we predict a negative correlation also between recombination rate and polygyny under the diversity hypothesis. Data from non-social species as well as from eusocial species with different colonial types (polygynous colonies, small colony sizes) are required for comparative tests of the adaptive role of recombination in advcanced eusocial insects.
Physical genome size
The size of the physical genome was determined from neural tissue of four Vespula vulgaris and seven V. germanica individuals (males and females), collected from nests that were dug up in Leuven in August 2009. Samples were prepared for flow cytometric genome size determination as described in . Neural tissue dissected from the head of each sampled specimen was placed along with a single head of a female Drosophila virilis into 1 ml of cold Galbraith buffer in a 2 ml Kontes Dounce tissue grinder, stroked 15 times with the "A" pestle and then filtered through 20 μm nylon mesh. Propidium iodide was added to each sample to a final concentration of 50 ppm and allowed to stain 40 minutes in the dark on ice. To determine the relative positions of sample and standard fluorescence peaks and test for possible artifacts of preparation , similar preparations were made with at least one insect of each species alone. The mean fluorescence of stained nuclei in replicate samples was quantified using a Partec Cyflow cytometer with a 100 mw green laser tuned at 532 nm. Fluorescence at >615 nm was detected by a photomultiplier screened by a long pass filter. To ensure that scoring included only intact nuclei free from cytoplasmic tags, counting was activated by red fluorescence, and only nuclei with low forward and low side scatter were included in the analysis. Samples were run to produce a total of at least 1000 nuclei under each scored peak. DNA content was determined from co-preparations of sample and standard by multiplying the ratio of the mean peak fluorescence of the 2C sample to the 2C mean fluorescent peak of D. virilis, times the genome size of the standard (328 Mbp for D. virilis, based on co-preparations with D. melanogaster, n = 30).
DNA extraction and quality, AFLP
Variable AFLP markers in haploid males were used to build linkage groups for V. vulgaris by using Mapmaker version 3.0 . This was done using the phase unknown procedure as described in [11, 40]. The Kosambi mapping function was used to transform the recombination fractions into centimorgans. This procedure reduces the map length by taking into account possible double recombination events. In two-point analysis, the preliminary linkage groups were formed with the thresholds of LOD score 5.0 and maximum marker distance 30 cM, and the order of markers within linkage groups was confirmed in multipoint analysis with the same criteria. When several other markers were already clustered and ordered in linkage groups, unlinked markers were added in the end parts of the groups with the criteria LOD 3.0 and distance 35 cM by using "near"-command that is based on two-point analyses . Finally, some of the existing groups were combined if the distal markers were within 35 cM distance and the statistical threshold LOD score was at least 3.0. The order within these newly formed groups was tested again and the final marker linkage groups were formed with the criteria LOD 3.0 and 35 cM. The final map was then run through a genotyping error detection process of Mapmaker 3.0, with the error detection option switched on with an a priori error probability of 1% to detect the presence of multiple crossing-overs in linkage groups. As a final check, the raw data were inspected manually, and arranged according to the order of markers to allow visual identification of possible falsely scored individuals and by rechecking the interpretation of marker peaks for the loci and individuals.
Acknowledgements and Funding
We thank Riitta Jokela for help with the laboratory analyses. The work has been supported by a grant from the Academy of Finland (122210 to P.P.) and by the Finnish Graduate School of Population Genetics (to A.S.) and by Finnish Cultural Foundation (to A.S.).
- Crozier RH, Pamilo P: Evolution of Social Insect Colonies. 1996, Oxford: Oxford University PressGoogle Scholar
- Barash DP, Holmes WG, Greene PJ: Exact versus probabilistic coefficient of relationship: some implications for sociobiology. Am Nat. 1978, 112: 355-363. 10.1086/283278.View ArticleGoogle Scholar
- Sherman PW: Insect chromosome numbers and eusociality. Am Nat. 1979, 113: 925-935. 10.1086/283445.View ArticleGoogle Scholar
- Templeton AR: Chromosome number, quantitative genetics and eusociality. Am Nat. 1979, 113: 937-941. 10.1086/283446.View ArticleGoogle Scholar
- Lorite P, Palomeque T: Karyotype evolution in ants (Hymenoptera: Formicidae) with a review of the known ant chromosome numbers. Myrmecol News. 2010, 13: 89-102.Google Scholar
- Schmid-Hempel P: Parasites in Social Insects. 1998, Princeton: Princeton University PressGoogle Scholar
- Hunt GJ, Page RE: Linkage map of the honey bee, Apis mellifera, based on RAPD markers. Genetics. 1995, 139: 1371-1382.PubMed CentralPubMedGoogle Scholar
- Solignac M, Vautrin D, Baudry E, Mougel F, Loiseau A, Cornuet JM: A microsatellite-based linkage map of the honeybee, Apis mellifera L. Genetics. 2004, 167: 253-262. 10.1534/genetics.167.1.253.PubMed CentralView ArticlePubMedGoogle Scholar
- Beye M, Gattermeier I, Hasselmann M, Gempe T, Schioett M, Baines JF: Exceptionally high levels of recombination across the honey bee genome. Genome Res. 2006, 16: 1339-1344. 10.1101/gr.5680406.PubMed CentralView ArticlePubMedGoogle Scholar
- Meznar ER, Gadau J, Koeniger N, Rueppell O: Comparative linkage mapping suggests a high recombination rate in all honeybees. J Hered. 2010, 101: S118-S126. 10.1093/jhered/esq002.View ArticlePubMedGoogle Scholar
- Sirviö A, Gadau J, Rueppell O, Lamatsch D, Boomsma JJ, Pamilo P: High recombination frequency creates genotypic diversity in colonies of the leaf-cutting ant Acromyrmex echinatior. J Evol Biol. 2006, 19: 1475-1485. 10.1111/j.1420-9101.2006.01131.x.View ArticlePubMedGoogle Scholar
- Sirviö A, Pamilo P, Johnson RA, Page RE, Gadau J: Origin and evolution of the dependent lineages in the genetic caste determination system of Pogonomyrmex ants. Evolution. 2011, 65: 869-884. 10.1111/j.1558-5646.2010.01170.x.View ArticlePubMedGoogle Scholar
- Gadau J, Page RE, Werren JH, Schmid-Hempel P: Genome organization and social evolution in Hymenoptera. Naturwiss. 2000, 87: 87-89. 10.1007/s001140050016.View ArticlePubMedGoogle Scholar
- Wilfert L, Gadau J, Schmid-Hempel P: Variation in genomic recombination rates among animal taxa and the case of social insects. Heredity. 2007, 98: 189-197. 10.1038/sj.hdy.6800950.View ArticlePubMedGoogle Scholar
- Robinson GE, Fernald RD, Clayton DF: Genes and social behavior. Science. 2008, 322: 896-900. 10.1126/science.1159277.PubMed CentralView ArticlePubMedGoogle Scholar
- Lincoln SE, Daly M, Lander ES: Constructing genetic maps with MAPMAKER/EXP 3.0. 1992, Cambridge MA: Whitehead Institute Technical ReportGoogle Scholar
- Hoshiba G, Matsuura M, Imai HT: Karyotype evolution in the social wasps (Hymenoptera, Vespidae). Jpn J Genet. 1989, 64: 209-222. 10.1266/jjg.64.209.View ArticleGoogle Scholar
- Butlin RK: Recombination and speciation. Mol Ecol. 2005, 14: 2621-2635. 10.1111/j.1365-294X.2005.02617.x.View ArticlePubMedGoogle Scholar
- Wilfert L, Gadau J, Schmid-Hempel P: A core linkage map of the bumblebee Bombus tererstris. Genome. 2006, 49: 1215-1226. 10.1139/g06-075.View ArticlePubMedGoogle Scholar
- Brady SG, Larkin L, Danforth BN: Bees, ants, and stinging wasps (Aculeata). The timetree of life. Edited by: Hedges SB, Kumar S. 2009, Oxford: Oxford University Press, 264-269.Google Scholar
- Brady SG, Schultz TR, Fisher BL, Ward PS: Evaluating alternative hypotheses for the early evolution and diversification of ants. Proc Natl Acad Sci USA. 2006, 103: 18172-18177. 10.1073/pnas.0605858103.PubMed CentralView ArticlePubMedGoogle Scholar
- Cardinal S, Straka J, Danforth BN: Comprehensive phylogeny of apid bees reveals the evolutionary origins and antiquity of cleptoparasitism. Proc Natl Acad Sci USA. 2010, 107: 16207-16211. 10.1073/pnas.1006299107.PubMed CentralView ArticlePubMedGoogle Scholar
- Paris M, Bonnes B, Ficetola GF, Poncet BN, Despres L: Amplified fragment length homoplasy: in silico analysis for model and non-model species. BMC Genomics. 2010, 11: 287-10.1186/1471-2164-11-287.PubMed CentralView ArticlePubMedGoogle Scholar
- Gort G, van Eeuwijk FA: Review and simulation of homoplasy and collision in AFLP. Euphytica. 2011Google Scholar
- Niehuis O, Gibson JD, Rosenberg MS, Pannebakker BA, Kovoets T: Recombination and its impact on the genome of the haplodiploid parasitoid wasp Nasonia. PLoS ONE. 2010, 5: e8597-10.1371/journal.pone.0008597.PubMed CentralView ArticlePubMedGoogle Scholar
- Gadau J, Werren JH, Page RE: Mapping of hybrid incompatibility loci in Nasonia. Genetics. 1999, 153: 1731-1741.PubMed CentralPubMedGoogle Scholar
- Miao XX, Xub SJ, Li MH, Li MW, Huang JH: Simple sequence repeat-based consensus linkage map of Bombyx mori. Proc Natl Acad Sci USA. 2006, 102: 16303-16308.View ArticleGoogle Scholar
- Yamamoto K, Nohata J, Kadono-Okuda K, Narukawa J, Sasanuma M: A BAC-based integrated linkage map of the silkworm Bomby mori. Genome Biology. 2008, 9: R21-10.1186/gb-2008-9-1-r21.PubMed CentralView ArticlePubMedGoogle Scholar
- Solignac M, Mougel F, Vautrin D, Monnerot M, Cornuet JM: A third-generation microsatellite-based linkage map of the honey bee, Apis mellifera, and its comparison with the sequence-based physical map. Genome Biology. 2007, 8: R66-10.1186/gb-2007-8-4-r66.PubMed CentralView ArticlePubMedGoogle Scholar
- Schmid-Hempel P: Mating, parasites and other trials of life in social insects. Microbes Infect. 2000, 2: 515-520. 10.1016/S1286-4579(00)00316-6.View ArticlePubMedGoogle Scholar
- Boomsma JJ, Kronauer DJC, Pedersen JS: The evolution of social insect mating systems. Organization of Insect Societies. Edited by: Gadau J, Fewell JH. 2009, Cambridge MA: Harvard University Press, 3-25.Google Scholar
- Hannonen M, Sundström L: Worker nepotism among polygynous ant. Nature. 2003, 421: 910-10.1038/421910a.View ArticlePubMedGoogle Scholar
- Wenseleers T: Nepotism absent in insect societies - or is it?. Mol Ecol. 2007, 16: 3063-3065. 10.1111/j.1365-294X.2007.03313.x.View ArticlePubMedGoogle Scholar
- Meunier J, Delaplace L, Chapuisat M: Reproductive conflicts and egg discrimination in a socially polymorphic ant. Behav Ecol Sociobiol. 2010, 64: 1655-1663. 10.1007/s00265-010-0979-z.View ArticleGoogle Scholar
- Hughes WOH, Ratnieks FLW, Oldroyd BP: Multiple paternity of multiple queens: two route to greater intracolonial genetic diversity in the eusocial Hymenoptera. J Evol Biol. 2008, 21: 1090-1095. 10.1111/j.1420-9101.2008.01532.x.View ArticlePubMedGoogle Scholar
- Johnston JS, Ross LD, Beani L, Hughes DP, Kathirithamby J: Tiny genomes and endoreduplication in Strepsiptera. Insect Mol Biol. 2004, 13: 581-585. 10.1111/j.0962-1075.2004.00514.x.View ArticlePubMedGoogle Scholar
- Bennett MD, Price HJ, Johnston JS: Anthocyanin inhibits propidium iodide DNA fluorescence in Euphorbia pulcherrima: Implications for genome size variation and flow cytometry. Ann Bot. 2008, 101: 777-790. 10.1093/aob/mcm303.PubMed CentralView ArticlePubMedGoogle Scholar
- Daly D, Archer ME, Watts PC, Speed MP, Hughes MR, Barker FS: Polymorphic microsatellite loci for eusocial wasps (Hymenoptera: Vespidae). Mol Ecol Notes. 2002, 2: 273-275.Google Scholar
- Vos P, Hogers R, Bleeker M, Reijans M, Van de Lee T, Hornes M: AFLP: a new technique for DNA fingerprinting. Nucleic Acids Res. 1995, 23: 4407-4414. 10.1093/nar/23.21.4407.PubMed CentralView ArticlePubMedGoogle Scholar
- Gadau J, Gerloff CU, Kruger N, Chan H, Schid-Hempel P, Wille A: A linkage analysis of sex determination in Bombus terrestris (L.) (Hymenoptera: Apidae). Heredity. 2001, 87: 234-242. 10.1046/j.1365-2540.2001.00919.x.View ArticlePubMedGoogle 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.