Search

A comparative analysis of 15 different ungulates species is undertaken to estimate the population recombination rate and to quantify levels of selection. Signatures of strong positive selection were observed and individual residues experiencing selection that were congruent with those constituting the peptide binding region of the human DRB gene were identified and in each species it was seen that recombination rates were significantly different from zero on the basis of likelihood-permutation tests.

Spatial pattern analysis of the variability of mitochondrial DNA and biparentally inherited genes of 443 chamois from Eastern Alps indicate that 67.09% of the variance is distributed among the four mitochondrial phylogenetic groups, while partioning among regional samples within the populations and among individuals within samples are 8.04% and 24.86%, respectively.

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 II B 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 II B 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. [ABSTRACT FROM AUTHOR]

Abstract Background In mammals, males typically have shorter lives than females. This difference is thought to be due to behavioural traits which enhance competitive abilities, and hence male reproductive success, but impair survival. Furthermore, in many species males usually show higher parasite burden than females. Consequently, the intensity of selection for genetic factors which reduce susceptibility to pathogens may differ between sexes. High variability at the major histocompatibility complex (MHC) genes is believed to be advantageous for detecting and combating the range of infectious agents present in the environment. Increased heterozygosity at these immune genes is expected to be important for individual longevity. However, whether males in natural populations benefit more from MHC heterozygosity than females has rarely been investigated. We investigated this question in a long-term study of free-living Alpine chamois (Rupicapra rupicapra), a polygynous mountain ungulate. Results Here we show that male chamois survive significantly (P = 0.022) longer if heterozygous at the MHC class II DRB locus, whereas females do not. Improved survival of males was not a result of heterozygote advantage per se, as background heterozygosity (estimated across twelve microsatellite loci) did not change significantly with age. Furthermore, reproductively active males depleted their body fat reserves earlier than females leading to significantly impaired survival rates in this sex (P < 0.008). This sex-difference was even more pronounced in areas affected by scabies, a severe parasitosis, as reproductively active males were less likely to survive than females. However, we did not find evidence for a survival advantage associated with specific MHC alleles in areas affected by scabies. Conclusions Increased MHC class II DRB heterozygosity with age in males, suggests that MHC heterozygous males survive longer than homozygotes. Reproductively active males appear to be less likely to survive than females most likely because of the energetic challenge of the winter rut, accompanied by earlier depletion of their body fat stores, and a generally higher parasite burden. This scenario renders the MHC-mediated immune response more important for males than for females, which implies a relatively stronger selection pressure on MHC genes in males than in females.

Abstract Background Mitotic recombination is important for inactivating tumour suppressor genes by copy-neutral loss of heterozygosity (LOH). Although meiotic recombination maps are plentiful, little is known about mitotic recombination. The APC gene (chr5q21) is mutated in most colorectal tumours and its usual mode of LOH is mitotic recombination. Methods We mapped mitotic recombination boundaries ("breakpoints") between the centromere (~50 Mb) and APC (~112 Mb) in early colorectal tumours. Results Breakpoints were non-random, with the highest frequency between 65 Mb and 75 Mb, close to a low copy number repeat region (68–71 Mb). There were, surprisingly, few breakpoints close to APC, contrary to expectations were there constraints on tumorigenesis caused by uncovering recessive lethal alleles or if mitotic recombination were mechanistically favoured by a longer residual chromosome arm. The locations of mitotic and meiotic recombination breakpoints were correlated, suggesting that the two types of recombination are influenced by similar processes, whether mutational or selective in origin. Breakpoints were also associated with higher local G+C content. The recombination and gain/deletion breakpoint maps on 5q were not, however, associated, perhaps owing to selective constraints on APC dosage in early colorectal tumours. Since polymorphisms within the region of frequent mitotic recombination on 5q might influence the frequency of LOH, we tested the 68–71 Mb low copy number repeat and nearby tagSNPs, but no associations with colorectal cancer risk were found. Conclusion LOH on 5q is non-random, but local factors do not greatly influence the rate of LOH at APC or explain inter differential susceptibility to colorectal tumours.

Asia is home to diverse hunter-gatherer populations characterized by significant morphological, anthropological, cultural, and linguistic diversity. Despite their importance in understanding ancestral human subsistence, little is known about the essential genetic adaptations of these groups. This study investigates the evolutionary pressures shaping the genome of the Maniq population, a nomadic hunter-gatherer group inhabiting the rainforests of southern Thailand. Using genome-wide approaches, including iHS, xp-EHH, PBE, and beta statistics, we identified signatures of positive and balancing selection. Genes under positive selection were enriched in pathways related to immunity, metabolic regulation, structural adaptation, cardiovascular performance, and neuromodulatory traits. Several genes associated with the Southeast Asian 'negrito-like' phenotype were also under positive selection. Balancing selection was primarily detected in immune-related genes, particularly within the HLA region, underscoring the critical role of genetic diversity in surviving pathogen-rich environments. Additionally, balancing selection in olfactory receptor genes highlights their importance in environmental sensing and adaptation. These results reveal the intricate interplay of positive and balancing selection in shaping the genetic landscape of the Maniq population and highlight their adaptations to the ecological and lifestyle challenges of life in the rainforest. This study contributes to our understanding of human evolutionary processes in tropical environments and hunter-gatherer societies., Competing Interests: Competing interests: The authors declare no competing interests. Ethics declarations: This study was approved by the Ethics Committee of the University of Vienna (reference number 00444) and the Human Research Ethics Committee of Khon Kaen University (reference number HE622223)., (© 2025. The Author(s).)

The Maniq of southern Thailand is one of the last remaining practicing hunter-gatherer communities in the world. However, our knowledge on their genetic origins and demographic history is still largely limited. We present here the genotype data covering ∼2.3 million single nucleotide polymorphisms of 11 unrelated Maniq individuals. Our analyses reveal the Maniq to be closely related to the Semang populations of Malaysia (Malay Negritos), who altogether carry an Andamanese-related ancestry linked to the ancient Hòabìnhian hunter-gatherers of Mainland Southeast Asia (MSEA). Moreover, the Maniq possess ∼35% East Asian-related ancestry, likely brought about by recent admixture with surrounding agriculturist communities in the region. In addition, the Maniq exhibit one of the highest levels of genetic differentiation found among living human populations, indicative of their small population size and historical practice of endogamy. Similar to other hunter-gatherer populations of MSEA, we also find the Maniq to possess low levels of Neanderthal ancestry and undetectable levels of Denisovan ancestry. Altogether, we reveal the Maniq to be a Semang group that experienced intense genetic drift and exhibits signs of ancient Hòabìnhian ancestry., (© The Author(s) 2022. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

ALDH2 is a key enzyme in alcohol metabolism that protects cells from acetaldehyde toxicity. Using iHS, iSAFE and F ST statistics, we identified regulatory acting variants affecting ALDH2 gene expression under positive selection in populations of European ancestry. Several SNPs (rs3184504, rs4766578, rs10774625, rs597808, rs653178, rs847892, rs2013002) that function as eQTLs for ALDH2 in various tissues showed evidence of strong positive selection. Very large pairwise F ST values indicated high genetic differentiation at these loci between populations of European ancestry and populations of other global ancestries. Estimating the timing of positive selection on the beneficial alleles suggests that these variants were recently adapted approximately 3000-3700 years ago. The derived beneficial alleles are in complete linkage disequilibrium with the derived ALDH2 promoter variant rs886205, which is associated with higher transcriptional activity. The SNPs rs4766578 and rs847892 are located in binding sequences for the transcription factor HNF4A, which is an important regulatory element of ALDH2 gene expression. In contrast to the missense variant ALDH2 rs671 (ALDH2*2), which is common only in East Asian populations and is associated with greatly reduced enzyme activity and alcohol intolerance, the beneficial alleles of the regulatory variants identified in this study are associated with increased expression of ALDH2. This suggests adaptation of Europeans to higher alcohol consumption., (© 2022. The Author(s).)

Enriching mitochondrial DNA (mtDNA) for sequencing entire mitochondrial genomes (mitogenomes) can be achieved by single long-range PCR. This avoids interference from the omnipresent nuclear mtDNA sequences (NUMTs). The approach is currently restricted to the use of samples collected from humans and ray-finned fishes. Here, we extended the use of single long-range PCR by introducing back-to-back oligonucleotides that target a sequence of extraordinary homology across vertebrates. The assay was applied to five hibernating rodents, namely alpine marmot, Arctic and European ground squirrels, and common and garden dormice, four of which have not been fully sequenced before. Analysis of the novel mitogenomes focussed on the prediction of mitochondrial-derived peptides (MDPs) providing another level of information encoded by mtDNA. The comparison of MOTS-c, SHLP4 and SHLP6 sequences across vertebrate species identified segments of high homology that argue for future experimentation. In addition, we evaluated four candidate polymorphisms replacing an amino acid in mitochondrially encoded subunits of the oxidative phosphorylation (OXPHOS) system that were reported in relation to cold-adaptation. No obvious pattern was found for the diverse sets of mammalian species that either apply daily or multiday torpor or otherwise cope with cold. In summary, our single long-range PCR assay applying a pair of back-to-back primers that target a consensus sequence motif of Vertebrata has potential to amplify (intact) mitochondrial rings present in templates from a taxonomically diverse range of vertebrates. It could be promising for studying novel mitogenomes, mitotypes of a population and mitochondrial heteroplasmy in a sensitive, straightforward and flexible manner., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Emser, Schaschl, Millesi and Steinborn.)

Five X-HxIP (Hx-amides) 6a-e, in which the N-terminus p-anisyl moiety is modified, were designed and synthesised with the purpose of optimising DNA binding, improving cellular uptake/nuclear penetration, and enhancing the modulation of the topoisomerase IIα (TOP2A) gene expression. The modifications include a fluorophenyl group and other heterocycles bearing different molecular shapes, size, and polarity. Like their parent compound HxIP 3, all five X-HxIP analogues bind preferentially to their cognate sequence 5'-TACGAT-3', which is found embedded on the 5' flank of the inverted CCAAT box-2 (ICB2) site in the TOP2A gene promoter, and inhibit protein complex binding. Interestingly, the 4-pyridyl analog 6a exhibits greater binding affinity for the target DNA sequence and abolishes the protein:ICB2 interaction in vitro, at a lower concentration, compared to the prototypical compound HxIP 3. Analogues 6b-e, display improved DNA sequence specificity, but reduced binding affinity for the cognate sequence, relative to the unmodified HxIP 3, with polyamides 6b and 6e being the most sequence selective. However, unlike 3 and 6b, 6a was unable to enter cells, access the nucleus and thereby affect TOP2A gene expression in confluent human lung cancer cells. These results show that while DNA binding affinity and sequence selectivity are important, consideration of cellular uptake and concentration in the nucleus are critical when exerting biological activity is the desired outcome. By characterising the DNA binding, cellular uptake and gene regulatory properties of these small molecules, we can elucidate the determinants of the elicited biological activity, which can be impacted by even small structural modifications in the polyamide molecular design., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

Background: In Europe, golden jackals (Canis aureus) have been expanding their range out of the southern and southeastern Balkans towards central Europe continually since the 1960s. Here, we investigated the level of functional diversity at the MHC class II DLA-DQA1 exon 2 in golden jackal populations from Bulgaria, Serbia, and Hungary. Specifically, we tested for positive selection on and geographic variation at that locus due to adaptation to supposedly regionally varying pathogenic landscapes. To test for potential fitness effects of different protein variants on individual body condition, we used linear modeling of individual body mass indexes (bmi) and accounted for possible age, sex, geographical, and climatic effects. The latter approach was performed, however, only on Serbian individuals with appropriate data., Results: Only three different DLA-DQA1 alleles were detected, all coding for different amino-acid sequences. The neutrality tests revealed no significant but positive values; there was no signal of spatial structuring and no deviation from the Hardy-Weinberg equilibrium across the studied range of expansion. However, we found a signal of trans-species polymorphism and significant test results for positive selection on three codons. Our information-theory based linear modeling results indicated an effect of ambient temperature on the occurrence of individual DLA-DQA1 genotypes in individuals from across the studied expansion range, independent from geographical position. Our linear modeling results of individual bmi values indicated that yearlings homozygous for DLA-DQA1*03001 reached values typical for adults contrary to yearlings carrying other genotypes (protein combinations). This suggested better growth rates and thus a possible fitness advantage of yearlings homozygous for DLA-DQA1*03001., Conclusions: Our results indicate a demographic (stochastic) signal of reduced DLA-DQA1 exon 2 variation, in line with the documented historical demographic bottleneck. At the same time, however, allelic variation was also affected by positive selection and adaptation to varying ambient temperature, supposedly reflecting geographic variation in the pathogenic landscape. Moreover, an allele effect on body mass index values of yearlings suggested differential fitness associated with growth rates. Overall, a combination of a stochastic effect and positive selection has shaped and is still shaping the variation at the studied MHC locus.

Background: Animal mitochondria play a central role in energy production in the cells through the oxidative phosphorylation (OXPHOS) pathway. Recent studies of selection on different mitochondrial OXPHOS genes have revealed the adaptive implications of amino acid changes in these subunits. In hares, climatic variation and/or introgression were suggested to be at the origin of such adaptation. Here we looked for evidence of positive selection in three mitochondrial OXPHOS genes, using tests of selection, protein structure modelling and effects of amino acid substitutions on the protein function and stability. We also used statistical models to test for climate and introgression effects on sites under positive selection., Results: Our results revealed seven sites under positive selection in ND4 and three sites in Cytb. However, no sites under positive selection were observed in the COX1 gene. All three subunits presented a high number of codons under negative selection. Sites under positive selection were mapped on the tridimensional structure of the predicted models for the respective mitochondrial subunit. Of the ten amino acid replacements inferred to have evolved under positive selection for both subunits, six were located in the transmembrane domain. On the other hand, three codons were identified as sites lining proton translocation channels. Furthermore, four codons were identified as destabilizing with a significant variation of Δ vibrational entropy energy between wild and mutant type. Moreover, our PROVEAN analysis suggested that among all positively selected sites two fixed amino acid replacements altered the protein functioning. Our statistical models indicated significant effects of climate on the presence of ND4 and Cytb protein variants, but no effect by trans-specific mitochondrial DNA introgression, which is not uncommon in a number of hare species., Conclusions: Positive selection was observed in several codons in two OXPHOS genes. We found that substitutions in the positively selected codons have structural and functional impacts on the encoded proteins. Our results are concordantly suggesting that adaptations have strongly affected the evolution of mtDNA of hare species with potential effects on the protein function. Environmental/climatic changes appear to be a major trigger of this adaptation, whereas trans-specific introgressive hybridization seems to play no major role for the occurrence of protein variants.

Background: Recent human transcriptomic analyses revealed a very large number of testis-enriched genes, many of which are involved in spermatogenesis. This comprehensive transcriptomic data lead us to the question whether positive selection was a decisive force influencing the evolution and variability of testis-enriched genes in humans. We used two methodological approaches to detect different levels of positive selection, namely episodic positive diversifying selection (i.e., past selection) in the human lineage within primate phylogeny, potentially driven by sperm competition, and recent positive directional selection in contemporary human populations, which would indicate adaptation to different environments., Results: In the human lineage (after correction for multiple testing) we found that only the gene TULP2, for which no functional data are yet available, is subject to episodic positive diversifying selection. Using less stringent statistical criteria (uncorrected p-values), also the gene SPATA16, which has a pivotal role in male fertility and for which episodes of adaptive evolution have been suggested, also displays a putative signal of diversifying selection in the human branch. At the same time, we found evidence for recent positive directional selection acting on several human testis-enriched genes (MORC1, SLC9B1, ROPN1L, DMRT1, PLCZ1, RNF17, FAM71D and WBP2NL) that play important roles in human spermatogenesis and fertilization. Most of these genes are population-specifically under positive selection., Conclusion: Episodic diversifying selection, possibly driven by sperm competition, was not an important force driving the evolution of testis-enriched genes in the human lineage. Population-specific, recent positive directional selection suggests an adaptation of male reproductive genes to different environmental conditions. Positive selection acts on eQTLS and sQTLs, indicating selective effects on important gene regulatory functions. In particular, the transcriptional diversity regulated by sQTLs in testis-enriched genes may be important for spermatocytes to respond to environmental and physiological stress.

Background: Recent 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., Results: We 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., Conclusions: positive 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.

Background: The evolutionary highly conserved neurohypophyseal hormones oxytocin and arginine vasopressin play key roles in regulating social cognition and behaviours. The effects of these two peptides are meditated by their specific receptors, which are encoded by the oxytocin receptor (OXTR) and arginine vasopressin receptor 1a genes (AVPR1A), respectively. In several species, polymorphisms in these genes have been linked to various behavioural traits. Little, however, is known about whether positive selection acts on sequence variants in genes influencing variation in human behaviours., Results: We identified, in both neuroreceptor genes, signatures of balancing selection in the cis-regulative acting sequences such as transcription factor binding and enhancer sequences, as well as in a transcriptional repressor sequence motif. Additionally, in the intron 3 of the OXTR gene, the SNP rs59190448 appears to be under positive directional selection. For rs59190448, only one phenotypical association is known so far, but it is in high LD' (>0.8) with loci of known association; i.e., variants associated with key pro-social behaviours and mental disorders in humans., Conclusions: Only for one SNP on the OXTR gene (rs59190448) was a sign of positive directional selection detected with all three methods of selection detection. For rs59190448, however, only one phenotypical association is known, but rs59190448 is in high LD' (>0.8), with variants associated with important pro-social behaviours and mental disorders in humans. We also detected various signatures of balancing selection on both neuroreceptor genes.

In lagomorphs, lymphocyte subset distributions and the importance of CD4(+) T cell levels has so far only been considered in the frame of rabbit disease models. In this study, the first assessment of CD4(+) T lymphocytes in peripheral blood cells in brown hares (Lepus europaeus L., 1758), a further leporid species using a cross-reactive rabbit anti-CD4 antibody in flow cytometry, is presented. In addition, the entire coding region of the hare CD4 gene (1380 bp) coding for a polypeptide of 459 amino acids has been sequenced. Using generalized least squares fitting by maximum likelihood (GLS) test, significantly (p=0.0095) higher CD4(+) T cell frequencies in males than in females and significantly (p=0.0001) higher frequencies for leverets (younger than 2 months of age) than for subadult and adult (older than 7 months of age) individuals were detected. No significant age influence, however, was found for subadult and adult hares. The study is particularly meant to provide a first step in establishing a toolbox for the assessment of the immune response in this leporid species., (Copyright © 2014 Elsevier B.V. All rights reserved.)

Facial asymmetries are commonly used as a proxy for human developmental imprecision resulting from inbreeding, and thus reduced genetic heterozygosity. Several environmental factors influence human facial asymmetry (e.g., health care, parasites), but the generalizability of findings on genetic stressors has been limited in humans by sample characteristics (island populations, endogamy) and indirect genetic assessment (inference from pedigrees). In a sample of 3215 adult humans from the Rotterdam Study, we therefore studied the relationship of facial asymmetry, estimated from nine mid-facial landmarks, with genetic variation at 102 single nucleotide polymorphism (SNP) loci recently associated with facial shape variation. We further tested whether the degree of individual heterozygosity is negatively correlated with facial asymmetry. An ANOVA tree regression did not identify any SNP relating to either fluctuating asymmetry or total asymmetry. In a general linear model, only age and sex--but neither heterozygosity nor any SNP previously reported to covary with facial shape--was significantly related to total or fluctuating asymmetry of the midface. Our study does not corroborate the common assumption in evolutionary and behavioral biology that morphological asymmetries reflect heterozygosity. Our results, however, may be affected by a relatively small degree of inbreeding, a relatively stable environment, and an advanced age in the Rotterdam sample. Further large-scale genetic studies, including gene expression studies, are necessary to validate the genetic and developmental origin of morphological asymmetries.

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., (Copyright © 2013 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

The genes of the major histocompatibility complex (MHC) are attractive candidates for investigating the link between adaptive variation and individual fitness. High levels of diversity at the MHC are thought to be the result of parasite-mediated selection and there is growing evidence to support this theory. Most studies, however, target just a single gene within the MHC and infer any evidence of selection to be representative of the entire gene region. Here we present data from three MHC class II beta genes (DPB, DQB, and DRB) for brown hares in two geographic regions and compare them against previous results from a class II alpha-chain gene (DQA). We report moderate levels of diversity and high levels of population differentiation in the DQB and DRB genes (Na  =  11, D (est) =  0.071 and Na  = 15, D (est) =  0.409, respectively), but not for the DPB gene (Na  =  4, D (est) =  0.00). We also detected evidence of positive selection within the peptide binding region of the DQB and DRB genes (95% CI, ω  >  1.0) but found no signature of selection for DPB. Mutation and recombination were both found to be important processes shaping the evolution of the class II genes. Our findings suggest that while diversifying selection is a significant contributor to the generally high levels of MHC diversity, it does not act in a uniform manner across the entire MHC class II region. The beta-chain genes that we have characterized provide a valuable set of MHC class II markers for future studies of the evolution of adaptive variation in Leporids.

The link between adaptive genetic variation, individual fitness and wildlife population dynamics is fundamental to the study of ecology and evolutionary biology. In this study, a Bayesian modelling approach was employed to examine whether individual variability at two major histocompatibility complex (MHC) class II loci (DQA and DRB) and eight neutral microsatellite loci explained variation in female reproductive success for wild populations of European brown hare (Lepus europaeus). We examined two aspects of reproduction: the ability to reproduce (sterility) and the number of offspring produced (fecundity). Samples were collected from eastern Austria, experiencing a sub-continental climatic regime, and from Belgium with a more Atlantic-influenced climate. As expected, reproductive success (both sterility and fecundity) was significantly influenced by age regardless of sampling locality. For Belgium, there was also a significant effect of DQA heterozygosity in determining whether females were able to reproduce (95% highest posterior density interval of the regression parameter [-3.64, -0.52]), but no corresponding effect was found for Austria. In neither region was reproduction significantly associated with heterozygosity at the DRB locus. DQA heterozygotes from both regions also showed a clear tendency, but not significantly so, to produce a larger number of offspring. Predictive simulations showed that, in Belgium, sub-populations of homozygotes will have higher rates of sterile individuals and lower average offspring numbers than heterozygotes. No similar effect is predicted for Austria. The mechanism for the spatial MHC effect is likely to be connected to mate choice for increased heterozygosity or to the linkage of certain MHC alleles with lethal recessives at other loci., (© 2010 Blackwell Publishing Ltd.)

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 Histocompatibility 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.

We surveyed the genetic diversity of the expressed major histocompatibility complex class II DQA locus in natural populations of European brown hares, Lepus europaeus, from Austria and Belgium (267 individuals in total). Based on cDNA sequences, we designed hare-specific primers to amplify the highly variable second exon of the DQA gene. Using cloning-sequencing methodology and capillary electrophoresis single-strand conformation polymorphism, we found ten alleles of the DQA exon 2 locus across these two European regions, of which eight are described for the first time. To search for signals of selection and recombination in the evolution of the DQA gene within the leporids, we augmented our sample with orthologous DQA alleles from the European rabbit, Oryctolagus cuniculus, in order to carry out a species level, species pairwise comparison. We found evidence of recombination in the history of the DQA sequences in leporids with some recombinant alleles bridging the species divide. In both species, selection on peptide binding site codons can be detected, though stronger for the rabbit. This result suggests that there may be a differential selection pressure in the deeper evolutionary history of these two species due to differences in several demographic and ecological traits likely subjecting them to differential selection by parasites. Finally, evolutionary relationships show a widespread and statistically significant intermingling of alleles from the two species. The many macroparasites shared between hares and rabbits may explain this pattern of trans-species polymorphism.

Genes of the major histocompatibility complex (MHC) encode molecules that control immune recognition and are highly polymorphic in most vertebrates. The remarkable polymorphisms at MHC loci may be maintained by selection from parasites, sexual selection, or both. If asexual species show equal (or higher) levels of polymorphisms at MHC loci as sexual ones, this would mean that sexual selection is not necessary to explain the high levels of diversity at MHC loci. In this study, we surveyed the MHC diversity of the asexual amazon molly (Poecilia formosa) and one of its sexual ancestors, the sailfin molly (P. latipinna), which lives in the same habitat. We found that the asexual molly has polymorphic MHC loci despite its clonal reproduction, yet not as polymorphic as the sexual species. Although the nucleotide diversity was similar between the asexual and sexual species, the sexual species exhibited a greater genotypic diversity compared to the asexual one from the same habitats. Within-genome diversity was similar for MHC class I loci, but for class IIB, the sexual species had higher diversity compared to the asexual--despite the hybrid origins and higher levels of heterozygosity at microsatellite loci in the asexual species. The level of positive selection appears to be similar between the two species, which suggests that these polymorphisms are maintained by selection. Thus, our findings do not allow us to rule out the sexual selection hypothesis for the evolution of MHC diversity, and although the sexual fish has higher levels of MHC-diversity compared to the asexual species, this may be due to differences in demography, parasites, or other factors, rather than sexual selection.

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.

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.

We examined the evolutionary processes contributing to genetic diversity at the major histocompatibility complex (MHC) class II DRB locus in chamois (Rupicapra spp., subfamily Caprinae). We characterised the pattern of intragenic recombination (or homologous gene conversion) and quantified the amount of recombination in the genealogical history of the two chamois species, Pyrenean chamois (Rupicapra pyrenaica) and Alpine chamois (Rupicapra rupicapra). We found evidence for intragenic recombination, and the estimated amount of population recombination suggests that recombination has been a significant process in generating DRB allelic diversity in the genealogical history of the genus Rupicapra. Moreover, positive selection appears to act on the same peptide-binding residues in both analysed chamois species, but not in identical intensity. Recombination coupled with positive selection drives the rapid evolution at the peptide-binding sites in the MHC class II DRB gene. Many chamois MHC class II DRB alleles are thus much younger than previously assumed.

Overdominant selection or heterozygote advantage can partly explain the extraordinary polymorphism found at classical major histocompatibility complex (MHC) loci. However, several studies employing only single infectious agents often failed to detect it. Here, we review recent studies suggesting that due to the dominant nature of MHC-mediated resistance, a heterozygote advantage is most likely to be detected in multiple pathogen challenges.

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.

This study characterises the expressed MHC class II DRB gene and its genetic variation in exon 2 among 59 Alpine chamois (Rupicapra r. rupicapra) from several locations in the Eastern Alps plus one population in the Western Alps. The expressed DRB gene contains all the expected characteristics including all of the conserved residues found in other mammalian species. The isolated expressed DRB gene from chamois can be considered orthologous to MHC class II DRB genes in other mammalian species. Cattle primers yielded exon 2 sequences of 284 bp length. Sequence analysis revealed 19 DRB exon 2 alleles among 59 chamois. Alleles exhibited a high degree of nucleotide and amino acid polymorphism, with most amino acid variations occurring at positions forming the peptide-binding sites. A significantly higher rate of non-synonymous (dN = 0.053+/-0.0132) than synonymous (dS = 0.007+/-0.003) substitutions indicated positive selection for diversity in the DRB locus.