Detecting the signature of selection on immune genes in highly structured populations of wild sheep (Ovis dalli)

Polymorphism of Selected Regions of Ovar-MHC and the Health Status of the Ovine Mammary Gland

Udder diseases (mastitis) are a serious cause of economic losses in sheep breeding as they have a negative impact on lamb rearing and the quality of dairy products. Thus far, progress in treatment and prevention of these diseases has been insufficient-giving ground for searching possibilities of using natural immunity to combat mastitis. This study aims to assess the relationship between the microsatellite polymorphism of selected Ovar-MHC genes and the health status of the mammary gland of sheep. The research is carried out on sheep of the Polish Heath and Polish Lowland breeds. In ovine milk, the SCC and the percentage of the lymphocyte subpopulation are assessed. Based on genomic DNA, molecular analysis of the Ovar-MHC gene fragments (OLADRB1, OLADRB2, OMHC1) polymorphism is performed. Significant differences are found in SCC level and the percentage of lymphocytes (CD4, CD8, CD19) in the milk depending on the alleles of the Ovar-MHC genes. Alleles of 488 bp (DRB1) and 284 bp (DRB2) are found more frequently in sheep, which milk contained <200 × 10³/mL SCC, while in carriers of the 508 bp (DRB1) and 272 bp (DRB2) alleles, SCC level in milk is significantly higher (>200 × 10³/mL). The obtained results justify the need for further research to better understand the genetic basis of mastitis, and to search for effective molecular markers that can be used in breeding practice.

Allometric relationships in morphological traits associated with foraging, swimming ability, and predator defense reveal adaptations toward brackish and freshwater environments in the threespine stickleback

Freshwater colonization by threespine stickleback has led to divergence in morphology between ancestral marine and derived freshwater populations, making them ideal for studying natural selection on phenotypes. In an open brackish-freshwater system, we previously discovered two genetically distinct stickleback populations that also differ in geometric shape: one mainly found in the brackish water lagoon and one throughout the freshwater system. As shape and size are not perfectly correlated, the aim of this study was to identify the morphological trait(s) that separated the populations in geometric shape. We measured 23 phenotypes likely to be important for foraging, swimming capacity, and defense against predation. The lateral plate morphs in freshwater displayed few significant changes in trait sizes, but the low plated expressed feeding traits more associated with benthic habitats. When comparing the completely plated genetically assigned populations, the freshwater, the hybrids, the migrants and the lagoon fish, many of the linear traits had different slopes and intercepts in trait-size regressions, precluding our ability to directly compare all traits simultaneously, which most likely results from low variation in body length for the lagoon and migrant population. We found the lagoon stickleback population to be more specialized toward the littoral zone, displaying benthic traits such as large, deep bodies with smaller eyes compared to the freshwater completely plated morph. Further, the lagoon and migrant fish had an overall higher body coverage of lateral plates compared to freshwater fish, and the dorsal and pelvic spines were longer. Evolutionary constraints due to allometric scaling relationships could explain the observed, overall restricted, differences in morphology between the sticklebacks in this study, as most traits have diversified in common allometric trajectories. The observed differences in foraging and antipredation traits between the fish with a lagoon and freshwater genetic signature are likely a result of genetic or plastic adaptations toward brackish and freshwater environments.

Heritability of Horn Size in Thinhorn Sheep

Understanding the genetic basis of fitness-related trait variation has long been of great interest to evolutionary biologists. Secondary sexual characteristics, such as horns in bovids, are particularly intriguing since they can be potentially affected by both natural and sexual selection. Until recently, however, the study of fitness-related quantitative trait variation in wild species has been hampered by a lack of genomic resources, pedigree, and/or phenotype data. Recent innovations in genomic technologies have enabled wildlife researchers to perform marker-based relatedness estimation and acquire adequate loci density, enabling both the “top-down” approach of quantitative genetics and the “bottom-up” approach of association studies to describe the genetic basis of fitness-related traits. Here we combine a cross species application of the OvineHD BeadChip and horn measurements (horn length, base circumference, and volume) from harvested thinhorn sheep to examine the heritability and to perform a genome-wide single-nucleotide polymorphism association study of horn size in the species. Thinhorn sheep are mountain ungulates that reside in the mountainous regions of northwestern North America. Thinhorn sheep males grow massive horns that determine the social rank and mating success. We found horn length, base circumference, and volume to be moderately heritable and two loci to be suggestively associated with horn length.

Developmental timing differences underlie armor loss across threespine stickleback populations

Comparing ontogenetic patterns within a well-described evolutionary context aids in inferring mechanisms of change, including heterochronies or deletion of developmental pathways. Because selection acts on phenotypes throughout ontogeny, any within-taxon developmental variation has implications for evolvability. We compare ontogenetic order and timing of locomotion and defensive traits in three populations of threespine stickleback that have evolutionarily divergent adult forms. This analysis adds to the growing understanding of developmental genetic mechanisms of adaptive change in this evolutionary model species by delineating when chondrogenesis and osteogenesis in two derived populations begin to deviate from the developmental pattern in their immediate ancestors. We found that differences in adult defensive morphologies arise through abolished or delayed initiation of these traits rather than via an overall heterochronic shift, that intra-population ontogenetic variation is increased for some derived traits, and that altered armor developmental timing differentiates the derived populations from each other despite parallels in adult lateral plate armor phenotypes. We found that changes in ossified elements of the pelvic armor are linked to delayed and incomplete development of an early-forming pelvic cartilage, and that this disruption likely presages the variable pelvic vestiges documented in many derived populations.

SNP discovery in candidate adaptive genes using exon capture in a free-ranging alpine ungulate

Identification of genes underlying genomic signatures of natural selection is key to understanding adaptation to local conditions. We used targeted resequencing to identify SNP markers in 5321 candidate adaptive genes associated with known immunological, metabolic and growth functions in ovids and other ungulates. We selectively targeted 8161 exons in protein-coding and nearby 5' and 3' untranslated regions of chosen candidate genes. Targeted sequences were taken from bighorn sheep (Ovis canadensis) exon capture data and directly from the domestic sheep genome (Ovis aries v. 3; oviAri3). The bighorn sheep sequences used in the Dall's sheep (Ovis dalli dalli) exon capture aligned to 2350 genes on the oviAri3 genome with an average of 2 exons each. We developed a microfluidic qPCR-based SNP chip to genotype 476 Dall's sheep from locations across their range and test for patterns of selection. Using multiple corroborating approaches (lositan and bayescan), we detected 28 SNP loci potentially under selection. We additionally identified candidate loci significantly associated with latitude, longitude, precipitation and temperature, suggesting local environmental adaptation. The three methods demonstrated consistent support for natural selection on nine genes with immune and disease-regulating functions (e.g. Ovar-DRA, APC, BATF2, MAGEB18), cell regulation signalling pathways (e.g. KRIT1, PI3K, ORRC3), and respiratory health (CYSLTR1). Characterizing adaptive allele distributions from novel genetic techniques will facilitate investigation of the influence of environmental variation on local adaptation of a northern alpine ungulate throughout its range. This research demonstrated the utility of exon capture for gene-targeted SNP discovery and subsequent SNP chip genotyping using low-quality samples in a nonmodel species.

Genetic variation at the MHC DRB1 locus is similar across Gunnison's prairie dog (Cynomys gunnisoni) colonies regardless of plague history

Yersinia pestis was introduced to North America around 1900 and leads to nearly 100% mortality in prairie dog (Cynomys spp.) colonies during epizootic events, which suggests this pathogen may exert a strong selective force. We characterized genetic diversity at an MHC class II locus (DRB1) in Gunnison's prairie dog (C. gunnisoni) and quantified population genetic structure at the DRB1 versus 12 microsatellite loci in three large Arizona colonies. Two colonies, Seligman (SE) and Espee Ranch (ES), have experienced multiple plague‐related die‐offs in recent years, whereas plague has never been documented at Aubrey Valley (AV). We found fairly low allelic diversity at the DRB1 locus, with one allele (DRB1*01) at high frequency (0.67–0.87) in all colonies. Two other DRB1 alleles appear to be trans‐species polymorphisms shared with the black‐tailed prairie dog (C. ludovicianus), indicating that these alleles have been maintained across evolutionary time frames. Estimates of genetic differentiation were generally lower at the MHC locus (F_ST = 0.033) than at microsatellite markers (F_ST = 0.098). The reduced differentiation at DRB1 may indicate that selection has been important for shaping variation at MHC loci, regardless of the presence or absence of plague in recent decades. However, genetic drift has probably also influenced the DRB1 locus because its level of differentiation was not different from that of microsatellites in an F_ST outlier analysis. We then compared specific MHC alleles to plague survivorship in 60 C. gunnisoni that had been experimentally infected with Y. pestis. We found that survival was greater in individuals that carried at least one copy of the most common allele (DRB1*01) compared to those that did not (60% vs. 20%). Although the sample sizes of these two groups were unbalanced, this result suggests the possibility that this MHC class II locus, or a nearby linked gene, could play a role in plague survival.

A comprehensive mapping of the structure and gene organisation in the sheep MHC class I region

Background

The major histocompatibility complex (MHC) is a chromosomal region that regulates immune responsiveness in vertebrates. This region is one of the most important for disease resistance because it has been associated with resistance or susceptibility to a wide variety of diseases and because the MHC often accounts for more of the variance than other loci. Selective breeding for disease resistance is becoming increasingly common in livestock industries, and it is important to determine how this will influence MHC polymorphism and resistance to diseases that are not targeted for selection. However, in sheep the order and sequence of the protein coding genes is controversial. Yet this information is needed to determine precisely how the MHC influences resistance and susceptibility to disease.

Methods

CHORI bacterial artificial chromosomes (BACs) known to contain sequences from the sheep MHC class I region were sub-cloned, and the clones partially sequenced. The resulting sequences were analysed and re-assembled to identify gene content and organisation within each BAC. The low resolution MHC class I physical map was then compared to the cattle reference genome, the Chinese Merino sheep MHC map published by Gao, et al. (2010) and the recently available sheep reference genome.

Results

Immune related class I genes are clustered into 3 blocks; beta, kappa and a novel block not previously identified in other organisms. The revised map is more similar to Bovidae maps than the previous sheep maps and also includes several genes previously not annotated in the Chinese Merino BAC assembly and others not currently annotated in the sheep reference chromosome 20. In particular, the organisation of nonclassical MHC class I genes is similar to that present in the cattle MHC. Sequence analysis and prediction of amino acid sequences of MHC class I classical and nonclassical genes was performed and it was observed that the map contained one classical and eight nonclassical genes together with three possible pseudogenes.

Conclusions

The comprehensive physical map of the sheep MHC class I region enhances our understanding of the genetic architecture of the class I MHC region in sheep and will facilitate future studies of MHC function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1992-4) contains supplementary material, which is available to authorized users.

Detection of Selection Signatures on the X Chromosome in Three Sheep Breeds

Artificial selection has played a critical role in animal breeding. Detection of artificial selection footprints in genomic regions can provide insights for understanding the function of specific phenotypic traits and better guide animal breeding. To more fully understand the relationship between genomic composition and phenotypic diversity arising from breed development, a genome-wide scan was conducted using an OvineSNP50 BeadChip and integrated haplotype score and fixation index analyses to detect selection signatures on the X chromosome in three sheep breeds. We identified 49, 34, and 55 candidate selection regions with lengths of 27.49, 16.47, and 25.42 Mb in German Mutton, Dorper, and Sunit sheep, respectively. Bioinformatics analysis showed that some of the genes in these regions with selection signatures, such as BMP15, were relevant to reproduction. We also identified some selection regions harboring genes that had human orthologs, including BKT, CENPI, GUCY2F, MSN, PCDH11X, PLP1, VSIG4, PAK3, WAS, PCDH19, PDHA1, and SRPX2. The VSIG4 and PCDH11X genes are associated with the immune system and disease, PDHA1 is associated with biosynthetic related pathways, and PCDH19 is expressed in the nervous system and skin. These genes may be useful as candidate genes for molecular breeding.

Molecular evolution of a polymorphic HSP40-like protein encoded in the histocompatibility locus of an invertebrate chordate

Allorecognition, the ability to distinguish self from non-self, occurs in most organisms. Despite the ubiquity of the allorecognition process, the genetic basis for allorecognition remains unexplored in most taxa outside vertebrates and flowering plants. The allorecognition system in the colonial ascidian Botryllus schlosseri is a notable exception. We have recently identified a polymorphic gene within the fuhc locus that may play a role in allorecognition. The encoded protein, called Hsp40-L, is a Type II member of the J-protein family which usually functions as a co-chaperone with Hsp70. While many of the residues that interact with Hsp70 are conserved in Hsp40-L, it may not be a housekeeping protein because it is surprisingly polymorphic and expressed in the ampullae, the site of allorecognition. While the majority of the Hsp40-L protein appears to evolve under purifying selection, a section of the C-terminal region likely experiences balancing/directional selection, characteristic of other allorecognition proteins.

Genomic scan of selective sweeps in thin and fat tail sheep breeds for identifying of candidate regions associated with fat deposition

BACKGROUND

Identification of genomic regions that have been targets of selection for phenotypic traits is one of the most important and challenging areas of research in animal genetics. However, currently there are relatively few genomic regions identified that have been subject to positive selection. In this study, a genome-wide scan using ~50,000 Single Nucleotide Polymorphisms (SNPs) was performed in an attempt to identify genomic regions associated with fat deposition in fat-tail breeds. This trait and its modification are very important in those countries grazing these breeds.

RESULTS

Two independent experiments using either Iranian or Ovine HapMap genotyping data contrasted thin and fat tail breeds. Population differentiation using FST in Iranian thin and fat tail breeds revealed seven genomic regions. Almost all of these regions overlapped with QTLs that had previously been identified as affecting fat and carcass yield traits in beef and dairy cattle. Study of selection sweep signatures using FST in thin and fat tail breeds sampled from the Ovine HapMap project confirmed three of these regions located on Chromosomes 5, 7 and X. We found increased homozygosity in these regions in favour of fat tail breeds on chromosome 5 and X and in favour of thin tail breeds on chromosome 7.

CONCLUSIONS

In this study, we were able to identify three novel regions associated with fat deposition in thin and fat tail sheep breeds. Two of these were associated with an increase of homozygosity in the fat tail breeds which would be consistent with selection for mutations affecting fat tail size several thousand years after domestication.

Evidence for selection at cytokine loci in a natural population of field voles (Microtus agrestis)

Individuals in natural populations are frequently exposed to a wide range of pathogens. Given the diverse profile of gene products involved in responses to different types of pathogen, this potentially results in complex pathogen-specific selection pressures acting on a broad spectrum of immune system genes in wild animals. Thus far, studies into the evolution of immune genes in natural populations have focused almost exclusively on the Major Histocompatibility Complex (MHC). However, the MHC represents only a fraction of the immune system and there is a need to broaden research in wild species to include other immune genes. Here, we examine the evidence for natural selection in a range of non-MHC genes in a natural population of field voles (Microtus agrestis). We concentrate primarily on genes encoding cytokines, signalling molecules critical in eliciting and mediating immune responses and identify signatures of natural selection acting on several of these genes. In particular, genetic diversity within Interleukin 1 beta and Interleukin 2 appears to have been maintained through balancing selection. Taken together with previous findings that polymorphism within these genes is associated with variation in resistance to multiple pathogens, this suggests that pathogen-mediated selection may be an important force driving genetic diversity at cytokine loci in voles and other natural populations. These results also suggest that, along with the MHC, preservation of genetic variation within cytokine genes should be a priority for the conservation genetics of threatened wildlife populations.

Characterization of the sheep Complement Factor B gene (CFB)

The Complement Factor B gene (CFB) of the alternative complement pathway has been identified in the sheep Major Histocompatibility Complex (MHC) and its genomic sequence determined. CFB is located approximately 600 bp upstream of the complement C2 gene, contains 18 exons, and manifests the domain signature characteristic of CFB protein. Thirteen single nucleotide polymorphisms were identified in merino sheep and interbreed variation was identified by comparison with International Sheep Genomics Consortium data. Two predicted non synonymous substitutions were observed and in-silico analysis indicates that these are likely to have a destabilizing effect on the protein structure. Sheep and cattle CFB were compared and shown to contain a common nine nucleotide deletion in exon 18 relative to human CFB. Predicted CFB amino acid sequences for these two species contain 761 aa relative to 764 aa in the human orthologue. Sequencing of the cosmid and BAC clones used in this study permitted the relative positions of three adjacent loci to be determined and showed that the previously described microsatellite locus (BfMs) is located within SKIV2L.

A single nomenclature and associated database for alleles at the major histocompatibility complex class II DRB1 locus of sheep

The development of standardised nomenclatures with associated databases containing reference sequences for alleles at polymorphic loci within the major histocompatibility complex (MHC) has been facilitated by the development of the immuno polymorphism database (IPD). Recently, included within IPD-MHC is information on allelic diversity within sheep species (IPD-MHC-OLA). Here, we present the first report of progress in populating the sheep IPD-MHC database with alleles at the class II MHC DRB1 locus. The sequence of 63 Ovar-DRB1 alleles within 24 allelic families is now held within the database, each meeting the minimum requirement of a complete second exon. These sequences are derived from a combination of genomic and cDNA-based approaches and represent the most extensive collection of validated alleles at the sheep DRB1 locus yet described. Although these 63 alleles probably represent only a fraction of the DRB1 allelic diversity in sheep species worldwide, we encourage the research community to use the official allelic nomenclature and to contribute allelic sequences to the database via its web-based submission tool. In time, the IPD-MHC-OLA resource will underpin population-based MHC genotyping studies and help to simplify meta-analyses of multi-source data from wild and domestic sheep populations.

Association of microsatellite polymorphisms with immune responses to a killed Mycobacterium avium subsp. paratuberculosis vaccine in Merino sheep

AIM

To study the association of polymorphisms at five microsatellite loci with immune responses to a killed Mycobacterium avium subsp. paratuberculosis (Map) vaccine.

METHODS

Merino sheep (504 vaccinates and 430 unvaccinated controls) from a long-term Johne's vaccine trial undertaken on three different properties in the Central Tablelands of New South Wales, Australia, were genotyped for five microsatellite markers located in three immunologically significant chromosome regions. The marker loci included three from the major histocompatibility complex (MHC), namely DYMS1, OLADRB and SMHCC1; and one each from the solute carrier family 11 member 1 (SLC11A1), OVINRA1, and the interferon-γ (IFN-γ), o(IFN)-γ, gene regions. Associations between immune responses and genetic polymorphisms at the marker loci were examined by analysing both allelic and genotypic effects.

RESULTS

The o(IFN)-γ locus had only two alleles, whereas the other four loci exhibited extensive polymorphism, with the number of alleles ranging from 10 (OVINRA1) to 21 (DYMS1), resulting in 30-92 genotypes per locus. Heterozygosities varied between 37% (o(IFN)-γ) and 87% (SMHCC1), while information on polymorphic contents ranged from 0.31 (o(IFN)-γ) to 0.87 (DYMS1). Each of the three properties exhibited unique allelic and genotypic frequencies. Analysis of immune response data revealed strong antibody and IFN-γ responses as early as 2 months post-vaccination. Immune responses in control animals on all three properties remained consistently low, except for slightly elevated IFN-γ responses at a few time-points on two properties, concomitant with exposure to natural infection. Genotype-phenotype association analyses revealed a number of marker genotypes/alleles to be significantly associated with antibody and IFN-γ responses. However, the effects of only five genotypes (one each at DYMS1, OLADRB, SMHCC1, OVINRA1 and o(IFN)-γ) and three alleles (one each at o(IFN)-γ, DYMS1and OLADRB) on IFN-γ responses were consistent across the three properties.

CONCLUSION

Considering the significance of IFN-γ responses in protection against Map, it is possible that the genotypes/alleles identified might have a role in protective immune responses to natural Map infections, and further studies are warranted to confirm this.

Identifying footprints of directional and balancing selection in marine and freshwater three-spined stickleback (Gasterosteus aculeatus) populations

Natural selection is expected to leave an imprint on the neutral polymorphisms at the adjacent genomic regions of a selected gene. While directional selection tends to reduce within-population genetic diversity and increase among-population differentiation, the reverse is expected under balancing selection. To identify targets of natural selection in the three-spined stickleback (Gasterosteus aculeatus) genome, 103 microsatellite and two indel markers including expressed sequence tags (EST) and quantitative trait loci (QTL)-associated loci, were genotyped in four freshwater and three marine populations. The results indicated that a high proportion of loci (14.7%) might be affected by balancing selection and a lower proportion (2.8%) by directional selection. The strongest signatures of directional selection were detected in a microsatellite locus and two indel markers located in the intronic regions of the Eda-gene coding for the number of lateral plates. Yet, other microsatellite loci previously found to be informative in QTL-mapping studies revealed no signatures of selection. Two novel microsatellite loci (Stn12 and Stn90) located in chromosomes I and VIII, respectively, showed signals of directional selection and might be linked to genomic regions containing gene(s) important for adaptive divergence. Although the coverage of the total genomic content was relatively low, the predominance of balancing selection signals is in agreement with the contention that balancing, rather than directional selection is the predominant mode of selection in the wild.

The genetic architecture of susceptibility to parasites

Background

The antagonistic co-evolution of hosts and their parasites is considered to be a potential driving force in maintaining host genetic variation including sexual reproduction and recombination. The examination of this hypothesis calls for information about the genetic basis of host-parasite interactions – such as how many genes are involved, how big an effect these genes have and whether there is epistasis between loci. We here examine the genetic architecture of quantitative resistance in animal and plant hosts by concatenating published studies that have identified quantitative trait loci (QTL) for host resistance in animals and plants.

Results

Collectively, these studies show that host resistance is affected by few loci. We particularly show that additional epistatic interactions, especially between loci on different chromosomes, explain a majority of the effects. Furthermore, we find that when experiments are repeated using different host or parasite genotypes under otherwise identical conditions, the underlying genetic architecture of host resistance can vary dramatically – that is, involves different QTLs and epistatic interactions. QTLs and epistatic loci vary much less when host and parasite types remain the same but experiments are repeated in different environments.

Conclusion

This pattern of variability of the genetic architecture is predicted by strong interactions between genotypes and corroborates the prevalence of varying host-parasite combinations over varying environmental conditions. Moreover, epistasis is a major determinant of phenotypic variance for host resistance. Because epistasis seems to occur predominantly between, rather than within, chromosomes, segregation and chromosome number rather than recombination via cross-over should be the major elements affecting adaptive change in host resistance.

Density-related changes in selection pattern for major histocompatibility complex genes in fluctuating populations of voles

Host-pathogen interactions are of particular interest in studies of the interplay between population dynamics and natural selection. The major histocompatibility complex (MHC) genes of demographically fluctuating species are highly suitable markers for such studies, because they are involved in initiating the immune response against pathogens and display a high level of adaptive genetic variation. We investigated whether two MHC class II genes (DQA1, DRB) were subjected to contemporary selection during increases in the density of fossorial water vole (Arvicola terrestris) populations, by comparing the neutral genetic structure of seven populations with that estimated from MHC genes. Tests for heterozygosity excess indicated that DQA1 was subject to intense balancing selection. No such selection operated on neutral markers. This pattern of selection became more marked with increasing abundance. In the low-abundance phase, when populations were geographically isolated, both overall differentiation and isolation-by-distance were more marked for MHC genes than for neutral markers. Model-based simulations identified DQA1 as an outlier (i.e. under selection) in a single population, suggesting the action of local selection in fragmented populations. The differences between MHC and neutral markers gradually disappeared with increasing effective migration between sites. In the high-abundance year, DQA1 displayed significantly lower levels of overall differentiation than the neutral markers. This gene therefore displayed stronger homogenization than observed under drift and migration alone. The observed signs of selection were much weaker for DRB. Spatial and temporal fluctuations in parasite pressure and locus-specific selection are probably the most plausible mechanisms underlying the observed changes in selection pattern during the demographic cycle.

Selection on MHC-linked microsatellite loci in sheep populations

Microsatellites within the major histocompatibility complex (MHC) region have received increasing attention as proxy measures of the level of polymorphism at the Mhc genes themselves. We assessed the diversity of microsatellite loci within or in close proximity of the Mhc genes in several breeds of domestic sheep (Ovis aries) and the wild Mouflon (Ovis orientalis musimon). This was compared to variation at other microsatellite loci scattered throughout the sheep genome. Significantly higher number of alleles were observed at the MHC microsatellites. The sheep breeds studied fell into high- and low-diversity group. This grouping is not related to the agricultural use of the breeds, whether for milk, meat or wool. It is, however, correlated with the geographic origins of the breeds. Southern breeds are genetically more diverse than northern breeds. The observed heterozygosity was in most cases lower than Hardy-Weinberg expectations. The potential impact of selective breeding by man on this is discussed. Neutrality tests indicated that for most of the breeds, the distribution of alleles at the MHC-linked microsatellites are more even than would be expected if the genes were neutral and sampled from populations under drift-mutation equilibrium. Hitchhiking due to tight linkage with alleles at the MHC loci that are under balancing selection is proposed as a possible explanation for this pattern.