Using a bacterial artificial chromosome (BAC) library, we analysed a 99.5 kb genomic segment containing the major histocompatibility class II genes of a teleost, the three-spined stickleback Gasterosteus aculeatus. Experiments with G. aculeatus have provided direct evidence for balancing selection by pathogens and mate choice driving MH class II beta polymorphism. Two sets of paralogous class II alpha genes and beta genes in a tandem arrangement were identified, designated Gaac-DAA/DAB and Gaac-DBA/DBB. Expression analysis of the beta genes using single-strand conformation polymorphism revealed that both gene copies are expressed. Based on an analysis of pairwise nucleotide polymorphisms, we estimate that the duplication into two paralogous class II loci occurred only 1.2-2.4 million years ago, 1-2 orders of magnitude more recently than in other fish, bird or mammalian species. At the 3'-direction of the classical MH loci, we identified another seven genes or gene fragments, two of which (small inducible cytokine, complement regulatory factor) are related to immune function in other vertebrates. None of these genes were associated with MH class II genes in zebrafish, suggesting a markedly different organisation of the MH class II region in sticklebacks, and presumably, across bony fishes. When the nucleotide substitution pattern of the novel class II beta genes was analysed together with a representative sequence sample isolated from fish in northern Germany (n=27), we found that the peptide binding region of the Gaac-DAB and Gaac-DBB loci had undergone an inter-locus gene conversion (P=0.007). In accordance, we found a 10- to 20-fold higher frequency of CpG-islands on the MH class II segment compared to other species, a feature that may be conducive for inter-locus recombination.

Major histocompatibility complex (MHC) antigen-presenting genes are the most variable loci in vertebrate genomes. Host-parasite co-evolution is assumed to maintain the excessive polymorphism in the MHC loci. However, the molecular mechanisms underlying the striking diversity in the MHC remain contentious. The extent to which recombination contributes to the diversity at MHC loci in natural populations is still controversial, and there have been only few comparative studies that make quantitative estimates of recombination rates. In this study, we performed a comparative analysis for 15 different ungulates species to estimate the population recombination rate, and to quantify levels of selection. As expected for all species, we observed signatures of strong positive selection, and identified individual residues experiencing selection that were congruent with those constituting the peptide-binding region of the human DRB gene. However, in addition for each species, we also observed recombination rates that were significantly different from zero on the basis of likelihood-permutation tests, and in other non-quantitative analyses. Patterns of synonymous and non-synonymous sequence diversity were consistent with differing demographic histories between species, but recent simulation studies by other authors suggest inference of selection and recombination is likely to be robust to such deviations from standard models. If high rates of recombination are common in MHC genes of other taxa, re-evaluation of many inference-based phylogenetic analyses of MHC loci, such as estimates of the divergence time of alleles and trans-specific polymorphism, may be required.

Parasites can strongly affect the evolution of their hosts, but their effects on host diversification are less clear. In theory, contrasting parasite communities in different foraging habitats could generate divergent selection on hosts and promote ecological speciation. Immune systems are costly to maintain, adaptable, and an important component of individual fitness. As a result, immune system genes, such as those of the Major Histocompatability Complex (MHC), can change rapidly in response to parasite-mediated selection. In threespine stickleback (Gasterosteus aculeatus), as well as in other vertebrates, MHC genes have been linked with female mating preference, suggesting that divergent selection acting on MHC genes might influence speciation. Here, we examined genetic variation at MHC Class II loci of sticklebacks from two lakes with a limnetic and benthic species pair, and two lakes with a single species. In both lakes with species pairs, limnetics and benthics differed in their composition of MHC alleles, and limnetics had fewer MHC alleles per individual than benthics. Similar to the limnetics, the allopatric population with a pelagic phenotype had few MHC alleles per individual, suggesting a correlation between MHC genotype and foraging habitat. Using a simulation model we show that the diversity and composition of MHC alleles in a sympatric species pair depends on the amount of assortative mating and on the strength of parasite-mediated selection in adjacent foraging habitats. Our results indicate parallel divergence in the number of MHC alleles between sympatric stickleback species, possibly resulting from the contrasting parasite communities in littoral and pelagic habitats of lakes.

Reduced FCGR3B copy number is associated with increased risk of systemic lupus erythematosus (SLE). The five FCGR2/FCGR3 genes are arranged across two highly paralogous genomic segments on chromosome 1q23. Previous studies have suggested mechanisms for structural rearrangements at the FCGR2/FCGR3 locus and have proposed mechanisms whereby altered FCGR3B copy number predisposes to autoimmunity, but the high degree of sequence similarity between paralogous segments has prevented precise definition of the molecular events and their functional consequences. To pursue the genomic pathology associated with FCGR3B copy-number variation, we integrated sequencing data from fosmid and bacterial artificial chromosome clones and sequence-captured DNA from FCGR3B-deleted genomes to establish a detailed map of allelic and paralogous sequence variation across the FCGR2/FCGR3 locus. This analysis identified two highly paralogous 24.5 kb blocks within the FCGR2C/FCGR3B/FCGR2B locus that are devoid of nonpolymorphic paralogous sequence variations and that define the limits of the genomic regions in which nonallelic homologous recombination leads to FCGR2C/FCGR3B copy-number variation. Further, the data showed evidence of swapping of haplotype blocks between these highly paralogous blocks that most likely arose from sequential ancestral recombination events across the region. Functionally, we found by flow cytometry, immunoblotting and cDNA sequencing that individuals with FCGR3B-deleted alleles show ectopic presence of FcγRIIb on natural killer (NK) cells. We conclude that FCGR3B deletion juxtaposes the 5'-regulatory sequences of FCGR2C with the coding sequence of FCGR2B, creating a chimeric gene that results in an ectopic accumulation of FcγRIIb on NK cells and provides an explanation for SLE risk associated with reduced FCGR3B gene copy number.

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.

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.

BackgroundRecent studies of selection on mitochondrial (mt) OXPHOS genes suggest adaptation due mainly to environmental variation. In this context, Tunisian hares that display several external phenotypes with phylogenetically rather homogenous gene pool and shallow population structure provide a good precondition to detect positive selection on mt genes related to environmental/climatic variation, specifically ambient temperature and precipitation.ResultsWe used codon-based methods along with population genetic data to test for positive selection on ATPase synthase 6 (ATP6) and NADH dehydrogenase 2 (ND2) of cape hares (Lepus capensis) collected along a steep ecological gradient in Tunisia. We found significantly higher differentiation at the ATP6 locus across Tunisia, with sub-humid Mediterranean, semi-arid, and arid Sahara climate than for fourteen unlinked supposedly neutrally evolving nuclear microsatellites and mt control region sequences. This suggested positive selection on ATP6 sequences, which was confirmed by several codon-based tests for one sequence site that together with a second site translated into four different amino acids. Positive selection on ND2 sequences was also confirmed by several codon-based tests. The corresponding frequencies of the two most prevalent variants at each locus varied significantly across climate regions, and our logistic general linear models of occurrence of those proteins indicated significant effects of mean annual temperature for ATP6 and mean minimum temperature of the coldest month of the year for ND2, independent of geographical location, annual precipitation, and the respective co-occurring protein at the second locus. Moreover, presence of the ancestral ATP6 protein, as inferred from phylogenetic networks, was positively affected by the simultaneous presence of the derived ND2 protein and vice versa, independent of temperature, precipitation, or geographic location. Finally, we obtained a significant coevolution signal for the ancestral ATP6 and derived ND2 sequences and vice versa. Conclusionspositive selection was strongly suggested by the population genetic approach and the codon-based tests in both mtDNA genes. Moreover, the two most prevalent proteins at the ATP6 locus were distributed at significantly varying frequencies across the study area with a significant effect of mean annual temperature on the occurrence of the ATP6 proteins independent of geographical coordinates and the co-occuring ND2 protein variant. For ND2, occurrence of the two most frequent protein variants was significantly influenced by the mean minimum temperature of the coldest month, independent of the co-occurring ATP6 protein variant and geographical coordinates. This strongly suggests direct involvement of ambient temperature in the adaptation of the studied mtOXPHOS genes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-017-0896-0) contains supplementary material, which is available to authorized users.

SummaryThe causes of autoimmune disease remain poorly defined. However, it is known that genetic factors contribute to disease susceptibility. Hitherto, studies have focused upon single nucleotide polymorphisms as both tools for mapping and as probable causal variants. Recent studies, using genome-wide analytical techniques, have revealed that, in the genome, segments of DNA ranging in size from kilobases to megabases can vary in copy number. These changes of DNA copy number represent an important element of genomic polymorphism in humans and in other species and may therefore make a substantial contribution to phenotypic variation and population differentiation. Furthermore, copy number variation (CNV) in genomic regions harbouring dosage-sensitive genes may cause or predispose to a variety of human genetic diseases. Several recent studies have reported an association between CNV and autoimmunity in humans such as systemic lupus, psoriasis, Crohn's disease, rheumatoid arthritis and type 1 diabetes. The use of novel analytical techniques facilitates the study of complex human genomic structures such as CNV, and allows new susceptibility loci for autoimmunity to be found that are not readily mappable by single nucleotide polymorphism-based association analyses alone.