Natural selection at major histocompatibility complex loci of vertebrates

[Cloning and sequencing HLA-A and -B genomic DNA and analyzing polymorphism in regulatory regions in Chinese Han individuals]

In the present study, a high-resolute method for cloning and sequencing genomic full-length HLA-A and -B using 20 Chinese Han individuals was established. We detected 10 HLA-A allele sequences 4.2 kb in length and 6 HLA-B allele sequences 3.7 kb in length, and the sequences included all exons, all introns, 5'promoter, and 3'UTR of the two genes. All sixteen sequences have been submitted to GenBank and IMGT/HLA database. A*1153 is a novel allele, and the introns of B*151101 are firstly reported here. The 5'promoter and 3'UTR sequences of 5 HLA-A alleles and 2 HLA-B alleles are also firstly disclosed, and all other alleles have extended the genomic full length sequences released in IMGT/HLA database. The polymorphic structures of upper 5'promoter and downstream 3'UTR, which were uncovered in IMGT/HLA database, are firstly depicted in Chinese Han individuals. Twenty-six single nucleotide polymorphisms (SNPs) and one 3 bp-insertion/deletion (Indel) were located in the upper 5'promoter and 14 SNPs were located in the 3'UTR of HLA-A. In addition, five SNPs and one 1 bp-indel were located in the upper 5'promoter and 5 SNPs were located in the 3'UTR of HLA-B. Through analyzing the phylogenetic trees of 5'promoter, exons and 3'UTR of the two genes, we found that the evolution history of regulatory regions and exons is different between the two genes. The regulatory regions are tightly linked with exons in most of HLA-A alleles excluding A*24020101. On the contrary, recombinant events may occur frequently between regulatory regions and exons in most HLA-B alleles.

Characterization of major histocompatibility complex (MHC) DRB exon 2 and DRA exon 3 fragments in a primary terrestrial rabies vector (Procyon lotor)

The major histocompatibility complex (MHC) presents a unique system to explore links between genetic diversity and pathogens, as diversity within MHC is maintained in part by pathogen driven selection. While the majority of wildlife MHC studies have investigated species that are of conservation concern, here we characterize MHC variation in a common and broadly distributed species, the North American raccoon (Procyon lotor). Raccoons host an array of broadly distributed wildlife diseases (e.g., canine distemper, parvovirus and raccoon rabies virus) and present important human health risks as they persist in high densities and in close proximity to humans and livestock. To further explore how genetic variation influences the spread and maintenance of disease in raccoons we characterized a fragment of MHC class II DRA exon 3 (250bp) and DRB exon 2 (228 bp). MHC DRA was found to be functionally monomorphic in the 32 individuals screened; whereas DRB exon 2 revealed 66 unique alleles among the 246 individuals screened. Between two and four alleles were observed in each individual suggesting we were amplifying a duplicated DRB locus. Nucleotide differences between DRB alleles ranged from 1 to 36 bp (0.4–15.8% divergence) and translated into 1 to 21 (1.3–27.6% divergence) amino acid differences. We detected a significant excess of nonsynonymous substitutions at the peptide binding region (P = 0.005), indicating that DRB exon 2 in raccoons has been influenced by positive selection. These data will form the basis of continued analyses into the spatial and temporal relationship of the raccoon rabies virus and the immunogenetic response in its primary host.

Antagonistic coevolution with parasites maintains host genetic diversity: an experimental test

Genetic variation in natural populations is a prime prerequisite allowing populations to respond to selection, but is under constant threat from forces that tend to reduce it, such as genetic drift and many types of selection. Haldane emphasized the potential importance of parasites as a driving force of genetic diversity. His theory has been taken for granted ever since, but despite numerous studies showing correlations between genetic diversity and parasitism, Haldane's hypothesis has rarely been tested experimentally for unambiguous support. We experimentally staged antagonistic coevolution between the host Tribolium castaneum and its natural microsporidian parasite, Nosema whitei, to test for the relative importance of two separate evolutionary forces (drift and parasite-induced selection) on the maintenance of genetic variation. Our results demonstrate that coevolution with parasites indeed counteracts drift as coevolving populations had significantly higher levels of heterozygosity and allelic diversity. Genetic drift remained a strong force, strongly reducing genetic variation and increasing genetic differentiation in small populations. To our surprise, differentiation between the evolving populations was smaller when they coevolved with parasites, suggesting parallel balancing selection. Hence, our results experimentally vindicate Haldane's original hypothesis 60 years after its conception.

A novel single cDNA amplicon pyrosequencing method for high-throughput, cost-effective sequence-based HLA class I genotyping

Human leukocyte antigen (HLA) genotype influences the immune response to pathogens and transplanted tissues; accurate HLA genotyping is critical for clinical and research applications. Sequence-based HLA typing is limited by the cost of Sanger sequencing genomic DNA (gDNA) and resolving cis/trans ambiguities, hindering both studies correlating high-resolution genotype with clinical outcomes, and population-specific allele frequency surveys. We present an assay for sequence-based HLA genotyping by titanium read length clonal Roche/454 pyrosequencing of a single, universally diagnostic polymerase chain reaction (PCR) amplicon from HLA class I cDNA that captures most of exons 2, 3, and 4 used for traditional sequence-based typing. The amplicon is predicted to unambiguously resolve 85% of known alleles. A panel of 48 previously HLA-typed samples was assayed with this method, demonstrating 100% non-null allele typing concordance. We show that this technique can multiplex at least 768 patients per sequencing run with multiplex identifier sequence bar-coding. Unprecedented typing throughput results from a novel single cDNA-PCR amplicon strategy requiring only 1 PCR amplification per sample. This method dramatically reduces cost for genotyping of large cohorts.

Characterisation of MHC class II DRB genes in the northern tree shrew (Tupaia belangeri)

Genes of the major histocompatibility complex (MHC) mainly code for proteins of the immune system of jawed vertebrates. In particular, MHC class I and II cell surface proteins are crucial for the self/non-self discrimination of the adaptive immune system and are the most polymorphic genes in vertebrates. Positive selection, gene duplications and pseudogenes shape the face of the MHC and reflect a highly dynamic evolution. Here, we present for the first time data of the highly polymorphic MHC class II DRB exon 2 of a representative of the mammalian order scandentia, the northern tree shrew Tupaia belangeri. We found up to eight different alleles per individual and determined haplotype constitution by intensively studying their inheritance. The alleles were assigned to four putative loci, all of which were polymorphic. Only the most polymorphic locus was subject to positive selection within the antigen binding sites and only alleles of this locus were transcribed.

Natural selection coupled with intragenic recombination shapes diversity patterns in the major histocompatibility complex class II genes of the giant panda

Ample variations of the major histocompatibility complex (MHC) genes are essential for vertebrates to adapt to various environmental conditions. In this study, we investigated the genetic variations and evolutionary patterns of seven functional MHC class II genes (one DRA, two DRB, two DQA, and two DQB) of the giant panda. The results showed the presence of two monomorphic loci (DRA and DQB2) and five polymorphic loci with different numbers of alleles (seven at DRB1, six at DRB3, seven at DQA1, four at DQA2, six at DQB1). The presence of balancing selection in the giant panda was supported by the following pieces of evidence: (1) The observed heterozygosity was higher than expected. (2) Amino acid heterozygosity was significantly higher at antigen-binding sites (ABS) compared with non-ABS sequences. (3) The selection parameter omega (d(N)/d(S)) was significantly higher at ABS compared with non-ABS sequences. (4) Approximately 95.45% of the positively selected codons (P>0.95) were located at or adjacent to an ABS. Furthermore, this study showed that (1) The Qinling subspecies exhibited high omega values across each locus (all >1), supporting its extensive positive selection. (2) The Sichuan subspecies displayed small omega at DRB1 (omega<0.72) and DQA2 (omega<0.48), suggesting that these sites underwent strong purifying selection. (3) Intragenic recombination was detected in DRB1, DQA1, and DQB1. The molecular diversity in classic Aime-MHC class II genes implies that the giant panda had evolved relatively abundant variations in its adaptive immunity along the history of host-pathogen co-evolution. Collectively, these findings indicate that natural selection accompanied by recombination drives the contrasting diversity patterns of the MHC class II genes between the two studied subspecies of giant panda.

Trans-species polymorphism and selection in the MHC class II DRA genes of domestic sheep

Highly polymorphic genes with central roles in lymphocyte mediated immune surveillance are grouped together in the major histocompatibility complex (MHC) in higher vertebrates. Generally, across vertebrate species the class II MHC DRA gene is highly conserved with only limited allelic variation. Here however, we provide evidence of trans-species polymorphism at the DRA locus in domestic sheep (Ovis aries). We describe variation at the Ovar-DRA locus that is far in excess of anything described in other vertebrate species. The divergent DRA allele (Ovar-DRA*0201) differs from the sheep reference sequences by 20 nucleotides, 12 of which appear non-synonymous. Furthermore, DRA*0201 is paired with an equally divergent DRB1 allele (Ovar-DRB1*0901), which is consistent with an independent evolutionary history for the DR sub-region within this MHC haplotype. No recombination was observed between the divergent DRA and B genes in a range of breeds and typical levels of MHC class II DR protein expression were detected at the surface of leukocyte populations obtained from animals homozygous for the DRA*0201, DRB1*0901 haplotype. Bayesian phylogenetic analysis groups Ovar-DRA*0201 with DRA sequences derived from species within the Oryx and Alcelaphus genera rather than clustering with other ovine and caprine DRA alleles. Tests for Darwinian selection identified 10 positively selected sites on the branch leading to Ovar-DRA*0201, three of which are predicted to be associated with the binding of peptide antigen. As the Ovis, Oryx and Alcelaphus genera have not shared a common ancestor for over 30 million years, the DRA*0201 and DRB1*0901 allelic pair is likely to be of ancient origin and present in the founding population from which all contemporary domestic sheep breeds are derived. The conservation of the integrity of this unusual DR allelic pair suggests some selective advantage which is likely to be associated with the presentation of pathogen antigen to T-cells and the induction of protective immunity.

Class II DRB polymorphism and sequence diversity in two vesper bats in the genus Myotis

Almost no studies have been done with respect to major histocompatibility complex (MHC) polymorphism and sequence diversity in bats, although they account for one in five living mammalian species. We analysed MHC Class II DRB polymorphism and sequence diversity in two Mexican verpertilionid bat species, the widespread continental species Myotis velifer and the narrowly distributed (and endangered) island endemic Myotis vivesi. We find extensive DRB polymorphism in the widespread M. velifer, similar to that commonly reported in other mammals. The geographically restricted M. vivesi by contrast shows only very limited polymorphism. We conclude that M. vivesi has undergone a dramatic loss of MHC polymorphism. The significance of this inference in light of other information on population structure and genetic diversity in this species is discussed.

MHC evolution in three salmonid species: a comparison between class II alpha and beta genes

The genes of the major histocompatibility complex (MHC) are amongst the most variable in vertebrates and represent some of the best candidates to study processes of adaptive evolution. However, despite the number of studies available, most of the information on the structure and function of these genes come from studies in mammals and birds in which the MHC class I and II genes are tightly linked and class II alpha exhibits low variability in many cases. Teleost fishes are among the most primitive vertebrates with MHC and represent good organisms for the study of MHC evolution because their class I and class II loci are not physically linked, allowing for independent evolution of both classes of genes. We have compared the diversity and molecular mechanisms of evolution of classical MH class II alpha and class II beta loci in farm populations of three salmonid species: Oncorhynchus kisutch, Oncorhynchus mykiss and Salmo salar. We found single classical class II loci and high polymorphism at both class II alpha and beta genes in the three species. Mechanisms of evolution were common for both class II genes, with recombination and point mutation involved in generating diversity and positive selection acting on the peptide-binding residues. These results suggest that the maintenance of variability at the class IIalpha gene could be a mechanism to increase diversity in the MHC class II in salmonids in order to compensate for the expression of one single classical locus and to respond to a wider array of parasites.

High-throughput high-resolution class I HLA genotyping in East Africa

HLA, the most genetically diverse loci in the human genome, play a crucial role in host-pathogen interaction by mediating innate and adaptive cellular immune responses. A vast number of infectious diseases affect East Africa, including HIV/AIDS, malaria, and tuberculosis, but the HLA genetic diversity in this region remains incompletely described. This is a major obstacle for the design and evaluation of preventive vaccines. Available HLA typing techniques, that provide the 4-digit level resolution needed to interpret immune responses, lack sufficient throughput for large immunoepidemiological studies. Here we present a novel HLA typing assay bridging the gap between high resolution and high throughput. The assay is based on real-time PCR using sequence-specific primers (SSP) and can genotype carriers of the 49 most common East African class I HLA-A, -B, and -C alleles, at the 4-digit level. Using a validation panel of 175 samples from Kampala, Uganda, previously defined by sequence-based typing, the new assay performed with 100% sensitivity and specificity. The assay was also implemented to define the HLA genetic complexity of a previously uncharacterized Tanzanian population, demonstrating its inclusion in the major East African genetic cluster. The availability of genotyping tools with this capacity will be extremely useful in the identification of correlates of immune protection and the evaluation of candidate vaccine efficacy.

[Single nucleotide polymorphisms in 22 HLA-Cw alleles in Chinese Han population]

To analyze the molecular genetic polymorphism of full-length HLA-Cw gene, a total of 28 samples with known genotypes from Chinese Han population were amplified by long-range PCR using high-fidelity Pfu polymerase. A fragment 4.5 kb in length of HLA-Cw gene was subjected to cloning and haplotype sequencing. The single nucleotide polymorphisms (SNPs) in all segments of the whole region of HLA-Cw gene were analyzed. As a result, we detected 22 different HLA-Cw alleles in 28 samples, all of which were submitted to GenBank and the IMGT/HLA Database. Among the 22 HLA-Cw alleles, the intronic sequences of Cw*030301, Cw*0706 and Cw*140201 were firstly elucidated. The novel intronic sequence and the SNPs information may help to design allele-specific primers for accurate sequence-based typing (SBT) and to avoid allele dropout events in SBT test. We aligned all the diploid sequences using ClustalX program and imported them into Dnasp4.0 to calculate polymorphism in all coding- and non-coding regions. We found 244 SNPs and 10 insertion/deletions (Indels). According to the analysis of polymorphism level, phylogenetic trees and frequency spectrum, we proposed that the evolution of intron 4 and exon 5 was under balancing selection. Selection on these segments indicated that they may be functionally important in evolution of HLA-Cw gene. The full-length sequences obtained and related SNPs information can be used as resources of markers for high-resolution typing, complex diseases association studies and human evolution.

A 454 multiplex sequencing method for rapid and reliable genotyping of highly polymorphic genes in large-scale studies

Background

High-throughput sequencing technologies offer new perspectives for biomedical, agronomical and evolutionary research. Promising progresses now concern the application of these technologies to large-scale studies of genetic variation. Such studies require the genotyping of high numbers of samples. This is theoretically possible using 454 pyrosequencing, which generates billions of base pairs of sequence data. However several challenges arise: first in the attribution of each read produced to its original sample, and second, in bioinformatic analyses to distinguish true from artifactual sequence variation. This pilot study proposes a new application for the 454 GS FLX platform, allowing the individual genotyping of thousands of samples in one run. A probabilistic model has been developed to demonstrate the reliability of this method.

Results

DNA amplicons from 1,710 rodent samples were individually barcoded using a combination of tags located in forward and reverse primers. Amplicons consisted in 222 bp fragments corresponding to DRB exon 2, a highly polymorphic gene in mammals. A total of 221,789 reads were obtained, of which 153,349 were finally assigned to original samples. Rules based on a probabilistic model and a four-step procedure, were developed to validate sequences and provide a confidence level for each genotype. The method gave promising results, with the genotyping of DRB exon 2 sequences for 1,407 samples from 24 different rodent species and the sequencing of 392 variants in one half of a 454 run. Using replicates, we estimated that the reproducibility of genotyping reached 95%.

Conclusions

This new approach is a promising alternative to classical methods involving electrophoresis-based techniques for variant separation and cloning-sequencing for sequence determination. The 454 system is less costly and time consuming and may enhance the reliability of genotypes obtained when high numbers of samples are studied. It opens up new perspectives for the study of evolutionary and functional genetics of highly polymorphic genes like major histocompatibility complex genes in vertebrates or loci regulating self-compatibility in plants. Important applications in biomedical research will include the detection of individual variation in disease susceptibility. Similarly, agronomy will benefit from this approach, through the study of genes implicated in productivity or disease susceptibility traits.

Evolution of an MHC class Ia gene fragment in four North American Morone species

A nucleotide sequence analysis of a fragment of a Morone MHC class Ia gene detected high levels of polymorphism in striped bass Morone saxatilis, white perch Morone americana and yellow bass Morone mississippiensis. Extremely low levels of MHC diversity, however, were detected in white bass Morone chrysops, suggesting the possibility of a severe population bottleneck for this species.

Ancestral polymorphism in exon 2 of bluethroat (Luscinia svecica) MHC class II B genes

The genes of the major histocompatibility complex (MHC) are important model genes for understanding selective forces in evolution. Here, we document, using a cloning and sequencing approach, high polymorphism at the exon 2 of the MHC class II B (MHCIIB) genes in the bluethroat (Luscinia svecica); a minimum of 61 unique alleles were detected in 20 individuals, and at least 11 functional loci. In addition, several pseudogenes were revealed. The specimens originated from three different bluethroat subspecies (azuricollis, cyanecula and svecica), and we also analysed four specimens of the closely related thrush nightingale (L. luscinia) for comparison. Phylogenetic analyses of the functional alleles revealed 258 equally parsimonious trees with poor statistical support for the majority of nodes. The distribution of the sequences in the trees point to an ancestral origin of the polymorphism in MHC class II B genes, a portion of which predated the phylogenetic split between the bluethroat and the thrush nightingale. Strong signatures of balancing selection were uncovered for the codons coding for the peptide-binding residues of the functional MHCIIB exon 2 alleles. Our results highlight the importance of duplication and recombination events for shaping passerine MHC and give insights in the evolutionary dynamics of MHC variation among closely related taxa.

Lineage pattern, trans-species polymorphism, and selection pressure among the major lineages of feline MHC-DRB peptide-binding region

The long-term evolution of major histocompatibility complex (MHC) involves the birth-and-death process and independent divergence of loci during episodes punctuated by natural selection. Here, we investigated the molecular signatures of natural selection at exon-2 of MHC class II DRB gene which includes a part of the peptide-binding region (PBR) in seven of eight putative extant Felidae lineages. The DRB alleles in felids can be mainly divided into five lineages. Signatures of trans-species polymorphism among major allelic lineages indicate that balancing selection has maintained the MHC polymorphism for a long evolutionary time. Analysis based on maximum likelihood models of codon substitution revealed overall purifying selection acting on the feline DRB. Sites that have undergone positive selection and those that are under divergent selective pressure among lineages were detected and found to fall within the putative PBR. This study increased our understanding of the nature of selective forces acting on DRB during feline radiation.

High MHC DQB variation and asymmetric allelic distribution in the endangered Yangtze finless porpoise, Neophocaena phocaenoides asiaeorientalis

The endangered Yangtze finless porpoise is found in the middle and lower reaches of the Yangtze River and its adjoining big lakes. To explore the major histocompatibility complex (MHC) genetic diversity and allelic distribution patterns across its range, we investigated variation at DQB exon 2. From 76 porpoises, we identified 18 DQB sequences. The freshwater Yangtze populations had much higher allelic diversity than marine populations. Among these freshwater populations, the middle-reach population had higher allelic diversity than the lower-reach population. The high DQB diversity level, relative to that of a neutral mtDNA locus, suggests that balancing selection is acting at the DQB gene and that rapid evolution and local positive selection play critical roles in generating and retaining high MHC diversity in the freshwater population. As the balancing selection might be driven by environmental pathogens, we suggest that maintaining MHC variation should be a high priority in the conservation and management of this endangered population, especially as an ex situ conservation strategy.

Targets of balancing selection in the human genome

Balancing selection is potentially an important biological force for maintaining advantageous genetic diversity in populations, including variation that is responsible for long-term adaptation to the environment. By serving as a means to maintain genetic variation, it may be particularly relevant to maintaining phenotypic variation in natural populations. Nevertheless, its prevalence and specific targets in the human genome remain largely unknown. We have analyzed the patterns of diversity and divergence of 13,400 genes in two human populations using an unbiased single-nucleotide polymorphism data set, a genome-wide approach, and a method that incorporates demography in neutrality tests. We identified an unbiased catalog of genes with signatures of long-term balancing selection, which includes immunity genes as well as genes encoding keratins and membrane channels; the catalog also shows enrichment in functional categories involved in cellular structure. Patterns are mostly concordant in the two populations, with a small fraction of genes showing population-specific signatures of selection. Power considerations indicate that our findings represent a subset of all targets in the genome, suggesting that although balancing selection may not have an obvious impact on a large proportion of human genes, it is a key force affecting the evolution of a number of genes in humans.

Assessing human variation data for signatures of natural selection

In this article, we highlight some of the different types of natural selection, their effects on patterns of DNA variation, and some of the statistical tests that are commonly used to detect such effects. We also explain some of the relative strengths and weaknesses of different strategies that can be used to detect signatures of natural selection at individual loci. These strategies are illustrated by their application to empirical data from gene variants that are often associated with differences in disease susceptibility. We briefly outline some of the methods proposed to scan the genome for evidence of selection. Finally, we discuss some of the problems associated with identifying signatures of selection and with making inferences about the nature of the selective process.

Mhc class II diversity and balancing selection in greater prairie-chickens

The major histocompatibility complex (Mhc) of domestic chickens has been characterized as small and relatively simple compared with that of mammals. However, there is growing evidence that the Mhc of many bird lineages may be more complex, even within the Order Galliformes. In this study, we measured genetic variation and balancing selection at Mhc loci in another galliform, the greater prairie-chicken. We cloned and sequenced a 239 bp fragment of Mhc Class II beta-chain (BLB) exon 2 in 14 individuals. There was a total of 10 unique sequences and a minimum of four BLB loci. The d(N)/d(S) ratio at peptide-binding codons was significantly greater than one, suggesting balancing selection is acting on the BLB. We also recovered two YLB sequences, which clustered tightly with YLB sequences from three other species: domestic chicken, black grouse and common quail. The relatively large number of loci revealed in our study suggests that even closely related galliforms differ in the level of Mhc variation and structure.

MHC genetic structure and divergence across populations of Chinook salmon (Oncorhynchus tshawytscha)

The major histocompatibility complex (MHC) is thought to be under strong selection pressure because of its integral role in pathogen recognition. Consequently, patterns of MHC genetic variation should reflect selection pressures across the landscape. We examined genetic variation and population genetic structure at the MHC class I-A1 and class II-B1 exons in five Chinook salmon (Oncorhynchus tshawytscha) populations from two geographic regions in British Columbia, Canada. We then compared estimates of population structure at the MHC genes with neutral estimates based on microsatellites to examine the potential for local adaptation at the MHC. Chinook salmon are in decline throughout much of their native range and understanding the degree of local adaptation exhibited by the MHC may be important in conservation planning. Comparisons among populations yielded higher G'(ST) estimates for the MHC class I than expected under neutrality based on the microsatellites. In contrast, the MHC class II tended to exhibit lower G'(ST) values than did the microsatellites. These results suggest that across populations unique selection pressures are driving allele frequency differences at the MHC class I but that the MHC class II may be the subject of homogenizing selection. Rates of nonsynonymous versus synonymous substitutions found in codons associated within the MHC class I and II peptide-binding regions provided strong evidence of positive selection. Together, these results support the hypothesis that selection is influencing genetic variation at the MHC, but suggest that selection pressures may vary at the two classes of loci both at the sequence and population levels.

MHC standing genetic variation and pathogen resistance in wild Atlantic salmon

Pathogens are increasingly emerging in human-altered environments as a serious threat to biodiversity. In this context of rapid environmental changes, improving our knowledge on the interaction between ecology and evolution is critical. The objective of this study was to evaluate the influence of an immunocompetence gene, the major histocompatibility complex (MHC) class IIbeta, on the pathogen infection levels in wild Atlantic salmon populations, Salmo salar, and identify selective agents involved in contemporary coevolution. MHC variability and bacterial infection rate were determined throughout the summer in juvenile salmon from six rivers belonging to different genetic and ecological regions in Québec, Canada. A total of 13 different pathogens were identified in kidney by DNA sequence analysis, including a predominant myxozoa, most probably recently introduced in North America. Infection rates were the highest in southern rivers at the beginning of the summer (average 47.6+/-6.3% infected fish). One MHC allele conferred a 2.9 times greater chance of being resistant to myxozoa, while another allele increased susceptibility by 3.4 times. The decrease in frequency of the susceptibility allele but not other MHC or microsatellite alleles during summer was suggestive of a mortality event from myxozoa infection. These results supported the hypothesis of pathogen-driven selection in the wild by means of frequency-dependent selection or change in selection through time and space rather than heterozygous advantage, and underline the importance of MHC standing genetic variation for facing pathogens in a changing environment.

Characterization and evolution of MHC class II B genes in Galápagos marine iguanas (Amblyrhynchus cristatus)

Major histocompatibility complex (MHC) class II molecules play a key role in the adaptive immune system of vertebrates. Class II B genes appear to evolve in a very different manner in mammals and birds. Orthology is commonly observed among mammal loci, while genes tend to cluster phylogenetically within bird species. Here we present class II B data from a representative of another major group of amniotes, the squamates (i.e. lizards, snakes, amphisbaenians), with the ultimate goal of placing mammalian and avian MHC evolution into a broader context. In this study, eight class II B cDNA sequences were obtained from the Galápagos marine iguana (Amblyrhynchus cristatus) which were divided into five locus groups, Amcr-DAB1 through -DAB5, based on similarities along most of the coding and noncoding portions of the transcribed gene. All marine iguana sequences were monophyletic with respect to class II genes from other vertebrates indicating that they originated from a common ancestral locus after squamates split from other reptiles. The beta-1 domain, which is involved in antigen binding, exhibited signatures of positive selection as well as interlocus gene conversion in both long and short tracts-a pattern also observed in birds and fish, but not in mammals. On the other hand, the beta-2 domain was divergent between gene groups, which is characteristic of mammals. Based on these results, we preliminarily show that squamate class II B genes have been shaped by a unique blend of evolutionary forces that have been observed in differing degrees in other vertebrates.

Correlated evolution of nucleotide substitution rates and allelic variation in Mhc-DRB lineages of primates

BACKGROUND

The major histocompatibility complex (MHC) is a key model of genetic polymorphism. Selection pressure by pathogens or other microevolutionary forces may result in a high rate of non-synonymous substitutions at the codons specifying the contact residues of the antigen binding sites (ABS), and the maintenance of extreme MHC allelic variation at the population/species level. Therefore, selection forces favouring MHC variability for any reason should cause a correlated evolution between substitution rates and allelic polymorphism. To investigate this prediction, we characterised nucleotide substitution rates and allelic polymorphism (i.e. the number of alleles detected in relation to the number of animals screened) of several Mhc class II DRB lineages in 46 primate species, and tested for a correlation between them.

RESULTS

First, we demonstrate that species-specific and lineage-specific evolutionary constraints favour species- and lineage-dependent substitution rate at the codons specifying the ABS contact residues (i.e. certain species and lineages can be characterised by high substitution rate, while others have low rate). Second, we show that although the degree of the non-synonymous substitution rate at the ABS contact residues was systematically higher than the degree of the synonymous substitution rate, these estimates were strongly correlated when we controlled for species-specific and lineage-specific effects, and also for the fact that different studies relied on different sample size. Such relationships between substitution rates of different types could even be extended to the non-contact residues of the molecule. Third, we provide statistical evidence that increased substitution rate along a MHC gene may lead to allelic variation, as a high substitution rate can be observed in those lineages in which many alleles are maintained. Fourth, we show that the detected patterns were independent of phylogenetic constraints. When we used phylogenetic models that control for similarity between species, due to common descent, and focused on variations within a single lineage (DRB1*03), the positive relationship between different substitution rates and allelic polymorphisms was still robust. Finally, we found the same effects to emerge in the analyses that eliminated within-species variation in MHC traits by using strictly single population-level studies. However, in a set of contrasting analyses, in which we focused on the non-functional DRB6 locus, the correlation between substitution rates and allelic variation was not prevalent.

CONCLUSION

Our results indicate that positive selection for the generation of allelic polymorphism acting on the functional part of the protein has consequences for the nucleotide substitution rate along the whole exon 2 sequence of the Mhc-DRB gene. Additionally, we proved that an increased substitution rate can promote allelic variation within lineages. Consequently, the evolution of different characteristics of genetic polymorphism is not independent.

Recent advances in understanding the biomolecular basis of chronic beryllium disease: a review

In this review we summarize the work conducted over the past decade that has advanced our knowledge of pulmonary diseases associated with exposure to beryllium that has provided a molecular-based understanding of the chemistry, immunopathology, and immunogenetics of beryllium toxicity. Beryllium is a strong and lightweight metal that generates and reflects neutrons, resists corrosion, is transparent to X-rays, and conducts electricity. Beryllium is one of the most toxic elements on the periodic table, eliciting in susceptible humans (a) an allergic immune response known as beryllium sensitization (BeS); (b) acute beryllium disease, an acutely toxic, pneumonitis-like lung condition resulting from exposure to high beryllium concentrations that are rarely seen in modern industry; and (c) chronic beryllium disease (CBD) following either high or very low levels of exposure. Because of its exceptional strength, stability, and heat-absorbing capability, beryllium is used in many important technologies in the modern world. In the early 1940s, beryllium was recognized as posing an occupational hazard in manufacturing and production settings. Although acute beryllium disease is now rare, beryllium is an insidious poison with a latent toxicity and the risk of developing CBD persists. Chronic beryllium disease-a systemic granulomatous lung disorder caused by a specific delayed immune response to beryllium within a few months to several decades after exposure-has been called the "unrecognized epidemic". Although not a disease in itself, BeS, the innate immune response to beryllium identified by an abnormal beryllium lymphocyte proliferation test result, is a population-based predictor of CBD. Genetic susceptibility to CBD is associated with alleles of the major histocompatibility gene, human leukocyte antigen DP (HLA-DP) containing glutamic acid at the 69th position of the beta chain (HLA-DPbeta-E69). Other genes are likely to be involved in the disease process, and research on this issue is in progress. The current Occupational Safety & Health Administration permissible exposure limit of 2 microg/m3 has failed to protect workers from BeS/CBD. As a safe exposure limit that will not lead to BeS or CBD has not yet been determined, the realization that the risk of CBD persists has led to a renaissance in research on the effects of the metal on human health. Current data support further reductions in exposure levels to help minimize the incidence of CBD. Steps that would directly impact both the power of epidemiologic studies and the cost of surveillance would be to develop and validate improved screening and diagnostic tests, and to identify more genetic factors that affect either sensitization or disease process. The major focus of this review is the recent research on the cellular and molecular basis of beryllium sensitization and disease, using a multidisciplinary approach of bioinorganic chemistry and immunology. First we present a historical background of beryllium exposure and disease, followed by occurrence of beryllium in the environment, toxicokinetics, biological effects, beryllium lung disease, and other human health effects.

Disentangling the effects of recombination, selection, and demography on the genetic variation at a major histocompatibility complex class II gene in the alpine chamois

The major histocompatibility complex (MHC) harbours some of the most polymorphic loci in vertebrate genomes. MHC genes are thought to be subject to some form of balancing selection, most likely pathogen-mediated selection. Hence, MHC genes are excellent candidates for exploring adaptive processes. In this study, we investigated the genetic variation at exon 2 of the DRB class II MHC locus in 191 alpine chamois (Rupicapra rupicapra) from 10 populations in the eastern Alps of Italy. In particular, we were interested in distinguishing and estimating the relative impact of selective and demographic factors, while taking into account the confounding effect of recombination. The extremely high d(n)/d(s) ratio and the presence of trans-species polymorphisms suggest that a strong long-term balancing selection effect has been operating at this locus throughout the evolutionary history of this species. We analysed patterns of genetic variation within and between populations, and the mitochondrial D-loop polymorphism patterns were analysed to provide a baseline indicator of the effects of demographic processes. These analyses showed that (i) the chamois experienced a demographic decline in the last 5000-30 000 years, most likely related to the postglacial elevation in temperature; (ii) this demographic process can explain the results of neutrality tests applied to MHC variation within populations, but cannot justify the much weaker divergence between populations implied by MHC as opposed to mitochondrial DNA; (iii) similar sets of divergent alleles are probably maintained with similar frequencies by balancing selection in different populations, and this mechanism is also operating in small isolated populations, which are strongly affected by drift.

Natural selection maintains a stable polymorphism at the circumsporozoite protein locus of Plasmodium falciparum in a low endemic area

Circumsporozoite protein gene sequences of Plasmodium falciparum were collected in 1996-1997 and in 2006-2007 from a single endemic area in Thailand. Repeat units were more similar within the same haplotype than between haplotypes, supporting the hypothesis that repeat arrays evolve by a process of concerted evolution. There was evidence that natural selection has favored amino acid changes in the Th2R and Th3R T-cell epitope regions. One haplotype in these epitopes, designated *5/*1, occurred in approximately 70% of sequences in both collection periods. The most common other haplotypes differed from *5/*1 by at least two amino acid replacements; and divergence in the epitopes was correlated with divergence in the repeats. These patterns are most consistent with balancing selection driven by interactions with the immune system of the vertebrate host, probably involving both T-cell recognition of the Th2R and Th3R epitopes and antibody responses to the repeats.

Evidence for directional selection at a novel major histocompatibility class I marker in wild common frogs (Rana temporaria) exposed to a viral pathogen (Ranavirus)

Whilst the Major Histocompatibility Complex (MHC) is well characterized in the anuran Xenopus, this region has not previously been studied in another popular model species, the common frog (Rana temporaria). Nor, to date, have there been any studies of MHC in wild amphibian host-pathogen systems. We characterise an MHC class I locus in the common frog, and present primers to amplify both the whole region, and specifically the antigen binding region. As no more than two expressed haplotypes were found in over 400 clones from 66 individuals, it is likely that there is a single class I locus in this species. This finding is consistent with the single class I locus in Xenopus, but contrasts with the multiple loci identified in axolotls, providing evidence that the diversification of MHC class I into multiple loci likely occurred after the Caudata/Anura divergence (approximately 350 million years ago) but before the Ranidae/Pipidae divergence (approximately 230 mya). We use this locus to compare wild populations of common frogs that have been infected with a viral pathogen (Ranavirus) with those that have no history of infection. We demonstrate that certain MHC supertypes are associated with infection status (even after accounting for shared ancestry), and that the diseased populations have more similar supertype frequencies (lower F_ST) than the uninfected. These patterns were not seen in a suite of putatively neutral microsatellite loci. We interpret this pattern at the MHC locus to indicate that the disease has imposed selection for particular haplotypes, and hence that common frogs may be adapting to the presence of Ranavirus, which currently kills tens of thousands of amphibians in the UK each year.

Molecular evolution of the antiretroviral TRIM5 gene

In 2004, the first report of TRIM5alpha as a cellular antiretroviral factor triggered intense interest among virologists, particularly because some primate orthologs of TRIM5alpha have activity against HIV. Since that time, a complex and eventful evolutionary history of the TRIM5 locus has emerged. A review of the TRIM5 literature constitutes a veritable compendium of evolutionary phenomena, including elevated rates of nonsynonymous substitution, divergence in subdomains due to short insertions and deletions, expansions and contractions in gene copy number, pseudogenization, balanced polymorphism, trans-species polymorphism, convergent evolution, and the acquisition of new domains by exon capture. Unlike most genes, whose history is dominated by long periods of purifying selection interspersed with rare instances of genetic innovation, analysis of restriction factor loci is likely to be complicated by the unpredictable and more-or-less constant influence of positive selection. In this regard, the molecular evolution and population genetics of restriction factor loci most closely resemble patterns that have been documented for immunity genes, such as class I and II MHC genes, whose products interact directly with microbial targets. While the antiretroviral activity encoded by TRIM5 provides plausible mechanistic hypotheses for these unusual evolutionary observations, evolutionary analyses have reciprocated by providing significant insights into the structure and function of the TRIM5alpha protein. Many of the lessons learned from TRIM5 should be applicable to the study of other restriction factor loci, and molecular evolutionary analysis could facilitate the discovery of new antiviral factors, particularly among the many TRIM genes whose functions remain as yet unidentified.

Structure and evolution of a new avian MHC class II B gene in a sub-Antarctic seabird, the thin-billed prion (Procellariiformes: Pachyptila belcheri)

The major histocompatibility complex encodes molecules that present foreign peptides to T cells of the immune system. The peptide binding region (PBR) of these molecules is among the most polymorphic regions found in vertebrate taxa. Genomic cloning approaches are improving our understanding of the evolution of this multigene family in nonmodel avian groups. By building a cosmid library, a new MHC class II B gene, Pabe-DAB1, was isolated and characterized at the genomic level in a sub-Antarctic seabird, the thin-billed prion (Pachyptila belcheri). Pabe-DAB1 exhibits the hallmark structural features of functional MHC class II loci. Direct sequencing of the PBR encoding exon in a panel of prions revealed significantly higher levels of genetic diversity compared to two noncoding neutral loci, with most alleles differing by at least one replacement substitution in the peptide binding codons. We estimated evolutionary dynamics for Pabe-DAB1 using a variety of Bayesian and other approaches. Evidence for balancing selection comes from a spatially variable ratio of nonsynonymous-to-synonymous substitutions (mean d (N)/d (S) = 2.87) in the PBR, with sites predicted to be functionally relevant exhibiting the highest omega values. We estimate the population recombination rate to be approximately 0.3 per site per generation, indicating an important role for recombination in generating polymorphism at this locus. Pabe-DAB1 is among the few avian class II loci characterized at the genomic level and with a known intron-exon structure, a feature that greatly facilitated the amplification and sequencing of a single MHC locus in this species.

Parental relatedness and major histocompatibility effects on early embryo survivorship in Atlantic salmon

Salmon have provided key insights into the relative influence of natural and sexual selection on major histocompatibility complex (MHC) variation. Natural selection on salmon MHC genes has been demonstrated in pathogen studies, and there is evidence of MHC-based mate choice (sexual selection). We tested whether parental MHC genes affect survivorship of juvenile Atlantic salmon (Salmo salar) by quantifying the influence of parental genome-wide relatedness and MHC genotype on survivorship to the swim-up stage. Thirteen microsatellite loci were used to estimate the influence of genome-wide relatedness between parents on offspring survivorship and MHC genotypes were determined by sequencing part of the class IIbeta gene. Our results revealed no significant relationship between early offspring survivorship and genome-wide relatedness, predicted MHC heterozygosity, or MHC allelic similarity. Overall, our data are consistent with the contention that excess MHC heterozygosity in Atlantic salmon juveniles is due to sexual selection as well as differential survival of offspring due to MHC genotype.

Spatio-temporal variation in the strength and mode of selection acting on major histocompatibility complex diversity in water vole (Arvicola terrestris) metapopulations

Patterns of spatio-temporal genetic variation at a class II major histocompatibility complex (MHC) locus and multiple microsatellite loci were analysed within and between three water vole metapopulations in Scotland, UK. Comparisons of MHC and microsatellite spatial genetic differentiation, based on standardised tests between two demographically asynchronous zones within a metapopulation, suggested that spatial MHC variation was affected by balancing selection, directional selection and random genetic drift, but that the relative effects of these microevolutionary forces vary temporally. At the metapopulation level, between-year differentiation for MHC loci was significantly correlated with that of microsatellites, signifying that neutral factors such as migration and drift were primarily responsible for overall temporal genetic change at the metapopulation scale. Between metapopulations, patterns of genetic differentiation implied that, at large spatial scales, MHC variation was primarily affected by directional selection and drift. Levels of MHC heterozygosity in excess of Hardy-Weinberg expectations were consistent with overdominant balancing selection operating on MHC variation within metapopulations. However, this effect was not constant among all samples, indicating temporal variation in the strength of selection relative to other factors. The results highlight the benefit of contrasting variation at MHC with neutral markers to separate the effects of stochastic and deterministic microevolutionary forces, and add to a growing body of evidence showing that the mode and relative strength of selection acting on MHC diversity varies both spatially and temporally.

Evolutionary history of an MHC gene in two leporid species: characterisation of Mhc-DQA in the European brown hare and comparison with the European rabbit

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

Review. Sympatric, parapatric or allopatric: the most important way to classify speciation?

The most common classification of modes of speciation begins with the spatial context in which divergence occurs: sympatric, parapatric or allopatric. This classification is unsatisfactory because it divides a continuum into discrete categories, concentrating attention on the extremes, and it subordinates other dimensions on which speciation processes vary, such as the forces driving differentiation and the genetic basis of reproductive isolation. It also ignores the fact that speciation is a prolonged process that commonly has phases in different spatial contexts. We use the example of local adaptation and partial reproductive isolation in the intertidal gastropod Littorina saxatilis to illustrate the inadequacy of the spatial classification of speciation modes. Parallel divergence in shell form in response to similar environmental gradients in England, Spain and Sweden makes this an excellent model system. However, attempts to demonstrate 'incipient' and 'sympatric' speciation involve speculation about the future and the past. We suggest that it is more productive to study the current balance between local adaptation and gene flow, the interaction between components of reproductive isolation and the genetic basis of differentiation.

Ancient trans-specific polymorphism at pheromone receptor genes in basidiomycetes

In the majority of sexual organisms, reproduction occurs almost exclusively through the combination of distinct and alternate forms, called sexes or mating types. In some fungi, there can be dozens to hundreds of alternate alleles that determine compatible mating types. Such extensive polymorphism is expected to be maintained by balancing selection, and in extreme cases may give rise to trans-specific polymorphism. Here, we analyzed sequences of two pheromone receptors in the Microbotryum fungal species complex (Basidiomycota), which has only two alternate mating types. Several lines of evidence strongly suggest that the pheromone receptors are two allelic sequences acting to determine the alternate A1 and A2 mating types required for mating in Microbotryum. Phylogenetic trees of pheromone receptors in the Microbotryum species complex indicated a trans-specific polymorphism: the Microbotryum sequences from a given mating type were all more similar to the pheromone receptors of distantly related classes of fungi than to the alternate pheromone receptor in the Microbotryum species. A phylogenetic tree built using other known pheromone receptors from basidiomycetes showed that trans-specific polymorphism is widespread. The pheromone receptor alleles from Microbotryum appeared as the oldest, being at least 370 million years old. This represents the oldest known trans-specific polymorphism known in any organism so far, which may be due to the existence of sex chromosomes, obligate sexuality, mitochondrial inheritance linked to the mating type, and a highly selfing mating system in Microbotryum.

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.