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

Population-level immunologic variation in wild threespine stickleback (Gasterosteusaculeatus)

Wild organisms are regularly exposed to a wide range of parasites, requiring the management of an effective immune response while avoiding immunopathology. Currently, our knowledge of immunoparasitology primarily derives from controlled laboratory studies, neglecting the genetic and environmental diversity that contribute to immune phenotypes observed in wild populations. To gain insight into the immunologic variability in natural settings, we examined differences in immune gene expression of two Alaskan stickleback (Gasterosteus aculeatus) populations with varying susceptibility to infection by the cestode Schistocephalus solidus. Between these two populations, we found distinct immune gene expression patterns at the population level in response to infection with fish from the high-infection population displaying signs of parasite-driven immune manipulation. Further, we found significant differences in baseline immune gene profiles between the populations, with uninfected low-infection population fish showing signatures of inflammation compared to uninfected high-infection population fish. These results shed light on divergent responses of wild populations to the same parasite, providing valuable insights into host-parasite interactions in natural ecosystems.

Polymorphism of MHC class IIB in an acheilognathid species, Rhodeus sinensis shaped by historical selection and recombination

BACKGROUND

Rhodeus sinensis is a bitterling species occurring throughout the numerous freshwater systems on the East Asia. Here, we analyzed the diversity of the MHC class IIB (DAB) genes from this species, which may offer meaningful insights into evolutionary processes in this species as well as other bitterlings.

RESULTS

Using cDNA and gDNA samples from 50 individuals, we discovered classical 140 allelic sequences that could be allocated into either DAB1 (Rhsi-DAB1) or DAB3 (Rhsi-DAB3). DAB sequences completely lacking the intron, but identical or similar to Rhsi-DAB1, were also discovered from our gDNA samples, and this intron loss likely originated from the retrotransposition events of processed mDNA. The β1 domain was the most polymorphic in both Rhsi-DAB1 and -DAB3. Putative peptide biding residues (PBRs) in Rhsi-DAB1, but not in Rhsi-DAB3, exhibited a significant dN/dS, presumably indicating that different selection pressures have acted on those two DABs. Recombination between different alleles seemed to have contributed to the increase of diversity in Rhsi-DABs. Upon phylogenetic analysis, Rhsi-DAB1 and -DAB3 formed independent clusters. Several alleles from other species of Cypriniformes were embedded in the clade of Rhsi-DAB1, whereas Rhsi-DAB3 clustered with alleles from the wider range of taxa (Cyprinodontiformes), indicating that these two Rhsi-DABs have taken different historical paths.

CONCLUSIONS

A great deal of MHC class IIB allelic diversity was found in R. sinensis, and gene duplication, selection and recombination may have contributed to this diversity. Based on our data, it is presumed that such historical processes have commonly or differently acted on the polymorphism of Rhsi-DAB1 and -DAB3.

Molecular and Structural Characterization of MHC Class II β Genes Reveals High Diversity in the Cold-Adapted Icefish Chionodraco hamatus

This study reports the presence of two distinct MHC class II β genes in the Antarctic icefish Chionodraco hamatus, belonging to the classical (ChhaDAB) and nonclassical (ChhaDBB) evolutionary lineages. By the application of targeted sequencing approach, a remarkable molecular diversity in the exon 2 sequence of the highly expressed gene ChhaDAB has been observed, resulting in an estimate of 92 different variants translated in 87 different peptides from 54 analysed icefish individuals. A highly conservative estimate, based on a 95% sequence identity threshold clustering, translate this variability in 41 different peptide clusters belonging to four different clades and showing the signature of different kinds of selection. In stark contrast, the poorly expressed ChhaDBB gene displayed a very low level of molecular diversity within exon 2, in agreement with expectations for a nonclassical MHC class II β gene.

Using Gene Transcription to Assess Ecological and Anthropological Stressors in Brown Bears

Increasingly, population- and ecosystem-level health assessments are performed using sophisticated molecular tools. Advances in molecular technology enable the identification of synergistic effects of multiple stressors on the individual physiology of different species. Brown bears (Ursus arctos) are an apex predator; thus, they are ideal candidates for detecting potentially ecosystem-level systemic perturbations using molecular-based tools. We used gene transcription to analyze 130 brown bear samples from three National Parks and Preserves in Alaska. Although the populations we studied are apparently stable in abundance and exist within protected and intact environments, differences in transcript profiles were noted. The most prevalent differences were among locations. The transcript patterns among groups reflect the influence of environmental factors, such as nutritional status, disease, and xenobiotic exposure. However, these profiles also likely represent baselines for each unique environment by which future measures can be made to identify early indication of population-level changes due to, for example, increasing Arctic temperatures. Some of those environmental changes are predicted to be potentially positive for brown bears, but other effects such as the manifestation of disease or indirect effects of oceanic acidification may produce negative impacts.

MHC class II DRB diversity predicts antigen recognition and is associated with disease severity in California sea lions naturally infected with Leptospira interrogans

We examined the associations between California sea lion MHC class II DRB (Zaca-DRB) configuration and diversity, and leptospirosis. As Zaca-DRB gene sequences are involved with antigen presentation of bacteria and other extracellular pathogens, we predicted that they would play a role in determining responses to these pathogenic spirochaetes. Specifically, we investigated whether Zaca-DRB diversity (number of genes) and configuration (presence of specific genes) explained differences in disease severity, and whether higher levels of Zaca-DRB diversity predicted the number of specific Leptospira interrogans serovars that a sea lion's serum would react against. We found that serum from diseased sea lions with more Zaca-DRB loci reacted against a wider array of serovars. Specific Zaca-DRB loci were linked to reactions with particular serovars. Interestingly, sea lions with clinical manifestation of leptospirosis that had higher numbers of Zaca-DRB loci were less likely to recover from disease than those with lower diversity, and those that harboured Zaca-DRB.C or -G were 4.5 to 5.3 times more likely to die from leptospirosis, regardless of the infective serovars. We propose that for leptospirosis, a disadvantage of having a wider range of antigen presentation might be increased disease severity due to immunopathology. Ours is the first study to examine the importance of Zaca-DRB diversity for antigen detection and disease severity following natural exposure to infective leptospires.

Gene Expression Contributes to the Recent Evolution of Host Resistance in a Model Host Parasite System

Heritable population differences in immune gene expression following infection can reveal mechanisms of host immune evolution. We compared gene expression in infected and uninfected threespine stickleback (Gasterosteus aculeatus) from two natural populations that differ in resistance to a native cestode parasite, Schistocephalus solidus. Genes in both the innate and adaptive immune system were differentially expressed as a function of host population, infection status, and their interaction. These genes were enriched for loci controlling immune functions known to differ between host populations or in response to infection. Coexpression network analysis identified two distinct processes contributing to resistance: parasite survival and suppression of growth. Comparing networks between populations showed resistant fish have a dynamic expression profile while susceptible fish are static. In summary, recent evolutionary divergence between two vertebrate populations has generated population-specific gene expression responses to parasite infection, affecting parasite establishment and growth.

Hook, Line and Infection: A Guide to Culturing Parasites, Establishing Infections and Assessing Immune Responses in the Three-Spined Stickleback

The three-spined stickleback (Gasterosteus aculeatus) is a model organism with an extremely well-characterized ecology, evolutionary history, behavioural repertoire and parasitology that is coupled with published genomic data. These small temperate zone fish therefore provide an ideal experimental system to study common diseases of coldwater fish, including those of aquacultural importance. However, detailed information on the culture of stickleback parasites, the establishment and maintenance of infections and the quantification of host responses is scattered between primary and grey literature resources, some of which is not readily accessible. Our aim is to lay out a framework of techniques based on our experience to inform new and established laboratories about culture techniques and recent advances in the field. Here, essential knowledge on the biology, capture and laboratory maintenance of sticklebacks, and their commonly studied parasites is drawn together, highlighting recent advances in our understanding of the associated immune responses. In compiling this guide on the maintenance of sticklebacks and a range of common, taxonomically diverse parasites in the laboratory, we aim to engage a broader interdisciplinary community to consider this highly tractable model when addressing pressing questions in evolution, infection and aquaculture.

No evidence of local adaptation of immune responses to Gyrodactylus in three-spined stickleback (Gasterosteus aculeatus)

Parasitism represents one of the most widespread lifestyles in the animal kingdom, with the potential to drive coevolutionary dynamics with their host population. Where hosts and parasites evolve together, we may find local adaptation. As one of the main host defences against infection, there is the potential for the immune response to be adapted to local parasites. In this study, we used the three-spined stickleback and its Gyrodactylus parasites to examine the extent of local adaptation of parasite infection dynamics and the immune response to infection. We took two geographically isolated host populations infected with two distinct Gyrodactylus species and performed a reciprocal cross-infection experiment in controlled laboratory conditions. Parasite burdens were monitored over the course of the infection, and individuals were sampled at multiple time points for immune gene expression analysis. We found large differences in virulence between parasite species, irrespective of host, and maladaptation of parasites to their sympatric host. The immune system responded to infection, with a decrease in expression of innate and Th1-type adaptive response genes in fish infected with the less virulent parasite, representing a marker of a possible resistance mechanism. There was no evidence of local adaptation in immune gene expression levels. Our results add to the growing understanding of the extent of host-parasite local adaptation, and demonstrate a systemic immune response during infection with a common ectoparasite. Further immunological studies using the stickleback-Gyrodactylus system can continue to contribute to our understanding of the function of the immune response in natural populations.

Spatio-temporal variation in parasite communities maintains diversity at the major histocompatibility complex class IIβ in the endangered Rio Grande silvery minnow

Climate change will strongly impact aquatic ecosystems particularly in arid and semi-arid regions. Fish-parasite interactions will also be affected by predicted altered flow and temperature regimes, and other environmental stressors. Hence, identifying environmental and genetic factors associated with maintaining diversity at immune genes is critical for understanding species' adaptive capacity. Here, we combine genetic (MHC class IIβ and microsatellites), parasitological and ecological data to explore the relationship between these factors in the remnant wild Rio Grande silvery minnow (Hybognathus amarus) population, an endangered species found in the southwestern United States. Infections with multiple parasites on the gills were observed and there was spatio-temporal variation in parasite communities and patterns of infection among individuals. Despite its highly endangered status and chronically low genetic effective size, Rio Grande silvery minnow had high allelic diversity at MHC class IIβ with more alleles recognized at the presumptive DAB1 locus compared to the DAB3 locus. We identified significant associations between specific parasites and MHC alleles against a backdrop of generalist parasite prevalence. We also found that individuals with higher individual neutral heterozygosity and higher amino acid divergence between MHC alleles had lower parasite abundance and diversity. Taken together, these results suggest a role for fluctuating selection imposed by spatio-temporal variation in pathogen communities and divergent allele advantage in maintenance of high MHC polymorphism. Understanding the complex interaction of habitat, pathogens and immunity in protected species will require integrated experimental, genetic and field studies.

Measuring the immune system of the three-spined stickleback - investigating natural variation by quantifying immune expression in the laboratory and the wild

Current understanding of the immune system comes primarily from laboratory-based studies. There has been substantial interest in examining how it functions in the wild, but studies have been limited by a lack of appropriate assays and study species. The three-spined stickleback (Gasterosteus aculeatus L.) provides an ideal system in which to advance the study of wild immunology, but requires the development of suitable immune assays. We demonstrate that meaningful variation in the immune response of stickleback can be measured using real-time PCR to quantify the expression of eight genes, representing the innate response and Th1-, Th2- and Treg-type adaptive responses. Assays are validated by comparing the immune expression profiles of wild and laboratory-raised stickleback, and by examining variation across populations on North Uist, Scotland. We also compare the immune response potential of laboratory-raised individuals from two Icelandic populations by stimulating cells in culture. Immune profiles of wild fish differed from laboratory-raised fish from the same parental population, with immune expression patterns in the wild converging relative to those in the laboratory. Innate measures differed between wild populations, whilst the adaptive response was associated with variation in age, relative size of fish, reproductive status and S. solidus infection levels. Laboratory-raised individuals from different populations showed markedly different innate immune response potential. The ability to combine studies in the laboratory and in the wild underlines the potential of this toolkit to advance our understanding of the ecological and evolutionary relevance of immune system variation in a natural setting.

Extensive copy-number variation of young genes across stickleback populations

Duplicate genes emerge as copy-number variations (CNVs) at the population level, and remain copy-number polymorphic until they are fixed or lost. The successful establishment of such structural polymorphisms in the genome plays an important role in evolution by promoting genetic diversity, complexity and innovation. To characterize the early evolutionary stages of duplicate genes and their potential adaptive benefits, we combine comparative genomics with population genomics analyses to evaluate the distribution and impact of CNVs across natural populations of an eco-genomic model, the three-spined stickleback. With whole genome sequences of 66 individuals from populations inhabiting three distinct habitats, we find that CNVs generally occur at low frequencies and are often only found in one of the 11 populations surveyed. A subset of CNVs, however, displays copy-number differentiation between populations, showing elevated within-population frequencies consistent with local adaptation. By comparing teleost genomes to identify lineage-specific genes and duplications in sticklebacks, we highlight rampant gene content differences among individuals in which over 30% of young duplicate genes are CNVs. These CNV genes are evolving rapidly at the molecular level and are enriched with functional categories associated with environmental interactions, depicting the dynamic early copy-number polymorphic stage of genes during population differentiation.

After a locus is duplicated in a genome, individuals from a population instantaneously differ in the number of copies of this locus producing a copy-number variation (CNV). Over time, the joint effects of selection and other evolutionary forces will act to either eliminate the extra genetic copy or retain it. Depending on this evolutionary interplay, young duplications, including newly duplicated genes, can persist for millions of years as CNVs. CNVs may especially be prevalent between populations that have colonized and adapted to disparate environments in which selective pressures differ. Using whole genome sequences from several populations of three-spined sticklebacks that inhabit different environments, we find that a third of young duplicated genes are CNVs. These young CNV genes are enriched with environmental response functions and evolving rapidly at the molecular level, making them promising candidates for a role in the rapid ecological adaptation to novel environments.

Absence of major histocompatibility complex class II mediated immunity in pipefish, Syngnathus typhle: evidence from deep transcriptome sequencing

The major histocompatibility complex (MHC)-mediated adaptive immune system is the hallmark of gnathostome immune defence. Recent work suggests that cod-like fishes (Gadidae) lack important components of the MHC class II mediated immunity. Here, we report a putative independent loss of functionality of this pathway in another species, the pipefish Syngnathus typhle, that belongs to a distantly related fish family (Syngnathidae). In a deep transcriptome sequencing approach comprising several independent normalized and non-normalized expressed sequence tag (EST) libraries with approximately 7.5 × 10(8) reads, sequenced with two next generation platforms (454 and Illumina), we were unable to identify MHC class IIα/β genes as well as genes encoding associated receptors. Along with the recent findings in cod, our results suggest that immune systems of the Euteleosts may be more variable than previously assumed.

The functional importance of sequence versus expression variability of MHC alleles in parasite resistance

Understanding selection processes driving the pronounced allelic polymorphism of the major histocompatibility complex (MHC) genes and its functional associations to parasite load have been the focus of many recent wildlife studies. Two main selection scenarios are currently debated which explain the susceptibility or resistance to parasite infections either by the effects of (1) specific MHC alleles which are selected frequency-dependent in space and time or (2) a heterozygote or divergent allele advantage. So far, most studies have focused only on structural variance in co-evolutionary processes although this might not be the only trait subject to natural selection. In the present study, we analysed structural variance stretching from exon1 through exon3 of MHC class II DRB genes as well as genotypic expression variance in relation to the gastrointestinal helminth prevalence and infection intensity in wild yellow-necked mice (Apodemus flavicollis). We found support for the functional importance of specific alleles both on the sequence and expression level. By resampling a previously investigated study population we identified specific MHC alleles affected by temporal shifts in parasite pressure and recorded associated changes in allele frequencies. The allele Apfl-DRB*23 was associated with resistance to infections by the oxyurid nematode Syphacia stroma and at the same time with susceptibility to cestode infection intensity. In line with our expectation, MHC mRNA transcript levels tended to be higher in cestode-infected animals carrying the allele Apfl-DRB*23. However, no support for a heterozygote or divergent allele advantage on the sequence or expression level was detected. The individual amino acid distance of genotypes did not explain individual differences in parasite loads and the genetic distance had no effect on MHC genotype expression. For ongoing studies on the functional importance of expression variance in parasite resistance, allele-specific expression data would be preferable.

Polymorphism and Balancing Selection of MHC Class II DAB Gene in 7 Selective Flounder (Paralichthys olivaceus) Families

In order to determine the genetic variation of the MHC class IIB exon2 allele in the offspring, 700 fry from seven families of Japanese flounder challenged with V. anguillarum were studied, and different mortality rates were found in those families. Five to ten surviving and dead fry from each of the seven families were selected to study the MHC class II B exon2 gene with PCR and a direct sequencing method. One hundred and sixteen different exon2 sequences were found and 116 different alleles were identified, while a minimum of four loci were revealed in the MHC class II B exon2 gene. The ratio (d_N/d_S) of nonsynonymous substitution (d_N) to synonymous substitutions (d_S) in the peptide-binding region (PBR) of the MHC class IIB gene was 6.234, which indicated that balancing selection is acting on the MHC class IIB genes. The MHC IIB alleles were thus being passed on to their progeny. Some alleles were significantly more frequent in surviving than dead individuals. All together our data suggested that the alleles Paol-DAB*4301, Paol-DAB*4601, Paol-DAB*4302, Paol-DAB*3803, and Paol-DAB*4101 were associated with resistance to V. anguillarum in flounder.

Distinct lineages of Schistocephalus parasites in threespine and ninespine stickleback hosts revealed by DNA sequence analysis

Parasitic interactions are often part of complex networks of interspecific relationships that have evolved in biological communities. Despite many years of work on the evolution of parasitism, the likelihood that sister taxa of parasites can co-evolve with their hosts to specifically infect two related lineages, even when those hosts occur sympatrically, is still unclear. Furthermore, when these specific interactions occur, the molecular and physiological basis of this specificity is still largely unknown. The presence of these specific parasitic relationships can now be tested using molecular markers such as DNA sequence variation. Here we test for specific parasitic relationships in an emerging host-parasite model, the stickleback-Schistocephalus system. Threespine and ninespine stickleback fish are intermediate hosts for Schistocephalus cestode parasites that are phenotypically very similar and have nearly identical life cycles through plankton, stickleback, and avian hosts. We analyzed over 2000 base pairs of COX1 and NADH1 mitochondrial DNA sequences in 48 Schistocephalus individuals collected from threespine and ninespine stickleback hosts from disparate geographic regions distributed across the Northern Hemisphere. Our data strongly support the presence of two distinct clades of Schistocephalus, each of which exclusively infects either threespine or ninespine stickleback. These clades most likely represent different species that diverged soon after the speciation of their stickleback hosts. In addition, genetic structuring exists among Schistocephalus taken from threespine stickleback hosts from Alaska, Oregon and Wales, although it is much less than the divergence between hosts. Our findings emphasize that biological communities may be even more complex than they first appear, and beg the question of what are the ecological, physiological, and genetic factors that maintain the specificity of the Schistocephalus parasites and their stickleback hosts.

Characterization and evolution of MHC class II B genes in Ardeid birds

Major histocompatibility complex (MHC) is a multi-gene family that is very suitable to investigate a wide range of open questions in evolutionary ecology. In this study, we characterized two expressed MHC class II B genes (DAB1 and DAB2) in the Grey Heron (Aves: Ardea cinerea). We further developed the primer pairs to amplify and sequence two MHC class II B loci in ten ardeid birds. Phylogenetic analysis revealed that different parts of the genes showed different evolutionary patterns. The exon 2 sequences tended to cluster two gene-specific lineages. In each lineage, exon 2 sequences from several species showed closer relationships than sequences within species, and two shared identical alleles were found between species (Egretta sacra and Nycticorax nycticorax; Egretta garzetta and Bubulcus ibis), supporting the hypothesis of trans-species polymorphism. In contrast, the species-specific intron 2 plus partial exon 3 tree suggested that DAB1 and DAB2 were subject to concerted evolution. GENECONV analyses showed the gene exchange played an important role in the ardeid MHC evolution.

Putative causes and consequences of MHC variation within and between locally adapted stickleback demes

Genes of the major histocompatibility complex (MHC) have been a source of considerable research interest, owing in large part to the growing body of evidence that they may be subject to both natural and sexual selection. However, much remains to be learned about the dynamics of MHC genes in subdivided populations, particularly those characterized by divergent ecological pressures. In this study, we attempt to disentangle the relative roles of both parasite-mediated selection and MHC-mediated mate choice in an open estuarine system inhabited by two parapatric, adaptively divergent threespine stickleback (Gasterosteus aculeatus) demes. We sequenced the putative peptide-binding region (PBR) of an estimated four Class IIβ loci from 127 individuals, identifying 329 sequence variants (276 translated amino acid sequences). Demes differed significantly both in the frequency of MHC alleles and in the communities of helminth parasites infecting resident sticklebacks. Strong signatures of natural selection were inferred from analyses of codon substitutions, particularly in the derived (freshwater) rather than the ancestral (marine) deme. Relationships between parasite load and MHC diversity were indicative of balancing selection, but only within the freshwater deme. Signals of MHC-mediated mate choice were weak and differed significantly between demes. Moreover, MHC-mediated mate choice was significantly influenced by environmental salinity and appeared of secondary importance to tendencies towards assortative mating. We discuss the implications of these findings in respect to ecological adaptation and the potential demographic consequences of possible outcomes of MHC-mediated mate choice.

Major histocompatibility complex polymorphism: dynamics and consequences of parasite-mediated local adaptation in fishes

Parasitism is a common form of life and represents a strong selective pressure for host organisms. In response to this evolutionary pressure, vertebrates have developed genetically coded defences such as the major histocompatibility complex (MHC). Mechanisms of parasite-mediated selection not only maintain outstanding polymorphism in these genes but have also been proposed to further promote host population divergence and ultimately speciation because it can drive evolution of local adaptation in which MHC genes play a crucial role. This review first highlights the dynamics and complexity of parasite-mediated selection in natural systems, which not only depends on dominating parasite strategies and on the taxonomic diversity of the parasite community but also includes the differences in parasite communities between habitats and niches, creating divergent selection on locally adapted populations. Then the different ways in which MHC genes potentially allow vertebrates to respond to these dynamics and to adapt locally are outlined. Finally, it is proposed that varying selection strength in time and space may lead to variation in the strength of precopulatory reproductive isolation which has evolved to maintain local adaptation.

How pathogens drive genetic diversity: MHC, mechanisms and misunderstandings

Major histocompatibility complex (MHC) genes have been put forward as a model for studying how genetic diversity is maintained in wild populations. Pathogen-mediated selection (PMS) is believed to generate the extraordinary levels of MHC diversity observed. However, establishing the relative importance of the three proposed mechanisms of PMS (heterozygote advantage, rare-allele advantage and fluctuating selection) has proved extremely difficult. Studies have attempted to differentiate between mechanisms of PMS using two approaches: (i) comparing MHC diversity with that expected under neutrality and (ii) relating MHC diversity to pathogen regime. Here, we show that in many cases the same predictions arise from the different mechanisms under these approaches, and that most studies that have inferred one mechanism of selection have not fully considered the alternative explanations. We argue that, while it may be possible to demonstrate that particular mechanisms of PMS are occurring, resolving their relative importance within a system is probably impossible. A more realistic target is to continue to demonstrate when and where the different mechanisms of PMS occur, with the aim of determining their relative importance across systems. We put forward what we believe to be the most promising approaches that will allow us to progress towards achieving this.

Does genetic diversity predict health in humans?

Genetic diversity, especially at genes important for immune functioning within the Major Histocompatibility Complex (MHC), has been associated with fitness-related traits, including disease resistance, in many species. Recently, genetic diversity has been associated with mate preferences in humans. Here we asked whether these preferences are adaptive in terms of obtaining healthier mates. We investigated whether genetic diversity (heterozygosity and standardized mean d²) at MHC and nonMHC microsatellite loci, predicted health in 153 individuals. Individuals with greater allelic diversity (d²) at nonMHC loci and at one MHC locus, linked to HLA-DRB1, reported fewer symptoms over a four-month period than individuals with lower d². In contrast, there were no associations between MHC or nonMHC heterozygosity and health. NonMHC-d² has previously been found to predict male preferences for female faces. Thus, the current findings suggest that nonMHC diversity may play a role in both natural and sexual selection acting on human populations.

Does intra-individual major histocompatibility complex diversity keep a golden mean?

An adaptive immune response is usually initiated only if a major histocompatibility complex (MHC) molecule presents pathogen-derived peptides to T-cells. Every MHC molecule can present only peptides that match its peptide-binding groove. Thus, it seems advantageous for an individual to express many different MHC molecules to be able to resist many different pathogens. However, although MHC genes are the most polymorphic genes of vertebrates, each individual has only a very small subset of the diversity at the population level. This is an evolutionary paradox. We provide an overview of the current data on infection studies and mate-choice experiments and conclude that overall evidence suggests that intermediate intra-individual MHC diversity is optimal. Selective forces that may set an upper limit to intra-individual MHC diversity are discussed. An updated mathematical model based on recent findings on T-cell selection can predict the natural range of intra-individual MHC diversity. Thus, the aim of our review is to evaluate whether the number of MHC alleles usually present in individuals may be optimal to balance the advantages of presenting an increased range of peptides versus the disadvantages of an increased loss of T-cells.

Cell markers and determinants in fish immunology

Despite the impressive increase in the cloning and expression of genes encoding fish immunoregulatory molecules, the knowledge on "in vivo" and "in vitro" functional immunology of the corresponding peptide products is still at an initial stage. This is partly due to the lacking of specific markers for immunoregulatory peptides, that represent an indispensible tool to dissect immune reactions and to trace the fate of cellular events downstream of the activation. In this review we summarise the available information on functional immune activities of some teleost species and discuss the obtained data in an evolutionary and applied context.

Mortality selection during the 2003 European heat wave in three-spined sticklebacks: effects of parasites and MHC genotype

Background

Ecological interaction strength may increase under environmental stress including temperature. How such stress enhances and interacts with parasite selection is almost unknown. We studied the importance of resistance genes of the major histocompatibility complex (MHC) class II in 14 families of three-spined sticklebacks Gasterosteus aculeatus exposed to their natural macroparasites in field enclosures in the extreme summer of 2003.

Results

After a mass die-off during the 2003-European heat wave killing 78% of 277 experimental fish, we found strong differences in survival among and within families. In families with higher average parasite load fewer individuals survived. Multivariate analysis revealed that the composition of the infecting parasite fauna was family specific. Within families, individuals with an intermediate number of MHC class IIB sequence variants survived best and had the lowest parasite load among survivors, suggesting a direct functional link between MHC diversity and fitness. The within family MHC effects were, however, small compared to between family effects, suggesting that other genetic components or non-genetic effects were also important.

Conclusion

The correlation between parasite load and mortality that we found at both individual and family level might have appeared only in the extraordinary heatwave of 2003. Due to global warming the frequency of extreme climatic events is predicted to increase, which might intensify costs of parasitism and enhance selection on immune genes.

Effect of multiple genetic polymorphisms on antigen presentation and susceptibility to Mycobacterium tuberculosis infection

Several molecules related to antigen presentation, including gamma interferon (IFN-gamma) and the major histocompatibility complex (MHC), are encoded by polymorphic genes. Some polymorphisms were found to affect susceptibility to tuberculosis (TB) when they were considered singly in epidemiological studies, but how multiple polymorphisms interact to determine susceptibility to TB in an individual remains an open question. We hypothesized that polymorphisms in some genes may counteract or intensify the effects of polymorphisms in other genes. For example, an increase in IFN-gamma expression may counteract the weak binding that a particular MHC variant displays for a peptide from Mycobacterium tuberculosis to establish the same T-cell response as another, more strongly binding MHC variant. To test this hypothesis, we developed a mathematical model of antigen presentation based on experimental data for the known effects of genetic polymorphisms and simulated time courses when multiple polymorphisms were present. We found that polymorphisms in different genes could affect antigen presentation to the same extent and therefore compensate for each other. Furthermore, we defined the conditions under which such relationships could exist. For example, increased IFN-gamma expression compensated for decreased peptide-MHC affinity in the model only above a certain threshold of expression. Below this threshold, changes in IFN-gamma expression were ineffectual compared to changes in peptide-MHC affinity. The finding that polymorphisms exhibit such relationships could explain discrepancies in the epidemiological literature, where some polymorphisms have been inconsistently associated with susceptibility to TB. Furthermore, the model allows polymorphisms to be ranked by effect, providing a new tool for designing association studies.

Housekeeping genes for quantitative expression studies in the three-spined stickleback Gasterosteus aculeatus

BACKGROUND

During the last years the quantification of immune response under immunological challenges, e.g. parasitation, has been a major focus of research. In this context, the expression of immune response genes in teleost fish has been surveyed for scientific and commercial purposes. Despite the fact that it was shown in teleostei and other taxa that the gene for beta-actin is not the most stably expressed housekeeping gene (HKG), depending on the tissue and experimental treatment, the gene has been used as a reference gene in such studies. In the three-spined stickleback, Gasterosteus aculeatus, other HKG than the one for beta-actin have not been established so far.

RESULTS

To establish a reliable method for the measurement of immune gene expression in Gasterosteus aculeatus, sequences from the now available genome database and an EST library of the same species were used to select oligonucleotide primers for HKG, in order to perform quantitative reverse-transcription (RT) PCR. The expression stability of ten candidate reference genes was evaluated in three different tissues, and in five parasite treatment groups, using the three algorithms BestKeeper, geNorm and NormFinder. Our results showed that in most of the tissues and treatments HKG that could not be used so far due to unknown sequences, proved to be more stably expressed than the one for beta-actin.

CONCLUSION

As they were the most stably expressed genes in all tissues examined, we suggest using the genes for the L13a ribosomal binding protein and ubiquitin as alternative or additional reference genes in expression analysis in Gasterosteus aculeatus.

Nutritional state influences shoaling preference for familiars

Preferences for grouping with familiar individuals are shown in many animal species, including the three-spined stickleback (Gasterosteus aculeatus). Shoaling with familiars is advantageous because of more precise anti-predator behaviours or more stable dominance hierarchies. Additionally, associations with familiar individuals facilitate the evolution of altruistic behaviour. Thus, in situations of increased competition one might expect an increased preference for familiar fish. We gave single juvenile sticklebacks of different nutritional state the choice between shoals composed either of familiar or unfamiliar individuals. Satiated fish preferred to shoal with familiar individuals. A comparative analysis of 8 stickleback studies with 15 different tests using familiars showed that all tests gave similar results, i.e. sticklebacks of all age classes preferred to shoal with familiars in a non-sexual context. In contrast, hungry test fish did not prefer to shoal with familiar fish, but even showed a preference for the unfamiliar group. Because sticklebacks use early-life familiarity to recognize kin, the results suggest the avoidance of competition with relatives. To our knowledge, this is the first study showing an impact of nutritional state on social interactions with familiar individuals.

Characterization of antisera raised against stickleback (Gasterosteus aculeatus) MHC class I and class II molecules

The three-spined stickleback (Gasterosteus aculeatus) is an important model organism for investigations on the maintenance of polymorphism of the major histocompatibility complex (MHC) of vertebrates. Analysis of functional aspects of MHC diversity in stickleback would benefit from the availability of MHC specific reagents. Here we characterize antisera raised against recombinant fusion proteins of stickleback MHC class I alpha and class II alpha and beta. Western blot analysis using recombinant proteins confirmed the specificity of the antisera. In brain and muscle preparations, neither of the MHC types was detectable. High levels of each MHC receptor type were observed in gills and spleen and lower levels in head kidneys. In histological sections of gills, epithelial cells of primary and secondary lamellae stained positive with MHC class I antiserum, while single, scattered cells stained positive for MHC class II. In sections of spleen and head kidney, considerable numbers of cells positive for either MHC type were detected. Molecular weight shift in SDS-PAGE after deglycosylation of MHC class I alpha and class II beta confirmed the predicted glyco-protein character of the molecules. The majority of MHC II alpha was not glycosylated; only a small fraction of MHC II alpha was susceptible to deglycosylation. This suggests differential expression of the two stickleback MHC II alpha genes (Gaac-DAA, Gaac-DBA) only one of which (Gaac-DBA) has a site for N-linked glycosylation.

Characterization, polymorphism, and evolution of MHC class II B genes in birds of prey

During the last decade, the major histocompatibility complex (MHC) has received much attention in the fields of evolutionary and conservation biology because of its potential implications in many biological processes. New insights into the gene structure and evolution of MHC genes can be gained through study of additional lineages of birds not yet investigated at the genomic level. In this study, we characterized MHC class II B genes in five families of birds of prey (Accipitridae, Pandionidae, Strigidae, Tytonidae, and Falconidae). Using PCR approaches, we isolated genomic MHC sequences up to 1300 bp spanning exons 1 to 3 in 26 representatives of each raptor lineage, finding no stop codons or frameshift mutations in any coding region. A survey of diversity across the entirety of exon 2 in the lesser kestrel Falco naumanni reported 26 alleles in 21 individuals. Bayesian analysis revealed 21 positively selected amino acid sites, which suggests that the MHC genes described here are functional and probably expressed. Finally, through interlocus comparisons and phylogenetic analysis, we also discuss genetic evidence for concerted and transspecies evolution in the raptor MHC.

Molecular cloning, differential expression and 3D structural analysis of the MHC class-II beta chain from sea bass (Dicentrarchus labrax L.)

The major histocompatibility complex class I and II molecules (MHC-I and MHC-II) play a pivotal role in vertebrate immune response to antigenic peptides. In this paper we report the cloning and sequencing of the MHC class II beta chain from sea bass (Dicentrarchus labrax L.). The six obtained cDNA sequences (designated as Dila-DAB) code for 250 amino acids, with a predicted 21 amino acid signal peptide and contain a 28bp 5'-UTR and a 478bp 3'-UTR. A multiple alignment of the predicted translation of the Dila-DAB sequences was assembled together with other fish and mammalian sequences and it showed the conservation of most amino acid residues characteristic of the MHC class II beta chain structure. The highest basal Dila-DAB expression was found in gills, followed by gut and thymus, lower mRNA levels were found in spleen, peripheral blood leucocytes (PBL) and liver. Stimulation of head kidney leukocytes with LPS for 4h showed very little difference in the Dila-DAB expression, but after 24h the Dila-DAB level decreased to a large extent and the difference was statistically significant. Stimulation of head kidney leukocytes with different concentrations of rIL-1beta (ranging from 0 to 100ng/ml) resulted in a dose-dependent reduction of the Dila-DAB expression. Moreover, two 3D Dila-DAB*0101 homology models were obtained based on crystallographic mouse MHC-II structures complexed with D10 T-cell antigen receptor or human CD4; features and differences between the models were evaluated and discussed. Taken together these results are of interest as MHC-II structure and function, molecular polymorphism and differential gene expression are in correlation with disease resistance to virus and bacteria in teleost fish.

Individual MHC class I and MHC class IIBdiversities 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.

Contrasting mode of evolution between the MHC class I genomic region and class II region in the three-spined stickleback (Gasterosteus aculeatus L.; Gasterosteidae: Teleostei)

Major histocompatibility complex (MHC) class I molecules display peptides on cell surfaces for subsequent T-cell recognition and are involved in the immune response against intracellular pathogens. In this study, a BAC library was created from a single three-spined stickleback and screened for clones containing MHC class I genes. In a 163.2-kb genomic sequence segment of a single clone, we identified three MHC class I genes in the same transcriptional orientation. Two class I genes are potentially expressed and functional. In one class I gene, the transmembrane region is missing and could therefore present a pseudogene. Alternatively, it presents a functional gene that encodes a soluble MHC class Ib molecule. Despite genomic similarities to the MHC class II region, which is characterized by interlocus recombination, we did not find any evidence for this kind of recombination in the class I genes. It thus seems that interlocus recombination may play a rather minor role in generating class I diversity in stickleback and that the class I region displays a higher genomic stability (i.e., lower local recombination rate). In addition, two non-MHC genes (Oct-2 beta and Na(+),K(+)-ATPasealpha3) have been identified in the analyzed class I region. The Oct-2 beta gene is a transcription factor that is expressed primarily in B lymphocytes, in activated T-cells, and in neuronal cells. The Na(+),K(+)-ATPasealpha3 gene is primarily expressed in the brain and heart and mediates catalytic activities. Both genes are located on the same linkage group together with the MHC class I genes in the zebra fish. In humans, however, homologues of Oct-2 beta and ATPasealpha3 lie outside the MHC region, which indicates that the concentration of immune genes found in mammalian genomes is a derived state.

MHC genes and oxidative stress in sticklebacks: an immuno-ecological approach

Individual variation in the susceptibility to infection may result from the varying ability of hosts to specifically recognize different parasite strains. Alternatively, there could be individual host differences in fitness costs of immune defence. Although, these two explanations are not mutually exclusive, they have so far been treated in separate experimental approaches. To analyse potential relationships, we studied body condition and oxidative stress, which may reflect costs of immunity, in three-spined sticklebacks that had been experimentally exposed to three species of naturally occurring parasite. These sticklebacks differed in a trait, which is crucial to specific parasite defence, i.e. individual genetic diversity at major histocompatibility complex (MHC) class IIB loci. Oxidative stress was quantified as tissue acrolein, a technique that has been applied to questions of immuno-ecology for the first time. We measured gene expression at the MHC and other estimates of immune activation. We found that fish with high levels of MHC expression had poor condition and elevated oxidative stress. These results indicate that MHC-based specific immunity is connected with oxidative stress. They could, thus, also be relevant in the broader context of the evolution of sexually selected signals that are based on carotenoids and are, thus supposed to reflect oxidative stress resistance.