Ancient polymorphism and functional variation in the primate MHC-DQA1 5' cis-regulatory region

Comparative three-dimensional genome architectures of adipose tissues provide insight into human-specific regulation of metabolic homeostasis

Elucidating the regulatory mechanisms of human adipose tissues (ATs) evolution is essential for understanding human-specific metabolic regulation, but the functional importance and evolutionary dynamics of three-dimensional (3D) genome organizations of ATs are not well defined. Here, we compared the 3D genome architectures of anatomically distinct ATs from humans and six representative mammalian models. We recognized evolutionarily conserved and human-specific chromatin conformation in ATs at multiple scales, including compartmentalization, topologically associating domain (TAD), and promoter-enhancer interactions (PEI), which have not been described previously. We found PEI are much more evolutionarily dynamic with respect to compartmentalization and topologically associating domain. Compared to conserved PEIs, human-specific PEIs are enriched for human-specific sequence, and the binding motifs of their potential mediators (transcription factors) are less conserved. Our data also demonstrated that genes involved in the evolutionary dynamics of chromatin organization have weaker transcriptional conservation than those associated with conserved chromatin organization. Furthermore, the genes involved in energy metabolism and the maintenance of metabolic homeostasis are enriched in human-specific chromatin organization, while housekeeping genes, health-related genes, and genetic variations are enriched in evolutionarily conserved compared to human-specific chromatin organization. Finally, we showed extensively divergent human-specific 3D genome organizations among one subcutaneous and three visceral ATs. Together, these findings provide a global overview of 3D genome architecture dynamics between ATs from human and mammalian models and new insights into understanding the regulatory evolution of human ATs.

HLA-DRB1 genes and the expression dynamics of HLA CIITA determine the susceptibility to T2DM

Type 2 diabetes mellitus (T2DM) is a disease with polygenic inheritance. The expression of major histocompatibility complex class II genes are regulated by several trans-activators. We have studied the expression of HLA-DRB1, RFX, CIITA-P1, PIV transactivators, immunophenotyping of cells, SNPs in CIITA-168 (A/G) and IFN-γ + 874 (T/A) in T2DM patients and controls (n = 201 each). We observed increased frequencies of DRB1*03, DRB1*04 and DRB1*07 and decreased frequencies of DRB1*10, DRB1*14, and DRB1*15 alleles among patients. Significant up-regulations of HLA-DRB1 genes were observed in patients (p < 0.0001). Down-regulated expressions were documented in DRB1*03-homo (p < 0.002) and DRB1*04-homo (p < 0.009) patients. No significant differences were observed for CIITA-P1 expression except DRB1*04-pooled (p < 0.0113). The CIITA-PIV was up-regulated in overall (p < 0.0001), DRB1*03-pooled (p < 0.0006), DRB1*03-hetero (p < 0.0006) and DRB1*03-homo (p < 0.001) T2DM patients. However, significant down-regulations were documented for DRB1*04-pooled (p < 0.040), DRB1*04-hetero (p < 0.060), and DRB1*04-homo (p < 0.027) combinations. Further, significant down-regulations of RFX5 were observed in overall (p < 0.0006), DRB1*04-pooled (p < 0.0022), and DRB1*04-hetero (p < 0.0004) combinations. Immunophenotyping studies revealed significant increase of CD45⁺ CD14^-, CD19⁺, CD14^- and CD8 cells and elevated level of expression of IFN-γ (p < 0.0001) in patients. A significant increase of TT (p < 3.35 × 10^-6) and decrease of TA (p < 4.57 × 10^-4) genotypes of IFN-γ + 874 (T/A) and an increase of GG (p < 0.001) and decrease of AG (p < 8.24 × 10^-5) genotypes of CIITA-168 A/G SNPs were observed. The combinatorial analysis revealed susceptible associations for DRB1*03 + AA, *03 + AG, *03 + GG and *04 + GG and protective associations for DRB1*10 + AG, *10 + GG, *15 + AG, and *14 + GG combinations. Thus, the present study corroborated the effect of differential expressions of promoters of risk alleles in the pathogenesis of T2DM.

Reorganization of 3D genome structure may contribute to gene regulatory evolution in primates

A growing body of evidence supports the notion that variation in gene regulation plays a crucial role in both speciation and adaptation. However, a comprehensive functional understanding of the mechanisms underlying regulatory evolution remains elusive. In primates, one of the crucial missing pieces of information towards a better understanding of regulatory evolution is a comparative annotation of interactions between distal regulatory elements and promoters. Chromatin conformation capture technologies have enabled genome-wide quantifications of such distal 3D interactions. However, relatively little comparative research in primates has been done using such technologies. To address this gap, we used Hi-C to characterize 3D chromatin interactions in induced pluripotent stem cells (iPSCs) from humans and chimpanzees. We also used RNA-seq to collect gene expression data from the same lines. We generally observed that lower-order, pairwise 3D genomic interactions are conserved in humans and chimpanzees, but higher order genomic structures, such as topologically associating domains (TADs), are not as conserved. Inter-species differences in 3D genomic interactions are often associated with gene expression differences between the species. To provide additional functional context to our observations, we considered previously published chromatin data from human stem cells. We found that inter-species differences in 3D genomic interactions, which are also associated with gene expression differences between the species, are enriched for both active and repressive marks. Overall, our data demonstrate that, as expected, an understanding of 3D genome reorganization is key to explaining regulatory evolution.

The way in which a genome folds affects the regulation of gene expression. This is often due to loops in the three-dimensional structure that bring linearly distant genes and regulatory elements into close proximity. Most studies examining three-dimensional structure genome-wide are limited to a single species. In this study, we compared three-dimensional structure in the genomes of induced pluripotent stem cells from humans and chimpanzees. We collected gene expression data from the same samples, which allowed us to assess the contribution of three-dimensional chromatin conformation to gene regulatory evolution in primates. Our results demonstrate that gene expression differences between the species may often be mediated by differences in three-dimensional genomic interactions. Our data also suggest that large-scale chromatin structures (i.e. topologically associating domains, TADs) are not well conserved in their placement across species. We hope the analytical paradigms we present here could serve as a basis for future comparative studies of three-dimensional genome organization, elucidating the putative functional regulatory loci driving speciation.

Negative Frequency-Dependent Selection Is Frequently Confounding

Persistent genetic variation within populations presents an evolutionary problem, as natural selection and genetic drift tend to erode genetic diversity. Models of balancing selection were developed to account for the maintenance of genetic variation observed in natural populations. Negative frequency-dependent selection is a powerful type of balancing selection that maintains many natural polymorphisms, but it is also commonly misinterpreted. This review aims to clarify the processes underlying negative frequency-dependent selection, describe classes of polymorphisms that can and cannot result from these processes, and discuss the empirical data needed to accurately identify processes that generate or maintain diversity in nature. Finally, the importance of accurately describing the processes affecting genetic diversity within populations as it relates to research progress is considered.

Cis-regulatory evolution in a wild primate: Infection-associated genetic variation drives differential expression of MHC-DQA1 in vitro

Few studies have combined genetic association analyses with functional characterization of infection-associated SNPs in natural populations of nonhuman primates. Here, we investigate the relationship between host genetic variation, parasitism and natural selection in a population of red colobus (Procolobus rufomitratus tephrosceles) in Kibale National Park, Uganda. We collected parasitological, cellular and genomic data to test the following hypotheses: (i) MHC-DQA1 regulatory genetic variation is associated with control of whipworm (Trichuris) infection in a natural population of red colobus; (ii) infection-associated SNPs are functional in driving differential gene expression in vitro; and (iii) balancing selection has shaped patterns of variation in the MHC-DQA1 promoter. We identified two SNPs in the MHC-DQA1 promoter, both in transcription factor binding sites, and both of which are associated with decreased control of Trichuris infection. We characterized the function of both SNPs by testing differences in gene expression between the two alleles of each SNP in two mammalian cell lines. Alleles of one of the SNPs drove differential gene expression in both cell lines, while the other SNP drove differences in expression in one of the cell lines. Additionally, we found evidence of balancing selection acting on the MHC-DQA1 promoter, including extensive trans-species polymorphisms between red colobus and other primates, and an excess of intermediate-frequency alleles relative to genome-wide, coding and noncoding RADseq data. Our data suggest that balancing selection provides adaptive regulatory flexibility that outweighs the consequences of increased parasite infection intensity in heterozygotes.

An evolutionary medicine perspective on Neandertal extinction

The Eurasian sympatry of Neandertals and anatomically modern humans - beginning at least 45,000 years ago and possibly lasting for more than 5000 years - has sparked immense anthropological interest into the factors that potentially contributed to Neandertal extinction. Among many different hypotheses, the "differential pathogen resistance" extinction model posits that Neandertals were disproportionately affected by exposure to novel infectious diseases that were transmitted during the period of spatiotemporal sympatry with modern humans. Comparisons of new archaic hominin paleogenome sequences with modern human genomes have confirmed a history of genetic admixture - and thus direct contact - between humans and Neandertals. Analyses of these data have also shown that Neandertal nuclear genome genetic diversity was likely considerably lower than that of the Eurasian anatomically modern humans with whom they came into contact, perhaps leaving Neandertal innate immune systems relatively more susceptible to novel pathogens. In this study, we compared levels of genetic diversity in genes for which genetic variation is hypothesized to benefit pathogen defense among Neandertals and African, European, and Asian modern humans, using available exome sequencing data (three individuals, or six chromosomes, per population). We observed that Neandertals had only 31-39% as many nonsynonymous (amino acid changing) polymorphisms across 73 innate immune system genes compared to modern human populations. We also found that Neandertal genetic diversity was relatively low in an unbiased set of balancing selection candidate genes for primates, those genes with the highest 1% genetic diversity genome-wide in non-human hominoids (apes). In contrast, Neandertals had similar or higher levels of genetic diversity than humans in 12 major histocompatibility complex (MHC) genes. Thus, while Neandertals may have been relatively more susceptible to some novel pathogens and differential pathogen resistance could be considered as one potential contributing factor in their extinction, the expectations of this model are not universally met.

Balancing selection on allorecognition genes in the colonial ascidian Botryllus schlosseri

Allorecognition is the capability of an organism to recognize its own or related tissues. The colonial ascidian Botryllus schlosseri, which comprises five genetically distinct and divergent species (Clades A-E), contains two adjacent genes that control allorecognition: fuhc^sec and fuhc^tm. These genes have been characterized extensively in Clade A and are highly polymorphic. Using alleles from 10 populations across the range of Clade A, we investigated the type and strength of selection maintaining this variation. Both fuhc genes exhibit higher within-population variation and lower population differentiation measures (F_ST) than neutral loci. The fuhc genes contain a substantial number of codons with >95% posterior probability of d_N/d_S > 1. fuhc^sec and fuhc^tm also have polymorphisms shared between Clade A and Clade E that were present prior to speciation (trans-species polymorphisms). These results provide robust evidence that the fuhc genes are evolving under balancing selection.

A Visual Interface for Querying Heterogeneous Phylogenetic Databases

Despite the recent growth in the number of phylogenetic databases, access to these wealth of resources remain largely tool or form-based interface driven. It is our thesis that the flexibility afforded by declarative query languages may offer the opportunity to access these repositories in a better way, and to use such a language to pose truly powerful queries in unprecedented ways. In this paper, we propose a substantially enhanced closed visual query language, called PhyQL, that can be used to query phylogenetic databases represented in a canonical form. The canonical representation presented helps capture most phylogenetic tree formats in a convenient way, and is used as the storage model for our PhyloBase database for which PhyQL serves as the query language. We have implemented a visual interface for the end users to pose PhyQL queries using visual icons, and drag and drop operations defined over them. Once a query is posed, the interface translates the visual query into a Datalog query for execution over the canonical database. Responses are returned as hyperlinks to phylogenies that can be viewed in several formats using the tree viewers supported by PhyloBase. Results cached in PhyQL buffer allows secondary querying on the computed results making it a truly powerful querying architecture.

A time transect of exomes from a Native American population before and after European contact

A major factor for the population decline of Native Americans after European contact has been attributed to infectious disease susceptibility. To investigate whether a pre-existing genetic component contributed to this phenomenon, here we analyse 50 exomes of a continuous population from the Northwest Coast of North America, dating from before and after European contact. We model the population collapse after European contact, inferring a 57% reduction in effective population size. We also identify signatures of positive selection on immune-related genes in the ancient but not the modern group, with the strongest signal deriving from the human leucocyte antigen (HLA) gene HLA-DQA1. The modern individuals show a marked frequency decrease in the same alleles, likely due to the environmental change associated with European colonization, whereby negative selection may have acted on the same gene after contact. The evident shift in selection pressures correlates to the regional European-borne epidemics of the 1800s.

A First Nation population declined after European contact, likely as a result of infectious disease. Here, researchers partner with indigenous communities to analyse ancient and modern Native American exomes, and find a shift in selection pressure on immune genes, correlated to European-borne epidemics.

Identification of Multiple Forms of RNA Transcripts Associated with Human-Specific Retrotransposed Gene Copies

The human genome contains thousands of retrocopies, mostly as processed pseudogenes, which were recently shown to be prevalently transcribed. In particular, those specifically acquired in the human lineage are able to modulate gene expression in a manner that contributed to the evolution of human-specific traits. Therefore, knowledge of the human-specific retrocopies that are transcribed or their full-length transcript structure contributes to better understand human genome evolution. In this study, we identified 16 human-specific retrocopies that harbor 5' CpG islands by in silico analysis and showed that 12 were transcribed in normal tissues and cancer cell lines with a variety of expression patterns, including cancer-specific expression. Determination of the structure of the transcripts associated with the retrocopies revealed that none were transcribed from their 5' CpG islands, but rather, from inside the 3' UTR and the nearby 5' flanking region of the retrocopies as well as the promoter of neighboring genes. The multiple forms of the transcripts, such as chimeric and individual transcripts in both the sense and antisense orientation, might have introduced novel post-transcriptional regulation into the genome during human evolution. These results shed light on the potential role of human-specific retrocopies in the evolution of gene regulation and genomic disorders.

Next-generation genotyping of hypervariable loci in many individuals of a non-model species: technical and theoretical implications

BACKGROUND

Across species, diversity at the Major Histocompatibility Complex (MHC) is critical to disease resistance and population health; however, use of MHC diversity to quantify the genetic health of populations has been hampered by the extreme variation found in MHC genes. Next generation sequencing (NGS) technology generates sufficient data to genotype even the most diverse species, but workflows for distinguishing artifacts from alleles are still under development. We used NGS to evaluate the MHC diversity of over 300 captive and wild ring-tailed lemurs (Lemur catta: Primates: Mammalia). We modified a published workflow to address errors that arise from deep sequencing individuals and tested for evidence of selection at the most diverse MHC genes.

RESULTS

In addition to evaluating the accuracy of 454 Titanium and Ion Torrent PGM for genotyping large populations at hypervariable genes, we suggested modifications to improve current methods of allele calling. Using these modifications, we genotyped 302 out of 319 individuals, obtaining an average sequencing depth of over 1000 reads per amplicon. We identified 55 MHC-DRB alleles, 51 of which were previously undescribed, and provide the first sequences of five additional MHC genes: DOA, DOB, DPA, DQA, and DRA. The additional five MHC genes had one or two alleles each with little sequence variation; however, the 55 MHC-DRB alleles showed a high dN/dS ratio and trans-species polymorphism, indicating a history of positive selection. Because each individual possessed 1-7 MHC-DRB alleles, we suggest that ring-tailed lemurs have four, putatively functional, MHC-DRB copies.

CONCLUSIONS

In the future, accurate genotyping methods for NGS data will be critical to assessing genetic variation in non-model species. We recommend that future NGS studies increase the proportion of replicated samples, both within and across platforms, particularly for hypervariable genes like the MHC. Quantifying MHC diversity within non-model species is the first step to assessing the relationship of genetic diversity at functional loci to individual fitness and population viability. Owing to MHC-DRB diversity and copy number, ring-tailed lemurs may serve as an ideal model for estimating the interaction between genetic diversity, fitness, and environment, especially regarding endangered species.

Recent Selection Changes in Human Genes under Long-Term Balancing Selection

Balancing selection is an important evolutionary force that maintains genetic and phenotypic diversity in populations. Most studies in humans have focused on long-standing balancing selection, which persists over long periods of time and is generally shared across populations. But balanced polymorphisms can also promote fast adaptation, especially when the environment changes. To better understand the role of previously balanced alleles in novel adaptations, we analyzed in detail four loci as case examples of this mechanism. These loci show hallmark signatures of long-term balancing selection in African populations, but not in Eurasian populations. The disparity between populations is due to changes in allele frequencies, with intermediate frequency alleles in Africans (likely due to balancing selection) segregating instead at low- or high-derived allele frequency in Eurasia. We explicitly tested the support for different evolutionary models with an approximate Bayesian computation approach and show that the patterns in PKDREJ, SDR39U1, and ZNF473 are best explained by recent changes in selective pressure in certain populations. Specifically, we infer that alleles previously under long-term balancing selection, or alleles linked to them, were recently targeted by positive selection in Eurasian populations. Balancing selection thus likely served as a source of functional alleles that mediated subsequent adaptations to novel environments.

Trans-species polymorphism in humans and the great apes is generally maintained by balancing selection that modulates the host immune response

Known examples of ancient identical-by-descent genetic variants being shared between evolutionarily related species, known as trans-species polymorphisms (TSPs), result from counterbalancing selective forces acting on target genes to confer resistance against infectious agents. To date, putative TSPs between humans and other primate species have been identified for the highly polymorphic major histocompatibility complex (MHC), the histo-blood ABO group, two antiviral genes (ZC3HAV1 and TRIM5), an autoimmunity-related gene LAD1 and several non-coding genomic segments with a putative regulatory role. Although the number of well-characterized TSPs under long-term balancing selection is still very small, these examples are connected by a common thread, namely that they involve genes with key roles in the immune system and, in heterozygosity, appear to confer genetic resistance to pathogens. Here, we review known cases of shared polymorphism that appear to be under long-term balancing selection in humans and the great apes. Although the specific selective agent(s) responsible are still unknown, these TSPs may nevertheless be seen as constituting important adaptive events that have occurred during the evolution of the primate immune system.

Epigenetic modifications are associated with inter-species gene expression variation in primates

BACKGROUND

Changes in gene regulation have long been thought to play an important role in evolution and speciation, especially in primates. Over the past decade, comparative genomic studies have revealed extensive inter-species differences in gene expression levels, yet we know much less about the extent to which regulatory mechanisms differ between species.

RESULTS

To begin addressing this gap, we perform a comparative epigenetic study in primate lymphoblastoid cell lines, to query the contribution of RNA polymerase II and four histone modifications, H3K4me1, H3K4me3, H3K27ac, and H3K27me3, to inter-species variation in gene expression levels. We find that inter-species differences in mark enrichment near transcription start sites are significantly more often associated with inter-species differences in the corresponding gene expression level than expected by chance alone. Interestingly, we also find that first-order interactions among the five marks, as well as chromatin states, do not markedly contribute to the degree of association between the marks and inter-species variation in gene expression levels, suggesting that the marginal effects of the five marks dominate this contribution.

CONCLUSIONS

Our observations suggest that epigenetic modifications are substantially associated with changes in gene expression levels among primates and may represent important molecular mechanisms in primate evolution.

Trans-Species Polymorphism in Immune Genes: General Pattern or MHC-Restricted Phenomenon?

Immunity exhibits extraordinarily high levels of variation. Evolution of the immune system in response to host-pathogen interactions in particular ecological contexts appears to be frequently associated with diversifying selection increasing the genetic variability. Many studies have documented that immunologically relevant polymorphism observed today may be tens of millions years old and may predate the emergence of present species. This pattern can be explained by the concept of trans-species polymorphism (TSP) predicting the maintenance and sharing of favourable functionally important alleles of immune-related genes between species due to ongoing balancing selection. Despite the generality of this concept explaining the long-lasting adaptive variation inherited from ancestors, current research in TSP has vastly focused only on major histocompatibility complex (MHC). In this review we summarise the evidence available on TSP in human and animal immune genes to reveal that TSP is not a MHC-specific evolutionary pattern. Further research should clearly pay more attention to the investigation of TSP in innate immune genes and especially pattern recognition receptors which are promising candidates for this type of evolution. More effort should also be made to distinguish TSP from convergent evolution and adaptive introgression. Identification of balanced TSP variants may represent an accurate approach in evolutionary medicine to recognise disease-resistance alleles.

Advantageous diversity maintained by balancing selection in humans

Most human polymorphisms are neutral or slightly deleterious, but some genetic variation is advantageous and maintained in populations by balancing selection. Considered a rarity and overlooked for years, balanced polymorphisms have recently received renewed attention with several lines of evidence showing their relevance in human evolution. From theoretical work on its role in adaptation to empirical studies that identify its targets, recent developments have showed that balancing selection is more prevalent than previously thought. Here we review these developments and discuss their implications in our understanding of the influence of balancing selection in human evolution. We also review existing evidence on the biological functions that benefit most from advantageous diversity, and the functional consequences of these variants. Overall, we argue that balancing selection must be considered an important selective force in human evolution.

The maintenance of mitochondrial genetic variation by negative frequency-dependent selection

Mitochondrial genes generally show high levels of standing genetic variation, which is puzzling given the accumulating evidence for phenotypic effects of mitochondrial genetic variation. Negative frequency-dependent selection, where the relative fitness of a genotype is inversely related to its frequency in a population, provides a potent and potentially general process that can maintain mitochondrial polymorphism. We assessed the change in mitochondrial haplotype frequencies over 10 generations of experimental evolution in 180 seed beetle populations in the laboratory, where haplotypes competed for propagation to subsequent generations. We found that haplotypes consistently increased in frequency when they were initially rare and decreased in frequency when initially common. Our results have important implications for the use of mtDNA haplotype frequency data to infer population level processes and they revive the general hypothesis that negative frequency-dependent selection, presumably caused by habitat heterogeneity, may commonly promote polymorphism in ecologically relevant life history genes.

Long-standing balancing selection in the THBS4 gene: influence on sex-specific brain expression and gray matter volumes in Alzheimer disease

The THBS4 gene encodes a glycoprotein involved in inflammatory responses and synaptogenesis. THBS4 is expressed at higher levels in the brain of humans compared with nonhuman primates, and the protein accumulates in β-amyloid plaques. We analyzed THBS4 genetic variability in humans and show that two haplotypes (hap1 and hap2) are maintained by balancing selection and modulate THBS4 expression in lymphocytes. Indeed, the balancing selection region covers a predicted transcriptional enhancer. In humans, but not in macaques and chimpanzees, THBS4 brain expression increases with age, and variants in the balancing selection region interact with sex in influencing THBS4 expression (pinteraction = 0.038), with hap1 homozygous females showing lowest expression. In Alzheimer disease (AD) patients, significant interactions between sex and THBS4 genotype were detected for peripheral gray matter (pinteraction = 0.014) and total gray matter (pinteraction = 0.012) volumes. Similarly to the gene expression results, the interaction is mainly mediated by hap1 homozygous AD females, who show reduced volumes. Thus, the balancing selection target in THBS4 is likely represented by one or more variants that regulate tissue-specific and sex-specific gene expression. The selection signature associated with THBS4 might not be related to AD pathogenesis, but rather to inflammatory responses.

The Loci of repeated evolution: a catalog of genetic hotspots of phenotypic variation

What is the nature of the genetic changes underlying phenotypic evolution? We have catalogued 1008 alleles described in the literature that cause phenotypic differences among animals, plants, and yeasts. Surprisingly, evolution of similar traits in distinct lineages often involves mutations in the same gene ("gene reuse"). This compilation yields three important qualitative implications about repeated evolution. First, the apparent evolution of similar traits by gene reuse can be traced back to two alternatives, either several independent causative mutations or a single original mutational event followed by sorting processes. Second, hotspots of evolution-defined as the repeated occurrence of de novo mutations at orthologous loci and causing similar phenotypic variation-are omnipresent in the literature with more than 100 examples covering various levels of analysis, including numerous gain-of-function events. Finally, several alleles of large effect have been shown to result from the aggregation of multiple small-effect mutations at the same hotspot locus, thus reconciling micromutationist theories of adaptation with the empirical observation of large-effect variants. Although data heterogeneity and experimental biases prevented us from extracting quantitative trends, our synthesis highlights the existence of genetic paths of least resistance leading to viable evolutionary change.

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

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

Adaptive diversity of innate immune receptor family short pentraxins in Murinae

The short pentraxins C-reactive protein (CRP) and serum amyloid P component (SAP) constitute a group of innate immune receptors that trigger immune activation by detecting molecules of the microbial cell wall. Here, we examined the evolution of short pentraxins in Murinae lineages. By molecular evolutionary analysis, CRP and SAP have been experiencing rapid diversification, driven by adaptive selection. Further, our protein modeling demonstrates that adaptively selected amino acids lie in the ligand-binding region and contact region between subunits. Our findings suggest that rapid diversification of these regions could contribute to the determinants of recognizing specificity and the interaction between subunits.

An evolutionary analysis of RAC2 identifies haplotypes associated with human autoimmune diseases

The human RAC2 gene encodes a small GTP-binding protein with a pivotal role in immune activation and in the induction of peripheral immune tolerance through restimulation-induced cell death (RICD). Different human pathogens target the protein product of RAC2, suggesting that the gene may be subject to natural selection, and that variants in RAC2 may affect immunological phenotypes in humans. We scanned the genomic region encompassing the entire transcription unit for the presence of putative noncoding regulatory elements conserved across mammals. This information was used to select two RAC2 gene regions and analyze their intraspecific genetic diversity. Results suggest that a region covering the 3' untranslated region has been a target of multiallelic balancing selection (or diversifying selection), and three major RAC2 haplogroups occur in human populations. Haplotypes belonging to one of these clades are associated with increased susceptibility to multiple sclerosis (P = 0.022) and earlier onset of disease symptoms (P = 0.025). This same haplogroup is significantly more common in patients with Crohn's disease compared with healthy controls (P = 0.048). These data reinforce recent evidences that susceptibility alleles/haplotypes are shared among multiple autoimmune disorders and support a causal "role for RAC2" variants in the pathogenesis of autoimmune diseases. Other genes with a role in RICD have previously been associated with autoimmunity in humans, suggesting that this pathway and RAC2 may represent novel therapeutic targets in autoimmune disorders.

HLA-DRB1 may be antagonistically regulated by the coordinately evolved promoter and 3'-UTR under stabilizing selection

HLA-DRB1 is the most polymorphic MHC (major histocompatibility complex) class II gene in human, and plays a crucial role in the development and function of the immune system. Extensive polymorphisms exist in the promoter and 3'-UTR of HLA-DRB1, especially a LTR (Long terminal repeat) element in the promoter, which may be involved in the expression regulation. However, it remains unknown how the polymorphisms in the whole promoter region and 3'-UTR to regulate the gene expression. In this study, we investigated the extensive polymorphisms in the HLA-DRB1 promoter and 3'-UTR, and how these polymorphisms affect the gene expression in both independent and jointly manners. It was observed that most of the haplotypes in the DRB1 promoter and 3'-UTR were clustered into 4 conserved lineages (H1, H2, H3 and H4), and showed high linkage disequilibrium. Compared with H1 and H2 lineage, a LTR element in the promoter of H3 and H4 lineage significantly suppressed the promoter activity, whereas the activity of the linked 3'-UTR increased, leading to no apparent difference in the final expression product between H1/H2 and H3/H4 lineage. Nevertheless, compared with the plasmid with a promoter and 3'-UTR from the same lineage, the recombinant plasmid with a promoter from H2 and a 3'-UTR from H3 showed about double fold increased luciferase activity, Conversely, the recombinant plasmid with a promoter from H3 and a 3'-UTR from H2 resulted in about 2-fold decreased luciferase activity. These results indicate that the promoter and 3'-UTR of HLA-DRB1 may antagonistically regulate the gene expression, which may be subjected to stabilizing selection. These findings may provide a novel insight into the mechanisms of the diseases associated with HLA-DRB1 genes.

Heligmosomoides polygyrus infection is associated with lower MHC class II gene expression in Apodemus flavicollis: indication for immune suppression?

Due to their key role in recognizing foreign antigens and triggering the subsequent immune response the genes of the major histocompatibility complex (MHC) provide a potential target for parasites to attack in order to evade detection and expulsion from the host. A diminished MHC gene expression results in less activated T cells and might serve as a gateway for pathogens and parasites. Some parasites are suspected to be immune suppressors and promote co-infections of other parasites even in other parts of the body. In our study we found indications that the gut dwelling nematode Heligmosomoides polygyrus might exert a systemic immunosuppressive effect in yellow-necked mice (Apodemus flavicollis). The amount of hepatic MHC class II DRB gene RNA transcripts in infected mice was negatively associated with infection intensity with H. polygyrus. The hepatic expression of immunosuppressive cytokines, such as transforming growth factor β and interleukin 10 was not associated with H. polygyrus infection. We did not find direct positive associations of H. polygyrus with other helminth species. But the prevalence and infection intensity of the nematodes Syphacia stroma and Trichuris muris were higher in multiple infected individuals. Furthermore, our data indicated antagonistic effects in the helminth community of A. flavicollis as cestode infection correlated negatively with H. polygyrus and helminth species richness. Our study shows that expression analyses of immune relevant genes can also be performed in wildlife, opening new aspects and possibilities for future ecological and evolutionary research.

Evidence that purifying selection acts on promoter sequences

We tested whether functionally important sites in bacterial, yeast, and animal promoters are more conserved than their neighbors. We found that substitutions are predominantly seen in less important sites and that those that occurred tended to have less impact on gene expression than possible alternatives. These results suggest that purifying selection operates on promoter sequences.

Balancing selection is common in the extended MHC region but most alleles with opposite risk profile for autoimmune diseases are neutrally evolving

Background

Several susceptibility genetic variants for autoimmune diseases have been identified. A subset of these polymorphisms displays an opposite risk profile in different autoimmune conditions. This observation open interesting questions on the evolutionary forces shaping the frequency of these alleles in human populations.

We aimed at testing the hypothesis whereby balancing selection has shaped the frequency of opposite risk alleles.

Results

Since balancing selection signatures are expected to extend over short genomic portions, we focused our analyses on 11 regions carrying putative functional polymorphisms that may represent the disease variants (and the selection targets). No exceptional nucleotide diversity was observed for ZSCAN23, HLA-DMB, VARS2, PTPN22, BAT3, C6orf47, and IL10; summary statistics were consistent with evolutionary neutrality for these gene regions. Conversely, CDSN/PSORS1C1, TRIM10/TRIM40, BTNL2, and TAP2 showed extremely high nucleotide diversity and most tests rejected neutrality, suggesting the action of balancing selection. For TAP2 and BTNL2 these signatures are not secondary to linkage disequilibrium with HLA class II genes. Nonetheless, with the exception of variants in TRIM40 and CDSN, our data suggest that opposite risk SNPs are not selection targets but rather have accumulated as neutral variants.

Conclusion

Data herein indicate that balancing selection is common within the extended MHC region and involves several non-HLA loci. Yet, the evolutionary history of most SNPs with an opposite effect for autoimmune diseases is consistent with evolutionary neutrality. We suggest that variants with an opposite effect on autoimmune diseases should not be considered a distinct class of disease alleles from the evolutionary perspective and, in a few cases, the opposite effect on distinct diseases may derive from complex haplotype structures in regions with high genetic diversity.

Gene expression differences among primates are associated with changes in a histone epigenetic modification

Changes in gene regulation are thought to play an important role in speciation and adaptation, especially in primates. However, we still know relatively little about the mechanisms underlying regulatory evolution. In particular, the extent to which epigenetic modifications underlie gene expression differences between primates is not yet known. Our study focuses on an epigenetic histone modification, H3K4me3, which is thought to promote transcription. To investigate the contribution of H3K4me3 to regulatory differences between species, we collected gene expression data and identified H3K4me3-associated genomic regions in lymphoblastoid cell lines (LCLs) from humans, chimpanzees, and rhesus macaques, using three cell lines from each species. We found strong evidence for conservation of H3K4me3 localization in primates. Moreover, regardless of species, H3K4me3 is consistently enriched near annotated transcription start sites (TSS), and highly expressed genes are more likely than lowly expressed genes to have the histone modification near their TSS. Interestingly, we observed an enrichment of interspecies differences in H3K4me3 at the TSS of genes that are differentially expressed between species. We estimate that as much as 7% of gene expression differences between the LCLs of humans, chimpanzees, and rhesus macaques may be explained, at least in part, by changes in the status of H3K4me3 histone modifications. Our results suggest a modest, yet important role for epigenetic changes in gene expression differences between primates.

A genome-wide study of DNA methylation patterns and gene expression levels in multiple human and chimpanzee tissues

The modification of DNA by methylation is an important epigenetic mechanism that affects the spatial and temporal regulation of gene expression. Methylation patterns have been described in many contexts within and across a range of species. However, the extent to which changes in methylation might underlie inter-species differences in gene regulation, in particular between humans and other primates, has not yet been studied. To this end, we studied DNA methylation patterns in livers, hearts, and kidneys from multiple humans and chimpanzees, using tissue samples for which genome-wide gene expression data were also available. Using the multi-species gene expression and methylation data for 7,723 genes, we were able to study the role of promoter DNA methylation in the evolution of gene regulation across tissues and species. We found that inter-tissue methylation patterns are often conserved between humans and chimpanzees. However, we also found a large number of gene expression differences between species that might be explained, at least in part, by corresponding differences in methylation levels. In particular, we estimate that, in the tissues we studied, inter-species differences in promoter methylation might underlie as much as 12%-18% of differences in gene expression levels between humans and chimpanzees.

Long-term balancing selection maintains trans-specific polymorphisms in the human TRIM5 gene

The human TRIM5 genes encodes a retroviral restriction factor (TRIM5α). Evolutionary analyses of this gene in mammals have revealed a complex and multifaceted scenario, suggesting that TRIM5 has been the target of exceptionally strong selective pressures, possibly exerted by recurrent waves of retroviral infections. TRIM5 displays inter-individual expression variability in humans and high levels of TRIM5 mRNA have been associated with a reduced risk of HIV-1 infection. We resequenced TRIM5 in chimpanzees and identified two polymorphisms in intron 1 that are shared with humans. Analysis of the gene region encompassing the two trans-specific variants in human populations identified exceptional nucleotide diversity levels and an excess of polymorphism compared to fixed divergence. Most tests rejected the null hypothesis of neutral evolution for this region and haplotype analysis revealed the presence of two deeply separated clades. Calculation of the time to the most recent common ancestor (TMRCA) for TRIM5 haplotypes yielded estimates ranging between 4 and 7 million years. Overall, these data indicate that long-term balancing selection, an extremely rare process outside MHC genes, has maintained trans-specific polymorphisms in the first intron of TRIM5. Bioinformatic analyses indicated that variants in intron 1 may affect transcription factor-binding sites and, therefore, TRIM5 transcriptional activity. Data herein confirm an extremely complex evolutionary history of TRIM5 genes in primates and open the possibility that regulatory variants in the gene modulate the susceptibility to HIV-1.

Evolutionary genetics in wild primates: combining genetic approaches with field studies of natural populations

Ecological and evolutionary studies of wild primates hold important keys to understanding both the shared characteristics of primate biology and the genetic and phenotypic differences that make specific lineages, including our own, unique. Although complementary genetic research on nonhuman primates has long been of interest, recent technological and methodological advances now enable functional and population genetic studies in an unprecedented manner. In the past several years, novel genetic data sets have revealed new information about the demographic history of primate populations and the genetics of adaptively important traits. In combination with the rich history of behavioral, ecological, and physiological work on natural primate populations, genetic approaches promise to provide a compelling picture of primate evolution in the past and in the present day.

Regulation of major histocompatibility complex class II gene expression, genetic variation and disease

Major histocompatibility complex (MHC) class II molecules are central to adaptive immune responses and maintenance of self-tolerance. Since the early 1970s, the MHC class II region at chromosome 6p21 has been shown to be associated with a remarkable number of autoimmune, inflammatory and infectious diseases. Given that a full explanation for most MHC class II disease associations has not been reached through analysis of structural variation alone, in this review we examine the role of genetic variation in modulating gene expression. We describe the intricate architecture of the MHC class II regulatory system, indicating how its unique characteristics may relate to observed associations with disease. There is evidence that haplotype-specific variation involving proximal promoter sequences can alter the level of gene expression, potentially modifying the emergence and expression of key phenotypic traits. Although much emphasis has been placed on cis-regulatory elements, we also examine the role of more distant enhancer elements together with the evidence of dynamic inter- and intra-chromosomal interactions and epigenetic processes. The role of genetic variation in such mechanisms may hold profound implications for susceptibility to common disease.

Genomic features that predict allelic imbalance in humans suggest patterns of constraint on gene expression variation

Variation in gene expression is an important contributor to phenotypic diversity within and between species. Although this variation often has a genetic component, identification of the genetic variants driving this relationship remains challenging. In particular, measurements of gene expression usually do not reveal whether the genetic basis for any observed variation lies in cis or in trans to the gene, a distinction that has direct relevance to the physical location of the underlying genetic variant, and which may also impact its evolutionary trajectory. Allelic imbalance measurements identify cis-acting genetic effects by assaying the relative contribution of the two alleles of a cis-regulatory region to gene expression within individuals. Identification of patterns that predict commonly imbalanced genes could therefore serve as a useful tool and also shed light on the evolution of cis-regulatory variation itself. Here, we show that sequence motifs, polymorphism levels, and divergence levels around a gene can be used to predict commonly imbalanced genes in a human data set. Reduction of this feature set to four factors revealed that only one factor significantly differentiated between commonly imbalanced and nonimbalanced genes. We demonstrate that these results are consistent between the original data set and a second published data set in humans obtained using different technical and statistical methods. Finally, we show that variation in the single allelic imbalance-associated factor is partially explained by the density of genes in the region of a target gene (allelic imbalance is less probable for genes in gene-dense regions), and, to a lesser extent, the evenness of expression of the gene across tissues and the magnitude of negative selection on putative regulatory regions of the gene. These results suggest that the genomic distribution of functional cis-regulatory variants in the human genome is nonrandom, perhaps due to local differences in evolutionary constraint.

Genomic analyses reveal a conserved glutathione homeostasis pathway in the invertebrate chordate Ciona intestinalis

The major thiol redox buffer glutathione (l-gamma-glutamyl-l-cysteinylglycine, GSH) is central to cell fate determination, and thus, associated metabolic and regulatory pathways are exquisitely sensitive to a wide range of environmental cues. An imbalance of cellular redox homeostasis has emerged as a pathologic hallmark of a diverse range of human gene-environment disorders. Despite the central importance of GSH in cellular homeostasis, underlying genetic regulatory pathways remain poorly defined. This report describes the annotation and expression analysis of genes contributing to GSH homeostasis in the invertebrate chordate Ciona intestinalis. A core pathway comprising 19 genes contributing to the biosynthesis of GSH and its use as both a redox buffer and a conjugate in phase II detoxification as well as known transcriptional regulators were analyzed. These genes exhibit a high level of sequence conservation with corresponding human, rat, and mouse homologs and were expressed constitutively in tissues of adult animals. The GSH biosynthetic genes Gclc and Gclm were also responsive to the prototypical antioxidant tert-butylhydroquinone. The present evidence of a conserved GSH homeostasis pathway in C. intestinalis together with its phylogenetic position as a basal chordate and lifestyle as a filter feeder constantly exposed to natural marine toxins introduces this species as an important animal model for defining molecular mechanisms that potentially underlie genetic susceptibility to environmentally associated stress.

Functional and geographical differentiation of candidate balanced polymorphisms in Arabidopsis thaliana

Molecular population genetic analysis of three chromosomal regions in Arabidopsis thaliana suggested that balancing selection might operate to maintain variation at three novel candidate adaptive trait genes, including SOLUBLE STARCH SYNTHASE I (SSI), PLASTID TRANSCRIPTIONALLY ACTIVE 7(PTAC7), and BELL-LIKE HOMEODOMAIN 10 (BLH10). If balanced polymorphisms are indeed maintained at these loci, then we would expect to observe functional variation underlying the previously detected signatures of selection. We observe multiple replacement polymorphisms within and in the 32 amino acids just upstream of the protein-protein interacting BELL domain at the BLH10 locus. While no clear protein sequence differences are found between allele types in SSI and PTAC7, these two genes show evidence for allele-specific variation in expression levels. Geographical patterns of allelic differentiation seem consistent with population stratification in this species and a significant longitudinal cline was observed at all three candidate loci. These data support a hypothesis of balancing selection at all three candidate loci and provide a basis for more detailed functional work by identifying possible functional differences that might be selectively maintained.

The signature of long-standing balancing selection at the human defensin beta-1 promoter

BACKGROUND

Defensins, small endogenous peptides with antimicrobial activity, are pivotal components of the innate immune response. A large cluster of defensin genes is located on human chromosome 8p; among them the beta defensin 1 (DEFB1) promoterhas been extensively studied since discovery that specific polymorphisms and haplotypes associate with asthma and atopy, susceptibility to severe sepsis, as well as HIV and Candida infection predisposition.

RESULTS

Here, we characterize the sequence variation and haplotype structure of the DEFB1 promoter region in six human populations. In all of them, we observed high levels of nucleotide variation, an excess of intermediate-frequency alleles, reduced population differentiation and a genealogy with common haplotypes separated by deep branches. Indeed, a significant departure from the expectation of evolutionary neutrality was observed in all populations and the possibility that this is due to demographic history alone was ruled out. Also, we verified that the selection signature is restricted to the promoter region and not due to a linked balanced polymorphism. A phylogeny-based estimation indicated that the two major haplotype clades separated around 4.5 million years ago, approximately the time when the human and chimpanzee lineages split.

CONCLUSION

Altogether, these features represent strong molecular signatures of long-term balancing selection, a process that is thought to be extremely rare outside major histocompatibility complex genes. Our data indicate that the DEFB1 promoter region carries functional variants and support previous hypotheses whereby alleles predisposing to atopic disorders are widespread in modern societies because they conferred resistance to pathogens in ancient settings.

Human MHC architecture and evolution: implications for disease association studies

Major histocompatibility complex (MHC) variation is a key determinant of susceptibility and resistance to a large number of infectious, autoimmune and other diseases. Identification of the MHC variants conferring susceptibility to disease is problematic, due to high levels of variation and linkage disequilibrium. Recent cataloguing and analysis of variation over the complete MHC has facilitated localization of susceptibility loci for autoimmune diseases, and provided insight into the MHC's evolution. This review considers how the unusual genetic characteristics of the MHC impact on strategies to identify variants causing, or contributing to, disease phenotypes. It also considers the MHC in relation to novel mechanisms influencing gene function and regulation, such as epistasis, epigenetics and microRNAs. These developments, along with recent technological advances, shed light on genetic association in complex disease.

Genetic evidence reveals temporal change in hybridization patterns in a wild baboon population

The process and consequences of hybridization are of interest to evolutionary biologists because of the importance of hybridization in understanding reproductive isolation, speciation, and the influence of introgression on population genetic structure. Recent studies of hybridization have been enhanced by the advent of sensitive, genetic marker-based techniques for inferring the degree of admixture occurring within individuals. Here we present a genetic marker-based analysis of hybridization in a large-bodied, long-lived mammal over multiple generations. We analysed patterns of hybridization between yellow baboons (Papio cynocephalus) and anubis baboons (Papio anubis) in a well-studied natural population in Amboseli National Park, Kenya, using genetic samples from 450 individuals born over the last 36 years. We assigned genetic hybrid scores based on genotypes at 14 microsatellite loci using the clustering algorithm implemented in STRUCTURE 2.0, and assessed the robustness of these scores by comparison to pedigree information and through simulation. The genetic hybrid scores showed generally good agreement with previous morphological assessments of hybridity, but suggest that genetic methods may be more sensitive for identification of low levels of hybridity. The results of our analysis indicate that the proportion of hybrids in the Amboseli population has grown over time, but that the average proportion of anubis ancestry within hybrids is gradually decreasing. We argue that these patterns are probably a result of both selective and nonselective processes, including differences in the timing of life-history events for hybrid males relative to yellow baboon males, and stochasticity in long-distance dispersal from the source anubis population into Amboseli.

Age at maturity in wild baboons: genetic, environmental and demographic influences

The timing of early life-history events, such as sexual maturation and first reproduction, can greatly influence variation in individual fitness. In this study, we analysed possible sources of variation underlying different measures of age at social and physical maturation in wild baboons in the Amboseli basin, Kenya. The Amboseli baboons are a natural population primarily comprised of yellow baboons (Papio cynocephalus) that occasionally hybridize with anubis baboons (Papio anubis) from outside the basin. We found that males and females differed in the extent to which various factors influenced their maturation. Surprisingly, we found that male maturation was most strongly related to the proportion of anubis ancestry revealed by their microsatellite genotypes: hybrid males matured earlier than yellow males. In contrast, although hybrid females reached menarche slightly earlier than yellow females, maternal rank and the presence of maternal relatives had the largest effects on female maturation, followed by more modest effects of group size and rainfall. Our results indicate that a complex combination of demographic, genetic, environmental, and maternal effects contribute to variation in the timing of these life-history milestones.

Genomic analysis of Ovis aries (Ovar) MHC class IIa loci

Determining the genomic organization of the Ovis aries (Ovar) major histocompatibility complex class IIa region is essential for future functional studies related to antigen presentation. In this study, a bacterial artificial chromosome (BAC) library of genomic DNA from peripheral blood leukocytes (PBL) of a Rambouillet ram was constructed, and BAC clone consisting of the major histocompatibility complex (MHC) class II DQB2, DQA2, DQB1, DQA1, and DRB1 loci was identified and completely sequenced. The BAC clone consists of 160,889 bp of finished sequence with the loci arranged in the following order: DQB2, DQA2, DQB1, DQA1, and DRB1 with 14.3, 25, 6.6, and 40.9 Kb spanning between the loci, respectively. All five of these loci were transcribed in the animal used to generate the MHC class II BAC clone. Repeat or retrotransposable elements along with MHC class II cis promoter elements consisting of S, X, and Y boxes were identified in the sequence. In addition, 16 non-coding conserved sequences amongst primates, carnivores, and ruminants were identified (p < 0.001). These conserved sequences include binding sites for transcription factors with known roles in immune cells, and they provide a basis for further functional investigation of the genes in this region. This is the first ruminant finished sequence of the DQB2-DRB1 region, and this sequence information will aid in whole genome and transcriptome analyses of MHC class II.

Parallel effects of genetic variation in ACE activity in baboons and humans

Like humans, savannah baboons (Papio sp.) show heritable interindividual variation in complex physiological phenotypes. One prominent example of such variation involves production of the homeostatic regulator protein angiotensin converting enzyme (ACE), which shows heritable variation in both baboons and humans. In humans, this phenotypic variation is associated with an Alu insertion-deletion polymorphism in the ACE gene, which explains approximately half of the variation in serum ACE activity. We identified a similar Alu insertion-deletion polymorphism in the baboon ACE homologue and measured its frequency in a wild population and a captive population of baboons. We also analyzed the contribution of ACE genotype at this indel to variation in serum ACE activity in the captive population. When conditioned on weight, a known factor affecting ACE activity in humans, age and ACE genotype both accounted for variance in ACE activity; in particular, we identified a significant nonadditive interaction between age and genotype. A model incorporating this interaction effect explained 21.6% of the variation in residual serum ACE activity. Individuals homozygous for the deletion mutation exhibited significantly higher levels of ACE activity than insertion-deletion heterozygotes at younger ages (10-14 years), but showed a trend towards lower levels of ACE activity compared with heterozygotes at older ages (> or =15 years). These results demonstrate an interesting parallel between the genetic architecture underlying ACE variation in humans and baboons, suggesting that further attention should be paid in humans to the relationship between ACE genetic variation and aging.

Linkage disequilibrium and recombination rate estimates in the self-incompatibility region of Arabidopsis lyrata

Genetic diversity is unusually high at loci in the S-locus region of the self-incompatible species of the flowering plant, Arabidopsis lyrata, not just in the S loci themselves, but also at two nearby loci. In a previous study of a single natural population from Iceland, we attributed this elevated polymorphism to linkage disequilibrium (LD) between variants at loci close to the S locus and the S alleles, which are maintained in the population by balancing selection. With the four S-flanking loci whose diversity we previously studied, we could not determine the extent of the region linked to the S loci in which neutral sites are affected. We also could not exclude the possibility of a population bottleneck, or of admixture, as causes of the LD. We have now studied four more distant loci flanking the S-locus region, and more populations, and we analyze the results using a theoretical model of the effect of balancing selection on diversity at linked neutral sites within and between different functional S-allelic classes. In the model, diversity is a function of the number of selectively maintained alleles and the recombination distances from the selectively maintained sites. We use the model to estimate the number of different functional S alleles, their turnover rate, and recombination rates between the S-locus region and other loci. Our estimates suggest that there is a small region of very low recombination surrounding the S-locus region.