MHC-mediated mate choice increases parasite resistance in salmon

Natural (parasite-driven) and sexual selection are thought to maintain high polymorphism in the genes of the major histocompatibility complex (MHC), but support for a link between mate choice, MHC variation and increased parasite resistance is circumstantial. We compared MHC diversity and Anisakis loads among anadromous Atlantic salmon (Salmo salar L.) returning to four rivers to spawn, which had originated from natural spawning (parents allowed to mate freely) or artificial crosses (parents deprived from the potential benefits of mate choice). We found that the offspring of artificially bred salmon had higher parasite loads and were almost four times more likely to be infected than free-mating salmon, despite having similar levels of MHC diversity. Moreover, the offspring of wild salmon were more MHC dissimilar than the offspring of artificially crossed salmon, and uninfected fish were more dissimilar for MHC than infected fish. Thus, our results suggest a link between disassortative mating and offspring benefits and indicate that MHC-mediated mate choice and natural (parasite-driven) selection act in combination to maintain MHC diversity, and hence fitness. Therefore, artificial breeding programmes that negate the potential genetic benefits of mate choice may result in inherently inferior offspring, regardless of population size, rearing conditions or genetic diversity.     

Genetic variation in MHC class II expression and interactions with MHC sequence polymorphism in three-spined sticklebacks

Genes of the major histocompatibility complex (MHC) have been studied for several decades because of their pronounced allelic polymorphism. Structural allelic polymorphism is, however, not the only source of variability subjected to natural selection. Genetic variation may also exist in gene expression patterns. Here, we show that in a natural population of three-spined sticklebacks (Gasterosteus aculeatus) the expression of MHC class IIB genes was positively correlated with parasite load, which indicates increased immune activation of the MHC when infections are frequent. To experimentally study MHC expression, we used laboratory-bred sticklebacks that were exposed to three naturally occurring species of parasite. We found strong differences in MHC class IIB expression patterns among fish families, which were consistent over two generations, thus demonstrating a genetic component. The average number of MHC class IIB sequence variants within families was negatively correlated to the MHC expression level suggesting compensatory up-regulation in fish with a low (i.e. suboptimal) MHC sequence variability. The observed differences among families and the negative correlation with individual sequence diversity imply that MHC expression is evolutionary relevant for the onset and control of the immune response in natural populations.     

TEMPORAL VARIATION AT THE MHC CLASS IIB IN WILD POPULATIONS OF THE GUPPY (POECILIA RETICULATA)

Understanding genetic diversity in natural populations is a fundamental objective of evolutionary biology. The immune genes of the major histocompatibility complex (MHC) are excellent candidates to study such diversity because they are highly polymorphic in populations. Although balancing selection may be responsible for maintaining diversity at these functionally important loci, temporal variation in selection pressure has rarely been examined. We examine temporal variation in MHC class IIB diversity in nine guppy (Poecilia reticulata) populations over two years. We found that five of the populations changed significantly more at the MHC than at neutral (microsatellite) loci as measured by F_ST, which suggests that the change at the MHC was due to selection and not neutral processes. Additionally, pairwise population differentiation measures at the MHC were higher in 2007 than in 2006, with the signature of selection changing from homogenizing to diversifying selection or neutral evolution. Interestingly, within the populations the magnitude of the change at the MHC between years was related to the change in the proportion of individuals infected by a common parasite, indicating a link between genetic structure and the parasite. Our data thereby implicate temporal variation in selective pressure as an important mechanism maintaining diversity at the MHC in wild populations.     

Balancing selection, random genetic drift, and genetic variation at the major histocompatibility complex in two wild populations of guppies (Poecilia reticulata)

Our understanding of the evolution of genes of the major histocompatibility complex (MHC) is rapidly increasing, but there are still enigmatic questions remaining, particularly regarding the maintenance of high levels of MHC polymorphisms in small, isolated populations. Here, we analyze the genetic variation at eight microsatellite loci and sequence variation at exon 2 of the MHC class IIB (DAB) genes in two wild populations of the Trinidadian guppy, Poecilia reticulata. We compare the genetic variation of a small (Ne, 100) and relatively isolated upland population to that of its much larger (Ne approximately 2400) downstream counterpart. As predicted, microsatellite diversity in the upland population is significantly lower and highly differentiated from the population further downstream. Surprisingly, however, these guppy populations are not differentiated by MHC genetic variation and show very similar levels of allelic richness. Computer simulations indicate that the observed level of genetic variation can be maintained with overdominant selection acting at three DAB loci. The selection coefficients differ dramatically between the upland (s > or = 0.2) and lowland (s < or = 0.01) populations. Parasitological analysis on wild-caught fish shows that parasite load is significantly higher on upland than on lowland fish, which suggests that large differences in selection intensity may indeed exist between populations. Based on the infection intensity, a substantial proportion of the upland fish would have suffered direct or indirect fitness consequences as a result of their high parasite loads. Selection by parasites plays a particularly important role in the evolution of guppies in the upland habitat, which has resulted in high levels of MHC diversity being maintained in this population despite considerable genetic drift.     

Individual MHC class I and MHC class IIB diversities are associated with male and female reproductive traits in the three-spined stickleback

Genes of the major histocompatibility complex (MHC) are indispensable for pathogen defence in vertebrates. With wild-caught three-spined sticklebacks (Gasterosteus aculeatus) we conducted the first study to relate individual reproductive parameters to both MHC class I and II diversities. An optimal MHC class IIB diversity was found for male nest quality. However, male breeding colouration was most intense at a maximal MHC class I diversity. One MHC class I allele was associated with male redness. Similarly, one MHC class IIB allele was associated with continuous rather than early female reproduction, possibly extending the reproductive period. Both alleles occurred more frequently with increasing individual allele diversity. We suggest that if an allele is currently not part of the optimum, it had not been propagated by choosy females. The parasite against which this allele provides resistance is therefore unlikely to have been predominant the previous year - a step to negative frequency-dependent selection.     

MHC variability, life-traits and parasite diversity of European cyprinid fish

Potential links between fish life history traits, an immune investment as measured by variability of the MHC genes and parasitism were analysed in 14 species of cyprinid fish. The hypothesis of the diversity of MHC genes being driven by high parasite diversity, i.e. species richness, was tested and a potential relationship between the MHC diversity and fish life-history traits including adult mortality rate, fecundity, longevity and maturity was investigated. Molecular techniques (SSCP and sequencing) were applied to analyse the MHC nucleotide diversity of the exons 2 and 3 of DAB genes belonging to the MHC class IIB. The comparative analyses, using phylogenetic independent contrasts, revealed a negative relationship between parasite species richness and adult fish mortality rate. We also found a positive relationship between nucleotide diversity of the exon 2 and parasite species richness. Our results suggest that fish species, of which populations are exposed to high parasite pressure, in terms of high parasite species richness, maintain a high genetic diversity of the exon 2 of the MHC genes (presenting the peptide binding regions), allowing them to decrease their natural mortality rate.     

Disentangling the role of MHC‐dependent ‘good genes' and ‘compatible genes' in mate‐choice decisions of three‐spined sticklebacks Gasterosteus aculeatus under semi‐natural conditions

To investigate and disentangle the role of major histocompatibility complex (MHC)‐based ‘good genes' and ‘compatible genes' in mate choice, three‐spined sticklebacks Gasterosteus aculeatus with specific MHC IIB genotypes were allowed to reproduce in an outdoor enclosure system. Here, fish were protected from predators but encountered their natural parasites. Mate choice for an intermediate genetic distance between parental MHC genotypes was observed, which would result in intermediate diversity in the offspring, but no mate choice based on good genes was found under the current semi‐natural conditions. Investigation of immunological variables revealed that the less‐specific innate immune system was more active in individuals with a genetically more divergent MHC allele repertoire. This suggests the need to compensate for an MHC‐diminished T‐cell repertoire and potentially explains the observed mate choice for intermediate MHC genetic distance. The present findings support a general pattern of mate choice for intermediate MHC diversity (i.e. compatible genes). In addition, the potentially dynamic role of MHC good genes in mate choice under different parasite pressures is discussed in the light of present and previous results.     

MHC class IIB alleles contribute to both additive and nonadditive genetic effects on survival in Chinook salmon

The genes of the major histocompatibility complex (MHC) are found in all vertebrates and are an important component of individual fitness through their role in disease and pathogen resistance. These genes are among the most polymorphic in genomes and the mechanism that maintains the diversity has been actively debated with arguments for natural selection centering on either additive or nonadditive genetic effects. Here, we use a quantitative genetics breeding design to examine the genetic effects of MHC class IIB alleles on offspring survivorship in Chinook salmon (Oncorhynchus tshawytscha). We develop a novel genetic algorithm that can be used to assign values to specific alleles or genotypes. We use this genetic algorithm to show simultaneous additive and nonadditive effects of specific MHC class IIB alleles and genotypes on offspring survivorship. The additive effect supports the rare-allele hypothesis as a potential mechanism for maintaining genetic diversity at the MHC. However, contrary to the overdominance hypothesis, the nonadditive effect led to underdominance at one heterozygous genotype, which could instead reduce variability at the MHC. Our algorithm is an advancement over traditional animal models that only partition variance in fitness to additive and nonadditive genetic effects, but do not allocate these effects to specific alleles and genotypes. Additionally, we found evidence of nonrandom segregation during meiosis in females that promotes an MHC allele that is associated with higher survivorship. Such nonrandom segregation could further reduce variability at the MHC and may explain why Chinook salmon has one of the lowest levels of MHC diversity of all vertebrates.     

Historical and contemporary selection of teleost MHC genes: did we leave the past behind?

The extreme polymorphism of antigen‐presenting genes of the major histocompatibility complex (MHC) has spurred intense research unparalleled for any other gene family. This applies also to teleosts where sequence information is available for 3559 MHC class I and class II allelic variants from 137 species. This review summarizes current knowledge on the origin and maintenance of diversity at classical MHC loci. Most studies identified positive selection (i.e. elevated rates of non‐synonymous over synonymous substitutions, dN/dS) as a sign of balancing selection. A meta‐analysis on nine species with sufficient numbers of class I and class II sequences revealed that recombination rate and intensity of positive selection were positively correlated, suggesting that recombination and gene conversion played a significant role in shaping the allelic repertoire. Processes that create diversity over long timescales need to be complemented by contemporary balancing selection, either through overdominance or frequency‐dependent selection, in order to explain the high allelic diversity observed today. While some evidence for overdominance exists for a few taxa (mainly salmonids) by correlating parasite infection data or survival to MHC genotypes, field or experimental data on negative frequency‐dependent selection are lacking altogether, even though some fish species are particularly suitable as model systems. Theoretical predictions suggest that negative frequency‐dependent selection is necessary to maintain the existing polymorphism. Hence, future empirical studies should focus on detecting signals that differentiate between mechanisms of contemporary selection rather than repeatedly showing historical selection events.     

Characterization of the major histocompatibility complex class II genes in miiuy croaker

Major histocompatibility complex (MHC) has a central role in the adaptive immune system by presenting foreign peptide to the T-cell receptor. In order to study the molecular function and genomic characteristic of class II genes in teleost, the full lengths of MHC class IIA and IIB cDNA and genomic sequence were cloned from miiuy croaker (Miichthys miiuy). As in other teleost, four exons and three introns were identified in miiuy croaker class IIA gene; but the difference is that six exons and five introns were identified in the miiuy croaker class IIB gene. The deduced amino acid sequence of class IIA and class IIB had 26.3-85.7% and 11.0-88.8% identity with those of mammal and teleost, respectively. Real-time quantitative RT-PCR demonstrated that the MHC class IIA and IIB were ubiquitously expressed in ten normal tissues; expression levels of MHC genes were found first upregulated and then downregulated, and finally by a recovery to normal level throughout the pathogenic bacteria infection process. In addition, we report on the underlying mechanism that maintains sequences diversity among many fish species. A series of site-model tests implemented in the CODEML program revealed that positive Darwinian selection is likely the cause of the molecular evolution in the fish MHC class II genes.     

Parasite diversity,patterns of MHC II variation and olfactory based mate choice in diverging three-spined stickleback ecotypes

Ecological speciation has been the subject of intense research in evolutionary biology but the genetic basis of the actual mechanism driving reproductive isolation has rarely been identified. The extreme polymorphism of the major histocompatibility complex (MHC), probably maintained by parasite-mediated selection, has been proposed as a potential driver of population divergence. We performed an integrative field and experimental study using three-spined stickleback river and lake ecotypes. We characterized their parasite load and variation at MHC class II loci. Fish from lakes and rivers harbor contrasting parasite communities and populations possess different MHC allele pools that could be the result of a combined action of genetic drift and parasite-mediated selection. We show that individual MHC class II diversity varies among populations and is lower in river ecotypes. Our results suggest the action of homogenizing selection within habitat type and diverging selection between habitat types. Finally, reproductive isolation was suggested by experimental evidence: in a flow channel design females preferred assortatively the odor of their sympatric male. This demonstrates the role of olfactory cues in maintaining reproductive isolation between diverging fish ecotypes.     

The role of major histocompatibility complex diversity in vigour of fish males (Abramis brama L.) and parasite selection

The major histocompatibility complex (MHC) presents a group of genes with highly polymorphic loci involved in specific immune responses. The factors maintaining extensive MHC polymorphism have been questioned, considering three possible hypotheses of parasite‐mediated selection driving an extensive MHC diversity (i.e. heterozygote advantage, rare‐allele advantage, and favouring optimal MHC diversity). The patterns of MHC diversity of class IIB genes were investigated following two noncontradicting hypotheses, parasite‐driven selection and MHC‐based mating preferences, using males of common bream collected in the spawning period. Two allelic groups DAB1 and DAB3 were recognized from the phylogenetic analyses. Individuals expressed one or two alleles of the same or different allelic groups. Several individuals shared identical alleles; however, the presence of parasite species was not associated with the occurrence of a particular allele. The presence of different allelic groups (only DAB1, only DAB3, or both DAB1 and DAB3) in individuals was not associated with parasite presence or diversity. The expression of two DAB1 alleles was associated with higher endoparasite abundance. Moreover, nucleotide diversity in individuals expressing a single type of alleles (DAB1 or DAB3) increased with the abundance of ectoparasitic Dactylogyrus spp. (Monogenea) and Ergasilus sp. (Crustacea). This suggests that the expression of two alleles of a single allelic type is related to high metazoan parasite infection whereas no significant influence of parasitism on the combined allelic form (the presence of both DAB1 and DAB3 alleles) was found. Moreover, the expression of two alleles of a single allelic type was related to decreased immunocompetence measured by spleen size. The condition factor was higher in fish expressing the combined allelic type. Thus, the presence of alleles of different lineages in individuals appears to be advantageous for individual male fitness. The expression of a single allelic type was related to higher sexual ornamentation, which could support the role of MHC in the hypothesis of the sexual selection of ‘good genes’. © 2007 The Linnean Society of London, Biological Journal of the Linnean Society, 2007, 90 , 525–538.     

Selection from parasites favours immunogenetic diversity but not divergence among locally adapted host populations

The unprecedented polymorphism in the major histocompatibility complex (MHC) genes is thought to be maintained by balancing selection from parasites. However, do parasites also drive divergence at MHC loci between host populations, or do the effects of balancing selection maintain similarities among populations? We examined MHC variation in populations of the livebearing fish Poecilia mexicana and characterized their parasite communities. Poecilia mexicana populations in the Cueva del Azufre system are locally adapted to darkness and the presence of toxic hydrogen sulphide, representing highly divergent ecotypes or incipient species. Parasite communities differed significantly across populations, and populations with higher parasite loads had higher levels of diversity at class II MHC genes. However, despite different parasite communities, marked divergence in adaptive traits and in neutral genetic markers, we found MHC alleles to be remarkably similar among host populations. Our findings indicate that balancing selection from parasites maintains immunogenetic diversity of hosts, but this process does not promote MHC divergence in this system. On the contrary, we suggest that balancing selection on immunogenetic loci may outweigh divergent selection causing divergence, thereby hindering host divergence and speciation. Our findings support the hypothesis that balancing selection maintains MHC similarities among lineages during and after speciation (trans‐species evolution).     

Comparative analysis of population structure across environments and geographical scales at major histocompatibility complex and microsatellite loci in Atlantic salmon (Salmo salar)

Evidence of selection acting on major histocompatibility complex (MHC) genes has been illustrated with the analysis of their nucleotide sequences and allele frequency distribution. Comparing the patterns of population differentiation at neutral markers and MHC genes in the wild may provide further insights about the relative role of selection and neutrality in shaping their diversity. In this study, we combine both methods to assess the role of selection on a MHC gene in Atlantic salmon. We compare variation at a MHC class II B locus and microsatellites among 14 samples from seven different rivers and seven subpopulations within a single river system covering a variety of habitats and different geographical scales. We show that diversifying selection is acting on the sites involved in antigen presentation and that balancing selection maintains a high level of polymorphism within populations. Despite important differences in habitat type, the comparison of the population structure at MHC and microsatellites on large geographical scales reveals a correlation between patterns of differentiation, indicating that drift and migration have been more important than selection in shaping population differentiation at the MHC locus. In contrast, strong discrepancies between patterns of population differentiation at the two types of markers provides support for the role of selection in shaping population structure within rivers. Together, these results confirm that natural selection is influencing MHC gene diversity in wild Atlantic salmon although neutral forces may also be important in their evolution.     

MHC-based mate choice combines good genes and maintenance of MHC polymorphism

Polymorphic genes of the major histocompatibility complex (MHC) are regarded as essential genes for individual fitness under conditions of natural and sexual selection. To test this hypothesis, we investigated the ultimate individual fitness trait — that of reproductive success. We used three-spined sticklebacks ( Gasterosteus aculeatus ) in seminatural enclosures, located in natural breeding areas where the experimental fish had been caught. During their reproductive period, fish were exposed continuously to their natural sympatric parasites. By genotyping almost 4000 eggs with nine microsatellites, we determined parenthood and inferred female mating decision. We found that with reference to their own MHC profile, female sticklebacks preferred to mate with males sharing an intermediate MHC diversity. In addition, males with a specific MHC haplotype were bigger and better at fighting a common parasite ( Gyrodactylus sp.). This translated directly into Darwinian fitness since fish harbouring this specific MHC haplotype were more likely to be chosen and had a higher reproductive output. We conclude that females also based their mating decision on a specific MHC haplotype conferring resistance against a common parasite. This identifies and supports 'good genes'. We argue that such an interaction between host and parasite driving assortative mating is not only a prerequisite for negative frequency-dependent selection — a potential mechanism to explain the maintenance of MHC polymorphism, but also potentially speciation.     

Variation in MHC genotypes in two populations of house sparrow (Passer domesticus) with different population histories

Small populations are likely to have a low genetic ability for disease resistance due to loss of genetic variation through inbreeding and genetic drift. In vertebrates, the highest genetic diversity of the immune system is located at genes within the major histocompatibility complex (MHC). Interestingly, parasite‐mediated selection is thought to potentially maintain variation at MHC loci even in populations that are monomorphic at other loci. Therefore, general loss of genetic variation in the genome may not necessarily be associated with low variation at MHC loci. We evaluated inter‐ and intrapopulation variation in MHC genotypes between an inbred (Aldra) and a relatively outbred population (Hestmannøy) of house sparrows (Passer domesticus) in a metapopulation at Helgeland, Norway. Genomic (gDNA) and transcribed (cDNA) alleles of functional MHC class I and IIB loci, along with neutral noncoding microsatellite markers, were analyzed to obtain relevant estimates of genetic variation. We found lower allelic richness in microsatellites in the inbred population, but high genetic variation in MHC class I and IIB loci in both populations. This suggests that also the inbred population could be under balancing selection to maintain genetic variation for pathogen resistance.     

Major histocompatibility complex diversity influences parasite resistance and innate immunity in sticklebacks

Proteins of the major histocompatibility complex (MHC) play a central role in the presentation of antigens to the adaptive immune system. The MHC also influences the odour-based choice of mates in humans and several animal taxa. It has recently been shown that female three-spined sticklebacks (Gasterosteus aculeatus) aim at a moderately high MHC diversity in their offspring when choosing a mate. Do they optimize the immune systems of their offspring? Using three-spined sticklebacks that varied in their individual numbers of MHC class IIB molecules, we tested, experimentally, whether allelic diversity at the MHC influences parasite resistance and immune parameters. We found that sticklebacks with low MHC diversity suffered more from parasite infection after experimental exposure to Schistocephalus solidus tapeworms and Glugea anomala microsporidians. They also showed the highest proportion of granulocytes and the strongest respiratory burst reaction, which are correlates of innate immunity. This indicates a strong activity of the innate immune system after challenge by parasites when MHC diversity is suboptimal. Individuals with very high allelic diversity at the MHC seemed inferior to those with moderately high diversity. Such a pattern is consistent with theoretical expectations of an optimal balance between the number of recognizable antigens and self-tolerance.     

Diversifying selection on MHC class I in the house sparrow (Passer domesticus)

Genes of the major histocompatibility complex (MHC) are the most polymorphic loci known in vertebrates. Two main hypotheses have been put forward to explain the maintenance of MHC diversity: pathogen-mediated selection and MHC-based mate choice. Host–parasite interactions can maintain MHC diversity via frequency-dependent selection, heterozygote advantage, and diversifying selection (spatially and/or temporally heterogeneous selection). In this study, we wished to investigate the nature of selection acting on the MHC class I across spatially structured populations of house sparrows ( Passer domesticus ) in France. To infer the nature of the selection, we compared patterns of population differentiation based on two types of molecular markers: MHC class I and microsatellites. This allowed us to test whether the observed differentiation at MHC genes merely reflects demographic and/or stochastic processes. At the global scale, diversifying selection seems to be the main factor maintaining MHC diversity in the house sparrow. We found that (i) overall population differentiation at MHC was stronger than for microsatellites, (ii) MHC marker showed significant isolation by distance. In addition, the slope of the regression of F _ST on geographical distance was significantly steeper for MHC than for microsatellites due to a stronger pairwise differentiation between populations located at large geographical distances. These results are in agreement with the hypothesis that spatially heterogeneous selective pressures maintain different MHC alleles at local scales, possibly resulting in local adaptation.     

EVALUATING PATTERNS OF CONVERGENT EVOLUTION AND TRANS‐SPECIES POLYMORPHISM AT MHC IMMUNOGENES IN TWO SYMPATRIC STICKLEBACK SPECIES

The immunologically important major histocompatibility complex (MHC) harbors some of the most polymorphic genes in vertebrates. These genes presumably evolve under parasite‐mediated selection and frequently show inconsistent allelic genealogies, where some alleles are more similar between species than within species. This phenomenon is thought to arise either from convergent evolution under parallel selection or from the preservation of ancient allelic lineages beyond speciation events (trans‐species polymorphism, TSP). Here, we examine natural populations of two sympatric stickleback species (Gasterosteus aculeatus and Pungitius pungitius) to investigate the contribution of these two mechanisms to the evolution of inconsistent allelic genealogies at the MHC. Overlapping parasite taxa between the two host species in three different habitats suggest contemporary parallel selection on the MHC genes. Accordingly, we detected a lack of species‐specific phylogenetic clustering in the immunologically relevant antigen‐binding residues of the MHC IIB genes which contrasted with the rest of the coding and noncoding sequence. However, clustering was not habitat‐specific and a codon‐usage analysis revealed patterns of similarity by descent. In this light, common descent via TSP, in combination with intraspecies gene conversion, rather than convergent evolution is the more strongly supported scenario for the inconsistent genealogy at the MHC.     

Neutral versus adaptive genetic variation in parasite resistance: importance of major histocompatibility complex supertypes in a free-ranging primate

Current discussions in evolutionary ecology and conservation genetics focus on the relative importance of using selective neutral markers or markers of coding genes to identify adaptive and evolutionary relevant processes. Genetic diversity might be particularly important in immune genes (e.g., in genes of the major histocompatibility complex, MHC), which are influencing pathogen and parasite resistance. We investigated the effects of neutral versus adaptive genetic variation in parasite resistance in a natural population of fat-tailed dwarf lemurs (Cheirogaleus medius). No association between neutral overall individual genetic diversity and parasite load could be detected. In 149 individuals, we identified 50 MHC class II alleles of the functionally important duplicated DRB locus. The investigation of the functional importance of immune gene (MHC) diversity and parasite selection in natural populations is often problematic due to extensive polymorphism in the MHC genes and restrictions in available sample sizes. Here, for the first time we applied an approach that has been developed in human medical studies. Eleven MHC class II supertypes were identified based on shared antigen-binding similarities. The number of individual MHC supertypes had no influence on the nematode burden. However, we found evidence for a specific MHC supertype (supertype 1) that was linked to infected individuals, a higher number of different nematode infections and high intensity of infection per individual. Moreover, one rare MHC supertype (supertype 7) was revealed to be advantageous with respect to parasite burden. Thus, our results add evidence to the small body of studies that show significant associations between specific MHC constitutions and naturally occurring parasites in the complexity of natural populations.