Diverse MHC IIB allele repertoire increases parasite resistance and body condition in the Long-tailed giant rat (Leopoldamys sabanus)
© Lenz et al; licensee BioMed Central Ltd. 2009
Received: 28 October 2008
Accepted: 23 November 2009
Published: 23 November 2009
Genes of the major histocompatibility complex (MHC) code for key functions in the adaptive immune response of vertebrates and most of them show exceptionally high polymorphism. This polymorphism has been associated with the selection by diverse and changing parasite communities. We analysed MHC class IIB diversity, gastrointestinal parasite load and body condition in the wild ranging tropical rat Leopoldamys sabanus (Thomas, 1887) under natural selection conditions in a highly variable rainforest environment in Borneo to explore the mechanisms that maintain these high levels of genetic polymorphism.
Allelic diversity was determined via SSCP and sequencing, and parasite screening was done through non-invasive faecal egg count. The detected alleles showed expected high levels of polymorphism and balancing selection. Besides a clear advantage for more diverse MHC genotypes in terms of number of alleles, reflected in better body condition and resistance against helminth infection, our data also suggested a positive effect of MHC allele divergence within an individual on these parameters.
In accordance with the heterozygote advantage hypothesis, this study provides evidence for an advantage of more diverse MHC genotypes. More specifically, the potential negative relation between individual allele divergence and number of parasite species is in line with the 'divergent allele advantage' hypothesis.
Genetic variation is the material for evolutionary processes. Its maintenance may to a large part be driven by the exceptional dynamics of host-parasite co-evolution [1, 2]. The major histocompatibility complex (MHC) represents a crucial part of the host side in these co-evolutionary processes and is one of the most polymorphic systems in the vertebrate genome. It has been studied extensively in model species under laboratory conditions , but also became the focus of an increasing number of studies on natural populations because of its functional importance for the immune system and mate choice [4–6]. Its genes code for cell surface molecules that present self and non-self antigens to T-cells . This function makes them a key factor of the adaptive immune system as it triggers a respective targeted response upon recognition of non-self antigens and therefore plays a vital role in the recognition of pathogens invading the body. The high polymorphism in this system, which is composed of a high total number of alleles as well as allelic divergence on the sequence level, is thought to facilitate populations to persist in and adapt to natural habitats [5, 8].
A number of hypotheses have been proposed to explain the forces that counteract the effects of random genetic drift and fixation [reviewed in [5, 6, 9–11]]. The 'negative frequency-dependent selection' hypothesis [12–14] grounds on the effects of host parasite co-evolution ['Red Queen hypothesis', ], stating that frequencies of alleles constantly change with the frequency of adapted and non-adapted pathogens. Rare alleles are being preserved within a population through such fluctuations or the emergence of new pathogens. The prerequisite for this co-evolution in terms of reciprocal adaptation is the direct interaction of individual alleles and distinct pathogenic morphotypes. Growing evidence comes from a variety of taxa that such an interaction exists [16–23]. Another hypothesis, the 'heterozygote advantage' has been proposed as a general mechanism to maintain high allele numbers in populations by an advantage for individuals carrying more than one allele, enabling them to present a broader range of antigens [24, 25]. There is experimental evidence that this is valid for animals exposed to multiple pathogens [, up to an optimum in ] but not to single pathogens . This seems plausible, as the advantage would most likely result from a combined resistance of individual alleles each presenting a certain range of parasites [29, 30]. Nevertheless, the exact mechanisms are still a matter of debate [31, 32]. De Boer et al.  used mathematical models to study the degree of MHC polymorphism arising when 'heterozygote advantage' is the only selection pressure. The simulations revealed that the advantage for an individual that simply carries two alleles instead of one is not sufficient to explain the high population diversity of the MHC. Wakeland et al.  proposed a special case of the heterozygote advantage, the 'divergent allele advantage', which refers to the genetic distance between alleles and expects an advantage for more divergent allele combinations, enabling the respective carrier to present more different antigens to its adaptive immune system. This mechanism may potentially explain a directional selection for both, a large number of alleles at the population level and a high divergence of alleles at the individual level.
Tropical rainforests are known for their exceptional species diversity  and this may also be true for the diversity of parasites an inhabitant has to face. Therefore studying a species from such diverse environment could provide a great opportunity to better understand the ongoing evolution of the MHC, as selection mechanisms should be pronounced. In this study we examine the variability in the exon 2 of MHC class IIB genes and its effect on parasite infection of the Long-tailed giant rat (Leopoldamys sabanus) to unravel the mechanisms of selection in a wild ranging species. This large nocturnal murid species lives semi-arboreal in tropical forested habitats and reaches an average body mass of 368 g , feeding on plant material and arthropods . This rodent was chosen, because it is a wide-spread species living in the diverse and dynamic environment of dipterocarp rainforests across Southeast Asia [38, 39], presumably conferring the need to constantly adapt to environmental changes such as the encounter of new parasites, but also to persist against ancient pathogen lines. The presumably high parasite pressure in these habitats could lead to pronounced and therefore detectable signs of elsewhere more subtle mechanisms of selection.
MHC II gene variability
In total we determined the MHC class IIB genotype of 46 individuals from the two different subpopulations. Overall 28 different sequence variants (alleles) could be distinguished via SSCP and confirmed by sequencing. The genomic origin of the alleles was investigated through a protein BLAST search on the NCBI database . Ten of the detected alleles showed unique DRB origin, blasting to different murid species, ten other alleles produced hits with DRB as well as DQB alleles from different murid species and the last eight alleles produced mainly DQB hits from several murid taxa. The same latter eight alleles showed a deletion of two amino acids, leading to a shortened amino acid sequence, but not to a frame shift in codon translation. One of the two deleted amino acid sites (position 65) is thought to be involved in antigen binding . The same deletion can also be found in some (but not all) published DQB alleles of the Brown rat (Rattus norvegicus; e.g. [GenBank:AY626188]) and House mouse (Mus musculus; e.g. [GenBank:AY740472]).
The two populations shared 82% of those alleles, which were present in more than two individuals (9 out of 11) and none of the private alleles, i.e. unique to one of the study sites, occurred in more than 3 individuals (Fig. 2). Genotypes ranged from one to four alleles per individual. We calculated the divergence of alleles for each individual, which ranged from 0 to 0.404 and correlated strongly with the number of alleles (nonparametric correlation, Spearman's rho = 0.58, p < 0.001).
Evidence for selection
Signatures of selection on the MHC class IIB loci.
# of aa sites
n = 28
0.44 ± 0.08
0.06 ± 0.03
0.13 ± 0.03
0.13 ± 0.04
0.20 ± 0.03
0.11 ± 0.03
Alleles with deletion
n = 8
0.19 ± 0.05
0.00 ± 0.00
0.09 ± 0.02
0.09 ± 0.04
0.11 ± 0.02
0.06 ± 0.03
Effects of MHC diversity on body condition and parasite load
We found a significant negative influence of allele divergence on the number of helminth morphotypes (parasite species richness) within an individual (GLM; X21,38 = 2.42, p = 0.042; Fig. 3b). ABS divergence showed a similar trend (GLM; X21,38 = 2.16, p = 0.056), while neither the number of alleles (GLM; X21,37 = 1.50, p = 0.11) nor body condition (GLM; X21,34 = 0.01, p = 0.97) had an effect on the number of helminth morphotypes. Again, a part of the association between allele divergence and parasite species richness was potentially due to the difference between individuals with one and those with more than one allele, since the latter showed a lower parasite species richness (Mann-Whitney U, p = 0.013). The models showed also that the number of helminth morphotypes differed significantly between the cohorts (GLM; X21,39 = 6.18, p = 0.014). The residuals of all four models were well normally distributed (Shapiro-Wilk test; all p ≥ 0.40).
When testing, whether the intensity of nematode infection (nematode EPG) was influenced by the presence of individual MHC-alleles, the parameters body condition, allele divergence and allele number were removed from the respective models according to AIC and therefore had no significant influence. In contrast, intensity of nematode infection differed between age classes, being higher in adult individuals, and was higher in Poring (negative binomial GLM; both p < 0.05). For statistical reasons we only tested the effect of those alleles that occurred at least in six individuals (6 alleles). Individuals with the second most common allele Lesa-MHCIIB*13 tended to show a decreased intensity of nematode infection (negative binomial GLM, X21,41 = 5.076, p = 0.024; not significant with Bonferroni corrected p = 0.008 for multiple testing). The intensity of nematode infection was independent of the sex (Mann-Whitney U, p = 0.46).
In this study, the variability of the MHC class IIB genes in the Long-tailed giant rat (Leopoldamys sabanus) was investigated and its effect on the parasite infection analysed. The aim was to further unravel the evolutionary mechanisms that lead to the observed high polymorphism in MHC genes by investigating a natural population in the diverse and dynamic environment of a tropical rainforest. As expected from studies on other small mammals, the MHC of the Long-tailed giant rat was highly polymorphic. This was shown in the diverse allele pool of the two subpopulations (28 distinct alleles in 46 individuals) as well as in the high percentage (70%) of polymorphic amino acid sites within the exon 2, the domain that is responsible for antigen binding .
The inconsistent BLAST results might indicate that MHC IIB locus differentiation is less consolidated than sometimes assumed. Mammals are the only order where the differentiation between DRB and DQB loci is actually being made, while in other vertebrate orders the different MHCIIB loci are usually indistinguishable. Most studies use the NCBI BLAST search for identification of locus origin and this might be sometimes misleading, as locus identification is commonly referred from distantly related taxa if no species specific sequence data is available. The so identified sequences then function again for later studies as reference and artificially reinforce locus determination. BLASTing for instance a published DRB1 allele of Canis familiaris [GenBank:AF016910] to the dog genome  brings as first hit the DLA DQB1 locus. Additionally there is support for the occurrence of interlocus exchange of sequence motifs [45, 46], which could also lead to misleading BLAST results.
The range of one to four alleles per individual suggests that the number of MHCIIB loci differs between haplotypes, albeit it is theoretically possible that the same allele occurs on several loci. The variation in MHCIIB loci number has already been resolved for several vertebrate species, including mammals [47–51] and seems to occur in the Long-tailed giant rat as well. Another explanation would be the presence of null alleles, meaning that we did not amplify all present alleles. However, the used primers are rather conservative, have been established before in several mammalian species and produced more than two alleles in some of them, indicating duplicated MHCIIB loci as well [22, 52–56]. The puzzling deletion in some of the alleles is not species specific and can be found in some (but not all) DQB alleles of the Brown rat  and the House mouse (e.g. [GenBank:AY740472]). However this sequence characteristic alone is not enough to refer locus origin, as it is not deterministic for one of the loci. The Neighbour-joining tree does not solve the question of locus origin either, as none of the Lesa allele groups gets significant bootstrap support. These results suggest that the divergence between the two detected MHC class II loci is either hindered by frequent inter-locus recombination or by the fact that both loci evolve under the same selective forces. Under these circumstances the apparent lack of traces of inter- and intra-locus recombination events therefore suggests that both loci code for molecules that provide the same function and are shaped by the same selective forces. As done in most studies that analyse MHC diversity in natural populations, we amplified only the exon 2 of MHC IIB genes, since this codes for the antigen binding region of the MHC molecule. Despite high functional and structural similarity of the molecules and presumable paralogous origin, the different MHC IIB loci might have diverged more strongly in other exons. Restricting the comparison to the exon 2 sequence only might therefore contribute to the dissatisfying locus resolution. However, a recent meta-analysis on MHC in several rodent taxa, some of them with duplicated loci, confirmed that the overall combined number of alleles from several class II loci is shaped by local parasite species richness .
The observed high ratio of non-synonymous to synonymous substitutions is a strong sign for balancing selection and strengthens the assumption that both loci are functional and currently involved in antigen presentation, which is not always the case as seen in a recent study on MHC genes in bank voles . This is supported by CODEML, which showed that the positive selection model M8 explained the obtained sequence data better than the model for neutrality M7. Overall, our results do not show any sign for non-functionality or pseudogenes and indicate that the analysed genes in this tropical rodent behave like classical MHC genes. We consequently assume that they are actively involved in antigen presentation, providing resistance against pathogens and confer a crucial adaptive trait for the Long-tailed giant rat. Therefore we analysed possible adaptive effects of the individual MHC-constitution on overall body condition as well as helminth parasite infestation to shed more light onto the selection mechanisms that lead to such high polymorphism.
Interestingly, the allele divergence (average genetic distance of alleles on the amino acid level within an individual) was correlated with a lower individual number of parasite morphotypes and a better overall body condition. This is in line with the 'divergent allele advantage' hypothesis , which assumes that individuals with more divergent alleles have an advantage, because their immune system can counteract against more different pathogens and therefore leads to 'overdominant selection' [29, 30]. Other studies have found an unexpected high divergence of MHC alleles in natural populations  and also signs for 'divergent allele advantage' in the context of mate choice [60–63]. Consuegra & Garcia de Leaniz  also found a positive effect of MHC dissimilarity on the prevalence of a single parasite. However, the current data set does not provide enough statistical power to disentangle the effects of allele divergence and pure heterozygote advantage.
The question remains whether the sequence divergence between alleles is a better functional predictor for immunocompetence than the plain number of alleles. It makes intuitively sense to give more importance to allele pairs that differ by 30 amino acids than to those that differ by only one amino acid. Studies that experimentally support the idea of overdominance in allele numbers [e.g. , up to an optimum in ], may in fact show indirectly the effect of allele divergence, which is supported by our observation that the genetic divergence within an individual correlates strongly with the number of alleles.
The genetic divergence increases strongest with the step from one allele to two alleles. Every additional allele thereafter only adds a fraction of its antigen binding potential, since alleles are known to overlap in their antigen repertoire  and each new allele only adds the fraction of antigens that are not yet bound by the already present alleles. The observed correlation between MHC allele divergence and parasite species richness might to a significant part be due to the observed difference between individuals with one allele and those with more alleles. This is, however, not the only reason for the correlation, since we would then also expect a stronger correlation of parasite species richness and mere allele number, which was not detected.
Our finding of a potential association between intensity of nematode infection and a certain allele, which looses significance after Bonferroni-correction, would support the 'negative frequency-dependent selection' hypothesis, which has already been shown in several other host parasite systems [16, 18, 20, 54, 65, 66]. The fact that the potentially 'resistant' allele Lesa-MHCIIB*13 is the second most common in our population, although a rather rare allele would be expected by this hypothesis, could be due to the sampling time point. The cycling pattern of frequencies of both alleles and parasites due to host-parasite coevolution can lead to a spread of an advantageous and formerly rare allele until the parasite evolves to avoid it . The 'negative frequency-dependent selection' hypothesis could also explain the maintenance of the observed high number of alleles that were found in one or two individuals only. It expects a potential advantage for individuals carrying rare alleles, which enable them to better resist new pathogen threats and therefore keeps those rare alleles in the gene pool instead of purging them.
The negative correlation of mean genetic distance between MHC IIB alleles and the number of parasite morphotypes could be interpreted in support of the 'divergent allele advantage', which has to our knowledge not yet been fully proven by empirical studies. However, the main effect originates likely from the difference between individuals with one and those with more alleles. In addition we found a substantial number of rare alleles and hints for an interaction between a specific allele and nematode infection intensity. We therefore conclude that both hypotheses, the heterozygosity advantage (and potentially the 'divergent allele advantage' as a special case thereof) as well as the 'negative frequency-dependent selection', work non-exclusively on the maintenance of a high variability in the MHC allele pool, as it has already been proposed in theory by Apanius et al. . More work has to be done on the molecular level to discover and unravel the interrelations of individual alleles and individual parasite genotypes to finally resolve the mystery of balancing selection on the MHC loci. A special focus should be put on the importance of individual allele divergence, e.g. by looking at a broader parasite range.
Study area and animal collection
The samples were collected in two lowland rainforest sites near Mount Kinabalu in the federal state of Sabah in Malaysia, North-western Borneo. While one site ('Poring') was located in an old-growth forest area within the Kinabalu National Park (6°22'N, 116°42'E), the other site ('Monggis') was located in a distance of ca. 21 km beeline in an adjacent logged forest area near the village Monggis (6°13'N, 116°45'E), which has been logged about 20 years ago. Both forest areas comprised more than 1000 hectares. Animals were captured in four trapping sessions (cohorts) between 2002 and 2004 with wire-mesh live traps placed along transects of 380 meter length and additional traps as described in Wells et al. .
Traps were baited with ripe banana and checked every morning. Captured individuals were weighed and first-capture individuals were anaesthetised with diethyl ether for further handling. We injected a subcutaneous transponder (AEG Identification Systems) for permanent identification. Then animals were aged, sexed and morphological measurements (weight, body length) were taken. We sampled ear tissue with a biopsy ear punch and stored it immediately in 95% ethanol. Faeces were collected from the bottom of the trap and stored in 2-3% formalin. Animals were handled and immediately released at the point of capture. The animal handling and sampling protocol followed guidelines of the American society of Mammalogists  and was approved by the Economic Planning Unit Malaysia and the Sabah Parks authority.
We focused on helminth parasite species because their prevalence and the intensity of infestation can be assessed non-invasively from faecal samples. We counted all eggs from nematodes, cestodes and trematodes (Plathelminthes) from faeces with a modified flotation and McMaster method [19, 68, 70]. This non-invasive technique has been shown to be accurate for quantification of helminth eggs (e.g. [71–73]). Samples of approximately 600 mg faeces were dissolved in 9 ml potassium iodide solution (specific gravity 1.5 g/ml), sieved to remove large debris and screened for helminth eggs by counting the content of two chambers of a McMaster slide. We calculated infection intensity as eggs per gram faeces (EPG). All eggs were photographed and measured (Zeiss, AxioCam and AxioVision software; 10-40× amplification). Images were assigned to operational taxonomic units (orders for nematodes, cestodes) based on features of egg shell and plasma, and further distinguished into morphotypes by size classes and shell thickness for strongyle nematodes [68, 74]. The overall number of helminth morphotypes within each rat was then used as a measure of parasite species richness.
DNA was extracted according to the mouse tail protocol using the DNeasy Tissue Kit (Qiagen, Hilden). We aimed at amplifying the highly variable exon 2 of the MHC class IIB genes. As there was no sequence data available for the Long-tailed giant rat, we used the well established primers JS1: 5'-GTGTCATTTCTACAACGGGACG-3'  and GH50: 5'-CTCCCCAACCCCGTAGTTGTGTCTGCA-3' , which have been established for DRB loci in small mammals. When BLASTed against the NCBI nr database, the first 1000 hits for both primers deliver only MHC II DRB genes from different species (searching for distant homologies and adjusting for short sequence queries). These primers lead to a fragment of 221 bp length, which spans most of the exon 2, including the most polymorphic part of the MHC IIB gene. Reactions of 20 μl contained 2-4 μl template DNA, 2 μl 10× incubation buffer incl. MgCl2, 0.175 mM of each dNTP, 1 unit taq polymerase (all QBioGene, Irvine) and 0.375 mM of each primer (MWG Biotech, Ebersberg). PCR was performed on a TGradient Thermocycler (Biometra, Göttingen) with an initial denaturing step of 96°C for 2 min, 34 cycles of 96°C for 30 s, 54.7°C for 60 s and 72°C for 60 s and a final elongation step of 72°C for 10 min. We checked the PCR product for quality and size on a 1.5% agarose gel and then purified using the QIAquick Gel Extraction Kit (Qiagen, Hilden) according to the manufacturers protocol, but eluting in dH2O. An additional ethanol precipitation was done to further purify the product.
Alleles were separated using Single Stranded Conformation Polymorphism (SSCP ). It is a sensitive method to distinguish even minimal allele differences [76, 77], which has been widely used in human genetics and became popular in population genetics and evolutionary ecology as well [19, 27, 62, 78–80].
Products were denatured at 95° for 5 min and immediately transferred to ice for snap-cooling to produce single-strands and hinder reannealing. The ssDNA was then mixed with loading dye and loaded on a non-denaturing polyacrylamide gel (CleanGel HP-15 from ETC, Kirchentellinsfurth) according to the manufacturer's protocol and with the following adaptations: 12°C running temperature, pre-run for 20 min at 200 V, 20 mA, 10 W and main run for 4.5 h at 450 V, 30 mA, 20 W.
Gels were silver stained, scanned and processed with Photoshop 7.0 (Adobe Systems Inc, Beaverton) to align individual's band patterns. We counted bands with the same mobility as same alleles and each of these alleles was sequenced at least once from both directions, if possible twice from different individuals, to confirm this assumption. Appropriate bands were cut from the polyacrylamide gel, eluted in TBE buffer, amplified with the same protocol as above but only 30 cycles and using 1-2 μl of eluate. The PCR product was gel purified, ethanol-precipitated as above and then sequenced using the chain termination method  and either the JS1 forward or the GH50 reverse primer.
Sequences were aligned manually and obvious PCR misincorporations (unreproducible SNPs) were omitted. Alleles were only counted when they could be supported by two independent sequences from one individual or, if possible, from two different individuals. Primer sequences were excluded from analysis.
Genetic analysis and statistical treatment
Sequence alignment was done manually with the sequence analysis software BioEdit . We used GENECONV , a powerful tool for the detection of sequence recombination , to test for traces of gene conversion events within and between loci. A phylogenetic tree was built with MEGA 4  based on neighbour-joining algorithm and nucleotide distance. The distance was estimated in MEGA 4 with a maximum composite likelihood method  based on a model by Tamura & Nei . This model corrects for multiple hits, takes into account the differences in substitution rates between nucleotides and the inequality of nucleotide frequencies and also distinguishes between transitional and transversional substitution rates.
We used the software MEGA 4 to calculate dN/dS ratios and the therein implemented Z-test to test for signs of selection. This was done for the full sequence and separately for those sites only that are directly involved in antigen binding (ABS) . To confirm this selection estimate we also performed a model test with CODEML. This tool is included in the PAML4 package  and calculates maximum likelihood estimates for the fit of models with different assumptions of selection patterns to the sequence data. We tested the models M7 (neutral) against M8 (positive selection). The models are described in detail in Yang et al.  and can be compared using the likelihood ratio test by calculating the likelihood difference 2Δl = 2(l1-l0) and compare it to a X2-distribution with the degree of freedom equal to the difference in the number of estimated parameters .
As a measure of MHC diversity besides the pure number of allelic variants, we determined allele divergence, i.e. the average genetic distance between all MHC IIB alleles of an individual. This distance was calculated in MEGA 4 as the sum of all pairwise amino acid p-distances between all alleles of an individual, divided by the number of possible pairings (2 alleles = 1, 3 alleles = 3, 4 alleles = 6). The p-distance was calculated according to the model of Tamura & Nei . The genetic distance for individuals with only one allele is zero. The same estimate for allele divergence was also calculated taking only the antigen binding sites into account (ABS divergence). Correlation between allele number and average distance was tested with Spearman rank correlation because of non-normality of the data.
Body condition as a measure of physical fitness is difficult to estimate in a living animal [90, 91]. Therefore we used an approximation index, taking the residuals from the regression of body mass against body length . To test effects of study site, cohort (time of sampling), age and MHC diversity on body condition, we used linear models (ANCOVA) and the step function in the statistical software R, version 2.6 , which selects the most appropriate model based on Akaike's information criterion (AIC). Three different estimates for MHC diversity were used independently because of co-linearity: allele number, allele divergence and ABS divergence. Normality of response variable and independence of included factors was verified. Sex could only be determined in a subsample of the individuals and was therefore left out of the models. However, we tested the effect of sex on body condition and parasite load (species richness and infection intensity) within this subsample using t-test or Mann-Whitney U tests. For analysis of infection parameters we included only parasitized individuals (n = 43). Since it is unlikely to find individuals without a single parasite egg (lowest EPG count in parasitized individuals = 50), we feared that the egg count was confounded in those three individuals and might deteriorate the results.
Parasite species richness within each individual, expressed as number of different helminth morphotypes, was not normally distributed (Shapiro-Wilk test, p < 0.001). Therefore we used a generalized linear model (GLM, quasipoisson distribution, log link). We analysed the effect of study site, cohort, age and MHC diversity on parasite species richness. Again we tested allele number, allele divergence and ABS divergence independently. The model was optimised manually, because the step function is not available for models with quasipoisson distribution. The effect of body condition was tested in a separate model without MHC diversity because of co-linearity. The validity of the chosen models was evaluated by testing the distribution of the residuals.
To analyse the infection intensity in more detail we focussed on nematodes, because of their high prevalence. Effects of study site, cohort, age and MHC diversity on nematode infection intensity (EPG of all nematode morphotypes pooled) were tested with a negative binomial GLM with log link (glm.nb function in the MASS package, ), because of its distribution. Again the step function was applied to find the best fitting model. The effect of allele number, allele divergence and body condition on infection intensity was tested in independent models because of co-linearity. To analyse the effect of individual alleles, they were included singly in the GLM with the other factors. A Bonferroni correction had to be applied for testing multiple alleles. All tests were performed two-tailed using R 2.6  or SPSS 13 (SPSS Inc., Chicago) and test assumptions were verified. All tests are also assigned in the results section.
We thank the Economic Planning Unit Malaysia for granting a research permit and various local authorities and people for support in the field. We are particularly indebted to Maklarin B. Lakim and Alim Biun (Sabah Parks) as well as Matamin and Awang Kiddra (Kampong Monggis) for support in the field. We also would like to thank Christophe Eizaguirre, two anonymous referees and the editor for comments on earlier versions of the manuscript and Irene Tomaschewski for help in the lab.
- Hamilton WD, Axelrod R, Tanese R: Sexual reproduction as an adaptation to resist parasites (a review). Proc Natl Acad Sci USA. 1990, 87 (9): 3566-3573. 10.1073/pnas.87.9.3566.PubMed CentralView ArticlePubMedGoogle Scholar
- Ebert D, Hamilton WD: Sex against virulence: the coevolution of parasitic diseases. Trends Ecol Evol. 1996, 11 (2): 79-82. 10.1016/0169-5347(96)81047-0.View ArticlePubMedGoogle Scholar
- Apanius V, Penn D, Slev PR, Ruff LR, Potts WK: The nature of selection on the major histocompatibility complex. Crit Rev Immunol. 1997, 17 (2): 179-224.View ArticlePubMedGoogle Scholar
- Bernatchez L, Landry C: MHC studies in nonmodel vertebrates: what have we learned about natural selection in 15 years. J Evol Biol. 2003, 16 (3): 363-377. 10.1046/j.1420-9101.2003.00531.x.View ArticlePubMedGoogle Scholar
- Sommer S: The importance of immune gene variability (MHC) in evolutionary ecology and conservation. Front Zool. 2005, 2 (1): 16-10.1186/1742-9994-2-16.PubMed CentralView ArticlePubMedGoogle Scholar
- Milinski M: The major histocompatibility complex, sexual selection, and mate choice. Annu Rev Ecol Evol Syst. 2006, 37 (1): 159-186. 10.1146/annurev.ecolsys.37.091305.110242.View ArticleGoogle Scholar
- Klein J: Natural history of the major histocompatibility complex. 1986Google Scholar
- Piertney SB, Oliver MK: The evolutionary ecology of the major histocompatibility complex. Heredity. 2005, 96 (1): 7-21.Google Scholar
- Potts WK, Wakeland EK: Evolution of diversity at the major histocompatibility complex. Trends Ecol Evol. 1990, 5 (6): 181-187. 10.1016/0169-5347(90)90207-T.View ArticlePubMedGoogle Scholar
- Edwards SV, Hedrick PW: Evolution and ecology of MHC molecules: From genomics to sexual selection. Trends Ecol Evol. 1998, 13 (8): 305-311. 10.1016/S0169-5347(98)01416-5.View ArticlePubMedGoogle Scholar
- Summers K, McKeon S, Sellars J, Keusenkothen M, Morris J, Gloeckner D, Pressley C, Price B, Snow H: Parasitic exploitation as an engine of diversity. Biol Rev. 2003, 78: 639-675. 10.1017/S146479310300616X.View ArticlePubMedGoogle Scholar
- Clarke B, Kirby DRS: Maintenance of histocompatibility polymorphisms. Nature. 1966, 211 (5052): 999-1000. 10.1038/211999a0.View ArticlePubMedGoogle Scholar
- Bodmer WF: Evolutionary significance of the HL-A system. Nature. 1972, 237 (5351): 139-183. 10.1038/237139a0.View ArticlePubMedGoogle Scholar
- Takahata N, Nei M: Allelic genealogy under overdominant and frequency-dependent selection and polymorphism of major histocompatibility complex loci. Genetics. 1990, 124 (4): 967-978.PubMed CentralPubMedGoogle Scholar
- Van Valen L: A new evolutionary law. Evol Theor. 1973, 1: 1-30.Google Scholar
- Hill AV, Allsopp CE, Kwiatkowski D, Anstey NM, Twumasi P, Rowe PA, Bennett S, Brewster D, McMichael AJ, Greenwood BM: Common west African HLA antigens are associated with protection from severe malaria. Nature. 1991, 352 (6336): 595-600. 10.1038/352595a0.View ArticlePubMedGoogle Scholar
- Thursz MR, Kwiatkowski D, Allsopp CE, Greenwood BM, Thomas HC, Hill AV: Association between an MHC class II allele and clearance of hepatitis B virus in the Gambia. N Engl J Med. 1995, 332 (16): 1065-1069. 10.1056/NEJM199504203321604.View ArticlePubMedGoogle Scholar
- Harf R, Sommer S: Association between major histocompatibility complex class II DRB alleles and parasite load in the hairy-footed gerbil, Gerbillurus paeba, in the southern Kalahari. Mol Ecol. 2005, 14: 85-91. 10.1111/j.1365-294X.2004.02402.x.View ArticlePubMedGoogle Scholar
- Meyer-Lucht Y, Sommer S: MHC diversity and the association to nematode parasitism in the yellow-necked mouse (Apodemus flavicollis). Mol Ecol. 2005, 14 (7): 2233-2243. 10.1111/j.1365-294X.2005.02557.x.View ArticlePubMedGoogle Scholar
- Westerdahl H, Waldenstrom J, Hansson B, Hasselquist D, von Schantz T, Bensch S: Associations between malaria and MHC genes in a migratory songbird. Proc R Soc Lond B Biol Sci. 2005, 272 (1571): 1511-1518. 10.1098/rspb.2005.3113.View ArticleGoogle Scholar
- Young K, Frodsham A, Doumbo OK, Gupta S, Dolo A, Hu JT, Robson KJH, Crisanti A, Hill AVS, Gilbert SC: Inverse associations of human leukocyte antigen and malaria parasite types in two West African populations. Infect Immun. 2005, 73 (2): 953-955. 10.1128/IAI.73.2.953-955.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Schwensow N, Fietz J, Dausmann KH, Sommer S: Neutral versus adaptive genetic variation in parasite resistance: importance of major histocompatibility complex supertypes in a free-ranging primate. Heredity. 2007, 99: 265-277. 10.1038/sj.hdy.6800993.View ArticlePubMedGoogle Scholar
- Loiseau C, Zoorob R, Garnier S, Birard J, Federici P, Julliard R, Sorci G: Antagonistic effects of a Mhc class I allele on malaria-infected house sparrows. Ecol Lett. 2008, 11 (3): 258-265. 10.1111/j.1461-0248.2007.01141.x.View ArticlePubMedGoogle Scholar
- Doherty PC, Zinkernagel RM: Enhanced immunological surveillance in mice heterozygous at H-2 gene complex. Nature. 1975, 256 (5512): 50-52. 10.1038/256050a0.View ArticlePubMedGoogle Scholar
- Klein J, Figueroa F: Evolution of the Major Histocompatibility Complex. Crit Rev Immunol. 1986, 6 (4): 295-386.PubMedGoogle Scholar
- Penn DJ, Damjanovich K, Potts WK: MHC heterozygosity confers a selective advantage against multiple-strain infections. Proc Natl Acad Sci USA. 2002, 99 (17): 11260-11264. 10.1073/pnas.162006499.PubMed CentralView ArticlePubMedGoogle Scholar
- Wegner KM, Kalbe M, Kurtz J, Reusch TBH, Milinski M: Parasite selection for immunogenetic optimality. Science. 2003, 301 (5638): 1343-1343. 10.1126/science.1088293.View ArticlePubMedGoogle Scholar
- Wedekind C, Walker M, Little TJ: The course of malaria in mice: Major histocompatibility complex (MHC) effects, but no general MHC heterozygote advantage in single-strain infections. Genetics. 2005, 170 (3): 1427-1430. 10.1534/genetics.105.040683.PubMed CentralView ArticlePubMedGoogle Scholar
- Hughes AL, Nei M: Pattern of nucleotide substitution at major histocompatibility complex class I loci reveals overdominant selection. Nature. 1988, 335 (6186): 167-170. 10.1038/335167a0.View ArticlePubMedGoogle Scholar
- Hughes AL, Nei M: Nucleotide substitution at major histocompatibility complex class-II loci - evidence for overdominant selection. Proc Natl Acad Sci USA. 1989, 86 (3): 958-962. 10.1073/pnas.86.3.958.PubMed CentralView ArticlePubMedGoogle Scholar
- Lipsitch M, Bergstrom CT, Antia R: Effect of human leukocyte antigen heterozygosity on infectious disease outcome: the need for allele-specific measures. BMC Med Genet. 2003, 4: 2-10.1186/1471-2350-4-2.PubMed CentralView ArticlePubMedGoogle Scholar
- Ilmonen P, Penn DJ, Damjanovich K, Morrison L, Ghotbi L, Potts WK: Major histocompatibility complex heterozygosity reduces fitness in experimentally infected mice. Genetics. 2007, 176 (4): 2501-2508. 10.1534/genetics.107.074815.PubMed CentralView ArticlePubMedGoogle Scholar
- De Boer RJ, Borghans JA, van Boven M, Kesmir C, Weissing FJ: Heterozygote advantage fails to explain the high degree of polymorphism of the MHC. Immunogenetics. 2004, 55 (11): 725-731. 10.1007/s00251-003-0629-y.View ArticlePubMedGoogle Scholar
- Wakeland EK, Boehme S, She JX, Lu CC, McIndoe RA, Cheng I, Ye Y, Potts WK: Ancestral polymorphisms of MHC class-II genes - divergent allele advantage. Immunol Res. 1990, 9 (2): 115-122. 10.1007/BF02918202.View ArticlePubMedGoogle Scholar
- Gaston KJ: Global patterns in biodiversity. Nature. 2000, 405 (6783): 220-227. 10.1038/35012228.View ArticlePubMedGoogle Scholar
- Wells K, Pfeiffer M, Lakim MB, Kalko EKV: Movement trajectories and habitat partitioning of small mammals in logged and unlogged rain forests on Borneo. J Anim Ecol. 2006, 75 (5): 1212-1223. 10.1111/j.1365-2656.2006.01144.x.View ArticlePubMedGoogle Scholar
- Wells K: Diverstiy, dynamics and spacing patterns of a small mammal community in a primary rainforest in Sabah, Malaysia. Comparing an arboreal and terrestrial habitat. MSc-thesis. 2002, Würzburg: Julius-Maximilian-UniversityGoogle Scholar
- Lim BL: Distribution, relative abundance, food habits, and parasite patterns of giant rats (Rattus) in West Malaysia. J Mammal. 1970, 51 (4): 730-740. 10.2307/1378298.View ArticlePubMedGoogle Scholar
- Wells K, Kalko EKV, Lakim MB, Pfeiffer M: Effects of rain forest logging on species richness and assemblage composition of small mammals in Southeast Asia. J Biogeogr. 2007, 34 (6): 1087-1099. 10.1111/j.1365-2699.2006.01677.x.View ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucl Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Brown JH, Jardetzky TS, Gorga JC, Stern LJ, Urban RG, Strominger JL, Wiley DC: Three-dimensional structure of the human class II histocompatibility antigen HLA-DR1. Nature. 1993, 364 (6432): 33-39. 10.1038/364033a0.View ArticlePubMedGoogle Scholar
- Nielsen R: Statistical tests of selective neutrality in the age of genomics. Heredity. 2001, 86: 641-647. 10.1046/j.1365-2540.2001.00895.x.View ArticlePubMedGoogle Scholar
- Hughes AL: Looking for Darwin in all the wrong places: the misguided quest for positive selection at the nucleotide sequence level. Heredity. 2007, 99 (4): 364-373. 10.1038/sj.hdy.6801031.View ArticlePubMedGoogle Scholar
- Ensembl: Dog Genome release 52. [http://dec2008.archive.ensembl.org/Canis_familiaris/]
- Go Y, Satta Y, Kawamoto Y, Rakotoarisoa G, Randrianjafy A, Koyama N, Hirai H: Frequent segmental sequence exchanges and rapid gene duplication characterize the MHC class I genes in lemurs. Immunogenetics. 2003, 55 (7): 450-461. 10.1007/s00251-003-0613-6.View ArticlePubMedGoogle Scholar
- Reusch TB, Schaschl H, Wegner KM: Recent duplication and inter-locus gene conversion in major histocompatibility class II genes in a teleost, the three-spined stickleback. Immunogenetics. 2004, 56 (6): 427-437. 10.1007/s00251-004-0704-z.View ArticlePubMedGoogle Scholar
- Malaga-Trillo E, Zaleska-Rutczynska Z, McAndrew B, Vincek V, Figueroa F, Sultmann H, Klein J: Linkage relationships and haplotype polymorphism among cichlid Mhc class II B loci. Genetics. 1998, 149 (3): 1527-1537.PubMed CentralPubMedGoogle Scholar
- Bontrop RE, Otting N, de Groot NG, Doxiadis GGM: Major histocompatibility complex class II polymorphisms in primates. Immunol Rev. 1999, 167 (1): 339-350. 10.1111/j.1600-065X.1999.tb01403.x.View ArticlePubMedGoogle Scholar
- Horton R, Gibson R, Coggill P, Miretti M, Allcock RJ, Almeida J, Forbes S, Gilbert JGR, Halls K, Harrow JL, et al: Variation analysis and gene annotation of eight MHC haplotypes: The MHC haplotype project. Immunogenetics. 2008, 60 (1): 1-18. 10.1007/s00251-007-0262-2.PubMed CentralView ArticlePubMedGoogle Scholar
- Bowen L, Aldridge BM, Gulland F, Bonn W, DeLong R, Melin S, Lowenstine LJ, Stott JL, Johnson ML: Class II multiformity generated by variable MHC-DRB region configurations in the California sea lion (Zalophus californianus). Immunogenetics. 2004, 56 (1): 12-27. 10.1007/s00251-004-0655-4.View ArticlePubMedGoogle Scholar
- Doxiadis GGM, Otting N, de Groot NG, Noort R, Bontrop RE: Unprecedented polymorphism of MHC-DRB region configurations in rhesus macaques. J Immunol. 2000, 164 (6): 3193-3199.View ArticlePubMedGoogle Scholar
- Erlich HA, Bugawan TL: HLA DNA Typing. PCR Protocols. Edited by: Innis MA, Gelfand DH, Sninsky JJ, White TJ. 1990, New York: Academic Press, 261-271.Google Scholar
- Schad J, Sommer S, Ganzhorn JU: MHC variability of a small lemur in the littoral forest fragments of Southeastern Madagascar. Conserv Genet. 2004, 5 (3): 299-309. 10.1023/B:COGE.0000031137.50239.d3.View ArticleGoogle Scholar
- Froeschke G, Sommer S: MHC class II DRB variability and parasite load in the striped mouse (Rhabdomys pumilio) in the southern Kalahari. Mol Biol Evol. 2005, 22 (5): 1254-1259. 10.1093/molbev/msi112.View ArticlePubMedGoogle Scholar
- Cutrera AP, Lacey EA: Major histocompatibility complex variation in Talas Tuco-Tucos: The influence of demography on selection. J Mammal. 2006, 87 (4): 706-716. 10.1644/05-MAMM-A-383R1.1.View ArticleGoogle Scholar
- Axtner J, Sommer S: Gene duplication, allelic diversity, selection processes and adaptive value of MHC class II DRB genes of the bank vole, Clethrionomys glareolus. Immunogenetics. 2007, 59 (5): 417-426. 10.1007/s00251-007-0205-y.View ArticlePubMedGoogle Scholar
- Hurt P, Walter L, Sudbrak R, Klages S, Muller I, Shiina T, Inoko H, Lehrach H, Gunther E, Reinhardt R, et al: The genomic sequence and comparative analysis of the rat major histocompatibility complex. Genome Res. 2004, 14 (4): 631-639. 10.1101/gr.1987704.PubMed CentralView ArticlePubMedGoogle Scholar
- Goüy De Bellocq J, Charbonnel N, Morand S: Coevolutionary relationship between helminth diversity and MHC class II polymorphism in rodents. J Evol Biol. 2008, 21 (4): 1144-1150. 10.1111/j.1420-9101.2008.01538.x.View ArticlePubMedGoogle Scholar
- Richman AD, Herrera LG, Nash D: MHC class II beta sequence diversity in the deer mouse (Peromyscus maniculatus): Implications for models of balancing selection. Mol Ecol. 2001, 10 (12): 2765-2773.View ArticlePubMedGoogle Scholar
- Landry C, Garant D, Duchesne P, Bernatchez L: 'Good genes as heterozygosity': The major histocompatibility complex and mate choice in Atlantic salmon (Salmo salar). Proc R Soc Lond B Biol Sci. 2001, 268 (1473): 1279-1285. 10.1098/rspb.2001.1659.View ArticleGoogle Scholar
- Consuegra S, Garcia de Leaniz C: MHC-mediated mate choice increases parasite resistance in salmon. Proc R Soc Lond B Biol Sci. 2008, 275 (1641): 1397-1403. 10.1098/rspb.2008.0066.View ArticleGoogle Scholar
- Schwensow N, Eberle M, Sommer S: Compatibility counts: MHC-associated mate choice in a wild promiscuous primate. Proc R Soc Lond B Biol Sci. 2008, 275 (1634): 555-564. 10.1098/rspb.2007.1433.View ArticleGoogle Scholar
- Neff BD, Garner SR, Heath JW, Heath DD: The MHC and non-random mating in a captive population of Chinook salmon. Heredity. 2008, 101 (2): 175-185. 10.1038/hdy.2008.43.View ArticlePubMedGoogle Scholar
- Sette A, Sidney J: HLA supertypes and supermotifs: a functional perspective on HLA polymorphism. Curr Opin Immunol. 1998, 10 (4): 478-482. 10.1016/S0952-7915(98)80124-6.View ArticlePubMedGoogle Scholar
- Langefors A, Lohm J, Grahn M, Andersen O, von Schantz T: Association between major histocompatibility complex class IIB alleles and resistance to Aeromonas salmonicida in Atlantic salmon. Proc R Soc Lond B Biol Sci. 2001, 268 (1466): 479-485. 10.1098/rspb.2000.1378.View ArticleGoogle Scholar
- Bonneaud C, Perez-Tris J, Federici P, Chastel O, Sorci G: Major histocompatibility alleles associated with local resistance to malaria in a passerine. Evolution. 2006, 60 (2): 383-389.View ArticlePubMedGoogle Scholar
- Borghans JA, Beltman JB, De Boer RJ: MHC polymorphism under host-pathogen coevolution. Immunogenetics. 2004, 55 (11): 732-739. 10.1007/s00251-003-0630-5.View ArticlePubMedGoogle Scholar
- Wells K, Smales LR, Kalko EKV, Pfeiffer M: Impact of rain-forest logging on helminth assemblages in small mammals (Muridae, Tupaiidae) from Borneo. J Trop Ecol. 2007, 23: 35-43. 10.1017/S0266467406003804.View ArticleGoogle Scholar
- Gannon WL, Sikes RS, The animal care and use committee of the American Society of Mammalogists: Guidelines of the American Socitey of Mammalogists for the use of wild animals in research. J Mammal. 2007, 88 (3): 809-823. 10.1644/06-MAMM-F-185R1.1.View ArticleGoogle Scholar
- Thienpont D, Rochette F, Vanparijs O: Diagnosing helminthiasis by coprological examination. 1986, Beerse, Belgium: Janssen Research FoundationGoogle Scholar
- Moss R, Trenholm IB, Watson A, Parr R: Parasitism, predation and survival of hen red grouse Lagopus lagopus scoticus in spring. J Anim Ecol. 1990, 59 (2): 631-642. 10.2307/4885.View ArticleGoogle Scholar
- Skorping A, Read AF, Keymer AE: Life-history covariation in intestinal nematodes of mammals. Oikos. 1991, 60 (3): 365-372. 10.2307/3545079.View ArticleGoogle Scholar
- Seivwright LJ, Redpath SM, Mougeot F, Watt L, Hudson PJ: Faecal egg counts provide a reliable measure of Trichostrongylus tenuis intensities in free-living red grouse Lagopus lagopus scoticus. J Helminthol. 2004, 78 (1): 69-76. 10.1079/JOH2003220.View ArticlePubMedGoogle Scholar
- Gillespie T: Noninvasive assessment of gastrointestinal parasite infections in free-ranging primates. Int J Primatol. 2006, 27 (4): 1129-1143. 10.1007/s10764-006-9064-x.View ArticleGoogle Scholar
- Hayashi K: PCR-SSCP: A method for detection of mutations. Genet Anal: Tech Appl. 1992, 9 (3): 73-79.View ArticleGoogle Scholar
- Lessa EP, Applebaum G: Screening techniques for detecting allelic variation in DNA sequences. Mol Ecol. 1993, 2 (2): 119-129. 10.1111/j.1365-294X.1993.tb00006.x.View ArticlePubMedGoogle Scholar
- Sunnucks P, Wilson ACC, Beheregaray LB, Zenger K, French J, Taylor AC: SSCP is not so difficult: the application and utility of single-stranded conformation polymorphism in evolutionary biology and molecular ecology. Mol Ecol. 2000, 9 (11): 1699-1710. 10.1046/j.1365-294x.2000.01084.x.View ArticlePubMedGoogle Scholar
- Sommer S, Tichy H: Major histocompatibility complex (MHC) class II polymorphism and paternity in the monogamous Hypogeomys antimena, the endangered, largest endemic Malagasy rodent. Mol Ecol. 1999, 8: 1259-1272. 10.1046/j.1365-294X.1999.00687.x.View ArticlePubMedGoogle Scholar
- Reusch TB, Haberli MA, Aeschlimann PB, Milinski M: Female sticklebacks count alleles in a strategy of sexual selection explaining MHC polymorphism. Nature. 2001, 414 (6861): 300-302. 10.1038/35104547.View ArticlePubMedGoogle Scholar
- Simkova A, Ottova E, Morand S: MHC variability, life-traits and parasite diversity of European cyprinid fish. Evol Ecol. 2006, 20 (5): 465-477.View ArticleGoogle Scholar
- Sanger F, Nicklen S, Coulson AR: DNA sequencing with chain-terminating inhibitors. Proc Natl Acad Sci USA. 1977, 74 (12): 5463-5467. 10.1073/pnas.74.12.5463.PubMed CentralView ArticlePubMedGoogle Scholar
- Hall T: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999, 41: 95-98.Google Scholar
- Sawyer S: Statistical tests for detecting gene conversion. Mol Biol Evol. 1989, 6 (5): 526-538.PubMedGoogle Scholar
- Posada D: Evaluation of methods for detecting recombination from DNA sequences: Empirical data. Mol Biol Evol. 2002, 19 (5): 708-717.View ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) Software Version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
- Tamura K, Nei M, Kumar S: Prospects of inferring very large phylogenies by using the neighbor-joining method. Proc Natl Acad Sci USA. 2004, 101: 11030-11035. 10.1073/pnas.0404206101.PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Nei M: Estimation of the number of nucleotide substitutions in the control region of mitochondrial-DNA in humans and chimpanzees. Mol Biol Evol. 1993, 10 (3): 512-526.PubMedGoogle Scholar
- Yang Z: PAML 4: Phylogenetic Analysis by Maximum Likelihood. Mol Biol Evol. 2007, 24 (8): 1586-1591. 10.1093/molbev/msm088.View ArticlePubMedGoogle Scholar
- Yang ZH, Bielawski JP: Statistical methods for detecting molecular adaptation. Trends Ecol Evol. 2000, 15 (12): 496-503. 10.1016/S0169-5347(00)01994-7.View ArticlePubMedGoogle Scholar
- Schulte-Hostedde AI, Zinner B, Millar JS, Hickling GJ: Restitution of mass-size residuals: Validating body condition indices. Ecology. 2005, 86 (1): 155-163. 10.1890/04-0232.View ArticleGoogle Scholar
- Stevenson RD, Woods WA: Condition indices for conservation: new uses for evolving tools. Integr Comp Biol. 2006, 46 (6): 1169-1190. 10.1093/icb/icl052.View ArticlePubMedGoogle Scholar
- R Development Core Team: R: A language and environment for statistical computing. 2007, Vienna, Austria: R Foundation for Statistical ComputingGoogle Scholar
- Venables W, Ripley B: Modern Applied Statistics with S. 2002, New York: Springer, 4View ArticleGoogle 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.